Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “resistive 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Influence of plasma sheet thickening toward the tail flanks on magnetotail stability and dynamics

On the basis of resistive MHD simulations we investigate the stability and dynamic evolution of the magnetotail for configurations with various degrees of plasma sheet thickening from midnight toward the tail flanks, associated with a similar increase of the B(sub z) component. This increase is varied by factors between 2 and 3.5, while the magnitude of B(sub z) at midnight is left unchanged. The increase of B(sub z) has a strong effect on the spatial extend of resistive tearing and reconnection across the tail and on the unstable evolution. For stronger thickening the cross-tail extent of the near-Earth neutral line, formed by reconnection, and correspondingly the amount of reconnected magnetic flux get reduced and the dipolarization effects earthward of the neutral line are more concentrated near midnight. An increase of B(sub z) from midnight toward the flanks by more than a factor of about 3 can possibly even suppress the resistive tearing instability. This indicates the possibility of a stability transition due to a reduction of B(sub z) in the flank regions in response to changes in the solar wind, even when the magnitude of B(sub z) near midnight and the fluctuation level in that region are unchanged (provided that a sufficient fluctuation level exists).

Birn, Joachim↗
NASA NTRS2016

Separator Reconnection at the Magnetopause for Predominantly Northward and Southward IMF: Techniques and Results

In this work, we demonstrate how to track magnetic separators in three-dimensional simulated magnetic fields with or without magnetic nulls, apply these techniques to enhance our understanding of reconnection at the magnetopause. We present three methods for locating magnetic separators and apply them to 3-D resistive MHD simulations of the Earth's magnetosphere using the Block-Adaptive-Tree Solar-wind Roe-type Upwind Scheme code. The techniques for finding separators and determining the reconnection rate are insensitive to interplanetary magnetic field (IMF) clock angle and can in principle be applied to any magnetospheric model. Moreover, the techniques have a number of advantages over prior separator finding techniques applied to the magnetosphere. The present work examines cases of high and low resistivity for two clock angles. We go beyond previous work examine the separator during Flux Transfer Events (FTEs). Our analysis of reconnection on the magnetopause yields a number of interesting conclusions: Reconnection occurs all along the separator even during predominately northward IMF cases. Multiple separators form in low-resistivity conditions, and in the region of an FTE the separator splits into distinct branches. Moreover, the local contribution to the reconnection rate, as determined by the local parallel electric field, drops in the vicinity of the FTE with respect to the value when there are none.

modeling↗
DOE OSTI2025

Exploring the transition from continuous turbulence fluctuations to bursting ELMs in high SOL density regimes

BOUT++ turbulence simulations of the DIII-D reveal that the density profile between the separatrix and pedestal plays a crucial role in the dynamics of edge localized modes (ELMs) and edge plasma turbulent transport. Nonlinear simulations demonstrate that small ELMs in the DIII-D hybrid scenario under high SOL density conditions are predominantly driven by local ballooning modes near the separatrix, stabilizing global instabilities while enhancing localized pressure fluctuations. A key control parameters for ELM dynamics is the separatrix-to-pedestal density ratio, n e,sep /n e,ped . A high ratio indicates a shallow gradient, favoring small ELMs, while a lower ratio signals a steep gradient, which increases the likelihood of large ELMs. Comprehensive parameter scans, including n e,sep /n e,ped , density gradient profiles near the separatrix, and resistivity, reveal the critical role of these parameters in shaping transitions between turbulence-driven transport and ELM bursting. The scans demonstrate that in high SOL density regimes, small ELMs can result from either global resistive MHD instabilities or local ballooning modes near the separatrix, depending on the steepness of the separatrix density gradient. These findings also highlight the transition from continuous turbulence to small ELMs. The post-crash peak in pressure fluctuations, δP rms serves as a critical metric for identifying transition from continuous turbulence fluctuations to ELM bursting. Larger δP rms values correlate with ELM bursts driven by local or global instabilities, whereas smaller values indicate turbulence-dominated transport. Drift-Alfvén and resistive ballooning turbulence enhance the entrainment of fluctuations from the pedestal to the SOL, contributing to the complex interplay of dynamics in this regime. These findings emphasize the importance of separatrix density shaping and pedestal gradient control for optimizing ELM behavior in ITER and future fusion devices.

