The set of transverse homoclinic intersections for a saddle-focus equilibrium in the planar equilateral restricted four body problem admits certain simple homoclinic orbits which form the skeleton of the complete homoclinic intersection, or homoclinic web. In this thesis, the planar restricted four body problem is viewed as an invariant subsystem of the spatial problem, and the influence of this planar homoclinic skeleton on the spatial dynamics is studied from a numerical point of view. Starting from the vertical Lyapunov families emanating from saddle focus equilibria, we compute the stable/unstable manifolds of these spatial periodic orbits and look for intersections between these manifolds near the fundamental planar homoclinics. In this way, we are able to continue all of the basic planar homoclinic motions into the spatial problem as homoclinics for appropriate vertical Lyapunov orbits which, by the Smale Tangle theorem, suggest the existence of chaotic motions in the spatial problem. While the saddle-focus equilibrium solutions in the planar problems occur only at a discrete set of energy levels, the cycle-to-cycle homoclinics in the spatial problem are robust with respect to small changes in energy. The method uses high order Fourier-Taylor and Chebyshev series approximations in conjunction with the parameterization method, a general functional analytic framework for invariant manifolds. Tools that admit a natural notion of a-posteriori error analysis. Finally, we develop and implement a validation algorithm which we later use to obtain Theorems confirming the existence of homoclinic dynamics. This approach, known as the Radii polynomial, is a contraction mapping argument which can be applied to both the parameterized manifold and the Chebyshev arcs. When the Theorem applies, it guarantees the existence of a true solution near the approximation and it provides an upper bound on the C0 norm of the truncation error. / Includes bibliography. / Dissertation (Ph.D.)--Florida Atlantic University, 2021. / FAU Electronic Theses and Dissertations Collection

This thesis is concerned with numerical techniques for resolving and continuing closed, compact invariant manifolds in parameter-dependent dynamical systems with specific emphasis on invariant tori under flows. In the first part, we review several numerical methods of continuing invariant tori and concentrate on one choice called the ``orthogonality condition'. We show that the orthogonality condition is equivalent to another condition on the smooth level and show that they both descend from the same geometrical relationship. Then we show that for hyperbolic, periodic orbits in the plane, the linearization of the orthogonality condition yields a scalar system whose characteristic multiplier is the same as the non-unity multiplier of the orbit. In the second part, we demonstrate that one class of discretizations of the orthogonality condition for periodic orbits represents a natural extension of collocation. Using this viewpoint, we give sufficient conditions for convergence of a periodic orbit. The stability argument does not extend to higher-dimensional tori, however, and we prove that the method is unconditionally unstable for some common types of two-tori embedded in R^3 with even numbers of points in both angular directions. In the third part, we develop several numerical examples and demonstrate that the convergence properties of the method and discretization can be quite complicated. In the fourth and final part, we extend the method to the general case of p-tori in R^n in a different way from previous implementations and solve the continuation problem for a three-torus embedded in R^8.

<p>The aim of this work has been to develop new knowledge about macromodelling of microsystems. Doing that, we have followed two different approaches for generating macromodels, namely model order reduction and lumped modelling. The latter is a rather mature method that has been widely recognized and used for a relatively long period of time. Model order reduction, on the other hand, is a relatively new area still in rapid development. Due to this, the parts considering reduced order modelling is strongly biased towards methodology and concepts, whereas parts on lumped modelling are biased towards systems and devices.</p><p>In the first part of this thesis, we focus on model order reduction. We introduce some approaches for reducing model order for linear systems, and we give an example related to squeeze-film damping. We then move on to investigate model order reduction of nonlinear systems, where we present and use the concept of invariant manifolds. While the concept of invariant manifolds is general, we utilize it for reducing models. An obvious advantage of using invariant manifold theory is that it offers a conceptually clear understanding of effects and behaviour of nonlinear system.</p><p>We exemplify and investigate the accuracy of one method for identifying invariant manifolds. The example is based on an industrialized dual-axis accelerometer.</p><p>A new geometrical interpretation of external forcing, relating to invariant manifolds, is presented. We show how this can be utilized to deal with external forcing in a manner consistent with the invariance property of the manifold. The interpretation also aids in reducing errors for reduce models.</p><p>We extend the asymptotic approach in a manner that makes it possible to create design-parameter sensitive models. We investigate an industrialized dual-axis accelerometer by means of the method and demonstrate capabilities of the method. We also discuss how manifolds for nonlinear dissipative systems can be found.</p><p>Focusing on lumped modelling, we analyse a microresonator. We also discuss the two analogies that can be used to build electrical equivalents of mechanical systems. It is shown how the f → V analogy, linking velocity to voltage, is the natural choice. General properties of lumped modelling are investigated using models with varying degrees of freedom.</p><p>Finally, we analyse an electromagnetic system, intended for levitating objects, and we demonstrate the scaling effects of the system. Furthermore, we prove the intrinsic stability of the system, although the floating disc will be slightly tilted. This is the first analysis done assessing the stability criterions of such a systems. The knowledge arising from the analysis gives strong indications on how such a system can be utilized, designed, and improved.</p>

