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 163 records · Page 9

Magnetohydrodynamic Wave Mode Conversion in a Pseudostreamer Structure

The interaction of a fast magnetoacoustic wave with a magnetic null point in a pseudostreamer magnetic configuration is modeled. The incident fast wave is excited by an impulsive velocity driver located in the corona outside the Alfvén-acoustic equipartition layer. The fast wave front experiences refraction and wraps around the null point, while the front segment approaching the null point along the magnetic bisector passes through it. In the equipartition layer surrounding the null point, part of the fast wave energy is converted into slow magnetoacoustic waves. Slow waves propagate outward from the null point along the magnetic separatrices. The fronts of slow waves exhibit a characteristic chevron shape. The fast-to-slow mode conversion efficiency increases with the width of the initial perturbation, following a power-law dependence with an index of approximately 1.06. For initial wave amplitudes below about 10% of the local Alfvén speed, the dependence of the conversion efficiency on the initial amplitude is weak. After interacting with the null point, the amplitudes of the fast waves decrease significantly owing to energy conversion and the spreading of wave energy over a broad spatial domain. In contrast, the energy of slow waves remains concentrated in the vicinity of the magnetic separatrices, resulting in less pronounced geometric attenuation of these waves. These findings are consistent with observational reports of the interaction of a coronal fast wave with a pseudostreamer.

