Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

At least 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The influence of boundary and edge-plasma modeling in computations of axisymmetric vertical displacement

A number of previously published studies apply extended-MHD computations to model macroscopic dynamics of tokamak vertical displacement events (VDEs). The characteristic plasma-wall contact and resistive-wall diffusion imply sensitivity to boundary conditions in a general sense, but which conditions have significant influence depends on how the plasma is modeled. This work examines these dependencies by comparing results of axisymmetric extended-MHD computations with different sets of boundary conditions, plasma thermal-conduction models, and wall resistivity values. The geometry and plasma parameters of the computations represent a modest-sized tokamak. The forced-displacement scenario and computational setup are similar to those used in previous 3D computations. The results show that for a given wall resistivity, the VDE time-scale is most sensitive to variations in the boundary and thermal-conduction parameters that affect electron thermal transport. The electrical conductivity depends on electron temperature, and the dependence on thermal transport stems from its influence on the electrical circuit that includes the open-field halo current. Lastly, conditions that lead to hotter, broader halo regions slow the evolution. Significant sensitivity to the boundary condition on plasma flow-velocity exists when electron thermal conduction is restricted and electron energy loss is convective, which is expected for conditions at the entrance of the magnetic presheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Gyrokinetic analysis of inter-edge localized mode transport mechanisms in a DIII-D pedestal

In this study, gyrokinetic simulations are used to study pedestal fluctuations for DIII-D discharge 174082 using the GENE code. Nonlinear local simulations indicate that electron heat flux has contributions from electron temperature gradient-driven transport but at levels insufficient to satisfy power balance. We show that microtearing modes (MTM) and neoclassical transport are likely to account for the remaining observed energy losses in the electron and ion channels, respectively. The MTM instabilities found in the simulations are consistent with the high-frequency fluctuations identified in the magnetic fluctuation data from Mirnov coils. The fluctuation data in this discharge also exhibit a low-frequency band of fluctuations. By modifying the equilibrium profiles and plasma β, simulations produce MHD modes, which may be responsible for these observed low-frequency fluctuations. We compare several metrics involving ratios of fluctuation amplitudes and transport quantities for both MTMs and MHD modes. This analysis suggests that the available data are consistent with the simultaneous activity of both MHD modes and MTMs provided that the former is limited largely to the particle transport channel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Pressure-driven tearing and thermal transport in finite-beta reversed field pinch computations

In this work, nonlinear resistive-magnetohydrodynamics (MHD) computation with heating and anisotropic transport is applied to examine the interaction between thermal energy and magnetic fluctuations in inductively driven reversed-field pinches (RFPs). The magnetic fluctuations underlie magnetic field reversal through dynamo-like correlations, and they enhance thermal energy transport through fluctuations of parallel heat flux density. With the unfavorable magnetic curvature that exists across the RFP profile, thermal energy also affects the magnetic fluctuations. Computations with the NIMROD code [Sovinec et al., J. Comput. Phys. 195, 355–386 (2004)] integrate nonlinear MHD dynamics with energy transport and reproduce an RFP state with experimentally relevant values of plasma-β. Equilibria constructed from results of the 3D computations are analyzed to assess the sources of free energy in the saturated nonlinear state. Linear computations for these profiles show unstable modes of tearing parity. Their eigenfunctions are used to evaluate and compare stabilizing and destabilizing contributions to the kinetic energy integral. An assessment of the drives in the integral reveals that the pressure gradient drive is of comparable magnitude to the parallel current drive, and only the sum of the two surpasses the stabilizing contributions. Correlation of magnetic and parallel heat flux density fluctuations in the nonlinear computations shows that fluctuation-induced thermal conduction is the dominant mode of energy loss, as expected from experimental evidence. Decomposition of the fluctuating heat flux density shows that second-order correlations, alone, do not explain the total energy transport. Higher-order correlations are also important.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

High resolution density pedestal measurements during edge localized modes by short-pulse reflectometry in the TCV tokamak

Here, this publication presents high spatio-temporal resolution (mm µs -1 ) density profile measurements of the pedestal top during type I, III, and small edge localized mode (ELM) H-mode plasmas in the tokamak à configuration variable (TCV). These measurements were performed using a novel short-pulse reflectometer. Average inter-ELM density profiles are obtained via conditional averaging using the D α trace as ELM indicator. Changes to the pedestal density profile gradients prior to type-III ELMs reveal unique pedestal dynamics leading to the ELM crash which can provide important experimental data for validation of non-linear MHD ELM simulations. The small-ELM scenario is found to feature a ~25–35 kHz quasi-coherent density fluctuation near the separatrix ρ ψ 0.993–1.05 not observed during a similar type-I ELM discharge. This oscillation is also found in low-field-side magnetic pick-up probes displaying a ballooning character and n = +1 toroidal mode number. This oscillation could help explain the markedly different pedestal dynamics observed in the small-ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Impact of impurities on peeling–ballooning modes and turbulence in tokamak plasmas

