SEARCH · Engineering Papers
Results for “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.
The virtual element method for a 2D incompressible MHD system
Not Available
MHD stability constraints on divertor heat flux width in DIII-D
Not Available
MHD stability and disruptions in the SPARC tokamak
SPARC is being designed to operate with a normalized beta of $β_N$ = 1.0, a normalized density of $n_G$ = 0.37 and a safety factor of $q_{95}$ ≈ 3.4, providing a comfortable margin to their respective disruption limits. Further, a low beta poloidal $β_p$ = 0.19 at the safety factor $\textit{q}$ = 2 surface reduces the drive for neoclassical tearing modes, which together with a frozen-in classically stable current profile might allow access to a robustly tearing-free operating space. Although the inherent stability is expected to reduce the frequency of disruptions, the disruption loading is comparable to and in some cases higher than that of ITER. The machine is being designed to withstand the predicted unmitigated axisymmetric halo current forces up to 50 MN and similarly large loads from eddy currents forced to flow poloidally in the vacuum vessel. Runaway electron (RE) simulations using GO+CODE show high flattop-to-RE current conversions in the absence of seed losses, although NIMROD modelling predicts losses of ~80 %; self-consistent modelling is ongoing. A passive RE mitigation coil designed to drive stochastic RE losses is being considered and COMSOL modelling predicts peak normalized fields at the plasma of order 10 –2 that rises linearly with a change in the plasma current. Massive material injection is planned to reduce the disruption loading. A data-driven approach to predict an oncoming disruption and trigger mitigation is discussed.
An adjoint-based method for optimising MHD equilibria against the infinite- n , ideal ballooning mode
We demonstrate a fast adjoint-based method to optimise tokamak and stellarator equilibria against a pressure-driven instability known as the infinite-n ideal ballooning mode. We present three finite-β (the ratio of thermal to magnetic pressure) equilibria: one tokamak equilibrium and two stellarator equilibria that are unstable against the ballooning mode. Using the self-adjoint property of ideal magnetohydrodynamics, we construct a technique to rapidly calculate the change in the eigenvalue, a measure of ideal ballooning instability. Using the SIMSOPT optimisation framework, we then implement our fast adjoint gradient-based optimiser to minimise the eigenvalue and find stable equilibria for each of the three originally unstable equilibria.
Exploring Localized Geomagnetic Disturbances in Global MHD: Physics and Numerics
One of the prominent effects of space weather is the formation of rapid geomagnetic field variations on Earth's surface driven by the magnetosphere-ionosphere system. These geomagnetic disturbances (GMDs) cause geomagnetically induced currents to run through ground conducting systems. In particular, localized GMDs (LGMDs) can be high amplitude and can have an effect on scale sizes less than 100 km, making them hazardous to power grids and difficult to predict. In this study, we examine the ability of the Space Weather Modeling Framework (SWMF) to reproduce LGMDs in the 7 September 2017 event using both existing and new metrics to quantify the success of the model against observation. We show that the high-resolution SWMF can reproduce LGMDs driven by ionospheric sources, but struggles to reproduce LGMDs driven by substorm effects. We calculate the global maxima of the magnetic fluctuations to show instances when the SWMF captures LGMDs at the correct times but not the correct locations. To remedy these shortcomings we suggest model developments that will directly impact the ability of the SWMF to reproduce LGMDs, most importantly updating the ionospheric conductance calculation from empirical to physics-based.
Planck intermediate results. XX. Comparison of polarized thermal emission from Galactic dust with simulations of MHD turbulence
Polarized emission observed by Planck HFI at 353 GHz towards a sample of nearby fields is presented, focusing on the statistics of polarization fractions p and angles ψ. The polarization fractions and column densities in these nearby fields are representative of the range of values obtained over the whole sky. We find that: (i) the largest polarization fractions are reached in the most diffuse fields; (ii) the maximum polarization fraction p max decreases with column density N H in the more opaque fields with N H > 10 21 cm -2 ; and (iii) the polarization fraction along a given line of sight is correlated with the local spatial coherence of the polarization angle. These observations are compared to polarized emission maps computed in simulations of anisotropic magnetohydrodynamical turbulence in which we assume a uniform intrinsic polarization fraction of the dust grains. We find that an estimate of this parameter may be recovered from the maximum polarization fraction p max in diffuse regions where the magnetic field is ordered on large scales and perpendicular to the line of sight. This emphasizes the impact of anisotropies of the magnetic field on the emerging polarization signal. The decrease of the maximum polarization fraction with column density in nearby molecular clouds is well reproduced in the simulations, indicating that it is essentially due to the turbulent structure of the magnetic field: an accumulation of variously polarized structures along the line of sight leads to such an anti-correlation. In the simulations, polarization fractions are also found to anti-correlate with the angle dispersion function S. However, the dispersion of the polarization angle for a given polarization fraction is found to be larger in the simulations than in the observations, suggesting a shortcoming in the physical content of these numerical models. In summary, we find that the turbulent structure of the magnetic field is able to reproduce the main statistical properties of the dust polarization as observed in a variety of nearby clouds, dense cores excluded, and that the large-scale field orientation with respect to the line of sight plays a major role in the quantitative analysis of these statistical properties.
Evidence for a resonant cyclotron interaction between runaway electrons and MHD modes in the experimental advanced superconducting tokamak
In the past, the resonant cyclotron interaction between runaway electrons and lower hybrid waves via anomalous Doppler broadening was experimentally investigated, and it was shown to be able to create a barrier to the energy that could be reached by the runaway electrons [E. Li et al., Nucl. Instrum. Methods Phys. Res. A 621, 566 (2010)]. In this paper, to our knowledge for the first time, experimental evidence will be provided for a resonant cyclotron interaction between runaway electrons and magnetohydrodynamics modes in a stochastic magnetic field in the experimental advanced superconducting tokamak (EAST), which has been theoretically proposed as a mechanism able to limit the maximum attainable energy by runaway electrons in tokamak plasmas [J. R. Martin-Solis and R. Sanchez, Phys. Plasmas 15, 112505 (2008)].
MHD flow of a micropolar fluid over a stretchable disk in a porous medium with heat and mass transfer
This article studies the simultaneous impacts of heat and mass transfer of an incompressible electrically conducting micropolar fluid generated by the stretchable disk in presence of porous medium. The thermal radiation effect is accounted via Rosseland’s approximation. The governing boundary layer equations are reduced into dimensionless form by employing the suitable similarity transformations. A finite difference base algorithm is utilized to obtain the solution expressions. The impacts of physical parameters on dimensionless axial velocity, radial velocity, micro-rotation, temperature and concentrations profiles are presented and examined carefully. Numerical computation is performed to compute shear stress, couple stress, heat and mass rate at the disk.
The direct criterion of Newcomb for the ideal MHD stability of an axisymmetric toroidal plasma
A method is presented for determining the ideal magnetohydrodynamic stability of an axisymmetric toroidal plasma, based on a toroidal generalization of the method developed by Newcomb for fixed-boundary modes in a cylindrical plasma. For toroidal mode number n≠0, the stability problem is reduced to the numerical integration of a high-order complex system of ordinary differential equations, the Euler-Lagrange equation for extremizing the potential energy, for the coupled amplitudes of poloidal harmonics m as a function of the radial coordinate ψ in a straight-fieldline flux coordinate system. Unlike the cylindrical case, different poloidal harmonics couple to each other, which introduces coupling between adjacent singular intervals. A boundary condition is used at each singular surface, where m = nq and q(ψ) is the safety factor, to cross the singular surface and continue the solutions beyond it. Fixed-boundary instability is indicated by the vanishing of a real determinant of a Hermitian complex matrix constructed from the fundamental matrix of solutions, the generalization of Newcomb's crossing criterion. In the absence of fixed-boundary instabilities, an M × M plasma response matrix W P with M the number of poloidal harmonics used, is constructed from the Euler-Lagrange solutions at the plasma-vacuum boundary. This is added to a vacuum response matrix W V to form a total response matrix W T . Finally, the existence of negative eigenvalues of W T indicates the presence of free-boundary instabilities. The method is implemented in the fast and accurate DCON code.
Study the effect of chemical reaction and variable viscosity on free convection MHD radiating flow over an inclined plate bounded by porous medium
An analysis is performed to study the free convection heat and mass transfer flow of an electrically conducting incompressible viscous fluid about a semi-infinite inclined porous plate under the action of radiation, chemical reaction in presence of magnetic field with variable viscosity. The dimensionless governing equations are steady, two-dimensional coupled and non-linear ordinary differential equation. Nachtsgeim-Swigert shooting iteration technique along with Runge-Kutta integration scheme is used to solve the non-dimensional governing equations. The effects of magnetic parameter, viscosity parameter and chemical reaction parameter on velocity, temperature and concentration profiles are discussed numerically and shown graphically. Therefore, the results of velocity profile decreases for increasing values of magnetic parameter and viscosity parameter but there is no effect for reaction parameter. The temperature profile decreases in presence of magnetic parameter, viscosity parameter and Prandtl number but increases for radiation parameter. Also, concentration profile decreases for the increasing values of magnetic parameter, viscosity parameter and reaction parameter. All numerical calculations are done with respect to salt water and fixed angle of inclination of the plate.
A conservative scheme of drift kinetic electrons for gyrokinetic simulation of kinetic-MHD processes in toroidal plasmas
A conservative scheme of drift kinetic electrons for gyrokinetic simulations of kinetic-magnetohydrodynamic processes in toroidal plasmas has been formulated and verified. Both vector potential and electron perturbed distribution function are decomposed into adiabatic part with analytic solution and non-adiabatic part solved numerically. The adiabatic parallel electric field is solved directly from the electron adiabatic response, resulting in a high degree of accuracy. Here, the consistency between electrostatic potential and parallel vector potential is enforced by using the electron continuity equation. Since particles are only used to calculate the non-adiabatic response, which is used to calculate the non-adiabatic vector potential through Ohm's law, the conservative scheme minimizes the electron particle noise and mitigates the cancellation problem. Linear dispersion relations of the kinetic Alfvén wave and the collisionless tearing mode in cylindrical geometry have been verified in gyrokinetic toroidal code simulations, which show that the perpendicular grid size can be larger than the electron collisionless skin depth when the mode wavelength is longer than the electron skin depth.
A robust solution for the resistive MHD toroidal Δ′ matrix in near real-time
Not Available
MHD stability trends and improved performance of LHD inward-shifted configurations: The role of the neutral beam current drive and thermal plasma density
The aim of the present study is to analyze the effect of the neutral beam current drive (NBCD), thermal plasma density, and NBI operational regime on the stability of pressure gradient-driven modes (PGDM) and Alfvén eigenmodes (AE) in LHD inward-shifted configurations. The stabilization of n/m = 1/2 PGDM (n toroidal mode and m poloidal mode) is observed in the discharge 167 800 during the co-NBCD phase. The iota profile evolution measured by motional stark effect diagnostic may indicate the iota profile up-shift caused by the co-NBCD can induce a non-resonant transition of the rational surface 1/2 before the mode stabilization. The evolution of the iota profile and continuum gaps in the discharge 167 805 during the ctr-NBCD phase leads to the stabilization of the AE, caused by the narrowing of the continuum gap as the iota profile down-shift. Opposite stability trends are identified for PGDM and AE stability with respect to the thermal plasma density. A larger thermal plasma density (larger thermal β) further enhances PGDM although the continuum gaps are narrower leading to configurations with stable AEs. The linear stability of AEs is analyzed using the gyro-fluid FAR3d code to reproduce the AE stability trends observed in the experiments with respect to the NBCD and thermal plasma density. The analysis of hypothetical scenarios dedicated to study different NBI operational regimes with respect to EP energy, and β and radial density profiles indicate off-axis NBI operation shows a higher EP β threshold to destabilize AEs compared to on-axis configuration. This is explained by the presence of a TAE gap in the inner plasma region, easily destabilized by an on-axis NBI injection. The control of the NBCD and thermal plasma in the discharge 167 800 shows a transitory stabilization of PGDM and AEs, as well as an improved discharge performance identified by an increment of the neutron fluxes.
Kelvin–Helmholtz instability under stabilizing parallel magnetic field in nonhomogeneous compressible MHD flows
We study the Kelvin–Helmholtz instability (KHI) for the general case of a compressible, nonhomogeneous, magnetized plasma flow. The study is limited to a vortex sheet interface with an imposed parallel magnetic field. We introduce a new formalism based on a convective Mach number M c , a convective Alfvénic Mach number M Ac , and a total convective Mach number that combines the two. We derive an analytic expression of the KHI growth rate for a homogeneous flow (i.e., zero Atwood number, A=0) that converges toward both the expression for unmagnetized compressible flow and Chandrasekhar's expression for magnetized incompressible flow. Otherwise, the dispersion relation is solved numerically and allows deriving general stability diagrams of magnetized KHI for the triplet (A, M c , β −plasma) parameters. We show these parameters uniquely define all configurations for a parallel magnetic field. We also construct diagrams with respect to the convective Alfvénic Mach number, the β − plasma parameter, or the magnetic field showing which magnetic field strength is required for stabilizing a given shear flow. The theoretical growth rates are compared with 18 simulations made with the GAMERA code, currently used for 3D magnetospheric simulations. Finally, we apply our results to the analysis of a past KHI experiment performed at the OMEGA laser facility, showing linear theory succeeds to provide accurate estimates of the growth rate at early times. We further discuss how our results can inform future experiments in the high-Mach magnetized regime at the National Ignition Facility. Possible limitations of the study due to resistive, mixing, or turbulence effects are discussed.
Orbit-averaging and deposition accuracy for runaway electron beams in hybrid kinetic-MHD simulations of the runaway plateau
We develop a new procedure that combines the kinetic orbit runaway electron code (KORC) and the NIMROD extended-magnetohydrodynamic code to simulate runaway electrons (REs) in the post-disruption plateau. KORC integrates guiding-center orbits, with a barycentric-based binary search strategy providing initial guesses for the Newton–Raphson logical-to-physical coordinate inversion, ensuring reliable particle-to-mesh mapping in NIMROD, whose fields remain static for the present study. Samples are drawn in accord with experimental parallel current profiles of RE beams during the plateau phase. Deposition in NIMROD is verified through comparison with a Python-based finite-element code that ensures periodicity in the poloidal direction and continuity at the magnetic axis. Accurate representation of near-axis fields requires finer mesh resolution to prevent under- and overshoots in current density from orbit inaccuracies. Yet, at a fixed particle count, increasing mesh resolution amplifies statistical noise in the deposited fields. An orbit-averaging method accumulates partial current deposits over multiple kinetic steps and reduces the statistical noise with little added computational cost. By coupling kinetic routines from KORC directly into the NIMROD codebase, these developments lay essential groundwork for future self-consistent KORC–NIMROD coupling.