SEARCH · Engineering Papers

Results for “N-body simulation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1979

Collapse and relaxation of rotating stellar systems

The evolution of initially spherical, uniform-density, and uniform-velocity-dispersion stellar systems is investigated by means of a three-dimensional N-body computer model. Some 100,000 simulation stars are used to follow the collapse and relaxation of the system for various amounts of angular momentum in solid-body rotation. For initially low values of the angular momentum satisfying the Ostriker-Pebbles stability criterion the systems quickly relax to an axisymmetric shape and resemble elliptical galaxies in appearance. The maximum flattening for these systems is equivalent only to an E2 system. For larger values of the initial angular momentum bars develop, and the systems undergo a much more drastic evolution. The apparent rotational and random velocities of the barred systems are very sensitive to the viewing direction. An additional complication is the frequent misalignment of the apparent major axis with the direction that reflects the maximum rotation.

Hohl, F.↗

NASA NTRS1994

Structure of merger remnants. 4: Isophotal shapes

This paper examines the shapes of isophotes of galaxy merger remnants. More specifically, we perform a series of numerical experiments to study galaxy mergers. The simulations explore a variety of encounter geometries, types of progenitor galaxies, and particle numbers. We 'observe' each of the remnants from 64 viewpoints to estimate how the isophotal shapes vary with the orientation of the remnant. Also, by comparing the results from the various simulations, we learn how encounter geometry and the structure of the progenitors can affect the shape of the remnants. The encounter geometry, the structure of the progenitors, and the orientation of the remnant play deciding roles in the shape of the 'observed' isophotes, so much so that these simulated merger remnants have both 'boxy' and 'disky' isophotes and ellipticities ranging from E0 to E7. Thus, the shapes seen in our merger simulation remnants span a similar range to observed isophotal shapes of real galaxies. Furthermore, to estimate the errors in this analysis, we introduce the statistical bootstrap. Bootstrapping is used rather widely in observational astronomy when dealing with small samples; however, its use is practically unheard of for tackling theoretical problems. We believe that bootstrapping can be extremely useful for dealing with the small samples found in numerical simulations; therefore, we present the rudiments and basis of the technique with emphasis on its use in N-body calculations. Also, by comparing the errors for the smaller simulations with those of the larger ones, we speculate on the number of particles required to accurately explore isophotal shapes in simulations.

Heyl, Jeremy S.↗

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

The power spectrum of IRAS galaxies

A technique is presented for removing the power spectrum in flux-limited galaxy catalogs is presented and tested using Monte Carlo simulations of the IRAS survey based on N-body realizations of a standard CDM universe. The technique is applied to the 1.2 Jy IRAS catalog, comparing the resulting IRAS power spectrum with previous results in the literature. A qualitative comparison of the IRAS power spectrum with various theoretical models is shown. An extensive set of Monte Carlo experiments designed to test quantitatively the validity of standard CDM on large scales is described.

Fisher, Karl B.↗

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 NTRS1994

The stability of perfect elliptic disks. 1: The maximum streaming case

Self-consistent distribution functions are constructed for two-dimensional perfect elliptic disks (for which the potential is exactly integrable) in the limit of maximum streaming; these are tested for stability by N-body integration. To obtain a discrete representation for each model, simulated annealing is used to choose a set of orbits which sample the distribution function and reproduce the required density profile while carrying the greatest possible amount of angular momentum. A quiet start technique is developed to place particles on these orbits uniformly in action-angle space, making the initial conditions as smooth as possible. The roundest models exhibit spiral instabilities similar to those of cold axisymmetric disks; the most elongated models show bending instabilities like those seen in prolate systems. Between these extremes, there is a range of axial ratios 0.25 approximately less than b/a approximately less than 0.6 within which these models appear to be stable. All the methods developed in this investigation can easily be extended to integrable potentials in three dimensions.

Levine, Stephen E.↗

NASA NTRS1993

Evolution of planetesimals. I - Dynamics: Relaxation in a thin disk. II - Numerical simulations

