Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrostaticity”

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 1,081 records · Page 60

Protostellar collapse in a self-gravitating sheet

We present preliminary calculations of protostellar cloud collapse starting from an isothermal, self-gravitating gaseous layer in hydrostatic equilibrium. This gravitationally unstable layer collapses into a flattened or toroidal density distribution, even in the absence of rotation or magnetic fields. We suggest that the flat infalling envelope recently observed in HL Tau by Hayashi et al.is the result of collapse from an initially nonspherical layer. We also speculate that the later evolution of such a flattened, collapsing envelope can produce a structure similar to the 'flared disk' invoked by Kenyon and Hartmann to explain the infrared excesses of many T Tauri stars.

Hartmann, Lee↗

The vertical structure and stability of accretion disks surrounding black holes and neutron stars

The structure and stability of the inner regions of accretion disks surrounding neutron stars and black holes have been investigated. Within the framework of the alpha viscosity prescription for optically thick disks, we assume the viscous stress scales with gas pressure only, and the alpha parameter, which is less than or equal to unity, is formulated as alpha(sub 0)(h/r)(exp n), where h is the local scale height and n and alpha(sub 0) are constants. We neglect advective energy transport associated with radial motions and construct the vertical structure of the disks by assuming a Keplerian rotation law and local hydrostatic and thermal equilibrium. The vertical structures have been calculated with and without convective energy transport, and it has been demonstrated that convection is important especially for mass accretion rates, M-dot, greater than about 0.1 times the Eddington value, M-dot(sub Edd). Although the efficiency of convection is not high, convection significantly modifies the vertical structure of the disk (as compared with a purely radiative model) and leads to lower temperatures at a given M-dot. The results show that the disk can be locally unstable and that for n greater than or = 0.75, an S-shaped relation can exist between M-dot and the column density, sigma, at a given radius. While the lower stable branch (derivative of M-dot/derivative of sigma greater than 0) and middle unstable branch (derivative of M-dot/derivative of sigma less than 0) represent structures for which the gas and radiation pressure dominate respectively, the stable upper branch (derivative of M-dot/derivative of sigma greater than 0) is a consequence of the saturation of alpha. This saturation of alpha can occur for large alpha(sub 0) and at M-dot less than or = M-dot(sub Edd). The instability is found to occur at higher mass accretion rates for neutron stars than for black holes. In particular, the disk is locally unstable for M-dot greater than or = 0.5 M-dot(sub Edd) for neutron stars and for M-dot greater than or = M-dot(sub Edd) for black holes for a viscosity prescription characterized by n = 1 and alpha(sub 0) = 10.

Milsom, J. A.↗

Observing stellar coronae with the Goddard High Resolution Spectrograph. 1: The dMe star AU microscopoii

We report on an observation of AU Mic taken with the Goddard High Resolution Spectrograph (GHRS) aboard the Hubble Space Telescope. The data consist of a rapid sequence of spectra covering the wavelength range 1345-1375 A with a spectral resolution of 10,000. The observations were originally intended to search for spectral variations during flares. No flares were detected during the 3.5 hr of monitoring. A method of reducing the noise while combining the individual spectra in the time series is described which resulted in the elimination of half of the noise while rejecting only a small fraction of the stellar signal. The resultant spectrum was of sufficient quality to allow the detection of emission lines with an integrated flux of 10(exp -15) ergs/sq cm(sec) or greater. Lines of C I, O I, O V, Cl I, and Fe XXI were detected. This is the first indisputable detection of the 1354 A Fe XXI line, formed at T approximately = 10(exp 7) K, on a star other than the Sun. The line was well resolved and displayed no significant bulk motions or profile asymmetry. From the upper limit on the observed line width, we derive an upper limit of 38 km/s for the turbulent velocity in the 10(exp 7) K plasma. An upper limit is derived for the flux of the 1349 A Fe XII line, formed at T approximately = 1.3 x 10(exp 6) K. These data are combined with contemporaneous GHRS and International Ultraviolet Explorer (IUE) data to derive the volume emission measure distribution of AU Mic over the temperature range 10(exp 4)-10(exp 7) K. Models of coronal loops in hydrostatic equilibrium are consistent with the observed volume emission measures of the coronal lines. The fraction of the stellar surface covered by the footprints of the loops depends upon the loop length and is less than 14% for lengths smaller than the stellar radius. From the upper limit to the estimated width of the Fe XXI line prfile we find that the we cannot rule out Alfven wave dissipation as a possible contributor to the required quiescent loop heating rate.

