Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magnetic diffusion”

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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Naval Research Laboratory's programs in advanced indium phosphide solar cell development

The Naval Research Laboratory (NRL) has been involved in the development of solar cells for space applications since the 1960s. It quickly became apparent in this work that radiation damage caused to solar cells by electrons and protons trapped by the earth's magnetic field would seriously degrade the power output of photovoltaic arrays in extended missions. Techniques were therefore developed to harden the cells by shielding them with coverglass, etc. Ultimately, however, there is a limit to such approaches, which is determined by the radiation response of the semiconductor material employed. A desire for high efficiency and radiation resistance led to the development of alternative cell technologies such as GaAs, which has since become the technology of choice for many applications. InP cells are currently the most radiation resistant, high efficiency, planar cells known. NRL first sponsored InP solar cell technology in 1986, when Arizona State University was contracted to grow p/n cells by liquid phase epitaxy. NRL's interest in InP cells was generated by the results presented by Yamaguchi and his co-workers in the early 1980s on the remarkable radiation resistance of cells grown by diffusion of S into Zn doped p-type InP substrates. These cells also had beginning of life (BOL) efficiencies approximately 16%(AM0). Related to the radiation resistance of the cells was the fact that radiation-induced damage could be optically annealed by sunlight. Relatively large quantities of 1 x 2 cm(exp 2) diffused junction cells were made and were used on the MUSES-A and the EXOS-D satellites. These cells were also available in the U.S. through NIMCO, and were studied at NRL and elsewhere. Workers at NASA Lewis became involved in research in InP cells about the same time as NRL.

Summers, Geoffrey P.↗
NASA NTRS2011

Ion-Neutral Coupling in Solar Prominences

Interactions between ions and neutrals in a partially ionized plasma are important throughout heliophysics, including near the solar surface in prominences. Understanding how ion-neutral coupling affects formation, support, structure, and dynamics of prominences will advance our physical understanding of magnetized systems involving a transition from a weakly ionized dense gas to a fully ionized tenuous plasma. We address the fundamental physics of prominence support, which is normally described in terms of a magnetic force on the prominence plasma that balances the solar gravitational force, and the implications for observations. Because the prominence plasma is only partially ionized, it is necessary to consider the support of the both the ionized and neutral components. Support of the neutrals is accomplished through a frictional interaction between the neutral and ionized components of the plasma, and its efficacy depends strongly on the degree of ionization of the plasma. More specifically, the frictional force is proportional to the relative flow of neutral and ion species, and for a sufficiently weakly ionized plasma, this flow must be relatively large to produce a frictional force that balances gravity. A large relative flow, of course, implies significant draining of neutral particles from the prominence. We evaluate the importance of this draining effect for a hydrogen-helium plasma, and consider the observational evidence for cross-field diffusion of neutral prominence material.

Gilbert, Holly↗
NASA NTRS1973

Drift shell splitting by internal geomagnetic multipoles.

Computations on an 80-coefficient model of the earth's field illustrate the 'topography' of the magnetic equatorial surface and the geometry of the drift shells of geomagnetically trapped particles. Individual terms in the spherical harmonic expansion of the geomagnetic scalar potential V(r, theta, phi) are either even or odd in cos theta, where theta = 90 deg denotes the dipole equator. Terms that are even in cos theta tend to 'warp' the equatorial surface, but do not (in first order) distort particle drift shells radially nor split the drift shells of particles having different equatorial pitch angles. Azimuthally asymmetric terms that are odd in cos theta do cause shell splitting in first order. Shell splitting at large L values (neglecting deformation of the earth's field by the solar wind) is found to be dominated by the geomagnetic octupole. At L approximately equal to 1, shell splitting is strongly enhanced by the South American and South African anomalies. When combined with pitch angle diffusion caused by atmospheric scattering, these results may be able to account for anomalous radial diffusion of inner zone electrons.

Roederer, J. G.↗
NASA NTRS2000

Vertical Shear of the Galactic Interstellar Medium

The detection of UV absorption, 21 cm, H alpha and other diffuse optical emission lines from gas up to ten kiloparsecs above the plane of the Milky Way and other galaxies provides the first, opportunity to probe the rotational properties of the ionized "atmospheres" of galaxies. This rotation has implications for our understanding of the Galactic gravitational potential, angular momentum transport in the Galactic disk, and the maintenance of a Galactic dynamo. The available evidence indicates that gas rotates nearly cylindrically up to a few kiloparsecs. This is in contrast to the expectation that there should be a significant gradient in rotation speed as a function of height assuming a reasonable mass model for the Galaxy. For example, for a vertical cut at galactocentric radius R = 5 kpc in NGC 891 by Rand, the rotation speed is observed to drop by approximately 30 kilometers per second from z = 1 to 5 kpc and is expected to drop by 80 kilometers per second. Magnetic tension forces may resolve this discrepancy. Other possibilities will be examined in the near future.

