Pellet ELM triggering in DIII-D with low collisionality, peeling-limited pedestals
Explore the source record for details and available documents.
Wilcox, Robert↗
SEARCH · Search NASA
Results for “collisionality”
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.
PREDICTION OF PELLET MASS THRESHOLDS FOR ELM TRIGGERING IN LOW-COLLISIONALITY, ITER-LIKE DISCHARGES
Explore the source record for details and available documents.
Wingen, Andreas↗
Generalized impurity pinch in partially magnetized multi-ion plasma
In a two-ion-species plasma with disparate ion masses, heavy ions tend to concentrate in the low-temperature region of collisionally magnetized plasma and in the high-temperature region of collisionally unmagnetized plasma, respectively. Moreover, collisional magnetization can be determined as the ratio of the light ion gyrofrequency to the collision frequency of light and heavy ion species, and the behavior of this effect in the intermediate regime of partially magnetized plasma is predominantly dependent on this Hall parameter. Multi-ion cross-field transport has been described before in the collisionally magnetized plasma regime, and generalized pinch relations, which describe densities of ion species in equilibrium in that plasma, are found in the literature. Here, the role of collisional magnetization and Larmor magnetization in multi-ion collisional transport is clarified, and generalized pinch relations are extended to the partially magnetized regime in which the ion Hall parameter may be small, as long as electrons remain collisionally magnetized. Equilibrium ion density profiles have the same dependence on external forces and on each other regardless of collisional magnetization of ions. The expansion of the range of validity of multi-ion collisional transport models makes them applicable to a wider range of laboratory plasma conditions. In particular, ion density profiles evolve sufficiently fast for radial impurity transport to be observable around stagnation on MagLIF, leading to expulsion of heavy ion impurities from the hotspot as long as plasma becomes sufficiently collisionally magnetized during the implosion.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
The Effect of Spatially Varying Collision Frequency on the Development of the Rayleigh–Taylor Instability
Abstract The Rayleigh–Taylor (RT) instability is ubiquitously observed, yet has traditionally been studied using ideal fluid models. Collisionality can vary strongly across the fluid interface, and previous work demonstrates the necessity of kinetic models to completely capture dynamics in certain collisional regimes. Where previous kinetic simulations used spatially and temporally constant collision frequency, this work presents five-dimensional (two spatial, three velocity dimensions) continuum-kinetic simulations of the RT instability using a more realistic spatially varying collision frequency. Three cases of collisional variation are explored for two Atwood numbers: low to intermediate, intermediate to high, and low to high. The low-to-intermediate case exhibits no RT instability growth, while the intermediate-to-high case is similar to a fluid-limit kinetic case with interface widening biased toward the lower-collisionality region. A novel contribution of this work is the low-to-high collisionality case that shows significantly altered instability growth through an upward movement of the interface and damped spike growth due to increased free-streaming particle diffusion in the lower region. Contributions to the energy flux from the non-Maxwellian portions of the distribution function are not accessible to fluid models and are greatest in magnitude in the spike and regions of low collisionality. Increasing the Atwood number results in greater RT instability growth and reduced upward interface movement. Deviation of the distribution function from Maxwellian is inversely proportional to collision frequency and concentrated around the fluid interface. The linear phase of RT instability growth is well described by theoretical linear growth rates accounting for viscosity and diffusion.
79 ASTRONOMY AND ASTROPHYSICS↗
Kinetic Landau-fluid closures of non-Maxwellian distributions
New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
The Implications of Collisions on the Spatial Profile of Electric Potential and the Space-Charge-Limited Current
The space-charge-limited current (SCLC) in a vacuum diode is given by the Child-Langmuir law (CLL), whose electric potential ϕ(x) ∝ (x/D) 4/3 , where x is the spatial coordinate across the gap and D is the gap separation distance. For a collisional diode, SCLC is given by the Mott-Gurney law (MGL) and ϕ(x) ∝ (x/D) 3/2 . Here, we apply a capacitance argument for SCLC and use the transit time from a recent exact solution for collisional SCLC to show that ϕ(x) ∝ (x/D)ξ for a general collisional gap, where 4/3 ≤ ξ ≤ 3/2 . Furthermore, ξ is strictly a function of νT, where ν is the collision frequency and T is the electron transit time. Using this definition of ξ, we estimate the spatial dependence of the electron velocity and use the gap capacitance to derive an analytic equation for collisional SCLC that agrees within ~4.5% of the exact solution that requires solving parametrically through T. This analytic equation for general ξ asymptotically recovers the CLL as ν → 0 and the MGL as ν → ∞. As a result, matching these limits shows that ξ ≈ 1.40 and V ∝ D 2 ν 2 at the transition from a vacuum to a collisional diode for any device condition.
Anodes↗
Anomalous transport in the H-mode pedestal of Alcator C-Mod discharges
Anomalous transport in the H-mode pedestal region of five Alcator C-Mod discharges, representing a collisionality scan is analyzed. The understanding of anomalous transport in the pedestal region is important for the development of a comprehensive model for the H-mode pedestal slope. In this research, a possible role of the drift resistive inertial ballooning modes (Rafiq et al 2010 Phys. Plasmas 17 082511) in the edge of Alcator C-Mod discharges is analyzed. The stability analysis, carried out using the TRANSP code, indicates that the DRIBM modes are strongly unstable in Alcator C-Mod discharges with large electron collisionality. An improved interpretive analysis of H-mode pedestal experimental data is carried out utilizing the additive flux minimization technique (Pankin et al 2013 Phys. Plasmas 20 102501) together with the guiding-center neoclassical kinetic XGC0 code. The neoclassical and neutral physics are simulated in the XGC0 code and the anomalous fluxes are computed using the additive flux minimization technique. The anomalous fluxes are reconstructed and compared with each other for the collisionality scan Alcator C-Mod discharges. It is found that the electron thermal anomalous diffusivities at the pedestal top increase with the electron collisionality. This dependence can also point to the drift resistive inertial ballooning modes as the modes that drive the anomalous transport in the plasma edge of highly collisional discharges.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Reimagining full wave rf quasilinear theory in a tokamak
The velocity dependent resonant interaction of particles with applied radiofrequency (rf) waves during heating and current drive in the presence of pitch angle scattering collisions gives rise to narrow collisional velocity space boundary layers that dramatically enhance the role of collisions as recently shown by Catto ( J. Plasma Phys. , vol. 86, 2020, 815860302). The behaviour is a generalization of the narrow collisional boundary layer that forms during Landau damping as found by Johnston ( Phys. Fluids , vol. 14, 1971, pp. 2719–2726) and Auerbach ( Phys. Fluids , vol. 20, 1977, pp. 1836–1844). For a wave of parallel wave number ${k_{||}}$ interacting with weakly collisional plasma species of collision frequency $\nu$ and thermal speed ${v_{\textrm{th}}}$ , the effective collision frequency becomes of order $\nu {({k_{||}}{v_{th}}/\nu )^{2/3}} \gg \nu $ . The narrow boundary layers that arise because of the diffusive nature of the collisions allow a physically meaningful wave–particle interaction time to be defined that is the inverse of this effective collision frequency. The collisionality implied by the narrow boundary layer results in changes in the standard quasilinear treatment of applied rf fields in tokamaks while remaining consistent with causality. These changes occur because successive poloidal interactions with the rf are correlated in tokamak geometry and because the resonant velocity space dependent interactions are controlled by the spatial and temporal behaviour of the perturbed full wave fields rather than just the spatially local Landau and Doppler shifted cyclotron wave–particle resonance condition associated with unperturbed motion of the particles. The correlation of successive poloidal circuits of the tokamak leads to the appearance in the quasilinear operator of transit averaged resonance conditions localized in velocity space boundary layers that maintain negative definite entropy production.
Physics↗
Finite orbit width effects in large aspect ratio stellarators
New orbit-averaged equations for low collisionality neoclassical fluxes in large aspect ratio stellarators with mirror ratios close to unity are derived. The equations retain finite orbit width effects by employing the second adiabatic invariant $J$ as a velocity-space coordinate and they have been implemented in the orbit-averaged neoclassical code KNOSOS (Velasco et al. , J. Comput. Phys. , vol. 418, 2020, 109512; Velasco et al. , Nucl. Fusion , vol. 61, 2021, 116013). The equations are used to study the $1/\nu$ regime and the lower collisionality regimes. For generic large aspect ratio stellarators with mirror ratios close to unity, as the collision frequency decreases, the $1/\nu$ regime transitions directly into the $\nu$ regime, without passing through a $\sqrt {\nu }$ regime. An explicit formula for the neoclassical fluxes in the $\nu$ regime is obtained. The formula includes the effect of particles that transition between different types of wells. While these transitions produce stochastic scattering independent of the value of the collision frequency in velocity space, the diffusion in real space remains proportional to the collision frequency. The $\sqrt {\nu }$ regime is only recovered in large aspect ratio stellarators close to omnigeneity: large aspect ratio stellarators with large mirror ratios and optimized large aspect ratio stellarators with mirror ratios close to unity. Neoclassical transport in large aspect ratio stellarators with large mirror ratios can be calculated with the orbit-averaged equations derived by Calvo et al. ( Plasma Phys. Control. Fusion , vol. 59, 2017, 055014). In these stellarators, the $\sqrt {\nu }$ regime exists in the collisionality interval $(a/R) \ln (R/a) \ll \nu _{ii} R a/\rho _i v_{ti} \ll R/a$ . In optimized large aspect ratio stellarators with mirror ratios close to unity, the $\sqrt {\nu }$ regime occurs in an interval of collisionality that depends on the deviation from omnigeneity $\delta$ : $\delta ^{2} |\ln \delta | \ll \nu _{ii} R a/\rho _i v_{ti} \ll 1$ . Here, $\nu _{ii}$ is the ion–ion collision frequency, $\rho _i$ and $v_{ti}$ are the ion gyroradius and thermal speed, and $a$ and $R$ are the minor and major radii.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Self-organization in collisionless, high- β turbulence
The magnetohydrodynamic (MHD) equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only possible but common. Recent studies have demonstrated that the kinetic physics inherent to this weakly collisional regime can fundamentally transform the evolution of such plasmas across a wide range of scales. Here, we explore the consequences of pressure anisotropy and Larmor-scale instabilities for collisionless, β $\gg$ 1, turbulence, focusing on the role of a self-organizational effect known as ‘magneto-immutability’. We describe this self-organization analytically through a high-β, reduced ordering of the Chew–Goldberger–Low-MHD (CGL-MHD) equations, finding that it is a robust inertial-range effect that dynamically suppresses magnetic-field-strength fluctuations, anisotropic-pressure stresses and dissipation due to heat fluxes. As a result, the turbulent cascade of Alfvénic fluctuations continues below the putative viscous scale to form a robust, nearly conservative, MHD-like inertial range. These findings are confirmed numerically via Landau-fluid CGL-MHD turbulence simulations that employ a collisional closure to mimic the effects of microinstabilities. We find that microinstabilities occupy a small (~5%) volume-filling fraction of the plasma, even when the pressure anisotropy is driven strongly towards its instability thresholds. We discuss these results in the context of recent predictions for ion-vs-electron heating in low-luminosity accretion flows and observations implying suppressed viscosity in ICM turbulence.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Ab initio quantum scattering calculations for the CO–O 2 system and a new CO–O 2 potential energy surface: O 2 and air broadening of the R(0) line in CO
We present ab initio calculations of the collisional broadening of the R(0) pure rotational line in CO (at 115 GHz) perturbed by O 2 . Our calculations are done in a fully quantum way by solving close-coupling quantum-scattering equations without any approximations. We also report a new, highly accurate CO–O 2 potential energy surface on which we did the quantum-scattering calculations. The calculated collisional broadening agrees with the available experimental data in a wide temperature range. The calculated collisional shift is negligible compared to the broadening, which is also consistent with the experimental data. We combine this result with our previous calculations for the same line in CO perturbed by N 2 [Józwiak et al., J. Chem. Phys. 154, 054314 (2021)]; the obtained air-perturbed broadening of the R(0) pure rotational line in CO and its temperature dependence perfectly agree with the HITRAN database. This result constitutes an important step toward developing a methodology for providing accurate ab initio reference data on spectroscopic collisional line-shape parameters for molecular systems relevant to the Earth’s atmosphere and for populating spectroscopic line-by-line databases.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Axisymmetric gyrokinetic simulation of ASDEX-Upgrade scrape-off layer using a conservative implicit BGK collision operator
Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become overly restrictive in an explicit time integration scheme, leading to high computational cost. With the aim of overcoming such restriction, we have implemented an implicit Bhatnagar–Gross–Krook (BGK) collision operator for use in the discontinuous Galerkin full-f gyrokinetic solver within the Gkeyll framework, which, when combined with Gkeyll's traditional explicit time integrator for collisionless advection, can significantly increase the time step in gyrokinetic simulations of highly collisional regimes. To ensure conservation of density, momentum, and energy, we utilize an iterative scheme to correct the discretized approximation to the equilibrium Maxwellian distribution to which the BGK collision operator relaxes. We have further generalized the BGK infrastructure, both the implicit scheme and the correction routine, to handle cross-species collisions. This improved implicit and conservative BGK operator is benchmarked against the more accurate but more computationally expensive Lenard–Bernstein–Dougherty (LBD) operator, which has been utilized in prior studies with Gkeyll. The implicit BGK operator enables 2D axisymmetric simulations of the ASDEX-Upgrade scrape-off layer to run 56 times faster to completion than the simulations with the LBD operator, because the BGK operator is more robust and converges at a lower resolution than is required by the LBD operator. Additionally, in this more collisional limit, we demonstrate that the results of our simulations utilizing the implicit BGK operator agreed well with simulations utilizing the more computationally expensive LBD operator.
Gyrokinetic simulations↗
Micro-tearing mode dominated electron heat transport in DIII-D H-mode pedestal
Abstract A new, comprehensive set of evidence reveals that Micro-Tearing Modes (MTMs) dominate pedestal electron heat transport in an H-mode experiment in the DIII-D tokamak. The experiment investigates the role of MTMs by scanning pedestal collisionality, a main drive of MTM instability, from 0.43 to 0.84 on the pedestal top. Broadband (150–800 kHz) magnetic and density fluctuations originating from the pedestal gradient region and highly consistent with MTMs are observed, with amplitude increasing during the scan. The higher magnetic fluctuation amplitude correlates with a lower pedestal electron temperature gradient, implying MTMs may regulate the pedestal electron heat transport. The collisionality scan results in profile and transport changes consistent with predicted transport capability of MTMs: (1) experimentally-determined electron heat diffusivity increases ∼40% at the location where the broadband density fluctuations peak; (2) ion heat diffusivity has less increase (<20%); and (3) a locally flattened region in the electron temperature pedestal is observed at high collisionality. A local, linear gyrokinetic simulation finds MTMs as the most unstable mode in the pedestal gradient region. In addition, local, nonlinear simulations suggest MTMs can dominate and drive experimentally-relevant, megawatt-level electron heat flux. This result establishes MTMs as an effective transport mechanism in the H-mode pedestal, in particular at high collisionality.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Linear simulations of a wide pedestal quiescent H-mode plasma with the extended-MHD code NIMROD
We present linear simulations of a Wide Pedestal QH (WPQH)-mode DIII-D plasma with the NIMROD code to address the role of ion multispecies collisionality on WPQH-modes. We use a full extended-MHD model that captures Alfvénic and magnetoacoustic waves, including two-fluid, ion gyroviscosity and cross heat flux. We show that two-fluid and ion gyroviscous effects destabilize electron-directed peeling–ballooning modes in the pedestal. The linear growth rates of the modes decrease when carbon (C) is added to a deuterium plasma because the multispecies collisional effects increase the plasma resistivity. When replacing C with tungsten (W), the multispecies collisional effects increase the resistivity further, as well as the stabilizing effect on the instabilities. The resistivity impact on mode stability is reversed when both two-fluid and ion gyroviscous effects are ignored, which is consistent with previous works, proving that multispecies collisional effects, together with two-fluid and ion gyroviscosity, should be considered to address the stability of plasmas with significant impurity content.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Thermodynamic Phase Transition in Magnetic Reconnection
By examining the entropy production in fully kinetic simulations of collisional plasmas, it is shown that the transition from collisional Sweet-Parker reconnection to collisionless Hall reconnection may be viewed as a thermodynamic phase transition. The phase transition occurs when the reconnection electric field satisfies E = $E_{D}\sqrt{m{e}/{m}_i}$ , where ${m{e}/{m}_i}$ is the electron-to-ion mass ratio and E D is the Dreicer electric field. This condition applies for all ${m{e}/{m}_i}$, including ${m{e}/{m}_i}$ = 1 , where the Hall regime vanishes and a direct phase transition from the collisional to the kinetic regime occurs. In the limit ${m{e}/{m}_i}$ → 0, this condition is equivalent to there being a critical electron temperature $Te ≈ m_iΩ^{2}_{i}δ^{2}$, where $Ω_i$ is the ion cyclotron frequency and δ is the current sheet half-thickness. The heat capacity of the current sheet changes discontinuously across the phase transition, and a critical power law is identified in an effective heat capacity. A model for the time-dependent evolution of an isolated current sheet in the collisional regime is derived.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Gaussian Mixture Model Solvers for the Boltzmann Equation
This report documents our experience constructing a numerical method for the collisional Boltzmann equation that is capable of accurately capturing the collisionless through strongly collisional limits. We explore three different functional representations and present a detailed account of a numerical method based on a spatially dependent Gaussian mixture model (GMM). The Kullback-Leibler divergence is used as a closeness measure and various expectation maximization (EM) solution algorithms are implemented to find a compact representation in velocity space for distribution functions that exhibit significant non-Maxwellian character. We discuss issues that appear with this representation over a range of Knudsen numbers for a prototypical test problem and demonstrate that the strongly collisional limit recovers a solution to Euler's equations. Looking forward, this approach is broadly applicable to the non-relativistic and relativistic collisional Vlasov equations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Reimagining full wave rf quasilinear theory in a tokamak
The velocity dependent resonant interaction of particles with applied radio frequency (rf) waves during heating and current drive in the presence of pitch angle scattering collisions gives rise to narrow collisional velocity space boundary layers that dramatically enhance the role of collisions as recently shown by Catto (J. Plasma Phys., vol. 86, 815860302, 2020). The behavior is a generalization of the narrow collisional boundary layer that forms during Landau damping as found by Johnston (Phys. Fluids, vol. 14, 1971, pp. 2719-2726) and Auerbach (Phys. Fluids, vol. 20, 1977, pp. 1836-1844). For a wave of parallel wave number k|| interacting with weakly collisional plasma species of collision frequency ν and thermal speed vth , the effective collision frequency becomes of order ν(k_||v_th /ν)^2/3>> ν . The narrow boundary layers that arise because of the diffusive nature of the collisions allows a physically meaningful wave-particle interaction time to be defined that is the inverse of this effective collision frequency. The collisionality implied by the narrow boundary layer results in changes in the standard quasilinear treatment of applied rf fields in tokamaks while remaining consistent with causality. These changes occur because successive poloidal interactions with the rf are correlated in tokamak geometry and because the resonant velocity space dependent interactions are controlled by the spatial and temporal behavior of the perturbed full wave fields rather than just the spatially local Landau and Doppler shifted cyclotron wave-particle resonance condition associated with unperturbed motion of the particles. The correlation of successive poloidal circuits of the tokamak leads to the appearance in the quasilinear operator of transit averaged resonance conditions localized in velocity space boundary layers that maintain negative definite entropy production.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Ionization dynamics of intense laser-produced argon plasmas revealed by NLTE modeling
The ionization dynamics and transient behavior of under-dense plasma irradiated by an intense laser are investigated. We report two significant effects in the ionization behavior: (1) a surprisingly large delay in ionization response and (2) a stepwise ionization process which involves collisional and laser-driven photoionization (LDP) processes. Ionization induced by intense lasers can exhibit delayed responses due to rapid changes in conditions, particularly when atomic transition processes occur more slowly than the relevant time scales. Furthermore, modeling reveals that the two-step ionization process—collisional excitation followed by LDP—plays an important role in this ionization delay, with collisional excitation acting as the bottleneck. Even low-energy photons (∼3.5 eV) can predominantly ionize plasmas, challenging the conventional belief that such energies are insufficient to overcome the binding energy of bound electrons. These findings underscore the necessity of including such processes into plasma simulations for various laser-plasma experiments.
Collisional excitation↗