This study investigates the impact of various impurity species on peeling–ballooning (PB) modes and microturbulence in tokamak plasmas through the extension of traditional two-fluid and gyro-landau-fluid (GLF) models. By incorporating finite Larmor radius (FLR) effects, the analysis provides a comprehensive understanding of impurity-driven impact and its interaction with plasma turbulence. Depending on charge state and local plasma conditions, heavy impurities may exhibit gyro-radii larger than those of main ions, which are captured in the extended GLF model presented. Following the presentation of modified two-fluid equations incorporating impurity effects, we systematically analyze the distinctions between impurity and main ion dynamics and their resultant feedback mechanisms on plasma behavior. Derivation of the linear dispersion relation enables quantification of impurity-mediated modifications to: plasma vorticity, diamagnetic drift and gyroviscous effects, electron Hall physics, and FLR effects. BOUT++ – based linear simulations corroborate this formalism, demonstrating systematic stabilization of PB modes upon impurity seeding. And then operational implications for practical impurity control strategies in tokamak devices are proposed. The results underscore the necessity of impurity management to maintain stability and optimize plasma confinement, with specific focus on how FLR effects contribute to transport dynamics. This work paves the way for enhanced modeling and simulation efforts, supporting the development of strategies to control impurity-induced turbulence and improve overall reactor performance.

BOUT++ simulation↗
DOE OSTI2022

Frequency multiplication with toroidal mode number of kink/fishbone modes on a static HL-2A-like tokamak

In the presence of energetic particles (EPs), the long-lived mode (LLM) frequency multiplication with n = 1, 2, 3, or higher is often observed on HL-2A, where n is the toroidal mode number. Hybrid kinetic-MHD model simulations of the energetic particle (EP) driven kink/fishbone modes on a static HL-2A-like tokamak using NIMROD code find that when the background plasma pressure is relatively high, and the EP pressure and the beam energy are relatively low, the mode frequency increases almost linearly with EP pressure, and the frequency is proportional to n (‘frequency multiplication’), even in the absence of any equilibrium plasma rotation. In addition, the frequency multiplication persists as the safety factor at the magnetic axis q 0 varies. In the absence of EPs, the growth rate of the 1/1 mode is the largest; however, as the EP pressure increases, the growth rate of 2/2 modes or 3/3 modes becomes dominant, suggesting that higher-n modes are more vulnerable to EPs. These results may shed light on the understanding of the toroidal mode number dependence of kink/fishbone modes in the advanced scenarios of tokamaks with weak or reversed central magnetic shear.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Kink-driven magnetic reconnection in relativistic jets: consequences for X-ray polarimetry of BL Lacs

ABSTRACT We investigate with relativistic MHD (magnetohydrodynamic) simulations the dissipation physics of BL Lac jets, by studying the synchrotron polarization signatures of particles accelerated by the kink instability in a magnetically dominated plasma column. The non-linear stage of the kink instability generates current sheets, where particles can be efficiently accelerated via magnetic reconnection. We identify current sheets as regions where s = Jδ/B is above some predefined threshold (where B is the field strength, J is the current density, and δ is the grid scale), and assume that the particle injection efficiency scales as ∝J2. X-ray emitting particles have short cooling times, so they only probe the field geometry of their injection sites. In contrast, particles emitting in the optical band, which we follow self-consistently as they propagate away from their injection sites while cooling, sample a larger volume, and so they may be expected to produce different polarimetric signatures. We find that the degree of polarization is roughly the same between X-ray and optical bands, because even the optical emitting particles do not travel far from the current sheet where they were injected, due to lack of sufficient kink-generated turbulence. The polarization angle shows a different temporal evolution between the two bands, due to the different regions probed by X-ray and optical emitting particles. In view of the upcoming IXPE satellite, our results can help constrain whether kink-induced reconnection (as opposed to shocks) can be the source of multiwavelength emission from BL Lacs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Modeling of resistive plasma response in toroidal geometry using an asymptotic matching approach