Zhong, Yu [Univ. of Warwick, Coventry (United King↗

Advancing our understanding of instabilities in high-energy-density systems through the study of benchmark datasets with advanced modeling tools

Numerous types of pulsed power driven inertial confinement fusion (ICF) and high energy density (HED) systems rely on implosion stability to achieve desired temperatures, pressures, and densities. Sandia National Laboratories Pulsed Power Sciences Center’s main ICF platform, Magnetized Liner Inertial Fusion (MagLIF), suffers from implosion instabilities which limit attainable fuel conditions and can compromise fuel confinement. This Truman Fellowship research primarily focused on computationally exploring (a) methods for improving our understanding of hydrodynamic and magnetohydrodynamic instabilities that form during cylindrical liner implosions, (b) methods for mitigating implosion instabilities, particularly those that degrade performance of MagLIF targets, and (c) novel MagLIF target designs intended to improve target performance primarily via enhanced implosion stability. Several multi-dimensional computational tools were used, including the magnetohydrodynamics code ALEGRA, the radiation-magnetohydrodynamics code HYDRA, and the magnetohydrodynamics code KRAKEN. This research succeeded in executing and analyzing simulations of automagnetizing liner implosions, shockless MagLIF implosions, dynamic screw pinch driven cylindrical liner implosions, and cylindrically convergent HED instability studies. The methods and tools explored and developed in this Truman Fellowship research have been published in several peer-reviewed journal articles and will serve as useful contributions to the fields of pulsed power science and engineering, particularly pertaining to pulsed power ICF and HED science.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Observation of magnetic islands in tokamak plasmas during the suppression of edge-localized modes

In tokamaks, a leading platform for fusion energy, periodic filamentary plasma eruptions known as edge-localized modes occur in plasmas with high-energy confinement and steep pressure profiles at the plasma edge. These edge-localized modes could damage the tokamak wall but can be suppressed using small three-dimensional magnetic perturbations. Here we demonstrate that these magnetic perturbations can change the magnetic topology just inside the steep gradient region of the plasma edge. We identify signatures of a magnetic island, and their observation is linked to the suppression of edge-localized modes. We compare high-resolution measurements of perturbed magnetic surfaces with predictions from ideal magnetohydrodynamic theory where the magnetic topology is preserved. Although ideal magnetohydrodynamics adequately describes the measurements in plasmas exhibiting edge-localized modes, it proves insufficient for plasmas where these modes are suppressed. Nonlinear resistive magnetohydrodynamic modelling supports this observation. Our study experimentally confirms the predicted role of magnetic islands in inhibiting the occurrence of edge-localized modes. This will be beneficial for physics-based predictions in future fusion devices to control these modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical study of $\delta$-function current sheets arising from resonant magnetic perturbations

We report general three-dimensional toroidal ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces are susceptible to forming singular current sheets when resonant perturbations are applied. The presence of singular current sheets indicates that, in the presence of non-zero resistivity, magnetic reconnection will ensue, leading to the formation of magnetic islands and potentially regions of stochastic field lines when islands overlap. Numerically resolving singular current sheets in the ideal magnetohydrodynamics (MHD) limit has been a significant challenge. This work presents numerical solutions of the Hahm–Kulsrud–Taylor (HKT) problem, which is a prototype for resonant singular current sheet formation. The HKT problem is solved by two codes: a Grad–Shafranov (GS) solver and the Stepped Pressure Equilibrium Code (SPEC) code. The GS solver has built-in nested flux surfaces with prescribed magnetic fluxes. The SPEC code implements multi-region relaxed magnetohydrodynamics (MRxMHD), whereby the solution relaxes to a Taylor state in each region while maintaining force balance across the interfaces between regions. As the number of regions increases, the MRxMHD solution appears to approach the ideal MHD solution assuming a continuum of nested flux surfaces. We demonstrate agreement between the numerical solutions obtained from the two codes through a convergence study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Direct prediction of nonlinear tearing mode saturation using a variational principle

It is shown that the variational principle of multi-region relaxed magnetohydrodynamics (MRxMHD) can be used to predict the stability and nonlinear saturation of tearing modes in strong guide field configurations without resolving the dynamics and without explicit dependence on the plasma resistivity. While the magnetic helicity is not a good invariant for tearing modes, we show that the saturated tearing mode can be obtained as an MRxMHD state of a priori unknown helicity by appropriately constraining the current profile. The predicted saturated island width in a tearing-unstable force-free slab equilibrium is shown to reproduce the theoretical scaling at small values of Δ' and the scaling obtained from resistive magnetohydrodynamics simulations at large Δ'.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Scaling theory of three-dimensional magnetic reconnection spreading

We develop a first-principles scaling theory of the spreading of three-dimensional (3D) magnetic reconnection of finite extent in the out of plane direction. This theory addresses systems with or without an out of plane (guide) magnetic field, and with or without Hall physics. The theory reproduces known spreading speeds and directions with and without guide fields, unifying previous knowledge in a single theory. New results include: (1) Reconnection spreads in a particular direction if an x-line is induced at the interface between reconnecting and non-reconnecting regions, which is controlled by the out of plane gradient of the electric field in the outflow direction. (2) The spreading mechanism for anti-parallel collisionless reconnection is convection, as is known, but for guide field reconnection it is magnetic field bending. We confirm the theory using 3D two-fluid and resistive-magnetohydrodynamics simulations. (3) The theory explains why anti-parallel reconnection in resistive-magnetohydrodynamics does not spread. (4) The simulation domain aspect ratio, associated with the free magnetic energy, influences whether reconnection spreads or convects with a fixed x-line length. (5) We perform a simulation initiating anti-parallel collisionless reconnection with a pressure pulse instead of a magnetic perturbation, finding spreading is unchanged rather than spreading at the magnetosonic speed as previously suggested. The results provide a theoretical framework for understanding spreading beyond systems studied here, and are important for applications including two-ribbon solar flares and reconnection in Earth’s magnetosphere.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Benchmarking NIMROD continuum kinetic formulations through the steady-state poloidal flow

In this work, continuum kinetic formulations are employed as a mechanism to include closure physics in an extended magnetohydrodynamics model. Two continuum kinetic approaches have been implemented in the plasma fluid code NIMROD [Sovinec et al., “Nonlinear magnetohydrodynamics with high-order finite elements,” J. Comput. Phys. 195, 355 (2004)] including a Chapman–Enskog-like (CEL) formulation and a more conventional df approach. Ion kinetic closure schemes are employed to describe the neoclassical flow properties in axisymmetric toroidal geometry. In particular, predictions for steady-state values of poloidal flow profiles in tokamak geometry are provided using both the df formulation and two different solution techniques for the CEL approach. These results are benchmarked against analytic theory predictions as well as results from the drift kinetic code DK4D. The continuum kinetic formulations employed here show agreement with both the analytic theory and DK4D results, and offer a novel velocity space representation involving higher-order finite elements in pitch angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Growing neoclassical tearing modes seeded via transient-induced-multimode interactions

Nonlinear extended magnetohydrodynamics simulations demonstrating seeding of neoclassical tearing modes (NTMs) via magnetohydrodynamics-transient-induced multimode interactions are presented. Simulations of NTMs are enabled by two recent NIMROD code developments: the implementation of heuristic neoclassical stresses and the application of transient magnetic perturbations (MPs) at the boundary. NTMs are driven unstable by the inherently pressure driven kinetic bootstrap current, which arises due to collisional viscosity between passing and trapped electrons. These simulations use heuristic closures that model the neoclassical electron and ion stresses. NTM growth requires a seed island, which is generated by a transiently applied MP in simulations. The capability is demonstrated using kinetic-based reconstructions with flow of a DIII-D ITER Baseline Scenario discharge (La Haye et al., in Proceedings IAEA FEC, 2020). The applied MP seeds a 2/1 NTM that grows in two phases: a slow growth phase followed by a faster robust growth phase. Additionally, an evolving sequence of higher order core modes are excited at first. Power transfer analysis shows that nonlinear interactions between the core modes and the 2/1 helps drive the initial slow growth. Once the induced 2/1 magnetic island reaches a critical width, the NTM transitions to faster robust growth, which is well described by the nonlinear modified Rutherford equation. Furthermore, this work highlights the role of nonlinear mode coupling in seeding NTMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hall instability driven seeding of helical magneto-Rayleigh–Taylor instabilities in axially premagnetized thin-foil liner Z-pinch implosions

Helical magneto-Rayleigh–Taylor instability (MRTI) structures have been observed in z-pinch-driven liner implosion experiments with a pre-imposed axial magnetic field. We show that the formation of these helical structures can be described by a Hall magnetohydrodynamical (HMHD) model. We used the 3D extended magnetohydrodynamics simulation code PERSEUS (which includes Hall physics) [Seyler and Martin, Phys. Plasmas 18, 012703 (2011)] to study these helical instabilities and show that a Hall interchange instability in low-density coronal plasma immediately surrounding the dense liner is responsible for producing helically oriented effects in the magnetic field and current density within the coronal layer. Here, this seeds the helical pitch angle of the MRTI even when other proposed helical seeding mechanisms are either not present in the experiments or not accounted for in the simulations. For example, this mechanism does not require low-density power-feed plasmas to be swept in from large radius or the development of electrothermal instabilities. The Hall Instability is, thus, a new, independent explanation for the origin of the helical instabilities observed in axially premagnetized liner experiments. Simulation results supporting this mechanism are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MHD-induced beta limits in the Large Helical Device

Using the extended-magnetohydrodynamics code, M3D-C1, we perform a systematic numerical study of the effect of externally applied heating on the achievable plasma beta in a ten field-period heliotron. Heat sources of varying intensity are applied to a vacuum magnetic field that is representative of the standard configuration of the Large Helical Device, with R 0 = 3.66 m, where R 0 is the radial position of the magnetic axis in vacuum. As the system is driven to a state that is unstable to low-n magnetohydrodynamic (MHD) modes, nonlinear mode interactions lead to the formation of chaotic magnetic fields. With sufficiently strong heating, a collapse of the electron temperature profile is observed. This demonstrates the necessity of simulating the self-consistent evolution of plasma profiles, without imposing assumptions on the structure of the magnetic field, to accurately determine transport properties in stellarator plasmas. It also highlights the value of these advanced simulation capabilities for accelerating the development of high-performance stellarator operating scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonresonant particle acceleration in strong turbulence: Comparison to kinetic and MHD simulations

Collisionless, magnetized turbulence offers a promising framework for the generation of nonthermal high-energy particles in various astrophysical sites. Yet, the detailed mechanism that governs particle acceleration has remained subject to debate. By means of 2D and 3D particle-in-cell, as well as 3D (incompressible) magnetohydrodynamic (MHD) simulations, we test here a recent model of nonresonant particle acceleration in strongly magnetized turbulence, which ascribes the energization of particles to their continuous interaction with the random velocity flow of the turbulence, in the spirit of the original Fermi model. To do so, we compare, for a large number of particles that were tracked in the simulations, the predicted and the observed histories of particles momenta. The predicted history is that derived from the model, after extracting from the simulations, at each point along the particle trajectory, the three force terms that control acceleration: the acceleration of the field line velocity projected along the field line direction, its shear projected along the same direction, and its transverse compressive part. Overall, we find a clear correlation between the model predictions and the numerical experiments, indicating that this nonresonant model can successfully account for the bulk of particle energization through Fermi-type processes in strongly magnetized turbulence. Additionally we also observe that the parallel shear contribution tends to dominate the physics of energization in the particle-in-cell simulations, while in the magnetohydrodynamic incompressible simulation, both the parallel shear and the transverse compressive term provide about equal contributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Numerical Analysis of Liquid Metal MHD Flow and Heat Transfer for Open-Surface Li Divertor in FNSF

Within the ongoing U.S.-based program on the development of liquid metal plasma-facing components, numerical simulations and analyses are performed to address the feasibility of the open-surface Li divertor. In the previous scoping studies (Smolentsev, 2021), heat-removal capabilities of the divertor were assessed using a simplified flow model for a slug-type velocity profile and constant flow thickness. Here, new analyses take into account forces acting on the flowing Li layer. Three reduced-order mathematical models are applied under the conditions of the U.S. Fusion Nuclear Science Facility (FNSF) to access magnetohydrodynamic (MHD) flow development effects, velocity distribution, and surface waves: 1) fully developed MHD flow; 2) quasi-2-D developing MHD flow; and 3) multiphase MHD flow. The obtained results for MHD flows and the surface heat flux computed with the plasma code scrape-off layer plasma simulation for ITER (SOLPS-ITER) are then used as input data to compute the temperature distribution in the divertor by solving the convection–diffusion energy equation.