Li, Nami [Lawrence Livermore National Laboratory (↗
DOE OSTI2024

Deep Koopman Neural Network for Analyzing High-Energy-Density Simulations of Electrical Wire Explosions

Megaampere-scale electrical wire experiments (EWEs) provide a platform for studying magnetohydrodynamic (MHD) instability growth in magneto-inertial fusion (MIF) devices. Even when nonlinear simulations of these experiments can digitally reproduce much of the experimentally observed instability growth, interpreting the results and understanding mode growth and evolution can be non-trivial. As a first step toward providing better interpretation of these simulation features, this work investigates the use of a deep neural network that uses Koopman operator theory to analyze the dynamics of pulsed-power-driven explosions of EWEs. This deep neural network is trained on 1-D resistive MHD simulations of EWEs. This neural network learns to transform the nonlinear data into a lower-dimensional representation where the time dynamics are linear. Layers of this neural network are shown to learn features of the simulations, including the locations of shock waves and different physical regimes of the simulation. Using the learned features, the network can compress a time state of the simulation consisting of 5120 data point into a 36-parameter lower-dimensional latent space embedding. Furthermore, these embeddings are shown to be clustered in the latent space by initial radius and time state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the frequency bifurcations of the MHD startup modes in NSTX

The observed bifurcations of the low frequency (<50 kHz) and low toroidal periodicity (n < 5) magnetohydrodynamic (MHD) activity often present in the initial part of the National Spherical Tokamak Experiment (NSTX) discharges can be explained by the evolution of the radial profile of the safety factor (q=rB φ /RB θ ) crossing multiple rational surfaces in the core. Important performance limiting instability mechanisms in the NSTX spherical tokamak are often linked to low frequency and low-n MHD activity. They are quite common in long-pulse NSTX plasmas. They can be present at the beginning of the plasma current flat-top, at the end of the discharge or during the whole duration, and they have been observed to deleteriously impact performance over a wide range of q 95 . An interesting feature observed in some NSTX discharges is the presence of a bifurcation in the frequency of the low n modes, as low as n = 1, that have frequencies comparable to the plasma core rotation divided by n. Equilibrium reconstructions constrained by magnetic diagnostics data and motional stark effect pitch angle radial profiles suggest that the observed bifurcations are linked to a fast evolving minimum value of q. Finally, 3D non-linear resistive MHD simulations show that these modes are ideal and exist as non-resonant before the correspondent rational surface enters the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Nonlinear error field response in the presence of plasma rotation and real frequencies due to favorable curvature

This works describes the interaction of an phenomenon known to arise in tokamak fusion experiments with the magnetic perturbations that can arise at the edge of such experiments. We present nonlinear NIMROD resistive MHD simulations of the response of a rotating plasma to an error field when the plasma has weakly damped linear tearing modes (TMs), stabilized by a pressure gradient and favorable curvature. The favorable curvature leads to the Glasser effect: the occurrence of real frequencies and stabilization with positive stability index Δ'.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Effects of external magnetic fields on the multi-mode Rayleigh–Taylor instability

The magneto-Rayleigh–Taylor instability (mRTI) is a key process in inertial confinement fusion and is thought to be widespread in the interstellar medium, where it can concentrate plasma into discrete structures. We present resistive MHD simulations of the nonlinear evolution of multi-mode mRTI in both two and three dimensions, examining the effects of uniform external magnetic fields oriented either parallel or perpendicular to the initial interface. In both 2D and 3D, weak parallel fields enhance mixing-zone growth, whereas stronger fields suppress it. For perpendicular fields, growth is initially inhibited but becomes enhanced at later times. These behaviors arise from magnetic tension, which modifies flow anisotropy, buoyancy, drag, and vortex dynamics. The interplay of these mechanisms governs the distinct ways in which magnetic fields influence mRTI evolution.

Flow instabilities↗
DOE OSTI2023

Species dependence of the impurity injection induced poloidal flow and magnetic island rotation in a tokamak

Abstract Recent experiments have demonstrated the species dependence of the impurity poloidal drift direction along with the magnetic island rotation in the poloidal plane. Our resistive MHD simulations have reproduced such a dependence of the impurity poloidal flow, which is found mainly determined by a local plasmoid formation due to the impurity injection. The synchronized magnetic island rotation is dominantly driven by the electromagnetic torque produced by the impurity radiation primarily through the modification to the axisymmetric components of current density.

Physics↗
DOE OSTI2020

Experimental studies of plasma-antenna coupling with the JET Alfvén Eigenmode Active Diagnostic

Here, this paper presents a dedicated study of plasma-antenna (PA) coupling with the Alfvén Eigenmode Active Diagnostic (AEAD) in JET. Stable AEs and their resonant frequencies f, damping rates γ < 0, and toroidal mode numbers n are measured for various PA separations and limiter versus X-point magnetic configurations. Two stable AEs are observed to be resonantly excited at distinct low and high frequencies in limiter plasmas. The values of f and n do not vary with PA separation. However, |γ| increases with PA separation for the low-f, but not high-f, mode, yet this may be due to slightly different edge conditions. The high-f AE is detected throughout the transition from limiter to X-point configuration, though its damping rate increases; the low-f mode, on the other hand, becomes unidentifiable. The linear, resistive MHD code CASTOR is used to simulate the frequency scan of an AEAD-like external antenna. For the limiter pulses, the high-f mode is determined to be an n = 0 GAE, while the low-f mode is likely an n = 2 TAE. During the transition from limiter to X-point configuration, CASTOR indicates that n = 1 and 2 EAEs are excited in the edge gap. These results extend previous experimental studies in JET and Alcator C-Mod; validate the computational work performed by Dvornova et al; and provide guidance for the optimization of PA coupling in upcoming JET energetic particle experiments, for which the AEAD will aim to identify the contribution of alpha particles to AE drive during the DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Beam modulation and bump-on-tail effects on Alfvén eigenmode stability in DIII-D

Beam modulation effects on Alfvén eigenmode stability have been investigated in a recent DIII-D experiment and show that variations in neutral beam modulation period can have an impact on the beam driven Alfvén eigenmode spectrum and resultant fast ion transport despite similar time-averaged input power. The experiment was carried out during the current ramp phase of L-mode discharges heated with sub-Alfvénic 50–80 kV deuterium neutral beams that drive a variety of Alfvén eigenmodes unstable. The modulation period of two interleaved beams with different tangency radii was varied from shot to shot in order to modify the relative time dependent mix of the beam pitch angle distribution as well as the persistence of a bump-on-tail feature near the injection energy (a feature confirmed by imaging neutral particle analyzer measurements). As the beam modulation period is varied from 7 ms to 30 ms on/off (typical full energy slowing down time of τ slow ≈ 50 ms at mid-radius), toroidicity-induced Alfvén eigenmodes (TAEs) located in the outer periphery of the plasma become intermittent and coincident with the more tangential beam. Core mode activity changes from reversed shear Alfvén eigenmodes (RSAEs) to a mix of RSAE and beta-induced Alfvén eigenmodes. Discharges with 30 ms on/off period do not have a persistent bump-on-tail feature, have the lowest average mode amplitude and least fast ion transport. Furthermore, detailed analysis of an individual TAE using TRANSP kick modeling (Monte Carlo evolution of the distribution function with probabilistic 'kicks' by the AEs) and the resistive MHD code with kinetic fast ions, MEGA, find no strong role of energy gradient drive due to bump-on-tail features. Instead, the observed TAE modulation with interleaved beams is likely a pitch angle dependent result combined with slowing down of the tangential beam between pulses. For the conditions investigated, bump-on-tail contributions to TAE drive were found to be 5% or less of the total drive at any given time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Controlling the size of non-axisymmetric magnetic footprints using resonant magnetic perturbations

The structure of the non-axisymmetric heat load distribution at the divertor plates is determined not only by the toroidal but also from the poloidal spectrum of non-axisymmetric field perturbations. Whether they are intrinsic, like error fields, or they are applied through 3D coils, the non-axisymmetric fields produce complex 3D edge magnetic topologies (footprints) that alter the properties of the heat and particle flux distributions on the divertor target plates. In this manuscript, a study of the impact of applied 3D field poloidal spectrum on the footprint size and structure is done for the DIII-D tokamak using the resistive MHD code M3D-C1 coupled with the field line tracing code TRIP3D. To resolve the impact of the poloidal spectrum of the magnetic perturbation, the relative phase of the two rows of in-vessel 3D coils used to produce both a n = 2 and a n = 3 perturbation is varied, where n is the toroidal harmonic of the magnetic perturbation. This shows that the largest footprint is predicted when the relative phase of the two rows is close to zero, which is also where the resonant coupling with the plasma is maximized. These results suggest that it will be challenging to decouple the footprint size from the requisite resonant coupling for RMP–ELM control. The correlation between the measured heat load and particle flux distributions at the outer divertor plates in DIII-D and the magnetic measurements is in good agreement with the predicted dependence of the magnetic footprint size on the amplitude of the resonant component of the plasma response.

3D fields↗
DOE OSTI2022

Toroidal modeling of 3D perturbations generated by current filaments in scrape-off layer in tokamak with biased divertor targets

The divertor biasing technique offers a promising alternative to control the edge localized mode (ELM) as well as the divertor heat load in tokamaks, as compared to the resonant magnetic perturbation (RMP) generated by magnetic coils. The linear resistive MHD code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) is employed to study the plasma response to the n = 2 (n is the toroidal mode number) field perturbation in HL-2A, produced by the divertor biasing current filaments flowing in the scrape-off layer region. The response field due to biasing currents is compared with the RMP field utilized for ELM control in HL-2A. The strength of the radial resonant field perturbation, produced by the biasing currents of 100 A level, is found to be comparable with RMP generated by several kA of ELM control coil currents for the reference plasma. The plasma normal displacement near the X-point and the associated neoclassical toroidal viscosity torque are also computed to be similar between these two techniques. The modeling results thus strongly suggest that the biasing technique can be applied to control ELMs. Moreover, the biasing currents produce field perturbations, including the plasma response, that are localized more near the plasma edge (compared to the RMP counterpart), thus reducing the chance of mode locking associated with core perturbations. Particle orbit tracing also reveals that the biasing current produced magnetic perturbation tends to widen the heat deposition region and induce the strike point splitting of the ion saturation flow on the outer divertor surface, consistent with experimental observations in HL-2A. These toroidal modeling results confirm the possibility of ELM control and plasma exhaust solution by the divertor biasing technique.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Magnetic island formation and rotation braking induced by low-Z impurity penetration in an EAST plasma

Abstract Recent observations of the successive formations of the 4 / 1 , 3 / 1 , and 2 / 1 magnetic islands as well as the subsequent braking of the 2 / 1 mode during a low- Z impurity penetration process in EAST experiments are well reproduced in our 3 D resistive MHD simulations. The enhanced parallel current perturbation induced by impurity radiation predominately contributes to the tearing mode growth, and the 2 / 1 island rotation is mainly damped by the impurity accumulation as results of the influence from high n modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Plasmoid-mediated reconnection during nonlinear peeling–ballooning edge-localized modes

Abstract Plasmoid-mediated reconnection is investigated for nonlinear peeling–ballooning edge-localized modes (P-B ELMs). The formation of current sheets and the transition to 3D current sheet instability is demonstrated through fully nonlinear resistive MHD simulations of P-B ELMs in DIII-D discharges. Large-scale axisymmetric current sheets, as well as small-scale poloidally extended current sheets, are formed as the coherent P-B ELM filaments nonlinearly evolve. It is observed that, at high Lundquist numbers, these current sheets break during a reconnection burst, i.e. a secondary exponential growth of intermediate modes followed by relaxation due to the suppression of P-B drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

3D modeling of n = 1 RMP driven heat fluxes on the SPARC tokamak PFCs using HEAT

3D heat flux calculations at the lower outer divertor plate of SPARC using the HEAT code show that 3D fields generated from error field correction coils can lead to enhanced peak heat fluxes up to 15 times larger compared to the axisymmetric case. Previously employed to simulate axisymmetric heat flux on 3D plasma facing components, the HEAT code can now predict 3D heat flux generated by non-axisymmetric plasmas. This is achieved via a new HEAT module which leverages the 3D field line tracing capabilities of MAFOT starting from an M3D-C1 (MHD resistive code) perturbed equilibria. The resulting heat flux is assigned using the magnetic footprint and the heat flux layer model, an extension of the 2D heat flux model also known as the Eich, to 3D non-axisymmetric plasmas. For SPARC, the new capabilities of HEAT are used to calculate the 3D heat loads resulting from n = 1 perturbation fields (with n indicating the toroidal periodicity) applied through a toroidal array of six picture frame coils with different amplitude. The comparison with the unperturbed case shows significant changes in shape and intensity of the heat flux profile. The results show that the application of n = 1 3D field leads to a localized enhancement of the heat flux peak, influenced by the wetted area impacted by the magnetic footprint, and the appearance of a secondary heat flux peak, whose intensity depends on amplitude of the applied 3D field and toroidal location.

3D heat flux↗
DOE OSTI2023

Contribution to: UNM-UNR-VT Project Research Performance Progress Report 'Seeding and Evolution of Magnetohydrodynamic Instabilities of a Metal Surface Driven by Intense Current' Year 4 (no-cost extension for FY 2022): Oct 2021 - Sept 2022

The multiphysics code Ares has been used to perform resistive MHD fluid simulations of the METI-II experiments (described in Section 1.1.2). The simulations include both 1D Lagrangian simulations and axisymmetric 2D ALE simulations. The results of these simulations are summarized here and a full manuscript is nearing completion for submission.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Collaborative Research: Unveiling the Interplay Between Magnetic Reconnection and Turbulence — Theory and Modeling

This is the final technical report for the project "Collaborative Research: Unveiling the Interplay Between Magnetic Reconnection and Turbulence — Theory and Modeling." This project investigated the interplay between two fundamental plasma processes — turbulence and magnetic reconnection — and the associated heating and particle acceleration. In particular, we focused on a new regime where the plasmoid instability mediates the turbulent energy cascade. In this new regime, reconnecting current sheets are disrupted by the growth of plasmoids/flux ropes on time scales shorter than the typical turnover times of eddies, speeding up the energy cascade and steepening the turbulence energy spectrum. This project addressed the following outstanding open questions: (1) What are the fundamental differences between plasmoid‐mediated turbulence cascade in three dimensions (3D) and two dimensions (2D)? (2) What are the essential features of plasmoid‐mediated turbulence in weakly collisional or collisionless regimes beyond resistive MHD? (3) How does the interplay of turbulence and reconnection regulate the energy release and dissipation in space and astrophysical plasmas? To address these questions, we have carried out three investigations: (1) Three-dimensional plasmoid-mediated turbulence energy cascade; (2) plasmoid-mediated reconnection and turbulence in three-dimensional Hall MHD; (3) reconnection and heating in coronal loops. This technical report summarizes the key findings of these investigations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