Search NASA⌕ Search

SEARCH · Search NASA

Results for “ideal MHD”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Asymptotic-preserving semi-implicit finite volume scheme for extended magnetohydrodynamics

A Finite Volume (FV) scheme is developed for solving the extended magnetohydrodynamic (XMHD) equations, yielding accurate results in the ideal, resistive, and Hall MHD limits. This is accomplished by first re-writing the XMHD equations such that it allows the algorithm to retain the use of ideal MHD Riemann solvers and the constrained transport method to preserve divergence-free magnetic fields. Incorporation of electron inertia and displacement current introduces additional numerical stiffness which motivates a semi-implicit FV scheme that re-formulates the XMHD model as a relaxation system. The equations are then advanced in time using an explicit 2nd-order Runge–Kutta scheme with operator splitting applied to the implicit source term updates at each sub-stage. For additional numerical stability, a density-dependent slope limiter is implemented to increase flux diffusivity at low density regions where non-ideal effects become significant. The algorithm is subsequently implemented in a scalable adaptive mesh refinement (AMR) framework. As the new algorithm retains many aspects of the ideal MHD formulations, it asymptotes naturally to the ideal MHD limit. Moreover, it shows promising results at the resistive and Hall MHD limits. This is verified against reference test problems for ideal, resistive and Hall MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1975

Acceleration of charged particles in oblique MHD shocks

The motion of high energy charged particles in ideal oblique MHD shocks has been studied extensively. The orbits of charged particles can be solved exactly from Lorentz force equation. It is found that a charged particle may cross the shock front many times before it leaves the shock and it is energized during crossings. The most energetic particles are those with initial pitch and phase near the boundary between crossing and reflected particles. The reflected particles cross the shock about twice as many times as the passing particles do. The energy gain increases for larger values of shock strength and speed. If reflection occurs, the energy gain is larger for smaller angles between the field and shock surface. The introduction of scattering results in a few more higher energy particles.

Chen, G.↗
DOE OSTI2020

Time-dependent relaxed magnetohydrodynamics: Inclusion of cross helicity constraint using phase-space action

In this work, a phase-space version of the ideal magnetohydrodynamic (MHD) Lagrangian is derived from first principles and shown to give a relabeling transformation when a cross-helicity constraint is added in Hamilton's Action Principle. A new formulation of time-dependent “relaxed” magnetohydrodynamics is derived using microscopic conservation of mass and macroscopic constraints on total magnetic helicity, cross helicity, and entropy under variations of density, pressure, fluid velocity, and magnetic vector potential. This gives Euler–Lagrange equations consistent with previous work on both ideal and relaxed MHD equilibria with flow, but generalizes the relaxation concept from statics to dynamics. The application of the new dynamical formalism is illustrated for short-wavelength linear waves, and the interface connection conditions for Multiregion Relaxed MHD (MRxMHD) are derived. The issue of whether E+u×B=0 should be a constraint is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

As a Matter of Tension: Kinetic Energy Spectra in MHD Turbulence

While magnetized turbulence is ubiquitous in many astrophysical and terrestrial systems, our understanding of even the simplest physical description of this phenomena, ideal magnetohydrodynamic (MHD) turbulence, remains substantially incomplete. As such, in this work we highlight the shortcomings of existing theoretical and phenomenological descriptions of MHD turbulence that focus on the joint (kinetic and magnetic) energy fluxes and spectra by demonstrating that treating these quantities separately enables fundamental insights into the dynamics of MHD turbulence. This is accomplished through the analysis of the scale-wise energy transfer over time within an implicit large eddy simulation of subsonic, super-Alfvénic MHD turbulence. Our key finding is that the kinetic energy spectrum develops a scaling of approximately k –4/3 in the stationary regime as magnetic tension mediates large-scale kinetic to magnetic energy conversion and significantly suppresses the kinetic energy cascade. This motivates a reevaluation of existing MHD turbulence theories with respect to a more differentiated modeling of the energy fluxes.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1999

On Europa's Magnetospheric Interaction: A MHD Simulation of the E4 Flyby

