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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Magnetic reconnection in 3D fusion devices: non-linear reduced equations and linear current-driven instabilities

Abstract Magnetic reconnection in 3D fusion devices is investigated. With the use of Boozer co-ordinates, we reduce the non-linear resistive magnetohydrodynamic equations in the limit of large aspect ratio and finite pressure fluctuations, to obtain a set of non-linear equations suitable for magnetic reconnection studies in stellarators. Magnetic flux unfreezing due to a finite electron mass is also considered. Equations that govern the linear regime and some of their general properties are given. We emphasise the role of magnetic geometry and identify how some aspects of stellarator optimisation could have an impact on reconnecting instabilities, in particular by exacerbating those enabled by electron inertia. The effect of 3D coupling on the linear reconnection rates and the mode structure is quantitatively addressed in the case in which the equilibrium rotational transform has one specific resonant location for which one mode can reconnect while coupled to an arbitrary number of non-resonant harmonics. The full problem is rigorously reduced to an equivalent cylindrical one, by introducing some geometrically modified plasma inertial and dissipative scales. The 3D scalings for the growth rates of reconnection instabilities and their destabilisation criteria are given.

Physics↗
DOE OSTI2023

Transition from retrograde to prograde drift instabilities in a magnetron microdischarge

Rotating plasma structures, or “spokes,” in magnetized discharges characterized by perpendicular electric and magnetic fields have been seen in a growing number of studies and are believed to be the result of gradient-driven drift instabilities. Under certain conditions, we have shown the spoke’s rotation to be opposite to the E B direction, i.e., retrograde in its expected direction. Recently, we have linked such counter-intuitive rotation to a local inversion of the electric field. Here, we give further experimental evidence for this inversion and attempt to provide an explanation for a relatively distinct transition seen between retrograde and positive (prograde) drift. In the experiments, a partially magnetized plasma forms inside a magnetron device of mm radius operated with argon. Discharge current–voltage measurements are acquired for a range of argon fill pressure and inter-electrode spacing. We find two branches of operation—a low current branch of negative resistance, coinciding with the retrograde spoke rotation, and a higher current branch of positive resistance, coincident with prograde spoke rotation. We postulate that at low discharge currents, high magnetic field confinement leads to a large density gradient, causing more electron transport to the anode than that demanded by the external circuit. At higher currents, anomalous axial electron transport (across the magnetic field lines) becomes dominant, eliminating the conditions for field inversion. Here, the current thresholds for the field inversion are found to be sensitive to the argon pressure and inter-electrode spacing.

42 ENGINEERING↗
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 OSTI2022

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The impact of magnetic fields on momentum transport and saturation of shear-flow instability by stable modes

The Kelvin–Helmholtz (KH) instability of a shear layer with an initially uniform magnetic field in the direction of flow is studied in the framework of 2D incompressible magnetohydrodynamics with finite resistivity and viscosity using direct numerical simulations. The shear layer evolves freely, with no external forcing, and thus broadens in time as turbulent stresses transport momentum across it. As with hydrodynamic KH, the instability here features a conjugate stable mode for every unstable mode in the absence of dissipation. Stable modes are shown to transport momentum up its gradient, shrinking the layer width whenever they exceed unstable modes in amplitude. In simulations with weak magnetic fields, the linear instability is minimally affected by the field, but enhanced small-scale fluctuations relative to the hydrodynamic case are observed. These enhanced fluctuations coincide with increased energy dissipation and faster layer broadening, with these features more pronounced in simulations with stronger fields. These trends result from the magnetic field reducing the effects of stable modes relative to the transfer of energy to small scales. As field strength increases, stable modes become less excited, thus transporting less momentum against its gradient. Furthermore, the energy that would otherwise transfer back to the driving shear because of the stable modes is instead allowed to cascade to small scales, where it is lost to dissipation. Approximations of the turbulent state in terms of a reduced set of modes are explored. While the Reynolds stress is well-described using just two modes per wavenumber at large scales, the Maxwell stress is not.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Local structure effects of carbon-doping on the phase change material Ge 2 Sb 2 Te 5

Ge 2 Sb 2 Te 5 is used in phase change memory, a nonvolatile memory technology, due to its phase change properties. The primary advantage of phase change memory over the state-of-the-art (flash memory) is its simple and small device geometry, which allows for denser nodes and lower power consumption. In phase change memory, resistive heating induces fast switching between the high resistance amorphous and low resistance crystalline phases, corresponding to storage of low and high digital states, respectively. However, the instability of the amorphous phase of Ge 2 Sb 2 Te 5 presents issues with processing and long-term data storage; such issues can be resolved by C doping, which stabilizes the amorphous phase and raises the crystallization temperature. To better understand the local structural effects of C doping on Ge 2 Sb 2 Te 5 , in situ Ge K-edge X-ray absorption spectroscopy measurements were taken during heating of films with various C doping concentrations. Here, the range of structural transformation temperatures derived from X-ray absorption near-edge structure analysis across the C doping series proved narrower than crystallization temperatures reported in similar in situ X-ray diffraction experiments, which may reflect changes in local structure that precede long-range ordering during crystallization. In addition, rigorous extended X-ray absorption fine structure fitting across and between temperature series revealed effects of C doping on the rigidity of Ge–Te bonds at low (2 at% and 4 at%) C concentrations.

