Search NASA⌕ Search

SEARCH · Search NASA

Results for “self-gravity”

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
All sourcesNASA NTRSDOE OSTI
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 NTRS1993

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

Numerical modeling of the interstellar medium in galactic disks

We have been developing detailed hydrodynamic models of the global interstellar medium in the hope of understanding the mass and volume occupied by various phases, as well as their structure and kinematics. In our model, the gas is modeled by one fluid while representative Pop 1 stars are modeled by a second fluid. The two fluids are coupled in that the gas forms into stars at a rate given by a Schmidt law while stellar mass loss returns matter into the gas phase (on a time scale of 100 Myr). Also, the stars heat the gas through stellar winds and the gas cools through optically thin radiation. The time behavior of these two fluids is studied in two spatial dimensions with the Eulerian finite difference numerical hydrodynamic code Zen. The two spatial dimensions are along the plane of a disk (x, total length of 2 kpc) and perpendicular to the disk (z, total height of +/- 15 kpc) and a galactic gravitational field in the z direction, typical of that at the solar circle, is imposed upon the simulation; self-gravity and rotation are absent. For the boundary conditions, outflow is permitted at the top and bottom of the grid (z = +/- 15 kpc) while periodic boundary conditions are imposed upon left and right sides of the grid. As initial conditions, we assumed a gaseous distribution like that seen for the H1 by earlier researchers, although the results are insensitive to the initial conditions. We have run simulations in which the heating due to stars, parameterized as a stellar wind velocity, a, is varied from low (a = 150 km/s), to intermediate (a = 300 km/s), to high (a = 600 km/s). Since the intermediate case is roughly equivalent to the Galactic energy injection rate from supernovae, this summary will concentrate on results from this simulation.

Rosen, A.↗
NASA NTRS1993

Spokes in ring galaxies

We examine the response of self-gravitating primary galaxies consisting of dark matter halos and disks containing both stars and gas to collisions with less massive companions. The primaries were constructed using a technique which makes it possible to realize multi-component systems that are stable and virtually in precise equilibrium. A total of 65,536 particles were employed to represent the primary and 4096 to represent the companion. Half the particles in the primary comprise its halo and the other half its disk. Gas makes up 10 percent of the disk mass and is represented by 8192 of the disk particles. A system of units is used where the gravitational constant, total disk mass, and disk exponential scale length are unity. The primary motivation of the present study is to determine whether effects associated with dissipation and self-gravity can account for the unusual morphology of the Cartwheel galaxy.

Hernquist, Lars↗
NASA NTRS1992

Numerical simulation of satellite-ring interactions - Resonances and satellite-ring torques

The Krook kinetic equation for planetary rings is numerically solved in two spatial dimensions and in time, with (1) interparticle collisions and (2) satellite-forcing, but (3) without self-gravity, for the case of a flattened planetary ring that undergoes gravitational perturbation by a nearby satellite. It is noted that the amplitude of wakes is limited by purely kinematic effects, even in the absence of collisions. Attention is given to the results of a simulation of an inner Lindblad-resonance location, as the distribution approaches steady state; these simulations do not show an increase in velocity dispersion in the resonance zone, obviating a net torque.

Brophy, Thomas G.↗
NASA NTRS1993

Nonlinear waves and solitons in molecular clouds

We begin a study of nonlinear wave phenomena in molecular clouds. These clouds exhibit highly nonlinear structure that is often described in terms of 'clumps' and 'filaments' which are bouncing around, twisting, and colliding within the cloud. These clouds are important because they ultimately produce the initial conditions for the star formation process. Our motivation is to explore the possibility that solitons (i.e., spatially localized, single-hump wave entities which often exhibit remarkable stability) can live in these molecular clouds and produce their observed structure. In this paper we focus on the case of one spatial dimension, and we show that a rich variety of nonlinear waves can exist in molecular cloud fluid systems (where self-gravity is included). We show that in the absence of magnetic fields no true soliton solutions are allowed, although highly nonlinear waves (whose crests become widely spaced and thus soliton-like) do exist. For clouds with embedded magnetic fields, we derive a model equation which describes the behavior of wave phenomena; this model equation allows solutions which correspond to nonlinear waves, solitons, and topological solitons. We briefly consider the stability of these wave entities and discuss the possible role they play in molecular cloud dynamics.