Maran, S. P.↗

Magnetic shaping of planetary nebulae and other stellar wind bubbles

As in the case of the solar wind, the magnetic field in the wind from a magnetized, rotating star becomes increasingly toroidal with distance from the star. The strength of the magnetic field can be characterized by sigma, the ratio of toroidal magnetic energy density to kinetic energy density in the equatorial plane of the wind. A fast wind shocks against the external medium and creates a bubble whose volume is dominated by shocked gas. The toroidal magnetic field increases in the shocked bubble and can dominate the thermal pressure. Because of the low velocities in the bubble, hydrostatic equilibrium is a good approximation and allows the calculation of the thermal and magnetic pressure in the bubble, as in the model of Begelman and Li (1992) for the Crab Nebula. The structure, which is axisymmetric and extended in the polar direction, depends on two parameters: sigma nu(sub w)/w(sub 0), where nu(sub w) is the wind velocity and w(sub 0) is the shell velocity in the polar direction, and lambda = nu(sub a)/w(sub 0), where nu(sub a) is the velocity of the slow wind. For small values of lambda, there is a cusp in the shell in the equatorial plane, i.e., there is an equatorial ring. For larger values of lambda, the maximum of the surface density moves away from the equator i.e., a double ring structure. Our models should apply to planetary nebulae, if their central stars are sufficiently magnetized; the calculated shapes do resemble the observed shapes of planetaries. In all cases, our model predicts that X-ray emission from the bubble is concentrated toward the polar axis. Finally, we briefly discuss the asymmetry of the Crab Nebula and 3C 58.

Chevalier, Roger A.↗

Inside the supernova: A powerful convective engine

We present an extensive study of the inception of supernova explosions by following the evolution of the cores of two massive stars (15 and 25 Solar mass) in multidimension. Our calculations begin at the onset of core collapse and stop several hundred milliseconds after the bounce, at which time successful explosions of the appropriate magnitude have been obtained. Similar to the classical delayed explosion mechanism of Wilson, the explosion is powered by the heating of the envelope due to neutrinos emitted by the protoneutron star as it radiates the gravitational energy liberated by the collapse. However, as was shown by Herant, Benz, & Colgate, this heating generates strong convection outside the neutrinosphere, which we demonstrate to be critical to the explosion. By breaking a purely stratified hydrostatic equilibrium, convection moves the nascent supernova away from a delicate radiative equilibrium between neutrino emission and absorption, Thus, unlike what has been observed in one-dimensional calculations, explosions are rendered quite insensitive to the details of the physical input parameters such as neutrino cross sections or nuclear equation of state parameters. As a confirmation, our comparative one-dimensional calculations with identical microphysics, but in which convection cannot occur, lead to dramatic failures. Guided by our numerical results, we have developed a paradigm for the supernova explosion mechanism. We view a supernova as an open cycle thermodynamic engine in which a reservoir of low-entropy matter (the envelope) is thermally coupled and physically connected to a hot bath (the protoneutron star) by a neutrino flux, and by hydrodynamic instabilities. This paradigm does not invoke new or modified physics over previous treatments, but relies on compellingly straightforward thermodynamic arguments. It provides a robust and self-regulated explosion mechanism to power supernovae that is effective under a wide range of physical parameters.

Herant, Marc↗

Loop models of low coronal structures observed by the Normal Incidence X-Ray Telescope (NIXT)

The X-ray pictures obtained with the Normal Incidence X-Ray Telescope (NIXT), apart from the ubiquitous coronal loops well known from previous X-ray observations, show a new and peculiar morphology: in many active regions there are wide and apparently low-lying areas of intense emission which resemble H alpha plages. By means of hydrostatic models of coronal arches, we analyze the distribution of temperature, density, emission measure, and plasma emissivity in the spectral band to which NIXT is sensitive, and we show that the above morphology can be explained by the characteristics of high pressure loops having a thin region of high surface brightness at the base. We therefore propose that this finding might help to identify high-pressure X-ray emitting coronal regions in NIXT images, and it is in principle applicable to any imaging instrument which has high sensitivity to 10(exp 4) - 10(exp 6) K plasma within a narrow coronal-temperature passband. As a more general result of this study, we propose that the comparison of NIXT observations with models of stationary loops might provide a new diagnostic: the determination of the loop plasma pressure from measurements of brightness distribution along the loop.