The aim of this work has been to develop new knowledge about macromodelling of microsystems. Doing that, we have followed two different approaches for generating macromodels, namely model order reduction and lumped modelling. The latter is a rather mature method that has been widely recognized and used for a relatively long period of time. Model order reduction, on the other hand, is a relatively new area still in rapid development. Due to this, the parts considering reduced order modelling is strongly biased towards methodology and concepts, whereas parts on lumped modelling are biased towards systems and devices. In the first part of this thesis, we focus on model order reduction. We introduce some approaches for reducing model order for linear systems, and we give an example related to squeeze-film damping. We then move on to investigate model order reduction of nonlinear systems, where we present and use the concept of invariant manifolds. While the concept of invariant manifolds is general, we utilize it for reducing models. An obvious advantage of using invariant manifold theory is that it offers a conceptually clear understanding of effects and behaviour of nonlinear system. We exemplify and investigate the accuracy of one method for identifying invariant manifolds. The example is based on an industrialized dual-axis accelerometer. A new geometrical interpretation of external forcing, relating to invariant manifolds, is presented. We show how this can be utilized to deal with external forcing in a manner consistent with the invariance property of the manifold. The interpretation also aids in reducing errors for reduce models. We extend the asymptotic approach in a manner that makes it possible to create design-parameter sensitive models. We investigate an industrialized dual-axis accelerometer by means of the method and demonstrate capabilities of the method. We also discuss how manifolds for nonlinear dissipative systems can be found. Focusing on lumped modelling, we analyse a microresonator. We also discuss the two analogies that can be used to build electrical equivalents of mechanical systems. It is shown how the f → V analogy, linking velocity to voltage, is the natural choice. General properties of lumped modelling are investigated using models with varying degrees of freedom. Finally, we analyse an electromagnetic system, intended for levitating objects, and we demonstrate the scaling effects of the system. Furthermore, we prove the intrinsic stability of the system, although the floating disc will be slightly tilted. This is the first analysis done assessing the stability criterions of such a systems. The knowledge arising from the analysis gives strong indications on how such a system can be utilized, designed, and improved.

