Search NASA⌕ Search

SEARCH · Search NASA

Results for “MHD modeling”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulation of MHD instabilities with fluid runaway electron model in M3D- C 1

Runaway electrons are generated in a tokamak during the start up, during normal operation and during a plasma disruption. During a disruption, runaway electrons can be accelerated to high energies, potentially damaging the first wall. To predict the consequences of runaway generation during a disruption, it is necessary to consider resonant interactions of runaways with the bulk plasma. Here we consider the interactions of runaways on low mode number tearing modes. We have developed a fluid runaway electron model for the 3D MHD code M3D-C1[Jardin,et al. J Comput. Sci Discovery 6 014002 (2012)]. To benchmark, we have reproduced the MHD linear tearing mode results (with runaway electrons) in a circular cylinder presented in previous analytic studies[[Helander, P., et al, Phys. Plasmas 14 144102 (2007)] and extended here with a numerical eigenvalue calculation. Furthermore, we find that the low mode number tearing mode has a rotation caused by the MHD - runaways interaction, and the toroidal current scale length is much smaller with runaways than that for without and decreases as the runaway speed increases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Resistive wall mode stability and resonant field amplification in MAST high beta plasma

The n=1 (n is the toroidal mode number) resistive wall mode (RWM) instability and the resonant field amplification (RFA) due to a stable RWM response are numerically investigated for a MAST high-pressure plasma scenario, utilizing the MARS-F (Liu Y.Q. et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503) codes, for the purpose of understanding the potential RWM behavior in MAST-U plasmas. Ideal kink stability analysis identifies a target plasma, with parameters similar to that of the reference equilibrium reconstructed from a MAST high-pressure discharge, that accesses the RWM regime. The unstable n=1 RWM for the target equilibrium is subject to strong damping by the plasma toroidal flow and/or the drift kinetic effects from thermal particles. As a result, the mode is found to be stable under the experimental flow conditions, or even without flow stabilization if drift kinetic stabilization is included. The stability prediction is robust against variation of the assumed resistive wall minor radius. Active magneto-hydrodynamic (MHD) spectroscopy modeling, using the magnetic coils designed for controlling the edge localized modes (ELMs) in MAST-U as the antenna, shows strong resonant field amplification due to a stable RWM response in the target plasma. Maximal amplification, of 6-7 times larger than the applied vacuum field within the plasma, is obtained assuming a coil phasing of 120-180 degrees between the upper and lower rows of the ELM control coils. The MHD-kinetic hybrid model predicts generally higher RFA than the fluid model, in particular near the low-field and high-field sides of the plasma boundary surface. Furthermore, the MARS-F/K modeling thus shows that the RWM in MAST-U is subject to strong flow and kinetic stabilization, but can nevertheless be detected by active MHD spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Physics basis for the reference flat-top plasma scenario in the ST–E1 fusion power plant

As part of the U.S. Department of Energy’s Milestone-Based Fusion Energy Development Program, Tokamak Energy has completed the pre-concept design of the ST–E1 fusion power plant. ST–E1 is envisaged to operate in two phases: a pilot plant phase, targeting sustained net power production of 300 - 500 MWe for a duration >1 hr, followed by a commercial power plant phase targeting steady-state operations and a normalised overnight capital cost of ⩽12 000 $\$$/kWe. The design process adopted was highly iterative, integrating all major plant systems and progressing in a phased fidelity approach. At the pre-conceptual stage, the emphasis has been on exploring the design space, identifying the main system-level trade-offs, and making the key decisions that define the overall plant concept, rather than optimising a single operating point. This paper, part of a focused collection detailing the ST–E1 pre-concept design, addresses the development of a series of reference flat-top plasma operating points for the pilot plant phase. A modelling workflow was established to develop and assess candidate plasma design points and explore key dependencies. The workflow includes integrated core plasma modelling, magnetohydrodynamic (MHD) stability assessment, equilibrium generation, scrape-off-layer and exhaust modelling, heating & current drive design and optimisation, and turbulent transport modelling. Using this framework, the impact of several key parameters on the flat-top operating space was investigated, including the density limit, core radiation fraction and divertor power loading, level of external heating and curent drive power and assumed pedestal characteristics. The MHD stability, controllability and micro-stability characteristics of these plasmas were also analysed. These investigations informed the definition of a set of fully non-inductive, flat-top reference operating points that satisfy the high-level ST–E1 mission, including a low and high density case, a case that is stable to resistive wall modes and a case with reduced divertor power loading.

ST–E1↗
DOE OSTI2023

Numerical study of magnetic island coalescence using magnetohydrodynamics with adaptively embedded particle-in-cell model

Collisionless magnetic reconnection typically requires kinetic treatment that is, in general, computationally expensive compared to fluid-based models. In this study, we use the magnetohydrodynamics with an adaptively embedded particle-in-cell (MHD-AEPIC) model to study the interaction of two magnetic flux ropes. This innovative model embeds one or more adaptive PIC regions into a global MHD simulation domain such that the kinetic treatment is only applied in regions where the kinetic physics is prominent. We compare the simulation results among three cases: (1) MHD with adaptively embedded PIC regions, (2) MHD with statically (or fixed) embedded PIC regions, and (3) a full PIC simulation. The comparison yields good agreement when analyzing their reconnection rates and magnetic island separations as well as the ion pressure tensor elements and ion agyrotropy. In order to reach good agreement among the three cases, large adaptive PIC regions are needed within the MHD domain, which indicates that the magnetic island coalescence problem is highly kinetic in nature, where the coupling between the macro-scale MHD and micro-scale kinetic physics is important.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Influence of initial conditions on data-driven model identification and information entropy for ideal mhd problems

Data-driven methods of model identification are able to discern governing dynamics of a system from data. Such methods are well suited to help us learn about systems with unpredictable evolution or systems with ambiguous governing dynamics given our current understanding. Many plasma problems of interest fall into these categories as there are a wide range of models that exist, however each model is only useful in a certain regime and often limited by computational complexity. To ensure data-driven methods align with theory, they must be consistent and predictable when acting on data whose governing dynamics are known. Weak Sparse Identification of Nonlinear Dynamics (WSINDy) is a recently developed data-driven method that has shown promise in learning governing dynamics from data with high noise levels [1]. This work examines how WSINDy acts on ideal MHD test problems as the initial conditions are varied and specifies limiting requirements for successful equation identification. Furthermore, it is hard to recover the governing dynamics from data that emphasize a single dominant behavior. In these low information cases, Shannon information entropy is able to pick up on the redundancies in the data that affect recoverability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Multiscale and multidimensional modeling of particle acceleration and transport in solar flares

Multi-messenger, multi-viewpoint, and time-resolved observations of solar flares are now providing unprecedented constraints on particle acceleration sites, energy conversion, and energy transport. The interpretation of current observations, including microwave imaging spectroscopy from EOVSA, hard x-ray (HXR) imaging from Solar Orbiter/STIX, gamma-ray diagnostics from Fermi, and in situ measurements from Parker Solar Probe and Solar Orbiter, collectively demands modeling frameworks that go beyond traditional spatially unresolved, one-zone models or single-mechanism descriptions. This review surveys multiscale and multidimensional modeling approaches, including kinetic, magnetohydrodynamic (MHD), and macroscopic particle models, that are being developed to meet the need. Kinetic simulations reveal that three-dimensional (3D) effects, including field-line chaos and self-generated turbulence, are essential for sustained power-law particle acceleration. MHD simulations now capture flux-rope eruptions, plasmoid-unstable current sheets, and turbulent flare regions in realistic magnetic topologies. Macroscopic models coupling MHD with energetic-particle models produce spatially resolved electron distributions and synthetic HXR and microwave emissions for direct comparison with observations. Despite these advances, outstanding challenges remain in bridging kinetic and global scales, improving MHD simulations and macroscopic particle models, and achieving quantitative model-observation closure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2025

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics↗
DOE OSTI2023

Magnetohydrodynamics with physics informed neural operators

Abstract The modeling of multi-scale and multi-physics complex systems typically involves the use of scientific software that can optimally leverage extreme scale computing. Despite major developments in recent years, these simulations continue to be computationally intensive and time consuming. Here we explore the use of AI to accelerate the modeling of complex systems at a fraction of the computational cost of classical methods, and present the first application of physics informed neural operators (NOs) (PINOs) to model 2D incompressible magnetohydrodynamics (MHD) simulations. Our AI models incorporate tensor Fourier NOs as their backbone, which we implemented with the TensorLY package. Our results indicate that PINOs can accurately capture the physics of MHD simulations that describe laminar flows with Reynolds numbers Re ⩽ 250 . We also explore the applicability of our AI surrogates for turbulent flows, and discuss a variety of methodologies that may be incorporated in future work to create AI models that provide a computationally efficient and high fidelity description of MHD simulations for a broad range of Reynolds numbers. The scientific software developed in this project is released with this manuscript.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Elimination of MHD current sheets by modifications to the plasma wall in a fixed boundary model

Models of magnetohydrodynamic (MHD) equilibia that for computational convenience assume the existence of a system of nested magnetic flux surfaces tend to exhibit singular current sheets. These sheets are located on resonant flux surfaces that are associated with rational values of the rotational transform. We study the possibility of eliminating these singularities by suitable modifications of the plasma boundary, which we prescribe in a fixed boundary setting. We find that relatively straightforward iterative procedures can be used to eliminate weak current sheets that are generated at resonant flux surfaces by the nonlinear interactions of resonating wall harmonics. These types of procedures may prove useful in the design of fusion devices with configurations that enjoy improved stability and transport properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Physics-constrained, low-dimensional models for magnetohydrodynamics: First-principles and data-driven approaches

Plasmas are highly nonlinear and multiscale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD), although these reduced models hold promise for understanding key physical mechanisms, efficient computation, and real-time optimization and control. Galerkin models, obtained by projection of the MHD equations onto a truncated modal basis, and data-driven models, obtained by modern machine learning and system identification, can furnish this gap in the lower levels of the model hierarchy. This work develops a reduced-order modeling framework for compressible plasmas, leveraging decades of progress in projection-based and data-driven modeling of fluids. We begin by formalizing projection-based model reduction for nonlinear MHD systems. To avoid separate modal decompositions for the magnetic, velocity, and pressure fields, we introduce an energy inner product to synthesize all of the fields into a dimensionally consistent, reduced-order basis. Next, we obtain an analytic model by Galerkin projection of the Hall-MHD equations onto these modes. We illustrate how global conservation laws constrain the model parameters, revealing symmetries that can be enforced in data-driven models, directly connecting these models to the underlying physics. We demonstrate the effectiveness of this approach on data from high-fidelity numerical simulations of a three-dimensional spheromak experiment. Finally, this manuscript builds a bridge to the extensive Galerkin literature in fluid mechanics and facilitates future principled development of projection-based and data-driven models for plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Impact of electron transport models on capillary discharge plasmas

