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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Forced magnetic reconnection in a plasma sheet with localized resistivity profile excited by lower hybrid drift type instability

A forced magnetic reconnection process with a temporal evolution of resistivity is studied for a plasma sheet with a nonuniform resistivity profile based on the nonlocal mode structure of the lower hybrid drift type instability. The growth rate of the mode found is almost independent of the resistivity at the neutral sheet, but depends on the resistivity of the region of maximum density gradient away from the neutral sheet. This is studied by using both a nonlinear numerical MHD simulation and a linear theory. The mode may be relevant to the prevalent theoretical concept of MHD reconnection and the localized anomalous resistivity profile based on the lower hybrid drift instability.

Hoshino, M.↗

DOE OSTI2022

Critical role of current-driven instabilities for ELMs in NSTX

The impact of different extended-magnetohydrodynamic (MHD) contributions on the stability thresholds of peeling-ballooning modes in ELMing and ELM-free plasmas in the spherical tokamak NSTX is investigated with the initial value code M3D-C1. We show that ELMing discharges in NSTX are limited by resistive current-driven peeling modes, whereas non-ELMing wide-pedestal H-mode discharges are located near the ideal pressure-driven ballooning threshold. It is demonstrated that extended-MHD can lead to more reliable edge stability predictions than existing ideal-MHD models. Resistive peeling-ballooning modes are found to exist well before the ideal stability threshold is met, and kink-peeling modes exhibit considerable sensitivity to plasma resistivity. Other effects not considered in ideal-MHD models affect PB modes in NSTX in a weaker way. Gyroviscous stress appears stabilizing such that the stability boundary lies closer to the experimental point. Equilibrium rotation can suppress ideal core modes and thus isolate edge modes. These results are important for the development of a predictive pedestal model for low-aspect ratio tokamaks.

pedestal↗

NASA NTRS2004

Adaptive Numerical Dissipative Control in High Order Schemes for Multi-D Non-Ideal MHD

The goal is to extend our adaptive numerical dissipation control in high order filter schemes and our new divergence-free methods for ideal MHD to non-ideal MHD that include viscosity and resistivity. The key idea consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free of numerical dissipation contamination. These scheme-independent detectors are capable of distinguishing shocks/shears, flame sheets, turbulent fluctuations and spurious high-frequency oscillations. The detection algorithm is based on an artificial compression method (ACM) (for shocks/shears), and redundant multi-resolution wavelets (WAV) (for the above types of flow feature). These filter approaches also provide a natural and efficient way for the minimization of Div(B) numerical error. The filter scheme consists of spatially sixth order or higher non-dissipative spatial difference operators as the base scheme for the inviscid flux derivatives. If necessary, a small amount of high order linear dissipation is used to remove spurious high frequency oscillations. For example, an eighth-order centered linear dissipation (AD8) might be included in conjunction with a spatially sixth-order base scheme. The inviscid difference operator is applied twice for the viscous flux derivatives. After the completion of a full time step of the base scheme step, the solution is adaptively filtered by the product of a 'flow detector' and the 'nonlinear dissipative portion' of a high-resolution shock-capturing scheme. In addition, the scheme independent wavelet flow detector can be used in conjunction with spatially compact, spectral or spectral element type of base schemes. The ACM and wavelet filter schemes using the dissipative portion of a second-order shock-capturing scheme with sixth-order spatial central base scheme for both the inviscid and viscous MHD flux derivatives and a fourth-order Runge-Kutta method are denoted.

Yee, H. C.↗

DOE OSTI2023

Tech-X contributions to the Center for Tokamak Transients (CTTS) Scientific Discovery through Advanced Computing (SciDAC) project (Final Report)