This dissertation presents geometric approaches of understanding chaotic transport in phase space that is fundamental across many disciplines in physical sciences and engineering. This approach is based on analyzing phase space transport using boundaries and regions inside these boundaries in presence of perturbation. We present a geometric view of defining such boundaries and study the transport that occurs by crossing such phase space structures. The structure in two dimensional non-autonomous system is the codimension 1 stable and unstable manifolds associated with the hyperbolic fixed points. The manifolds separate regions with varied dynamical fates and their time evolution encodes how the initial conditions in a given region of phase space get transported to other regions. In the context of four dimensional autonomous systems, the corresponding structure is the stable and unstable manifolds of unstable periodic orbits which reside in the bottlenecks of energy surface. The total energy and the cylindrical (or tube) manifolds form the necessary and sufficient condition for global transport between regions of phase space. Furthermore, we adopt the geometric view to define escaping zones for avoiding transition/escape from a potential well using partial control. In this approach, the objective is two fold: finding the minimum control that is required for avoiding escape and obtaining discrete representation called disturbance of continuous noise that is present in physical sciences and engineering. In the former scenario, along with avoiding escape, the control is constrained to be smaller than the disturbance so that it can not exactly cancel out the disturbances. / Ph. D. / The prediction and control of critical events in engineering systems has been a major objective of scientific research in recent years. The multifaceted problems facing the modern society includes critical events such as spread of pathogens and pollutants in atmosphere and ocean, capsize of boats and cruise ships, space exploration and asteroid collision, to name but a few. Although, at first glance they seem to be disconnected problems in different areas of engineering and science, however, they have certain features that are inherently common. This can be studied using the abstraction of phase space which can be thought of as the universe where all possible solutions of the governing equations, derived using principles of physics, live and evolve in time. The <i>phase space</i> can be just 2D, 3D or even infinite dimensional but the critical events manifest themselves as volumes of phase space, which represent solutions at a given instant of time, get transported from one region to another due to the underlying dynamics. This mathematical abstraction is called phase space transport and studied under the umbrella of dynamical systems theory. The geometric view of the solutions that live in the phase space provides insight into the mechanisms of how the critical events occur, and the understanding of these mechanisms is useful in deciding about control strategies. A slightly different view for understanding critical events is to consider a thought experiment where a ball is rolling on a multi-well surface or potential well. As the time evolves, the ball will escape from its initial well and roll into another well, and eventually start exploring all the wells in a seemingly unpredictable way. However, these unpredictable escape/transition can be studied systematically using methods of chaos and dynamical systems. The escape/transition in a potential well implies a dramatic change in the behavior of the system, and hence the significance in prediction and control of <i>escaping dynamics</i>. The control aspect becomes more challenging due to inherent disturbance in the system that is difficult to model and we may not have the equal or more control authority to cancel those disturbances. However, we can usually estimate the maximum values of the disturbance, and try to avoid escaping from the potential well while using a smaller control. This idea is called <i>partial control of escaping dynamics</i> and can guarantee avoidance of escape for <i>ad infinitum</i>. In this doctoral research, we focus on the two mechanisms, phase space transport and escaping dynamics, by considering problems from fluid dynamics and capsize of a ship. The applications are used for numerical demonstration and evidence of the general approach in studying a large class of problems in classical physics.

This dissertation expands across three topics the geometric theory of phase space transit in the circular restricted three-body problem (CR3BP) and its generalizations. The first topic generalizes the low energy transport theory that relies on linearizing the Lagrange points in the CR3BP to time-periodic perturbations of the CR3BP, such as the bicircular problem (BCP) and the elliptic restricted three-body problem (ER3BP). The Lagrange points are no longer invariant under perturbation and are replaced by periodic orbits, which we call Lagrange periodic orbits. Calculating the monodromy matrix of the Lagrange periodic orbit and transforming into eigenbasis coordinates reveals that the transport geometry is a discrete analogue of the continuous transport geometry in the unperturbed problem. The second topic extends the theory of low energy phase space transit in periodically perturbed models using a nonlinear analysis of the geometry. This nonlinear analysis relies on calculating the monodromy tensors, which generalize monodromy matrices in order to encode higher order behavior, about the Lagrange periodic orbit. A nonlinear approximate map can be obtained which can be used to iterate initial conditions within the linear eigenbasis, providing a computationally efficient means of distinguishing transit and nontransit orbits that improves upon the predictions of the linear framework. The third topic demonstrates that the recently-discovered "arches of chaos" that stretch through the solar system, causing substantial phase space divergence for high energy particles, may be identified with the stable and unstable manifolds to the singularities of the CR3BP. We also study the arches in terms of particle orbital elements and demonstrate that the arches correspond to gravity assists in the two-body limit. / Doctor of Philosophy / Suppose that we have a spacecraft and we want to model its motion under gravity. Depending upon what trade-offs we are willing to make between accuracy and complexity, we have several options at our disposal. For example, the restricted three-body problem (R3BP) and its generalizations prove useful in many real-world situations and are rich in theoretical power despite seeming mathematically simple. The simplest restricted three-body problem is the circular restricted three-body problem (CR3BP). In the CR3BP, two masses (like a star and a planet or a planet and a moon) orbit their common center of gravity in circular orbits, while a much smaller body (like a spacecraft) moves freely, influenced by the gravitational fields that the two masses create. If we add in an extra force that acts on the spacecraft in a periodic, cycling way, the regular CR3BP becomes a periodically-perturbed CR3BP. Examples of periodically-perturbed CR3BP's include the bicircular problem (BCP), which adds in a third mass that appears to orbit the center of the system from a distance, and the elliptic restricted three-body problem (ER3BP), which allows the two masses to orbit more realistically as ellipses rather than circles. The purpose of this dissertation is to determine how to select trajectories that move spacecraft between places of interest in restricted three-body models. We generalize existing theories of CR3BP spacecraft motion to periodically-perturbed CR3BP's in the first two topics, and then we investigate some new areas of research in the unperturbed CR3BP in the third topic. We utilize numerical computations and mathematical methods to perform these analyses.