Peres, G.↗

Dynamics of emerging active region flux loops

The buoyant rise of a magnetic flux loop arising from a single perturbed segment of a toroidal flux ring lying slightly beneath the base of the convection zone is studied by way of numerical simulations. We have considered flux loop evolution assuming both solid-body rotation, and differential rotation consistent with recent results from helioseismology. Our major results are presented, and we offer some speculations on the decay of active regions, based on the results of our studies. We speculate that as plasma in the tube attempts to establish hydrostatic equilibrium along the field lines after the flux emergence has taken place, the tube field strength at some intermediate depths below the surface becomes sufficiently small at the surface portions of the tube (which have cooled and undergone convective collapse) become dynamically disconnected from those portions near the base of the convection zone. The surface proportions of the emerged flux tubes are then transported by motions near the photosphere, such as supergranular convection and meridional flow.

Fan, Y.↗

The r-process and neutrino-heated supernova ejecta

As a neutron star is formed by the collapse of the iron core of a massive star, its Kelvin-Helmholtz evolution is characterized by the release of gravitational binding energy as neutrinos. The interaction of these neutrinos with heated material above the neutron star generates a hot bubble in an atmosphere that is nearly in hydrostatic equilibrium and heated, after approximately 10 s, to an entropy of S/N(sub AS)k greater than or approximately = 400. The neutron-to-proton ratio for material moving outward through this bubble is set by the balance between neutrino and antineutrino capture on nucleons. Because the electron antineutrino spectrum at this time is hotter than the electron neutrino spectrum, the bubble is neutron-rich (0.38 less than or approximately = Y(sub e) less than or approximately = 0.47). Previous work using a schematic model has shown that these conditions are well suited to the production of heavy elements by the r-process. In this paper we have advanced the numerical modeling of a 20 solar mass 'delayed' supernova explosion to the point that we can follow the detailed evolution of material moving through the bubble at the late times appropiate to r-process nucleosynthesis. The supernova model predicts a final kinetic energy for the ejecta of 1.5 x 10(exp 51) ergs and leaves behind a remnant with a baryon mass of 1.50 solar mass (and a gravitational mass of 1.445 solar mass). We follow the thermodynamic and compositional evolution of 40 trajectories in rho(t), T(t), Y(sub e)(t) for a logarithmic grid of mass elements for the last approximately = 0.03 solar mass to be ejected by the proto-neutron star down to the last less than 10(exp -6) solar mass of material expelled at up to approximately = 18 s after core collapse. We find that an excellent fit to the solar r-process abundance distribution is obtained with no adjustable parameters in the nucleosynthesis calculations. Moreover, the abundances are produced in the quantities required to account for the present Galactic abundances. However, at earlier times, this one-dimensional model ejects too much material with entropies S/N(sub A)k approximately 50 and Y(sub e) approximately 0.46. This leads to an acceptable over production of N = 50 nuclei, particularly Sr-88, Y-89, and Zr-90, relative to their solar abundances. We speculate on various means to avoid the early overproduction and/or ejection of N = 50 isotonic nuclei while still producing and ejecting the correct amount of r-process material.

Woosley, S. E.↗

Super-Eddington radiation transfer in soft gamma repeaters

Bursts from soft gamma repeaters (SGRs) have been shown to be super-Eddington by a factor of 1000 and have been persuasively associated with compact objects. Super-Eddington radiation transfer on the surface of a strongly magnetic (greater than or equal to 10(exp 13) G) neutron star is studied and related to the observational constraints on SGRs. In strong magnetic fields, Thompson scattering is suppressed in one polarization state, so super-Eddington fluxes can be radiated while the plasma remains in hydrostatic equilibrium. We discuss a model which offers a somewhat natural explanation for the observation that the energy spectra of bursts with varying intensity are similar. The radiation produced is found to be linearly polarized to one part in 1000 in a direction determined by the local magnetic field, and intensity variations between bursts are understood as a change in the radiating area on the source. The net polarization is inversely correlated with burst intensity. Further, it is shown that for radiation transfer calculations in limit of superstrong magnetic fields, it is sufficient to solve the radiation transfer for the low opacity state rather than the coupled equations for both. With this approximation, standard stellar atmosphere techniques are utilized to calculate the model energy spectrum.

