Search NASA⌕ Search

SEARCH · Search NASA

Results for “MAGNETOHYDRODYNAMICS”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

How a Realistic Magnetosphere Alters the Polarizations of Surface, Fast Magnetosonic, and Alfvén Waves

System-scale magnetohydrodynamic (MHD) waves within Earth's magnetosphere are often understood theoretically using box models. While these have been highly instructive in understanding many fundamental features of the various wave modes present, they neglect the complexities of geospace such as the inhomogeneities and curvilinear geometries present. Here, we show global MHD simulations of resonant waves impulsively excited by a solar wind pressure pulse. Although many aspects of the surface, fast magnetosonic (cavity/waveguide), and Alfvén modes present agree with the box and axially symmetric dipole models, we find some predictions for large-scale waves are significantly altered in a realistic magnetosphere. The radial ordering of fast mode turning points and Alfvén resonant locations may be reversed even with monotonic wave speeds. Additional nodes along field lines that are not present in the displacement/velocity occur in both the perpendicular and compressional components of the magnetic field. Close to the magnetopause, the perpendicular oscillations of the magnetic field have the opposite handedness to the velocity. Finally, widely used detection techniques for standing waves, both across and along the field, can fail to identify their presence. We explain how all these features arise from the MHD equations when accounting for a non-uniform background field and propose modified methods that might be applied to spacecraft observations.

ULF Waves↗
NASA NTRS2022

Variability of the Reconnection Guide Field in Solar Flares

Solar flares may be the best-known examples of the explosive conversion of magnetic energy into bulk motion, plasma heating, and particle acceleration via magnetic re-connection. The energy source for all flares is the highly sheared magnetic field of a filament channel above a polarity inversion line (PIL). During the flare, this shear field becomes the so-called reconnection guide field (i.e., the non-reconnecting component), which has been shown to play a major role in determining key properties of the recon-nection including the efficiency of particle acceleration. We present new high-resolution, three-dimensional, magnetohydrodynamics simulations that reveal the detailed evolution of the magnetic shear/guide field throughout an eruptive flare. The magnetic shear evolves in three distinct phases: shear first builds up in a narrow region about the PIL, then expands outward to form a thin vertical current sheet, and finally is transferred by flare reconnection into an arcade of sheared flare loops and an erupting flux rope. We demonstrate how the guide field may be inferred from observations of the sheared flare loops. Our results indicate that initially the guide field is larger by about a factor of 5 than the reconnecting component, but it weakens by more than an order of magnitude over the course of the flare. Instantaneously, the guide field also varies spatially over a similar range along the three-dimensional current sheet. We discuss the implications of the remarkable variability of the guide field for the timing and localization of efficient particle acceleration in flares.

Solar magnetic reconnection↗
NASA NTRS2023

Understanding the Deflection of the “Cartwheel CME”: Data Analysis and Modeling

We study the low corona evolution of the “Cartwheel” coronal mass ejection (CME; 2008 April 9) by- reconstructing its three-dimensional path and modeling it with magnetohydrodynamic simulations. This event exhibited a double deflection that has been reported and analyzed in previous works but whose underlying cause remained unclear. The Cartwheel CME traveled toward a coronal hole (CH) and against the magnetic gradients. Using a high-cadence, full-trajectory reconstruction, we accurately determine the location of the magnetic flux rope (MFR) and, consequently, the magnetic environment in which it is immersed. We find a pseudostreamer (PS) structure whose null point may be responsible for the complex evolution of the MFR at the initial phase. From the preeruptive magnetic field reconstruction, we estimate the dynamic forces acting on the MFR and provide a new physical insight into the motion exhibited by the 2008 April 9 event. By setting up a similar magnetic configuration in a 2.5D numerical simulation we are able to reproduce the observed behavior, confirming the importance of the PS null point. We find that the magnetic forces directed toward the null point cause the first deflection, directing the MFR toward the CH. Later, the magnetic pressure gradient of the CH produces the reversal motion of the MFR.

