Search NASA⌕ Search

Engineering topics

Bodenheimer, Peter

Publications and source records attributed to Bodenheimer, Peter.

17 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Enhancement of the Accretion of Jupiters Core by a Voluminous Low-Mass Envelope

We present calculations of the early stages of the formation of Jupiter via core nucleated accretion and gas capture. The core begins as a seed body of about 350 kilometers in radius and orbits in a swarm of planetesimals whose initial radii range from 15 meters to 100 kilometers. We follow the evolution of the swarm by accounting for growth and fragmentation, viscous and gravitational stirring, and for drag-induced migration and velocity damping. Gas capture by the core substantially enhances the cross-section of the planet for accretion of small planetesimals. The dust opacity within the atmosphere surrounding the planetary core is computed self-consistently, accounting for coagulation and sedimentation of dust particles released in the envelope as passing planetesimals are ablated. The calculation is carried out at an orbital semi-major axis of 5.2 AU and an initial solids' surface density of 10/g/cm^2 at that distance. The results give a core mass of 7 Earth masses and an envelope mass of approximately 0.1 Earth mass after 500,000 years, at which point the envelope growth rate surpasses that of the core. The same calculation without the envelope gives a core mass of only 4 Earth masses.

Jupiter↗
NASA NTRS2007

Theoretical Developments in Understanding Massive Star Formation

Except under special circumstances massive stars in galactic disks will form through accretion. The gravitational collapse of a molecular cloud core will initially produce one or more low mass quasi-hydrostatic objects of a few Jupiter masses. Through subsequent accretion the masses of these cores grow as they simultaneously evolve toward hydrogen burning central densities and temperatures. We review the evolution of accreting (proto-)stars, including new results calculated with a publicly available stellar evolution code written by the authors.

stellar evolution↗
NASA NTRS1997

Effects of FU Orionis Outbursts on Protoplanetary Disks

In the early stages of work under this grant, we developed simulations to match the light curves of the three best studied systems: FU Ori, V1515 Cyg, & V1057 Cyg (Bell et al. 1995). We compared the details of model results to observations to test the validity of the thermal ionization instability model for outburst. In this paper, we were able to answer several of the key objections to the accretion disk outburst model for the FU Orionis phenomenon (eg. Herbig 1989). The declines in line width and reddening observed in V1057 Cyg following peak light had been used as arguments against the disk instability model. We showed these effects to be natural consequences of the slow outward progression and limited radial excursion of the ionization front during outburst. By the end of the grant period, we had begun combining the inner and outer disk models to derive the radiation expected in the planet forming region of the disk. A crucial step in this was the development of a radiative transfer model of the complete inner and outer disk system (Turner, Bodenheimer, & Bell 1997). In this work, wavelength dependent opacities were used to synthesize images and spectral energy distributions of FU Ori objects. New detailed opacities provided by Alexander (1995) allowed the resolution of coarse features such as silicate emission lines. Data for the fits were taken from the time dependent simulations in Bell et al. (1995) to which was added the effect of disk to disk or "self"-reprocessing which accounts for the illumination of the outer disk by the inner disk (Bell 1998). Through the course of work on this grant we have made considerable progress in computing detailed models of both the active outburst region of the disk and the outer quiescent disk. We have begun an investigation into the effects of a cocooning envelope. Direct comparisons with observations lend confidence that the basic underlying assumptions of the outburst model are reasonable. We are now ready to build upon these results to investigate the effects of outbursts on the central star and on the accumulation of small particles in the planet-forming regions.

Bodenheimer, Peter↗
NASA NTRS1995

The formation of protostellar disks. 2: Disks around intermediate-mass stars