Magnetohydrodynamics (MHD) can be used to model capillary discharge waveguides in laser-wakefield accelerators. However, the predictive capability of MHD can suffer due to poor microscopic closure models. Here, we study the impact of electron heating and thermal conduction on the capillary waveguide performance as part of an effort to understand and quantify uncertainties in modeling and designing next-generation plasma accelerators. To do so, we perform two-dimensional high-resolution MHD simulations using an argon-filled capillary discharge waveguide with three different electron transport coefficients models. The models tested include (i) Davies et al., (ii) Spitzer, and (iii) Epperlein–Haines (EH). We found that the EH model overestimates the electron temperature inside the channel by over 20% while predicting a lower azimuthal magnetic field. Furthermore, the Spitzer model, often used in MHD simulations for plasma-based accelerators, predicts a significantly higher electron temperature than the other models suggest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling Electron Acceleration and Transport in the Early Impulsive Phase of the 2017 September 10th Solar Flare

The X8.2-class limb flare on 2017 September 10 is among the best studied solar flare events owing to its great similarity to the standard flare model and the broad coverage by multiple spacecraft and ground-based observations. These multiwavelength observations indicate that electron acceleration and transport are efficient in the reconnection and flare looptop regions. However, there lacks a comprehensive model for explaining and interpreting the multi-faceted observations. In this work, we model the electron acceleration and transport in the early impulsive phase of this flare. We solve the Parker transport equation that includes the primary acceleration mechanism during magnetic reconnection in the large-scale flare region modeled by MHD simulations. We find that electrons are accelerated up to several MeV and fill a large volume of the reconnection region, similar to the observations shown in microwaves. The electron spatial distribution and spectral shape in the looptop region agree well with those derived from the microwave and hard X-ray emissions before magnetic islands grow large and dominate the acceleration. Future emission modelings using the electron maps will enable direct comparison with microwave and hard X-ray observations. These results shed new light on the electron acceleration and transport in a broad region of solar flares within a data-constrained realistic flare geometry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Impact of equilibrium radial electric field on energy loss process after pedestal collapse