Solar prominences↗
NASA NTRS2023

Investigating Potential Causes for the Prediction of Spurious Magnetopause Crossings at Geosynchronous Orbit in MHD Simulations

During intense geomagnetic storms, the magnetopause can move in as far as geosynchronous orbit, leaving the satellites in that orbit out in the magnetosheath. Spacecraft operators turn to numerical models to predict the response of the magnetopause to solar wind conditions, but the predictions of the models are not always accurate. This study investigates four storms with a magnetopause crossing by at least one GOES satellite, using four magnetohydrodynamic models at NASA's Community Coordinated Modeling Center to simulate the events, and analyzes the results to investigate the reasons for errors in the predictions. Two main reasons can explain most of the erroneous predictions. First, the solar wind input to the simulations often contains features measured near the L1 point that did not eventually arrive at Earth; incorrect predictions during such periods are due to the solar wind input rather than to the models themselves. Second, while the models do well when the primary driver of magnetopause motion is a variation in the solar wind density, they tend to overpredict or underpredict the integrated Birkeland currents and their effects during times of strong negative interplanetary magnetic field (IMF) Bz, leading to poorer prediction capability. Coupling the MHD codes to a ring current model, when such a coupling is available, generally will improve the predictions but will not always entirely correct them. More work is needed to fully characterize the response of each code under strong southward IMF conditions as it relates to prediction of magnetopause location.

magnetopause↗
NASA NTRS2023

Comparing Magnetopause Predictions From Two MHD Models During A Geomagnetic Storm and A Quiet Period

Magnetopause location is an important prediction of numerical simulations of the magnetosphere, yet the models can err, either under-predicting or over-predicting the motion of the boundary. This study compares results from two of the most widely used magnetohydrodynamic (MHD) models, the Lyon–Fedder–Mobarry (LFM) model and the Space Weather Modeling Framework (SWMF), to data from the GOES 13 and 15 satellites during the geomagnetic storm on 22 June 2015, and to THEMIS A, D, and E during a quiet period on 31 January 2013. The models not only reproduce the magnetopause crossings of the spacecraft during the storm, but they also predict spurious magnetopause motion after the crossings seen in the GOES data. We investigate the possible causes of the over-predictions during the storm and find the following. First, using different ionospheric conductance models does not significantly alter predictions of the magnetopause location. Second, coupling the Rice Convection Model (RCM) to the MHD codes improves the SWMF magnetopause predictions more than it does for the LFM predictions. Third, the SWMF produces a stronger ring current than LFM, both with and without the RCM and regardless of the LFM spatial resolution. During the non-storm event, LFM predicts the THEMIS magnetopause crossings due to the southward interplanetary magnetic field better than the SWMF. Additionally, increasing the LFM spatial grid resolution improves the THEMIS predictions, while increasing the SWMF grid resolutions does not.

magnetohydrodynamics↗
DOE OSTI2024

HL-2A's ELM cycle simulations by integrating BOUT++'s drift MHD and transport code

A new integrating model has been developed to couple tokamak edge multiscale magnetohydrodynamic (MHD) events and transport simulations, such as edge-localized mode (ELM) cycles. As a proof of principle, we first start from a set of three-field two-fluid model equations, which includes the pressure, current, and vorticity. Here, the equations are separated into the slowly evolving part of the axisymmetric component by taking a time average of the axisymmetric component. The time-averaged fluxes, which are quadratic in fluctuating quantities, act as driven terms for the time-averaged axisymmetric quantities that determine the plasma transport, and therefore the large-scale evolution of the plasma profiles. Then the HL-2A's ELM cycles are simulated using the model. Good agreements of ELM size and pedestal recovery time have been achieved for the solutions obtained from the coupled simulation compared with experiment. For one ELM cycle simulation, the coupled code can achieve a speedup of a factor of up to 30 over standalone code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Magnetic Reconnection and Associated Particle Acceleration in High-Energy Astrophysics