Benjamin, Robert A.↗
NASA NTRS1999

Low-Frequency Oscillations and Transport Processes Induced by Multiscale Transverse Structures in the Polar Wind Outflow: A Three-Dimensional Simulation

Multiscale transverse structures in the magnetic-field-aligned flows have been frequently observed in the auroral region by FAST and Freja satellites. A number of multiscale processes, such as broadband low-frequency oscillations and various cross-field transport effects are well correlated with these structures. To study these effects, we have used our three-dimensional multifluid model with multiscale transverse inhomogeneities in the initial velocity profile. Self-consistent-frequency mode driven by local transverse gradients in the generation of the low field-aligned ion flow and associated transport processes were simulated. Effects of particle interaction with the self-consistent time-dependent three-dimensional wave potential have been modeled using a distribution of test particles. For typical polar wind conditions it has been found that even large-scale (approximately 50 - 100 km) transverse inhomogeneities in the flow can generate low-frequency oscillations that lead to significant flow modifications, cross-field particle diffusion, and other transport effects. It has also been shown that even small-amplitude (approximately 10 - 20%) short-scale (approximately 10 km) modulations of the original large-scale flow profile significantly increases low-frequency mode generation and associated cross-field transport, not only at the local spatial scales imposed by the modulations but also on global scales. Note that this wave-induced cross-field transport is not included in any of the global numerical models of the ionosphere, ionosphere-thermosphere, or ionosphere-polar wind. The simulation results indicate that the wave-induced cross-field transport not only affects the ion outflow rates but also leads to a significant broadening of particle phase-space distribution and transverse particle diffusion.

Ganguli, Supriya B.↗
NASA NTRS1992

The lunar environment and its effect on optical astronomy

