Primary crystallization and partial remelting of chondrules in the protoplanetary disk: Petrographic, mineralogical and chemical constraints recorded in zoned type-I chondrules
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Solar system origin model successfully predicts relative terrestrial abundances of noble gases by taking into account ambipolar diffusion
Explore the source record for details and available documents.
Models of giant protoplanets computed by DeCampli and Cameron (1979) show that Fe and other minerals in the planet interior are in a liquid state during one of the stages of protoplanet evolution. A model of coalescence of liquid drops was developed using the 'stochastic' collection equation of Slattery (1978); the growth times to droplets was much shorter than the period during which the drops are in a liquid state. Brownian collection quickly coalesced the tiny droplets to a radius of 0.005 cm; gravitational collection was required to form droplets of radii greater than 0.005 cm.
The 'core instability' mechanism presently postulated for the interaction of planetesimals with the envelopes of growing giant planets involves core growth through the accretion of solid bodies in the solar nebula; the core thereby becomes more effective in the gravitational attraction of gas from the surrounding nebula into an envelope, until a runaway accretion stage is reached. Two-body trajectory curves of planetesimals through this envelope are computed which demonstrate the depth of penetration and the cross-sections for capture of icy and rocky planetesimals.
A computational method for treating the generation of dynamo magnetic fields in astrophysical disks is presented. The numerical difficulty of handling the boundary condition at infinity in the cylindrical disk geometry is overcome by embedding the disk in a spherical computational space and matching the solutions to analytically tractable spherical functions in the surrounding space. The lowest lying dynamo normal modes for a 'thick' astrophysical disk are calculated. The generated modes found are all oscillatory and spatially localized. Tha potential implications of the results for the properties of dynamo magnetic fields in real astrophysical disks are discussed.
The following subject areas are covered: (1) magnetization of meteorites; (2) the possible origin of a nebular magnetic field; (3) the possibility of magnetic flares; (4) external manifestations; (5) dynamical effects of the magnetic fields; and (6) the problem of ionization.
Accretion disks evolve by transferring angular momentum outwards, which allows mass to flow inwards onto the central object. The radial structure and time-dependent behavior of low-mass, Keplerian, viscous disks is analyzed. Potential mechanisms for the angular momentum transport in disks are reviewed. An evolutionary model of a low mass disk around a pulsar is presented.
The spectral energy distributions (SED's) of all the known Class I (protostellar) sources in the Taurus molecular cloud were modelled. The Tereby, Shu, & Cassen (1984, TSC) density distribution for a rotating, infalling envelope was adopted. The radiative equilibrium temperature distribution from the spherical average of the TSC density distribution was calculated. The resulting spherically-symmetric temperature distribution then provides the source function to obtain the emergent spectrum at a given inclination angle i from the formal solution of the transfer equation, using the exact density (opacity) distribution. Results showed that the SED's of the protostar candidates in Taurus can be reproduced with TSC models having infall rates close to the values predicted by the theory of isothermal cloud collapse. Flat Spectrum T Tauri Stars were studied. The mid- to far-infrared fluxes of 'flat spectrum' T Tauri stars can be explained by radiative equilibrium emission from infalling dusty envelopes. This explanation indicated that models employing 'active' disks, in which the temperature distribution is a parameterized power law, should be invoked with caution. Infall also naturally explains the scattered light nebulae detected around many flat-spectrum sources. Results showed that the SED's of the protostar candidates in Taurus can be reproduced with TSC models having infall rates close to the values predicted by the theory of isothermal cloud collapse. Flattened infalling envelope models are also being studied.
Turbulence is known to possess structure on many scales. The largest or integral scale L is usually taken to be the largest dimension of the system. Two different turbulent regimes are likely to be of importance. Early stage convective turbulence in a hot nebula probably extends over the entire nebula scale height, and probably has typical eddy frequency comparable to the orbit frequency. At a later stage, midplane shear layer turbulence can be driven by the presence of a differently rotating, densely settled particle layer, with different length and timescales.
Many primitive meteorites contain high concentrations of chondrules - millimeter-sized beads of glassy silicate - in a matrix of essentially unprocessed material. These chondrules were apparently melted by some transient high-energy event, then cooled rapidly on timescales of minutes to hours. Among the phenomena that have been suggested as driving mechanisms for the heating event are: (1) high-velocity impacts; (2) magnetic flares; (3) aerodynamic drag heating; and (4) lightning. Although nebular lightning has been made of the efficiency of such lightning in melting chondrule precursors, a quantitative treatment of processes that might produce nebular lightning has not been attempted. We have developed a one-dimensional nebular lightning model that utilizes a spectrum of particle sizes. This model has also been applied to the Earth as a test of its predictive ability. Particles grow to an initial size distribution, then undergo elastic charge-separating collisions. Charge separation is opposed by the mobile free charge in the nebula, which can either accumulate on and neutralize the charged ice particles or move to cancel the large scale electric field. Since the free charge density depends on the grain size distribution, it is determined self-consistently as the particle distribution changes.
It is thought that planets form from solid particles in a flattened, rotating, 99% gaseous nebula. These grains gradually coagulate into millimeter-to-meter sized aggregates which settle toward the midplane of the nebula. It is widely believed that the resulting dense layer eventually becomes gravitationally unstable and collapses into 'planetesimals.' A new numerical model is presented to simulate the predominant processes (gravitation, vertical convection, and shear-driven turbulence) during the stage while the particulate material is still dispersed about the midplane of the nebula. In our previous work, particles were assumed to be spheres of a single radius; in the present work, particles are spheres of different radii. Results indicate that neither a broad nor a narrow distribution of particle sizes is likely to become gravitationally unstable.
In this paper, we study the stage of planetary formation during which the particulate material is still dispersed as centimeter-to-meter sized primordial aggregates. During this stage, particles are able to settle toward the midplane into a layer of mass density comparable to or much greater than that of the gas. Nonlinear, coupled interactions between the particles and the nebula gas become significant and ultimately determine the vertical profiles of the particle density and the mean velocities of the particles and the nearby gas. This is the environment in which the earliest planetesimals probably form. Our numerical models rely on the Reynolds averaged Navier-Stokes equations for the gas and particles and are fully viscous and turbulent (and, for the particles. compressible). Our turbulence modeling uses a Prandtl local shear parametrization, validated by laboratory experiments. We have developed a new model for particle diffusion involving the profile of the nebula gas turbulence and the particle Schmidt number, which is a function of particle size and density. We model a cool, quiescent nebula at 1 AU (280 K) and 10 AU (90 K) and a higher temperature stage (1000 K) at 1 AU. Our main results include: (a) rapid accretion of planetesimals by gravitationally unstable fragmentation on an orbital time scale (the 'Goldreich-Ward instability') is unlikely to occur until objects have already accreted by some other process to the mass of the largest known meteorite samples, if at all; (b) from 'seeds' as small as 10 m, growth of 10- to 100-km planetesimals can proceed rapidly by drift-augmented accretion in the particle-rich midplane with orbital decay of about 1% for the growing planetesimals; (c) outward transport of vapor and small entrained chips can account for significant radial compositional and mineralogical mixing in primitive meteorite parent bodies. Other implications may be drawn concerning particle random collisional velocities and the time scales of midplane turbulence.
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.
Significant results were obtained and published in the literature. The first optical detection of a circumbinary disk was reported in the ApJ at millimetric wavelengths. The size and inclination of this disk were found to be consistent with millimetric observations. Evidence was found for a cavity inside the disk as theory predicts from dust clearing by the stellar companion.
Our research focused on the structure and composition of the disks in the T Tau and Haro 6-10 multiple star systems and on the nature of the 'infrared companions' that these systems contain. This work has resulted in two papers, one on T Tau, presently under review at the Astrophysical Journal, and the other, on Haro 6-10, about to be submitted to Astronomy and Astrophysics. In the paper 'Evidence for Extinction and Accretion Variability in T Tau S' by Tracy L. Beck, L. Prato, and M. Simon, we present angularly resolved spectra of T Tau N (the visible star) and T Tau S (the infrared companion, IRC) in the three micrometer water ice-feature and the K-band. Most of the water ice absorption lies along the line of sight to T Tau S, confirming it is seen through stronger obscuration. A decrease in the ice-band absorption toward T Tau S between 12/98 and 1/00, significant, at the two-sigma level, was associated with an increase in its near IR flux. Bracket gamma emission is detected in T Tau N and S, and H2 (2.12 micrometer) emission only toward T Tau S, consistent with previous studies of IRCs. Our results suggest that the near IR variability of T Tau S is probably caused by both variations in accretion rate and extinction. Our paper on Haro 6-10, 'The Near IR and Ice Band Variability of Haro, 6-10' by Ch. Leinert, T.L. Beck, S. Ligori, M. Simon, J. Woitas, and R.R. Howell, represents a fusion of originally independent efforts at the Max Planck Institut fur Astronomie (Heidelberg) and Stony Brook. Our combined observations demonstrate that both Haro 6-10 S (the visible star) and Haro 6-10 N (the IRC) vary significantly in near IR flux on time scales as short as a month. The substantial decrease of Haro 6-10 S in the last four years carries the photometric signature of increased extinction. However, a comparable K-band flux increase of the IRC is associated with a dimming at H so cannot be explained by lower extinction. Absorption in the 3.1 micrometer ice band was always greater toward the IRC during our observations indicating a large amount of obscuring material. We detect variable extinction in the ice-band toward Haro 6-10 S and the IRC, significant at the 3.5-sigma and 2.0-sigma levels, respectively.