Engineering Papers⌕ Search

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.

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

At least 73 records · Page 4

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

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

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

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

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

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

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

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↗

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

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↗

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↗

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

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

Gas distribution and starbursts in shell galaxies

Detailed maps of most elliptical galaxies reveal that, whereas the greatest part of their luminous mass originates from a smooth distribution with a surface brightness approximated by a de Vaucouleurs law, a small percentage of their light is contributed by low surface brightness distortions termed 'fine structures'. The sharp-edged features called 'shells' are successfully reproduced by merger and infall models involving accretion from less massive companions. In this context, dwarf spheroidal and compact disk galaxies are likely progenitors of these stellar phenomena. However, it is probable that the sources of shell-forming material also contain significant amounts of gas. This component may play an important role in constraining the formation and evolution of shell galaxies. To investigate the effects of the gaseous component, numerical simulations were performed to study the tidal disruption of dwarf galaxies containing both gas and stars by more massive primaries, and the evolution of the ensuing debris. The calculations were performed with a hybrid N-body/hydrodynamics code. Collisionless matter is evolved using a conventional N-body technique and gas is treated using smoothed particle hydrodynamics in which self-gravitating fluid elements are represented as particles evolving according to Lagrangian hydrodynamic equations. An isothermal equation of state is employed so the gas remains at a temperature 104 K. Owing to the large mass ratio between the primary and companion, the primary is modeled as a rigid potential and the self-gravity of both galaxies is neglected.

Weil, Melinda L.↗

The morphology of and locations of star formation in impact induced ring galaxies

Observed ring galaxies appear to fall into two major types. The first tends to consist of isolated galaxies which display a smooth, apparently circular ring and a central nucleus. These have been variously classified as R(S) by de Vaucouleurs (1959) and as type O by Few and Madore (1986). The second class of ring galaxy nearly always has a close companion of comparable size (no less than about one tenth that of the ring galaxy). In these objects the ring is knotty in appearance, is usually elliptical, even when deprojected on the sky, and is often open on one side, having a 'horse shoe' or 'banana' shape. The nucleus does not usually appear at the center of the ring and is sometimes apparently absent, giving rise to an 'empty ring' galaxy. deVaucouleurs et al. (1976) designated this second type as RING, while Few and Madore (1986) have classified similar galaxies as P type. These galaxies have elevated far IR emission, bright HII regions, and blue spectral colors. The different environments of the two types or ring galaxy, together with their overall morphological and spectral differences suggest that the R(S)/O type are most probably the result of an instability that occurs in isolated galaxies, whereas the RING/P type appears to be the result of a recent collision between two roughly equal mass objects, at least one of which is a disk galaxy. Theys and Spiegel (1976) studied a sample of this latter type and identified three subclasses: RE: galaxies with crisp, empty rings; RN: galaxies like those of RE but with off-center nuclei; RK: galaxies having single dominant knots or condensations in the rings. A presentation of a preliminary understanding of the connections between these different observed forms in terms of parameters which are intrinsic to the galaxy system, such as time since collision and impact parameter, and in terms of our line of sight view is the purpose of this paper. Here we report results we have obtained from three dimensional computer simulations of collisions between equal mass galaxies, one of which is a rotating, disk galaxy containing both gas and stars and the other is an elliptical containing stars only. We have used a combined n-body/SPH program (see Balsara, 1990) to model fully self consistent models in which the halo mass is 2.5 times that of the disk and gas comprises ten percent of the disk mass.

Lamb, Susan A.↗

Cloud-particle galactic gas dynamics and star formation

Galactic gas dynamics, spiral structure, and star formation are discussed in relation to N-body computational studies based on a cloud-particle model of the interstellar medium. On the small scale, the interstellar medium is seen as cloud-dominated and supernova-perturbed. It is noted that the cloud-particle model simulates cloud-cloud collisions, the formation of stellar associations, and supernova explosions as dominant local processes. On the large scale, in response to a spiral galactic gravitational field, global density waves and galactic shocks develop having large-scale characteristics similar to those found in continuum gas dynamical studies. Both the system of gas clouds and the system of young stellar associations forming from the clouds figure in the global spiral structure. However, with the attributes of neither assuming a continuum of gas (as in continuum gas dynamical studies) or requiring a prescribed equation of state (such as the isothermal condition), the cloud-particle picture retains much of the detail lost in earlier work. By detail is meant the small-scale features and structures so important in understanding the local, turbulent state of the interstellar medium as well as the degree of raggedness often seen to be superposed on the global spiral structure.

Roberts, W. W., Jr.↗

Collisional dynamics of perturbed particle disks in the solar system

Theoretical and computational studies were carried out on galactic and planetary disks. With the goal of addressing important open questions centered on galactic structure, the cloudy interstellar medium, giant molecular clouds, and star formation in galactic disks and the collisional dynamics of perturbed particulate matter in planetary disks, focus was largely on the self-gravitational effects, dissipative effects, and collisional dynamics of cloud-particle disks. N-body, 'cloud-particle' computational algorithms were developed for the purpose of isolating the role of gaseous self gravity from the roles of other dominant physical mechanisms and dynamical processes, e.g. the collisional dynamics and dissipative processes. The efforts focused largely on galactic disks show that observational constraints provide stringent tests of the numerical simulation techniques developed. Self gravitational effects of the galactic interstellar medium's gas clouds are included by means of Fourier Transform technique.

Roberts, William W., Jr.↗

Relaxation time in disk galaxy simulations.

An initially stationary and stable axisymmetric disk of stars with a mass spectrum of stars is used to determine the collisional relaxation time. The relaxation time as determined from the rate of energy equipartition was found to be 560 rotation periods for stars with 10 times the mean star mass and was 1700 rotation periods for stars with 0.55 times the mean star mass. These times are in general agreement with theoretical predictions for the relaxation time of the two mass groups. The results show that the model used for the large-scale gravitational N-body calculations is indeed 'collisionless.'

Hohl, F.↗