Abstract Magnetic reconnection occurs ubiquitously in the universe and is often invoked to explain fast energy release and particle acceleration in high-energy astrophysics. The study of relativistic magnetic reconnection in the magnetically dominated regime has surged over the past two decades, revealing the physics of fast magnetic reconnection and nonthermal particle acceleration. Here we review these recent progresses, including the magnetohydrodynamic and collisionless reconnection dynamics as well as particle energization. The insights in astrophysical reconnection strongly connect to the development of magnetic reconnection in other areas, and further communication is greatly desired. We also provide a summary and discussion of key physics processes and frontier problems, toward a better understanding of the roles of magnetic reconnection in high-energy astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Ohm’s Law, the Reconnection Rate, and Energy Conversion in Collisionless Magnetic Reconnection

Magnetic reconnection is a ubiquitous plasma process that transforms magnetic energy into particle energy during eruptive events throughout the universe. Reconnection not only converts energy during solar flares and geomagnetic substorms that drive space weather near Earth, but it may also play critical roles in the high energy emissions from the magnetospheres of neutron stars and black holes. In this review article, we focus on collisionless plasmas that are most relevant to reconnection in many space and astrophysical plasmas. Guided by first-principles kinetic simulations and spaceborne in-situ observations, we highlight the most recent progress in understanding this fundamental plasma process. We start by discussing the non-ideal electric field in the generalized Ohm’s law that breaks the frozen-in flux condition in ideal magnetohydrodynamics and allows magnetic reconnection to occur. We point out that this same reconnection electric field also plays an important role in sustaining the current and pressure in the current sheet and then discuss the determination of its magnitude (i.e., the reconnection rate), based on force balance and energy conservation. This approach to determining the reconnection rate is applied to kinetic current sheets with a wide variety of magnetic geometries, parameters, and background conditions. We also briefly review the key diagnostics and modeling of energy conversion around the reconnection diffusion region, seeking insights from recently developed theories. Finally, future prospects and open questions are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Simulation Models for Exploring Magnetic Reconnection

Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This article reviews many of these methods, laying out critical assumptions, numerical techniques, and giving examples of scientific results. Plasma models described include magnetohydrodynamics (MHD), Hall MHD, Hybrid, kinetic particle-in-cell (PIC), kinetic Vlasov, Fluid models with embedded PIC, Fluid models with direct feedback from energetic populations, and the Rice Convection Model (RCM).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Three-step LES-C models for flows at high Reynolds numbers

Here, we investigate the need for the second correction step in the recently proposed LES-C models for fluid flows at high Reynolds numbers. These models use a predictor-corrector idea to enhance the efficiency of the existing Large Eddy Simulation models. Different three-step(one defect step, two corrections) LES-C models, based on the Leray-α, ADM and NS-ω LES models, are tested in three different situations. The new Leray-α-C2 model (C2 stands for two correction steps) is applied to the Navier–Stokes equations; the ADC2 is applied to the MagnetoHydroDynamic flow; and the NS-ω-C2 is used in the fluid-fluid interaction problem. We evaluate the effectiveness of the second correction step in all these settings, using qualitative and quantitative numerical tests.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A review of high order strong stability preserving two-derivative explicit, implicit, and IMEX methods