An impact of the equilibrium radial electric field on energy loss processes after pedestal collapse is numerically investigated using the BOUT++ framework. Using linear stability analysis, the resistive ballooning mode is shown to be stabilized by the radial shear of the equilibrium radial electric field. On the other hand, the energy loss level after the pedestal collapse increases if the equilibrium radial electric field is taken into account. The spatio-temporal and phase diagram analyses show that the equilibrium radial electric field partially cancels the fluctuation-driven toroidally axisymmetric radial electric field and weakens the E × B shearing rate after pedestal collapse, weakening the turbulence suppression by vortex shearing. Finally, the equilibrium radial electric field therefore increases turbulence intensity in nonlinear cyclic oscillations among pressure gradient, E × B shearing rate, and turbulence intensity, which gives rise to subsequent bursts of turbulent transport and increases the energy loss level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Early Solar Wind and Dynamo Magnetic Field Topology Predictions for (16) Psyche and Other Asteroids

Abstract Asteroid (16) Psyche is a metal‐rich body that might record an ancient coherent magnetization if some relict crust or mantle is preserved. Herein, we use magnetohydrodynamic simulations to predict (16) Psyche's field topology for several distinct pathways: (i) an early solar wind‐induced magnetization imparted after a larger body was impacted, forming the present‐day asteroid, (ii) a core dynamo magnetization imparted in an asteroid that is either presently largely intact or was a rubble pile, and (iii) magnetization in the turbulent solar nebula disk. For pathway (i) we find the field to be predominantly dipolar and spin axis‐aligned. For pathway (ii) we find the field to be either dipolar and spin axis‐misaligned, or highly multipolar. We also find that (iii) a field produced earlier before the solar nebula cleared, would be highly multipolar. In cases (i) and (ii) we also place constraints on the field strength. Simple detection of a magnetic field without constraining its topology and temporal variability would be insufficient to confirm a remanent source, due to the influence of the present‐day solar wind, electromagnetic induction, and (16) Psyche's high obliquity. For sufficiently strong fields however, the field topology and orientation may reveal key observable consequences of the nature and history of (16) Psyche. Our framework is also broadly applicable to the study of magnetic fields from other asteroids.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Resistive contributions to the stability of DIII-D ITER baseline demonstration discharges