The Center for Tokamak Transients Simulations (CTTS) consists of computational plasma physicists, computer scientists, and applied mathematicians who are together developing and applying simulation software to improve our understanding of and ability to predict and control transient events in tokamak discharges, particularly those that lead to violent disruptions. We have major simulation activities concerning disruptions that are caused by violating ideal MHD stability criteria, vertical displacement events, resistive wall modes, and tearing modes. We are performing comprehensive modeling of disruption mitigation techniques, particularly by shattered pellet injection (SPI). The Center is built around two mature extended MHD codes: NIMROD and M3D-C1. These codes are being extended in a number of ways as required by the demands of disruption modeling, and are being optimized to work well for larger problem size and modern computer architectures as required to model tokamak plasmas with realistic parameters over the required timescales. The codes will be coupled to the FronTier code for state-of-the-art modeling of SPI. Validation with experimental data on DIII-D, NSTX, and JET is being facilitated by rigorous uncertainty quantification (UQ). The high-level objectives of the center are to: (1) Develop better understanding and improved predictive capability of when crossing a linear stability boundary will likely lead to a disruption, (2) Provide a tool for calculating vessel forces, both axisymmetric and non-axisymmetric, for worst case VDE and other disruptions in tokamaks, (3) Develop increased understanding and improved predictive capability for locked mode disruptions and how best to avoid them, (4) Develop, verify, and validate 3D models for disruption mitigation by shattered pellets that can be used to design and optimize a system for future tokamaks, and (5) Increase the efficiency and scalability of the M3D-C1 and NIMROD codes so that they can model disruptions in ITER scale devices using the most powerful DOE HPC hardware, and (6) Demonstrate code-coupling via FronTier coupling with NIMROD and M3D-C1 and by interfacing these codes within a WDM framework.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1990

Investigation of turbulent processes in magnetospheric boundary layers

A self-consistent non-evolving two dimensional slab model of a viscous low-latitude boundary layer (LLBL) coupled to the ionosphere was developed by Phan, et al., (1989). Numerical results from the model and possible use of observations to determine the model parameters are discussed. The dynamical model developed by Lotko, et al., (1987) was used by Lotko and Shen (1991) to examine dynamical processes relevant to the LLBL with particular application to post-noon auroral shear layers. Initial results from a magnetohydrodynamic study of flank-side mangetopause boundary configuration are described. Effects of compressibility, scalar viscosity, and electrical resistivity are included in the MHD equations.

Lotko, William↗

NASA NTRS2005

Adaptive Numerical Dissipation Control in High Order Schemes for Multi-D Non-Ideal MHD

The required type and amount of numerical dissipation/filter to accurately resolve all relevant multiscales of complex MHD unsteady high-speed shock/shear/turbulence/combustion problems are not only physical problem dependent, but also vary from one flow region to another. In addition, proper and efficient control of the divergence of the magnetic field (Div(B)) numerical error for high order shock-capturing methods poses extra requirements for the considered type of CPU intensive computations. The goal is to extend our adaptive numerical dissipation control in high order filter schemes and our new divergence-free methods for ideal MHD to non-ideal MHD that include viscosity and resistivity. The key idea consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free from numerical dissipation contamination. These scheme-independent detectors are capable of distinguishing shocks/shears, flame sheets, turbulent fluctuations and spurious high-frequency oscillations. The detection algorithm is based on an artificial compression method (ACM) (for shocks/shears), and redundant multiresolution wavelets (WAV) (for the above types of flow feature). These filters also provide a natural and efficient way for the minimization of Div(B) numerical error.

Yee, H. C.↗

NASA NTRS2012

Comparison of Computed and Measured Performance of a Pulsed Inductive Thruster Operating on Argon Propellant

Pulsed inductive plasma accelerators are electrodeless space propulsion devices where a capacitor is charged to an initial voltage and then discharged through a coil as a high-current pulse that inductively couples energy into the propellant. The field produced by this pulse ionizes the propellant, producing a plasma near the face of the coil. Once a plasma is formed if can be accelerated and expelled at a high exhaust velocity by the Lorentz force arising from the interaction of an induced plasma current and the magnetic field. A recent review of the developmental history of planar-geometry pulsed inductive thrusters, where the coil take the shape of a flat spiral, can be found in Ref. [1]. Two concepts that have employed this geometry are the Pulsed Inductive Thruster (PIT)[2, 3] and the Faraday Accelerator with Radio-frequency Assisted Discharge (FARAD)[4]. There exists a 1-D pulsed inductive acceleration model that employs a set of circuit equations coupled to a one-dimensional momentum equation. The model was originally developed and used by Lovberg and Dailey[2, 3] and has since been nondimensionalized and used by Polzin et al.[5, 6] to define a set of scaling parameters and gain general insight into their effect on thruster performance. The circuit presented in Fig. 1 provides a description of the electrical coupling between the current flowing in the thruster I1 and the plasma current I2. Recently, the model was upgraded to include an equation governing the deposition of energy into various modes present in a pulsed inductive thruster system (acceleration, magnetic flux generation, resistive heating, etc.)[7]. An MHD description of the plasma energy density evolution was tailored to the thruster geometry by assuming only one-dimensional motion and averaging the plasma properties over the spatial dimensions of the current sheet to obtain an equation for the time-evolution of the total energy. The equation set governing the dynamics of the coupled electrodynamic-current sheet system is composed of first-order, coupled ordinary differential equations that can be easily solved numerically without having to resort to much more complex 2-D finite element plasma simulations.

