Search NASA⌕ Search

SEARCH · Search NASA

Results for “microturbulence”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Impact of impurities on peeling–ballooning modes and turbulence in tokamak plasmas

This study investigates the impact of various impurity species on peeling–ballooning (PB) modes and microturbulence in tokamak plasmas through the extension of traditional two-fluid and gyro-landau-fluid (GLF) models. By incorporating finite Larmor radius (FLR) effects, the analysis provides a comprehensive understanding of impurity-driven impact and its interaction with plasma turbulence. Depending on charge state and local plasma conditions, heavy impurities may exhibit gyro-radii larger than those of main ions, which are captured in the extended GLF model presented. Following the presentation of modified two-fluid equations incorporating impurity effects, we systematically analyze the distinctions between impurity and main ion dynamics and their resultant feedback mechanisms on plasma behavior. Derivation of the linear dispersion relation enables quantification of impurity-mediated modifications to: plasma vorticity, diamagnetic drift and gyroviscous effects, electron Hall physics, and FLR effects. BOUT++ – based linear simulations corroborate this formalism, demonstrating systematic stabilization of PB modes upon impurity seeding. And then operational implications for practical impurity control strategies in tokamak devices are proposed. The results underscore the necessity of impurity management to maintain stability and optimize plasma confinement, with specific focus on how FLR effects contribute to transport dynamics. This work paves the way for enhanced modeling and simulation efforts, supporting the development of strategies to control impurity-induced turbulence and improve overall reactor performance.

BOUT++ simulation↗
DOE OSTI2021

An experimental characterization of core turbulence regimes in Wendelstein 7-X

First results from the optimized helias Wendelstein 7-X stellarator (W7-X) have shown that core transport is no longer mostly neoclassical, as is the case in previous kinds of stellarators. Instead, power balance analysis has shown that turbulent transport poses a serious limitation to the global performance of the machine. Furthermore, several studies have found this particularly relevant for ion transport, with core ion temperatures becoming clamped at relatively low values of T i ≃ 1.7 keV, except in the few scenarios in which turbulence can be suppressed. In order to understand the precise turbulent mechanisms at play and thus design improved performance scenarios, it is important to have a clear understanding of the parametric dependencies of turbulent fluctuations, and the relation between them and turbulent transport. As a first step in this direction, in this work we use Doppler reflectometry measurements carried out during a number of relevant operational scenarios to provide a systematic characterization of ion-scale ( k ⊥ ρ i ≃ 1) density fluctuations in the core of W7-X. Then, we study the relation between fluctuation amplitude and plasma profiles and show how distinct regimes can be defined for the former, depending on normalized gradients and . Furthermore, we discuss the importance of other potentially relevant parameters such as T e / T i , E r or collisionality. Comparing the different regimes, we find that turbulence amplitude depends generally on the gradient ratio η i = , as would be expected for ITG modes, with the exception of a range of discharges, for which turbulence suppression may be better explained by an ITG to TEM transition triggered by a drop in collisionality. Finally, we show a number of scenarios under which T i,core > 1.7 keV is achieved and how core fluctuations are suppressed in all of them, thus providing experimental evidence of microturbulence being the main responsible for the limited ion confinement in W7-X.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effects of negative triangularity shaping on energetic particle driven Alfvén eigenmodes in DIII-D

Shape variations from positive to negative triangularity may lead to improved performance regimes in a tokamak with a reduction in microturbulence as well as improved thermal confinement. The motivation of this investigation is to model and compare the neutral beam driven Alfvén eigenmode (AE) activity in two DIII-D discharges with positive and negative triangularity shaping of plasma. The simulations are performed using the linear version of the FAR3d code, which solves the reduced MHD equations for thermal plasma with addition of moment equations for the energetic ion density and parallel velocity with appropriate Landau closure relations. Overall, our numerical results indicate that for similar physical parameters, the unstable AEs observed in the negative triangularity case have lower growth rates as compared to the positive triangularity regime. Our findings may be useful to analyze the influence of the reverse-D like geometry on the AE instabilities in DIII-D and may lead to better configurations for minimizing fast ion losses in a tokamak device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Density wavenumber spectrum measurements, synthetic diagnostic development, and tests of quasilinear turbulence modeling in the core of electron-heated DIII-D H-mode plasmas

