Search NASA⌕ Search

SEARCH · Search NASA

Results for “Instability waves”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Momentum transfer and plasma rotation caused by destabilized eigenmodes in tokamaks

In this work, the influence of magnetohydrodynamic eigenmodes destabilized by energetic ions on the momentum of these ions and concomitant sheared plasma rotation are studied. Two mechanisms affecting rotation are revealed: (i) spatial channelling (SC) – radially separated emission and absorption of the momentum; (ii) mode induced redistribution (MIR) across the magnetic field of the momentum of energetic ions. Forces arising during SC and MIR produce both toroidal and poloidal rotations. In addition, the momentum emission during SC leads to a radial flux of fast ions and generation of a radial electric field. Using the developed theory, estimates were made for the ITER (International Thermonuclear Experimental Reactor) 15 MA baseline scenario. They show that a global toroidicity-induced Alfvén eigenmode destabilized by alpha particles and neutral beam injection can result in significant radial electric field and forces applied to plasma. However, available data are not sufficient for a reliable prediction of the effects of SC and MIR in ITER. In general, one can expect that sheared rotation arising after destabilization of Alfvén modes and fast magnetoacoustic modes by energetic ions will tend to suppress the turbulence and improve plasma performance. The importance of plasma rotation is supported, in particular, by the fact that during the JET DTE1 experimental campaign the best parameters were achieved in a deuterium–tritium discharge where the rotation frequency was largest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Radiative and continuum dampings of reversed shear Alfvén eigenmodes and perturbative analysis limitations for tokamaks

A careful theoretical analysis of the excitation of Alfvén eigenmodes (AEs), such as TAE (toroidicity-induced AE) and RSAE (reversed shear AE), by superalfvenic energetic particles is required for reliable predictions of energetic ion relaxation in present day fusion experiments. This includes the evaluation of different AE damping mechanisms including radiative and continuum dampings which are the focus of this study. A recent comprehensive benchmark of different eigenmode solvers including gyrokinetic, gyrofluid and hybrid magenetohydrodynamics (MHD) has shown that employed models may have deficiencies when addressing some of them (Taimourzadeh et al., Nucl. Fusion, vol. 59, 2019, 066006). Here, in this paper, we are studying the radiative and continuum dampings of RSAEs in details which were missing in hybrid NOVA/NOVA-C calculations to prepare a NOVA-C package with a substantial upgrade. Both dampings require the finite Larmor radius (FLR) corrections to AE mode structures to be accounted for. Accurately calculating different damping rates and understanding their parametric dependencies, we resolve the limitation coming out of the perturbative approach. In particular, here, the radiative damping is included perturbatively, whereas the continuum damping is computed non-perturbatively. Our comparison leads to the conclusion that the non-perturbative treatment of the unstable RSAE modes is needed to find the agreement with the gyrokinetic calculations. We expect that the RSAE mode structure modification plays a dominant role in determining the RSAE stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Evaluation of Implicit‐Explicit Additive Runge‐Kutta Integrators for the HOMME‐NH Dynamical Core

Abstract The nonhydrostatic High‐Order Method Modeling Environment (HOMME‐NH) atmospheric dynamical core supports acoustic waves that propagate significantly faster than the advective wind speed, thus greatly limiting the time step size that can be used with standard explicit time integration methods. Resolving acoustic waves is unnecessary for accurate climate and weather prediction. This numerical stiffness is addressed herein by considering implicit‐explicit additive Runge‐Kutta (ARK IMEX) methods that can treat the acoustic waves in a stable manner without requiring implicit treatment of nonstiff modes. Various ARK IMEX methods are evaluated for their efficiency in producing accurate solutions, ability to take large time step sizes, and sensitivity to grid cell length ratio. Both the gravity wave test and baroclinic instability test from the 2012 Dynamical Core Model Intercomparison Project are used to recommend 5 of the 27 ARK IMEX methods tested for use in HOMME‐NH.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Precursor of pair-density wave in doping Kitaev spin liquid on the honeycomb lattice

Abstract We study the effects of doping the Kitaev model on the honeycomb lattice where the spins interact via the bond-directional interaction J K , which is known to have a quantum spin liquid as its exact ground state. The effect of hole doping is studied within the t - J K model on a three-leg cylinder using density-matrix renormalization group. Upon light doping, we find that the ground state of the system has a dominant quasi-long-range charge-density-wave correlations but short-range single-particle correlations. In the pairing channel, the even-parity superconducting correlation is dominant with d -wave-like symmetry, which oscillates in sign as a function of separation with a period equal to that of the spin-density wave and two times the charge-density wave. Although these correlations fall rapidly (possibly exponentially) at long distances, this is never-the-less the example where a pair-density wave is the leading instability in the pairing channel on the honeycomb lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Hot-electron generation at direct-drive ignition-relevant plasma conditions at the National Ignition Facility