Adams, Fred C.↗
NASA NTRS1993

Gravitational tides in the outer planets. I - Implications of classical tidal theory. II - Interior calculations and estimation of the tidal dissipation factor

Classical tidal theory is applied to the atmospheres of the outer planets. The tidal geopotential due to satellites of the outer planets is discussed, and the solution of Laplace's tidal equation for Hough modes appropriate to tides on the outer planets is examined. The vertical structure of tidal modes is described, noting that only relatively high-order meridional mode numbers can propagate vertically with growing amplitude. Expected magnitudes for tides in the visible atmosphere of Jupiter are discussed. The classical theory is extended to planetary interiors taking the effects of spherically and self-gravity into account. The thermodynamic structure of Jupiter is described and the WKB theory of the vertical structure equation is presented. The regions for which inertial, gravity, and acoustic oscillations are possible are delineated. The case of a planet with a neutral interior is treated, discussing the various atmospheric boundary conditions and showing that the tidal response is small.

Ioannou, Petros J.↗
NASA NTRS1993

Ejecting basaltic achondrites from Vesta: Hydrodynamical impact models

Vesta is a large asteroid whose crust is mostly basaltic. Spectral heterogeneity suggests a sizable olivine feature which may be explained as impact excavation. The spectral data probably show a localized feature approximately 200 km in diameter or a diffuse feature approximately 400 km in diameter. Lightcurve irregularities suggest heterogeneity on a similar scale. This heterogeneity may represent the crater bowl, the extent of its ejecta deposit, or indeed something unrelated to cratering. In any case, drawing direct inferences about the state of Vesta's surface on the basis of these observations involves substantial speculation. We are still far from understanding the boundary between cratering and catastrophic disruption, particularly on targets for which strength and self-gravity both matter. But we are now able to model the underlying physical process -- dynamic fragmentation -- accurately with fragmentation hydrocodes such as SALE 2D and SPH3D. We present several impact scenarios for Vesta; our study is similar to a previous impact model for the formation of Stickney crater on Phobos. We illustrate the effects of impactors of different sizes and velocities, and the effect of gravity and structural inhomogeneity.

Asphaug, E.↗
NASA NTRS1994

A model for collisionally induced disturbed structure in disk galaxies

We derive analytic expressions, using the impulse and epicycle approximations, which describe the kinematic response of a disk galaxy following a collision with a second spherical galaxy which collides perpendicular to, but not through the center of, the disk. This model can reporduce the morphologies found in n-body experiments in which distant encounters produce two-armed spiral patterns and more central collisions produce rings in the disk galaxy, thereby confirming that simple kinematics can be used to describe the early evolution of these systems. Application of this procedure provides a convenient method with which to conduct parameter studies of these collisions. Comparison of the kinematic description with a fully self-gravitating, three-dimensional n-body/gasdynamics computer model shows that the disk galaxy's response is initially well represented by the kinematic model but that the self-gravity of the disk becomes important at longer times after the collision. The flows of gas and stars decouple from one another where stellar orbits cross, leaving regions of elevated gas density behind as the stars move freely past each other. If star formation rates are enhanced in these regions of high gas density, active star formation could be taking place where there is no corresponding dense feature in the old stellar population.

Gerber, Richard A.↗
NASA NTRS1994

Dynamics of binary-disk interaction. 1: Resonances and disk gap sizes

