Source author record

Alessandra Celletti

Alessandra Celletti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2021arXiv

Proper Elements for Space Debris

Proper elements are quasi-invariants of a Hamiltonian system, obtained through a normalization procedure. Proper elements have been successfully used to identify families of asteroids, sharing the same dynamical properties. We show that proper elements can also be used within space debris dynamics to identify groups of fragments associated to the same break-up event. The proposed method allows to reconstruct the evolutionary history and possibly to associate the fragments to a parent body. The procedure relies on different steps: (i) the development of a model for an approximate, though accurate, description of the dynamics of the space debris; (ii) the construction of a normalization procedure to determine the proper elements; (iii) the production of fragments through a simulated break-up event. We consider a model that includes the Keplerian part, an approximation of the geopotential, and the gravitational influence of Sun and Moon. We also evaluate the contribution of Solar radiation pressure and the effect of noise on the orbital elements. We implement a Lie series normalization procedure to compute the proper elements associated to semi-major axis, eccentricity and inclination. Based upon a wide range of samples, we conclude that the distribution of the proper elements in simulated break-up events (either collisions and explosions) shows an impressive connection with the dynamics observed immediately after the catastrophic event. The results are corroborated by a statistical data analysis based on the check of the Kolmogorov-Smirnov test and the computation of the Pearson correlation coefficient.

preprint2020arXiv

KAM estimates for the dissipative standard map

From the beginning of KAM theory, it was realized that its applicability to realistic problems depended on developing quantitative estimates on the sizes of the perturbations allowed. In this paper we present results on the existence of quasi-periodic solutions for conformally symplectic systems in non-perturbative regimes. We recall that, for conformally symplectic systems, finding the solution requires also to find a "drift parameter". We present a proof on the existence of solutions for values of the parameters which agree with more than three figures with the numerically conjectured optimal values. The verification of the estimates presented here is not completely rigorous since we do not control the round-off error. Nevertheless, running with different precision shows very little difference in the results. Given the high precision of the calculation and the simplicity of the estimates, this does not seem to affect the results. A full verification should be done implementing interval arithmetic. We make available the approximate solutions, the highly efficient algorithms to generate them (incorporating high precision based on the MPFR library) and the routines used to verify the applicability of the theorem.

preprint2020arXiv

KAM theory for some dissipative systems

Dissipative systems play a very important role in several physical models, most notably in Celestial Mechanics, where the dissipation drives the motion of natural and artificial satellites, leading them to migration of orbits, resonant states, etc. Hence the need to develop theories that ensure the existence of structures such as invariant tori or periodic orbits and device efficient computational methods. In this work we concentrate on the existence of invariant tori for the specific case of dissipative systems known as "conformally symplectic" systems, which have the property that they transform the symplectic form into a multiple of itself. To give explicit examples of conformally symplectic systems, we will present two different models: a discrete system, the standard map, and a continuous system, the spin-orbit problem. In both cases we will consider the conservative and dissipative versions, that will help to highlight the differences between the symplectic and conformally symplectic dynamics. For such dissipative systems we will present a KAM theorem in an a-posteriori format. The method of proof is based on extending geometric identities originally developed in [39] for the symplectic case. Besides leading to streamlined proofs of KAM theorem, this method provides a very efficient algorithm which has been implemented. Coupling an efficient numerical algorithm with an a-posteriori theorem, we have a very efficient way to provide rigorous estimates close to optimal. Indeed, the method gives a criterion (the Sobolev blow up criterion) that allows to compute numerically the breakdown. We will review this method as well as an extension of J. Greene's method and present the results in the conservative and dissipative standard maps. Computing close to the breakdown, allows to discover new mathematical phenomena such as the "bundle collapse mechanism".

preprint2019arXiv

Existence of whiskered KAM tori of conformally symplectic systems