Laser–plasma interaction instabilities can be detrimental for direct-drive inertial confinement fusion by generating high-energy electrons that preheat the target. Moreover, an experimental platform has been developed and fielded on the National Ignition Facility to investigate hot-electron production from laser–plasma instabilities at direct-drive ignition-relevant conditions. The radiation-hydrodynamic code DRACO has been used to design planar-target experiments that generate plasma and interaction conditions comparable to direct-drive ignition designs: I L ~10 15 W/cm 2 , T e > 3 keV, and density-gradient scale lengths of L n ~ 600 μ m in the quarter-critical density region. The hot-electron properties were inferred by comparing the experimentally observed hard x-ray spectra to Monte Carlo simulations of hard x-ray emission from hot electrons depositing energy in the target. Hot-electron temperatures of ~40 keV to 60 keV and the fraction of laser energy converted to hot electrons of ~0.5% to 5% were inferred in plastic targets for laser intensities at the quarter-critical density surface of (~4 to 14) × 10 14 W/cm 2 . The use of silicon ablators was found to mitigate the hot-electron preheat by increasing the threshold laser intensity for hot-electron generation from ~3.5 × 10 14 W/cm 2 in plastic to ~6 × 10 14 W/cm 2 in silicon. The overall hot-electron production is further reduced in silicon ablators when the intensity threshold is exceeded.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

On the local measurement of electric currents and magnetic fields using Thomson scattering in Weibel-unstable plasmas

Here, we demonstrate the capability of the Thomson Scattering (TS) diagnostic to measure locally the microscopic electron and ion currents in counter-streaming plasmas unstable to the Weibel or current-filamentation instability. Synthetic TS spectra are calculated with particle distribution functions obtained from particle-in-cell simulations and used to accurately reproduce the simulated currents. We show that this technique allows accurate local measurements of the magnetic field, thus opening the way for the complete experimental characterization of the growth rate, saturation, and nonlinear dynamics of electromagnetic instabilities in plasmas. We illustrate the application of this diagnostic to experimental TS data, which yields local measurements of the magnetic field in Weibel-unstable plasmas and indicates that the magnetic energy density reaches ~1% of the kinetic energy density of the flows, in agreement with previous numerical studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An improved ten-moment closure for reconnection and instabilities

The integration of kinetic effects in fluid models is important for global simulations of Earth's magnetosphere. The use of the two-fluid ten-moment model, which includes the pressure tensor for both species, has had some success in simulating Ganymede and Mercury with a simple closure model. Here, we discuss a heat flux closure which accounts for some limitations of the earlier work while remaining computationally tractable. Comparisons with kinetic simulations for magnetic reconnection and lower-hybrid drift instabilities show good agreement with kinetic results and improvements on previous closure models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The effect of laser entrance hole foil thickness on MagLIF-relevant laser preheat

The magnetized liner inertial fusion (MagLIF) scheme relies on coupling laser energy into an underdense fuel raising the fuel adiabat at the start of the implosion. To deposit energy into the fuel, the laser must first penetrate a laser entrance hole (LEH) foil which can be a significant energy sink and introduce mix. In this paper, we report on experiments investigating laser energy coupling into MagLIF-relevant gas cell targets with LEH foil thicknesses varying from 0.5 μm to 3 μm. Two-dimensional (2D) axisymmetric simulations match the experimental results well for 0.5 μm and 1 μm thick LEH foils but exhibit whole-beam self-focusing and excessive penetration of the laser into the gas for 2 μm and 3 μm thick LEH foils. Additionally, better agreement for the 2 μm-thick foil is found when using a different thermal conductivity model in 2D simulations, while only 3D Cartesian simulations come close to matching the 3 μm-thick foil experiments. The study suggests that simulations may over-predict the tendency for the laser to self-focus during MagLIF preheat when thicker LEH foils are used. This effect is pronounced with 2D simulations where the azimuthally symmetric density channel effectively self-focuses the rays that are forced to traverse the center of the plasma. The extra degree of freedom in 3D simulations significantly reduces this effect. The experiments and simulations also suggest that, in this study, the amount of energy coupled into the gas is highly correlated with the laser propagation length regardless of the LEH foil thickness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Suppression of power losses during laser pulse propagation in underdense plasma slab

