Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Symplectic integration”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

16 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

A mapping method for the gravitational few-body problem with dissipation

Recently a new class of numerical integration methods - 'mixed variable symplectic integrators' - has been introduced for studying long-term evolution in the conservative gravitational few-body problem. These integrators are an order of magnitude faster than conventional ordinary differential equations (ODE) integration methods. Here we present a simple modification of this method to include small non-gravitational forces. The new scheme provides a similar advantage of computational speed for a larger class of problems in Solar System dynamics.

Malhotra, Renu↗
NASA NTRS1997

Modelling the Centers of Galaxies

The key to studying central regions by means of nobody numerical experiments is to concentrate on the central few parsecs of a galaxy, replacing the remainder of the galaxy by a suitable boundary condition, rather after the manner in which stellar interiors can be studied without a detailed stellar atmosphere by replacing the atmosphere with a boundary condition. Replacements must be carefully designed because the long range gravitational force means that the core region is sensitive to mass outside that region and because particles can exchange between the outer galaxy and the core region. We use periodic boundary conditions, coupled with an iterative procedure to generate initial particle loads in isothermal equilibrium. Angular momentum conservation is ensured for problems including systematic rotation by a circular reflecting boundary and by integrating in a frame that rotates with the mean flow. Mass beyond the boundary contributes to the gravitational potential, but does not participate in the dynamics. A symplectic integration scheme has been developed for rotating coordinate systems. This combination works well, leading to robust configurations. Some preliminary results with this combination show that: (1) Rotating systems are extremely sensitive to non-axisymmetric external potentials, and (2) that a second core, orbiting near the main core (like the M31 second core system), shows extremely rapid orbital decay. The experimental setups will be discussed, along with preliminary results.

Smith, B. F.↗
NASA NTRS1991

Symplectic maps for the n-body problem

The present study generalizes the mapping method of Wisdom (1982) to encompass all gravitational n-body problems with a dominant central mass. The rationale for the generalized mapping method is discussed as well as details for the mapping for the n-body problem. Some refinements of the method are considered, and the relationship of the mapping method to other symplectic integration methods is shown. The method is used to compute the evolution of the outer planets for a billion years. The resulting evolution is compared to the 845 million year evolution of the outer planets performed on the Digital Orerry using standard numerical integration techniques. This calculation provides independent numerical confirmation of the result of Sussman and Wisdom (1988) that the motion of the planet Pluto is chaotic.

Wisdom, Jack↗
NASA NTRS1992

Experimenting with galaxies

A study to demonstrate how the dynamics of galaxies may be investigated through the creation of galaxies within a computer model is presented. The numerical technique for simulating galaxies is shown to be both highly efficient and highly robust. Consideration is given to the anatomy of a galaxy, the gravitational N-body problem, numerical approaches to the N-body problem, use of the Poisson equation, and the symplectic integrator.

Miller, Richard H.↗
NASA NTRS2003

Terrestrial Planet Formation in Binary Star Systems

Most stars reside in multiple star systems; however, virtually all models of planetary growth have assumed an isolated single star. Numerical simulations of the collapse of molecular cloud cores to form binary stars suggest that disks will form within such systems. Observations indirectly suggest disk material around one or both components within young binary star systems. If planets form at the right places within such circumstellar disks, they can remain in stable orbits within the binary star systems for eons. We are simulating the late stages of growth of terrestrial planets within binary star systems, using a new, ultrafast, symplectic integrator that we have developed for this purpose. We show that the late stages of terrestrial planet formation can indeed take place in a wide variety of binary systems and we have begun to delineate the range of parameter space for which this statement is true. Results of our initial simulations of planetary growth around each star in the alpha Centauri system and other 'wide' binary systems, as well as around both stars in very close binary systems, will be presented.

Lissauer, Jack J.↗
NASA NTRS2003

Terrestrial Planet Formation Around Close Binary Stars