We study the existence of whiskered tori in a family $f_μ$ of conformally symplectic maps depending on parameters $μ$. Whiskered tori are tori on which the motion is a rotation, but they have as many expanding/contracting directions as allowed by the preservation of the geometric structure. Our main result is formulated in an "a-posteriori" format. We fix $ω$ satisfying Diophantine conditions. We assume that we are given 1) a value of the parameter $μ_0$, 2) an embedding of the torus $K_0$ into the phase space, approximately invariant under $f_{μ_0}$ in the sense that $f_{μ_0} \circ K_0 - K_0 \circ T_ω$ is small, 3) a splitting of the tangent space at the range of $K_0$, into three bundles which are approximately invariant under $D f_{μ_0}$ and such that the derivative satisfies "rate conditions" on each of the components. Then, if some non-degeneracy conditions are satisfied, we show that there is another parameter $μ_\infty$, an embedding $K_\infty$ and splittings close to the original ones which are invariant under $f_{μ_\infty}$. We also bound $|μ_\infty - μ_0|$, $\|K_\infty - K_0 \|$ and the distance of the initial and final splittings in terms of the initial error. The proof of the main theorem consists in describing an iterative process that takes advantage of cancellations coming from the geometry. Then, we show that the process converges to a true solution when started from an approximate enough solution. The iterative process leads to an efficient algorithm that is quite practical to implement. As an application, we study the singular problem of effects of small dissipation on whiskered tori. We develop formal expansions in the perturbative parameter and use them as input for the a-posteriori theorem. This allows to obtain lower bounds for the domain of analyticity of the tori as function of the perturbative parameter.

preprint2019arXiv

Resonances in the Earth's Space Environment

We study the presence of resonances in the region of space around the Earth. We consider a massless body (e.g, a dust particle or a small space debris) subject to different forces: the gravitational attraction of the geopotential, the effects of Sun and Moon. We distinguish different types of resonances: tesseral resonances are due to a commensurability involving the revolution of the particle and the rotation of the Earth, semi-secular resonances include the rates of variation of the mean anomalies of Moon and Sun, while secular resonances just depend on the rates of variation of the arguments of perigee and the longitudes of the ascending nodes of the perturbing bodies. We characterize such resonances, giving precise statements on the regions where the resonances can be found and provide examples of some specific commensurability relations.

preprint2016arXiv

A study of the lunisolar secular resonance $2\dotω+\dotΩ=0$

The dynamics of small bodies around the Earth has gained a renewed interest, since the awareness of the problems that space debris can cause in the nearby future. A relevant role in space debris is played by lunisolar secular resonances, which might contribute to an increase of the orbital elements, typically of the eccentricity. We concentrate our attention on the lunisolar secular resonance described by the relation $2\dotω+\dotΩ=0$, where $ω$ and $Ω$ denote the argument of perigee and the longitude of the ascending node of the space debris. We introduce three different models with increasing complexity. We show that the growth in eccentricity, as observed in space debris located in the MEO region at the inclination about equal to $56^\circ$, can be explained as a natural effect of the secular resonance $2\dotω+\dotΩ=0$, while the chaotic variations of the orbital parameters are the result of interaction and overlapping of nearby resonances.

preprint2016arXiv

Analytical development of the lunisolar disturbing function and the critical inclination secular resonance

We provide a detailed derivation of the analytical expansion of the lunar and solar disturbing functions. Although there exist several papers on this topic, many derivations contain mistakes in the final expansion or rather (just) in the proof, thereby necessitating a recasting and correction of the original derivation. In this work, we provide a self-consistent and definite form of the lunisolar expansion. We start with Kaula's expansion of the disturbing function in terms of the equatorial elements of both the perturbed and perturbing bodies. Then we give a detailed proof of Lane's expansion, in which the elements of the Moon are referred to the ecliptic plane. Using this approach the inclination of the Moon becomes nearly constant, while the argument of perihelion, the longitude of the ascending node, and the mean anomaly vary linearly with time. We make a comparison between the different expansions and we profit from such discussion to point out some mistakes in the existing literature, which might compromise the correctness of the results. As an application, we analyze the long--term motion of the highly elliptical and critically--inclined Molniya orbits subject to quadrupolar gravitational interactions. The analytical expansions presented herein are very powerful with respect to dynamical studies based on Cartesian equations, because they quickly allow for a more holistic and intuitively understandable picture of the dynamics.