Polzin, Kurt A.↗

NASA NTRS1985

Resistive decay of Alfven waves in a non-uniform plasma

The effect of resistive dissipation on the propagation of an MHD disturbance in a nonuniform plasma is examined. The present analysis, based on a boundary-layer technique, shows the existence of resistive normal modes with complex eigenfrequencies. The real part of the eigenfrequency is associated with an oscillatory behavior and defines the location in space of the layer where resistivity is important. The dissipation mechanism is responsible for the damping of the wave, in contrast with previous works in which the ideal MHD theory was used.

Mok, Y.↗

NASA NTRS1986

Magnetic energy storage and conversion in the solar atmosphere

According to the approach employed in this investigation, particularly important simple configurations of magnetic field and plasma are identified, and it is attempted to achieve an understanding of the large-scale dynamic processes and transformations which these systems can undergo. Fundamental concepts are discussed, taking into account aspects of magnetic energy generation, ideal MHD theory, non-MHD properties, the concept of 'anomalous' resistivity, and global electrodynamic coupling. Questions of magnetically controlled energy conversion are examined, giving attention to magnetic modifications of plasma transport, the transition region structure and flows, channeling and acceleration of plasma, channeling and dissipation of MHD waves, and anomalous dissipation of field-aligned currents. A description of the characteristics of magnetohydrodynamic energy conversion is also provided, and outstanding questions are discussed.

Spicer, D. S.↗

NASA NTRS1989

The stability analysis of magnetohydrodynamic equilibria - Comparing the thermodynamic approach with the energy principle

This paper is a contribution to the stability analysis of current-carrying plasmas, i.e., plasma systems that are forced by external mchanisms to carry a nonrelaxing electrical current. Under restriction to translationally invariant configurations, the thermodynamic stability criterion for a multicomponent plasma is rederived within the framework of nonideal MHD. The chosen dynamics neglects scalar resistivity, but allows for other types of dissipation effects both in Ohm's law and in the equation of motion. In the second section of the paper, the thermodynamic stability criterion is compared with the ideal MHD based energy principle of Bernstein et al. With the help of Schwarz's inequality, it is shown that the former criterion is always more 'pessimistic' than the latter, i.e., that thermodynamic stability implies stability according to the MHD principle, but not vice versa. This reuslt confirms the physical plausible idea that dissipational effects tend to weaken the stability properties of current-carrying plasma equilibria by breaking the constraints of ideal MHD and allowing for possibly destabilizing effects such as magnetic field line reconfiguration.

Brinkmann, R. P.↗

NASA NTRS2012

Multi-Scale Modeling of Magnetospheric Dynamics

Magnetic reconnection is a key element in many phenomena in space plasma, e.g. Coronal mass Ejections, Magnetosphere substorms. One of the major challenges in modeling the dynamics of large-scale systems involving magnetic reconnection is to quantifY the interaction between global evolution of the magnetosphere and microphysical kinetic processes in diffusion regions near reconnection sites. Recent advances in small-scale kinetic modeling of magnetic reconnection significantly improved our understanding of physical mechanisms controlling the dissipation in the vicinity of the reconnection site in collisionless plasma. However the progress in studies of small-scale geometries was not very helpful for large scale simulations. Global magnetosphere simulations usually include non-ideal processes in terms of numerical dissipation and/or ad hoc anomalous resistivity. Comparative studies of magnetic reconnection in small scale geometries demonstrated that MHD simulations that included non-ideal processes in terms of a resistive term 11 J did not produce fast reconnection rates observed in kinetic simulations. In collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is nongyrotropic pressure effects with spatial scales comparable with the particle Larmor radius. We utilize the global MHD code BATSRUS and replace ad hoc parameters such as "critical current density" and "anomalous resistivity" with a physically motivated model of dissipation. The primary mechanism controlling the dissipation in the vicinity of the reconnection site in incorporated into MHD description in terms of non-gyrotropic corrections to the induction equation. We will demonstrate that kinetic nongyrotropic effects can significantly alter the global magnetosphere evolution. Our approach allowed for the first time to model loading/unloading cycle in response to steady southward IMF driving. The role of solar wind parameters and ionospheric conductance on the dynamics of the loading/unloading cycle will be discussed.

Kuznetsova, M. M.↗

NASA NTRS1995

Heating of the solar middle chromosphere by large-scale electric currents