Ulmer, Andrew↗

Disruption of a helmet streamer by photospheric shear

Helmet streamers on the Sun have been observed to be the site of coronal mass ejections, dynamic events that eject coronal plasma and magnetic fields into the solar wind. We develop a two-dimensional (azimuthally symmetric) helmet streamer configuration by computing solutions of the time-dependent magnetohydrodynamic (MHD) equations, and we investigate the evolution of the configuration when photospheric shearing motions are imposed. We find that the configuration disrupts when a critical shear is exceeded, ejecting a plasmoid into the solar wind. The results are similar to the case of a sheared dipole magnetic field in a hydrostatic atmosphere (Mikic & Linker 1994). However, the presence of the outflowing solar wind makes the disruption significantly more energetic when a helmet streamer is sheared. Our resutls suggest that shearing of helmet streamers may initiate coronal mass ejections.

Linker, Jon A.↗

HST eclipse mapping of dwarf nova OY Carinae in quiescence: An 'Fe II curtain' with Mach approx. = 6 velocity dispersion veils the white dwarf

Hubble Space Telescope (HST) observations of the eclipsing dwarf nova OY Car in its quiescent state are used to isolate the ultraviolet spectrum (1150-2500 A at 9.2 A Full Width at Half Maximum (FWHM) resolution) of the white dwarf, the accretion disk, and the bright spot. The white dwarf spectrum has a Stark-broadened photospheric L(alpha) absorption, but is veiled by a forest of blended Fe II features that we attribute to absorption by intervening disk material. A fit gives T(sub w) approx. = 16.5 x 10(exp 3) K for the white dwarf with a solar-abundance, log g = 8 model atmosphere, and T approx. = 10(exp 4) K, n(sub e) approx. = 10(exp 13)/cu cm, N(sub H) approx. = 10(exp 22) sq cm, and velocity dispersion delta V approx. = 60 km/s for the veil of homogeneous solar-abundance local thermodynamic equilibrium (LTE) gas. The veil parameters probably measure characteristic physical conditions in the quiescent accretion disk or its chromosphere. The large velocity dispersion is essential for a good fit; it lowers (chi square)/778 from 22 to 4. Keplerian shear can produce the velocity dispersion if the veiling gas is located at R approx. = 5 R(sub W) with (delta R)/R approx. = 0.3, but this model leaves an unobscured view to the upper hemisphere of the white dwarf, incompatible with absorptions that are up to 80% deep. The veiling gas may be in the upper atmosphere of the disk near its outer rim, but we then require supersonic (Mach approx. = 6) but sub-Keplerian (delta V/V(sub Kep) approx. = 0.07) velocity disturbances in this region to produce both the observed radial velocity dispersion and vertical motions sufficient to elevate the gas to z/R = cos i = 0.12. Such motions might be driven by the gas stream, since it may take several Kepler periods to reestablish the disk's vertical hydrostatic equilibrium. The temperature and column density of the gas we see as Fe II absorption in the ultraviolet are similar to what is required to produce the strong Balmer jump and line emissions seen in optical spectra of OY Car and similar quiescent dwarf novae. The outer accretion disk is detected at mid-eclipse with a spectrum that rises from 0.05 to 0.3 mJy between 2000 and 2500 A, consistent with combinations of cool blackbodies, blended Fe II emission lines, and Balmer continuum emission. The total disk flux density is 0.5 mJy at 2500 A, and this shallow disk eclipse implies a roughly flat surface brightness distribution. The bright spot, somewhat bluer than the disk, has a flux density rising from 0.05 to 0.15 mJy between 1600 and 2500 A. The C IV emission line has a broad shallow eclipse, but the radial velocity variations observed during the eclipse do not clearly distinguish between a disk or wind origin. The only possible indications of boundary layer emission are fast UV flares that appear to arise from near the central object -- not from the bright spot.