Smolentsev, Sergey↗

Proton Acceleration in Low- β Magnetic Reconnection with Energetic Particle Feedback

Magnetic reconnection regions in space and astrophysics are known as active particle acceleration sites. There is ample evidence showing that energetic particles can take a substantial amount of converted energy during magnetic reconnection. However, there has been a lack of studies understanding the backreaction of energetic particles at magnetohydrodynamical scales in magnetic reconnection. To address this, we have developed a new computational method to explore the feedback by nonthermal energetic particles. This approach considers the backreaction from these energetic particles by incorporating their pressure into magnetohydrodynamics (MHD) equations. The pressure of the energetic particles is evaluated from their distribution evolved through Parker's transport equation, solved using stochastic differential equations (SDEs), so we coin the name MHD-SDE. Applying this method to low-β magnetic reconnection simulations, we find that reconnection is capable of accelerating a large fraction of energetic particles that contain a substantial amount of energy. When the feedback from these particles is included, their pressure suppresses the compression structures generated by magnetic reconnection, thereby mediating particle energization. Consequently, the feedback from energetic particles results in a steeper power-law energy spectrum. These findings suggest that feedback from nonthermal energetic particles plays a crucial role in magnetic reconnection and particle acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗

Demonstration of RELAP5-3D for transient analysis of a dual coolant lead lithium fusion blanket concept

