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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Non‐Linear Kinetics of The Lithium Metal Anode on Li 6 PS 5 Cl at High Current Density: Dendrite Growth and the Role of Lithium Microstructure on Creep

Abstract Interfacial instability, viz., pore formation in the lithium metal anode (LMA) during discharge leading to high impedance, current focusing induced solid–electrolyte (SE) fracture during charging, and formation/behaviour of the solid–electrolyte interphase (SEI), at the anode, is one of the major hurdles in the development of solid‐state batteries (SSBs). Also, understanding cell polarization behaviour at high current density is critical to achieving the goal of fast‐charging battery and electric vehicle. Herein, via in situ electrochemical scanning electron microscopy (SEM) measurements, performed with freshly deposited lithium microelectrodes on transgranularly fractured fresh Li6PS5Cl (LPSCl), the LiǀLPSCl interface kinetics are investigated beyond the linear regime. Even at relatively small overvoltages of a few mV, the LiǀLPSCl interface shows non‐linear kinetics. The interface kinetics possibly involve multiple rate‐limiting processes, i.e., ion transport across the SEI and SE|SEI interfaces, as well as charge transfer across the LiǀSEI interface. The total polarization resistance R P of the microelectrode interface is determined to be ≈ 0.8 Ω cm 2 . It is further shown that the nanocrystalline lithium microstructure can lead to a stable LiǀSE interface via Coble creep along with uniform stripping. Also, spatially resolved lithium deposition, i.e., at grain surface flaws, grain boundaries, and flaw‐free surfaces, indicates exceptionally high mechanical endurance of flaw‐free surfaces toward cathodic load (>150 mA cm −2 ). This highlights the prominent role of surface defects in dendrite growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Wide-Bandgap Semiconductor Amplifiers for Fusion Plasma Heating and Control

This paper discusses power electronics developed under the ARPA-E GAMOW program to support nuclear fusion power production. The goal of this project was to develop and assess the potential for wide-bandgap (WBG) semiconductor devices in power electronics to enable high-efficiency and high-voltage solid-state systems for fusion plasma generation, heating, and control. The power electronics use an architecture in which multiple high-power boards can be combined to produce megawatt-level power, where using multiple boards provides high reliability. Two main areas of power electronics boards are developed in this project for fusion plasma heating and control applications: (1) pulse generation and control and (2) radiofrequency generation. The first area is for boards capable of driving high-voltage millisecond pulses at high duty cycles. The envisioned application of these pulses is in plasma control of magnetohydrodynamic instabilities, plasma position, and edge-localized modes. Pulse-width modulation allows for the implementation of a wide variety of linear and nonlinear control systems. The boards developed for this project could actuate control coils based on digital input signals and can be parallelized to provide megawatts of output power. The design of the pulse generator is a low-side load switch. A load switch was designed and constructed that utilized 2-kV-rated field-effect transistor (FET)-based cascodes developed by Qorvo under this project to perform initial testing of these cascodes. The second area is being implemented using class E amplifiers with WBG devices and a reactance steering network to handle inductive or capacitive plasma loads. Applications include ion cyclotron resonance heating (ICRH) and high-harmonic fast-wave (HHFW) heating. A class E reactance steering network is demonstrated in modeling and experiment with a resistive-inductive load that models an inductively-coupled plasma. Power combining of boards with class E reactance steering networks is also simulated and demonstrated experimentally, to enable scaling up to high power. Modeling of high-power-density cooling and remaining useful life is conducted to enable reliable, effectively cooled high-power electronics for fusion applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Local measurements of the spatial magnetic field distribution in a z-pinch plasma during and near stagnation using polarization spectroscopy