The global three-dimensional interaction of Europa with the Jovian magnetosphere is modeled by using a complete set of ideal magnetohydrodynamic (MHD) equations. The model accounts for exospheric mass loading, ion-neutral charge exchange, recombination, and a possible intrinsic dipole magnetic field of Europa. The single-fluid MHD equations are solved by using a modem, finite volume, higher-order, Godunov-type method on an adoptively refined unstructured grid, which allows detailed modeling of the region near Europa while still resolving both the upstream region and the satellite's wake. The magnetic field and plasma density measured during Galileo's EGA flyby of December 19, 1996, are reproduced reasonably well in the simulation. We find the agreement between the data and our model particularly convincing if we assume that the plasma velocity during the EGA flyby deviated from the nominal corotation direction by approximately 20 deg. Evidence from the Galileo energetic particle detector also supports this assumption. In this case, we can fit the data using a dipole with orientation close to that of an induced dipole arising from the interaction of a hypothetical conducting subsurface layer on Europa with the periodically changing magnetic field of Jupiter. However, the magnitude of the dipole in our model is somewhat smaller (70%) than that suggested by Khurana et al. The total mass loading and ion-neutral charge exchange rates are consistent with the estimates of Europa's atmosphere and ionosphere.

Kabin, K.↗
NASA NTRS2000

Two-Species, 3D, MHD Simulation of Europa's Interaction with Jupiter's Magnetosphere

The interaction of Europa with the Jovian a magnetosphere has been studied by using a two species in ideal magnetohydrodynamic (MHD) numerical model. This model considers the upstream plasma in the Jovian magnetosphere and the molecular oxygen ions in the ionosphere of Europa, separately. We present results a from simulation studies, which take into account impact ionization, recombination, and the effect of a possible induced dipole magnetic field of Europa. The total mass loading of the magnetospheric flow and the ionization frequency used in the model are consistent with the estimates of Europa's ionosphere and atmosphere. The multi-species MHD equations are solved by using a finite volume, high-order, Godunov-type method on an adoptively refined unstructured grid, which allows a detailed modeling of the region near Europa's surface, while still resolving both the upstream region and the satellite's wake. We have paid special attention to the wake of Europa, in order to be able to make comparisons with the Galileo's E4 flyby observations, as well as other model calculations. The calculated escape flux of a O2+ down the tail was found to be about 5.6 x 10(exp 25) s(sup -1).

Liu, Yifan↗
NASA NTRS1990

Magnetohydrodynamic modeling of the solar corona

The ideal and resistive magnetohydrodynamic (MHD) model is used to examine the dynamics and structure of the solar corona. When the coronal magnetic field is deformed by photospheric flow it can evolve to states that become unstable to ideal MHD modes. The nonlinear evolution of these instabilities can lead to the generation of current sheets, field line reconnection, and energy release. The disruption of an arcade field and the kinking of coronal loops is described. The braiding of the large-scale coronal field by convective photospheric motions develops fine-scale structure in the magnetic field and leads to the development of intense current filaments. The resistive dissipation of these currents can provide an efficient coronal heating mechanism.

Mikic, Zoran↗
DOE OSTI2020

Edge ICRF simulations in 3D geometry: From MHD equilibrium to coupling determination

We present in this work a consistent numerical scheme that allows the computation of 3D magnetic fields a nd 3D density profiles and their usage in ion cyclotron range of frequencies (ICRF) coupling simulations. We first utilize the PARVMEC code to compute the 3D free-boundary plasma equilibrium in the ideal magnetohydrodynamic (MHD) approximation. Since the PARVMEC solution is only defined within the last closed flux surface (LCFS), the magnetic field domain is extended to the scrape-off layer (SOL) via the BMW code, which computes a divergence-free magnetic field solution arising from the external conductors’ vacuum field a nd t he PARVMEC flux surface currents. This magnetic reconstruction is then used in the EMC3-EIRENE transport code in order to compute 3D density profiles. In the last step, the RAPLICASOL code is utilized to compute the ICRF antenna S-matrices resulting from the 3D density profiles. We exemplify this scheme for the ASDEX Upgrade tokamak. A new implementation of a curved model for the ASDEX Upgrade ICRF 2-strap antenna in RAPLICASOL allows simulations in realistic geometry, without any coordinate transformations.

