Search NASA⌕ Search

SEARCH · Search NASA

Results for “AMBIPOLAR 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Formation of magnetized prestellar cores with ambipolar diffusion and turbulence

We investigate the roles of magnetic fields and ambipolar diffusion during prestellar core formation in turbulent giant molecular clouds, using three-dimensional numerical simulations. Our simulations focus on the shocked layer produced by a converging large-scale flow and survey varying ionization and the angle between the upstream flow and magnetic field. We also include ideal magnetohydrodynamic (MHD) and hydrodynamic models. From our simulations, we identify hundreds of self-gravitating cores that form within 1 Myr, with masses M ∼ 0.04-2.5 M {sub ☉} and sizes L ∼ 0.015-0.07 pc, consistent with observations of the peak of the core mass function. Median values are M = 0.47 M {sub ☉} and L = 0.03 pc. Core masses and sizes do not depend on either the ionization or upstream magnetic field direction. In contrast, the mass-to-flux ratio does increase with lower ionization, from twice to four times the critical value. The higher mass-to-flux ratio for low ionization is the result of enhanced transient ambipolar diffusion when the shocked layer first forms. However, ambipolar diffusion is not necessary to form low-mass supercritical cores. For ideal MHD, we find similar masses to other cases. These masses are one to two orders of magnitude lower than the value M {sub mag,} {sub sph} = 0.007B {sup 3}/(G {sup 3/2}ρ{sup 2}) that defines a magnetically supercritical sphere under post-shock ambient conditions. This discrepancy is the result of anisotropic contraction along field lines, which is clearly evident in both ideal MHD and diffusive simulations. We interpret our numerical findings using a simple scaling argument that suggests that gravitationally critical core masses will depend on the sound speed and mean turbulent pressure in a cloud, regardless of magnetic effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

GMC Collisions as Triggers of Star Formation. IV. The Role of Ambipolar Diffusion

We investigate the role of ambipolar diffusion (AD) in collisions between magnetized giant molecular clouds (GMCs), which may be an important mechanism for triggering star cluster formation. Three-dimensional simulations of GMC collisions are performed using a version of the Enzo magnetohydrodynamics code that has been extended to include AD. The resistivities are calculated using the 31-species chemical model of Wu et al. (2015). We find that in the weak-field, 10 μG case, AD has only a modest effect on the dynamical evolution during the collision. However, for the stronger-field, 30 μG case involving near-critical clouds, AD results in the formation of dense cores in regions where collapse is otherwise inhibited. The overall efficiency of formation of cores with n{sub H}⩾10{sup 6} cm{sup −3} in these simulations is increases from about 0.2% to 2% once AD is included, comparable to observed values in star-forming GMCs. The gas around these cores typically has relatively slow infall at speeds that are a modest fraction of the free-fall speed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Global simulations of protoplanetary disks with OHMIC resistivity and ambipolar diffusion

Protoplanetary disks (PPDs) are believed to accrete onto their central T Tauri star because of magnetic stresses. Recently published shearing box simulations indicate that Ohmic resistivity, ambipolar diffusion (AD) and the Hall effect all play important roles in disk evolution. In the presence of a vertical magnetic field, the disk remains laminar between 1–5 AU, and a magnetocentrifugal disk wind forms that provides an important mechanism for removing angular momentum. Questions remain, however, about the establishment of a true physical wind solution in the shearing box simulations because of the symmetries inherent in the local approximation. We present global MHD simulations of PPDs that include Ohmic resistivity and AD, where the time-dependent gas-phase electron and ion fractions are computed under FUV and X-ray ionization with a simplified recombination chemistry. Our results show that the disk remains laminar, and that a physical wind solution arises naturally in global disk models. The wind is sufficiently efficient to explain the observed accretion rates. Furthermore, the ionization fraction at intermediate disk heights is large enough for magneto-rotational channel modes to grow and subsequently develop into belts of horizontal field. Depending on the ionization fraction, these can remain quasi-global, or break-up into discrete islands of coherent field polarity. The disk models we present here show a dramatic departure from our earlier models including Ohmic resistivity only. It will be important to examine how the Hall effect modifies the evolution, and to explore the influence this has on the observational appearance of such systems, and on planet formation and migration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Dust transport in MRI turbulent disks: Ideal and non-ideal MHD with ambipolar diffusion