Abstract Measurements of the turbulent density wavenumber spectrum, δ n ˆ e ( k ⊥ ) , using the Doppler Back-Scattering (DBS) diagnostic are reported from DIII-D H-mode plasmas with electron cyclotron heating as the only auxiliary heating method. These electron-heated plasmas have low collisionality, ν e ∗ < 1 , T e / T i > 1 , and zero injected torque—a regime expected to be relevant for future fusion devices. We probe density fluctuations in the core ( ρ ≈ 0.7) over a broad wavenumber range, 0.5 ⩽ k ⊥ ⩽ 16 cm −1 ( 0.1 ⩽ k ⊥ ρ s ⩽ 5 ), to characterize plasma instabilities and compare with theoretical predictions. We present a novel synthetic DBS diagnostic to relate the back-scattered power spectrum, P s ( k ⊥ ) —which is directly measured by DBS—to the underlying electron density fluctuation spectrum, δ n ˆ e ( k ⊥ ) . The synthetic DBS P s ( k ⊥ ) spectrum is calculated by combining the SCOTTY beam-tracing code with a model δ n ˆ e ( k ⊥ ) predicted either analytically or numerically. In this work we use the quasi-linear code Trapped Gyro-Landau Fluid (TGLF) to approximate the δ n ˆ e ( k ⊥ ) spectrum. We find that TGLF, using the experimental profiles, is capable of closely reproducing the DBS measurements. Both the DBS measurements and the TGLF-DBS synthetic diagnostic show a wavenumber spectrum with variable decay. The measurements show weak decay ( k −0.6 ) for k < 3.5 cm −1 , with k −2.6 at intermediate- k ( 3.5 ⩽ k ⩽ 8.5 cm −1 ), and rapid decay ( k −9.4 ) for k > 8.5 cm −1 . Scans of physics parameters using TGLF suggest that the normalized ∇ T e scale-length, R / L T e , is an important factor for distinguishing microturbulence regimes in these plasmas. A combination of DBS observations and TGLF simulations indicate that fluctuations remain peaked at ITG-scales (low k ) while R / L T e -driven TEM/ETG-type modes (intermediate/high k ) are marginally sub-dominant.

synthetic diagnostics↗
DOE OSTI2024

Global linear and nonlinear gyrokinetic simulations of tearing modes

Abstract To better understand multi-scale interactions between global tearing modes and microturbulence in the Madison Symmetric Torus (MST) reversed-field pinch (RFP), the global gyrokinetic code Gene is modified to describe global tearing mode instability via a shifted Maxwellian (SM) distribution consistent with experimental equilibria. The implementation of the SM is tested and benchmarked by comparisons with different codes and models. Good agreement is obtained in code-code and code-theory comparisons. Linear stability of tearing modes of a non-reversed MST discharge is studied. A collisionality scan is performed to the lowest order unstable modes ( n = 5, n = 6) and shown to behave consistently with theoretical scaling. The nonlinear evolution is simulated, and saturation is found to arise from mode coupling and transfer of energy from the most unstable tearing mode to small-scale stable modes mediated by the m = 2 tearing mode. The work described herein lays the foundation for nonlinear simulation and analysis of the interaction of tearing modes and gyroradius-scale instabilities in RFP plasmas.

Physics↗
DOE OSTI2024

Fast ion relaxation in ITER mediated by Alfvén instabilities

Abstract We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The role of ion-scale micro-turbulence in pedestal width of the DIII-D wide-pedestal QH mode