We present here the detailed measurements of radial distribution of the magnetic field in a gas-puff z-pinch plasma at the final stages of the implosion phase and at stagnation. While the measurements are chordal, the radial distribution of different charge states was utilized to measure the magnetic field locally for certain radii, so that unlike chordal measurements in general, the magnetic field radial distribution was obtained with no need for the Abel inversion of the data. The distribution was measured using the Zeeman effect via a novel spectroscopic technique, at several axial locations, and demonstrates striking features such as the peak field remaining at a radius much larger than the stagnation radius at all times. Furthermore, while the distribution observed is sometimes monotonic with respect to the radius, it is often not, a behavior that can be linked to 2D features in the plasma column resulting from the Rayleigh–Taylor instability. The current flowing through the stagnating plasma was found to be a small fraction of the total current, resulting in clearly insufficient magnetic pressure to balance the plasma pressure at stagnation. The magnetic field data, taken over several axial positions, are used to obtain the true inductance in the imploding plasma for the first time; it is found that the data cannot explain the current turnover at stagnation. A simulation with the MACH2-Tabular Collisional-Radiative Equilibrium magnetohydrodynamics code in the r–z plane shows that the peak of the magnetic field pinches to a much smaller radius than is observed in the spectroscopic data. Furthermore, the depth of the computed current turnover at stagnation is smaller than the measured one. The two observed features of a radially extended magnetic field at stagnation together with a deep current turnover are a challenge to match in simulations. Various calculations and estimates of the inductive and resistive load voltages are examined to ascertain if they are responsible for the observed current notch. Lastly, the results demonstrate that the knowledge of the true inductance in the driven load requires such magnetic-field-distribution measurements and that imaging data or electrical measurements are insufficient.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Mn(Pt 1–x Pd x ) 5 P: Isovalent tuning of Mn-sublattice magnetic order

Here, e report the growth and characterization of MnPd 5 P, a rare-earth-free ferromagnet, with T C ≈ 295 K and planar anisotropy, and conduct a substitutional study with its antiferromagnetic analogue MnPt 5 P. We provide a solution route to grow large single crystals of MnPd 5 P and the series Mn(Pt 1–x Pd x ) 5 P by adding Mn into Pd-P and (Pt 1–x Pd x )-P based melts. All compounds in the family adopt the layered anti-CeCoIn5 type structure with the space group P4/mmm, and EDS and X-ray diffraction results indicate that MnPt 5 P and MnPd 5 P form a complete solid solution. Based on measurements of the temperature- and field-dependent magnetization and resistance, we construct a temperature-composition (T-x) phase diagram for Mn(Pt 1–x Pd x ) 5 P and demonstrate that the initial antiferromagnetic order found in MnPt 5 P is extraordinarily sensitive to Pd substitution. At low Pd fractions (x < 0.010), the single antiferromagnetic transition in pure MnPt 5 P splits into a higher temperature ferromagnetic transition followed first, upon cooling, by a lower temperature ferromagnetic to antiferromagnetic transition and then by a re-entrant antiferromagnetic to ferromagnetic transition at even lower temperatures. The antiferromagnetic region makes up a bubble phase that persists up to x ≈ 0.008-0.009 for T ≈ 150 K, with all samples x < 0.008 recovering their initial ferromagnetic state upon further cooling to base temperature. Over the same low substitution range we find a non-monotonic change in the room temperature value of the unit cell volume, further suggesting that pure MnPt 5 P is very close to an instability. Once x > 0.010, Mn(Pt 1–x Pd x ) 5 P undergoes a only single transition into the ferromagnetic phase. The Curie temperature initially increases rapidly with x, rising from T C ≈ 197 K at x = 0.013 to a maximum of T C ≈ 312 K for x ≈ 0.62, and then falling back to T C ≈ 295 K for pure MnPd 5 P (x = 1.00). Given that Pt and Pd are isoelectronic, this work raises questions as to the origin of the extreme sensitivity of the magnetic ground state in MnPt 5 P upon introducing Pd.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Advanced Nanocarbons for Enhanced Performance and Durability of Platinum Catalysts in Proton Exchange Membrane Fuel Cells

Insufficient stability of current carbon supported Pt and Pt alloy catalysts is a significant barrier for proton-exchange membrane fuel cells (PEMFCs). As a primary degradation cause to trigger Pt nanoparticle migration, dissolution, and aggregation, carbon corrosion remains a significant challenge. Compared with enhancing Pt and PtM alloy particle stability, improving support stability is rather challenging due to carbon's thermodynamic instability under fuel cell operation. In recent years, significant efforts have been made to develop highly durable carbon-based supports concerning innovative nanostructure design and synthesis along with mechanistic understanding. Here, this review critically discusses recent progress in developing carbon-based materials for Pt catalysts and provides synthesis–structure–performance correlations to elucidate underlying stability enhancement mechanisms. The mechanisms and impacts of carbon support degradation on Pt catalyst performance are first discussed. The general strategies are summarized to tailor the carbon structures and strengthen the metal–support interactions, followed by discussions on how these designs lead to enhanced support stability. Based on current experimental and theoretical studies, the critical features of carbon supports are analyzed concerning their impacts on the performance and durability of Pt catalysts in fuel cells. Finally, the perspectives are shared on future directions to develop advanced carbon materials with favorable morphologies and nanostructures to increase Pt utilization, strengthen metal-support interactions, facilitate mass/charge transfer, and enhance corrosion resistance.

