Search NASA⌕ Search

SEARCH · Search NASA

Results for “resistive instabilities”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Exploring the transition from continuous turbulence fluctuations to bursting ELMs in high SOL density regimes

BOUT++ turbulence simulations of the DIII-D reveal that the density profile between the separatrix and pedestal plays a crucial role in the dynamics of edge localized modes (ELMs) and edge plasma turbulent transport. Nonlinear simulations demonstrate that small ELMs in the DIII-D hybrid scenario under high SOL density conditions are predominantly driven by local ballooning modes near the separatrix, stabilizing global instabilities while enhancing localized pressure fluctuations. A key control parameters for ELM dynamics is the separatrix-to-pedestal density ratio, n e,sep /n e,ped . A high ratio indicates a shallow gradient, favoring small ELMs, while a lower ratio signals a steep gradient, which increases the likelihood of large ELMs. Comprehensive parameter scans, including n e,sep /n e,ped , density gradient profiles near the separatrix, and resistivity, reveal the critical role of these parameters in shaping transitions between turbulence-driven transport and ELM bursting. The scans demonstrate that in high SOL density regimes, small ELMs can result from either global resistive MHD instabilities or local ballooning modes near the separatrix, depending on the steepness of the separatrix density gradient. These findings also highlight the transition from continuous turbulence to small ELMs. The post-crash peak in pressure fluctuations, δP rms serves as a critical metric for identifying transition from continuous turbulence fluctuations to ELM bursting. Larger δP rms values correlate with ELM bursts driven by local or global instabilities, whereas smaller values indicate turbulence-dominated transport. Drift-Alfvén and resistive ballooning turbulence enhance the entrainment of fluctuations from the pedestal to the SOL, contributing to the complex interplay of dynamics in this regime. These findings emphasize the importance of separatrix density shaping and pedestal gradient control for optimizing ELM behavior in ITER and future fusion devices.

