Search NASA⌕ Search

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.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

On the collapse and violent relaxation of protoglobular clusters

During the formation of stellar systems such as globular clusters, low-mass subcondensations which eventually form stars must retain a geometric size throughout the collapse process that is small compared to the characteristic distance separating them. If the local velocity dispersion of the subcondensations is small, the overall dimension of the system can decrease substantially before reaching a dynamical equilibrium state. The maximum collapse factor is deduced by examining the growth of the velocity dispersion and the spread in arrival times at the origin caused by local and global fluctuations. It is shown, analytically as well as in a series of N-body simulations, that the maximum reduction in the characteristic dimension of a system of N fragments with an initial homogeneous distribution subject to N exp 1/2 fluctuations is proportional to N exp 1/3. Direct physical collisions between low-mass subcondensations are therefore unlikely to occur in protoglobular clusters. The results are discussed in the context of fragmentation and violent relaxation.

Aarseth, S. J.↗
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 NTRS1984

Numerical experiments on galaxy clustering in open universes

Numerical studies were performed on the dynamical effects on the evolution of clumps, filaments, voids, and galaxy clusters by various final Omega values. The final Omega values examined ranged from 0.03-1, and attention was given to defining observations of superclusters which would aid in determining an actual value for Omega. The numerical trials consisted of n-body integration programs governed by the total expansion and final Omega value and included comparisons between results for open and closed universes. All runs started from the same initial conditions. The dispersion trajectories of particles and the final galactic cluster forms were found to be equivalent, regardless of the final Omega value. The possibility of deriving a value for Omega from velocity data on galaxies which have not yet joined clusters is discussed.

Miller, R. H.↗
NASA NTRS1993

The spatial correlation properties of dark galaxy halos in a CDM universe

We use the Hierarchical Particle Mesh (HPM) N-body code written by J. V. Villumsen (Villumsen, 1989) to investigate the two-point spatial correlation function, xi(r), of dark galaxy halos as a function of halo mass and local environment (i.e. high, low, or average mass density). We assume a standard cold dark matter (CDM) universe (omega = 1, delta = 0, H sub 0 = 50,km/sec/Mpc). Because of the large dynamic ranges in mass and length that can be obtained with the HPM code, it is well-suited to an investigation of this sort.

Brainerd, Tereasa G.↗
NASA NTRS1983

Numerical simulations of the decay of satellite galaxy orbits

A multiple three-body technique is used to study the orbital evolution of satellite galaxies which is similar to the N-body method but neglects two-body forces between stars in the halo of the parent galaxy. It is found that, for satellites orbiting within the halo, Chandrasekhar's (1960) dynamical friction formula accurately describes the orbital decay rate, including its variation with satellite mass and size and with the number density and mass of halo stars. Significant frictional forces are present even outside the halo, and the orbital decay rate, instead of depending on the procedure used to place the satellite in its orbit, is determined only by the current orbital parameters. This semirestricted N-body method is sufficiently fast to have permitted the running of 200 simulations to date, many more than would have been possible by means of the conventional N-body technique.

Lin, D. N. C.↗
NASA NTRS1993

A test of the adhesion approximation for gravitational clustering

We quantitatively compare a particle implementation of the adhesion approximation to fully non-linear, numerical 'N-body' simulations. Our primary tool, cross-correlation of N-body simulations with the adhesion approximation, indicates good agreement, better than that found by the same test performed with the Zel-dovich approximation (hereafter ZA). However, the cross-correlation is not as good as that of the truncated Zel-dovich approximation (TZA), obtained by applying the Zel'dovich approximation after smoothing the initial density field with a Gaussian filter. We confirm that the adhesion approximation produces an excessively filamentary distribution. Relative to the N-body results, we also find that: (a) the power spectrum obtained from the adhesion approximation is more accurate than that from ZA or TZA, (b) the error in the phase angle of Fourier components is worse than that from TZA, and (c) the mass distribution function is more accurate than that from ZA or TZA. It appears that adhesion performs well statistically, but that TZA is more accurate dynamically, in the sense of moving mass to the right place.

Melott, Adrian L.↗
NASA NTRS1994

A test of the adhesion approximation for gravitational clustering

We quantitatively compare a particle implementation of the adhesion approximation to fully nonlinear, numerical 'N-body' simulations. Our primary tool, cross-correlation of N-body simulations with the adhesion approximation, indicates good agreement, better than that found by the same test performed with the Zel'dovich approximation (hereafter ZA). However, the cross-correlation is not as good as that of the truncated Zel'dovich approximation (TZA), obtained by applying the Zel'dovich approximation after smoothing the initial density field with a Gaussian filter. We confirm that the adhesion approximation produces an excessively filamentary distribution. Relative to the N-body results, we also find that: (a) the power spectrum obtained from the adhesion approximation is more accurate that that from ZA to TZA, (b) the error in the phase angle of Fourier components is worse that that from TZA, and (c) the mass distribution function is more accurate than that from ZA or TZA. It appears that adhesion performs well statistically, but that TZA is more accurate dynamically, in the sense of moving mass to the right place.

Melott, Adrian L.↗
NASA NTRS1993

Effects of gas on the global stability of galactic disks - Radial flows

We study numerically the effect of gas on the global stability of a two-component self-gravitating galactic disk embedded in a live halo. The stars are evolved by using a 3D collisionless N-body code, and the gas is represented by an ensemble of finite size inelastic particles. The gravitational interaction of stars and gas is calculated using a TREE method. We find that the evolution of the gaseous distribution in the globally unstable disks can be described by two different regimes. When the gas mass fraction is less than about 10 percent, the gas is channeled toward the galactic center by a growing stellar bar. For higher gas fractions, the gas becomes highly inhomogeneous, and the bar instability in the disk is heavily damped. The gas falls toward the inner kpc due to dynamical friction. Domains of both regimes depend on the efficiency of dissipation in the gas. We also discuss the relevance of the Jeans instability and give an empirical criterion for the global bar instability in a two-component self-gravitating disk.

Shlosman, Isaac↗
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↗
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.↗
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.↗