A global resistive, two-dimensional, time-dependent magnetohydrodynamic (MHD) model is used to introduce and support the hypothesis that the quiet solar middle chromosphere is heated by resistive dissipation of large-scale electric currents which fill most of its volume. The scale height and maximum magnitude of the current density are 400 km and 31.3 m/sq m, respectively. The associated magnetic field is almost horizontal, has the same scale height as the current density, and has a maximum magnitude of 153 G. The current is carried by electrons flowing across magnetic field lines at 1 m/s. The resistivity is the electron contribution to the Pedersen resitivity for a weakly ionized, strongly magnetized, hydrogen gas. The model does not include a driving mechanism. Most of the physical quantities in the model decrease exponentially with time on a resistive timescale of 41.3 minutes. However, the initial values and spatial; dependence of these quantities are expected to be essentially the same as they would be if the correct driving mechanism were included in a more general model. The heating rate per unit mass is found to be 4.5 x 10(exp 9) ergs/g/s, independent of height and latitude. The electron density scale height is found to be 800 km. The model predicts that 90% of the thermal energy required to heat the middle chromosphere is deposited in the height range 300-760 km above the temperature minimum. It is shown to be consistent to assume that the radiation rate per unit volume is proportional to the magnetic energy density, and then it follows that the heating rate per unit volume is also proportional to the energy from the photosphere into the overlying chromosphere are briefly discussed as possible driving mechanisms for establishing and maintaining the current system. The case in which part of or all of the current is carried by protons and metal ions, and the contribution of electron-proton scattering to the current are also considered, with the conclusion that these effects do not change the qualitative prediction of the model, but probably change the quantitative predictions slightly, mainly by increasing the maximum magntiude of the current density and magnetic field to at most approximately 100 mA/m and approximately 484 G, respectively. The heating rate per unit mass, current density scale height, magnetic field scale height, temperatures, and pressures are unchanged or are only slightly changed by including these additional effects due to protons and ions.

Goodman, M. L.↗

NASA NTRS2007

Multi-Scale Modeling of Magnetospheric Reconnection

One of the major challenges in modeling the magnetospheric magnetic reconnection is to quantify the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions near reconnection sites. There is still considerable debate as to what degree microphysical processes on kinetic scales affect the global evolution and how important it is to substitute numerical dissipation and/or ad hoc anomalous resistivity by a physically motivated model of dissipation. Comparative studies of magnetic reconnection in small scale geometries demonstrated that MHD simulations that included non-ideal processes in terms of a resistive term $\eta J$ did not produce the fast reconnection rates observed in kinetic simulations. For a broad range of physical parameters in collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. We utilize the global MHD code BATSRUS and incorporate nongyrotropic effects in diffusion regions in terms of corrections to the induction equation. We developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated selfconsistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites is updated during the simulations. To clarify the role of nongyrotropic effects in diffusion region on the global magnetospheric dynamic we perform simulations with steady southward IMF driving of the magnetosphere. Ideal MHD simulations with magnetic reconnection supported by numerical resistivity produce steady configuration with almost stationary near-earth neutral line (NENL). Simulations with non-gyrotropic corrections demonstrate dynamic quasi-periodic response to the steady driving condition. The loading/unloading cycle in non-gyrotropic MHD results has a non-stationary reconnection site in the magnetotail, with the retreating during the stretching phase and then a new NENL forming in the resulting thin plasma sheet. We expect that this model will lead to improved representations of space weather event in the magnetosphere.

Kuznetsova, M. M.↗

DOE OSTI2025

Simulations of saturated MHD activity in the HBT-EP tokamak

The NIMROD code is used to perform simulations of MHD activity in the HBT-EP tokamak, including the effect of a resistive wall. Linear simulations are conducted to scan current and pressure profiles to establish self-consistent baseline equilibria that result in MHD instabilities with low error between experimental diagnostics and synthetic equivalents. A family of such equilibria is found with relatively peaked current and pressure profiles. Nonlinear simulations are performed using the optimal equilibrium, which is linearly unstable to a Resistive Wall Mode—stable in the ideal wall limit and with an Alfvénic growth rate in the no-wall limit. Using resistive wall boundary conditions, nonlinear Hall MHD simulations yield a rotating, saturated n=1 mode similar to the type observed experimentally in HBT-EP discharges. An n=1 mode around the q=2 surface mediates a localized axisymmetric perturbation that flattens the plasma current from a linearly unstable state to a 3D, rotating, stable state. During the saturated period, periodic changes in inter-diagnostic phase lag allow for qualitative estimations of effective, macroscopic transport coefficients used in MHD closure models and could be extended to model post-disruption MHD structure rotation after thermal quenches seen in HBT-EP.