The low-edge rotation, intrinsically ELM-free, and improved confinement wide-pedestal quiescent H-mode (QH-mode), discovered in DIII-D tokamak, has pedestal widths exceeding the EPED-kinetic-ballooning mode (KBM) model scaling typically by at least 25%. Ion-scale ($k_yp_s$ < 1) microturbulence and its role in setting the pedestal structure is investigated using the radially local δ$f$ gyrokinetic code CGYRO. The electromagnetic trapped electron mode (TEM) is unstable at the pedestal top, while plasma beta (β$_e$) is ∼60% below the KBM onset threshold and the electron temperature gradient mode is found to be unstable in the peak gradient region. Nonlinear simulation reveals that the ion-scale turbulence could produce electron energy flux consistent with the flux inferred from power balance at the pedestal top, with a reasonable variation of the local shearing rate; and the local neoclassical transport from NEO is dominant over the simulated turbulent transport in the ion energy flux channel. The simulated ion-scale turbulence produces much lower electron energy flux than inferred from experiment in the pedestal peak gradient region. A correction to the EPED-KBM pedestal width scaling is obtained based on the two-dimensional scan of pedestal top plasma beta (β$_e$) and normalized electron density and temperature scale lengths,${a}$/${L_n}_e$, ${a}$/${L_T}_n$using CGYRO linear simulations. Mode transitions among TEM, micro-tearing mode, ion-temperature gradient mode and KBM, are observed in the 2D scan at the pedestal top. A fixed normalized growth rate for these drift-type modes is taken to determine the pedestal width scaling, which shows good consistency with the QH experimental database on pedestal heights and widths. The onset of KBM instabilities and the local E × B shear suppression criterion set the lower and upper limit for the pedestal width of standard QH-mode, wide-pedestal QH-mode and type-I ELMy H mode. A potentially higher and wider pedestal is expected from the new scaling of pedestal width. This work presents an improved understanding of the ion-scale micro-turbulence of wide-pedestal QH-mode and sheds light on a promising scenario for future reactors, including ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Perturbative model for the saturation of energetic-particle-driven modes limited by self-generated zonal modes

We present a simplified energy-conserving approach to incorporate wave–wave nonlinear effects within the framework commonly used to describe wave–particle nonlinearities. In particular, the effects of zonal mode (ZM) generation on the determination of the saturation amplitude of energetic particle (EP)-driven Alfvénic instabilities is studied. The model assumes that the zonal perturbations grow at a rate twice that of the original (pump) wave, consistent with a beat-driven (or force-driven) generation mechanism. The evolution and saturation of the mode amplitude are investigated both analytically and numerically within our reduced model assumptions, in both the collisionless and scattering-dominated regimes. These studies underscore the crucial role of sources and sinks in capturing the impact and the role of beat-driven zonal perturbations on mode evolution. In the realistic case of saturation set by sources and sinks, we discuss the role of a finite amplitude ZM in reducing microturbulent particle scattering, thus limiting the energy source for the resonant mode. We then discuss comparisons between the model’s predictions and simulation results. The model reproduces key features observed in gyrokinetic simulations as the reduction in saturated mode amplitude and the onset of wave–wave nonlinear effects as functions of mode growth rate and amplitude. Thanks to its simplicity, it can be readily implemented into codes based on reduced models, thereby improving their predictive capability for strongly driven instabilities.

Alfvén eigenmodes↗
DOE OSTI2026

First observation of currents induced by Alfven eigenmodes in a magnetic confinement device

Instabilities driven by energetic particles normally reduce fusion power but recent studies find Alfven eigenmode instabilities can drive ``zonal" flows and currents that reduce thermal transport, improving overall performance. Motional Stark effect (MSE) data from the DIII-D tokamak reveal changes in magnetic pitch angle during the phase when Alfven-eigenmode nonlinearity drives zonal flows, suppresses microturbulence, and increases the electron and ion temperatures. Finally, measured changes in safety factor of ∆q/q ≃ 5% occur in roughly 20 ms, consistent with a rough estimate based on the measured zonal flow.