Our research aims to build a multiphysics framework for transient analysis of the Dual- Coolant Lead-Lithium (DCLL) blanket design of the proposed Fusion Nuclear Science Facility (FNSF) for all Onticipated Operational Occurrences (AOOs) using Reduced Order Models (ROMs). This paper presents an assessment of using RELAP5-3D for transient thermal-hydraulic blanket analysis of the DCLL design. Preliminary validation studies for the RELAP5-3D properties were conducted using a representative vertical flow loop model and a comparison study using experimental data from the thermal convection corrosion loops at ORNL. These studies inspected the thermal hydraulic response of Lead-Lithium Eutectic (PbLi) systems within the RELAP5-3D code; mainly pressure drop and heat transfer, as compared to analytical calculations and experimental data. The Magnetohydrodynamic (MHD) pressure drop effect for liquid metal under the effects of a uniform magnetic field was implemented into our models using an equivalent forms loss coefficient based on the phasic momentum equations. We verified our implementation of MHD pressure drop within RELAP5-3D using a systematic approach consisting of two studies. Models were developed to perform comparison studies based on experimental data from the Argonne Liquid metal EXperiment (ALEX) facility and the Magnetohydrodynamic PbLi Experiment (MaPLE) facility. Using the developed validation basis, we built a simplified DCLL channel model within RELAP5-3D, inspected MHD pressure drop and performed transient temperature analysis. The DCLL channel model used radial heating data from the International Thermonuclear Experimental Reactor (ITER) which was normalized to the FNSF design. Here, a representative startup transient was performed using a relative power curve from the initial ramp up of a long-term power pulse of the Tore Supra reactor. This study demonstrated that RELAP5-3D ROMs are capable of accurate transient analysis of the DCLL blanket. We have since utilized MCNP to develop heating profiles based on the FNSF design and are working toward implementing them into our DCLL model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mechanisms of electromagnetic field control on mineral scaling in brackish water reverse osmosis: Combined homogenous and heterogeneous nucleation