Aquesta tesi doctoral està emmarcada en el camp de l'astrodinàmica. Presenta solucions a problemes identificats en el disseny de missions que utilitzen òrbites entorn dels punts de libració, fent servir la teoria de sistemes dinàmics.El problema restringit de tres cossos és un model per estudiar el moviment d'un cos de massa infinitessimal sota l'atracció gravitatòria de dos cossos molt massius. Els cinc punts d'equilibri d'aquest model, en especial L1 i L2, han estat motiu de nombrosos estudis per aplicacions pràctiques en les últimes dècades (SOHO, Genesis...). Genèricament, qualsevol missió en òrbita al voltant del punt L2 del sistema Terra-Sol es veu afectat per ocultacions degudes a l'ombra de la Terra. Si l'òrbita és al voltant de L1, els eclipsis són deguts a la forta influència electromagnètica del Sol. D'entre els diferents tipus d'òrbites de libració, les òrbites de Lissajous resulten de la combinació de dues oscil.lacions perpendiculars. El seu principal avantatge és que les amplituds de les oscil.lacions poden ser escollides independentment i això les fa adapatables als requeriments de cada missió. La necessitat d'estratègies per evitar eclipsis en òrbites de Lissajous entorn dels punts L1 i L2 motivaren la primera part de la tesi. En aquesta part es presenta una eina per la planificació de maniobres en òrbites de Lissajous que no només serveix per solucionar el problema d'evitar els eclipsis, sinó també per trobar trajectòries de transferència entre òrbites d'amplituds diferents i planificar rendez-vous. Per altra banda, existeixen canals de baix cost que uneixen els punts L1 i L2 d'un sistema donat i representen una manera natural de transferir d'una regió de libració a l'altra. Gràcies al seu caràcter hiperbòlic, una òrbita de libració té uns objectes invariants associats: les varietats estable i inestable. Si tenim present que la varietat estable està formada per trajectòries que tendeixen cap a l'òrbita a la qual estan associades quan el temps avança, i que la varietat inestable fa el mateix però enrera en el temps, una intersecció entre una varietat estable i una d'inestable proporciona un camí asimptòtic entre les òrbites corresponents. Un mètode per trobar connexions d'aquest tipus entre òrbites planes entorn de L1 i L2 es presenta a la segona part de la tesi, i s'hi inclouen els resultats d'aplicar aquest mètode als casos dels problemes restringits Sol Terra i Terra-Lluna.La idea d'intersecar varietats hiperbòliques es pot aplicar també en la cerca de camins de baix cost entre les regions de libració del sistema Sol-Terra i Terra-Lluna. Si existissin camins naturals de les òrbites de libració solars cap a les lunars, s'obtindria una manera barata d'anar a la Lluna fent servir varietats invariants, cosa que no es pot fer de manera directa. I a l'inversa, un camí de les regions de libració lunars cap a les solars permetria, per exemple, que una estació fos col.locada en òrbita entorn del punt L2 lunar i servís com a base per donar servei a les missions que operen en òrbites de libració del sistema Sol-Terra. A la tercera part de la tesi es presenten mètodes per trobar trajectòries de baix cost que uneixen la regió L2 del sistema Terra-Lluna amb la regió L2 del sistema Sol-Terra, primer per òrbites planes i més endavant per òrbites de Lissajous, fent servir dos problemes de tres cossos acoblats. Un cop trobades les trajectòries en aquest model simplificat, convé refinar-les per fer-les més realistes. Una metodologia per obtenir trajectòries en efemèrides reals JPL a partir de les trobades entre òrbites de Lissajous en el model acoblat es presenta a la part final de la tesi. Aquestes trajectòries necessiten una maniobra en el punt d'acoblament, que és reduïda en el procés de refinat, arribant a obtenir trajectòries de cost zero quan això és possible. / This PhD. thesis lies within the field of astrodynamics. It provides solutions to problems which have been identified in mission design near libration points, by using dynamical systems theory. The restricted three body problem is a well known model to study the motion of an infinitesimal mass under the gravitational attraction of two massive bodies. Its five equilibrium points, specially L1 and L2, have been the object of several studies aimed at practical applications in the last decades (SOHO, Genesis...). In general, any mission in orbit around L2 of the Sun-Earth system is affected by occultations due to the shadow of the Earth. When the orbit is around L1, the eclipses are caused by the strong electromagnetic influence of the Sun. Among all different types of libration orbits, Lissajous type ones are the combination of two perpendicular oscillations. Its main advantage is that the amplitudes of the oscillations can be chosen independently and this fact makes Lissajous orbits more adaptable to the requirements of each particular mission than other kinds of libration motions. The need for eclipse avoidance strategies in Lissajous orbits around L1 and L2 motivated the first part of the thesis. It is in this part where a tool for planning maneuvers in Lissajous orbits is presented, which not only solves the eclipse avoidance problem, but can also be used for transferring between orbits having different amplitudes and for planning rendez-vous strategies.On the other hand, there exist low cost channels joining the L1 and L2 points of a given sistem, which represent a natural way of transferring from one libration region to the other one. Furthermore, there exist hyperbolic invariant objects, called stable and unstable manifolds, which are associated with libration orbits due to their hyperbolic character. If we bear in mind that the stable manifold of a libration orbit consists of trajectories which tend to the orbit as time goes by, and that the unstable manifold does so but backwards in time, any intersection between a stable and an unstable manifold will provide an asymptotic path between the corresponding libration orbits. A methodology for finding such asymptotic connecting paths between planar orbits around L1 and L2 is presented in the second part of the dissertation, including results for the particular cases of the Sun-Earth and Earth-Moon problems. Moreover, the idea of intersecting hyperbolic manifolds can be applied in the search for low cost paths joining the libration regions of different problems, such as the Sun-Earth and the Earth-Moon ones. If natural paths from the solar libration regions to the lunar ones was found, it would provide a cheap way of transferring to the Moon from the vicinity of the Earth, which is not possible in a direct way using invariant manifolds. And the other way round, paths from the lunar libration regions to the solar ones would allow for the placement of a station in orbit around the lunar L2, providing services to solar libration missions, for instance. In the third part of the thesis, a methodology for finding low cost trajectories joining the lunar L2 region and the solar L2 region is presented. This methodology was developed in a first step for planar orbits and in a further step for Lissajous type orbits, using in both cases two coupled restricted three body problems to model the Sun-Earth-Moon spacecraft four body problem. Once trajectories have been found in this simplified model, it is convenient to refine them to more realistic models. A methodology for obtaining JPL real ephemeris trajectories from the initial ones found in the coupled models is presented in the last part of the dissertation. These trajectories need a maneuver at the coupling point, which can be reduced in the refinement process until low cost connecting trajectories in real ephemeris are obtained (even zero cost, when possible).