We investigate the gravitational interaction of a generally eccentric binary star system with circumbinary and circumstellar gaseous disks. The disks are assumed to be coplanar with the binary, geometrically thin, and primarily governed by gas pressure and (turbulent) viscosity but not self-gravity. Both ordinary and eccentric Lindblad resonances are primarily responsible for truncating the disks in binaries with arbitrary eccentricity and nonextreme mass ratio. Starting from a smooth disk configuration, after the gravitational field of the binary truncates the disk on the dynamical timescale, a quasi-equilibrium is achieved, in which the resonant and viscous torques balance each other and any changes in the structure of the disk (e.g., due to global viscous evolution) occur slowly, preserving the average size of the gap. We analytically compute the approximate sizes of disks (or disk gaps) as a function of binary mass ratio and eccentricity in this quasi-equilibrium. Comparing the gap sizes with results of direct simulations using the smoothed particle hydrodynamics (SPH), we obtain a good agreement. As a by-product of the computations, we verify that standard SPH codes can adequately represent the dynamics of disks with moderate viscosity, Reynolds number R approximately 10(exp 3). For typical viscous disk parameters, and with a denoting the binary semimajor axis, the inner edge location of a circumbinary disk varies from 1.8a to 2.6a with binary eccentricity increasing from 0 to 0.25. For eccentricities 0 less than e less than 0.75, the minimum separation between a component star and the circumbinary disk inner edge is greater than a. Our calculations are relevant, among others, to protobinary stars and the recently discovered T Tau pre-main-sequence binaries. We briefly examine the case of a pre-main-sequence spectroscopic binary GW Ori and conclude that circumbinary disk truncation to the size required by one proposed spectroscopic model cannot be due to Linblad resonances, even if the disk is nonviscous.

Artymowicz, Pawel↗
NASA NTRS1995

0n the dynamics of warped accretion disks

We present some analyses of warping modes in inviscid near Keplerian disks taking the three-dimensional structure fully into account. The results of this investigation verify the validity of a vertical averaging approximation for thin disks when the radial wavelength is significantly longer than the disk thickness. They also indicate that long wavelength disturbances may persist for long times. Shorter wavelength disturbances in non-self gravitating disks are found to propagate with little dispersion at a speed related to the sound speed. When self-gravity becomes important, fast and slow waves are found which also propagate with little dispersion. When a small viscosity is included, the evolution of the disturbances becomes more diffusive in character.

Papaloizou, J. C. B.↗
NASA NTRS1994

A quantitative analysis of IRAS maps of molecular clouds

We present an analysis of IRAS maps of five molecular clouds: Orion, Ophiuchus, Perseus, Taurus, and Lupus. For the classification and description of these astrophysical maps, we use a newly developed technique which considers all maps of a given type to be elements of a pseudometric space. For each physical characteristic of interest, this formal system assigns a distance function (a pseudometric) to the space of all maps: this procedure allows us to measure quantitatively the difference between any two maps and to order the space of all maps. We thus obtain a quantitative classification scheme for molecular clouds. In this present study we use the IRAS continuum maps at 100 and 60 micrometer(s) to produce column density (or optical depth) maps for the five molecular cloud regions given above. For this sample of clouds, we compute the 'output' functions which measure the distribution of density, the distribution of topological components, the self-gravity, and the filamentary nature of the clouds. The results of this work provide a quantitative description of the structure in these molecular cloud regions. We then order the clouds according to the overall environmental 'complexity' of these star-forming regions. Finally, we compare our results with the observed populations of young stellar objects in these clouds and discuss the possible environmental effects on the star-formation process. Our results are consistent with the recently stated conjecture that more massive stars tend to form in more 'complex' environments.

Wiseman, Jennifer J.↗
NASA NTRS1995

A new criterion for bar-forming instability in rapidly rotating gaseous and stellar systems. 1: Axisymmetric form