Li, Nami [Lawrence Livermore National Laboratory (↗
NASA NTRS1992

MHD modeling of magnetotail instability for anisotropic pressure

Results of three-dimensional MHD simulations of magnetotail dynamics with anisotropic pressure are presented. The pressure tensor is assumed to be gyrotropic, satisfying a modified double adiabatic approximation including Ohmic heating. It is found that these constraints tend to stabilize the tail with respect to the resistive tearing instability. Including an increasing level of additional isotropization results in a gradual transition to the fast instability found in an isotropic model, which is discussed for comparison. Possible isotropization mechanisms yielding magnetotail instability are discussed also. The stabilizing anisotropies are found primarily in the plasma sheet boundary region. This fact indicates a possibly important role of isotropization mechanisms operating in this region in destabilizing the magnetotail and initiating a tearing instability without the necessity of changes in fluctuation levels in the neutral sheet.

Hesse, Michael↗
DOE OSTI2021

Computation of linear MHD instabilities with the multi-region relaxed MHD energy principle

Over the last decade, a variational principle based on a generalisation of Taylor's relaxation, referred to as multi-region relaxed magnetohydrodynamics (MRxMHDs) has been developed. The numerical solutions of the MRxMHD equilibria have been constructed using the Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502). In principle, SPEC could also be established to describe the MRxMHD stability of an equilibrium. In this work, a theoretical framework is developed to relate the second variation of the energy functional to the so-called Hessian matrix, enabling the prediction of MHD linear instabilities of cylindrical plasmas, and is implemented in SPEC. The negative and positive eigenvalues of the Hessian matrix predict the stability of an equilibrium. Verification studies of SPEC stability results with the M3D-C 1 code and the tearing mode $\Delta^{^{\prime}}$ criterion have been conducted for ideal and resistive MHD instabilities, respectively, in a pressureless cylindrical tokamak, and have shown good agreement. Overall, our stability analysis is a critical step towards understanding the MHD stability of three-dimensional MHDs where nested flux surfaces, magnetic islands and stochastic regions co-exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Resistive Wall Simulations with LAMDA

The prospect for a new downstream transport (DST) beam line for multipulse target testing (MPTL) has raised questions about the resistive-wall instability, which has been theoretically shown to couple between the multiple pulses in a pulse train. Recent upgrades to computer operating systems and the LAMDA beam dynamics code mandate that the code be revalidated before using to explore MPTL parameters. This note is a brief account of the revalidation of LAMDA for resistive wall problems.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Linear and nonlinear benchmarks between the CLT code and the M3D-C1 code for the 2/1 resistive tearing mode and the 1/1 resistive kink mode

Here, the linear and nonlinear benchmarks between the CLT code and the M3D-C1 code for the 2/1 resistive tearing mode and the 1/1 resistive kink mode are presented. CLT is an explicit finite difference code, while M3D-C1 is an implicit finite element code. Although the implementations of CLT and M3D-C1 are totally different, we find that the simulation results of the resistive-kink mode and the m/n=2/1 tearing mode from M3D-C1 and CLT are almost the same, including the linear and nonlinear growth rates, the mode structures, the nonlinear saturation levels, the Poincare plots, and the scaling laws. This confirms that the nonlinear results for the 1/1 resistive-kink mode and 2/1 tearing mode are accurate and reliable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The Impact of Anomalous Resistivity in Vacuum Contaminant Plasmas on the Electrothermal Instability

This manuscript presents an assessment of the electrothermal instability (ETI) in the presence of anomalous resistivity (AR) in vacuum contaminant plasmas (VCP) when applied to a magnetized liner inertial fusion (MagLIF)-like load.Pulsed-power driven dielectrically coated metallic liners, like in MagLIF, experience the current-driven electrothermal instability which occurs when a material’s resistivity changes with temperature and is subject to ohmic heating. Large scale pulsed-power facilities that use magnetically insulated transmission lines (MITL) have been shown to generate low-density plasma which enters the target chamber and coalesces around the load. The low-density high-temperature vacuum contaminant plasmas (VCP) can parasitically divert current from the load through causing a short in the anode-cathode gap inside the target chamber. Resistive magnetohydrodynamic (MHD) simulations of these VCP experience unphysical runaway ohmic heating due to under predicting the resistivity by using a purely collisional resistivity model.AR provides a physics-based way to address this runaway heating through increasing the resistivity in a proportional way with the drift speed. In this work, 1D simulations probe the effect that AR in VCP has on the magnetic diffusion rate, and 2D simulations show how this effect manifests in the nonlinear striation form of the ETI for a MagLIF-like load. Beryllium and aluminum dielectrically coated liners are used for the 1D and 2D simulations in this work. The1D simulations show that a VCP causes a delay in the current delivery to the load by upwards of 8 ns at 60 ns into atypical current pulse for an 1 × 10 7 A scale pulsed-power accelerator. The 2D simulations show the delay observed in the 1D simulations is reduced substantially (4 ns) by 2D turbulence that disrupts the VCP layer early in time. Here, the 2Dsimulated ETI growth varies across AR models, more so for the beryllium liner than the aluminum liner because the beryllium liner shows an enhanced rate of penetration for the magnetic diffusion wave in comparison to the aluminum liner. The 2D simulations show the bulk dielectric thickness varies across AR models with the Davidson AR model beingthe largest and the Buneman AR model being the smallest, and in connection with the thickness the Rayleigh-Taylor bubble-spike distances varies correspondingly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Toroidal modeling of runaway electron loss due to 3D fields in ITER

Mitigation of runaway electrons (REs) by three-dimensional (3D) magnetic field perturbations is numerically investigated for the ITER 15 MA baseline D–T scenario, utilizing the MARS-F code (Liu et al Phys. Plasmas 7 3681) with a drift orbit test particle tracing module. Considered are two types of 3D fields: the n = 3 (n is the toroidal mode number) resonant magnetic perturbation (RMP) utilized for the purpose of controlling the edge localized modes in ITER, and perturbations generated by the n = 1 magneto-hydrodynamic (MHD) instabilities in a post-disruption plasma. The RMP field, applied to a pre-disruption plasma, is found to be moderately effective in mitigating the RE seeds in ITER when vacuum field model is assumed. Up to ~40% loss fraction is possible at 90 kA-turn coil current. The mitigation efficiency is however substantially reduced, down to less than 5%, when the plasma response is taken into account. This is due to strong screening of the resonant magnetic field components by the plasma response resulting in much less field line stochasticity. On the other hand, the MARS-F modeling, based on the DINA-simulated post-disruption equilibria, shows that the n = 1 resistive kink instabilities develop in these plasmas, as the edge safety factor qa evolves and drops below integer numbers. RE mitigation by these MHD instabilities is sensitive to the eigenmode structure. The best mitigation is achieved as q a drops below 3, when a global kink instability occurs that encompasses both internal and external components. This global instability is found to be capable of mitigating over 80% MeV-level passing RE orbits at a field perturbation |δB|/B 0 that is comparable to that observed in DIII-D experiments, and full mitigation if the perturbation amplitude is doubled. The 'wetted' area on the ITER limiting surface, due to MHD instability induced RE loss, generally increases with the perturbation amplitude (together with increasing loss fraction). At the highest perturbation level assumed in this study, the wetted area reaches ~60% of the total limiting surface area. Lastly, the lost RE orbits mainly strike the outer divertor region of the limiting surface, with some fraction also hitting a wide area along the inboard side of the surface.

3D fields↗
DOE OSTI2020

Toroidal modeling of runaway electron loss due to 3-D fields in DIII-D and COMPASS

The 3-D field induced relativistic runaway electron (RE) loss has been simulated for DIII-D, COMPASS and ITER plasmas, utilizing the MARS-F code incorporated with the recently developed and updated RE orbit module (REORBIT). Modeling shows effectively 100% loss of a post-disruption, high-current runaway beam in DIII-D, due to the 1 kG level of magnetic field perturbation produced by a fast growing n = 1 resistive kink instability. The RE loss is shown to be independent of the particle energy or the initial location of particles in the configuration space. Applied resonant magnetic perturbation (RMP) fields from in-vessel coils are not effective for RE beam mitigation in DIII-D, but do produce finite (>10%) RE loss in COMPASS, consistent with experimental observations in above two devices. The major reasons for this difference in RE control by RMP between these two devices are (i) the coil proximity to the RE beam and (ii) the effective coil current scaling versus the machine size and the toroidal magnetic field. About 10% RE loss fraction is also predicted for an ITER 15 MA scenario with pre-disruption plasma, highlighting the role of the plasma response. Up to 30% loss is computed, however, by artificially scaling the equilibrium pressure to zero. This is due to the more resistive plasma response and stronger resulting field line stochasticity. Distributions of the lost REs to the limiting surface show poloidally peaked profile near the high-field-side in both DIII-D and COMPASS, covering about 100o poloidal angle. Higher perturbation field level and/or higher particle energy also result in REs being lost to the low-field-side of the limiting surface of these two devices, increasing the effective wetted area. Finally, for the modeled ITER plasmas, the applied RMP field, with optimal poloidal spectrum that maximizes the edge localized mode control, leads to REs being lost to the lower divertor region of the limiting surface within a narrow poloidal band.

3-D field↗
DOE OSTI2024

Electro-Chemo-Mechanical Evolution at the Garnet Solid Electrolyte–Cathode Interface

Solid-state batteries promise higher energy density and improved safety compared with lithium-ion batteries. However, electro-chemomechanical instabilities at the solid electrolyte interface with the cathode and the anode hinder their large scale implementation. Here, in this study, we focus on resolving electro-chemo-mechanical instability mechanisms and their onset conditions between a state-of-the-art cathode, LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC622), and the garnet Li 7 La 3 Zr 2 O 12 (LLZO) solid electrolyte. We used thin-film NMC622 on LLZO pellets to place the interfacial region within the detection depth of the X-ray characterization techniques. The experimental probes of the near-interface region included in operando X-ray absorption spectroscopy and ex situ focused ion beam scanning electron microscopy. Electrochemical degradation was not observable during cycling at room temperature with 4.3 V versus Li/Li + charge voltage cutoff, or with stepwise potentiostatic hold up to 4.1 V versus Li/Li + . In contrast, secondary phases including reduced transition metal species (Ni 2+ , Co 2+ ) were found after cycling up to 4.3 V versus Li/Li + at 80 °C and during potentiostatic hold at 4.3 V versus Li/Li + (Ni 2+ ). Intergranular cracks between NMC622 grains and delamination at the NMC622|LLZO interface occurred readily after the first charge. These interface reaction products and mechanical failure lowered the capacity and cell efficiency due to partial loss of the NMC622 phase, partial loss of contact at the interface, and a higher polarization resistance. Electrochemical instability between delithiated NMC622 and LLZO could be mitigated by using a low charge voltage cutoff or cycling at lower temperature. Ways to engineer the mechanical properties to avoid crack deflection and delamination at the interface are also discussed for enhancing mechanical stability.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Impurity radiation seeding of neoclassical tearing mode growth

Abstract The physics of neoclassical tearing modes (NTMs) is of great concern to tokamak plasma stability and performance, especially in the burning plasma regime. Whereas in many situations the seed events can be clearly identified, such as sawteeth and edge localized modes, the potential seeding mechanism of NTMs due to the resistive tearing instability driven by the impurity radiation cooling still needs more study. Recent NIMROD simulations have demonstrated that local impurity radiation cooling can drive the seed island growth and trigger the subsequent onset of NTM instability. The seed island is mainly driven by the local helical perturbation of the diamagnetic current induced by the perturbed pressure gradient as a result of the impurity radiative cooling on the rational surface. A heuristic closure for the neoclassical viscosity is adopted, and the seed island is further driven by the perturbed bootstrap current induced from the neoclassical electron viscous stress in the extended Ohm’s law. The growth rate of the NTM in simulations is found proportional to the electron neoclassical viscosity, and a theoretical neoclassical driving term is adopted to account for the nonlinear neoclassical island growth in the simulations.

Zeng, Shiyong (ORCID:000000017098818X)↗
DOE OSTI2022

Long-Term Stability Metrics of Photoelectrochemical Water Splitting

Photoelectrochemical (PEC) water splitting, one of the most promising technologies for clean hydrogen generation, has drawn considerable attention over the past few decades. Achieving simultaneous highly efficient and stable unassisted PEC water splitting has been the “holy grail” in clean and renewable fuel generation. State-of-the-art photoelectrodes have shown relatively high efficiencies (∼10–20%). Still, their stability is limited due to photoelectrode chemical instability, electrolyte resistance, mass transfer issues, and an often unoptimized experimental setup. In this work, we present a framework and a set of protocols for conducting long-term stability experiments and further provide details on several critical factors such as light source calibration, choosing the right counter electrode, the configuration of the PEC cell, and photoelectrode sample preparation.

08 HYDROGEN↗
NASA NTRS1994

Influence of plasma sheet thickening toward the tail flanks on magnetotail stability and dynamics

On the basis of resistive MHD simulations we investigate the stability and dynamic evolution of the magnetotail for configurations with various degrees of plasma sheet thickening from midnight toward the tail flanks, associated with a similar increase of the B(sub z) component. This increase is varied by factors between 2 and 3.5, while the magnitude of B(sub z) at midnight is left unchanged. The increase of B(sub z) has a strong effect on the spatial extend of resistive tearing and reconnection across the tail and on the unstable evolution. For stronger thickening the cross-tail extent of the near-Earth neutral line, formed by reconnection, and correspondingly the amount of reconnected magnetic flux get reduced and the dipolarization effects earthward of the neutral line are more concentrated near midnight. An increase of B(sub z) from midnight toward the flanks by more than a factor of about 3 can possibly even suppress the resistive tearing instability. This indicates the possibility of a stability transition due to a reduction of B(sub z) in the flank regions in response to changes in the solar wind, even when the magnitude of B(sub z) near midnight and the fluctuation level in that region are unchanged (provided that a sufficient fluctuation level exists).

Birn, Joachim↗
NASA NTRS2001

Plasma Waves Observed in the Cusp Turbulent Boundary Layer: An Analysis of High Time Resolution Wave and Particle Measurements from the Polar Spacecraft

The boundary layer located in the cusp and adjacent to the magnetopause is a region that is quite turbulent and abundant with waves. The Polar spacecraft's orbit and sophisticated instrumentation are ideal for studying this region of space. Our analysis of the waveform data obtained in this turbulent boundary layer shows broadband magnetic noise extending up to a few kilohertz (but less than the electron cyclotron frequency); sinusoidal bursts (a few tenths of a second) of whistler mode waves at around a few tens of hertz, a few hundreds of hertz, and just below the electron cyclotron frequency; and bipolar pulses, interpreted as electron phase-space holes. In addition, bursts of electron cyclotron harmonic waves are occasionally observed with magnetic components. We show evidence of broadband electrostatic bursts covering a range of approx. 3 to approx. 25 kHz (near but less than the plasma frequency) occurring in packets modulated at the frequency of some of the whistler mode waves. On the basis of high time resolution particle data from the Polar HYDRA instrument, we show that these bursts are consistent with generation by the resistive medium instability. The most likely source of the whistler mode waves is the magnetic reconnection site closest to the spacecraft, since the waves are observed propagating both toward and away from the Earth, are bursty, which is often the case with reconnection, and do not fit on the theoretical cold plasma dispersion relation curve.

Pickett, J. S.↗
NASA NTRS2018

Generation Mechanism for Interlinked Flux Tubes on the Magnetopause

We use a global magnetohydrodynamics simulation to analyze transient magnetic reconnection processes at the magnetopause. The solar wind conditions have been kept constant, and an interplanetary magnetic field with large duskward BY and southward BZ components has been imposed. Five flux transfer events (FTEs) with clear bipolar magnetic field signatures have been observed. We observed a peculiar structure defined as interlinked flux tubes (IFTs) in the first and fourth FTE, which had very different generation mechanisms. The first FTE originates as an IFTs and remains with this configuration until its final moment. However, the fourth FTE develops as a classical flux rope but changes its 3-D magnetic configuration to that of IFTs. This work studies the mechanism for generating IFTs. The growth of the resistive tearing instability has been identified as the cause for the first IFTs formation. We believe that the instability has been triggered by the accumulation of interplanetary magnetic field at the subsolar point where the grid resolution is very high. The evidence shows that two new reconnection lines form northward and southward of the subsolar region. The IFTs have been generated with all the classical signatures of a single flux rope. The other IFTs detected in the fourth FTE developed as a result of magnetic reconnection inside its complex and twisted magnetic fields, which leads to a change in the magnetic configuration from a flux rope of twisted magnetic field lines to IFTs.

Farinas Perez, G.↗
NASA NTRSDate not supplied

Three Dimensional Instantaneous Spray Measurements of Additively Manufactured Candidate Injector Schemes at Varying Pressure Drops

Additively manufactured (AM) injectors have the potential to improve upon the combustion performance achieved using traditionally manufactured (TM) injectors. Traditionally manufactured injector combustion performance decreases greatly in off-design conditions. Additively manufactured injectors hold the potential to sustain high combustion performance over a broader range of operating conditions than currently seen in TM injectors. Additively manufactured injectors promise to offer higher resistances to instability and achieve better overall propellant mixing and atomization for all types of combustion devices. This study expands upon previous studies that used 2-D imaging techniques to characterize the injector spray pattern of candidate AM injector elements. In this study, 3-D instantaneous measurements of the injector spray pattern are taken, allowing for deeper insight into the liquid spray breakup, an inherently 3-D phenomena. Stereo imaging is used for three dimensional analysis, which utilizes two cameras at different angles to reconstruct a three dimensional reconstruction. With this information, AM injectors can be designed to improve propellant mixing and atomization and reduce breakup length. As part of this work, a high-pressure, cold flow experimental spray facility was modified to measure each injector’s spray characteristics in terms of 𝑅𝑒 𝐿 and 𝑊𝑒 𝑔 . This facility is operated to capture instantaneous stereo photographic images, mass flow rates, and injector element thrust for pressure drop across the injector face ranging from 100 to 200 psig at a stiffness ratio of 1.0.

Murphy Mitchell↗
DOE OSTI2020

Stable Glasses of Organic Semiconductor Resist Crystallization

The instability of glassy solids poses a key limitation to their use in several technological applications. Well-packed organic glasses, prepared by physical vapor deposition (PVD), have drawn attention recently because they can exhibit significantly higher thermal and chemical stability than glasses prepared from more traditional routes. We show here that PVD glasses can also show enhanced resistance to crystallization. By controlling the deposition temperature, resistance toward crystallization can be enhanced by at least a factor of ten in PVD glasses of the model organic semiconductor Alq3 (tris(8-hydroxyquinolinato) aluminum). PVD glasses of Alq3 first transform into a supercooled liquid before crystallizing. By controlling the deposition temperature, we increase the glass → liquid transformation time thereby also increasing the overall time for crystallization. Furthermore, we thus demonstrate a new strategy to stabilize glasses of organic semiconductors against crystallization, which is a common failure mechanism in organic light emitting diode devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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.↗
NASA NTRS1974

Solar Physics - Plasma Physics Workshop

A summary of the proceedings of a conference whose purpose was to explore plasma physics problems which arise in the study of solar physics is provided. Sessions were concerned with specific questions including the following: (1) whether the solar plasma is thermal or non-themal; (2) what spectroscopic data is required; (3) what types of magnetic field structures exist; (4) whether magnetohydrodynamic instabilities occur; (5) whether resistive or non-magnetohydrodynamic instabilities occur; (6) what mechanisms of particle acceleration have been proposed; and (7) what information is available concerning shock waves. Very few questions were answered categorically but, for each question, there was discussion concerning the observational evidence, theoretical analyses, and existing or potential laboratory and numerical experiments.

Baum, P. J.↗