In this dissertation, we demonstrate the importance of invariant objects in many areas of applied research. The areas of application we consider are chemistry, celestial mechanics and aerospace engineering, plasma physics, and coupled map lattices. In the context of chemical reactions, stable and unstable manifolds of fixed points separate regions of phase space that lead to a certain outcome of the reaction. We study how these regions change under the influence of exposing the molecules to a laser. In celestial mechanics and aerospace engineering, we compute periodic orbits and their stable and unstable manifolds for a object of negligible mass (e.g. a satellite or spacecraft) under the presence of Jupiter and two of its moons, Europa and Ganymede. The periodic orbits serve as convenient spot to place a satellite for observation purposes, and computing their stable and unstable manifolds have been used in constructing low-energy transfers between the two moons. In plasma physics, an important and practical problem is to study barriers for heat transport in magnetically confined plasma undergoing fusion. We compute barriers for which heat cannot pass through. However, such barriers break down and lead to robust partial barriers. In this latter case, heat can flow across the barrier, but at a very slow rate. Finally, infinite dimensional coupled map lattice systems are considered in a wide variety of areas, most notably in statistical mechanics, neuroscience, and in the discretization of PDEs. We assume that the interaction amont the lattice sites decays with the distance of the sites, and assume the existence of an invariant whiskered torus that is localized near a collection of lattice sites. We prove that the torus has invariant stable and unstable manifolds that are also localized near the torus. This is an important step in understanding the global dynamics of such systems and opens the door to new possible results, most notably studying the problem of energy transfer between the sites. / text