The Moon's geologic environment features: (1) gravity field one-sixth that of Earth; (2) sidereal rotation period of 27.3 days; (3) surface with greater curvature than Earth's surface (a chord along a 10 km baseline would have a bulge of 7.2 m); (4) seismically and tidally stable platform on which to make astronomical observations (most moonquakes have magnitudes of 1 to 2 on the Richter scale, within the earth's seismic noise, resulting in ground motions only 1 nm); (5) tenuous atmosphere (the total mass at night is only 10(exp 4) kg) that has an optical depth of 10(exp -6) and does not cause wind induced stresses and vibrations on structures; (6) large diurnal temperature variation (100 to 385 K in equatorial regions), which telescopes must be designed to withstand; (7) weak magnetic field, ranging from 3 to 330 x 10(exp -9) T, compared to 3 x 10(exp -5) T on Earth at the equator; (8) surface exposed to radiation, the most dangerous of which are high energy (1 to 100 Mev) particles resulting from solar flares; (9) high flux of micrometeorites which are not slowed down from their cosmic velocities because of the lack of air (data indicate that microcraters greater than 10 microns across will form at the rate of 3000/sq m/yr); (10) regolith 2 to 30 m thick which blankets the entire lunar surface (this layer is fine-grained (average grain sizes range from 40 to 268 microns), has a low density (800 to 1000 kg/cu m in the upper few mm, rising to 1500 to 1800 kg/cu m at depths of 10 to 20 cm), is porous (35 to 45 pct), cohesive (0.1 to 1.0 kN/sq m), and has a low thermal diffusivity (0.7 to 1.0 x 110-8 sq m/sec); about 29 pct of the regolith is less than 20 micron in size (this dust could pose a hazard to optical telescopes); (11) rubbly upper several hundred meters in which intact bedrock is uncommon, especially in the lunar highlands; and (12) craters with diameter-to-depth ratios of 5 if fresh and less than km across (larger and eroded craters have diameter-to-depth ratios greater than 5).

Taylor, G. Jeffrey↗
NASA NTRS2013

High-Energy X-rays from J174545.5-285829, the Cannonball: a Candidate Pulsar Wind Nebula Associated with Sgr a East

We report the unambiguous detection of non-thermal X-ray emission up to 30 keV from the Cannonball, a few arcsecond long diffuse X-ray feature near the Galactic Center, using the NuSTAR X-ray observatory. The Cannonball is a high-velocity (v(proj) approximately 500 km s(exp −1)) pulsar candidate with a cometary pulsar wind nebula (PWN) located approximately 2' north-east from Sgr A*, just outside the radio shell of the supernova remnant Sagittarius A (Sgr A) East. Its non-thermal X-ray spectrum, measured up to 30 keV, is well characterized by a Gamma is approximately 1.6 power law, typical of a PWN, and has an X-ray luminosity of L(3-30 keV) = 1.3 × 10(exp 34) erg s(exp −1). The spectral and spatial results derived from X-ray and radio data strongly suggest a runaway neutron star born in the Sgr A East supernova event. We do not find any pulsed signal from the Cannonball. The NuSTAR observations allow us to deduce the PWN magnetic field and show that it is consistent with the lower limit obtained from radio observations.

stars: neutron↗
NASA NTRS1987

AMPTE IRM observations of waves associated with flux transfer events in the magnetosphere

The AMPTE IRM wave instrument has been applied to the study of flux transfer events (FTE's). This initial investigation concentrates on FTE's observed in the magnetosphere during the fall of 1984. The wave morphology consists of four significant features: at frequencies below the ion gyrofrequency, magnetic fluctuations occur with amplitudes of the order of 1 nT; at frequencies from a few hertz to a few hundred hertz, electric field fluctuations are observed which have a broadband amplitude of a few millivolts per meter, which are perpendicularly polarized at the lowest frequencies, and which are partly electrostatic and partly electromagnetic; in the frequency range from 100 Hz to 10 kHz, broadband spikelike waves occur with a time-averaged amplitude of about 0.1 mV/m; and near the electron plasma frequency, bursts of waves occur at the boundaries of FTE's. In none of the 25 events surveyed does the total broadband amplitude of all the waves exceed a few millivolts per meter, an amplitude far too small to provide the dissipation required by 'traditional' reconnection models in which the dissipation region has a thickness the order of an ion gyroradius. Thus, either all of the FTE's are observed at some distance from the diffusion region, or the observable waves play no significant role in the diffusion process, or 'traditional' reconnection models do not apply to reconnection in FTE's.

Labelle, J.↗
NASA NTRS1993

The origin of chromitic chondrules and the volatility of Cr under a range of nebular conditions

We characterize ten chromatic chondrules, two spinelian chondrules andd one spinel-bearing chondrule and summarize data for 120 chromitic inclusions discovered in an extensive survey of ordinary chondrites. Compositional and petrographic evidence suggests that chromitic chondrules and inclusions are closely related. The Cr/(Cr + Al) ratios in the spinal of these objects range from 0.5 to 0.9 and bulk Al2O3 contents are uniformly high (greater than 10 wt%, except for one with 8 wt%). No other elements having comparable solar abundances are so stongly enriched, and alkali feldspar and merrillite are more common than in normal chondrules. The Cr/Mg ratios in chromitic chondrules are 180-750 times the ratios in the bulk chondrite. With the possible exception of magnetic clumping of chromite in the presolar cloud, mechanical processes cannot account for this enrichment. Examination of nebular equilibrium processes shows that 50%-condensation temperatures of Cr at pH2/pH2O of 1500 are several tens of degrees below those of Mg as Mg2SiO4; the condensation of Cr is primarily as MgCr2O4 dissolved in MgAl2O4 at nebular pressures of 10(exp -4) atm or below. At pH2 = 10(exp -3) atm condesation as Cr in Fe-Ni is favored. Making the nebula much more oxidizing reduces the difference in condensation temperatures but Mg remains more refractory. We conclude that nebular equilibrium processes are not responsible for the enhanced Cr/Mg ratios. We propose that both Cr and Al became enriched in residues formed by incomplete evaporation of presolar lumps. We suggest that spinals remained as solid phases when the bulk of the silicates were incorporated into the evaporating melt; vaporization of Al and Cr were inhibited by the slow kinetics of diffusion. Subsequent melting and crystallization of these residues fractionated Cr from Al. The resulting materials constituted major components in the precursors of chromitic chondrules. Our model implies that chromitic chondrules and inclusions preserve the Cr isotopic record of presolar sources.

Krot, Alexander↗
NASA NTRS2002

X-ray Spectra and Pulse Frequency Changes in SAX J2103.5+4545

The November 1999 outburst of the transient pulsar SAX J2103.5+4545 was monitored with the large area detectors of the Rossi X-Ray Timing Explorer until the pulsar faded after a year. The 358 s pulsar was spun up for 150 days, at which point the flux dropped quickly by a factor of approximately 7, the frequency saturated and, as the flux continued to decline, a weak spin-down began. The pulses remained strong during the decay and the spin-up/flux correlation can be fit to the Ghosh and Lamb derivations for the spin-up caused by accretion from a thin, pressure-dominated disk, for a distance approximately 3.2 kpc and a surface magnetic field approximately 1.2 x 10(exp 13) Gauss. During the bright spin-up part of the outburst, the flux was subject to strong orbital modulation, peaking approximately 3 days after periastron of the eccentric 12.68 day orbit, while during the faint part, there was little orbital modulation. The X-ray spectra were typical of accreting pulsars, describable by a cut-off power-law, with an emission line near the 6.4 keV of Kappa(sub alpha) fluorescence from cool iron. The equivalent width of this emission did not share the orbital modulation, but nearly doubled during the faint phase, despite little change in the column density. The outburst could have been caused by an episode of increased wind from a Be star, such that a small accretion disk is formed during each periastron passage. A change in the wind and disk structure apparently occurred after 5 months such that the accretion rate was no longer modulated or the diffusion time was longer. The distance estimate implies the X-ray luminosity observed was between 1 X 10(exp 36) ergs s(exp -1) and 6 x 10(exp 34) ergs s(exp -1), with a small but definite correlation of the intrinsic power-law spectral index.

Baykal, A.↗
NASA NTRS2019

Energetic Particles of Cosmic Accelerators I: Galactic Accelerators

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Our pursuit of more than a century to uncover the origins and fate of these cosmic energetic particles has given rise to some of the most interesting and challenging questions in astrophysics. Energetic particles in our own galaxy, galactic cosmic rays (GCRs), engage in a complex interplay with the interstellar medium and magnetic fields in the galaxy, giving rise to many of its key characteristics. For instance, GCRs act in concert with galactic magnetic fields to support its disk against its own weight. GCR ionization and heating are essential ingredients in promoting and regulating the formation of stars and protostellar disks. GCR ionization also drives astrochemistry, leading to the build up of complex molecules in the interstellar medium. GCR transport throughout the galaxy generates and maintains turbulence in the interstellar medium, alters its multi-phase structure, and amplifies magnetic fields. GCRs could even launch galactic winds that enrich the circumgalactic medium and alter the structure and evolution of galactic disks. As crucial as they are for many of the varied phenomena in our galaxy, there is still much we do not understand about GCRs. While they have been linked to supernova remnants (SNRs), it remains unclear whether these objects can fully account for their entire population, particularly at the lower (approximately less than 1 GeV per nucleon) and higher (~PeV) ends of the spectrum. In fact, it is entirely possible that the SNRs that have been found to accelerate CRs merely re-accelerate them, leaving the origins of the original GCRs a mystery. The conditions for particle acceleration that make SNRs compelling source candidates are also likely to be present in sources such as protostellar jets, superbubbles, and colliding wind binaries (CWBs), but we have yet to ascertain their roles in producing GCRs. For that matter, key details of diffusive shock acceleration (DSA) have yet to be revealed, and it remains to be seen whether DSA can adequately explain particle acceleration in the cosmos. This White Paper is the first of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. For the case of GCRs, MeV astronomy will: 1) Search for fresh acceleration of GCRs in SNRs; 2) Test the DSA process, particularly in SNRs and CWBs; 3) Search for signs of CR acceleration in protostellar jets and superbubbles.