25 ENERGY STORAGE↗
DOE OSTI2020

Fluctuation dynamo in a weakly collisional plasma

The turbulent amplification of cosmic magnetic fields depends upon the material properties of the host plasma. In many hot, dilute astrophysical systems, such as the intracluster medium (ICM) of galaxy clusters, the rarity of particle–particle collisions allows departures from local thermodynamic equilibrium. These departures – pressure anisotropies – exert anisotropic viscous stresses on the plasma motions that inhibit their ability to stretch magnetic-field lines. We present an extensive numerical study of the fluctuation dynamo in a weakly collisional plasma using magnetohydrodynamic (MHD) equations endowed with a field-parallel viscous (Braginskii) stress. When the stress is limited to values consistent with a pressure anisotropy regulated by firehose and mirror instabilities, the Braginskii-MHD dynamo largely resembles its MHD counterpart, particularly when the magnetic field is dynamically weak. If instead the parallel viscous stress is left unabated – a situation relevant to recent kinetic simulations of the fluctuation dynamo and, we argue, to the early stages of the dynamo in a magnetized ICM – the dynamo changes its character, amplifying the magnetic field while exhibiting many characteristics reminiscent of the saturated state of the large-Prandtl-number ( ${Pm}\gtrsim {1}$ ) MHD dynamo. We construct an analytic model for the Braginskii-MHD dynamo in this regime, which successfully matches simulated dynamo growth rates and magnetic-energy spectra. A prediction of this model, confirmed by our numerical simulations, is that a Braginskii-MHD plasma without pressure-anisotropy limiters will not support a dynamo if the ratio of perpendicular and parallel viscosities is too small. This ratio reflects the relative allowed rates of field-line stretching and mixing, the latter of which promotes resistive dissipation of the magnetic field. Finally, in all cases that do exhibit a viable dynamo, the generated magnetic field is organized into folds that persist into the saturated state and bias the chaotic flow to acquire a scale-dependent spectral anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulations of radiatively cooled magnetic reconnection driven by pulsed power

Magnetic reconnection is an important process in astrophysical environments, as it reconfigures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the reconnection layer. In this paper, we perform two- and three-dimensional simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively cooled regime. These simulations are performed in GORGON, an Eulerian two-temperature resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and P 1/3 multi-group radiation transport. The load produces highly collisional, super-Alfvénic (Alfvén Mach number M A ≈ 1.5), supersonic (Sonic Mach number M S ≈ 4–5) strongly driven plasma flows which generate an elongated reconnection layer (Aspect Ratio L/δ≈100 , Lundquist number S L ≈ 400). The reconnection layer undergoes radiative collapse when the radiative losses exceed the rates of ohmic and compressional heating (cooling rate/hydrodynamic transit rate = τ$^{–1}_{cool}$/τ$^{–1}_{H}$ ≈100); this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Lastly, the current sheet is also unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Observations of heteroclinic bifurcations in resistive magnetohydrodynamic simulations of the plasma response to resonant magnetic perturbations

A class of topological magnetic island bifurcations that has not previously been observed in toroidal plasmas is described. Increasing an externally applied three-dimensional magnetic field in resistive magnetohydrodynamic simulations results in the asymmetric elongation of resonant island flux surfaces followed by a sequence of heteroclinic bifurcations. Here, these bifurcations produce new sets of hyperbolic-elliptic fixed points as predicted by the Poincaré-Birkoff fixed point theorem. Field line calculations verify that the new fixed points do not connect to those of the prebifurcated islands as required for heteroclinic bifurcations on a torus with winding numbers composed of common integer factors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Probing the Formation of Cathode-Electrolyte Interphase on Lithium Iron Phosphate Cathodes via Operando Mechanical Measurements