The study examines the effects of density inhomogeneity and differential rotation as well as inelastic collisions on the dynamical evolution of planetesimals. Consideration is given to a three-step analysis: the dynamical evolution of the planetesimals, collisions and mass accumulation, and interaction with gas. It is shown that the velocity dispersion of a cold system of planetesimals increases rapidly due to elastic gravitational scattering. When the dispersion in the epicycle amplitude becomes comparable to the planetesimals' Roche radius, energy is transferred from the systematic Keplerian shear to the dispersive motion. With a numerical N-body scheme, gravitational scattering and physical collisions among a system of planetesimals is simulated. It is shown that dynamical equilibrium is attained with a velocity dispersion comparable to the surface escape velocity of those planetesimals which contribute most of the system mass.

Palmer, P. L.↗

NASA NTRS1993

A new method to simulate vertical and horizontal structure in galactic disks

We have modified the particles in an N-body treecode to have different softening lengths in the horizontal and vertical directions. This allows us to simultaneously have thin enough particles to resolve the vertical structure in galactic disks, and horizontally large enough particles to suppress the vertical heating due to two-body effects.

Dalcanton, Julianne J.↗

NASA NTRS1991

Orbital decay in aspherical galaxies. II - Triaxial systems

A previous analysis of orbital decay of small satellite galaxies into moderately oblate, axisymmetrical cannibal galaxies with Staeckel potentials is generalized here to triaxial systems with arbitary anisotropy and internal streaming. The orbital evolution to be expected in systems with differing degrees of anisotropy or rotation is briefly discussed. The evolution of some representative orbits is presented in detail using the orbit-averaged anisotropic Chandrasekhar drag and the distribution functions for perfect ellipsoids. The validity of orbit averaging is assessed by comparing the results with ensembles of single-particle orbits integrated directly in the same potential and the same drag prescription. An attempt is made to gauge the validity of the Chandrasekhar formula by comparison with restricted N-body solutions. It is found that the orbit-averaged evolution is consistent with such simulations.

Statler, Thomas S.↗

NASA NTRS1978

Computer experiments on the effect of retrograde stars in disk galaxies

Using large-scale N-body calculations for flat disk galaxies, we examine the effect of reversing the angular momentum for various fractions of the stars upon the global bar-forming mode. The initial conditions for these simulations are based on stationary states of two classes of models: the isochrones studied recently by Kalnajs by means of linear theory, and a model resembling the Schmidt model of our own Galaxy. In both cases, as the fraction of retrograde stars is increased, the growth of the bar-forming mode is inhibited (although not eliminated). These N-body results for the isochrones agree with the predictions of linear theory, quantitatively as well as qualitatively.

Zang, T. A.↗

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 NTRS1985

Navigation systems

The elements of the measurement and communications network comprising the global deep space navigation system (DSN) for NASA missions are described. Among the measurement systems discussed are: VLBI, two-way Doppler and range measurements, and optical measurements carried out on board the spacecraft. Processing of navigation measurement is carried out using two modules: an N-body numerical integration of the trajectory (and state transition partial derivatives) based on pre-guessed initial conditions; and partial derivatives of simulated observables corresponding to each actual observation. Calculations of velocity correction parameters is performed by precise modelling of all physical phenomena influencing the observational measurements, including: planetary motions; tracking station locations, gravity field structure, and transmission media effects. Some of the contributions to earth-relative orbit estimate errors for the Doppler/range system on board Voyager are discussed in detail. A line drawing of the DSN navigation system is provided.

Jordan, J. F.↗

NASA NTRS1990

The spiral-compact galaxy pair AM 2208-251: Computer simulations versus observations