The method of solving linear resistive plasma response, based on the asymptotic matching approach, is developed for full toroidal tokamaks by upgrading the Resistive DCON code [A.H. Glasser, Z.R. Wang and J.-K. Park, Physics of Plasmas, \textbf{23}, 112506 (2016)]. The derived matching matrix, asymptotically matching the outer and inner regions, indicates that the applied three dimension (3-D) magnetic perturbations contribute additional small solutions at each resonant surface due to the toroidal coupling of poloidal modes. In contrast, the resonant harmonic only affects the corresponding resonant surface in the cylindrical plasma. Since the solution of ideal outer region is critical to the asymptotic matching and is challenging to be solved in the toroidal geometry due to the singular power series solution at the resonant surfaces, systematic verification of the outer region $\Delta^\prime$ matrix is made by reproducing the well known analytical $\Delta^{\prime}$ result in [H.P. Furth, P.H. Rutherford and H. Selberg, The Physics of Fluids, \textbf{16}, 1054-1063 (1073)] as well as by making a quantitative benchmark with the PEST3 code [A. Pletzer and R.L. Dewar, J. Plasma Physics, \textbf{45}, 427-451 (1991)]. Finally, the reconstructed numerical solution of resistive plasma response from the toroidal matching matrix is presented. Comparing with the ideal plasma response, the global structure of the response can be affected by the small finite island at the resonant surfaces.

MHD simulation↗
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 OSTI2023

Setting up resistive wall modeling for the GA Compact Advanced Tokamak (GA FPP) (Final Report)

Resistive wall implementation for tokamaks has been achieved by both a 'thin wall' approximation in which time-varying fields are matched across a boundary for all modes; and, by meshing the region from the first wall to the outer boundary and setting certain fields to mimic a resistive wall. Both implementations were tested on existing tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Magnetic Field Strength from Turbulence Theory. I. Using Differential Measure Approach

The mean plane-of-sky magnetic field strength is traditionally obtained from the combination of polarization and spectroscopic data using the Davis–Chandrasekhar–Fermi (DCF) technique. However, we identify the major problem of the DCF technique to be its disregard of the anisotropic character of MHD turbulence. On the basis of the modern MHD turbulence theory we introduce a new way of obtaining magnetic field strength from observations. Unlike the DCF technique, the new technique uses not the dispersion of the polarization angle and line-of-sight velocities, but increments of these quantities given by the structure functions. To address the variety of astrophysical conditions for which our technique can be applied, we consider turbulence in both media with magnetic pressure higher than the gas pressure, corresponding, e.g., to molecular clouds, and media with gas pressure higher than the magnetic pressure, corresponding to the warm neutral medium. We provide general expressions for arbitrary admixtures of Alfvén, slow, and fast modes in these media and consider in detail particular cases relevant to diffuse media and molecular clouds. We successfully test our results using synthetic observations obtained from MHD turbulence simulations. We demonstrate that our differential measure approach, unlike the DCF technique, can be used to measure the distribution of magnetic field strengths, can provide magnetic field measurements with limited data, and is much more stable in the presence of induced large-scale variations of nonturbulent nature. Furthermore, our study uncovers the deficiencies of earlier DCF research.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The Role of Magnetic Reconnection in Energizing Protons and Heavier Ions at the Heliospheric Current Sheet

During near-Sun crossings of the heliospheric current sheet (HCS), Parker Solar Probe (PSP) observed populations of high-energy protons and heavier ions, indicating possible energization by magnetic reconnection up to 10 s–100 s keV nucleon −1 . Here we study ion acceleration by magnetic reconnection at the HCS. To estimate ion energization, we solve the Parker transport equation coupled to a large-scale 2D MHD reconnection simulation. We find that multiple ion species develop power-law distributions with both spectral index and high-energy cutoff E max consistent with in situ data. By accounting for the injection physics determined by kinetic simulations, we confirm that the charge-to-mass ratio scales as E max ∝ (Q/M) α with α ∼ 0.8–1.1, approximately consistent with PSP measurements in the broader range α ∼ 0.6–1.7. In the limit where ions are injected at the same energy per nucleon, α can be as low as ∼0.3. These findings further support the role of magnetic reconnection in producing high-energy heavy ions at the HCS.