Heidbrink, W. W. [Univ. of California, Irvine, CA ↗
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

Weak collisionless shocks mediated by ion gyroviscosity

Collisionless shocks are ubiquitous in space and astrophysical plasmas, and they are essential dynamical features of these systems. Lacking Coulomb collisions, these shocks are mediated by the anomalous dissipation provided by nonlinear plasma instabilities. By numerically resolving the structure of a steady-state, ion gyroviscous shock, we show that ion gyroviscosity, alone, can produce weak (M≲1.1, where M is the sonic Mach number) shocks in a collisionless, magnetized plasma. We emphasize that this effect does not require an appeal to plasma microturbulence. Moreover, while most collisionless systems may be unsuitable to support purely gyroviscous shocks, we argue that gyro-viscous heating may be an overlooked mechanism, generally; and it may be a key driver within magnetohydrodynamic shocks at large. In conclusion, representative examples include the plasma environments produced on the plasma liner experiment and the magnetized liner inertial fusion platforms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Overview of research completed with funding from the SciDAC CSPM project (Final Technical Report)

The Center for the Study of Plasma Microturbulence was a SciDAC project managed by the Office of Fusion Energy Sciences in the US Department of Energy. Funding for the University of Maryland was initially modest (40k/year) for the first two years. The remaining four years were funded at 201k/year, with a single $50k supplemental in 2015, and a no-cost extension that ran through March, 2018. The CSPM involved several institutions whose efforts are documented separately. At UMD, the CSPM supported two graduate students (Kanekar and Wilkie), two post-doctoral fellows (Landreman and TenBarge), and summer salary for the UMD PI (Dorland). The supplemental in 2015 provided additional salary for a sabbatical by Dorland in 2015/16. Landreman and TenBarge both transitioned to separately funded research staff positions at Maryland at the ends of their post-doctoral, CSPM-supported, positions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Visualization techniques for the gyrokinetic tokamak simulation code

Gyrokinetic simulations of plasma microturbulence in tokamaks are challenging to visualize because the compute grid follows the magnetic field lines that spiral around the torus. We have overcome this challenge by developing three new approaches that improve visualization of gyrokinetics. Our techniques work directly with the topology of magnetic flux surfaces where the simulation stores variables in concentric rings on poloidal planes (vertical cross sections of the torus). Our visualization preview step triangulates each consecutive pair of rings to display the data on a poloidal plane. The second visualization technique follows spiral field lines around the torus and constructs polygons to visualize a flux surface. Third, the poloidal triangles are connected between planes to form prisms that compose a 3-D model of the entire torus. The visualization workflow produces detailed geometry that matches the high resolution, irregular compute grid for every time step. The surface and solid models are displayed in scientific visualization programs to effectively explore and communicate the results, including fluctuation of electron density, ion temperature, and electrostatic potential. Highly detailed renderings verify plasma behavior along magnetic field lines over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

SDSS-IV MaStar: Data-driven Parameter Derivation for the MaStar Stellar Library

The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) Stellar Library (MaStar) is a large collection of high-quality empirical stellar spectra designed to cover all spectral types and ideal for use in the stellar population analysis of galaxies observed in the MaNGA survey. The library contains 59,266 spectra of 24,130 unique stars with spectral resolution R ~ 1800 and covering a wavelength range of 3622–10,354 Å. In this work, we derive five physical parameters for each spectrum in the library: effective temperature (T eff ), surface gravity ($\mathrm{log} g$), metallicity ([Fe/H]), microturbulent velocity ($\mathrm{log}({v}_{\mathrm{micro}})$), and alpha-element abundance ([α/Fe]). These parameters are derived with a flexible data-driven algorithm that uses a neural network model. We train a neural network using the subset of 1675 MaStar targets that have also been observed in the Apache Point Observatory Galactic Evolution Experiment (APOGEE), adopting the independently-derived APOGEE Stellar Parameter and Chemical Abundance Pipeline parameters for this reference set. For the regions of parameter space not well represented by the APOGEE training set (7000 ≤ T ≤ 30,000 K), we supplement with theoretical model spectra. We present our derived parameters along with an analysis of the uncertainties and comparisons to other analyses from the literature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Microphysics of Relativistic Collisionless Electron-ion-positron Shocks