High order strong stability preserving time discretizations ensure the nonlinear non-inner-product strong stability properties of spatial discretizations suited for the stable simulation of hyperbolic PDEs in a wide variety of application areas including fluid dynamics, magnetohydrodynamics, semiconductor devices, electromagnetics, and astrophysics. Over the past decade multiderivative time-stepping have been increasingly used for the time-evolution hyperbolic PDEs, so that the strong stability properties of these methods have become important. In this work we review sufficient conditions for a two-derivative multistage method to preserve the strong stability properties of spatial discretizations in a forward Euler and different conditions on the second derivative. In particular we present the strong stability preserving theory for explicit and implicit two-derivative Runge–Kutta schemes, including a special condition on the second derivative under which these implicit methods may be unconditionally strong stability preserving. This special condition is natural for the stiff component of wide range of plasma physics problems, and can be useful in the context of strong stability preserving implicit-explicit multi-derivative Runge–Kutta schemes, where the time-step restriction is then independent of the stiff term. Lastly, we present the strong stability preserving theory for implicit-explicit multi-derivative general linear methods, and some novel second and third order methods where the time-step restriction is independent of the stiff term.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following Du Toit et al. (2018), the scheme transforms mixed-derivative diffusion fluxes (which are responsible for the lack of a discrete maximum principle) into nonlinear advective fluxes, amenable to nonlinear-solver-friendly monotonicity-preserving limiters. The scheme enables accurate multi-dimensional heat transport simulations with up to seven orders of magnitude of heat-transport-coefficient anisotropies with low cross-field numerical error pollution and excellent algorithmic performance, with the number of linear iterations scaling very weakly with grid resolution and grid anisotropy, and scaling with the square-root of the implicit timestep. We propose a multigrid preconditioning strategy based on a lower-order approximation that renders the scheme efficient and scalable under grid refinement. Several numerical tests are presented that display the expected spatial convergence rates and strong algorithmic performance, including fully nonlinear magnetohydrodynamics simulations of kink instabilities in a Bennett pinch in 2D helical geometry and of ITER in 3D toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2026

Fast solvers for tokamak fluid models with PETSc