Arnold, David A. (ORCID:0009000059535761)↗

DOE OSTI2020

Low- n global ideal MHD instabilities in the CFETR baseline scenario

This article reports an evaluation on the linear ideal magnetohydrodynamic (MHD) stability of the China Fusion Engineering Test Reactor (CFETR) baseline scenario for various first-wall locations. The initial-value code NIMROD and eigen-value code AEGIS are employed in this analysis. Despite the distinctly different approaches in modeling the scrape off layer region, the dominant growth in each of the low-n (n=1–10) modes are consistent between the two codes. We report the higher-n modes are dominated by ballooning modes and localized in the pedestal region, while the lower-n modes have more prominent external kink components and broader mode profiles. The influence of the plasma-vacuum profile and wall shape are also examined using NIMROD. In the presence of a resistive wall, the low-n ideal MHD instabilities are further studied using AEGIS. For the designed first-wall location, the n=1 resistive wall mode is found to be unstable, which could be fully stabilized by uniform toroidal rotation above the 2.9% core Alfven speed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1986

The effect of different resistivity models on magnetotail dynamics

Two-dimensional, time-dependent MHD simulations of the interaction between the solar wind and the earth's magnetosphere have been performed to study magnetotail dynamics with varying forms of anomalous resistivity. In general, the resulting models conform to the neutral line model proposed for substorms with near-earth x points and high-speed tailward flows occurring in the magnetotail. However, in the case where Joule heating is included in the physical description of the system, the near-earth x point never moves far down tail, and high-speed tailward flows never stop. Only in the case where there is no Joule heating does the x point move down tail. Simultaneously, the high-speed tailward flows cease. These results indicate that the mechanism of energy dissipation can have an important effect on reconnection processes and the global magnetospheric dynamics.

Lyon, J. G.↗

NASA NTRS1985

Resistive Alfven normal modes in a non-uniform plasma

Resistive normal mode solutions of the MHD equations are found numerically in a smooth, nonuniform, magnetic field. The (alpha, S) boundary within which normal mode solutions exist is explicitly computed, where alpha is the normalized wavenumber and S the Lundquist number. As an extension of previous analytic results of Mok and Einaudi (1985), the damping rate of these modes is computed to a higher accuracy, and is found to have an a + bS exp -1/3 dependence, where a and b are independent of S.

Einaudi, G.↗

DOE OSTI2020

The role of edge resonant magnetic perturbations in edge-localized-mode suppression and density pump-out in low-collisionality DIII-D plasmas

Two-fluid nonlinear MHD simulations using the TM1 code demonstrate that the formation of magnetic islands at the top and bottom of the H-mode pedestal, together with the strong screening of resonant fields in the gradient region of the pedestal, can account for ELM suppression and density pump-out by n = 2 resonant magnetic perturbations (RMPs) in low-collisionality DIII-D ITER similar shape (ISS) plasmas. Using experimentally relevant transport coefficients, neoclassical resistivity, electron collisionality, and RMP amplitudes, nonlinear MHD simulations reproduce the observed level of density reduction (density pump-out) in DIII-D due to the formation of narrow magnetic islands in the resistive foot of pedestal. For large amplitude RMPs ( B r / B t > 1 × 10 -4 ) simulations predict resonant field penetration and resulting in pressure reduction at the top of the pedestal, consistent with experimental observations at the onset of ELM suppression. The predicted reduction in the height and width of the pedestal by magnetic island enhanced collisional transport provides a quantitative mechanism for the stabilization of the peeling-ballooning mode (PBM). Importantly, these simulations predict strong screening of resonant fields in the steep gradient region of the pedestal due to strong E × B rotation and diamagnetic flows. Yet, if the plasma resistivity is made artificially larger (~ 10X) than neoclassical, the simulations predict magnetic stochasticity throughout the plasma edge and the collapse of the pedestal due to the reduction in the penetration threshold with increasing resistivity. A scaling relation for resonant field penetration at the pedestal top, using several hundred nonlinear simulations, reproduces the density and E × B rotation dependence of the ELM suppression threshold observed in DIII-D. Simulations using ITER model equilibria indicate that the penetration threshold at the top of the ITER pedestal will be lower than in DIII-D due to the anticipated lower perpendicular flow velocities in ITER, resulting in the weaker screening of resonant fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