Engineering PapersSearch

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 19 records

On percolation as a cosmological test

Difficulties in the use of percolation as a complementary statistic for the galaxy clustering pattern are pointed out by studying simple toy models and dynamical N-body models that represent the competing clustering scenarios. The percolation properties are found not to be very sensitive to the presence of pancakes and strings once they are clumpy, and hence they do not distinguish properly between models that are very different. In the case of very smooth pancakes, the ability to percolate depends on sampling parameters, such as the mean number density and the volume, in a way which is unknown a priori because it depends on the same properties that the test ought to measure. This problem could, in principle, be eased by using volume-limited samples of high mean number density (an order of magnitude denser than the CfA redshift survey volume limited at 4000 km/s) and by comparing to models of identical number density and volume. An alternative approach, based on the sampling effects themselves, may provide a qualitative test for pancakes in samples of lower densities.

Dekel, A.

A unified N-body and statistical treatment of stellar dynamics. I - The hybrid code. II - Applications to globular cluster cores

A unified N-body and statistical treatment of stellar dynamics is developed and applied to the late stages of core collapse and early stages of post collapse evolution in globular clusters. A 'hybrid' computer code is joined to a direct N-body code which is used to calculate exactly the behavior of particles in the inner spatial region, and the combination is used to follow particles statistically in the outer spatial region. A transition zone allows the exchange of particles and energy between the two regions. The main application results include: formation of a hard central binary system, reversal of core collapse and expansion due to the heat input from this binary, ejection of the binary from the core, and recollapse of the core; density profiles that form a one-parameter sequence during the core oscillations; and indications that these oscillations will eventually cease.

Mcmillan, S. L. W.

A unified N-body and statistical treatment of stellar dynamics. III - Early postcollapse evolution of globular clusters

The period immediately following the core collapse phase in the evolution of a globular cluster is studied using a hybrid N-body/Fokker-Planck stellar dynamical code. Several core oscillations of the type predicted in earlier work are seen. The oscillations are driven by the formation, hardening, and ejection of binaries by three-body processes, and appear to decay on a timescale of about 10 to the 7th yr, for the choice of 'typical' cluster parameters made here. There is no evidence that they are gravothermal in nature. The mechanisms responsible for the decay are discussed in some detail. The distribution of hard binaries produced by the oscillations is compared with theoretical expectations and the longer term evolution of the system is considered.

Mcmillan, S. L. W.

A unified n-body and statistical treatment of stellar dynamics

The methods of a new 'hybrid' computed code for stellar dynamics are summarized. All particles in the inner spatial region are followed exactly via a direct N-body code and all particles in the outer spatial region are treated statistically via a distribution function and Fokker-Planck type methods. An intermediate region, with features of both, allows exchange of particles and energy between the outer and inner regions. The code is applied to the period just before core collapse and just after and the results are summarized.

Lightman, A. P.

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.

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.

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.

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.

A numerical experiment on the equilibrium and stability of a rotating galactic bar

A self-consistent, three-dimensional numerical experiment is performed on an N-body system whose initial state is a realization of a certain theoretical model of a rotating triaxial galaxy. The model is a stellar-dynamical counterpart of a uniformly rotating polytrope of index equal to 0.5. The aim of the experiment is to study the equilibrium of the system and, in particular, to test its stability. The experimental system behaves in the mean like a realization of the theoretical model for at least seven crossing times. The principal departure of the system from equilibrium is an oscillation which is identified as a radial pulsation. There is no indication in its behavior that the system is unstable with respect to anu mode with an e-folding time shorter than or of the order of two crossing times. Certain changes that occur in the state of the system are interpreted, with the aid of the theoretical model, as secular changes which result from a slight failure of our numerical methods to conserve the mass, energy, and angular momentum of the system; these effects are small enough that they do not vitiate the experiment on a dynamical time scale.

Miller, R. H.

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.

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.

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.