We study dust transport in turbulent protoplanetary disks using three-dimensional global unstratified MHD simulations including Lagrangian dust particles. The turbulence is driven by the magnetorotational instability (MRI) with either ideal or non-ideal MHD that includes ambipolar diffusion (AD). In ideal MHD simulations, the surface density evolution (except for dust that drifts fastest), turbulent diffusion, and vertical scale height of dust can all be reproduced by simple one-dimensoinal and/or analytical models. However, in AD dominated simulations which simulate protoplanetary disks beyond 10s of AU, the vertical scale height of dust is larger than previously predicted. To understand this anomaly in more detail, we carry out both unstratified and stratified local shearing box simulations with Lagrangian particles, and find that turbulence in AD dominated disks has very different properties (e.g., temporal autocorrelation functions and power spectra) than turbulence in ideal MHD disks, which leads to quite different particle diffusion efficiency. For example, MRI turbulence with AD has a longer correlation time for the vertical velocity, which causes significant vertical particle diffusion and large dust scale height. In ideal MHD the Schmidt numbers (Sc) for radial and vertical turbulent diffusion are Sc{sub r}∼1 and Sc{sub z}≳3, but in the AD dominated regime both Sc{sub r} and Sc{sub z} are ≲1. Particle concentration in pressure bumps induced by MRI turbulence has also been studied. Since non-ideal MHD effects dominate most regions in protoplanetary disks, our study suggests that modeling dust transport in turbulence driven by MRI with non-ideal MHD effects is important for understanding dust transport in realistic protoplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Alfvénic turbulence beyond the ambipolar diffusion scale

We investigate the nature of the Alfvénic turbulence cascade in two-fluid magnetohydrodynamic (MHD) simulations in order to determine if turbulence is damped once the ion and neutral species become decoupled at a critical scale called the ambipolar diffusion scale (L{sub AD}). Using mode decomposition to separate the three classical MHD modes, we study the second-order structure functions of the Alfvén mode velocity field of both neutrals and ions in the reference frame of the local magnetic field. On scales greater than L{sub AD} we confirm that two-fluid turbulence strongly resembles single-fluid MHD turbulence. Our simulations show that the behavior of two-fluid turbulence becomes more complex on scales less than L{sub AD}. We find that Alfvénic turbulence can exist past L{sub AD} when the turbulence is globally super-Alfvénic, with the ions and neutrals forming separate cascades once decoupling has taken place. When turbulence is globally sub-Alfvénic and hence strongly anisotropic, with a large separation between the parallel and perpendicular decoupling scales, turbulence is damped at L{sub AD}. We also find that the power spectrum of the kinetic energy in the damped regime is consistent with a k{sup −4} scaling (in agreement with the predictions of Lazarian et al.).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Ambipolar and non-ambipolar diffusion in an rf plasma source containing a magnetic filter

By placing a magnetic filter across a rectangular plasma source (closed at one end with a ceramic plate and an rf antenna, and terminated at the opposite end by a grounded grid), we experimentally investigate the effect of conducting and insulating source walls on the nature of the plasma diffusion phenomena. The use of a magnetic filter creates a unique plasma, characterized by a high upstream electron temperature (T{sub e{sub u}}∼5 eV) near the rf antenna and a low downstream electron temperature (T{sub e{sub d}}∼1 eV) near the grid, which more clearly demonstrates the role of the source wall materials. For conducting walls a net ion current to ground is measured on the grid, and the plasma potential is determined by a mean electron temperature within the source. For insulating walls the plasma potential is determined by the downstream electron temperature (i.e., V{sub p}∼5.2T{sub e{sub d}} in argon), and the net current to the grid is exactly zero. Furthermore, by inserting a small additional upstream conductor (that can be made floating or grounded through an external circuit switch), we demonstrate that the plasma potential can be controlled and set to a low (V{sub p}∼5.2T{sub e{sub d}}), or high (V{sub p}∼5.2T{sub e{sub u}}) value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Two-dimensional Radiative Magnetohydrodynamic Simulations of Partial Ionization in the Chromosphere. II. Dynamics and Energetics of the Low Solar Atmosphere