Abstract We perform particle-in-cell simulations to elucidate the microphysics of relativistic weakly magnetized shocks loaded with electron-positron pairs. Various external magnetizations σ ≲ 10 −4 and pair-loading factors Z ± ≲ 10 are studied, where Z ± is the number of loaded electrons and positrons per ion. We find the following: (1) The shock becomes mediated by the ion Larmor gyration in the mean field when σ exceeds a critical value σ L that decreases with Z ± . At σ ≲ σ L the shock is mediated by particle scattering in the self-generated microturbulent fields, the strength and scale of which decrease with Z ± , leading to lower σ L . (2) The energy fraction carried by the post-shock pairs is robustly in the range between 20% and 50% of the upstream ion energy. The mean energy per post-shock electron scales as E ¯ e ∝ Z ± + 1 − 1 . (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low as σ ≈ 5 × 10 −6 . The ions then become essentially thermal with mean energy E ¯ i , while electrons form a nonthermal tail, extending from E ∼ Z ± + 1 − 1 E ¯ i to E ¯ i . When σ = 0, particle acceleration is enhanced by the formation of intense magnetic cavities that populate the precursor during the late stages of shock evolution. Here, the maximum energy of the nonthermal ions and electrons keeps growing over the duration of the simulation. Alongside the simulations, we develop theoretical estimates consistent with the numerical results. Our findings have important implications for models of early gamma-ray burst afterglows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Determining Stellar Elemental Abundances from DESI Spectra with the Data-driven Payne

Abstract Stellar abundances for a large number of stars provide key information for the study of Galactic formation history. Large spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI) and LAMOST take median-to-low-resolution (R≲ 5000) spectra in the full optical wavelength range for millions of stars. However, the line-blending effect in these spectra causes great challenges for elemental abundance determination. Here we employDD-Payne, a data-driven method regularized by differential spectra from stellar physical models, to the DESI early data release spectra for stellar abundance determination. Our implementation delivers 15 labels, including effective temperatureT eff , surface gravity log g , microturbulence velocityv mic , and the abundances for 12 individual elements, namely C, N, O, Mg, Al, Si, Ca, Ti, Cr, Mn, Fe, and Ni. Given a spectral signal-to-noise ratio of 100 per pixel, the internal precisions of the label estimates are about 20 K forT eff , 0.05 dex for log g , and 0.05 dex for most elemental abundances. These results agree with the theoretical limits from the Crámer–Rao bound calculation within a factor of 2. The majority of the accreted halo stars contributed by the Gaia–Enceladus–Sausage are discernible from the disk and in situ halo populations in the resultant [Mg/Fe]–[Fe/H] and [Al/Fe]–[Fe/H] abundance spaces. We also provide distance and orbital parameters for the sample stars, which spread over a distance out to ∼100 kpc. The DESI sample has a significantly higher fraction of distant (or metal-poor) stars than the other existing spectroscopic surveys, making it a powerful data set for studying the Galactic outskirts. The catalog is publicly available.

Astronomy & Astrophysics↗
DOE OSTI2022

Origin of Intense Electron Heating in Relativistic Blast Waves

The modeling of gamma-ray burst afterglow emission bears witness to strong electron heating in the precursor of Weibel-mediated, relativistic collisionless shock waves propagating in unmagnetized electron–ion plasmas. In this Letter, we propose a theoretical model, which describes electron heating via a Joule-like process caused by pitch-angle scattering in the decelerating, self-induced microturbulence and the coherent charge-separation field induced by the difference in inertia between electrons and ions. The emergence of this electric field across the precursor of electron–ion shocks is confirmed by large-scale particle-in-cell (PIC) simulations. Integrating the model using a Monte Carlo-Poisson method, we compare the main observables to the PIC simulations to conclude that the above mechanism can indeed account for the bulk of electron heating.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Effect of radial pressure corrugations and profile shearing on turbulence in Fusion plasmas

Microturbulence can produce stationary fine-scale radial corrugations on the plasma density and temperature gradients in magnetic confinement fusion devices. We show that these structures play a significant role in regulating turbulent transport. We focus on the pedestal, studying electron-temperature-gradient (ETG) mode destabilisation and saturation in the presence of radial corrugations on the electron temperature gradient that could result from microtearing turbulence. A linear dispersion relation is derived for a shearless slab case, which indicates that in the presence of a sinusoidal background corrugation, each ETG mode splits into three distinct eigenvalues, with one being the original, one being more unstable and one being less unstable. However, despite the presence of more unstable linear modes, nonlinear gyrokinetic simulations of ETG with corrugated background electron temperature show a reduction of fluxes. Our investigation reveals a radial variation of the phase velocity of the modes that is proportional to the diamagnetic drift velocity and the local pressure gradient. The associated profile shearing breaks the turbulent eddies apart, reducing the transport level. This profile shearing resulting from fine-scale pressure corrugations could be a ubiquitous turbulence saturation mechanism not just in Fusion plasmas, but in Astrophysics and other areas.

FOS: Physical sciences↗