Optimal hybrid propulsion systems thrust trajectories in a patched-conic N-body force field
Euler-Lagrange equations for approximating optimal thrust trajectories for limited hybrid propulsion system in N-body dynamic environment
Pines, S.↗
SEARCH · Search NASA
Results for “n-body dynamics”
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.
One Million Open-source Cislunar Orbits
Cislunar space, encompassing the region from geosynchronous orbit to beyond the Moon, is poised to become a cornerstone for future exploration, scientific discovery, and national security. Missions in this region, spanning durations from weeks to decades, require robust infrastructure and reliable transit capabilities. The complex gravitational influences of the Moon, Sun, and planets, along with thermal radiation from Earth and the Sun, lead to significant trajectory deviations, resulting in kilometer-scale errors within days. Leveraging the high-performance computing resources at Lawrence Livermore National Laboratory (LLNL), we have simulated one million high-fidelity cislunar trajectories, now publicly available via LLNL’s Green Data Oasis and the Unified Data Library. Generated using the open-source Space Situational Awareness Python package, these trajectories match the precision of commercial tools such as AGI’s Systems Tool Kit and NASA’s General Mission Analysis Tool. This data set is a valuable resource for reference, statistical analysis of cislunar orbit populations, and training machine learning models for rapid orbit classification with minimal observational input. Preliminary analysis reveals stable bands in Keplerian element space, particularly around five geosynchronous radii across a range of inclinations and eccentricities. Beyond this threshold, the Moon’s influence disrupts most unassisted orbits, though co-orbiting L4/L5 Lunar Trojans persist throughout the six-year simulation.
Astronomy and AstroPhysics↗
Pseudoatom molecular dynamics plasma microfields
Spectral lines are powerful diagnostic tools for both laboratory and astrophysical plasmas, as their shape is sensitive to the plasma environment. The low-frequency component of the electric microfield is an important input for semi-analytic line broadening codes. Here, in this paper, we detail a new method of calculating plasma microfields using configuration-resolved pseudoatom molecular dynamics. This approach accounts for both quantum atomic structure and N-body effects, similar to density functional theory molecular dynamics, but with less computational cost. We present pseudoatom microfields at plasma conditions relevant for recent high energy density laboratory astrophysics experiments conducted at the Sandia Z-Machine, National Ignition Facility, and Linac Coherent Light Source. Compared to established microfield codes we find moderate deviations at solid density conditions and strong agreement at lower plasma densities.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Equations of motion for coupled n-body systems
Computer program, developed to analyze spacecraft attitude dynamics, can be applied to large class of problems involving objects that can be simplified into component parts. Systems of coupled rigid bodies, point masses, symmetric wheels, and elastically flexible bodies can be analyzed. Program derives complete set of non-linear equations of motion in vectordyadic format. Numerical solutions may be printed out. Program is in FORTRAN IV for batch execution and has been implemented on IBM 360.
Frisch, H. P.↗
An Open Benchmark of One Million High-Fidelity Cislunar Trajectories
Cislunar space spans from geosynchronous altitudes to beyond the Moon and will underpin future exploration, science, and security operations. We describe and release an open dataset of one million numerically propagated cislunar trajectories generated with the open-source Space Situational Awareness Python package (SSAPy). The model includes high-degree Earth/Moon gravity, solar gravity, and Earth/Sun radiation pressure; other planetary gravities are omitted by design for computational efficiency. Initial conditions uniformly sample commonly used osculating-element ranges, and each trajectory is propagated for up to six years under a single, fixed start epoch. The dataset is intended as a reusable benchmark for method development (e.g., space domain awareness, navigation, and machine-learning pipelines), a reference library for statistical studies of orbit families, and a starting point for community-driven extensions (e.g., alternative epochs). We report empirically observed stability trends (e.g., a band near ~5 GEO and persistence of some co-orbital classes including L4/L5 librators) as dataset descriptors rather than new dynamical results. The chief contribution is the scale, fidelity, organization (CSV/HDF5 with full state time series and metadata), and open availability, which together lower the barrier to comparative and data-driven studies in the cislunar regime.
79 ASTRONOMY AND ASTROPHYSICS↗
A suboptimal stochastic controller for an N-body spacecraft
Considerable attention, in the open literature, is being focused on the problem of developing a suitable set of deterministic dynamical equations for a complex spacecraft. This paper considers the problem of determining a stochastic optimal controller for an n-body spacecraft. The approach used in obtaining the stochastic controller involves the application, interpretation, and combination of advanced dynamical principles and the theoretical aspects of modern control theory. The stochastic controller obtained herein for a complicated model of a spacecraft uses sensor angular measurements associated with the base body to obtain smoothed estimates of the entire state vector. It can be easily implemented, and it enables system performance to be significantly improved.
Larson, V.↗
State equations for an n-body spacecraft
Considerable attention, in the open literature, is being focused on the problem of developing a suitable set of deterministic dynamical equations for a complex spacecraft. The present paper addresses the problem of determining a set of state equations for an n-body spacecraft. The approach used in obtaining the state equations involves the application and interpretation of advanced dynamical principles. This set of state equations can be effectively used in the development of a stochastic controller for the spacecraft (in this latter development, the deterministic model developed in the present paper will be appropriately corrupted by plant noise). The major effort in the paper revolves around the determination of the plant matrices and the specification of the state vector and the control vector.
Larson, V.↗
Numerical experiments on galaxy formation. I - Introduction and first results
The large, fully three-dimensional n-body programs designed for numerical experiments on the dynamics of galaxies have been used for investigations of protogalactic collapses. The dynamical aspects of galaxy formation dominate in this approach. The present experiments start from a large rotating gaseous mass that is made up of many gas clouds. The large mass collapses because it is out of equilibrium. Stars form as the collapse proceeds. The problem is formulated, limiting-case examples are described, and results from simple 'first cut' models are given. Results with no star formation are contrasted with results for very rapid star formation. All models with gas formed thin disklike structures; some were circular, while others were oval. The principal result from the simple models is that stars continue to form rapidly even after the gas has been depleted to unacceptably low levels (1% of the galactic mass in gas).
Miller, R. H.↗
Six collapses
The self-consistent dynamical development of six stellar systems, started from rotating spherical configurations, has been studied by means of a fully three-dimensional n-body integration. The six examples had different initial angular velocities and velocity dispersions. All settled down into prolate bars rotating about a short axis within two initial rotation periods. The bars are long-lived, robust, and stable. Bars are the natural form toward which rapidly rotating stellar dynamical systems develop, instead of the flattened axisymmetric disks that had been expected. The early stages of each collapse are reasonably well described by a theoretical model according to which a collapse passes through a sequence of rigidly rotating, uniform-density spheroids. The first significant departures from spheroidal form were axisymmetric in all cases. Rings formed in some examples, sheets in others, with transition cases between these extremes. Nonaxisymmetric forms developed from these intermediate stages.
Miller, R. H.↗
Suboptimal stochastic controller for an n-body spacecraft
The problem is studied of determining a stochastic optimal controller for an n-body spacecraft. The approach used in obtaining the stochastic controller involves the application, interpretation, and combination of advanced dynamical principles and the theoretical aspects of modern control theory. The stochastic controller obtained for a complicated model of a spacecraft uses sensor angular measurements associated with the base body to obtain smoothed estimates of the entire state vector, can be easily implemented, and enables system performance to be significantly improved.
Larson, V.↗
On a Critical Acceleration Scale of Dark Matter in ΛCDM and Dynamical Dark Energy
Abstract Universal acceleration a 0 emerges in various empirical laws, yet its fundamental nature remains unclear. Using Illustris and Virgo N -body simulations, we focus on the velocity and acceleration fluctuations in collisionless dark matter involving long-range gravity. For comparison, in the kinetic theory of gases, molecules undergo random elastic collisions involving short-range interactions, where only velocity fluctuations are relevant. Hierarchical structure formation proceeds through the merging of smaller halos to form larger halos, which facilitates a continuous energy cascade from small to large halos at a constant rate ε u ≈ −10 −7 m 2 s −3 . Velocity fluctuations involve a critical velocity u c ∝ (1 + z ) −3/4 . Acceleration fluctuations involve a critical acceleration a c ∝ (1 + z ) 3/4 . Two critical quantities are related by the rate of energy cascade ε u ≈ − a c u c /[2(3 π ) 2 ], where factor 3 π is from the angle of incidence during merging. With critical velocity u c on the order of 300 km s −1 at z = 0, the critical acceleration is determined to be a c 0 ≡ a c ( z = 0) ≈ 10 −10 m s −2 , suggesting a c might explain the universal acceleration a 0 ≈ 10 −10 m s −2 in the empirical Tully–Fisher relation or modified Newtonian dynamics. The redshift evolution a c ∝ (1 + z ) 3/4 is in good agreement with Magneticum and EAGLE simulations and in reasonable agreement with limited observations. This suggests a larger a 0 at a higher redshift such that galaxies of fixed mass rotate faster at a higher redshift. Note that for dark energy (DE) density ρ DE 0 ≈ a c 0 2 / G = 1 0 − 10 J m −3 , we postulate an entropic origin of the DE from acceleration fluctuations of dark matter, analogous to the gas pressure from velocity fluctuations. This leads to a dynamical DE coupled to the structure evolution involving a relatively constant DE density followed by a slow weakening phase, suggesting possible deviations from the standard ΛCDM paradigm.
N-body simulations↗
ELG×LRG Distribution through Dark Matter Halo Dynamics
We investigate the clustering and halo occupation distribution (HOD) of DESI Y1 emission-line (ELGs) and luminous red (LRGs) galaxies at 0.8 < z < 1.1, including their cross-correlation (ELG×LRG), using the A BACUS S UMMIT suite and a new Halo Occupation Model (H OME ) for galaxy multitracers. This integrates intrahalo dynamics, halo exclusion, and quenching, bridging insights from hydrodynamical, HOD, abundance-matching, and semianalytic studies. Leveraging full phase-space information from the Uchuu N-body simulation, and sampling satellites from dark-matter particle positions via physically motivated prescriptions, Home reproduces the anisotropic clustering down to s = 200 h −1 kpc with unprecedented accuracy. Model parameters are inferred solely from two-point statistics using a two-level Bayesian framework, yielding high-fidelity ELG, LRG, and cross-reference catalogs. We find that satellite ELGs behave as incoherent flows within their parent halos, dominating the clustering below 4 h −1 Mpc. The HOD from the best-fit Home has the following properties: (i) 90.50% (85.91%) of ELGs (LRGs) are central galaxies without satellites, residing in halos of M vir ∼ 6.6 × 10 11 (1.2 × 10 13 ) h −1 M ⊙ ; (ii) the ELG×LRG cross-correlation is governed by central-central pairs and shaped by halo exclusion on 2–5 h −1 Mpc scales; (iii) 9.50% (14.09%) of ELGs (LRGs) are satellites, of which 1.09% (3.52%) inhabit halos with a central galaxy of the same species in a maximally conformal configuration, 7.02% (0.005%) orbit complementary hosts in a minimally conformal state, and 0.58% (10.57%) are orphans. The high sensitivity of Home precisely captures the dynamics of satellites in different host environments, opening a promising avenue for understanding systematics and the dynamical nature of dark matter, potentially distinguishing gravity models.
Favole, Ginevra [Universidad de La Laguna (Spain);↗
Reconstructing the Stripping History of the Sagittarius Stream with Neural Networks
The Sagittarius (Sgr) Stream is produced by the ongoing disruption of the Sgr dwarf spheroidal (dSph) galaxy and is thought to contain multiple wraps that were stripped during different pericentric passages. In this study, we introduce a neural-network–based method trained on N-body simulations to infer the stripping time of Sgr Stream stars directly from their phase-space coordinates. We combine spectroscopic data from SEGUE, APOGEE DR17, and LAMOST DR7 low-resolution spectroscopic (LRS) survey with Gaia EDR3 astrometry and distance estimates from the latest StarHorse catalog to identify high-quality Sgr Stream members. Applying our method to these stars, we measure a clear metallicity gradient with stripping time, well described by a linear relation with slope ∼0.3 dex Gyr −1 . We further predict the stripping times of globular clusters previously suggested to originate from the Sgr dSph. M 54, Terzan 7, Terzan 8, and Arp 2 exhibit stripping times consistent with being currently bound to the Sgr remnant. Pal 12, Whiting 1, and NGC 2419 are inferred to have been stripped 0.9 ± 0.1, 1.1 ± 0.2, and 2.1 ± 0.2 Gyr ago, respectively. For NGC 4147 and NGC 5634, whose membership in the Sgr system remains uncertain, our analysis suggests stripping times of 1.1 ± 0.4 and 1.1 ± 0.1 Gyr, respectively, if they are ultimately confirmed as genuine Sgr members. These results demonstrate that data-driven models of dynamical stripping histories offer a promising approach for reconstructing the formation and chemical evolution of the Sgr Stream.
79 ASTRONOMY AND ASTROPHYSICS↗
Effects of Dark Matter Self-interactions on Sagittarius and its Stream
This work explores how assumptions regarding the particle-physics nature of dark matter can alter the evolution of the Sagittarius (Sgr) dwarf spheroidal galaxy and its expansive stellar stream. We run a large suite of N-body simulations to model the infall of a spherically symmetric Sgr-like dwarf, exploring how the presence of dark matter self-interactions impacts its evolution. For a scattering cross section of σ/m χ = 30 cm 2 g −1 (at orbital velocity scales), these interactions result in significantly less stellar mass and little to no dark matter bound to the progenitor at the present day. To isolate the cause of this mass loss, we introduce a novel technique for controlling which pairs of dark matter simulation particles can interact. This enables us to identify ram-pressure evaporation—the scattering of satellite and host dark matter particles—as the primary source of the enhanced mass loss. The rapid disintegration of the Sgr progenitor when self-interactions are allowed alters some key properties of the resulting stellar stream, most dramatically suppressing the presence of a “spur” on the apocenter of the trailing stream arm that correlates with the mass of the satellite at last pericenter. We demonstrate how the effects on the Sgr system scale with the particular choice of self-interaction cross section, which affects the degree of ram-pressure evaporation. These findings generalize beyond the Sgr system, underscoring that dwarf stellar streams and dwarf galaxies with close passages may serve as sensitive probes for dark matter self-interactions.
Milky Way dark matter halo↗
Generalized representative structures for atomistic systems
A new method is presented to generate atomic structures that reproduce the essential characteristics of arbitrary material systems, phases, or ensembles. Previous methods allow one to reproduce the essential characteristics (e.g. the chemical disorder) of a large random alloy within a small crystal structure. The ability to generate small representations of random alloys, along with the restriction to crystal systems, results from using the fixed-lattice cluster correlations to describe structural characteristics. A more general description of the structural characteristics of atomic systems is obtained using complete sets of atomic environment descriptors. These are used within for generating representative atomic structures without restriction to fixed lattices. A general data-driven approach is provided here utilizing the atomic cluster expansion (ACE) basis. The N-body ACE descriptors are a complete set of atomic environment descriptors that span both chemical and spatial degrees of freedom and are used within for describing atomic structures. The generalized representative structure (GRS) method presented within generates small atomic structures that reproduce ACE descriptor distributions corresponding to arbitrary structural and chemical complexity. It is shown that systematically improvable representations of crystalline systems on fixed parent lattices, amorphous materials, liquids, and ensembles of atomic structures may be produced efficiently through optimization algorithms. With the GRS method, we highlight reduced representations of atomistic machine-learning training datasets that contain similar amounts of information and small 40–72 atom representations of liquid phases. The ability to use GRS methodology as a driver for informed novel structure generation is also demonstrated. The advantages over other data-driven methods and state-of-the-art methods restricted to high-symmetry systems are highlighted.
atomic cluster expansion↗
Dynamical models for a spacecraft idealized as a set of multi-hinged rigid bodies
A brief description is presented of a canonical set of equations which governs the behavior of an n-body spacecraft. General results are given for the case in which the spacecraft is modeled in terms of n rigid bodies connected by dissipative elastic joints. The final equations are free from constraint torques and involve only r variables (r is the number of degrees of freedom of the system). An advantage which accompanies the elimination of the constraint torques is a decrease in the computer run time (especially when n is large).
Larson, V.↗
Hybrid-coordinate spacecraft dynamics using large-deformation modal coordinates
For a vehicle amenable to idealization as a set of n + 1 rigid bodies interconnected by n line hinges (implying tree topology), a set of minimum-dimension matrix equations of motion is recorded in terms of variables including the relative rotations of contiguous bodies. When a substructure or portion of the vehicle is subjected to large gross deformations but is restricted to 'small' strains or local deformations, the n-body model is (for sufficiently large n) assumed to experience 'small' relative rotations of contiguous bodies within the substructure, and the equations of motion are linearized in these variables. The application of coordinate transformations and truncations to the linearized subset of the system equations of motion then produces a hybrid-coordinate formulation of the problem involving both discrete and distributed coordinates.
Likins, P. W.↗
Unified connected theory of few-body reaction mechanisms in N-body scattering theory
A unified treatment of different reaction mechanisms in nonrelativistic N-body scattering is presented. The theory is based on connected kernel integral equations that are expected to become compact for reasonable constraints on the potentials. The operators T/sub +-//sup ab/(A) are approximate transition operators that describe the scattering proceeding through an arbitrary reaction mechanism A. These operators are uniquely determined by a connected kernel equation and satisfy an optical theorem consistent with the choice of reaction mechanism. Connected kernel equations relating T/sub +-//sup ab/(A) to the full T/sub +-//sup ab/ allow correction of the approximate solutions for any ignored process to any order. This theory gives a unified treatment of all few-body reaction mechanisms with the same dynamic simplicity of a model calculation, but can include complicated reaction mechanisms involving overlapping configurations where it is difficult to formulate models.
Polyzou, W. N.↗