We analyze previous results on the stability of uniformly and differentialy rotating, self-gravitating, gaseous and stellar, axisymmetric systems to derive a new stability criterion for the appearance of torodial, m = 2 intermediate or I-modes and bar modes. In the process, we demonstrate that the bar modes in stellar systems and the m = 2 I-modes in gaseous systems have many common physical characteristics and only one substantial difference: because of the anisotropy of the stress tensor, dynamical instability sets in at lower rotation in stellar systems. This difference is reflected also in the new stability criterion. The new stability parameter alpha equals (T(sub J))/(absolute value of W) is formulated first for uniformly rotating systems and is based on the angular momentum content rather than on the energy content of a system. (T(sub J) is defined as ((L)(Omega(sub J)))/2; L is the total angular momentum; Omega(sub J) is the Jeans frequency introduced by self-gravity; and W is the total gravitational potential energy.) For stability of stellar systems alpha less than or equal to 0.254-0.258 while alpha less than or equal to 0.341-0.354 for stability of gaseous systems. For uniform rotation, one can write alpha = ((ft)/2)(exp 1/2), where t is defined as T/(absolute value of W), T is the total kinetic energy due to rotation, and f is a function characteristic of the topology/connectedness and the geometric shape of a system. Equivalently, alpha equals t/(chi), where chi is defined as Omega/Omega(sub J) and Omega is the rotation frequency. Using these forms, alpha can be extended to and calculated for a variety of differentially rotating, gaseous and stellar, axisymmetric disk and spheroidal models whose equilibrium structures and stability characteristics are known. In this paper, we also estimate alpha for gaseous torodial models and for stellar disk systems embedded in an inert or responsive 'halo.' We find that the new stability criterion holds equally well for all these previously published axisymmetric models.

Christodoulou, Dimitris M.↗
NASA NTRS1995

Star formation in the Galactic Center GMC cores: Sagittarius B2 and the dust ridge

The total far-infrared luminosity and the ionizing flux inferred from radio continuum observations of the Galactic center region imply a rate of star formation per unit mass of molecular material comparable to that in the Galactic disk. However, H2O and OH masers commonly found in sites of high-mass star formation are relatively rare in the nuclear disk. Far-infrared studies suggest that the formation rate of stars with masses greater than approximately 20 Solar Mass is reduced in the central region compared to the Galactic disk. Star formation might be suppressed currently in the central region as a result of the different geometry and strength of the magnetic fields there, which arguably might tend to inhibit cloud collapse. High gas pressures implied by observations of the diffuse X-ray emission suggest that giant molecular clouds (GMCs) in the nuclear disk may be held together by external pressure rather than self-gravity. The gravitational collapse leading to the formation of high density cores may thus be suppressed in all but the most massive clouds.

Lis, D. C.↗
NASA NTRS1993

A Simulation Atlas of Tidal Features in Galaxies

Detailed simulations of tidally induced structure in disk galaxies have either concentrated on specific systems or consisted of a few encounters with relatively small numbers of particles and no self-gravity. Observers need a 'dictionary' of simulations that covers many encounter parameters with fine morphological resolution and includes effects of self-gravitation. Observers can then search the dictionary for the parameters that best match a particular observed morphology. Alternatively, the dictionary can be used with observational samples for statistical studies of system parameters. To fill this need, we present a survey of model tidal encounters using a self-gravitating, 180,000 particle, two-component ('stars' and 'gas') disk. A wide variety of fascinating morphologies results. There are 86 different encounters that vary orbit tilt, perigalacticon distance, galaxy to companion mass ratio, and the amount of halo dark matter relative to the disk. For morphological comparisons, over 1700 images of the entire survey are available in video form. While there is a rich variety of tidal structure covering much of this parameter space, some general patterns may be remarked. There is a strong orbital inclination dependence of the symmetry of tidal patterns, most symmetric for planar orbits and nearly one-sided for polar encounters. Retrograde encounters produce only broad fanlike global patterns, but rich small-scale internal structure. In both kinds of encounter, our numerical resolution allows us to track internal spiral structure driven by the outer material arms, especially in the lighter halo simulations. We note also that polar encounters generate series of expanding, essentially non-rotating loops resembling shell structures in some respects.

Howard, Sethanne↗
NASA NTRS1996

Cometary Nuclei and Tidal Disruption: The Geologic Record of Crater Chains on Callisto and Ganymede