Most stars reside in multiple star systems; however, virtually all models of planetary growth have assumed an isolated single star. Numerical simulations of the collapse of molecular cloud cores to form binary stars suggest that disks will form within such systems. Observations indirectly suggest disk material around one or both components within young binary star systems. If planets form at the right places within such circumstellar disks, they can remain in stable orbits within the binary star systems for eons. We are simulating the late stages of growth of terrestrial planets around close binary stars, using a new, ultrafast, symplectic integrator that we have developed for this purpose. The sum of the masses of the two stars is one solar mass, and the initial disk of planetary embryos is the same as that used for simulating the late stages of terrestrial planet growth within our Solar System and in the Alpha Centauri wide binary star system. Giant planets &are included in the simulations, as they are in most simulations of the late stages of terrestrial planet accumulation in our Solar System. When the stars travel on a circular orbit with semimajor axis of up to 0.1 AU about their mutual center of mass, the planetary embryos grow into a system of terrestrial planets that is statistically identical to those formed about single stars, but a larger semimajor axis and/or a significantly eccentric binary orbit can lead to significantly more dynamically hot terrestrial planet systems.

Lissauer, Jack J.↗
NASA NTRS2008

The Detectability of Exo-Earths and Super-Earths via Resonant Signatures in Exozodiacal Clouds

Directly imaging extrasolar terrestrial planets necessarily means contending with the astrophysical noise of exozodiacal dust and the resonant structures created by these planets in exozodiacal clouds. Using a custom tailored hybrid symplectic integrator we have constructed 120 models of resonant structures created by exo-Earths and super-Earths on circular orbits interacting with collisionless steady-state dust clouds around a Sun-like star. Our models include enough particles to overcome the limitations of previous simulations that were often dominated by a handful of long-lived particles, allowing us to quantitatively study the contrast of the resulting ring structures. We found that in the case of a planet on a circular orbit, for a given star and dust source distribution, the morphology and contrast of the resonant structures depend on only two parameters: planet mass and (square root)ap/Beta, where ap is the planet's semi-major axis and Beta is the ratio of radiation pressure force to gravitational force on a grain. We constructed multiple-grain-size models of 25,000 particles each and showed that in a collisionless cloud, a Dohnanyi crushing law yields a resonant ring whose optical depth is dominated by the largest grains in the distribution, not the smallest. We used these models to estimate the mass of the lowest-mass planet that can be detected through observations of a resonant ring for a variety of assumptions about the dust cloud and the planet's orbit. Our simulations suggest that planets with mass as small as a few times Mars' mass may produce detectable signatures in debris disks at ap greater than or approximately equal to 10 AU.

Stark, Christopher C.↗
NASA NTRS2012

A Survey of Symplectic and Collocation Integration Methods for Orbit Propagation

Demands on numerical integration algorithms for astrodynamics applications continue to increase. Common methods, like explicit Runge-Kutta, meet the orbit propagation needs of most scenarios, but more specialized scenarios require new techniques to meet both computational efficiency and accuracy needs. This paper provides an extensive survey on the application of symplectic and collocation methods to astrodynamics. Both of these methods benefit from relatively recent theoretical developments, which improve their applicability to artificial satellite orbit propagation. This paper also details their implementation, with several tests demonstrating their advantages and disadvantages.

integration↗
NASA NTRS1999

Dynamical Chaos in the Wisdom-Holman Integrator: Origins and Solutions

We examine the nonlinear stability of the Wisdom-Holman (WH) symplectic mapping applied to the integration of perturbed, highly eccentric (e-0.9) two-body orbits. We find that the method is unstable and introduces artificial chaos into the computed trajectories for this class of problems, unless the step size chosen 1s small enough that PeriaPse is always resolved, in which case the method is generically stable. This 'radial orbit instability' persists even for weakly perturbed systems. Using the Stark problem as a fiducial test case, we investigate the dynamical origin of this instability and argue that the numerical chaos results from the overlap of step-size resonances; interestingly, for the Stark-problem many of these resonances appear to be absolutely stable. We similarly examine the robustness of several alternative integration methods: a time-regularized version of the WH mapping suggested by Mikkola; the potential-splitting (PS) method of Duncan, Levison, Lee; and two original methods incorporating approximations based on Stark motion instead of Keplerian motion. The two fixed point problem and a related, more general problem are used to conduct a comparative test of the various methods for several types of motion. Among the algorithms tested, the time-transformed WH mapping is clearly the most efficient and stable method of integrating eccentric, nearly Keplerian orbits in the absence of close encounters. For test particles subject to both high eccentricities and very close encounters, we find an enhanced version of the PS method-incorporating time regularization, force-center switching, and an improved kernel function-to be both economical and highly versatile. We conclude that Stark-based methods are of marginal utility in N-body type integrations. Additional implications for the symplectic integration of N-body systems are discussed.