preprint2016arXiv

Bifurcations of lunisolar secular resonances for space debris orbits

Using bifurcation theory, we study the secular resonances induced by the Sun and Moon on space debris orbits around the Earth. In particular, we concentrate on a special class of secular resonances, which depend only on the debris' orbital inclination. This class is typically subdivided into three distinct types of secular resonances: those occurring at the critical inclination, those corresponding to polar orbits, and a third type resulting from a linear combination of the rates of variation of the argument of perigee and the longitude of the ascending node. The model describing the dynamics of space debris includes the effects of the geopotential, as well as the Sun's and Moon's attractions, and it is defined in terms of suitable action-angle variables. We consider the system averaged over both the mean anomaly of the debris and those of the Sun and Moon. Such multiply-averaged Hamiltonian is used to study the lunisolar resonances which depend just on the inclination. Borrowing the technique from the theory of bifurcations of Hamiltonian normal forms, we study the birth of periodic orbits and we determine the energy thresholds at which the bifurcations of lunisolar secular resonances take place. This approach gives us physically relevant information on the existence and location of the equilibria, which help us to identify stable and unstable regions in the phase space. Besides their physical interest, the study of inclination dependent resonances offers interesting insights from the dynamical point of view, since it sheds light on different phenomena related to bifurcation theory.

preprint2016arXiv

Geostationary secular dynamics revisited: application to high area-to-mass ratio objects

The long-term dynamics of the geostationary Earth orbits (GEO) is revisited through the application of canonical perturbation theory. We consider a Hamiltonian model accounting for all major perturbations: geopotential at order and degree two, lunisolar perturbations with a realistic model for the Sun and Moon orbits, and solar radiation pressure. The long-term dynamics of the GEO region has been studied both numerically and analytically, in view of the relevance of such studies to the issue of space debris or to the disposal of GEO satellites. Past studies focused on the orbital evolution of objects around a nominal solution, hereafter called the forced equilibrium solution, which shows a particularly strong dependence on the area-to-mass ratio. Here, we i) give theoretical estimates for the long-term behavior of such orbits, and ii) we examine the nature of the forced equilibrium itself. In the lowest approximation, the forced equilibrium implies motion with a constant non-zero average `forced eccentricity', as well as a constant non-zero average inclination, otherwise known in satellite dynamics as the inclination of the invariant `Laplace plane'. Using a higher order normal form, we demonstrate that this equilibrium actually represents not a point in phase space, but a trajectory taking place on a lower-dimensional torus. We give analytical expressions for this special trajectory, and we compare our results to those found by numerical orbit propagation. We finally discuss the use of proper elements, i.e., approximate integrals of motion for the GEO orbits.

preprint2016arXiv

Poynting-Robertson drag and solar wind in the space debris problem

We analyze the combined effect of Poynting-Robertson and solar wind drag on space debris. We derive a model within Cartesian, Gaussian and Hamiltonian frameworks. We focus on the geosynchronous resonance, although the results can be easily generalized to any resonance. By numerical and analytical techniques, we compute the drift in semi-major axis due to Poynting-Robertson and solar wind drag. After a linear stability analysis of the equilibria, we combine a careful investigation of the regular, resonant, chaotic behavior of the phase space with a long-term propagation of a sample of initial conditions. The results strongly depend on the value of the area-to-mass ratio of the debris, which might show different dynamical behaviors: temporary capture or escape from the geosynchronous resonance, as well as temporary capture or escape from secondary resonances involving the rate of variation of the longitude of the Sun. Such analysis shows that Poynting-Robertson and solar wind drag must be taken into account, when looking at the long-term behavior of space debris. Trapping or escape from the resonance can be used to place the debris in convenient regions of the phase space.

preprint2016arXiv

The theory of secondary resonances in the spin-orbit problem