Prominent crater chains on Ganymede and Callisto are most likely the impact scars of comets tidally disrupted by Jupiter and are not secondary crater chains. We have examined the morphology of these chains in detail in order to place constraints on the properties of the comets that formed them and the disruption process. In these chains, intercrater spacing varies by no more than a factor of 2 and the craters within a given chain show almost no deviation from linearity (although the chains themselves are on gently curved small circles). All of these crater chains occur on or very near the Jupiter-facing hemisphere. For a given chain, the estimated masses of the fragments that formed each crater vary by no more than an order of magnitude. The mean fragment masses for all the chains vary by over four orders of magnitude (W. B. McKinnon and P. M. Schenk 1995, Geophys. Res. Lett. 13, 1829-1832), however. The mass of the parent comet for each crater chain is not correlated with the number of fragments produced during disruption but is correlated with the mean mass of the fragments produced in a given disruption event. Also, the larger fragments are located near the center of each chain. All of these characteristics are consistent with those predicted by disruption simulations based on the rubble pile cometary nucleus model (in which nuclei are composed on numerous small fragments weakly bound by self-gravity), and with those observed in Comet D/Shoemaker-Levy 9. Similar crater chains have not been found on the other icy satellites, but the impact record of disrupted comets on Callisto and Ganymede indicates that disruption events occur within the Jupiter system roughly once every 200 to 400 years.

Schenk, Paul M.↗
NASA NTRS1994

The Disruption and Demise of Periodic Comet Shoemaker-Levy 9

The impact of the fragmented comet Shoemaker-Levy 9 (SL9) into Jupiter this July promises to change our understanding of the outer solar system. More than twenty mountain-sized conglomerates of ice and rock will hit the atmosphere at approx. 50 km/s over the course of a week beginning July 16, releasing approx. 10(exp 4) to 10(exp8) megatons of energy per burst, and providing unique and perhaps pivotal clues to the properties of comets and the physics of massive atmospheres. Because the fragments will strike the far side of Jupiter, data acquisition, analysis and interpretation will be quite sensitive to the actual size and energy of the fragments. We therefore examine an event which took place two summers ago, unnoticed and unobserved: the disruption of SL9 into a "string of pearls' as it passed within the Roche limit at perijove. We first demonstrate, on the basis of timescales of tidal interaction, that the comet could not have broken into 20+ fragments through a hierarchy of brittle fracture events. Next, noting that the tidal stress was too weak to have even fragmented an uncompressed mass of freshly fallen snow, we run models for a strengthless comet held together only by self-gravity. We explore the initial size, density, and rotation. We conclude that a 4 km diameter comet (smaller if a prograde rotator) of density approx. 0.5 g/cu cm disrupts and disperses into a chain of fragments similar to Shoemaker-Levy 9, whether we begin with 21, 85, 169, 700 or 2000 sub-grains. Gravitational reaccumulation is evidently the answer, and there is no need to invoke the presence of 21 "cometesimals" as the subscale of the comet. To explain how a comet can be weaker than uncompacted snow, we show that the ring-plane crossing prior to perijove could have caused total damage. Finally, we compute the tidal stress on impactors as they approach Jupiter this July. Objects of various density are moderately distorted but not disrupted by the time they strike the planet.

Asphaug, Erik↗
NASA NTRS1998

Tidal Disruption of Strengthless Rubble Piles: A Dimensional Analysis

A relatively simple prescription for estimating the number of debris clumps (n) that form after a catastrophic tidal disruption event is presented. Following the breakup event, it is assumed that the individual debris particles follow keplerian orbits about the planet until the debris' gravitational contraction timescale (t(sub c)) becomes shorter than its orbital spreading timescale (t(sub s)). When the two timescales become comparable, self-gravity breaks up the debris train into n = L/D clumps, which is the debris length/diameter ratio at that instant. The clumps subsequently orbit the planet independent of each other. The predicted number of clumps n is in good agreement with more sophisticated N-body treatments of tidal breakup for parabolic encounters, and the dependence of n upon the progenitor's density as well as its orbit is also mapped out for hyperbolic encounters. These findings may be used to further constrain both the orbits and densities of the tidally disrupted bodies that struck Callisto and Ganymede. A cursory analysis shows that the Gomul and Gipul crater chains on Callisto, which have the greatest number of craters among the known chains, were formed by projectiles having comet-like densities estimated at rho(sub o) < 1 gm/cc.

Hahn, Joseph M.↗