Hydrodynamical calculations of the evolution of a collapsing, rotating axisymmetric 10 solar masses molecular clump, including the effects of radiative acceleration but without magnetic fields, are represented. The initial cloud is assumed to be uniformly rotating, centrally condensed sphere with rho is proportional to r(exp -2). Several cases are considered, in which both the overall clump size and the total amount of angular momentum are varied. The calculations show how a warm, quasi-hydrostatic disk surrounding a central unresolved core of only a few solar masses forms and grows in size and mass. The disk is encased in two distinct accretion shock fronts, both of which are several scale heights above the equatorial plane. At the end of the calculation of our standard case, the central unresolved region is found to have a mass of 2.7 solar masses and a ratio of rotational to gravitational energy of approximately 0.45, sufficiently large to be unstable to nonaxisymmetric perturbations. In addition, the inner portions of the disk containing most of the mass are unstable according to the local Toomre criterion, implying that also in this region nonaxisymmetric perturbations will lead to rapid evolution. Under the assumption that gravitational torques would transport angular momentum out of this region, a central core of less than or approximately 8 solar masses with a stable disk of greater than or approximately = 2 solar masses should result. Frequency-dependent radiative transfer calculations of the standard case at selected ages show how the continuum spectrum of the structure depends on the disk's orientation and age and how the observed isophotal contours vary with wavelength. Because of the strong dependence on viewing angle, continuum spectra alone should not be used to estimate the evolutionary stage of development of these objects. Comparable results were obtained for the other cases considered.

Yorke, Harold W.↗
NASA NTRS1995

Formation Of the Giant Planets By Concurrent Accretion Of Solids And Gas

New numerical simulations of the formation of the giant planets are presented, in which for the first time both the gas and planetesimal accretion rates are calculated in a self-consistent, interactive fashion. The simulations combine three elements: 1) three-body accretion cross-sections of solids onto an isolated planetary embryo, 2) a stellar evolution code for the planet's gaseous envelope, and 3) a planetesimal dissolution code within the envelope, used to evaluate the planet's effective capture radius and the energy deposition profile of accreted material. Major assumptions include: The planet is embedded in a disk of gas and small planetesimals with locally uniform initial surface mass density, and planetesimals are not allowed to migrate into or out of the planet's feeding zone. All simulations are characterized by three major phases. During the first phase, the planet's mass consists primarily of solid material. The planetesimal accretion rate, which dominates that of gas, rapidly increases owing to runaway accretion, then decreases as the planet's feeding zone is depleted. During the second phase, both solid and gas accretion rates are small and nearly independent of time. The third phase, marked by runaway gas accretion, starts when the solid and gas masses are about equal. It is engendered by a strong positive feedback on the gas accretion rates, driven by the rapid contraction of the gaseous envelope and the rapid expansion of the outer boundary, which depends on the planet's total mass. The overall evolutionary time scale is generally determined by the length of the second phase. The actual rates at which the giant planets accreted small planetesimals is probably intermediate between the constant rates assumed in most previous studies and the highly variable rates that we have used. Within the context, of the adopted model of planetesimal accretion, the joint constraints of the time scale for dissipation of the solar nebula and the current high-Z masses of the giant planets lead to estimates of the initial surface density (sigma(sub init)) of planetesimals in the outer region of the solar nebula. The results show sigma(sub init) approx. = 10 g/sq cm near Jupiter's orbit and that sigma(sub init) proportional to alpha(sup -2), where alpha is the distance from the Sun. These values are a factor of 3 - 4 times as high as that of the "minimum mass" solar nebula at Jupiter's distance and a factor of 2 - 3 times as high it Saturn's distance. Our estimates for the formation time of Jupiter and Saturn are 1 - 10 million years while those for Uranus fall in the range of 2 - 16 million years. These estimates follow from the properties of our Solar System and do not necessarily apply to giant planets in other planetary systems.

Pollack, James B.↗
NASA NTRS1994

Nonaxisymmetric evolution in protostellar disks

We present a two-dimensional, multigridded hydrodynamical simulation of the collapse of an axisymmetric, rotating, 1 solar mass protostellar cloud, which forms a resolved, hydrotastic disk. The code includes the effects of physical viscosity, radiative transfer and radiative acceleration but not magnetic fields. We examine how the disk is affected by the inclusion of turbulent viscosity by comparing a viscous simulation with an inviscid model evolved from the same initial conditions, and we derive a disk evolutionary timescale on the order of 300,000 years if alpha = 0.01. Effects arising from non-axisymmetric gravitational instabilities in the protostellar disk are followed with a three-dimensional SPH code, starting from the two-dimensional structure. We find that the disk is prone to a series of spiral instabilities with primary azimulthal mode number m = 1 and m = 2. The torques induced by these nonaxisymmetric structures elicit material transport of angular momentum and mass through the disk, readjusting the surface density profile toward more stable configurations. We present a series of analyses which characterize both the development and the likely source of the instabilities. We speculate that an evolving disk which maintains a minimum Toomre Q-value approximately 1.4 will have a total evolutionary span of several times 10(exp 5) years, comparable to, but somewhat shorter than the evolutionary timescale resulting from viscous turbulence alone. We compare the evolution resulting from nonaxisymmetric instabilities with solutions of a one-dimensional viscous diffusion equation applied to the initial surface density and temperature profile. We find that an effective alpha-value of 0.03 is a good fit to the results of the simulation. However, the effective alpha will depend on the minimum Q in the disk at the time the instability is activated. We argue that the major fraction of the transport characterized by the value of alpha is due to the action of gravitational torques, and does not arise from inherent viscosity within the smoothed particle hydrodynamics method.

