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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2013

MHD Simulations of the Plasma Flow in the Magnetic Nozzle

The magnetohydrodynamic (MHD) flow of plasma through a magnetic nozzle is simulated by solving the governing equations for the plasma flow in the presence of an static magnetic field representing the applied nozzle. This work will numerically investigate the flow and behavior of the plasma as the inlet plasma conditions and magnetic nozzle field strength are varied. The MHD simulations are useful for addressing issues such as plasma detachment and to can be used to gain insight into the physical processes present in plasma flows found in thrusters that use magnetic nozzles. In the model, the MHD equations for a plasma, with separate temperatures calculated for the electrons and ions, are integrated over a finite cell volume with flux through each face computed for each of the conserved variables (mass, momentum, magnetic flux, energy) [1]. Stokes theorem is used to convert the area integrals over the faces of each cell into line integrals around the boundaries of each face. The state of the plasma is described using models of the ionization level, ratio of specific heats, thermal conductivity, and plasma resistivity. Anisotropies in current conduction due to Hall effect are included, and the system is closed using a real-gas equation of state to describe the relationship between the plasma density, temperature, and pressure.A separate magnetostatic solver is used to calculate the applied magnetic field, which is assumed constant for these calculations. The total magnetic field is obtained through superposition of the solution for the applied magnetic field and the self-consistently computed induced magnetic fields that arise as the flowing plasma reacts to the presence of the applied field. A solution for the applied magnetic field is represented in Fig. 1 (from Ref. [2]), exhibiting the classic converging-diverging field pattern. Previous research was able to demonstrate effects such as back-emf at a super-Alfvenic flow, which significantly alters the shape of the magnetic field in both the near- and far-field regions. However, in that work the downstream domain was constrained to a channel of constant cross-sectional area. In the present work we seek to address this issue by modeling the downstream region with a domain that permits free expansion of the plasma, permitting a better evaluation of the downstream effects the applied field has on the plasma. The inlet boundary conditions and applied magnetic field values will also be varied to determine the effect the initial plasma energy content and applied magnetic field energy density have on the near- and far-field plasma properties on the MHD code. This will determine the effect of inlet boundary conditions on the results downstream and address issues related to the restrictive numerical domain previously used.

Smith, T. E. R.↗
DOE OSTI2026

Linear simulations of a wide pedestal quiescent H-mode plasma with the extended-MHD code NIMROD

We present linear simulations of a Wide Pedestal QH (WPQH)-mode DIII-D plasma with the NIMROD code to address the role of ion multispecies collisionality on WPQH-modes. We use a full extended-MHD model that captures Alfvénic and magnetoacoustic waves, including two-fluid, ion gyroviscosity and cross heat flux. We show that two-fluid and ion gyroviscous effects destabilize electron-directed peeling–ballooning modes in the pedestal. The linear growth rates of the modes decrease when carbon (C) is added to a deuterium plasma because the multispecies collisional effects increase the plasma resistivity. When replacing C with tungsten (W), the multispecies collisional effects increase the resistivity further, as well as the stabilizing effect on the instabilities. The resistivity impact on mode stability is reversed when both two-fluid and ion gyroviscous effects are ignored, which is consistent with previous works, proving that multispecies collisional effects, together with two-fluid and ion gyroviscosity, should be considered to address the stability of plasmas with significant impurity content.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

Numerical analysis method for linear induction machines.

A numerical analysis method has been developed for linear induction machines such as liquid metal MHD pumps and generators and linear motors. Arbitrary phase currents or voltages can be specified and the moving conductor can have arbitrary velocity and conductivity variations from point to point. The moving conductor is divided into a mesh and coefficients are calculated for the voltage induced at each mesh point by unit current at every other mesh point. Combining the coefficients with the mesh resistances yields a set of simultaneous equations which are solved for the unknown currents.

Elliott, D. G.↗
DOE OSTI2024

Boundary condition effects on runaway electron mitigation coil modeling for the SPARC and DIII-D tokamaks