We study the resonant dynamics in a simple one degree of freedom, time dependent Hamiltonian model describing spin-orbit interactions. The equations of motion admit periodic solutions associated with resonant motions, the most important being the synchronous one in which most evolved satellites of the Solar system, including the Moon, are observed. Such primary resonances can be surrounded by a chain of smaller islands which one refers to as secondary resonances. Here, we propose a novel canonical normalization procedure allowing to obtain a higher order normal form, by which we obtain analytical results on the stability of the primary resonances as well as on the bifurcation thresholds of the secondary resonances. The procedure makes use of the expansion in a parameter, called the detuning, measuring the shift from the exact secondary resonance. Also, we implement the so-called `book-keeping' method, i.e., the introduction of a suitable separation of the terms in orders of smallness in the normal form construction, which deals simultaneously with all the small parameters of the problem. Our analytical computation of the bifurcation curves is in excellent agreement with the results obtained by a numerical integration of the equations of motion, thus providing relevant information on the parameter regions where satellites can be found in a stable configuration.

preprint2015arXiv

A study of the main resonances outside the geostationary ring

We investigate the dynamics of satellites and space debris in external resonances, namely in the region outside the geostationary ring. Precisely, we focus on the 1:2, 1:3, 2:3 resonances, which are located at about 66 931.4 km, 87 705.0 km, 55 250.7 km, respectively. Some of these resonances have been already exploited in space missions, like XMM-Newton and Integral. Our study is mainly based on a Hamiltonian approach, which allows us to get fast and reliable information on the dynamics in the resonant regions. Significative results are obtained even by considering just the effect of the geopotential in the Hamiltonian formulation. For objects (typically space debris) with high area-to-mass ratio the Hamiltonian includes also the effect of the solar radiation pressure. In addition, we perform a comparison with the numerical integration in Cartesian variables, including the geopotential, the gravitational attraction of Sun and Moon, and the solar radiation pressure. We implement some simple mathematical tools that allows us to get information on the terms which are dominant in the Fourier series expansion of the Hamiltonian around a given resonance, on the amplitude of the resonant islands and on the location of the equilibrium points. We also compute the Fast Lyapunov Indicators, which provide a cartography of the resonant regions, yielding the main dynamical features associated to the external resonances. We apply these techniques to analyze the 1:2, 1:3, 2:3 resonances; we consider also the case of objects with large area-to-mass ratio and we provide an application to the case studies given by XMM-Newton and Integral.

preprint2015arXiv

Domains of analyticity of Lindstedt expansions of KAM tori in dissipative perturbations of Hamiltonian systems

Many problems in Physics are described by dynamical systems that are conformally symplectic (e.g., mechanical systems with a friction proportional to the velocity, variational problems with a small discount or thermostated systems). Conformally symplectic systems are characterized by the property that they transform a symplectic form into a multiple of itself. The limit of small dissipation, which is the object of the present study, is particularly interesting. We provide all details for maps, but we present also the modifications needed to obtain a direct proof for the case of differential equations. We consider a family of conformally symplectic maps $f_{μ, ε}$ defined on a $2d$-dimensional symplectic manifold $\mathcal M$ with exact symplectic form $Ω$; we assume that $f_{μ,ε}$ satisfies $f_{μ,ε}^*Ω=λ(ε) Ω$. We assume that the family depends on a $d$-dimensional parameter $μ$ (called drift) and also on a small scalar parameter $ε$. Furthermore, we assume that the conformal factor $λ$ depends on $ε$, in such a way that for $ε=0$ we have $λ(0)=1$ (the symplectic case). We study the domains of analyticity in $ε$ near $ε=0$ of perturbative expansions (Lindstedt series) of the parameterization of the quasi--periodic orbits of frequency $ω$ (assumed to be Diophantine) and of the parameter $μ$. Notice that this is a singular perturbation, since any friction (no matter how small) reduces the set of quasi-periodic solutions in the system. We prove that the Lindstedt series are analytic in a domain in the complex $ε$ plane, which is obtained by taking from a ball centered at zero a sequence of smaller balls with center along smooth lines going through the origin. The radii of the excluded balls decrease faster than any power of the distance of the center to the origin.

preprint2015arXiv

Dynamical investigation of minor resonances for space debris

We study the dynamics of the space debris in regions corresponding to minor resonances; precisely, we consider the resonances 3:1, 3:2, 4:1, 4:3, 5:1, 5:2, 5:3, 5:4, where a j:l resonance (with j, l integers) means that the periods of revolution of the debris and of rotation of the Earth are in the ratio j/l. We consider a Hamiltonian function describing the effect of the geopotential and we use suitable finite expansions of the Hamiltonian for the description of the different resonances. In particular, we determine the leading terms which dominate in a specific orbital region, thus limiting our computation to very few harmonics. Taking advantage from the pendulum-like structure associated to each term of the expansion, we are able to determine the amplitude of the islands corresponding to the different harmonics. By means of simple mathematical formulae, we can predict the occurrence of splitting or overlapping of the resonant islands for different values of the parameters. We also find several cases which exhibit a transcritical bifurcation as the inclination is varied. These results, which are based on a careful mathematical analysis of the Hamiltonian expansion, are confirmed by a numerical study of the dynamical behavior obtained by computing the so-called Fast Laypunov Indicators. Since the Hamiltonian approach includes just the effect of the geopotential, we validate our results by performing a numerical integration in Cartesian variables of a more complete model including the gravitational attraction of Sun and Moon, as well as the solar radiation pressure.

preprint2015arXiv

Halo orbits around the collinear points of the restricted three-body problem

We perform an analytical study of the bifurcation of the halo orbits around the collinear points $L_1$, $L_2$, $L_3$ for the circular, spatial, restricted three--body problem. Following a standard procedure, we reduce to the center manifold constructing a normal form adapted to the synchronous resonance. Introducing a detuning, which measures the displacement from the resonance and expanding the energy in series of the detuning, we are able to evaluate the energy level at which the bifurcation takes place for arbitrary values of the mass ratio. In most cases, the analytical results thus obtained are in very good agreement with the numerical expectations, providing the bifurcation threshold with good accuracy. Care must be taken when dealing with $L_3$ for small values of the mass-ratio between the primaries; in that case, the model of the system is a singular perturbation problem and the normal form method is not particularly suited to evaluate the bifurcation threshold.

preprint2015arXiv

Qualitative and analytical results of the bifurcation thresholds to halo orbits

We study the dynamics in the neighborhood of the collinear Lagrangian points in the spatial, circular, restricted three--body problem. We consider the case in which one of the primaries is a radiating body and the other is oblate (although the latter is a minor effect). Beside having an intrinsic mathematical interest, this model is particularly suited for the description of a mission of a spacecraft (e.g., a solar sail) to an asteroid. The aim of our study is to investigate the occurrence of bifurcations to halo orbits, which take place as the energy level is varied. The estimate of the bifurcation thresholds is performed by analytical and numerical methods: we find a remarkable agreement between the two approaches. As a side result, we also evaluate the influence of the different parameters, most notably the solar radiation pressure coefficient, on the dynamical behavior of the model. To perform the analytical and numerical computations, we start by implementing a center manifold reduction. Next, we estimate the bifurcation values using qualitative techniques (e.g. Poincaré surfaces, frequency analysis, FLIs). Concerning the analytical approach, following \cite{CPS} we implement a resonant normal form, we transform to suitable action-angle variables and we introduce a detuning parameter measuring the displacement from the synchronous resonance. The bifurcation thresholds are then determined as series expansions in the detuning. Three concrete examples are considered and we find in all cases a very good agreement between the analytical and numerical results.

preprint2015arXiv

Response solutions for quasi-periodically forced, dissipative wave equations

We consider several models of nonlinear wave equations subject to very strong damping and quasi-periodic external forcing. This is a singular perturbation, since the damping is not the highest order term. We study the existence of response solutions (i.e., quasi-periodic solutions with the same frequency as the forcing). Under very general non-resonance conditions on the frequency, we show the existence of asymptotic expansions of the response solution; moreover, we prove that the response solution indeed exists and depends analytically on $\varepsilon$ (where $\varepsilon$ is the inverse of the coefficient multiplying the damping) for $\varepsilon$ in a complex domain, which in some cases includes disks tangent to the imaginary axis at the origin. In other models, we prove analyticity in cones of aperture $π/2$ and we conjecture it is optimal. These results have consequences for the asymptotic expansions of the response solutions considered in the literature. The proof of our results relies on reformulating the problem as a fixed point problem, constructing an approximate solution and studying the properties of iterations that converge to the solutions of the fixed point problem.

preprint2014arXiv

On the dynamics of space debris: 1:1 and 2:1 resonances

We study the dynamics of the space debris in the 1:1 and 2:1 resonances, where geosynchronous and GPS satellites are located. By using Hamiltonian formalism, we consider a model including the geopotential contribution for which we compute the secular and resonant expansions of the Hamiltonian. Within such model we are able to detect the equilibria and to study the main features of the resonances in a very effective way. In particular, we analyze the regular and chaotic behavior of the 1:1 and 2:1 resonant regions by analytical methods and by computing the Fast Lyapunov Indicators, which provide a cartography of the resonances. This approach allows us to detect easily the location of the equilibria, the amplitudes of the libration islands and the main dynamical stability features of the resonances, thus providing an overview of the 1:1 and 2:1 resonant domains under the effect of Earth's oblateness. The results are validated by a comparison with a model developed in Cartesian coordinates, including the geopotential, the gravitational attraction of Sun and Moon and the solar radiation pressure.

preprint2014arXiv

The effect of Poynting-Robertson drag on the triangular Lagrangian points

We investigate the stability of motion close to the Lagrangian equilibrium points L4 and L5 in the framework of the spatial, elliptic, restricted three- body problem, subject to the radial component of Poynting-Robertson drag. For this reason we develop a simplified resonant model, that is based on averaging theory, i.e. averaged over the mean anomaly of the perturbing planet. We find temporary stability of particles displaying a tadpole motion in the 1:1 resonance. From the linear stability study of the averaged simplified resonant model, we find that the time of temporary stability is proportional to beta a1 n1 , where beta is the ratio of the solar radiation over the gravitational force, and a1, n1 are the semi-major axis and the mean motion of the perturbing planet, respectively. We extend previous results (Murray (1994)) on the asymmetry of the stability indices of L4 and L5 to a more realistic force model. Our analytical results are supported by means of numerical simulations. We implement our study to Jupiter-like perturbing planets, that are also found in extra-solar planetary systems.

preprint2012arXiv

Stability estimates of nearly-integrable systems with dissipation and non-resonant frequency

We consider a dissipative vector field which is represented by a nearly-integrable Hamiltonian flow to which a non symplectic force is added, so that the phase space volume is not preserved. The vector field depends upon two parameters, namely the perturbing and dissipative parameters, and by a drift function. We study the general case of an l-dimensional, time-dependent vector field. Assuming to start with non-resonant initial conditions, we prove the stability of the variables which are actions of the conservative system (namely, when the dissipative parameter is set to zero) for exponentially long times. In order to construct the normal form, a suitable choice of the drift function must be performed. We also provide some simple examples in which we construct explicitly the normal form, we make a comparison with a numerical integration and we compute theoretical bounds on the parameters as well as we give explicit stability estimates.

preprint2012arXiv

Stability of nearly-integrable systems with dissipation

We study the stability of a vector field associated to a nearly-integrable Hamiltonian dynamical system to which a dissipation is added. Such a system is governed by two parameters, named the perturbing and dissipative parameters, and it depends on a drift function. Assuming that the frequency of motion satisfies some resonance assumption, we investigate the stability of the dynamics, and precisely the variation of the action variables associated to the conservative model. According to the structure of the vector field, one can find linear and exponential stability times, which are established under smallness con- ditions on the parameters. We also provide some applications to concrete examples, which exhibit a linear or exponential stability behavior.