Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Space Interferometry Mission: Dynamical Observations of Galaxies (SIMDOG)

Space Interferometry Mission (SIM) will be used to obtain proper motions for a sample of 27 galaxies; the first proper motion measurements of galaxies beyond the satellite system of the Milky Way. SIM measurements lead to knowledge of the full 6-dimensional position and velocity vector of each galaxy. In conjunction with new gravitational flow models, the result will be the first total mass measurements of individual galaxies. The project, includes developnient of powerful theoretical methods for orbital calculations. This SIM study will lead to vastly improved determinations of individual galaxy masses, halo sizes, and the fractional contribution of dark matter. Astronomers have struggled to calculate the orbits of galaxies with only position and redshift information. Traditional N-body techniques are unsuitable for an analysis backward in time from a present distribution if any components of velocity or position are not very precisely known.

Shaya, Edward J.↗
NASA NTRS1986

What have we learned about galaxies through large computations?

Numerical models permit modelling processes within galaxies that are seemingly beyond the reach of physical contact or close-up viewing from different angles. N-body equations permit simple computations of the interactions of individual stars within galaxies, but the motions of these objects feature nonlinear characteristics because of the forces exerted by the complex as a whole which are not amenable to examination with Newtonian laws of force. The models have, however, revealed that thin disk galactic are not stable forms, that barlike structures rotating end-over-end are more stable structures, and that galactic contraction occurs during collisions between galaxies. Internal effects have been successfully predicted for galaxies which are part of clusters, growth rate disturbances due to the expansion of the Universe have been quantified as 2-3 percent and dynamical friction has been shown to damage galaxies which are sliding into background material.

Miller, R. H.↗
NASA NTRS1990

A simulation survey of galaxy interactions

Many carefully selected samples of interacting galaxies have been observed extensively in attempts to clarify whether interaction produces activity in galaxies. Because the sample members represent a wide range of encounter parameters and times, one can then study whether there are correlations between observable encounter features and, for example, Seyfert activity. On the other hand, in theoretical studies, simulations typically deal with either time-consuming detailed modelling of single galaxy pairs or tracing a few model encounters over time. The authors extend the observational survey approach by combining it with a simulation survey. The authors are conducting a survey of model encounters, covering the most important encounter parameters over a wide range. Some parameters, such as companion structure and initial velocity, are demonstratably less important and can be ignored in a first pass. The parameter range must be richly enough sampled so that the authors can evaluate the uniqueness of the observable morphology and velocity structure of the resulting simulated pairs to diagnose unobservable companion orbit parameters. They are using a self-gravitating polar n-body code run on the Cray X-MP at the Alabama Supercomputer Network. For each simulation, the authors have stellar and gas distributions predicted over, typically, a billion years, along with information on gas motions within the disk and any material captured by the companion or lost to the system. Features of disturbed spiral galaxies are sensitive enough to time and encounter parameters so that a match of the simulation survey results to observations can be applied as starting points to infer unobservable orbital or system parameters in actual sample members. This should enable them to examine whether interesting observed properties (Seyfert activity, nuclear star-formation rate) are functions of unobservable dynamical properties which characterize each encounter. Any correlations (or lack of some expected ones) will provide strong clues as to how or whether these phenomena are related to interactions. Aside from its use with such observed samples, this survey should greatly speed determination of initial orbital parameters for more detailed subsequent simulations of individual systems.

Byrd, Gene G.↗
NASA NTRS1993

Gas dynamic simulations of galaxy formation

Results are presented from a simulation modeling the formation of a group of galaxies in a 'standard' cold, dark matter universe with delta = 1, h sub 0 = 50 km/(s(Mpc)), baryon fraction omega sub b = 0.1 and spectrum normalization sigma sub 8 = 0.6 (bias parameter b = 1.7). Initial conditions are generated within a periodic box with comoving length 16 Mpc in a manner constrained to produce a small cluster of total mass approximately 10 exp 14 solar mass. Two sets of 643 particles are used to model the dark matter and baryon fluids. Each gas particle represents 1.08 x 10 exp -8 solar mass, implying an L* galaxy is resolved by approximately 1000 particles. The system is evolved self-consistently in three dimensions using the combined N-body/hydrodynamic scheme P3MSPH up to a final redshift z = 1. Evolving to the present is prohibited by the fact that the mean density in the simulated volume is above critical and the entire volume would be going nonlinear beyond this point, We are currently analyzing another run with somewhat poorer mass resolution which was evolved to the present.