Laughlin, Gregory↗
NASA NTRS1993

The formation of protostellar disks. I - 1 M(solar)

Hydrodynamical calculations of the collapse of an axisymmetric, rotating one solar mass protostellar cloud, including the effects of radiative transfer and radiative acceleration but without magnetic fields, are presented. The results include calculations of infrared protostellar spectra as a function of time and viewing angle. A numerical algorithm involving explicit nested grids is used to resolve the region of initial disk formation and at the same time to include the outer regions in the calculation. The central part of the protostar is modeled approximately. Initial conditions are systematically varied to investigate their influence on the evolution and final configuration of central star plus circumstellar disk. The initial state for the standard case is a centrally condensed molecular cloud core of one solar mass with a mean density of 8 x 10 exp -18 g/cu cm and a specific angular momentum at the outer edge of 7 x 10 exp 20 sq cm/s. The collapse is followed for 8 x 10 exp 4 yr, at which point 0.45 solar mass is contained in a rapidly rotating central object and most of the remainder in a surrounding equilibrium disk. The stability of this final structure is qualitatively analyzed.

Yorke, Harold W.↗
NASA NTRS1993

Luminosity functions for very low mass stars and brown dwarfs

A theoretical investigation of the luminosity function for low-mass objects to constrain the stellar initial mass function at the low-mass end is reported. The ways in which luminosity functions for low-mass stars are affected by star formation histories, brown dwarf and premain-sequence cooling rates and main-sequence mass luminosity relations, and the IMF are examined. Cooling rates and the mass-luminosity relation are determined through a new series of evolutionary calculations for very low mass stars and brown dwarfs in the range 0.05-0.50 solar mass. Model luminosity functions are constructed for specific comparison with the results of four recent observational surveys. The likelihood that the stellar mass function in the solar neighborhood is increasing at masses near the bottom of the main sequence and perhaps at lower masses is confirmed. In the most optimistic case, brown dwarfs contribute half of the local missing disk mass. The actual contribution is likely to be considerably less.

Laughlin, Gregory↗
NASA NTRS1993

Planets around pulsars - Implications for planetary formation

Data on planets around pulsars are summarized, and different models intended to explain the formation mechanism are described. Both theoretical and observational evidence suggest that very special circumstances are required for the formation of planetary systems around pulsars, namely, the prior presence of a millisecond pulsar with a close binary companion, probably a low mass main-sequence star. It is concluded that the discovery of two planets around PSR 1257+12 is important for better understanding the problems of dynamics and stellar evolution. The process of planetary formation should be learned through intensive studies of the properties of disks near young objects and application of techniques for detection of planets around main-sequence solar-type stars.

Bodenheimer, Peter↗
NASA NTRS1993

Stellar multiple systems - Constraints on the mechanism of origin

The discovery rate of premain-sequence binary and multiple systems has increased dramatically over the last few years, providing us with our first insights into the nature of such systems at young ages. The observational data derived from spectroscopic, occultation, speckle, visual, and other techniques are reviewed. Where possible, the premain-sequence binary population is compared with main-sequence systems. These results are interpreted in terms of binary formation mechanisms. While wide binaries (P is greater than 100 yr) are ordinarily explained by a fragmentation process during the collapse of a rotating protostar, the origin of some systems by capture is a possibility. The origin of close systems is not well understood: various processes, such as fragmentation during late stages of collapse, gravitational instabilities in disks, or orbital decay from a longer-period system are examined.

Bodenheimer, Peter↗
NASA NTRS1990

The formation phase of the solar nebula