We investigate the effects of interactions between ions and neutrals on the chromosphere and overlying corona using 2.5D radiative MHD simulations with the Bifrost code. We have extended the code capabilities implementing ion–neutral interaction effects using the generalized Ohm’s law, i.e., we include the Hall term and the ambipolar diffusion (Pedersen dissipation) in the induction equation. Our models span from the upper convection zone to the corona, with the photosphere, chromosphere, and transition region partially ionized. Our simulations reveal that the interactions between ionized particles and neutral particles have important consequences for the magnetothermodynamics of these modeled layers: (1) ambipolar diffusion increases the temperature in the chromosphere; (2) sporadically the horizontal magnetic field in the photosphere is diffused into the chromosphere, due to the large ambipolar diffusion; (3) ambipolar diffusion concentrates electrical currents, leading to more violent jets and reconnection processes, resulting in (3a) the formation of longer and faster spicules, (3b) heating of plasma during the spicule evolution, and (3c) decoupling of the plasma and magnetic field in spicules. Our results indicate that ambipolar diffusion is a critical ingredient for understanding the magnetothermodynamic properties in the chromosphere and transition region. The numerical simulations have been made publicly available, similar to previous Bifrost simulations. This will allow the community to study realistic numerical simulations with a wider range of magnetic field configurations and physics modules than previously possible.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Generation of magnetic structures on the solar photosphere

The lower solar atmosphere is a partially ionized plasma consisting of electrons, ions, and neutral atoms. In this, which is essentially a three-fluid system, the Hall effect arises from the treatment of the electrons and ions as two separate fluids and the ambipolar diffusion arises from the inclusion of neutrals as the third fluid. The Hall effect and ambipolar diffusion have been shown to be operational in a region beginning from near the photosphere up to the chromosphere. In a partially ionized plasma, the magnetic induction is subjected to ambipolar diffusion and the Hall drift in addition to the usual resistive dissipation. These nonlinear effects create sharp magnetic structures which then submit themselves to various relaxation mechanisms. A first-principles derivation of these effects in a three-fluid system and an analytic solution to the magnetic induction equation in a stationary state are presented, which in the general case includes the Hall effect, ambipolar diffusion, and ohmic dissipation. The temporal evolution of the magnetic field is then investigated under the combined as well as the individual effects of the Hall drift and ambipolar diffusion to demonstrate the formation of steep magnetic structures and the resultant current sheet formation. These structures have just the right features for the release of magnetic energy into the solar atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

ON THE MISALIGNMENT BETWEEN CHROMOSPHERIC FEATURES AND THE MAGNETIC FIELD ON THE SUN

Observations of the upper chromosphere show an enormous amount of intricate fine structure. Much of this comes in the form of linear features, which are most often assumed to be well aligned with the direction of the magnetic field in the low plasma β regime that is thought to dominate the upper chromosphere. We use advanced radiative magnetohydrodynamic simulations, including the effects of ion-neutral interactions (using the generalized Ohm’s law) in the partially ionized chromosphere, to show that the magnetic field is often not well aligned with chromospheric features. This occurs where the ambipolar diffusion is large, i.e., ions and neutral populations decouple as the ion-neutral collision frequency drops, allowing the field to slip through the neutral population; where currents perpendicular to the field are strong; and where thermodynamic timescales are longer than or similar to those of ambipolar diffusion. We find this often happens in dynamic spicule or fibril-like features at the top of the chromosphere. This has important consequences for field extrapolation methods, which increasingly use such upper chromospheric features to help constrain the chromospheric magnetic field: our results invalidate the underlying assumption that these features are aligned with the field. In addition, our results cast doubt on results from 1D hydrodynamic models, which assume that plasma remains on the same field lines. Finally, our simulations show that ambipolar diffusion significantly alters the amount of free energy available in the coronal part of our simulated volume, which is likely to have consequences for studies of flare initiation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

RADIATION MAGNETOHYDRODYNAMIC SIMULATIONS OF PROTOSTELLAR COLLAPSE: NONIDEAL MAGNETOHYDRODYNAMIC EFFECTS AND EARLY FORMATION OF CIRCUMSTELLAR DISKS

The transport of angular momentum by magnetic fields is a crucial physical process in the formation and evolution of stars and disks. Because the ionization degree in star-forming clouds is extremely low, nonideal magnetohydrodynamic (MHD) effects such as ambipolar diffusion and ohmic dissipation work strongly during protostellar collapse. These effects have significant impacts in the early phase of star formation as they redistribute magnetic flux and suppress angular momentum transport by magnetic fields. We perform three-dimensional nested-grid radiation magnetohydrodynamic simulations including ohmic dissipation and ambipolar diffusion. Without these effects, magnetic fields transport angular momentum so efficiently that no rotationally supported disk is formed even after the second collapse. Ohmic dissipation works only in a relatively high density region within the first core and suppresses angular momentum transport, enabling formation of a very small rotationally supported disk after the second collapse. With both ohmic dissipation and ambipolar diffusion, these effects work effectively in almost the entire region within the first core and significant magnetic flux loss occurs. As a result, a rotationally supported disk is formed even before a protostellar core forms. The size of the disk is still small, about 5 AU at the end of the first core phase, but this disk will grow later as gas accretion continues. Thus, the nonideal MHD effects can resolve the so-called magnetic braking catastrophe while keeping the disk size small in the early phase, which is implied from recent interferometric observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Induced core formation time in subcritical magnetic clouds by large-scale trans-Alfvénic flows