Multigrid (MG) is widely recognized as a highly effective solver for the model problem, the Laplacian, but textbook MG fails on most problems of interest. MG methods have been applied to complex, real-world applications with careful consideration of the physical model and discretization. In this work we develop the first step in applying MG methods to science and engineering relevant magnetohydrodynamics (MHD) tokamak models in the M3D-C1 (https://m3dc1.pppl.gov) fusion energy science code. The semi-implicit time integrator in M3D-C1 is composed of many linear solves. The implicit advance of the momentum equation is the most challenging and is the focus of this work. The current production solver in M3D-C1 is a block Jacobi (BJ) preconditioner within a Krylov solver, where blocks group degrees of freedom on planes of constant toroidal coordinate. BJ convergence degrades as the number of planes increases due to the spectral properties of the matrix preconditioned with BJ. The partially magnetic field-aligned, regular toroidal grid structure in M3D-C1 is amenable to semi-coarsening geometric MG in the toroidal direction. This paper develops such a solver and demonstrates competitive performance on a runaway electron model of a SPARC (https://cfs.energy/technology/sparc) disruption, and superior robustness on a stellarator model on which the BJ solver fails to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nonlinear control of the minimum safety factor in tokamaks by optimal allocation of spatially moving electron cyclotron current drive

The minimum value of the safety factor profile is related to the magnetohydrodynamic (MHD) stability of the plasma confined in a tokamak. Therefore, active control of the minimum safety factor may mitigate MHD instabilities that can degrade or even terminate plasma confinement. Typically, in most tokamak scenarios, the minimum safety factor evolves spatially with time, i.e., the location at which the safety factor achieves the minimum value changes with time. In addition to the inherent nonlinearities in the minimum safety factor evolution, its spatial variation makes the control design challenging. In particular, complexity in control design may arise from the need for time-dependent nonlinear models that account for spatial variation of the minimum safety factor. Furthermore, the minimum safety factor may drift to locations where the actuator authority is low. The problem of minimum safety factor control with target location tracking and moving electron cyclotron current drive (ECCD) is addressed in this work. A nonlinear time-dependent model that incorporates the spatial variation of the minimum safety factor is presented. A nonlinear controller based on optimal feedback linearization is developed to track a target minimum safety factor. The proposed controller treats the ECCD position as a controllable variable. In other words, the controller prescribes the ECCD position (in addition to the non-inductive powers) in real time based on an optimal criterion that is defined a priori. This work also presents the steps necessary to integrate the minimum safety factor controller with a total energy controller to achieve multiple control objectives simultaneously. In conclusion, the proposed integrated control algorithm is tested using nonlinear simulations in the Control Oriented Transport SIMulator (COTSIM) for a DIII-D tokamak scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Liquid lithium divertor analysis using coupled plasma material interaction model

A liquid lithium divertor can improve performance of future fusion devices by creating efficient power exhaust and improving the energy confinement via pumping of the hydrogen isotopes. In addition, significantly higher heat fluxes can be handled if controlled vapor shielding is used to redistribute the divertor heat flux over a wider area. Design and optimization of such a system calls for an analysis model which includes a strong two-way coupling between the plasma and divertor material. The incoming plasma heat and particle flux will affect the divertor surface temperature, which is a defining factor of the lithium evaporative and sputtered flux going into the plasma. Results of the coupled model based on the plasma edge code SOLPS-ITER and the computational fluid dynamics (CFD) code ANSYS-CFX will be presented for different configurations. An analytical slab flow model is used as a heat transfer boundary condition for SOLPS, defining particle flux from the wall via calculation of the surface temperature. At the final step, results of the SOLPS analysis are verified using a 3D CFD magnetohydrodynamics (MHD) analysis which uses heat and particle flux from SOLPS as a boundary condition. In addition to plasma heat flux, both analytical and CFD temperature models include several plasma material interaction effects, such as lithium evaporation, condensation and sputtering based on deuterium target flux. New adatom sputtering model based on the available experimental data is presented. Analytical model is expanded to include free surface axisymmetric configurations. Results of parametric studies of the divertor configurations with different lithium inlet temperature and velocity will be presented leading to the optimal design resulting in the lowest possible lithium contamination in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Magneto-Stokes flow in a shallow free-surface annulus

In this study, we analyse ‘magneto-Stokes’ flow, a fundamental magnetohydrodynamic (MHD) flow that shares the cylindrical-annular geometry of the Taylor–Couette cell but uses applied electromagnetic forces to circulate a free-surface layer of electrolyte at low Reynolds numbers. The first complete, analytical solution for time-dependent magneto-Stokes flow is presented and validated with coupled laboratory and numerical experiments. Three regimes are distinguished (shallow-layer, transitional and deep-layer flow regimes), and their influence on the efficiency of microscale mixing is clarified. The solution in the shallow-layer limit belongs to a newly identified class of MHD potential flows, and thus induces mixing without the aid of axial vorticity. We show that these shallow-layer magneto-Stokes flows can still augment mixing in distinct Taylor dispersion and advection-dominated mixing regimes. The existence of enhanced mixing across all three distinguished flow regimes is predicted by asymptotic scaling laws and supported by three-dimensional numerical simulations. Mixing enhancement is initiated with the least electromagnetic forcing in channels with order-unity depth-to-gap-width ratios. If the strength of the electromagnetic forcing is not a constraint, then shallow-layer flows can still yield the shortest mixing times in the advection-dominated limit. Our robust description of momentum evolution and mixing of passive tracers makes the annular magneto-Stokes system fit for use as an MHD reference flow.

58 GEOSCIENCES↗
DOE OSTI2024

Simulations of radiatively cooled magnetic reconnection driven by pulsed power

Magnetic reconnection is an important process in astrophysical environments, as it reconfigures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the reconnection layer. In this paper, we perform two- and three-dimensional simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively cooled regime. These simulations are performed in GORGON, an Eulerian two-temperature resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and P 1/3 multi-group radiation transport. The load produces highly collisional, super-Alfvénic (Alfvén Mach number M A ≈ 1.5), supersonic (Sonic Mach number M S ≈ 4–5) strongly driven plasma flows which generate an elongated reconnection layer (Aspect Ratio L/δ≈100 , Lundquist number S L ≈ 400). The reconnection layer undergoes radiative collapse when the radiative losses exceed the rates of ohmic and compressional heating (cooling rate/hydrodynamic transit rate = τ$^{–1}_{cool}$/τ$^{–1}_{H}$ ≈100); this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Lastly, the current sheet is also unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