Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Plasma 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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Computational general relativistic force-free electrodynamics: II. Characterization of numerical diffusivity

Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe’s most extreme flows of energy. The discretization of Maxwell’s equations – needed to make highly magnetized (astro)physical plasma amenable to its numerical modeling – introduces numerical diffusion. It acts as a source of dissipation independent of the system’s physical constituents. Understanding the numerical diffusion of scientific codes is the key to classifying their reliability. It gives specific limits in which the results of numerical experiments are physical. We aim at quantifying and characterizing the numerical diffusion properties of our recently developed numerical tool for the simulation of general relativistic force-free electrodynamics by calibrating and comparing it with other strategies found in the literature. Our code correctly models smooth waves of highly magnetized plasma. We evaluate the limits of general relativistic force-free electrodynamics in the context of current sheets and tearing mode instabilities. We identify that the current parallel to the magnetic field (j ∥ ), in combination with the breakdown of general relativistic force-free electrodynamics across current sheets, impairs the physical modeling of resistive instabilities. We find that at least eight numerical cells per characteristic size of interest (e.g., the wavelength in plasma waves or the transverse width of a current sheet) are needed to find consistency between resistivity of numerical and of physical origins. High-order discretization of the force-free current allows us to provide almost ideal orders of convergence for (smooth) plasma wave dynamics. The physical modeling of resistive layers requires suitable current prescriptions or a sub-grid modeling for the evolution of j ∥ .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Experimental evidence of runaway electron tail generation via localized helical structure in pellet-triggered tokamak disruptions