Interfacial instabilities in electrodes control the performance and lifetime of Li-ion batteries. While the formation of solid-electrolyte interface (SEI) on anodes has received much attention, there is still lack of understanding about the formation of cathode-electrolyte interface (CEI) on the cathodes. To fill this gap, we report on dynamic deformations on lithium iron phosphate, LiFePO 4 cathodes during charge / discharge by utilizing in-operando digital image correlation, impedance spectroscopy and Cryo X-ray photoelectron spectroscopy. LiFePO4 cathodes were cycled in either LiPF 6 , LiClO 4 or LiTFSI- containing organic liquid electrolytes. Beyond the first cycle, Li-ion intercalation results in a nearly linear correlation between electrochemical strains and the state of (dis)-charge, regardless of the electrolyte chemistry. However, during the first charge in LiPF 6 - containing electrolyte, there is a distinct irreversible positive strain evolution at the onset of anodic current rise as well as current decay at around 4.0V. Impedance studies show the increase in surface resistance in the same potential window, suggesting the formation of CEI layers on the cathode. The chemistry of the CEI layer was characterized by X-ray photoelectron spectroscopy. LiF is detected in CEI layer starting as early as 3.4 V and Li $x$ PO $y$ F $z$ appeared at voltages higher than 4.0 V during the first charge. In conclusion, our approach offers new insights into the formation mechanism of CEI layers on the cathode electrodes, which is crucial for the development of robust cathode and electrolyte chemistries for higher performance batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

The Madison plasma dynamo experiment: A facility for studying laboratory plasma astrophysics

The Madison plasma dynamo experiment (MPDX) is a novel, versatile, basic plasma research device designed to investigate flow driven magnetohydrodynamic instabilities and other high-β phenomena with astrophysically relevant parameters. A 3 m diameter vacuum vessel is lined with 36 rings of alternately oriented 4000 G samarium cobalt magnets, which create an axisymmetric multicusp that contains ∼14 m3of nearly magnetic field free plasma that is well confined and highly ionized (>50%). At present, 8 lanthanum hexaboride (LaB6) cathodes and 10 molybdenum anodes are inserted into the vessel and biased up to 500 V, drawing 40 A each cathode, ionizing a low pressure Ar or He fill gas and heating it. Up to 100 kW of electron cyclotron heating power is planned for additional electron heating. The LaB6cathodes are positioned in the magnetized edge to drive toroidal rotation through J × B torques that propagate into the unmagnetized core plasma. Dynamo studies on MPDX require a high magnetic Reynolds number Rm > 1000, and an adjustable fluid Reynolds number 10 1). Initial results from MPDX are presented along with a 0-dimensional power and particle balance model to predict the viscosity and resistivity to achieve dynamo action.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Drift kinetic effects on plasma response to resonant magnetic perturbation for EU DEMO design

Abstract A systematic investigation of the plasma response to an externally applied resonant magnetic perturbation (RMP) field, for the purpose of controlling edge localized modes, is carried out for an EU DEMO reference plasma. Particular emphasis is placed on the role of kinetic effects associated with both thermal particles and fusion-born alphas. The single fluid, resistive model predicts a large peak amplification of the n = 1 ( n is the toroidal mode number) plasma response to the target equilibrium, which is found to be close to the Troyon no-wall limit. A more advanced response model, including kinetic resonances between the RMP perturbation and drift motions of thermal and energetic particles, on the other hand, finds a strong suppression of the n = 1 field amplification. A major role is played by the precessional drift resonance of fusion-born alphas. A strong parallel sound wave damping model is found to well reproduce the full kinetic response results for the DEMO plasma, in terms of both the resonant field response amplitude and the plasma displacement. Finally, both fluid and kinetic models produce similar responses for the n = 2 and 3 RMP fields for the considered DEMO plasma, whilst kinetic effects again become important for the n = 4 RMP due to proximity of the reference plasma to the no-wall limit for the n = 4 ideal kink instability.

Physics↗
DOE OSTI2024

Extended-MHD simulations of disruption mitigation via massive gas injection in SPARC