The system AM2208-251 is a roughly edge-on spiral extending east-west with a smaller round compact E system about 60 arcsec east of the spiral nucleus along the major axis of the spiral. Bertola, Huchtmeier, and Zeilinger (1990) have presented optical spectroscopic as well as single dish 21 cm observations of this system. Their spectroscopic data show, via emission lines lambda lambda 3727-29A, a rising rotation curve near the nucleus. These spectroscopic observations may indicate a tidal interaction in the system. In order to learn more about such pairs, the authors simulated the interaction using the computer model developed by Miller (1976 a,b, 1978) and modified by the authors (Byrd 1986, 1987, 1988). To do the simulation they need an idea of the mutual orbits of the two galaxies. Their computer model is a two-dimensional polar N-body program. It consists of a self-gravitating disk of particles, within an inert axially symmetric stabilizing halo potential. The particles are distributed in a 24(radial) by 36(azimuthal) polar grid. Self consistent calculations can be done only within the grid area. The disk is modeled with a finite Mestel disk, where all the particles initially move in circular orbits with constant tangential velocities (Mestel 1963), resulting in a flat rotation curve. The gas particles in the spiral's disk, which make up 30 percent of its mass, collide in the following manner. The number of particles in each bin of the polar grid is counted every time step. If it is greater than a given critical density, all the particles in the bin collide, obtaining in the result the same velocities, equal to the average for the bin. This process produces clumps of gas particles-the star formation sites. The authors suppress the collision in the inner part of the disk (within the circle r = 6) to represent the hole seen in the gas in the nuclear bulge of spirals. They thus avoid spurious effects due to collisions in that region. They also varied the size of the collisional bins, which did not affect their conclusions.

Klaric, Mario↗

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 NTRS1992

Starbursts in the nuclei of shell galaxies

Simulations of shell-forming galaxy mergers are reported in which gas and stars are tracked separately to study the implications of observations that indicate fine structures in elliptical galaxies. The hybrid N-body/hydrodynamics code by Hernquist and Katz (1989) is employed neglecting the self-gravity of galaxies in tidal disruptions by massive primaries of small dwarf galaxies. The evolution of the galactic debris is examined showing that the gas and stars are effectively segregated with the stars forming sharp-edged features and the gas forming disks in the nucleus of the primary. Star formation is theorized to occur in the nuclei of shell galaxies due to shocks in the gas, and the added presence of a black hole can lead to an active galactic nucleus. The simulation is of interest to the interpretation of such observations as those of NGC1275 in which dormant elliptical galaxies are activated by mergers between large galaxies and gas-rich satellites.

Hernquist, Lars↗

NASA NTRS1971

Space Trajectories Error Analysis (STEAP) Programs. Volume 1: Analytic manual, update

Manual revisions are presented for the modified and expanded STEAP series. The STEAP 2 is composed of three independent but related programs: NOMAL for the generation of n-body nominal trajectories performing a number of deterministic guidance events; ERRAN for the linear error analysis and generalized covariance analysis along specific targeted trajectories; and SIMUL for testing the mathematical models used in the navigation and guidance process. The analytic manual provides general problem description, formulation, and solution and the detailed analysis of subroutines. The programmers' manual gives descriptions of the overall structure of the programs as well as the computational flow and analysis of the individual subroutines. The user's manual provides information on the input and output quantities of the programs. These are updates to N69-36472 and N69-36473.

Source record↗

NASA NTRS1975

Spiral structure produced by tidal interaction of galaxies

Simulation of a close passage of two galaxies where the self-gravity of the disk stars is taken into account was performed with the aid of a large-scale N-body model (Hohl and Hockney, 1969; Hohl, 1972). The calculation was performed by perturbing an initially stable axisymmetric disk consisting of 100,000 stars by the passage of a companion galaxy having one-fourth the mass of the primary galaxy. Parameters appropriate for M 51 and its companion NGC 5195 were used. The results show the evolution of a weak two arm spiral. The effect of large random velocities on the formation of spiral structure was studied by plotting the spatial distribution of stars for various velocity intervals. The spiral structure is quite pronounced for stars with the lower velocity dispersion, but it can hardly be detected for stars with higher velocities. The tidal interaction causes only slight changes in the azimuthally averaged density distribution.

Hohl, F.↗

NASA NTRS1971

Simulated trajectories error analysis program, version 2. Volume 2: Programmer's manual

A series of three computer programs for the mathematical analysis of navigation and guidance of lunar and interplanetary trajectories was developed. All three programs require the integration of n-body trajectories for both interplanetary and lunar missions. The virutal mass technique is used in all three programs. The user's manual contains the information necessary to operate the programs. The input and output quantities of the programs are described. Sample cases are given and discussed.

Vogt, E. D.↗