Electromagnetic field (EMF) treatment has emerged as a promising approach for scaling control due to its cost-effectiveness, simplicity, and low energy consumption. However, there is a limited understanding of the mechanisms by which applied EMF impacts mineral scaling in reverse osmosis (RO) systems. This has led to inconclusive and varied results and uncertainties regarding its effectiveness. This study elucidates the impacts of EMF on homogenous and heterogeneous nucleation and membrane performance during RO desalination of different feedwaters. Our results reveal that EMF exhibits greater efficacy in treating near-saturated water (SI∼0), especially when coupled with extended hydraulic flushing (HF). For saturated brackish water desalination, heterogeneous scaling predominantly occurs on membrane surfaces, with the effectiveness of EMF in inhibiting scaling primarily attributed to the hydration effect. In supersaturated solutions, EMF promotes bulk precipitation due to the magnetohydrodynamic effect, quickly blocking membrane pores. Thus, when the saturation reaches a certain high level during RO desalination, magnetohydrodynamic EMF effects can accelerate flux decline caused by homogeneous scaling. In conclusion, this work provides an efficient method for predicting EMF efficiency, emphasizing the importance of saturation conditions and HF cleaning duration in determining membrane performance, suggesting these show promise for improving undersaturated or near-saturated feedwater desalination via RO.

42 ENGINEERING↗

Multi-scale Interaction Mechanism for Edge-Localized-Mode Suppression in the Tokamak Edge

A central challenge in fusion energy is reconciling the high-confinement mode required for reactor performance with the intense intermittent relaxation events it produces, known as edge-localized modes. These instabilities arise in the steep pressure pedestal at the plasma edge when magnetohydrodynamic thresholds are crossed, inflicting damaging heat loads on reactor components. Here, we show that multiscale interactions between microscopic turbulence and macroscopic magnetohydrodynamic modes provide encouraging prospects for self-organized edge-localized modes regulation. Using direct quantitative measurements of multiscale modes, eddy dynamics, and turbulent flux, we show that small-scale electron drift wave turbulence actively scatters the large-scale peeling-ballooning modes. This scattering decorrelates the pressure and velocity fields of the instability, so arresting its growth. Our modeling and theoretical analysis confirm this suppression mechanism is effective even when conventional linear stability thresholds are exceeded. This work establishes a nonlinear principle for edge-localized modes stability, revealing how ambient micro-turbulence can be leveraged to maintain a macro-stable, high-performance pedestal for future fusion reactors.

experimental nuclear physics↗