Murtas, Giulia [West Virginia Univ., Morgantown, W↗
DOE OSTI2022

High resolution density pedestal measurements during edge localized modes by short-pulse reflectometry in the TCV tokamak

This publication presents high spatio-temporal resolution (mm/μs) density profile measurements of the pedestal top during type I, III, and small edge localized mode (ELM) H-mode plasmas in the Tokamak à Configuration Variable (TCV). These measurements were performed using a novel short-pulse reflectometer. Average inter-ELM density profiles are obtained via conditional averaging using the Dα trace as ELM indicator. Changes to the pedestal density profile gradients prior to type-III ELMs reveal unique pedestal dynamics leading to the ELM crash which can provide important experimental data for validation of non-linear MHD ELM simulations. The small-ELM scenario is found to feature a ∼25-35 kHz quasi-coherent density fluctuation near the separatrix rho_psi ∼0.993-1.05 not observed during a similar type-I ELM discharge. This oscillation is also found in low-field-side magnetic pick- up probes displaying a ballooning character and n=+1 toroidal mode number. This oscillation could help explain the markedly different pedestal dynamics observed in the small-ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

On the formation of coronal cavities

The formation of a coronal cavity and its relation to a quiescent prominence is studied theoretically. The stability of condensation modes of a plasma in the coronal streamer model (Steinolfson et al., 1982) is considered using a two-dimensional time-dependent ideal MHD numerical simulation. It is found that a plasma with beta = 0.5 is unstable but one with beta = 4 is stable because the density enhancement of the plasma trapped by the closed fields increases with the strength of the magnetic field. The means by which condensation modes can produce a coronal cavity and/or initiate the formation of a prominence (depending on the field configuration) are discussed. It is argued that prominence and cavity material is all supplied from the chromospheric level in the form of spicules.

An, C.-H.↗
NASA NTRS1992

Constraints on cosmic ray propagation in the galaxy

The goal was to derive a more detailed picture of magnetohydrodynamic turbulence in the interstellar medium and its effects on cosmic ray propagation. To do so, radio astronomical observations (scattering and Faraday rotation) were combined with knowledge of solar system spacecraft observations of MHD turbulence, simulations of wave propagation, and modeling of the galactic distribution to improve the knowledge. A more sophisticated model was developed for the galactic distribution of electron density turbulence. Faraday rotation measure data was analyzed to constrain magnetic field fluctuations in the ISM. VLBI observations were acquired of compact sources behind the supernova remnant CTA1. Simple calculations were made about the energies of the turbulence assuming a direct link between electron density and magnetic field variations. A simulation is outlined of cosmic ray propagation through the galaxy using the above results.

Cordes, James M.↗
NASA NTRS1993

Magnetospheric response to solar wind dynamic pressure variations: Interaction of interplanetary tangential discontinuities with the bow shock

Some magnetic impulse events observed in the polar region are related to vortices associated with plasma convection in the ionosphere. Recent analyses of satellite and ground data suggest that the interaction of solar wind dynamic pressure pulses and the magnetosphere may lead to the formation of velocity vortices in the magnetopause boundary layer region. This can in turn lead to the presence of vortices in the polar ionosphere. However, before reaching the Earth's magnetopause, these interplanetary pressure pulses must interact with and pass through the bow shock. A variation of the solar wind dynamic pressure may be associated with shocks, magnetic holes, or tangential discontinuities (TDs) in the interplanetary medium. We study the interaction of interplanetary TDs with the Earth's bow shock (BS) using both theoretical analysis and MHD computer simulations. It is found that as a result of the collision between a TD and the BS, the jump in the solar wind dynamic pressure associated with the TD is significantly modified, the bow shock moves, and a new fast shock or fast rarefaction wave, which propagates in the downstream direction, is excited. Our theoretical analysis shows that the change in the plasma density across the interplanetary TD plays the most important role in the collision process. In the case with an enhanced dynamic pressure behind the interplanetary TD, the bow shock is intensified in strength and moves in the earthward direction. The dynamic pressure jump associated with the transmitted TD is generally reduced from the value before the interaction. A fast compressional shock is excited ahead of the transmitted TD and propagates toward the Earth's magnetosphere. For the case in which the dynamic pressure is reduced behind the interplanetary TD, the pressure jump across the transmitted TD is substantially weakened, the bow shock moves in the sunward direction, and a rarefaction wave which propagates downstream is excited. We also simulate and discuss the interaction of a pair of tangential discontinuities, which may correspond to a magnetic hole, with the BS.

Wu, Bor-Han↗
NASA NTRS1995

Origin of coronal mass ejection and magnetic cloud: Thermal or magnetic driven?

A fundamental problem in Solar-Terrestrial Physics is the origin of the solar transient plasma output, which includes the coronal mass ejection and its interplanetary manifestation, e.g. the magnetic cloud. The traditional blast wave model resulted from solar thermal pressure impulse has faced with challenge during recent years. In the MHD numerical simulation study of CME, the authors find that the basic feature of the asymmetrical event on 18 August 1980 can be reproduced neither by a thermal pressure nor by a speed increment. Also, the thermal pressure model fails in simulating the interplanetary structure with low thermal pressure and strong magnetic field strength, representative of a typical magnetic cloud. Instead, the numerical simulation results are in favor of the magnetic field expansion as the likely mechanism for both the asymmetrical CME event and magnetic cloud.

Zhang, Gong-Liang↗