A novel detector, using stacked BGO crystals, is developed for runaway electron (REs) studies in the DIII-D tokamak. It is able to resolve fast dynamics of high-energy tail formation of REs with an ultra-high time resolution of ~1 μs. As a cost, the detector estimates the `effective' energy of a given shape of γ-ray spectra and sacrifices the energy resolution. In aid of the new measurement capability, a rapid, inhomogeneous growth of RE tail is observed in detail during a major disruption triggered by an argon pellet. It is found that both the population and energy of a well-confined RE tail significantly oscillate at the early period of the growth. The oscillation phase is locked to a slow rotating magnetohydrodynamic instability, which is briefly destabilized for only ~1 ms at the early period of the current quench. The oscillation ceases promptly, when the mode disappears. As a result, the data suggests that the high-energy RE tail is well-confined and accelerated via a localized helical structure in the plasma core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Integrated simulations of premagnetized and self-magnetizing dynamic screw pinch-driven MagLIF

Magnetically driven implosions such as in magnetized liner inertial fusion (MagLIF) on the Z accelerator suffer from magneto-Rayleigh–Taylor instabilities (MRTI) that dynamically redistribute liner mass during implosion, limiting fusion fuel compression and confinement, which ultimately degrades performance. Driving the implosion with an initially helical drive field that dynamically shifts the direction of the magnetic field surrounding the liner (i.e., a dynamic screw pinch, DSP) is a method proposed to mitigate MRTI in-flight and improve target performance. In DSPs, the axial drive magnetic field component implodes the liner and diffuses through the shocked, melted liner material into the fuel throughout the implosion. Liners can be designed to enable enough axial magnetic flux to diffuse through the liner material to effectively magnetize the fuel region without the need of an initial axial magnetic field (i.e., from external field coils). We present results from three-dimensional radiation-magnetohydrodynamic simulations of MagLIF implosions employing drive magnetic fields composed of axial and azimuthal components (a helical drive field). These simulated DSP-driven MagLIF targets demonstrate improved fuel conditions and thermonuclear yield compared to a traditional MagLIF target implosion. Synthetic x-ray radiography of the imploding liner material and x-ray emission images of the fuel region at the time of peak neutron yield rate indicate superior implosion morphology for DSP-MagLIF implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Case Study of Nonresonant Mode 3-s ULF Waves Observed by MMS

The nature of the 3-s ultralow frequency (ULF) wave in the Earth's foreshock region and the associated wave-particle interaction are not yet well understood. In this work, we investigate the 3-s ULF waves using Magnetospheric Multiscale (MMS) observations. By combining the plasma rest frame wave properties obtained from multiple methods with the instability analysis based on the velocity distribution in the linear wave stage, the ULF wave is determined to be due to the ion/ion nonresonant mode instability. The interaction between the wave and ions is analyzed using the phase relationship between the transverse wave fields and ion velocities and using the longitudinal momentum equation. During the stage when ULF waves have sinusoidal waveforms up to |dB|/|B 0 | ~ 3, where dB is the wave magnetic field and B0 is the background magnetic field, the wave electric fields perpendicular to B 0 do negative work to solar wind ions; along B 0 , a longitudinal electric field develops, but the V × B force is stronger and leads to solar wind ion deceleration. During the same wave stage, the backstreaming beam ions gain energy from the transverse wave fields and get deceleration along B 0 by the longitudinal electric field. The ULF wave leads to electron heating, preferentially in the direction perpendicular to the local magnetic field. Secondary waves are generated within the ULF waveforms, including whistler waves near half of the electron cyclotron frequency, high-frequency electrostatic waves, and magnetosonic whistler waves. The work improves the understanding of the nature of 3-s ULF waves and the associated wave-particle interaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Instability seeding mechanisms due to internal defects in inertial confinement fusion targets

Performance degradation in laser-driven inertial confinement fusion (ICF) implosions is caused by several effects, one of which is Rayleigh–Taylor instability growth. Defects in ICF targets, such as internal voids and surface roughness, create instability seeds in the shell as shocks propagate through the target. A comprehensive understanding of seeding mechanisms is essential to characterize the impact of target defects on inflight shell integrity and mass injection into the central, lower-density vapor region. An analysis of early-time behavior of both single-mode shell mass modulations and isolated voids is performed by examining the evolution of the acoustic waves launched by these target imperfections. A systematic study of localized perturbation growth as a function of defect placement and size is presented. Furthermore, the use of low-density ablator materials (such as foams) is suggested as a potential mitigation strategy to improve target robustness against the impact of defect-initiated growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Deceleration-stage Rayleigh–Taylor growth in a background magnetic field studied in cylindrical and Cartesian geometries

Experiments have identified the Rayleigh–Taylor (RT) instability as one of the greatest obstacles to achieving inertial confinement fusion. Consequently, mitigation strategies to reduce RT growth and fuel–ablator mixing in the hotspot during the deceleration phase of the implosion are of great interest. In this work, the effect of seed magnetic fields on deceleration-phase RT growth are studied in planar and cylindrical geometries under conditions relevant to the National Ignition Facility (NIF) and Omega experiments. The magnetohydrodynamic (MHD) and resistive-MHD capabilities of the FLASH code are used to model imploding cylinders and planar blast-wave-driven targets. Realistic target and laser parameters are presented that suggest the occurrence of morphological differences in late-time RT evolution in the cylindrical NIF case and a measurable difference in spike height of single-mode growth in the planar NIF case. The results of this study indicate the need for target designs to utilize an RT-unstable foam–foam interface in order to achieve sufficient magnetic field amplification to alter RT evolution. Benchmarked FLASH simulations are used to study these magnetic field effects in both resistive and ideal MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Evolution of highly multimodal Rayleigh–Taylor instabilities

Rayleigh–Taylor (RT) instabilities are important fluid instabilities that arise in inertial confinement fusion (ICF) capsule implosions, and many other contexts. Multi-mode coupling is observed in experiments and plays a substantial role in material mix from RT instabilities. In this work, we study the evolution of highly multimodal perturbations (power law distribution) that approximate those found at manufactured material interfaces. We use simulations of over 2000 different perturbations in the LANL code xRAGE to identify distinct phases in the processes of bubble growth and bubble merger which can be visualized in a 2D phase portrait with clear regimes of mode growth and decay. Our results show that the dynamic evolution of the instability strongly depends on the mode of the perturbations and mode interactions. The merger process accelerates bubble growth. A non-Markovian region and a transition of the instability from: (1) initial exponential growth to (2) linear growth and to (3) quadratic growth and asymptotic behavior, are clearly captured in the phase space. We have developed a quantitative model of bubble growth that reproduces the dynamic behavior of ensembles of perturbations. Implications for ICF capsules designed for robustness against instabilities are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

3D Numerical Simulation of Kink-driven Rayleigh–Taylor Instability Leading to Fast Magnetic Reconnection

Fast magnetic reconnection involving non-MHD microscale physics is thought to underlie both solar eruptions and laboratory plasma current disruptions. While there is extensive research on both the MHD macroscale physics and the non-MHD microscale physics, the process by which large-scale MHD couples to the microscale physics is not well understood. An MHD instability cascade from a kink to a secondary Rayleigh–Taylor instability in the Caltech astrophysical jet laboratory experiment provides new insights into this coupling and motivates a 3D numerical simulation of this transition from large to small scale. A critical finding from the simulation is that the axial magnetic field inside the current-carrying dense plasma must exceed the field outside. In addition, the simulation verifies a theoretical prediction and experimental observation that, depending on the strength of the effective gravity produced by the primary kink instability, the secondary instability can be Rayleigh–Taylor or mini-kink. Finally, it is shown that the kink-driven Rayleigh–Taylor instability generates a localized electric field sufficiently strong to accelerate electrons to very high energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Liner implosion experiments driven by a dynamic screw pinch

This paper expands upon recent experimental results [Campbell et al., Phys. Rev. Lett. 125, 035001 (2020)], where thin-foil liner implosions were driven by a dynamic screw pinch (DSP) and found to have magneto-Rayleigh–Taylor instability (MRTI) amplitudes up to three times smaller than in implosions driven by a standard z-pinch (SZP). Here, the expanded discussion presented herein includes: (1) a detailed comparison of the MRTI growth measured in the experiment with that calculated from theory; (2) measurements of axial magnetic field injection into the liner interior prior to the implosion, as well as the subsequent compression of this field during the implosion; (3) an in-depth description of how the helical geometry of the DSP can result in earlier implosion and stagnation times relative to the SZP; and (4) particle-in-cell simulations showing different electron drift behavior in the anode–cathode gap of the DSP relative to the SZP, and how this difference may be related to the different current waveforms recorded during the experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Visualizing plasmons and ultrafast kinetic instabilities in laser-driven solids using X-ray scattering

Ultra-intense lasers that ionize atoms and accelerate electrons in solids to near the speed of light can lead to kinetic instabilities that alter the laser absorption and subsequent electron transport, isochoric heating, and ion acceleration. These instabilities can be difficult to characterize, but X-ray scattering at keV photon energies allows for their visualization with femtosecond temporal resolution on the few nanometer mesoscale. Here, we perform such experiment on laser-driven flat silicon membranes that shows the development of structure with a dominant scale of 60 nm in the plane of the laser axis and laser polarization, and 95 nm in the vertical direction with a growth rate faster than 0.1 fs -1 . Combining the XFEL experiments with simulations provides a complete picture of the structural evolution of ultra-fast laser-induced plasma density development, indicating the excitation of plasmons and a filamentation instability. Particle-in-cell simulations confirm that these signals are due to an oblique two-stream filamentation instability. These findings provide new insight into ultra-fast instability and heating processes in solids under extreme conditions at the nanometer level with possible implications for laser particle acceleration, inertial confinement fusion, and laboratory astrophysics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Time-step dependent force interpolation scheme for suppressing numerical Cherenkov instability in relativistic particle-in-cell simulations

The WT scheme, a piecewise polynomial force interpolation scheme with time-step dependency, is proposed in this paper for relativistic particle-in-cell (PIC) simulations. The WT scheme removes the lowest order numerical Cherenkov instability (NCI) growth rate for arbitrary time steps allowed by the Courant condition. While NCI from higher order resonances is still present, the numerical tests show that for smaller time steps, the numerical instability grows much slower than using the optimal time step found in previous studies. The WT scheme is efficient for improving the quality and flexibility of relativistic PIC simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A two-dimensional numerical study of the magneto-Rayleigh–Taylor instability with FLASH: Application to the staged Z-pinch concept

Magnetically driven implosions involving a liner collapsing onto a target are inherently vulnerable to the magneto-Rayleigh–Taylor instability (MRTI). Among the various approaches proposed to achieve fusion conditions within the target, the staged Z-pinch (SZP) concept employs a high-Z liner, the advantages of which remain an active area of investigation. Consequently, ongoing design optimization efforts are essential, while critical physical processes such as magnetic field (B-field) diffusion and radiation transport increase the complexity of required simulations. In this study, we utilize the new capabilities of the FLASH code to simulate in 2D, for the first time using this code, staged Z-pinch configurations (designated SZP1 and SZP1*), focusing on their stability. First, a comparison of simulation results with theoretical predictions of MRTI growth provide new insights into the often-overlooked influences of high mesh resolution and initial perturbation seeding on instability dynamics. These findings then lead to a discussion on potential improvements for future SZP experiments. These include the use of an axial B-field stabilization and the optimization of radiation transport processes during the implosion. The results presented herein establish the framework for simulating multi-dimensional Z pinches using the FLASH code and pave the way for the development of innovative experimental configurations leveraging its advanced simulation capabilities.

Adaptive mesh refinement↗
DOE OSTI2025

Hot-spot model for inertial confinement fusion implosions with an applied magnetic field

Imposing a magnetic field on inertial confinement fusion implosions magnetizes the electrons in the compressed fuel; this suppresses thermal losses, which increases temperature and fusion yield. Indirect-drive experiments at the National Ignition Facility with 12 and 26 T applied magnetic fields demonstrate up to 40% increase in temperature, 3× increase in fusion yield, and indicate that magnetization alters the radial temperature profile [Moody et al., Phys. Rev. Lett. 129, 195002 (2022); Lahmann et al., APS DPP (2022)]. In this work, we develop a semi-analytic hot-spot model, which accounts for the two-dimensional (2D) Braginskii anisotropic heat flow due to an applied axial magnetic field. First, we show that hot-spot magnetization alters the radial temperature profile, increasing the central peakedness, which is most pronounced for moderately magnetized implosions (with 8–14 T applied field), compared to both unmagnetized (with no applied field) and highly magnetized (with 26 T or higher applied field) implosions. This model explains the trend in the experimental data, which finds a similarly altered temperature profile in the 12 T experiment. Next, we derive the hot-spot model for gas-filled (Symcap) implosions, accounting for the effects of magnetization on the thermal conduction and in changing the radial temperature (and density) profiles. Using this model, we compute predicted central temperature amplification and yield enhancement scaling with the applied magnetic field. The central temperature fits the experimental data accurately, and the discrepancy in the yield suggests a systematic (independent of applied field) degradation, such as mix, and additional degradation in the reference unmagnetized shot, such as reduced laser drive, increased implosion asymmetry, or the magnetic field suppressing ablator mixing into the hot-spot.

Alpha particles↗
DOE OSTI2021

A comparative study of internal kink stability in EU DEMO designs with negative and positive triangularity

Abstract Internal kink (IK) instability is investigated for European demonstration fusion reactor (EU DEMO) plasmas in both negative triangularity (NT) and positive triangularity (PT) configurations. For NT plasmas, the IK becomes more unstable as an ideal conformal wall moves away from the plasma boundary, with the mode growth rate saturating at the wall radial location of about b / a = 1.5 , where a is the plasma minor radius and b the wall radial location. The plasma resistivity destabilizes the IK mode. The effect of sub-sonic toroidal plasma flow is sufficiently weak and can thus be ignored for these EU DEMO equilibria. These results are consistent with those for PT plasmas, albeit with larger mode growth rate in the NT configuration. Both perturbative and self-consistent magneto-hydrodynamic (MHD)-kinetic hybrid calculations predict (partial) stabilization of the IK modes in both NT and PT configurations, with inclusion of various kinetic contributions. Precessional drift motion of trapped fusion-born alphas in EU DEMO produces weak stabilization to the IK mode. Stronger stabilization occurs with the toroidal precession of trapped thermal particles (ions and electrons) and the bounce-transit motion of thermal ions. The stabilization is similar between the NT and PT configurations, due to the similarity of the mode eigenfunction (occupying a nearly circular region in the plasma core) despite the sign difference in the triangularity. The non-perturbative MHD-kinetic hybrid model predicts much less stabilization of the mode than the perturbative model, primarily due to the self-consistent determination of the mode eigenvalue in the former. Generally, no significant difference in the IK mode stability is found between the NT and PT plasmas in EU DEMO.

Physics↗
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 OSTI2022

Harmonic Generation and Inverse Cascade in the z-Pinch Driven, Preseeded Multimode, Magneto-Rayleigh-Taylor Instability

We report the magneto-Rayleigh-Taylor instability (MRTI) plays an essential role in astrophysical systems and in magneto-inertial fusion, where it is known to be an important degradation mechanism of confinement and target performance. In this Letter, we show for the first time experimental evidence of mode mixing and the onset of an inverse-cascade process resulting from the nonlinear coupling of two discrete preseeded axial modes (400- and 550 - μm wavelengths) on an Al liner that is magnetically imploded using the 20-MA, 100-ns rise-time Z Machine at Sandia National Laboratories. Four radiographs captured the temporal evolution of the MRTI. We introduce a novel unfold technique to analyze the experimental radiographs and compare the results to simulations and to a weakly nonlinear model. We find good quantitative agreement with simulations using the radiation magnetohydrodynamics code hydra. Spectral analysis of the MRTI time evolution obtained from the simulations shows evidence of harmonic generation, mode coupling, and the onset of an inverse-cascade process. The experiments provide a benchmark for future work on the MRTI and motivate the development of new analytical theories to better understand this instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

AC compensation of 3D magnetic diagnostic signals in DIII-D and National Spherical Torus Experiment-Upgrade (NSTX-U) for real-time application

A time domain algorithm has been developed to remove the vacuum pickup generated by both coil current (DC) and induced vessel current (AC) in real time from three dimensional (3D) magnetic diagnostic signals in the National Spherical Torus Experiment-Upgrade (NSTX-U) and DIII-D tokamaks. The possibility of detecting 3D plasma perturbations in real time is essential in modern and future tokamaks to avoid and control MHD instabilities. The presence of vacuum field pickup, due to toroidally asymmetric (3D) coils or to misalignment between sensors and axisymmetric (2D) coils, pollutes the measured plasma 3D field, making the detection of the magnetic field produced by the plasma challenging. Although the DC coupling between coils and sensors can be easily calculated and removed, the AC part is more difficult. Here, an algorithm based on a layered low-pass filter approach for the AC compensation and its application for DIII-D and NSTX-U data is presented, showing that this method reduces the vacuum pickup to the noise level. Comparison of plasma response measurements with and without vacuum compensation shows that accurate mode locking detection and plasma response identification require precise AC and DC compensations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Statistical model of the stimulated forward Brillouin scattering driven by a randomized laser beam in plasma

The modeling of a spatially incoherent laser beam remains a central problem of the parametric instabilities in the context of inertial confinement fusion. This letter gives a simplified and comprehensive overview of the recent analytical developments regarding the modeling of these laser beams and a comparison with a dedicated experiment. Our model accounts for the first time for the statistical standard deviation of the gain and accurately captures the entanglement between wave mixing processes and the speckle correlations thus resolving the longstanding contradictions between the random phase approximation and the model of independent speckles. It is successfully compared to a recent laser beam spray experiment and the associated paraxial simulations, demonstrating that backscattering predictions require accounting for the beam spray. Furthermore, our framework thus provides a way to evaluate and guide the analysis of parametric instabilities in high laser energy experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