Recent developments to the M3D-C1 code enable higher fidelity modeling of disruptions, and can be applied in the design verification of reactor-scale tokamaks. Among these new capabilities is a method to mesh conducting vessel structures such as coils and passive plates, packing of the toroidal mesh around gas injectors, as well as anisotropic resistivity inside the vessel structures. We present extended-magnetohydrodynamic (MHD) simulations of disruption mitigation via massive gas injection (MGI) in SPARC. The goal of this study is to inform the disruption mitigation layout of SPARC and aid in the design of an effective gas injector configuration. Fully three-dimensional simulations with M3D-C1 are carried out for various injector configurations with the primary goal of determining the effect of different MGI parameters on heat loads and vessel forces. The simulations include a model for impurity ionization, recombination, advection and radiation, as well as spatially resolved conducting structures around the plasma. A localized mixture of deuterium and neon with a small toroidal and poloidal width is injected in up to six locations. We demonstrate that M3D-C1 can model a rapid shutdown via MGI using narrow and more realistic gas plumes than in previous simulations. As a result of the q = 1 surface in the SPARC baseline case a sawtooth is observed early in the simulations. Despite the sawtooth and the onset of edge MHD instabilities, the impurity distribution remains localized around the injector locations, but enables a radiative shutdown of the plasma. We find that using the maximum of six gas injectors results in a lower peaking factor and leads to a more even distribution of radiation toroidally than using two injectors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A novel path to runaway electron mitigation via deuterium injection and current-driven MHD instability

Relativistic electron (RE) beams at high current density (low safety factor, qa) yet very low free-electron density accessed with D2 secondary injection in the DIII-D and JET tokamak are found to exhibit large-scale MHD instabilities that benignly terminate the RE beam. In JET, this technique has enabled termination of MA-level RE currents without measurable first-wall heating. This scenario thus offers an unexpected alternate pathway to achieve RE mitigation without collisional dissipation. Benign termination is explained by two synergistic effects. First, during the MHD-driven RE loss events both experiment and MHD orbit-loss modeling supports a significant increase in the wetted area of the RE loss. Second, as previously identified at JET and DIII-D, the fast kink loss timescale precludes RE beam regeneration and the resulting dangerous conversion of magnetic to RE kinetic energy. During the termination, the RE kinetic energy is lost to the wall, but the current fully transfers to the cold bulk thus enabling benign Ohmic dissipation of the magnetic energy on longer timescales via a conventional current quench. Hydrogenic (D2) secondary injection is found to be the only injected species that enables access to the benign termination. D2 injection: 1) facilitates access to low qa in existing devices (via reduced collisionality & resistivity), 2) minimizes the RE avalanche by ‘purging’ the high-Z atoms from the RE beam, 3) drives recombination of the background plasma, reducing the density and Alfven time, thus accelerating the MHD growth. Furthermore, this phenomenon is found to be accessible when crossing the low qa stability boundary with rising current, falling toroidal field, or contracting minor radius - the latter being the expected scenario for vertically unstable RE beams in ITER. While unexpected, this path scales favorably to fusion-grade tokamaks and offers a novel RE mitigation scenario in principle accessible with the day-one disruption mitigation system (DMS) of ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Superconductivity at carrier density 10 17 cm – 3 in quasi-one-dimensional Li 0.9 Mo 6 O 17

Superconductivity is a fascinating phenomenon that involves an attractive interaction by which electrons are paired and can move without electrical resistance below a critical temperature T c . A paradigm for understanding superconductive pairing, the successful “BCS” theory, developed more than 65 years ago, is challenged by certain classes of materials for which parameters (e.g. the value of T c , the density of electrons, etc.) fall outside the range for which BCS theory is applicable. The study of such materials is important because new physics may be revealed and because superconductivity has many potential applications in sensing and electronics. Superconductivity at very low carrier density is a rare and fascinating phenomenon that is of considerable importance and interest in condensed matter physics. The conditions of high density of states and effective screening of Coulomb repulsion, both key ingredients of stable Cooper pairing in the BCS theory, may not be met at low carrier density. Thus, such systems may also harbor unconventional electron pairing (e.g. non-phononic) or pairing without phase coherence as a precursor to Bose-Einstein condensation. An intriguing route toward extreme low-density (and possibly unconventional) superconductivity may be found in quasi-one-dimensional metals for which density-wave instabilities compete for the ground state and can lead to gapping of substantial portions of the Fermi surface (FS), leaving a residual (reconstructed) FS comprising a very small carrier density. Our work highlights a novel material in this unusual category: Li 0.9 Mo 6 O 17 (T c =2 K) (also known as "lithium purple bronze" or "LiPB"). We present both electrical and thermoelectric magnetotransport measurements that reveal superconductivity to occur at carrier density 2 x10 17 cm -3 , among the lowest known for any superconductor, and four orders of magnitude lower than is expected from its chemical valence and band structure. With its exceptionally large ratio of superconducting to Fermi temperature, T c /TF ≳ 0.1, LiPB is properly placed at the top of a very small group of superconductors (including cuprates, heavy-fermion and Fe-based superconductors) for which this ratio falls in the BCS-BEC crossover regime and for which pairing is believed to be spin mediated.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effects of secondary void-initiating particles on the steady-state crack growth resistance of high-strength steel