For current state-of-the-art terawatt lasers, the primary laser scattering mechanisms in plasma include Forward Raman Scattering (FRS), excitation of plasma waves, and the self-modulational instability (SMI). Using 2D PIC simulations, we demonstrate the dominance of the FRS in the regime with medium-to-low density plasma and non-relativistic laser fields. However, the use of multi-colored lasers with frequency detuning exceeding the plasma frequency, Δω > ω pe , suppresses the FRS. Finally, the laser power can then be transmitted efficiently.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Thresholds of absolute two-plasmon-decay and stimulated Raman scattering instabilities driven by multiple broadband lasers

Here, thresholds for the absolute stimulated Raman scattering (SRS) and two-plasma decay (TPD) instabilities driven by multiple broadband laser beams are evaluated using 3-D simulations at conditions relevant to inertial confinement fusion experiments. Multibeam TPD and SRS backscatter are found to be easier to mitigate with bandwidth than the corresponding single-beam instabilities. A relative bandwidth of 1% increases the threshold for absolute SRS backscatter by a factor of 4 at conditions relevant to ongoing National Ignition Facility experiments and should be sufficient to keep all of the absolute instabilities below threshold in experiments with similar conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Upgrades to the ion cyclotron emission diagnostic on the DIII-D tokamak

The ion cyclotron emission diagnostic on the DIII-D tokamak comprises seven single-turn loops that measure high-frequency (1–100 MHz) magnetic field fluctuations that are often excited by energetic particles in the plasma. The raw voltage signals induced in the loops in response to these fluctuations travel through a series of cables, isolation transformer DC blocks, low-pass filters, and finally a digitizer before being analyzed in frequency space. The diagnostic has been recently upgraded, most notably to include four additional graphite tile loops and a new eight-channel digitizer. The previous three loops are all on the low-field side of the tokamak. The measurement capabilities of the system have been expanded by the addition of a new horizontally oriented loop on the low-field side, an additional toroidal loop on the low-field side, and two toroidal loops on the high-field side. These loops will be used to provide approximate mode polarization, improved toroidal mode number calculations, and information on modes in inward-shifted plasmas, respectively.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Dissipation in the magnetic turbulence of reversed field pinch plasmas

Reversed field pinch plasmas are subject to tearing instability that creates a broad spectrum of magnetic fluctuations. The dominant fluctuations have poloidal and toroidal mode numbers (m,n) = (1,6 - 10) and can grow to 2 - 3% of the mean magnetic field. Through non-linear coupling this growth culminates in a strong reconnection event and broadening of the magnetic spectrum extending to the ion gyroradius scale. Multiple developments occur during the reconnection stage: ions and electrons are energized, magnetic fluctuation amplitudes increase, plasma flow is halted, and the toroidal magnetic flux increases in a sawtooth-like fashion as the RFP dynamo becomes stronger. Magnetic fluctuations are measured in the plasma edge at multiple radial locations from r/a = 0.75 to 0.96 to assess and characterize the magnetic turbulence. The measured spectrum perpendicular to the mean field, S(k ⊥ ), can be fit to a model spectrum consisting of power law and exponential component with one free parameter that characterizes dissipation. The measured dissipation is much larger than estimated from classical viscous or resistive dissipation, but it is consistent with a flow damping measurement of anomalous viscosity. The measurements show an evolution of the spectrum during which fluctuation power builds up in the smallest wavenumbers and cascades to the larger wavenumber due to the non-linear coupling between the linear (m,n) = (1, > 6) and the non-linear (m,n) = (0, 1) tearing mo

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Development of five-moment two-fluid modeling for Z-pinch physics