Langhout, John D.↗
DOE OSTI2023

Ideal MHD induced temperature flattening in spherical tokamaks

This paper extends the analysis first presented in Jardin et al. [Phys. Rev. Lett. 128, 245001 (2022)] to more thoroughly examine the stability of spherical torus equilibrium to ideal magnetohydrodynamic (MHD) infernal modes and their nonlinear consequences. Here we demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation of a NSTX discharge, anomalous transport can occur due to these instabilities. Here we generate a family of equilibrium of differing β and use this to show that these instabilities could explain the experimentally observed flattening of the electron temperature profile at modest β. The modes studied in this paper are found to occur with poloidal mode number m and toroidal mode number n when the ratio m/ n is in the range of 1.2–1.5, when the central safety factor is in this range or slightly lower, and when the central region has very low magnetic shear. Our analysis gives some insight as to why the unstable linear growth rates are oscillatory functions of the toroidal mode number n. We present a simulation of an initially stable configuration that passes through a stability boundary at a critical β as it is heated. We also show that a particular NSTX discharge is unstable to these modes over a timescale of several hundred ms. We conclude that these modes must be taken into account when performing predictive modeling. An appendix shows that similar modes can be found in R / a = 4 tokamaks for certain q-profiles and β values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Instability and turbulent relaxation in a stochastic magnetic field

An analysis of instability dynamics in a stochastic magnetic field is presented for the tractable case of the resistive interchange. Externally prescribed static magnetic perturbations convert the eigenmode problem to a stochastic differential equation, which is solved by the method of averaging. The dynamics are rendered multi-scale, due to the size disparity between the test mode and magnetic perturbations. Maintaining quasi-neutrality at all orders requires that small-scale convective cell turbulence be driven by disparate scale interaction. Here, the cells in turn produce turbulent mixing of vorticity and pressure, which is calculated by fluctuation-dissipation type analyses, and are relevant to pump-out phenomena. The development of correlation between the ambient magnetic perturbations and the cells is demonstrated, showing that turbulence will 'lock on' to ambient stochasticity. Magnetic perturbations are shown to produce a magnetic braking effect on vorticity generation at large scale. Detailed testable predictions are presented. The relations of these findings to the results of available simulations and recent experiments are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Lattice instabilities and phonon thermal transport in TlBr

Materials with simple crystal structure usually exhibit substantial lattice thermal conductivity ( κ latt ) due to a scarcity of strong phonon-phonon scattering channels. Yet here we find simple CsCl-structure TlBr to have an extremely low room temperature κ latt ~ 0.5 W / m-K , while at low temperature we find a high κ latt of 155 W / m-K . This extreme range of conductivity behavior results from lattice-instability-related anharmonicity, effective only at temperatures above ~ 30 K . Manipulating these instabilities (e.g., pressure, strain) may yield a pathway to ultralow thermal conductivity via intrinsic phonon resistance in similar materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Time-resolved X-ray imaging of the current filamentation instability in solid-density plasmas

The streaming of energetic charged particles can magnetize astrophysical and laboratory plasmas via the current filamentation instability. Despite its importance, the experimental characterization of this instability has remained a challenge. Here, we report an experiment combining a high-intensity optical laser with a high-brightness X-ray free electron laser that successfully images the instability in solid-density plasmas with 200 nm spatial and 50 fs temporal resolution. We characterize the development of μm-scale filamentary structures and their evolution over tens of picoseconds through a non-linear merging process. The measured plasma density modulations and long merging time reveal the critical importance of space-charge effects and ion motion on this electron-driven instability. Supporting theoretical analysis and kinetic simulations help distinguish the relative role of space-charge and resistive effects. Our findings indicate that magnetic fields on the order of 10 megagauss are produced, with important implications for transport and radiation emission of energetic particles in plasmas.