In this investigation, we focus on the dependence of stable crack growth resistance, as measured by the short-rod fracture toughness test, on secondary void-initiating particles such as AlN and Ti-based grain-refining precipitates in a variety of high-strength steels with tempered martensitic microstructures. The model developed by Ritchie and Thompson is modified to illustrate the significant amounts of toughening that can result from the refinement of secondary particles. Analysis of the data suggests that material strength is a predominant factor in increasing the short-rod fracture toughness relative to linear-elastic measures of initiation fracture toughness, but the extent of toughening is limited by the size and number density of secondary particles in the microstructure. The variation in estimates of secondary microvoid initiation and growth strains with precipitate size reinforce the notion that primary fracture at non-metallic inclusions and secondary fracture at smaller particles occur as sequential processes with a degree of concurrence that is dependent on the state of precipitation in both particle dispersions. Toughening in this connection is maximized by increases in microvoid growth strain that result from decreases in the size and areal number density of secondary void-initiating particles. Finally, the occurrence of transient instabilities during crack extension in short-rod specimens is explained with a phenomenological model that relates crack growth stability to natural variations in the dispersion of secondary void-initiating particles in the microstructure.

42 ENGINEERING↗
DOE OSTI2020

Axisymmetric simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

A benchmark exercise for the modeling of vertical displacement events (VDEs) is presented and applied to the 3D nonlinear magnetohydrodynamic codes M3D-C 1 , JOREK, and NIMROD. The simulations are based on a vertically unstable NSTX equilibrium enclosed by an axisymmetric resistive wall with a rectangular cross section. A linear dependence of the linear VDE growth rates on the resistivity of the wall is recovered for sufficiently large wall conductivity and small temperatures in the open field line region. We report the benchmark results show good agreement between the VDE growth rates obtained from linear NIMROD and M3D-C 1 simulations and from the linear phase of axisymmetric nonlinear JOREK, NIMROD, and M3D-C 1 simulations. Axisymmetric nonlinear simulations of a full VDE performed with the three codes are compared, and an excellent agreement is found regarding the plasma location and plasma currents, as well as eddy and halo currents in the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Exploration of Rayleigh-Taylor seeding mechanisms in laser- and pulsed-power-driven inertial fusion

The ubiquitous Rayleigh-Taylor instability (RTI) has been observed in fluids and plasmas in a wide range of parameter regimes from astrophysical to laboratory plasmas. In high-energy-density (HED) laboratory plasmas, such as laser-driven inertial confinement fusion (ICF) plasmas, the RTI can have a detrimental impact on achieving fusion ignition through the generation of hydrodynamic mix. Hence, an understanding of the seeding mechanisms that produce the RTI and identifying ways to mitigate hydrodynamic mix is of critical importance. Appropriately aligned magnetic fields have been known to stabilize short-wavelength RTI. A number of numerical and experimental studies have demonstrated the benefits of using imposed magnetic fields in laser-driven ICF to achieve higher ion temeratures, higher neutron yields, and a relaxation of the ignition criteria. This work addresses seeding mechanisms in laser-driven implosions to understand critical early-stage physics that ultimately leads to substantial growth of the RTI along with mechanisms for mitigation of this growth. Surface perturbations due to machining tolerances and single-feature seeds (for example due to fill tubes) can produce substantial RTI growth in the ignition-relevant high-convergence ratio targets for laser- driven ICF implosions. RTI growth from single-feature seeding in laser-based implosions has the potential to be mitigated through appropriately aligned externally applied magnetic fields and the goal of this work has been to quantify that numerically and experimentally. Resistive-magnetohydrodynamic (MHD) simu- lations are used to study the seeding and evolution of the RTI leveraging previous support from the DOE HEDLP program. Furthermore, experimental data has been obtained and applied towards code validation of unmagnetized and magnetized evolution of single-feature seeded RTI growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