Rauch, Kevin P.↗
NASA NTRS2002

Comet and Asteroid Hazard to the Terrestrial Planets

We made computer simulations of orbital evolution for intervals of at least 5-10 Myr of N=2000 Jupiter-crossing objects (JCOs) with initial orbits close to those of real comets with period P less than 10 yr, 500 objects with orbits close to that of Comet 10P, and the asteroids initially located at the 3:1 and 5:2 resonances with Jupiter at initial eccentricity e(sub 0)=0.15 and initial inclination i(sub 0)=10(sup 0). The gravitational influence of all planets, except for Mercury and Pluto, was taken into account (without dissipative factors). We calculated the probabilities of collisions of bodies with the terrestrial planets, using orbital elements obtained with a step equal to 500 yr, and then summarized the results for all bodies, obtaining, the total probability Psigma of collisions with a planet and the total time interval Tsigma during which perihelion distance q of bodies was less than a semimajor axis of the planet. The values of p(sub r) =10(exp 6)Psigma/N and T(sub r)=T/1000 yr (where T=Tsigma/N) are presented in a table together with the ratio r of the total time interval when orbits were of Apollo type (at a greater than 1 AU, q less than 1.017 AU, e less than 0.999) to that of Amor type (1.017 less than q less than 1.33 AU), r(sub 2) is the same as r but for Apollo objects with e less than 0.9. For asteroids we present only results obtained by direct integration, as a symplectic method can give large errors for these resonances.

Ipatov, S. I.↗
NASA NTRS1991

Mapping the stability field of Jupiter Trojans

Jupiter Trojans are a remnant of outer solar system planetesimals captured into stable or quasistable libration about the 1:1 resonance with the mean motion of Jupiter. The observed swarms of Trojans may provide insight into the original mass of condensed solids in the zone from which the Jovian planets accumulated, provided that the mechanisms of capture can be understood. As the first step toward this understanding, the stability field of Trojans were mapped in the coordinate proper eccentricity, e(sub p), and libration amplitude, D. To accomplish this mapping, the orbits of 100 particles with e(sub p) in the range of 0 to 0.8 and D in the range 0 to 140 deg were numerically integrated. Orbits of the Sun, the four Jovian planets, and the massless particles were integrated as a full N-body system, in a barycentric frame using fourth order symplectic scheme.

Levison, H. F.↗
NASA NTRS2002

Solar System Dynamics

In these 18 years, the research has touched every major dynamical problem in the solar system, including: the effect of chaotic zones on the distribution of asteroids, the delivery of meteorites along chaotic pathways, the chaotic motion of Pluto, the chaotic motion of the outer planets and that of the whole solar system, the delivery of short period comets from the Kuiper belt, the tidal evolution of the Uranian arid Galilean satellites, the chaotic tumbling of Hyperion and other irregular satellites, the large chaotic variations of the obliquity of Mars, the evolution of the Earth-Moon system, and the resonant core- mantle dynamics of Earth and Venus. It has introduced new analytical and numerical tools that are in widespread use. Today, nearly every long-term integration of our solar system, its subsystems, and other solar systems uses algorithms that was invented. This research has all been primarily Supported by this sequence of PGG NASA grants. During this period published major investigations of tidal evolution of the Earth-Moon system and of the passage of the Earth and Venus through non-linear core-mantle resonances were completed. It has published a major innovation in symplectic algorithms: the symplectic corrector. A paper was completed on non-perturbative hydrostatic equilibrium.

Wisdom, Jack↗
NASA NTRS1992

Polynomial approximation of Poincare maps for Hamiltonian system

Different methods are proposed and tested for transforming a non-linear differential system, and more particularly a Hamiltonian one, into a map without integrating the whole orbit as in the well-known Poincare return map technique. We construct piecewise polynomial maps by coarse-graining the phase-space surface of section into parallelograms and using either only values of the Poincare maps at the vertices or also the gradient information at the nearest neighbors to define a polynomial approximation within each cell. The numerical experiments are in good agreement with both the real symplectic and Poincare maps.

Froeschle, Claude↗