Here, measurements of the plasma response to applied n = 1 perturbations from DIII-D ITER baseline scenario (IBS) demonstration discharges exhibit increasing amplitude prior to disruption inducing mode locking events. Simulations reveal that changes in ideal and resistive stability impact the response in the DIII-D IBS regime, well below the pressure limit of the external kink mode. The dependencies of the response measurements on the plasma normalized internal inductance ℓ i and beta β N are qualitatively consistent with ideal MHD, although in most cases the amplitude of the measurements exceeds predictions, indicating that the experimental discharges are less stable than expected. Resistive MHD simulations of the response show improved compatibility with the measurements, and the closest agreement is obtained by including the experimental plasma rotation in the simulations. Although the input neutral beam torque is near zero, the simulations show that the corresponding level of rotation can lead to nearly complete screening of the pitch-resonant field component at the q = 2 surface. However, the simulation results exhibit a significant weakening of this screening effect at the lowest rotation levels in the data set, concurrent with predictions of elevated response amplitudes that are consistent with measurements. These comparisons with DIII-D data provide a qualitative validation of a linear, resistive MHD response model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Two-stage crash process in resistive drift ballooning mode driven ELM crash

We report a two-stage crash process in edge localized mode (ELM) driven by resistive drift-ballooning modes (RDBMs) numerically simulated in a full annular torus domain with a scale-separated four-field reduced MHD (RMHD) model using the BOUT++ framework. In the early nonlinear phase, the small first crash is triggered by linearly unstable RDBMs, and m/n=2/1 magnetic islands are nonlinearly excited by nonlinear coupling of RDBMs as well as their higher harmonics. Here, m is the poloidal mode number, n is the toroidal mode number, the q = 2 rational surface exists near the pressure gradient peak, and q is the safety factor. Simultaneously, middle-n RDBM turbulence develops but is poloidally localized around X-points of the magnetic islands, leading to the small energy loss. The second large crash occurs in the late nonlinear phase. Higher harmonics of m/n=2/1 magnetic islands well develop around the q = 2 surface via nonlinear coupling and make the magnetic field stochastic by magnetic island overlapping. Turbulence heat transport develops at X-points of higher harmonics of m/n=2/1 magnetic islands, resulting in the turbulence spreading in the poloidal direction. The large second crash is triggered when the turbulence covers the whole poloidal region so that the magnetic island generation and magnetic field stochastization before the large crash can be interpreted as ELM precursors. It is concluded that the ELM trigger is attributed to the turbulent spreading in the poloidal direction in synchronization with the magnetic field stochastization and the crash is driven by E × B convection rather than the conventional Rechester–Rosenbluth anomalous electron heat transport.

Seto, H. (ORCID:0000000259278643)↗
DOE OSTI2021

Stabilization of kink/peeling modes by coupled rotation and ion diamagnetic drift effects in quiescent H-mode plasmas in DIII-D and JT-60U

Magnetohydrodynamic stability at the edge pedestal in several quiescent H-mode (QH-mode) plasmas in DIII-D and JT-60U experiments was analyzed by considering plasma rotation and ion diamagnetic drift effects. It was identified that a kink/peeling mode, which is a prime candidate for a trigger of edge harmonic oscillation in QH-mode, is stabilized by plasma rotation when considering the ion diamagnetic drift simultaneously in both experiments. Furthermore, the stabilizing effect by rotation becomes more effective in case using the rotation profile of the main ion species evaluated by assuming radial force balance. In addition, when inverting the rotation direction, it was found that the kink/peeling mode is more stabilized when considering the rotation of the main ion species, though the mode is less stabilized by the rotation of impurity ion species. The result implies that the kink/peeling mode stability in QH-mode plasmas is sensitive to how the rotation profile is evaluated, but it is shown that a qualitative trend stabilizing the kink/peeling mode by rotation can be reproduced with the poloidal rotation profile of an impurity predicted numerically based on the neoclassical theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