Abstract Extended-MHD modeling of planned Runaway Electron Mitigation Coils (REMC) for SPARC and DIII-D is performed with the NIMROD code. A coil has been designed for each machine, with the two differing in shape and location, but both having n = 1 symmetry (with n the toroidal mode number). Compared to previous modeling efforts, three improvements are made to the simulations boundary conditions. First a resistive wall model is used in place of an ideal wall. Second, the ThinCurr code is used to compute the time-dependent 3D fields used as magnetic boundary conditions for the simulations. Third, the simulation boundary is moved from the first-wall location to the Vacuum Vessel (VV), which extends the boundary past the location of the internal REMC. To remove the 3D coil from the simulation domain, an equivalent set of 3D fields is calculated at the VV boundary that produce approximately the same field distribution at the last closed flux surface assuming vacuum between the two. Each of these three boundary condition improvements leads to an improvement in the predicted performance of the REMC for both machines. The resistive wall alone primarily effects the resonance of the coil with the plasma after the TQ, affecting the q-profile evolution in the SPARC modeling, and allowing the applied spectrum to be modified in response to the plasma in the DIII-D modeling. The movement of the simulation boundary has the most significant effect on the RE confinement overall, including in the early stages, particularly for a DIII-D inner wall limited equilibrium, where the RE loss fraction increases from 90% to > 99%, with SPARC RE losses also occurring much earlier when the boundary is placed at the VV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Viscous normal modes on coronal inhomogeneities and their role as a heating mechanism

Viscous damping of Alfven surface waves is examined both analytically and numerically using incompressible MHD. Normal modes are shown to exist on discontinuous as well as continuously varying interfaces in Alfven speed. The waves experience negligible decay below the transition zone. High-frequency waves damp just above the transition region, while those of lower frequency lose energy further out. A comparison of dissipative decay rates shows that wave damping by viscosity proceeds approximately two orders of magnitude faster than by resistivity.

Steinolfson, R. S.↗
NASA NTRS1979

Global magnetohydrodynamic simulation of the two-dimensional magnetosphere

The time-dependent magnetohydrodynamic interaction of the solar wind with a two-dimensional dipole magnetic field has been simulated using a novel Lagrangian particle type of MHD code that can treat local low density or vacuum regions without numerical instability. This enables one to simulate the time-dependent magnetic tail. When the solar wind field is southward, a magnetic field line topology consistent with Dungey's model emerges in steady state. The tail, however, is short, and the x-points are only slightly shifted from their vacuum locations, because of strong numerical resistivity. Different configurations resulting from different relative orientations of the solar wind magnetic field and dipole axis are also presented. While the magnetic field is relatively steady, the density and flow in the magnetosheath are turbulent, as are the bow shock and magnetopause; the Kelvin-Helmholtz instability may account for these phenomena. We also model a 'substorm' as the passage of a rotational discontinuity in the solar wind over the dipole. Both 90 and 180 deg shifts to a southward solar wind field cause a violent readjustment of the magnetic tail which eventually settles down to the Dungey configuration.

Leboeuf, J. N.↗
DOE OSTI2023

Stability and Control of Burning Tokamak Plasmas with Resistive Walls (Final Report)

This research has focused on quantitative prediction of the stability, control, and equilibrium state evolution in toroidal burning plasmas. The stability of long pulse burning plasmas is highly sensitive to the physics of resonant layers in the plasma, sources of momentum and flow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal flow equilibrium state, sustained primarily by energetic alpha particles from fusion reactions, will require the consideration of all of these key elements to predict quantitatively the stability and evolution. The principal investigators on this proposal are leading experts in the relevant theoretical and computational areas, and aimed to perform computations guided by analytic modeling, to address this physics in realistic configurations. The overall goal is to understand the key physics mechanisms that describe resistive toroidal burning plasmas, surrounded by a resistive wall, under active feedback control. With the physics of the energetic ions, resonant layers, resistive wall, and toroidal momentum transport included, this study will extend from recent publications in theory and simulation of individual effects and move toward predictive modeling for burning plasmas.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Potential Vorticity Mixing in a Tangled Magnetic Field

A theory of potential vorticity (PV) mixing in a disordered (tangled) magnetic field is introduced. The analysis is in the context of β -plane MHD, with a special focus on the physics of momentum transport in the stably stratified, quasi-2D solar tachocline. A physical picture of mean PV evolution by vorticity advection and tilting of magnetic fields is proposed. In the case of weak field perturbations, quasi-linear theory predicts that the Reynolds and magnetic stresses balance as turbulence Alfvénizes for a larger mean magnetic field. Jet formation is explored quantitatively in the mean field–resistivity parameter space. However, since even a modest mean magnetic field leads to large magnetic perturbations for large magnetic Reynolds number, the physically relevant case is that of a strong but disordered field. We show that numerical calculations indicate that the Reynolds stress is modified well before Alfvénization—i.e., before fluid and magnetic energies balance. To understand these trends, a double-average model of PV mixing in a stochastic magnetic field is developed. Calculations indicate that mean-square fields strongly modify Reynolds stress phase coherence and also induce a magnetic drag on zonal flows. The physics of transport reduction by tangled fields is elucidated and linked to the related quench of turbulent resistivity. We propose here a physical picture of the system as a resisto-elastic medium threaded by a tangled magnetic network. Applications of the theory to momentum transport in the tachocline and other systems are discussed in detail.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