Hydrodynamical calculations of the collapse of an axisymmetric, rotating protostellar cloud, with radiation transport and without magnetic fields, are presented. The collapse is assumed to start from a centrally condensed sphere of radius 5 x 10 to the 15th cm, a mean density of 4 x 10 to the -15th g/cu cm, a total mass of 1 solar mass, and a total angular momentum of 10 to the 53rd g sq cm per sec. The numerical grid is chosen to resolve the region of disk formation between 1 and 60 AU from the center. Frequency-dependent radiative transfer calculations show how the emergent spectrum of the structure depends upon viewing angle with respect to the rotation axis and how the observed isophotal contours should depend on wavelength and viewing angle. The central part of the protostar, interior to 1 AU, is not resolved numerically but is modeled approximately. At the end of the calculation, this region is found to have a mass of 0.6 solar mass and a ratio of rotational to gravitational energy of about 0.4, sufficiently large to be unstable to nonaxisymmetric perturbations. Although the disk is gravitationally stable according to the local Toomre criterion, the nonaxisymmetric structure in the center is likely to lead to angular momentum transport.

Bodenheimer, Peter↗
NASA NTRS1990

Brown dwarfs as close companions to white dwarfs

The influence of the radiation flux emitted by a white dwarf primary on the evolution of a closely orbiting brown dwarf (BD) companion is investigated. Full stellar evolutionary calculations are presented for both isolated and thermal bath cases, including effects of large variations in the atmospheric grain opacities. High grain opacities significantly increase the radii of the BDs, but the thermal bath does not. The major influence of the thermal bath is to increase substantially the surface temperature and luminosity of the BD at a given age. These results are compared with the observational properties of the possible BD companion of the white dwarf G29-38. Inclusion of both physical effects, high grain opacities and thermal bath, increases the mass range (0.034-0.063 solar masses) of viable models significantly, yet the final determination of whether the object is indeed a BD requires improvements in the observations of the system's properties.

Stringfellow, Guy S.↗
NASA NTRS1990

On the formation and expansion of H II regions

The evolution of H II regions in spherical clouds with small, constant-density cores and power-law density distributions r exp -w outside the core is described analytically. It is found that there is a critical exponent above which the cloud becomes completely ionized. Its value in the formation phase depends on the initial conditions, but it has a well-defined value w(crit) = 3/2 during the expansion phase. For w less than w(crit), the radius of the H II region grows at a given rate, while neutral mass accumulates in the interphase between the ionization and shock fronts. For w = w(crit), the fronts move together without mass accumulation. Cases with w greater than w(crit) lead to the champagne phase: once the cloud is fully ionized, the expansion becomes supersonic. For self-gravitating disks without magnetic fields, the main features include a new 'variable-size' stage. The initial shape of the H II region has a critical point beyond which the disk becomes completely ionized.

Franco, Jose↗
NASA NTRS1989

Structure and evolution of gaseous protoplanets

The evolutionary stages of a giant planet are reviewed, starting with the formation of the solar nebula and ending at the present time. The general mechanism for formation of a giant planet is assumed to involve the accumulation of solid material to form a core which grows to a critical mass, after which rapid accretion of the gaseous envelope occurs. The object then evolves at constant mass to its present state. Detailed numerical simulations of these evolutionary phases, under the assumption of spherical symmetry are reviewed, with particular attention to formation time scales, chemical composition, rotation, and satellite formation.

Bodenheimer, Peter↗
NASA NTRS1989

The impact of stellar evolution on planetary system development

The connection between stellar evolution and planet formation is investigated. Particular attention is given to the problem posed by the fact that the formation of Jupiter occurred before the formation of Mars and that the formation of the solid core of Saturn was completed before the dissipation of the gas in the nebula. Several possible solutions to this problem are suggested.

Bodenheimer, Peter↗
NASA NTRS1987

Numerical studies of collapsing interstellar clouds

Numerical simulation of the structure and evolution of interstellar clouds was conducted. Steps were taken toward an integrated treatment of the dynamical, thermal, and chemical processes entering model calculations, and a detailed study was made of radiative transfer in molecular lines to allow model predictions to be tested against empirical data. It is shown that the shapes of molecular lines are sensitive to details of the cloud structure and evolutionary state and are thus useful in inferring the cloud density, temperature, chemical composition, age, and initial conditions. The calculations have successfully reproduced and explained several observed cloud properties, including abundances of complex molecular species and the apparent depletion of CO in dense cores.

Bodenheimer, Peter↗