We clarify the mechanism of accelerated core formation by large-scale nonlinear flows in subcritical magnetic clouds by finding a semi-analytical formula for the core formation time and describing the physical processes that lead to them. Recent numerical simulations show that nonlinear flows induce rapid ambipolar diffusion that leads to localized supercritical regions that can collapse. Here, we employ non-ideal magnetohydrodynamic simulations including ambipolar diffusion for gravitationally stratified sheets threaded by vertical magnetic fields. One of the horizontal dimensions is eliminated, resulting in a simpler two-dimensional simulation that can clarify the basic process of accelerated core formation. A parameter study of simulations shows that the core formation time is inversely proportional to the square of the flow speed when the flow speed is greater than the Alfvén speed. We find a semi-analytical formula that explains this numerical result. The formula also predicts that the core formation time is about three times shorter than that with no turbulence, when the turbulent speed is comparable to the Alfvén speed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Observation and interpretation of energy efficient, diffuse direct current glow discharge at atmospheric pressure

A diffuse direct-current glow discharge was realized with low energy consumption and high energy utilization efficiency at atmospheric pressure. The formation of diffuse discharge was demonstrated by examining and comparing the electrical properties and optical emissions of plasmas. In combination with theoretical derivation and calculation, we draw guidelines that appearance of nitrogen ions at low electron density is crucial to enhance the ambipolar diffusion for the expansion of discharge channel and the increasing ambipolar diffusion near the cathode plays a key role in the onset of diffuse discharge. An individual-discharge-channel expansion model is proposed to explain the diffuse discharge formation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

HALL EFFECT CONTROLLED GAS DYNAMICS IN PROTOPLANETARY DISKS. II. FULL 3D SIMULATIONS TOWARD THE OUTER DISK

We perform three-dimensional stratified shearing-box magnetohydrodynamic (MHD) simulations on the gas dynamics of protoplanetary disks with a net vertical magnetic flux of B {sub z0}. All three nonideal MHD effects, Ohmic resistivity, the Hall effect, and ambipolar diffusion, are included in a self-consistent manner based on equilibrium chemistry. We focus on regions toward outer disk radii, from 5 to 60 AU, where Ohmic resistivity tends to become negligible, ambipolar diffusion dominates over an extended region across the disk height, and the Hall effect largely controls the dynamics near the disk midplane. We find that at around R = 5 AU the system launches a laminar or weakly turbulent magnetocentrifugal wind when the net vertical field B {sub z0} is not too weak. Moreover, the wind is able to achieve and maintain a configuration with reflection symmetry at the disk midplane. The case with anti-aligned field polarity (Ω⋅B{sub z0}<0) is more susceptible to the magnetorotational instability (MRI) when B {sub z0} decreases, leading to an outflow oscillating in radial directions and very inefficient angular momentum transport. At the outer disk around and beyond R = 30 AU, the system shows vigorous MRI turbulence in the surface layer due to far-UV ionization, which efficiently drives disk accretion. The Hall effect affects the stability of the midplane region to the MRI, leading to strong/weak Maxwell stress for aligned/anti-aligned field polarities. Nevertheless, the midplane region is only very weakly turbulent in both cases. Overall, the basic picture is analogous to the conventional layered accretion scenario applied to the outer disk. In addition, we find that the vertical magnetic flux is strongly concentrated into thin, azimuthally extended shells in most of our simulations beyond 15 AU, leading to enhanced radial density variations know as zonal flows. Theoretical implications and observational consequences are briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Gravito-turbulence and dynamo in poorly ionized protostellar discs – I. Zero-net-flux case

In their early stages, protoplanetary discs are sufficiently massive to undergo gravitational instability (GI). This instability is thought to be involved in mass accretion, planet formation via gas fragmentation, the generation of spiral density waves, and outbursts. A key and very recent area of research is the interaction between the GI and magnetic fields in young protoplanetary discs, in particular whether this instability is able to sustain a magnetic field via a dynamo. We conduct 3D, stratified shearing-box simulations using two independent codes, PLUTO and Athena++, to characterize the GI dynamo in poorly ionized protostellar discs subject to ambipolar diffusion. We find that the dynamo operates across a large range of ambipolar Elssaser number Am (which characterizes the strength of ambipolar diffusion) and is particularly strong in the regime Am = 10–100, with typical magnetic to thermal energy ratios of order unity. The dynamo is only weakly dependent on resolution (at least for Am ≲ 100), box size, and cooling law. The magnetic field is produced by the combination of differential rotation and large-scale vertical roll motions associated with spiral density waves. Our results have direct implications for the dynamo process in young protoplanetary discs and possibly some regions of active galactic nucleus discs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