Venters, Tonia M.↗
NASA NTRS1994

Time-dependent simulation of cosmic-ray shocks, including Alfven transport

Time evolution of plane, cosmic-ray modified shocks has been simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven-wave energy transport equations. Nonlinear feedback from the cosmic rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier, Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the timescale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗
NASA NTRS1994

The roots of coronal structure in the Sun's surface

We have compared the structures seen on X-ray images obtained by a flight of the NIXT sounding rocket payload on July 11, 1991 with near-simultaneous photospheric and chromospheric structures and magnetic fields observed at Big Bear. The X-ray images reflect emission of both Mg X and Fe XVI, formed at 1 x 10(exp 6) K and 3 x 10(exp 6) K, respectively. The brightest H(alpha) sources correspond to a dying sub-flare and other active region components, all of which reveal coronal enhancements situated spatially well above the H(alpha) emission. The arches appear to lie in a small range of angle in the meridian plane connecting their footpoints. Sunspots are dark on the surface and in the corona. For the first time we see an emerging flux region in X-rays and find the emission extends twice as high as the H(alpha) arches. Many features which we believe to correspond to `X-ray bright points' (XBPs) were observed. Whether by resolution or spectral band, the number detected greatly exceeds that from previous work. All of the brighter XBPs correspond to bipolar H(alpha) features, while unipolar H(alpha) bright points are the base of more diffuse comet-like coronal arches, generally vertical. These diverge from individual features by less than 30 deg, and give a good measure of what the `canopies' must do. The H(alpha) data shows that all the H(alpha) features were present present the entire day, so they are not clearly disappearing or reappearing. We find a new class of XBPs which we call `satellite points', elements of opposite polarity linked to nearby umbrae by invisible field lines. The satellite points change rapidly in X-ray brightness during the flight. An M1.9 flare occurred four hours after the flight; examination of the pre-flare structures reveals nothing unusual.