The Z-pinch m = 0 instability as well as its stabilization by radially sheared axial flow is studied using the nonlinear ideal five-moment twofluid (5M2F) model with an extension of that model to include Braginskii heat and momentum transport. Using the ideal 5M2F model, linear growth rate results are compared with prior work using MHD and Hall MHD. At small normalized wavenumber, 1 < k a < 4 , where a is the effective pinch radius, 5M2F results agree with Hall MHD within 20% in scenarios without radially sheared axial flow. With the sheared flow and focusing on ka = 10/3, agreement with Hall MHD is excellent. In the limit of small ion inertial length, results also match with MHD. A comparison with PIC modeling of shear-free m = 0 stability focuses on a plasma scenario based on recent experimental results. In a scan of mode wavenumber, ideal 5M2F results are qualitatively similar to PIC: the growth rate rises to a peak at a moderate wavenumber and declines at a large wavenumber in contrast to MHD results, which show the saturation of the growth rate with the increasing wavenumber rather than a decline. The peak normalized 5M2F growth rate is &#x3B3; &#x3C4; A = 1.5 , where sA is the Alfven transit time across the pinch. The peak occurs at normalized wavenumber ka = 10. For comparison, PIC results have a peak growth of &#x3B3; &#x3C4; A = 0.77 at ka = 5. Including Braginskiibased closure of the 5M2F model does not qualitatively change the ideal results in this particular case. Nonlinear saturation is studied using the 5M2F model with the dissipative Braginskii-based closure in cases with pinch-edge sheared-flow speed equal to half the Alfven speed. Nonlinear mixing due to the sheared flow yields a quasi-steady state after modest losses of pinch ion inventory and pinch thermal energy, approximately 30% and 10%, respectively. 5M2F modeling captures the essential physics of m = 0 instability and offers a computationally tractable route to high-fidelity modeling of 3D Z-pinch behavior, including m = 1 instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Understanding the negative triangularity ELM trigger and ELM free state on DIII-D with ECE-imaging

Here, the Electron Cyclotron Emission Imaging (ECEI) diagnostic was used to observe a finite-n interchange mode structure in the edge of negative triangularity shaped plasmas on DIII-D. At a small negative triangularity (δ u = -0.2), the plasma is in the H-mode with ELMs that are triggered by a low-n interchange mode. At a larger negative triangularity (δ u = -0.4) and low NBI power (2 MW), a dithering oscillation is observed that is triggered by a low-n interchange mode, whereas at higher NBI power (>2 MW), the edge reverts to L-mode and the low-n interchange mode is present continuously. In all cases, the edge pressure gradient is clamped when the interchange mode is present. It is concluded that the low-n interchange mode prevents the plasma from transitioning to H-mode at a large negative triangularity. This agrees with linear BOUT++ simulations which suggest that the interchange-type MHD can be a resistive ballooning mode whereby resistivity can significantly increase the finite-n ballooning mode growth rate. The absence of H-mode at large negative triangularity can, thus, be explained by the excitation of low-n pressure driven resistive ballooning modes in the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Close coupling of the peeling ballooning and external kink modes in high poloidal-beta discharges with strong shape and large pedestal width

In plasmas with strong shape and large pedestal width, like DIII-D high poloidal-beta discharges, the eigenvalue spectrum of the peeling ballooning mode (PBM) is shown to shift to the very low-n regime near one, where n is the toroidal mode number. This peeling-type eigenvalue spectrum has a further shift to n = 1, thus being smoothly connected to the n = 1 external kink mode (EKM), as the normalized beta (β N ) increases. Once this connection occurs, the mode takes a mixture form of the PBM and EKM with its mode structure varying from the PBM-like to the EKM-like one as β N increases. The mode stability also becomes sensitive to both the local pedestal gradient and global β N , thus allowing an anti-correlation between the two driving forces. Further, these results appear to provide a qualitative explanation of the two unexpected features observed in the DIII-D high poloidal-beta discharges, that is, the dominance of the n = 1 mode in the edge-localized-modes and the negative correlation between the pedestal height and the internal transport barrier strength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ∼1/Np with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly non-equilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Radiation and heat transport in divergent shock–bubble interactions

Shock–bubble interactions (SBIs) are important across a wide range of physical systems. In inertial confinement fusion, interactions between laser-driven shocks and micro-voids in both ablators and foam targets generate instabilities that are a major obstacle in achieving ignition. Experiments imaging the collapse of such voids at high energy densities (HED) are constrained by spatial and temporal resolution, making simulations a vital tool in understanding these systems. In this study, we benchmark several radiation and thermal transport models in the xRAGE hydrodynamic code against experimental images of a collapsing mesoscale void during the passage of a 300 GPa shock. We also quantitatively examine the role of transport physics in the evolution of the SBI. This allows us to understand the dynamics of the interaction at timescales shorter than experimental imaging framerates. We find that all radiation models examined reproduce empirical shock velocities within experimental error. Radiation transport is found to reduce shock pressures by providing an additional energy pathway in the ablation region, but this effect is small (∼1% of total shock pressure). Employing a flux-limited Spitzer model for heat conduction, we find that flux limiters between 0.03 and 0.10 produce agreement with experimental velocities, suggesting that the system is well-within the Spitzer regime. Higher heat conduction is found to lower temperatures in the ablated plasma and to prevent secondary shocks at the ablation front, resulting in weaker primary shocks. Finally, we confirm that the SBI-driven instabilities observed in the HED regime are baroclinically driven, as in the low energy case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β⁠. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