Block Preconditioning for Magnetic Confinement Fusion Relevant Resistive MHD Simulations.
Abstract not provided.
Ohm, Peter↗
SEARCH · Search NASA
Results for “MHD simulation”
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.
The effect of time-dependent level kinetics in a 1-D MHD simulation of an argon gas z-pinch implosion*
Explore the source record for details and available documents.
Schaeuble, Marc-Andre↗
3D MHD Simulations for Jet Experiments: DOE INFUSE [Slides]
Abstract not provided.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Understanding Electrothermal Instability Growth by Comparing Z-Pinches with Engineered Defects to 3D-MHD Simulations.
Abstract not provided.
Awe, Thomas James↗
Pedestal collapse by resonant magnetic perturbations
Pedestal collapse (i.e., the complete loss of the edge transport barrier (ETB)) in DIII-D H-mode plasmas occurs when resonant magnetic perturbations (RMPs) penetrate the steep gradient region at the plasma edge. Normally, RMP driven magnetic islands can occur at the top and bottom of the H-mode pedestal and these islands generate conditions consistent with edge-localized-mode (ELM) suppression and density pump-out, respectively, based on nonlinear two-fluid MHD simulations. In contrast, MHD simulations show that the steep pressure gradient region between the top and bottom of the DIII-D pedestal is generally immune to resonant field penetration due to large local E × B and diamagnetic flows. By this fortuitous circumstance, the edge-transport-barrier and H-mode confinement can be maintained while achieving ELM suppression. However, pedestal collapse can occur in DIII-D when the screening flows are inadequate to prevent field penetration in the steep gradient region of the pedestal. Non-linear two-fluid MHD simulations support the role of resonant field penetration in pedestal collapse for DIII-D H-mode plasmas with weak edge E × B and diamagnetic screening flows. ITER will likely have weaker edge screening flows than present experiments due to its much larger size, making it more susceptible to resonant field penetration in the steep gradient region of the pedestal. Analysis of model ITER equilibria demonstrates that resonant field penetration in the steep pressure gradient region is possible for RMP levels of the order required for ELM suppression. The effect of such penetration on the ITER pedestal will depend sensitively on the resulting degree of island overlap.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Role of edge-localized neoclassical tearing modes in quiescent H-mode plasmas in the DIII-D tokamak
The edge-harmonic oscillations (EHOs) in standard quiescent H-mode (QH-mode) plasmas in DIII-D are consistent with edge-localized neoclassical tearing modes (NTMs) based on nonlinear two-fluid MHD simulations. Using kinetic equilibria constrained by edge profile measurements, the MHD simulations show that the n = 1 NTM and its harmonics can be destabilized at the pedestal top of QH-mode plasma by the edge bootstrap current. The simulations further show that the unstable NTMs can saturate either at small (<2% ψ N ) or large (>4% ψ N ) island width depending on the magnitude of the edge bootstrap current, where ψ N is the normalized radius in poloidal flux. The onset of the EHO also results in a prompt decrease in the pedestal width and height, consistent with simulation results for the onset of the NTM at the top of the QH-mode pedestal. This suggests that the avoidance of edge-localized modes (ELMs) in QH-mode can be attributed to the enhanced local transport induced by the NTM that is sufficient to prevent the expansion of the pedestal to an unstable width, analogous to the mechanism explored for ELM suppression by resonant magnetic perturbations. Nonlinear MHD simulations scanning the E × B frequency and the ratio of parallel and perpendicular thermal diffusivity (χ | /χ ⊥ ) at the pedestal top show that edge-localized NTMs are destabilized for conditions of high E × B frequency, high pedestal temperature, and low pedestal density, qualitatively consistent with experimental conditions required for observing the EHO.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Computational study of runaway electrons in MST tokamak discharges with applied resonant magnetic perturbation
A numerical study of magnetohydrodynamics (MHD) and tracer-particle evolution investigates the effects of resonant magnetic perturbations (RMPs) on the confinement of runaway electrons (REs) in tokamak discharges conducted in the Madison Symmetric Torus. In computational results of applying RMPs having a broad toroidal spectrum but a single poloidal harmonic, m = 1 RMP does not suppress REs, whereas m = 3 RMP achieves significant deconfinement but not the complete suppression obtained in the experiment. MHD simulations with the NIMROD code produce sawtooth oscillations, and the associated magnetic reconnection can affect the trajectory of REs starting in the core region. Simulations with m = 3 RMP produce chaotic magnetic topology over the outer region, but the m = 1 RMP produces negligible changes in field topology, relative to applying no RMP. Using snapshots of the MHD simulation fields, full-orbit relativistic electron test particle computations with KORC show ≈50% loss from the m = 3 RMP compared to the 10%–15% loss from the m = 1 RMP. Here, test particle computations of the m = 3 RMP in the time-evolving MHD simulation fields show correlation between MHD activity and late-time particle losses, but total electron confinement is similar to computations using magnetic-field snapshots.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
Transition in particle transport under resonant magnetic perturbations in a tokamak
Abstract Nonlinear 3D MHD simulations and validations reveal that the hybrid particle-MHD transport is a key process for driving the pump-out in the presence of Resonant Magnetic Perturbations (RMPs) in the KSTAR tokamak. Particle transport and the resulting density pump-out by RMPs are shown to be composed of not only the classical flow convection near magnetic islands due to polarization but also the neoclassical ion diffusion across perturbed magnetic surfaces. The latter is known as the Neoclassical Toroidal Viscosity (NTV) and is integrated into nonlinear MHD simulations here for the first time, revealing that the two-stage pump-outs observed in KSTAR experiments are reproduced only with such integrated nonlinear MHD and transport evolution. Near-resonant responses, which have received less attention than the resonant response, play distinct roles in the pump-out along with the island formation. In addition, this modeling is used to investigate the pump-outs in double-null-like plasmas and numerically capture the effect of the double-null shape on the pump-outs, which may explain the difficulty of Edge Localized Mode (ELM) suppression access in double-like plasmas. This reveals new aspects of the impact toroidal geometry and mode coupling have on 3D physics and reveals the importance of near-resonant components in suppressing ELMs.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Assessment of dynamic-screw-pinch-driven, current-scaled MagLIF target implosion performance using 3D magnetohydrodynamic simulations
Analytic studies and two-dimensional “clean” radiation-magnetohydrodynamic (rad-MHD) simulations employing dynamical similarity driver-target scaling prescriptions [Ruiz et al., Phys. Plasmas 30, 032708 (2023)] suggest that Magnetized Liner Inertial Fusion (MagLIF) target implosions can scale to > 10 MJ DT fusion yields when peak drive current is increased beyond 60 MA. We present results from three-dimensional (3D) rad-MHD simulations of similarity-scaled MagLIF target implosions at peak drive currents ranging from 15 to 40 MA. Simulations in this study suggest that magneto-Rayleigh–Taylor instability (MRTI) growth and feedthrough to the fuel region are more severe at higher drive current scales, which reduces the fusion yield compared to prior analytic and 2D clean simulation predictions. In contrast to standard MagLIF, simulations of current-scaled MagLIF target implosions driven by a dynamic screw pinch (DSP) demonstrate reduced MRTI feedthrough and greater fuel magnetization, resulting in improved thermonuclear performance and enhanced performance scaling with peak drive current. DSP drive enables additional scaling of the liner mass to increase liner radius but maintain implosion time, resulting in higher implosion velocities at the expense of increased susceptibility to MRTI. We present a current- and mass-scaled simulated DSP-MagLIF target implosion at the ∼ 40 MA peak current level that produces ignition scale performance, demonstrating a burn-averaged Lawson ignition parameter above unity and DT fusion yield above 1 MJ.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Sensitivity of magnetohydrodynamic simulations of Joule-heated conductors to the vaporization curve in equations of state
Magnetohydrodynamic (MHD) simulations of electrically exploded aluminum and copper rods demonstrate a technique to validate equations of state (EOS) for rapidly Joule-heated conductors. The balance of internal and magnetic forces at the conductor-insulator interface drives the metal there along the vaporization phase boundary. Variations between critical points and vaporization curves in existing models predict differing densities and temperatures in MHD simulations for these models. Here, the inclusion of Maxwell constructs in the liquid-vapor biphase region of the EOS caused the rod surface to vaporize earlier in time than unmodified tables with van der Waals loops. Velocimetry of recent experiments is used to validate the location of the vaporization curve in existing EOS models and differentiate between the vapor dome treatments. Dielectric coatings applied to the metal surface restricted the conductor’s expansion and diverted the metal into the warm dense matter regime.
36 MATERIALS SCIENCE↗
MHD modeling of shattered pellet injection in JET
Abstract Nonlinear 3D MHD simulations of shattered-pellet injection (SPI) in JET show prototypical SPI-driven disruptions using the M3D-C1 and NIMROD extended-MHD codes. Initially, radiation-driven thermal quenches are accelerated by MHD activity as the pellet crosses rational surfaces, leading to a radiation spike, global stochasticization of the magnetic field, and a complete thermal quench. Eventually, current quenches, preceded by a current spike are seen as the Ohmic heating becomes equal to the radiative cooling. The results are qualitatively similar for both a single monolithic pellet, pencil-beam model, and a realistic shatter to represent the SPI plume. A scan in viscosity from 500 to 2000 m 2 s −1 for MHD simulations finds that reducing viscosity increases MHD activity and decreases thermal quench time slightly. A realistic cloud of fragments modeling shows that mixed-D–Ne pellet travels deeper into the plasma core before the thermal quench. At the slow pellet speeds, the pellet is found to be moving slowly enough inward that even the 5% neon in the mixed pellet is enough to effectively radiate the thermal energy available. Radiation toroidal peaking is predicted to be at levels consistent with experimental observations and reduced as the pellet travels deeper into the plasma. These simulations lay the ground work for more-sophisticated validative and predictive modeling of SPI in JET using both M3D-C1 and NIMROD.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A finite-volume scheme for modeling compressible magnetohydrodynamic flows at low Mach numbers in stellar interiors
Fully compressible magnetohydrodynamic (MHD) simulations are a fundamental tool for investigating the role of dynamo amplification in the generation of magnetic fields in deep convective layers of stars. The flows that arise in such environments are characterized by low (sonic) Mach numbers ($\mathscr{M}$ son ≲ 10 –2 ). In these regimes, conventional MHD codes typically show excessive dissipation and tend to be inefficient as the Courant–Friedrichs–Lewy (CFL) constraint on the time step becomes too strict. In this work we present a new method for efficiently simulating MHD flows at low Mach numbers in a space-dependent gravitational potential while still retaining all effects of compressibility. The proposed scheme is implemented in the finite-volume SEVEN-LEAGUE HYDRO (SLH) code, and it makes use of a low-Mach version of the five-wave Harten–Lax–van Leer discontinuities (HLLD) solver to reduce numerical dissipation, an implicit–explicit time discretization technique based on Strang splitting to overcome the overly strict CFL constraint, and a well-balancing method that dramatically reduces the magnitude of spatial discretization errors in strongly stratified setups. The solenoidal constraint on the magnetic field is enforced by using a constrained transport method on a staggered grid. We carry out five verification tests, including the simulation of a small-scale dynamo in a star-like environment at $\mathscr{M}$ son ~ 10 –3 . We demonstrate that the proposed scheme can be used to accurately simulate compressible MHD flows in regimes of low Mach numbers and strongly stratified setups even with moderately coarse grids.
79 ASTRONOMY AND ASTROPHYSICS↗
Long-term general relativistic magnetohydrodynamics simulations of magnetic field in isolated neutron stars
ABSTRACT Strong magnetic fields play an important role in powering the emission of neutron stars. Nevertheless, a full understanding of the interior configuration of the field remains elusive. In this work, we present general relativistic magnetohydrodynamics (MHD) simulations of the magnetic field evolution in neutron stars lasting ${\sim } {880}\,$ms (∼6.5 Alfvén crossing periods) and up to resolutions of $0.1155\,$km using Athena++. We explore two different initial conditions, one with purely poloidal magnetic field and the other with a dominant toroidal component, and study the poloidal and toroidal field energies, the growth times of the various instability-driven oscillation modes, and turbulence. We find that the purely poloidal setup generates a toroidal field, which later decays exponentially reaching $1{{\ \rm per\ cent}}$ of the total magnetic energy, showing no evidence of reaching equilibrium. The initially stronger toroidal field setup, on the other hand, loses up to 20 per cent of toroidal energy and maintains this state till the end of our simulation. We also explore the hypothesis, drawn from previous MHD simulations, that turbulence plays an important role in the quasi-equilibrium state. An analysis of the spectra in our higher resolution setups reveals, however, that in most cases we are not observing turbulence at small scales, but rather a noisy velocity field inside the star. We also observe that the majority of the magnetic energy gets dissipated as heat increasing the internal energy of the star, while a small fraction gets radiated away as electromagnetic radiation.
79 ASTRONOMY AND ASTROPHYSICS↗
Nature of transonic sub-Alfvénic turbulence and density fluctuations in the near-Sun solar wind
Context. Recent Parker Solar Probe (PSP) measurements have revealed that solar wind (SW) turbulence transits from a subsonic to a transonic regime near the Sun, while remaining sub-Alfvénic. These observations call for a revision of the existing SW models, where turbulence is considered to be both subsonic and sub-Alfvénic. Aims. In this work, we introduce a new magnetohydrodynamic (MHD) model of transonic sub-Alfvénic turbulence (TsAT). Methods. We used 3D MHD simulations initialized with parameters measured by PSP to investigate the properties of the new near-Sun SW transonic turbulent regime. We then derived a reduced set of MHD equations in the transonic sub-Alfvénic limit to interpret our numerical results. Results. Our TsAT model shows that turbulence is effectively nearly incompressible (NI) and has a 2D + slab (quasi-2D) geometry not only in the subsonic limit, but also in the transonic regime, as long as it remains sub-Alfvénic, a condition essentially enforced everywhere in the heliosphere by the strong local magnetic field. These predictions are consistent with 3D MHD simulations, showing that transonic turbulence is dominated by low-frequency quasi-2D incompressible structures, while compressible fluctuations are a minor component corresponding to low-frequency slow modes and high-frequency fast modes. Conclusions. Our new TsAT model extends existing NI theories of turbulence, and is potentially relevant for the theoretical and numerical modeling of space and astrophysical plasmas, including the near-Sun SW, the solar corona, and the interstellar medium.
79 ASTRONOMY AND ASTROPHYSICS↗
UNR Grant Proposal Contribution
The work of Carrier et. al. (citation needed for Matt’s effect of surface roughness on phase transitions/ETI paper) demonstrates that 1D and 2D resistive magnetohydrodynamic (MHD) simulations can reliably model exploding aluminum rods driven by megaampere currents. Close agreement between simulation results and photonic Doppler velocimetry data from the Mykonos electrothermal instability (METI-II) campaign builds confidence in predictive modelling capabilities for pulsed-power HED experiments. Furthermore, 2D MHD simulations show how machined features and micro-scale surface roughness can seed the electrothermal instability (ETI). Surface roughness was observed to reduce the time of melt by 19%, with ETI growth driving the enhanced heating of the rod surface.
36 MATERIALS SCIENCE↗
Burn propagation in magnetized high-yield inertial fusion
Recent experiments at the National Ignition Facility (NIF) have demonstrated ignition for the first time in an inertial confinement fusion (ICF) experiment, a major milestone allowing the possibility of high energy gain through burn propagation. Use of external magnetic fields, applied primarily to reduce thermal losses, could increase hotspot temperature and ease requirements for ignition, opening up the capsule design space for high energy gain. However, this same restriction of thermal transport has the potential to inhibit burn propagation, which is vital in the attainment of high gain. In this work, radiation-magnetohydrodynamics (MHD) simulations carried out using the code Chimera are used to investigate the effect of a pre-imposed magnetic field on ignition and burn propagation. This paper studies the propagation of burn using both an idealized planar model and in fully integrated 2D MHD simulations of an igniting NIF capsule. A study of magnetized burn propagation in the idealized planar model identifies three regimes of magnetized burn propagation: (1) thermal conduction driven; (2) alpha transport driven; and (3) fully suppressed burn. Simulations of NIF shot N210808 with an applied 40 T axial field show clear indication of burn suppression perpendicular to field lines, with rapid burn observed along field lines. Implosion shape is altered by the field, and anisotropic conduction causes significant modification to the rate of ablation during stagnation. These results highlight the fundamental changes to implosion dynamics in high-yield magnetized ICF and motivate further study to better optimize future magnetized target designs for high gain.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Quantification of MagLIF morphology using the Mallat scattering transformation
The morphology of the stagnated plasma resulting from magnetized liner inertial fusion is measured by imaging the self-emission x rays coming from the multi-keV plasma. Equivalent diagnostic responses can be generated by integrated radiation-magnetohydrodynamic (rad-MHD) simulations from programs such as HYDRA and GORGON. There have been only limited quantitative ways to compare the image morphology, that is the texture, of simulations and experiments. We have developed a metric of image morphology based on the Mallat scattering transformation (MST), a transformation that has proved to be effective at distinguishing textures, sounds, and written characters. This metric is designed, demonstrated, and refined by classifying ensembles (i.e., classes) of synthetic stagnation images and by regressing an ensemble of synthetic stagnation images to the morphology (i.e., model) parameters used to generate the synthetic images. We use this metric to quantitatively compare simulations to experimental images, experimental images to each other, and to estimate the morphological parameters of the experimental images with uncertainty. This coordinate space has proved to be very adept at doing a sophisticated relative background subtraction in the MST space. This was needed to compare the experimental self-emission images to the rad-MHD simulation images.
Glinsky, Michael E. (ORCID:0000000324933326)↗