What Causes the High Apparent Speeds in Chromospheric and Transition Region Spicules on the Sun?

Spicules are the most ubuiquitous type of jets in the solar atmosphere. The advent of high-resolution imaging and spectroscopy from the Interface Region Imaging Spectrograph ( IRIS ) and ground-based observatories has revealed the presence of very high apparent motions of order 100–300 km s{sup −1} in spicules, as measured in the plane of the sky. However, line of sight measurements of such high speeds have been difficult to obtain, with values deduced from Doppler shifts in spectral lines typically of order 30–70 km s{sup −1}. In this work, we resolve this long-standing discrepancy using recent 2.5D radiative MHD simulations. This simulation has revealed a novel driving mechanism for spicules in which ambipolar diffusion resulting from ion-neutral interactions plays a key role. In our simulation, we often see that the upward propagation of magnetic waves and electrical currents from the low chromosphere into already existing spicules can lead to rapid heating when the currents are rapidly dissipated by ambipolar diffusion. The combination of rapid heating and the propagation of these currents at Alfvénic speeds in excess of 100 km s{sup −1} leads to the very rapid apparent motions, and often wholesale appearance, of spicules at chromospheric and transition region temperatures. In our simulation, the observed fast apparent motions in such jets are actually a signature of a heating front, and much higher than the mass flows, which are of order 30–70 km s{sup −1}. Our results can explain the behavior of transition region “network jets” and the very high apparent speeds reported for some chromospheric spicules.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Effect of Radiation on Chromospheric Magnetic Reconnection: Reactive and Collisional Multi-fluid Simulations

We study magnetic reconnection under chromospheric conditions in five different ionization levels from 0.5% to 50% using a self-consistent two-fluid (ions + neutrals) model that accounts for compressibility, collisional effects, chemical inequilibrium, and anisotropic heat conduction. Results with and without radiation are compared, using two models for the radiative losses: an optically thin radiation loss function, and an approximation of the radiative losses of a plasma with photospheric abundances. The results without radiation show that reconnection occurs faster for the weakly ionized cases as a result of the effect of ambipolar diffusion and fast recombination. The tearing mode instability appears earlier in the low ionized cases and grows rapidly. We find that radiative losses have a stronger effect than was found in previous results as the cooling changes the plasma pressure and the concentration of ions inside the current sheet. This affects the ambipolar diffusion and the chemical equilibrium, resulting in thin current sheets and enhanced reconnection. The results quantify this complex nonlinear interaction by showing that a strong cooling produces faster reconnections than have been found in models without radiation. The results accounting for radiation show timescales and outflows comparable to spicules and chromospheric jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Magnetically self-regulated formation of early protoplanetary disks

The formation of protoplanetary disks during the collapse of molecular dense cores is significantly influenced by angular momentum transport, notably by the magnetic torque. In turn, the evolution of the magnetic field is determined by dynamical processes and non-ideal MHD effects such as ambipolar diffusion. Considering simple relations between various timescales characteristic of the magnetized collapse, we derive an expression for the early disk radius, r≃18 au (η{sub AD}/0.1s){sup 2/9}(B{sub z}/0.1G){sup −4/9}(M/0.1M{sub ⊙}){sup 1/3}, where M is the total disk plus protostar mass, η{sub AD} is the ambipolar diffusion coefficient, and B {sub z} is the magnetic field in the inner part of the core. This is significantly smaller than the disks that would form if angular momentum was conserved. The analytical predictions are confronted against a large sample of 3D, non-ideal MHD collapse calculations covering variations of a factor 100 in core mass, a factor 10 in the level of turbulence, a factor 5 in rotation, and magnetic mass-to-flux over critical mass-to-flux ratios 2 and 5. The disk radius estimates are found to agree with the numerical simulations within less than a factor 2. A striking prediction of our analysis is the weak dependence of circumstellar disk radii upon the various relevant quantities, suggesting weak variations among class-0 disk sizes. In some cases, we note the onset of large spiral arms beyond this radius.

79 ASTRONOMY AND ASTROPHYSICS↗