Local gamma-band (~30-100Hz) oscillations in the brain, produced by feedback inhibition on a characteristic timescale, appear in multiple areas of the brain and are associated with a wide range of cognitive functions. Some regions producing gamma also receive gamma-rhythmic input, and the interaction and coordination of these rhythms has been hypothesized to serve various functional roles. This thesis consists of three stand-alone chapters, each of which considers the response of a gamma-rhythmic neuronal circuit to input in an analytical framework. In the first, we demonstrate that several related models of a gamma-generating circuit under periodic forcing are asymptotically drawn onto an attracting invariant torus due to the convergence of inhibition trajectories at spikes and the convergence of voltage trajectories during sustained inhibition, and therefore display a restricted range of dynamics. In the second, we show that a model of a gamma-generating circuit under forcing by square pulses cannot maintain multiple stably phase-locked solutions. In the third, we show that a separation of time scales of membrane potential dynamics and synaptic decay causes the gamma model to phase align its spiking such that periodic forcing pulses arrive under minimal inhibition. When two of these models are mutually coupled, the same effect causes excitatory pulses from the faster oscillator to arrive at the slower under minimal inhibition, while pulses from the slower to the faster arrive under maximal inhibition. We also show that such a time scale separation allows the model to respond sensitively to input pulse coherence to an extent that is not possible for a simple one-dimensional oscillator. We draw on a wide range of mathematical tools and structures including return maps, saltation matrices, contraction methods, phase response formalism, and singular perturbation theory in order to show that the neuronal mechanism of gamma oscillations is uniquely suited to reliably phase lock across brain regions and facilitate the selective transmission of information.

The study of the long time behavior of nonlinear systems is not effortless, but it is very rewarding. The computation of invariant objects, in particular manifolds provide the scientist with the ability to make predictions at the frontiers of science. However, due to the presence of strong nonlinearities in many important applications, understanding the propagation of errors becomes necessary in order to quantify the reliability of these predictions, and to build sound foundations for future discoveries. This dissertation develops methods for the accurate computation of high-order polynomial approximations of stable/unstable manifolds attached to long periodic orbits in discrete time dynamical systems. For this purpose a multiple shooting scheme is applied to invariance equations for the manifolds obtained using the Parameterization Method developed by Xavier Cabre, Ernest Fontich and Rafael De La Llave in [CFdlL03a, CFdlL03b, CFdlL05]. / Includes bibliography. / Dissertation (Ph.D.)--Florida Atlantic University, 2020. / FAU Electronic Theses and Dissertations Collection