Evrard, August E.↗
NASA NTRS1994

The formation of compact groups of galaxies. I: Optical properties

The small crossing time of compact groups of galaxies (t(sub cr)H(sub 0) approximately less than 0.02) makes it hard to understand why they are observable at all. Our dissipationless N-body simulations show that within a single rich collapsing group compact groups of galaxies continually form. The mean lifetime of a particular compact configuration if approximately 1 Gyr. On this time scale, members may merge and/or other galaxies in the loose group may join the compact configuration. In other words, compact configurations are continually replaced by new systems. The frequency of this process explains the observability of compact groups. Our model produces compact configurations (compact groups (CG's) with optical properties remarkably similar to Hickson's (1982) compact groups (HCG's): (1) CG's have a frequency distribution of members similar to that of HCG's; (2) CG's are approximately equals 10 times as dense as loose groups; (3) CG's have dynamical properties remarkably similar to those of HCG's; (4) most of the galaxy members of CG's are not merger remnants. The crucial aspect of the model is the relationship between CG's and the surrounding rich loose group. Our model predicts the frequency of occurrence of CG's. A preliminary analysis of 18 rich loose groups is consistent with the model prediction. We suggest further observational tests of the model.

Diaferio, Antonaldo↗
NASA NTRS1993

Cluster evolution as a probe of primordial density fluctuations

Although COBE's detection of large angle microwave background anisotropies fixes the amplitude of density fluctuations on length scales k exp -1 approximately = (300-6000) h(exp -1)Mpc, what is crucial for the level of large scale clustering is the amplitude of density fluctuations on scales (5-50) h(exp -1)Mpc. The level of dynamical clustering is usually parameterized by the size of the mass fluctuations in 8 h exp -1 Mpc spheres, sigma sub 8. For the cold dark matter model, COBE gives sigma sub 8 approximately = 1, while models with extra large scale power give sigma sub 8 approximately = 1/2. The most massive clusters of galaxies (greater than or approximately = 10 exp 15 solar mass) form from rare 'peak patches' found in the initial mass density distribution. Their abundance as a function of redshift is a sensitive probe of the wave number band k(exp -1) approx. (3-8) h(exp -1)Mpc, hence of sigma sub 8, and so cluster evolution can discriminate among models allowed by the COBE results. We use our Hierarchical Peaks Method, which accurately reproduces the results of P3M N-body simulations, to calculate the evolution of cluster x-ray flux counts, luminosity, and temperature functions as a function of sigma sub 8 for CDM models and those with more large scale power. We find that the EMSS and Edge et al. cluster samples support sigma sub 8 in the range from approx. 0.6-0.9, and that models with more large scale power (and hence flatter fluctuation spectra in the cluster regime) fit the x-ray bright end better.

Bond, J. Richard↗
NASA NTRS1994

Star cluster evolution with primordial binaries. 3: Effect of the Galactic tidal field

We present the results of N-body simulations of tidally limited star clusters with an initial population of 0%-20% binaries. We find that (1) if enough binaries are initially present, the binary fraction may fall to a minimum value, then increase at late times; (2) the cluster evaporation timescale is quite insensitive to the details of the initial binary distribution; (3) the cluster core radius stabilizes at a few percent of the half-mass radius when binaries are present, just as in the case of isolated clusters; and (4) there may be a marked difference between the spatial distribution of low-energy and high-energy binaries as the cluster evolves. Specifically, the spatial distribution of the lower energy systems is often substantially more extended than that of the more tightly bound pairs. At no time are our simulated clusters well described by simple dynamical models that neglect the close coupling between the binding energies and the center-of-mass energies of the binaries they contain.

Mcmillan, Steve↗
DOE OSTI2026

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↗
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 NTRS1994

A simulation of the intracluster medium with feedback from cluster galaxies

We detail method and report first results from a three-dimensional hydrodynamical and N-body simulation of the formation and evolution of a Coma-sized cluster of galaxies, with the intent of studying the history of the hot, X-ray emitting intracluster medium. Cluster gas, galaxies, and dark matter are included in the model. The galaxies and dark matter fell gravitational forces; the cluster gas also undergoes hydrodynamical effects such as shock heating and PdV work. For the first time in three dimensions, we include modeling of ejection of processed gas from the simulated galaxies by winds, including heating and heavy element enrichment. For comparison, we employ a `pure infall' simulation using the same initial conditions but with no galaxies or winds. We employ an extreme ejection history for galactic feedback in order to define the boundary of likely models. As expected, feedback raises the entropy of the intracluster gas, preventing it from collapsing to densities as high as those attained in the infall model. The effect is more pronounced in subclusters formed at high redshift. The cluster with feedback is always less X-ray luminous, but experiences more rapid luminosity evolution, than the pure infall cluster. Even employing an extreme ejection model, the final gas temperature is only approximately 15% larger than in the infall model. The radial temperature profile is very nearly isothermal within 1.5 Mpc. The cluster galaxies in the feedback model have a velocity dispersion approximately 15% lower than the dark matter. This results in the true ratio of specific energies in galaxies to gas being less than one, beta(sub spec) approximately 0.7. The infall model predicts beta(sub spec) approximately 1.2. Large excursions in these values occur over time, following the complex dynamical history of the cluster. The morphology of the X-ray emission is little affected by feedback. The emission profiles of both clusters are well described by the standard beta-model with beta(sub fit) approximately equal to 0.7 - 0.9. X-ray mass estimates based on the assumptions of hydrostatic equilibrium and the applicability of the beta-model are quite accurate in both cases. A strong, radial iron abundance gradient is present, which develops as a consequence of the steepening of the galaxy density profile over time. Spectroscopic observations using nonimaging detectors with wide (approximately 45 min) fields of view dramatically smear the gradient. Observations with arcminute resolution, made available with the ASCA satellite, would readily resolve the gradient.

Metzler, Christopher A.↗
DOE OSTI2025

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↗
DOE OSTI2024

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↗
NASA NTRS1973

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.↗
NASA NTRS1973

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.↗
NASA NTRS1978

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.↗
NASA NTRS1982

Galaxy clusters with multiple components. I - The dynamics of Abell 98

The nature and evolution of rich galaxy clusters containing multiple condensations is studied. Velocities for 13 members of A98 are presented and the determination of dynamical properties of the individual subclusters is substantially improved. The application of techniques for quantifying gravitational scale lengths are stressed without invocation of the symmetry arguments demanded by the usual parameter fitting procedures. The entire A98 cluster is modeled as a two-body system. The probability that A98 is bound is found to be 98 percent. From this model and the results of previous N-body simulations, it is shown that the most likely description of the system is that it reached maximum expansion about 3.5 billion years ago and is presently in a state of collapse. The subclusters will merge in another three billion years.

Beers, T. C.↗
NASA NTRS2007

Newton-Euler Dynamic Equations of Motion for a Multi-body Spacecraft

The Magnetospheric MultiScale (MMS) mission employs a formation of spinning spacecraft with several flexible appendages and thruster-based control. To understand the complex dynamic interaction of thruster actuation, appendage motion, and spin dynamics, each spacecraft is modeled as a tree of rigid bodies connected by spherical or gimballed joints. The method presented facilitates assembling by inspection the exact, nonlinear dynamic equations of motion for a multibody spacecraft suitable for solution by numerical integration. The building block equations are derived by applying Newton's and Euler's equations of motion to an "element" consisting of two bodies and one joint (spherical and gimballed joints are considered separately). Patterns in the "mass" and L'force" matrices guide assembly by inspection of a general N-body tree-topology system. Straightforward linear algebra operations are employed to eliminate extraneous constraint equations, resulting in a minimum-dimension system of equations to solve. This method thus combines a straightforward, easily-extendable, easily-mechanized formulation with an efficient computer implementation.

Stoneking, Eric↗
NASA NTRS2015

The 2014 KCG Meteor Outburst: Clues to a Parent Body

The Kappa Cygnid (KCG) meteor shower exhibited unusually high activity in 2014, producing ten times the typical number of meteors. The shower was detected in both radar and optical systems and meteoroids associated with the outburst spanned at least five decades in mass. In total, the Canadian Meteor Orbit Radar, European Network, and NASA All Sky and Southern Ontario Meteor Network produced thousands of KCG meteor trajectories. Using these data, we have undertaken a new and improved characterization of the dynamics of this little-studied, variable meteor shower. The Cygnids have a di use radiant and a significant spread in orbital characteristics, with multiple resonances appearing to play a role in the shower dynamics. We conducted a new search for parent bodies and found that several known asteroids are orbitally similar to the KCGs. N-body simulations show that the two best parent body candidates readily transfer meteoroids to the Earth in recent centuries, but neither produces an exact match to the KCG radiant, velocity, and solar longitude. We nevertheless identify asteroid 2001 MG1 as a promising parent body candidate.

Moorhead, Althea V.↗