Lopez, Graham↗
NASA NTRS1997

A Self-Consistent Numerical Magnetohydrodynamic (MHD) Model of Helmet Streamer and Flux-Rope Interactions: Initiation and Propagation of Coronal Mass Ejections (CMEs)

We present results for an investigation of the interaction of a helmet streamer arcade and a helical flux-rope emerging from the sub-photosphere. These results are obtained by using a three-dimensional axisymmetric, time-dependent ideal magnetohydrodynamic (MHD) model. Because of the physical nature of the flux-rope, we investigate two types of flux-ropes; (1) high density flux-rope (i.e. flux-rope without cavity), and (2) low density flux rope (i.e. flux-rope with cavity). When the streamer is disrupted by the flux-rope, it will evolve into a configuration resembling the typical observed loop-like Coronal Mass Ejection (CMES) for both cases. The streamer-flux rope system with cavity is easier to be disrupted and the propagation speed of the CME is faster than the streamer-flux rope system without cavity. Our results demonstrate that magnetic buoyancy force plays an important role in disrupting the streamer.

Wu, S. T.↗
NASA NTRS1997

Self-Consistent and Time-Dependent Solar Wind Models

We describe the first results from a self-consistent study of Alfven waves for the time-dependent, single-fluid magnetohydrodynamic (MHD) solar wind equations, using a modified version of the ZEUS MHD code. The wind models we examine are radially symmetrical and magnetized; the initial outflow is described by the standard Parker wind solution. Our study focuses on the effects of Alfven waves on the outflow and is based on solving the full set of the ideal nonlinear MHD equations. In contrast to previous studies, no assumptions regarding wave linearity, wave damping, and wave-flow interaction are made; thus, the models naturally account for the back-reaction of the wind on the waves, as well as for the nonlinear interaction between different types of MHD waves. Our results clearly demonstrate when momentum deposition by Alfven waves in the solar wind can be sufficient to explain the origin of fast streams in solar coronal holes; we discuss the range of wave amplitudes required to obtained such fast stream solutions.

Ong, K. K.↗
DOE OSTI2023

Can we observe the ion-neutral drift velocity in prestellar cores?

Given the low-ionization fraction of molecular clouds, ambipolar diffusion is thought to be an integral process in star formation. However, chemical and radiative-transfer effects, observational challenges, and the fact that the ion-neutral drift velocity is inherently very small render a definite detection of ambipolar diffusion extremely non-trivial. Here, we study the ion-neutral drift velocity in a suite of chemodynamical, non-ideal magnetohydrodynamic (MHD), two-dimensional axisymmetric simulations of prestellar cores where we alter the temperature, cosmic-ray ionization rate, visual extinction, mass-to-flux ratio, and chemical evolution. Subsequently, we perform a number of non-local thermodynamic equilibrium (non-LTE) radiative-transfer calculations considering various idealized and non-idealized scenarios in order to assess which factor (chemistry, radiative transfer, and/or observational difficulties) is the most challenging to overcome in our efforts to detect the ion-neutral drift velocity. We find that temperature has a significant effect in the amplitude of the drift velocity with the coldest modelled cores (T = 6 K) exhibiting drift velocities comparable to the sound speed. Against expectations, we find that in idealized scenarios (where two species are perfectly chemically co-evolving) the drift velocity ‘survives’ radiative-transfer effects and can in principle be observed. However, we find that observational challenges and chemical effects can significantly hinder our view of the ion-neutral drift velocity. Finally, we propose that HCN and HCNH + ⁠, being chemically co-evolving, could be used in future observational studies aiming to measure the ion-neutral drift velocity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Asymptotic quasisymmetric high-beta three-dimensional magnetohydrodynamic equilibria near axisymmetry

Quasisymmetry (QS), a hidden symmetry of the magnetic field strength, is known to support nested flux surfaces and provide superior particle confinement in stellarators. In this work, we study the ideal magnetohydrodynamic (MHD) equilibrium and stability of high-beta plasma in a large-aspect-ratio stellarator. In particular, we show that the lowest-order description of a near-axisymmetric equilibrium vastly simplifies the problem of three-dimensional quasisymmetric MHD equilibria, which can be reduced to a standard elliptic Grad–Shafranov equation for the flux function. We show that any large-aspect-ratio tokamak, deformed periodically in the vertical direction, is a stellarator with approximate volumetric QS. We discuss exact analytical solutions and numerical benchmarks. Finally, we discuss the ideal ballooning and interchange stability of some of our equilibrium configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Determination of a macro- to micro-scale progression leading to a magnetized plasma disruption