A numerical simulation of magnetic reconnection and radiative cooling in line-tied current sheets

Radiative MHD equations are used for an optically thin plasma to carry out a numerical experiment related to the formation of 'postflare' loops. The numerical experiment starts with a current sheet that is in mechanical and thermal equilibrium but is unstable to both tearing-mode and thermal-condensation instabilities. The current sheet is line-tied at one end to a photospheric-like boundary and evolves asymmetrically. The effects of thermal conduction, resistivity variation, and gravity are ignored. In general, reconnection in the nonlinear stage of the tearing-mode instability can strongly affect the onset of condensations unless the radiative-cooling time scale is much smaller than the tearing-mode time scale. When the ambient plasma is less than 0.2, the reconnection enters a regime where the outflow from the reconnection region is supermagnetosonic with respect to the fast-mode wave speed. In the supermagnetosonic regime the most rapidly condensing regions occur downstream of a fast-mode shock that forms where the outflow impinges on closed loops attached to the photospheric-like boundary. A similar shock-induced condensation might occur during the formation of 'postflare' loops.

Forbes, T. G.↗
DOE OSTI2021

Suppression of MHD modes with active phase-control of probe-injected currents

Abstract Active phase-control of probe-injected current is shown to both suppress and amplify long-wavelength rotating magnetohydrodynamic instabilities in the HBT-EP tokamak. Four probes are connected in quadrature and energized to drive non-axisymmetric currents through the edge of the tokamak, creating magnetic perturbations comparable to previously-studied saturated kink modes or resonant magnetic perturbations that are generated by an external control coil array. Measurements of the magnetic perturbations from the probe-injected currents determine a set of current-carrying helical filaments used to model active feedback control of resistive wall modes. These experiments suggest current-injection feedback may be an effective alternative to external control coils for control of RWMs and other long-wavelength kink-like modes at the edge of tokamaks.

Physics↗
NASA NTRS1984

Externally driven magnetic reconnection

A model is presented in which externaly driven reconnection is simulated by solving the MHD equations in an initially plane current sheet. Magnetic reconnection in the earth's magnetotail is widely believed to be the direct cause of magnetospheric substorms. Both 2-D and 3-D versions of the model have been developed. It is postulated that connection in the tail is triggered by a local compression of the plasma sheet which results from an invasion of the solar wind into the magnetotail. Thus, the simulation is started by introducing flow from the lobes normal to the plasma sheet. When resistivity is generated in a local region of the neutral sheet, reconnection develops and magnetic energy is converted into plasma bulk flow. Although the driven reconnection model is highly simplified, it can aid in understanding many features of substorms in the tail; in particular, results show that rapid flows both earthward and tailward of the neutral line and the nightside substorm current system are natural consequences of driven magnetic reconnection.

Walker, R. J.↗
NASA NTRS1987

Equilibrium and tearing stability of thin current layers in magnetic arcades

The MHD equilibrium of a thin, force-free current layer in a magnetic arcade is solved analytically. Various approximations are made in order to achieve a lowest order description that is physically relevant as well as mathematically tractable: The arcade is assumed to emanate from feet that are well localized, the current emanates from a localized sheet within the feet and is relatively weak. The resulting expressions for B (sup arrow) are relatively simple and natural flux coordinates are identifiable. The stability of such a current sheet to resistive filamentation is then investigated.

Hassam, A. B.↗
NASA NTRS1990

The kink instability in infinite cylindrical flux tubes - Eigenvalues for power-law twist profiles

Simple, accurate methods of calculating ideal MHD instability eigenvalues for infinitely long cylindrical tubes with twist function T(r) are developed. The results show that the most rapidly growing and energetic instabilities occur in the Gold-Hoyle v = 0 field, with the instability progressively weakening with increasing v. However, the maximum force eigenvalue is always small, so that even in the Gold-Hoyle case only a small proportion of the available magnetic energy can be released in the linear phase. The results also confirm that the linear pinch is remarkably weak yet relatively resistant to line-tying. It is shown that the weakness of the force eigenvalue implies that the influence of uniform gas pressure on stability is negligible. Implications for the energy-release mechanism in solar flares are discussed.