Schoenwaelder, Christopher [SLAC National Accelera↗
DOE OSTI2026

Time-resolved X-ray imaging of the current filamentation instability in solid-density plasmas

The streaming of energetic charged particles can magnetize astrophysical and laboratory plasmas via the current filamentation instability. Despite its importance, the experimental characterization of this instability has remained a challenge. Here, we report an experiment combining a high-intensity optical laser with a high-brightness X-ray free electron laser that successfully images the instability in solid-density plasmas with 200nm spatial and 50fs temporal resolution. We characterize the development of µm-scale filamentary structures and their evolution over tens of picoseconds through a non-linear merging process. The measured plasma density modulations and long merging time reveal the critical importance of space-charge effects and ion motion on this electron-driven instability. Supporting theoretical analysis and kinetic simulations help distinguish the relative role of space-charge and resistive effects. Our findings indicate that magnetic fields on the order of 10megagauss are produced, with important implications for transport and radiation emission of energetic particles in plasmas.

Schoenwaelder, Christopher↗
DOE OSTI2020

Time- and strain-dependent nanoscale structural degradation in phase change epitaxial strontium ferrite films

Topotactic phase transition between metallic, perovskite SrFeO 3 and insulating, Brownmillerite SrFeO 2.5 has been extensively studied due to the potential applications in resistive switching devices for neuromorphic computing. However, its practical utilization as memristors has been hindered by the structural instability of SrFeO 3 , which is often ascribed to the generation of oxygen vacancies to form SrFeO 3-δ . Here we reveal the dominating defects generated in SrFeO 3 epitaxial thin films are atomic scale gaps as a result of interfacial strain. Our correlated time- and strain-dependent measurements show that tensile strained SrFeO 3 films form vertical, nanoscale gaps that are SrO-rich, which are accountable for the observed metal-to-insulator transition over time. While compressively strained or small lattice mismatched SrFeO 3 films mainly yield horizontal gaps with a smaller impact on the in-plane transport. The atomic scale origin of such defects and their impact on device performance need to be further understood in order to integrate phase change materials in oxide electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High-Performance Layered Ni-Rich Cathode Materials Enabled by Stress-Resistant Nanosheets

Here, layered O3-type transition metal oxides are promising cathode candidates for high-energy-density Li-ion batteries. However, the structural instability at the highly delithiated state and low kinetics at the fully lithiated state are arduous challenges to overcome. Here, a facile approach is developed to make secondary particles of Ni-rich materials with nanosheet primary grains. Because the alignment of the primary grains reduces internal stress buildup within the particle during charge-discharge and provides straightforward paths for Li transport, the as-synthesized Ni-rich materials do not undergo cracking upon cycling with higher overall Li + ion diffusion rates. Specifically, a LiNi 0.75 Co 0.14 Mn 0.11 O 2 cathode with nanosheet grains delivers a high reversible capacity of 206 mAh g -1 and shows ultrahigh cycling stability, e.g., 98% capacity retention over 500 cycles in a full cell with a graphite anode.

25 ENERGY STORAGE↗
DOE OSTI2021

Simulations of heat fluxes in an ELMy H-mode discharge on HL-2A

In order to study the distribution and evolution of the transient heat flux on HL-2A during edge-localized-mode (ELM) bursts, the BOUT++ electromagnetic six-field two-fluid model is used to simulate the pedestal collapse under the lower single-null divertor geometry. The equilibrium profiles of HL-2A ELMy H-mode discharge No. 24 953 are adopted as the initial condition in the original case. In this instance, linear analysis shows that the resistive ballooning mode (RBM) and drift-Alfven wave are unstable to this equilibrium, and RBM is the dominant instability. The evolutions of the radial heat fluxes at the outer mid-plane and heat fluxes to the inner and outer targets during the ELM event are presented. Six more equilibria are constructed based on the original case to find out the influence of the pedestal profiles on the peak electron heat flux. The results indicate that the heat flux increases with temperature and/or density, and the theoretical analysis and simulation results consistently show that the heat flux q ∥e is proportional to n e0,SEP $T^{\frac32}_{e0,SEP}$.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ultrathin stable Ohmic contacts for high-temperature operation of β -Ga2O3 devices

Beta gallium oxide (β-Ga2O3) shows significant promise in high-temperature, high-power, and sensing electronics applications. However, long-term stable metallization layers for Ohmic contacts at high temperatures present unique thermodynamic challenges. The current most common Ohmic contact design based on 20 nm of Ti has been repeatedly demonstrated to fail at even moderately elevated temperatures (300–400 °C) due to a combination of nonstoichiometric Ti/Ga2O3 interfacial reactions and kinetically favored Ti diffusion processes. Here, we demonstrate stable Ohmic contacts for Ga2O3 devices operating up to 500–600 °C using ultrathin Ti layers with a self-limiting interfacial reaction. The ultrathin Ti layer in the 5 nm Ti/100 nm Au contact stack is designed to fully oxidize while forming an Ohmic contact, thereby limiting both thermodynamic and kinetic instability. This novel contact design strategy results in an epitaxial conductive anatase titanium oxide interface layer that enables low-resistance Ohmic contacts that are stable both under long-term continuous operation (>500 h) at 600 °C in vacuum (≤10−4 Torr), as well as after repeated thermal cycling (15 times) between room temperature and 550 °C in flowing N2. This stable Ohmic contact design will accelerate the development of high-temperature devices by enabling research focus to shift toward rectifying interfaces and other interfacial layers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

W7-X and the sawtooth instability: towards realistic simulations of current-driven magnetic reconnection

Magnetic reconnection in W7-X is studied by means of global numerical simulations in a series of models of increasing complexity. The magnetic geometry ranges from that of a cylinder to the full three-dimensional field of W7-X, and the equations solved range from ideal magnetohydrodynamics (MHD) to gyrokinetics. We simulate plasmas from the first operation phase with electron cyclotron current drive (ECCD). These are characterized by an equilibrium magnetic field featuring an ECCD-distorted 'humped' profile of the rotational transform ι, with ι = 1 in two radial locations. Such plasmas generally show sawtooth activity, hence motivating the present study. We pay particular attention to the role of equilibrium current density gradients in the destabilization of reconnecting modes. When the equilibrium temperature and density gradients are artificially suppressed (to eliminate the pressure gradient drive), the perturbed electrostatic potential is radially localized between the locations at which ι = 1. This is shown with a purely collisionless gyrokinetic model, in cylindrical geometry. In the real toroidal geometry of W7-X, for a non-ideal MHD model including a uniform resistivity, electron inertia and (numerical) viscosity, the same qualitative behaviour is observed. In particular, even if a resonant (m, n) = (1, -1) perturbation is initialized, the most unstable mode is the (m, n) = (-4, 4), where m and n are the poloidal and toroidal mode numbers, respectively. Other modes are destabilized due to geometric coupling. The growth rate of this instability scales as η 1/3 , where η is the plasma resistivity, thus suggesting that ECCD drives ideal MHD stable W7-X plasmas towards non-ideal marginality. An ideal magnetohydrodynamic analysis confirms the result. A fluid-kinetic hybrid version of the EUTERPE code shows that gyrokinetic ions have a stabilizing effect on these modes. For W7-X relevant collisionalities, the growth rate scales linearly with the electron skin depth, d e . Implications of our results for sawtoothing W7-X operation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Multifunctional Effect of Fe Substitution in Na Layered Cathode Materials for Enhanced Storage Stability

In this study, developing stable cathode materials that are resistant to storage degradation is essential for practical development and industrial processing of Na-ion batteries, as many sodium layered oxide materials are susceptible to hygroscopicity and instability when exposed to ambient air. Among the various layered compounds, Fe-substituted O 3 -type Na(Ni 1/2 Mn 1/2 ) 1-x Fe x O 2 materials have emerged as a promising option for high-performance and low-cost cathodes. While previous reports have noted the decent air-storage stability of these materials, the role and origin of Fe substitution in improving storage stability remain unclear. In this study, we investigate the air resistant effect of Fe substitution in O 3 -Na(Ni 1/2 Mn 1/2 ) 1-x Fe x O 2 cathode materials by performing systematic surface and structural characterizations. We find that the improved storage stability can be attributed to the multifunctional effect of Fe substitution, which forms a surface protective layer containing an Fe-incorporated spinel phase and decreases the thermodynamical driving force for bulk chemical sodium extraction. With these mechanisms, the Fe-containing cathodes can suppress the cascades of cathode degradation processes and better retain the electrochemical performance after air storage.

25 ENERGY STORAGE↗
DOE OSTI2026

Experimental considerations and applications for time-resolved impedance on silicon-anode battery systems

Dynamic electrochemical impedance spectroscopy (dEIS) utilizes a superimposed multi-sinusoidal waveform and enables temporally resolved investigations on various electrochemical processes; however, the signal-to-noise ratio must be maximized while retaining linearity and stationarity within the system. In the present work, we probe the impact of the waveform’s total amplitude as well as the waveform and phase profiles. The equal amplitude multi-sinusoidal waveform with 35 mV total amplitude (<1 mV per frequency) and random phase profile exhibits the lowest impedance noise. This input signal is sufficient in avoiding the high frequency signal attenuation from the potentiostat’s low-pass filter. Additionally, we demonstrate the optimized dEIS waveform’s utility in investigating surface passivation within thin-film and composite silicon (Si) anodes. Within a two-electrode Si-lithium metal coin cell, the charge transfer resistance (RCT) associated with alloying kinetics dominates the overall impedance. RCT increases with state of charge indicating a kinetic bottleneck at higher lithiation states. To explore this further, we utilize a three-electrode Si-LiNi0.8Mn0.1Co0.1O2 pouch cell to deconvolute the Si impedance contributions. The evolution of the solid electrolyte interphase (SEI) and charge transfer resistances correlate to specific Li–Si alloy phases. We attribute this to the influence of Si particle volume expansion on SEI instabilities and a kinetic bottleneck at higher lithiation states.

Lam, Steven [ORNL]↗