We report here the observations of a plasma jet evolving through a macro- to micro-scale progression sequence. This leads to a fast magnetic reconnection that results in the jet breaking off from its originating electrode and forming a force-free state. A sausage-like pinching occurs first and squeezes an initially fat, short magnetized jet so that it becomes thin. The thin jet then becomes kink unstable. The lengthening of the jet by the kinking thins the jet even more since the kink is an incompressible instability. When the jet radius becomes comparable to the ion-skin depth, Hall and electron inertial physics become important and establish the environment for a fast magnetic reconnection. This fast reconnection occurs, disrupting the jet and establishing a force-free state. X-ray bursts and whistler waves, evidence of a magnetic reconnection, are observed when the plasma jet breaks off from the electrode. This experimentally confirmed sequence of successive thinning from pinching followed by kinking is reproduced in a three-dimensional ideal Magnetohydrodynamic (MHD) numerical simulation. The results of the experiment and the numerical simulation, together demonstrate a viable path from macro-scale MHD physics to micro-scale non-MHD physics where fast reconnection occurs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effects of zero and reversed magnetic shear on resistive wall modes in a straight limiter tokamak with zero beta

Advanced tokamak scenarios often feature equilibriums with zero and reversed magnetic shear. To isolate and investigate their impacts on the resistive wall mode (RWM) instability analytically, we construct a series of cylindrical limiter equilibriums with different magnetic shear in the core and toward plasma edge, as a prototype of the configurations in advanced tokamak scenarios. Equilibriums with zero β are considered so that the analysis focuses on the current-driven RWMs. Analytical solutions for the n=1 resistive wall mode, derived from the reduced ideal magnetohydrodynamics (MHD) equations, consistently demonstrate key insights for both the core and the edge: the growth rate of RWMs is inversely proportional to the averaged magnetic shear when the safety factor at the edge (qa) is held constant. Specifically, an increase in positive averaged magnetic shear results in a decreased growth rate, whereas an increase in negative averaged magnetic shear leads to an increased growth rate. Furthermore, the influence of nqa on the mode's stability significantly surpasses the effects of magnetic shear changes previously discussed. These analytical calculation results provide a verification of NIMROD simulations.

Wan, Sui (ORCID:0009000464121677)↗
DOE OSTI2021

Time-resolved turbulent dynamo in a laser plasma

Understanding magnetic-field generation and amplification in turbulent plasma is essential to account for observations of magnetic fields in the universe. A theoretical framework attributing the origin and sustainment of these fields to the so-called fluctuation dynamo was recently validated by experiments on laser facilities in low-magnetic-Prandtl-number plasmas (Pm < 1). However, the same framework proposes that the fluctuation dynamo should operate differently when Pm ≳ 1, the regime relevant to many astrophysical environments such as the intracluster medium of galaxy clusters. This paper reports an experiment that creates a laboratory Pm ≳ 1 plasma dynamo. We provide a time-resolved characterization of the plasma’s evolution, measuring temperatures, densities, flow velocities, and magnetic fields, which allows us to explore various stages of the fluctuation dynamo’s operation on seed magnetic fields generated by the action of the Biermann-battery mechanism during the initial drive-laser target interaction. The magnetic energy in structures with characteristic scales close to the driving scale of the stochastic motions is found to increase by almost three orders of magnitude and saturate dynamically. It is shown that the initial growth of these fields occurs at a much greater rate than the turnover rate of the driving-scale stochastic motions. Here, our results point to the possibility that plasma turbulence produced by strong shear can generate fields more efficiently at the driving scale than anticipated by idealized magnetohydrodynamics (MHD) simulations of the nonhelical fluctuation dynamo; this finding could help explain the large-scale fields inferred from observations of astrophysical systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Structure of pressure-gradient-driven current singularity in ideal magnetohydrodynamic equilibrium