Craig, I. J. D.↗
NASA NTRS1995

Electric conductivity of plasma in solar wind

One of the most important parameters in MHD description of the solar wind is the electric conductivity of plasma. There exist now two quite different approaches to the evaluation of this parameter. In the first one a value of conductivity taken from the most elaborated current theory of plasma should be used in calculations. The second one deals with the empirical, phenomenological value of conductivity. E.g.: configuration of interplanetary magnetic field, stretched by the expanding corona, depends on the magnitude of electrical conductivity of plasma in the solar wind. Knowing the main empirical features of the field configuration, one may estimate the apparent phenomenological value of resistance. The estimations show that the electrical conductivity should be approximately 10(exp 13) times smaller than that calculated by Spitzer. It must be noted that the empirical value should be treated with caution. Due to the method of its obtaining it may be used only for 'large-scale' description of slow processes like coronal expansion. It cannot be valid for 'quick' processes, changing the state of plasma, like collisions with obstacles, e.g., planets and vehicles. The second approach is well known in large-scale planetary hydrodynamics, stemming from the ideas of phenomenological thermodynamics. It could formulate real problems which should be solved by modern plasma physics, oriented to be adequate for complicated processes in space.

Chertkov, A. D.↗
NASA NTRS2008

The Adaptively-Refined MHD Solver (ARMS)

ARMS is a robust code for solving the single-fluid, fully compressible equations of magnetohydrodynamics. In this presentation I will describe how to access and run ARMS and how to use HelioSpace to analyze the results. The main features of the code are the following: 1. Equations are solved in either Cartesian or Spherical coordinates. For the latter, either a linear or a Log(r) coordinate can be used. 2. The latest flux-corrected-transport schemes developed at NRZ, are used to advance the convective equations. All waves and shocks are accurately captured. 3. A fully adaptive grid can be used. The code incorporates the PARAMESH technology for block-adaptation. Arbitrary adaptation criteria can be specified. 4. The code maintains a divergence-free magnetic field to machine accuracy throughout all refinement and de-refinement. 5. A uniform resistivity and thermal conduction can be included. Solution is either explicit or implicit. 6. Arbitrary boundary and initial conditions can be used in either Cartesian or spherical coordinates. 7. An accompanying visualization package, HelioSpace, is available for viewing and analyzing the native ARMS output, (without need for interpolating onto a structured grid).

Antichos, S. K.↗
DOE OSTI2020

Thermomagnetic instability of plasma composition gradients

In this paper, we show that, under Braginskii magneto-hydrodynamics, anti-parallel gradients in an average ion charge state and electron temperature can be unstable to the growth of self-generated magnetic fields. The instability is analogous to the field-generating thermomagnetic instability, although it is driven by the collisional thermal force magnetic source term rather than the Biermann battery term. The gradient in ion charge state causes a gradient in collisionality, which couples with temperature perturbations to create a self-generated magnetic field. This magnetic field deflects the electron heat flux in a way that reinforces the temperature perturbation. The derived linearized growth rate, typically on hydrodynamic timescales, includes resistive and thermal smoothing. It increases with large ion composition gradients and electron heat flux, conditions typical of the hohlraum walls or contaminant mix jets in inertial confinement fusion implosions. However, extended magneto-hydrodynamic simulations indicate that the instability is usually dominated and stabilized by nonlinear Nernst advection, in a similar manner to the standard thermomagnetic instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

The substorm current wedge and field-aligned currents in MHD simulations of magnetotail reconnection

Results are reported from a three-dimensional MHD simulation of magnetotail reconnection, which, in contrast to earlier ones, starts from an equilibrium configuration that contains a region 1 type field-aligned current system at the plasma sheet/lobe boundary. This current system is found from equilibrium theory as the consequence of boundary conditions, which require that field lines at large distances from the earth become aligned-with the tail axis, consistent with observations and the expected consequences of the interaction with the solar wind. As in earlier simulations, the dynamic evolution develops from a slow diffusion to the fast growth of a three-dimensional tearing mode, initiated by the sudden occurrence or increase of resistivity, leading again to the formation and subsequent tailward ejection of a plasmoid, associated with fast plasma flows. In addition there are changes of the electric current system which are consistent with the substorm current wedge picture inferred from observations.

Birn, J.↗