Horne, Keith↗

X-ray-emitting gas surrounding the spiral galaxy NGC 891

We observed the edge-on spiral galaxy NGC 891 with the Position Sensitive Proportional Counter (PSPC) on Roentgen Satellite (ROSAT) to search for how extraplanar gas expected in the galactic fountain model. Diffuse X-ray emission surrounds the disk with a Half Width at Half Maximum (HWHM) for the surface brightness perpendicular to the disk of 50 sec (2.4 kpc) and a radial extent of approximately 6.5 kpc, both of which are similar in extent to the extended H(alpha) and radio halo component; the implied density scale height for the hot gas is 7 kpc. The spectrum is best fitted with a hard stellar component and a soft diffuse gas component of temperature 3.6 x 10(exp 6) K. The density of this gas is 2 x 10(exp -3)/cu cm, the luminosity is 4.4 x 10(exp 39) ergs/s, the mass is 1 x 10(exp 8) solar mass, and the pressure (P/k) is 1.4 10(exp 4) K/cu cm. These data are consistent with this gas participating in a galactic fountain, where the material approaches hydrostatic equilibrium before cooling at a rate of 0.12 solar mass/yr. The cooled material may be responsible for some of the H(alpha) emission.

Bregman, Joel N.↗

Diabatic forcing and intialization with assimilation of cloud water and rainwater in a forecast model

In this study, diabatic forcing, and liquid water assimilation techniques are tested in a semi-implicit hydrostatic regional forecast model containing explicit representations of grid-scale cloud water and rainwater. Diabatic forcing, in conjunction with diabatic contributions in the initialization, is found to help the forecast retain the diabatic signal found in the liquid water or heating rate data, consequently reducing the spinup time associated with grid-scale precipitation processes. Both observational Special Sensor Microwave/Imager (SSM/I) and model-generated data are used. A physical retrieval method incorporating SSM/I radiance data is utilized to estimate the 3D distribution of precipitating storms. In the retrieval method the relationship between precipitation distributions and upwelling microwave radiances is parameterized, based upon cloud ensemble-radiative model simulations. Regression formulae relating vertically integrated liquid and ice-phase precipitation amounts to latent heating rates are also derived from the cloud ensemble simulations. Thus, retrieved SSM/I precipitation structures can be used in conjunction with the regression-formulas to infer the 3D distribution of latent heating rates. These heating rates are used directly in the forecast model to help initiate Tropical Storm Emily (21 September 1987). The 14-h forecast of Emily's development yields atmospheric precipitation water contents that compare favorably with coincident SSM/I estimates.

Raymond, William H.↗

A simulation of the intracluster medium with feedback from cluster galaxies