Abstract Singular currents typically appear on rational surfaces in non-axisymmetric ideal magnetohydrodynamic (MHD) equilibria with a continuum of nested flux surfaces and a continuous rotational transition. These currents have two components: a surface current (Dirac δ -function in flux surface labeling) that prevents the formation of magnetic islands, and an algebraically divergent Pfirsch–Schlüter current density when a pressure gradient is present across the rational surface. On flux surfaces adjacent to the rational surface, the traditional treatment gives the Pfirsch–Schlüter current density scaling as J ∼ 1 / Δ ι , where Δ ι is the difference of the rotational transform relative to the rational surface. If the distance s between flux surfaces is proportional to Δ ι , the scaling relation J ∼ 1 / Δ ι ∼ 1 / s will lead to a paradox that the Pfirsch–Schlüter current is not integrable. In this work, we investigate this issue by considering the pressure-gradient-driven singular current in the Hahm–Kulsrud–Taylor problem, which is a prototype for singular currents arising from resonant magnetic perturbations. We show that not only the Pfirsch–Schlüter current density but also the diamagnetic current density are divergent as ∼ 1 / Δ ι . However, due to the formation of a Dirac δ -function current sheet at the rational surface, the neighboring flux surfaces are strongly packed with s ∼ ( Δ ι ) 2 . Consequently, the singular current density J ∼ 1 / s , making the total current finite, thus resolving the paradox. Furthermore, the strong packing of flux surfaces causes a steepening of the pressure gradient near the rational surface, with ∇ p ∼ d p / d s ∼ 1 / s . In general non-axisymmetric MHD equilibrium, contrary to Grad’s conjecture that the pressure profile is flat around densely distributed rational surfaces, our result suggests a pressure profile that densely steepens around them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Kinetic-ballooning-limited pedestals in spherical tokamak plasmas

Abstract A theoretical model is presented that for the first time matches experimental measurements of the pedestal width-height Diallo scaling in the low-aspect-ratio high- β tokamak NSTX. Combining linear gyrokinetics with self-consistent pedestal equilibrium variation, kinetic-ballooning, rather than ideal-ballooning plasma instability, is shown to limit achievable confinement in spherical tokamak pedestals. Simulations are used to find the novel Gyrokinetic Critical Pedestal constraint, which determines the steepest pressure profile a pedestal can sustain subject to gyrokinetic instability. Gyrokinetic width-height scaling expressions for NSTX pedestals with varying density and temperature profiles are obtained. These scalings for STs depart significantly from that of conventional aspect ratio tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Formation and evolution of protostellar accretion discs – I. Angular-momentum budget, gravitational self-regulation, and numerical convergence

Here, we investigate the formation and early evolution of a protostellar disc from a magnetized prestellar core using non-ideal magnetohydrodynamic (MHD) simulations including ambipolar diffusion and Ohmic dissipation. The dynamical contraction of the prestellar core ultimately leads to the formation of a first hydrostatic core, after ambipolar diffusion decouples the magnetic field from the predominantly neutral gas. The hydrostatic core accumulates angular momentum from the infalling material, evolving into a rotationally supported torus; this ‘first hydrostatic torus’ then forms an accreting protostar and a rotationally supported disc. The disc spreads out by gravitational instability, reaching ~30 au in diameter at ~3 kyr after protostar formation. The total mass and angular momentum of the protostar–disc system are determined mainly by accretion of gas from an infalling pseudo-disc, which has low specific angular momentum because of magnetic braking; their removal from the protostar–disc system by outflow and disc magnetic braking are negligible, in part because the magnetic field is poorly coupled there. The redistribution of angular momentum within the protostar–disc system is facilitated mainly by gravitational instability; this allows formation of relatively large discs even when the specific angular momentum of infalling material is low. We argue that such discs should remain marginally unstable as they grow (with Toomre Q ~ 1–2), an idea that is broadly consistent with recent observational estimates for Class 0/I discs. We discuss the numerical convergence of our results, and show that properly treating the inner boundary condition is crucial for achieving convergence at an acceptable computational cost.

79 ASTRONOMY AND ASTROPHYSICS↗