Golub, Leon↗
NASA NTRS1997

Numerical Simulation of Rotation-Driven Plasma Transport In the Jovian Magnetosphere

A Jupiter version of the Rice Convection Model (RCM-J) was developed with support of an earlier NASA SR&T grant. The conversion from Earth to Jupiter included adding currents driven by centrifugal force, reversing the planetary magnetic field, and rescaling various parameters. A series of informative runs was carried out, all of them solving initial value problems. The simulations followed an initial plasma torus configuration as it fell apart by interchange instability. Some conclusions from the simulations were the following: 1. We confirmed that, for conventional values of the torus density and ionospheric conductance, the torus disintegrates by interchange instability on a time scale of approx. one day, which is 1-2 orders of magnitude shorter than the best estimates of the average residence time of plasma in the torus. 2. In the model, the instability could be slowed to an arbitrary degree by the addition of sufficient impounding energetic particles, as suggested earlier by Siscoe et al (1981). However, the observed energetic particles do not seem sufficient to guarantee impoundment (e.g., Mauk et al., 1996). 3. Whether inhibited by impoundment or not, the interchange was found to proceed by the formation of long fingers, which get thinner as they get longer. This picture differed dramatically from the conventional radial-diffusion picture (e.g., Siscoe and Summers (1981)), more superficially with the outward-moving-blob picture (Pontius and Hill, 1989). The obvious limitation of the original RCM-J was that it could not represent a plasma source. We could represent the decay of a pre-existing torus, but we could not represent the way ionization of material from Io continually replenishes the plasma. We consequently were precluded from studying a whole set of fundamental issues of torus theory, including whether the system can come to a steady state.

Wolf, Richard A.↗
NASA NTRS1993

Time dependent simulation of cosmic-ray shocks including Alfven transport

Time evolution of plane, cosmic-ray modified shocks was simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven wave energy transport equations. Nonlinear feedback from the cosmic-rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the time scale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic-rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗
NASA NTRS2017

Kinetic Theory and Fast Wind Observations of the Electron Strahl

We develop a model for the strahl population in the solar wind - a narrow, low-density and high-energy electron beam centred on the magnetic field direction. Our model is based on the solution of the electron drift-kinetic equation at heliospheric distances where the plasma density, temperature and the magnetic field strength decline as power laws of the distance along a magnetic flux tube. Our solution for the strahl depends on a number of parameters that, in the absence of the analytic solution for the full electron velocity distribution function (eVDF), cannot be derived from the theory. We however demonstrate that these parameters can be efficiently found from matching our solution with observations of the eVDF made by the Wind satellite's SWE strahl detector. The model is successful at predicting the angular width (FWHM) of the strahl for the Wind data at 1 au, in particular by predicting how this width scales with particle energy and background density. We find that the strahl distribution is largely determined by the local temperature Knudsen number γ ∼ |T dT/dx|/n, which parametrizes solar wind collisionality. We compute averaged strahl distributions for typical Knudsen numbers observed in the solar wind, and fit our model to these data. The model can be matched quite closely to the eVDFs at 1 au; however, it then overestimates the strahl amplitude at larger heliocentric distances. This indicates that our model may be improved through the inclusion of additional physics, possibly through the introduction of 'anomalous diffusion' of the strahl electrons.

Horaites, Konstantinos↗
NASA NTRS2015

The North Galactic Pole Rift and the Local Hot Bubble

The North Galactic Pole Rift (NGPR) is one of the few distinct neutral hydrogen clouds at high Galactic latitudes that have well-defined distances. It is located at the edge of the Local Cavity (LC) and provides an important test case for understanding the Local Hot Bubble (LHB), the presumed location for the hot diffuse plasma responsible for much of the observed 1/4 keV emission originating in the solar neighborhood. Using data from the ROSAT All- Sky Survey and the Planck reddening map, we find the path length within the LC (LHB plus Complex of Local Interstellar Clouds) to be 98 plus or minus 27 pc, in excellent agreement with the distance to the NGPR of 98 +/- 6 pc. In addition, we examine another 14 directions that are distributed over the sky where the LC wall is apparently optically thick at 1/4 keV. We find that the data in these directions are also consistent with the LHB model and a uniform emissivity plasma filling most of the LC.

ISM: bubbles aEuro" ISM: clouds aEuro" ISM: magnet↗
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↗