We detail method and report first results from a three-dimensional hydrodynamical and N-body simulation of the formation and evolution of a Coma-sized cluster of galaxies, with the intent of studying the history of the hot, X-ray emitting intracluster medium. Cluster gas, galaxies, and dark matter are included in the model. The galaxies and dark matter fell gravitational forces; the cluster gas also undergoes hydrodynamical effects such as shock heating and PdV work. For the first time in three dimensions, we include modeling of ejection of processed gas from the simulated galaxies by winds, including heating and heavy element enrichment. For comparison, we employ a `pure infall' simulation using the same initial conditions but with no galaxies or winds. We employ an extreme ejection history for galactic feedback in order to define the boundary of likely models. As expected, feedback raises the entropy of the intracluster gas, preventing it from collapsing to densities as high as those attained in the infall model. The effect is more pronounced in subclusters formed at high redshift. The cluster with feedback is always less X-ray luminous, but experiences more rapid luminosity evolution, than the pure infall cluster. Even employing an extreme ejection model, the final gas temperature is only approximately 15% larger than in the infall model. The radial temperature profile is very nearly isothermal within 1.5 Mpc. The cluster galaxies in the feedback model have a velocity dispersion approximately 15% lower than the dark matter. This results in the true ratio of specific energies in galaxies to gas being less than one, beta(sub spec) approximately 0.7. The infall model predicts beta(sub spec) approximately 1.2. Large excursions in these values occur over time, following the complex dynamical history of the cluster. The morphology of the X-ray emission is little affected by feedback. The emission profiles of both clusters are well described by the standard beta-model with beta(sub fit) approximately equal to 0.7 - 0.9. X-ray mass estimates based on the assumptions of hydrostatic equilibrium and the applicability of the beta-model are quite accurate in both cases. A strong, radial iron abundance gradient is present, which develops as a consequence of the steepening of the galaxy density profile over time. Spectroscopic observations using nonimaging detectors with wide (approximately 45 min) fields of view dramatically smear the gradient. Observations with arcminute resolution, made available with the ASCA satellite, would readily resolve the gradient.

Metzler, Christopher A.↗

Interactions between gravity waves and cold air outflows in a stably stratified uniform flow

Interactions between gravity waves and cold air outflows in a stably stratified uniform flow forced by various combinations of prescribed heat sinks and sources are studied using a hydrostatic two-dimensional nonlinear numerical model. The formation time for the development of a stagnation point or reversed flow at the surface is not always directly proportional to the Froude number when wave reflections exist from upper levels. A density current is able to form by the wave-otuflow interaction, even though the Froude number is greater than a critical value. This is the result of the wave-outflow interaction shifting the flow response to a different location in the characteristic parameter space. A density current is able to form or be destroyed due to the wave-outflow interaction between a traveling gravity wave and cold air outflow. This is proved by performing experiments with a steady-state heat sink and an additional transient heat source. In a quiescent fluid, a region of cold air, convergence, and upward motion is formed after the collision between two outflows produced by two prescribed heat sinks. After the collision, the individual cold air outflows lose their own identity and merge into a single, stationary, cold air outflow region. Gravity waves tend to suppress this new stationary cold air outflow after the collision. The region of upward motion associated with the collision is confined to a very shallow layer. In a moving airstream, a density current produced by a heat sink may be suppressed or enhanced nonlinearly by an adjacent heat sink due to the wave-outflow interaction.

Lin, Yuh-Lang↗

ROSAT PSPC observations of the early-type galaxies NGC 507 and NGC 499: Central cooling and mass determination

We present the results of a deep observation of NGC 507 and NGC 499 with the ROSAT Position Sensitive Proportional Counter (PSPC). The X-ray emission of NGC 507 is extended at least out to 1000 sec (458 kpc at a distance of 94.5 Mpc). The radial profile of X-ray surface brightness goes as Sigma(sub x) is approximately r(exp -1.8) outside the core region. The radial profile is a function of energy such that the softer X-rays have a smaller core radius and a flatter slope. Spectral analysis reveals that the emission temperature, with an average of 1 keV, peaks at an intermediate radius of 2-3 min and falls toward the center (possibly decreases outward as well). The absorption column density is consistent with the Galactic line-of-sight value. The X-ray emission of NGC 499 is extended to 300 sec and suggests a similarly cooler core. The cooler cores of NGC 507 and NGC 499 are strong evidence of the presence of cooling flows in these galaxies. Assuming hydrostatic equilibrium outside the cooling radius, the estimated mass-to-light ratio of NGC 507 is 97 +/- 16 within 458 kpc, indicative of the presence of a heavy halo. Similarly, the mass-to-light ratio of NGC 499 is 89 +/- 14 within 137 kpc. Near the edge of the X-ray-emitting region of NGC 507 we detect 19 soft, unresolved sources. These sources do not have optical counterparts and are significantly in excess of the expected number of background serendipitous sources. We speculate that they may represent cooling clumps in the halo of NGC 507. If there are many undetected cooling clumps distributed at large radii, then the radial profile of the X-ray surface brightness does not directly reflect the potential, adding uncertainty to the measurement of the binding mass; the gas mass could also be overestimated.

Kim, Dong-Woo↗

An X-ray temperature map of Abell 754: A major merger

We present the first two-dimensional X-ray temperature map of the rich cluster of galaxies Abell 754. We also present an X-ray surface brightness map with improved spatial resolution and sensitivity compared with previous maps. Both the temperature map and the surface brightness map show that A754 is in the throes of a violent merger; it is probably far from hydrostatic equilibrium.

Henry, J. Patrick↗

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