A first-principles self-consistent model of plasma turbulence and kinetic neutral dynamics in the tokamak scrape-off layer

A first-principles self-consistent model that couples plasma and neutral physics suitable for the simulation of turbulent plasma behavior in the tokamak SOL is presented. While the plasma is modeled by the drift-reduced two fluid Braginskii equations, a kinetic model for the neutrals is developed, valid in short and in long mean free path scenarios. The model includes ionization, charge-exchange, recombination, and elastic collisional processes. The solution of the neutral kinetic equation is implemented within the GBS plasma turbulence code (Ricci et al 2012 Plasma Phys. Control. Fusion 54 124047) and it is performed by using the method of characteristics. The details of the numerical implementation are discussed. Finally, we show initial results of the first self-consistent simulations of plasma turbulence and neutral dynamics

Impact of neutral density fluctuations on gas puff imaging diagnostics

A three-dimensional turbulence simulation of the SOL and edge regions of a toroidally limited tokamak is carried out. The simulation couples self-consistently the drift-reduced two-fluid Braginskii equations to a kinetic equation for neutral atoms. A diagnostic neutral gas puff on the low-field side midplane is included and the impact of neutral density fluctuations on D alpha light emission investigated. We find that neutral density fluctuations affect the D alpha emission. In particular, at a radial distance from the gas puff smaller than the neutral mean free path, neutral density fluctuations are anti-correlated with plasma density, electron temperature, and D alpha fluctuations. It follows that the neutral fluctuations reduce the D alpha emission in most of the observed region and, therefore, have to be taken into account when interpreting the amplitude of the D alpha emission. On the other hand, higher order statistical moments (skewness, kurtosis) and turbulence characteristics (such as correlation length, or the autocorrelation time) are not significantly affected by the neutral fluctuations. At distances from the gas puff larger than the neutral mean free path, a non-local shadowing effect influences the neutral density fluctuations. There, the D alpha fluctuations are correlated with the neutral density fluctuations, and the high-order statistical moments and measurements of other turbulence properties are strongly affected by the neutral density fluctuations

Seasonal phenology and starch allocation patterns in populations of Oxycaryum cubense f. cubense and paraguayense in Mississippi, Louisiana, and Florida

Phenological studies for Cuban bulrush [Oxycaryum cubense (Poepp. & Kunth) Lye] have been limited to the monocephalous form in Lake Columbus (Mississippi). Accordingly, there is little available information on potential phenological differences among O. cubense forms (monocephalous vs. polycephalous) and populations in other geographic locations in the United States. Therefore, seasonal patterns of biomass and starch allocation in O. cubense were quantified from two populations in Lake Columbus on the Tennessee-Tombigbee Waterway in Mississippi (monocephalous), two populations from Lake Martin in Louisiana (polycephalous), and two populations from Orange Lake in Florida (polycephalous). Monthly samples of O. cubense inflorescence, emergent, and submersed tissue were harvested from two plots per state from October 2021 to September 2022. During monthly data collection, air temperature and photoperiod were recorded. Starch allocation patterns were similar among all sites, with starch storage being less than 1.5% dry weight for all plant tissues. Biomass was greatest in Lake Columbus (monocephalous; 600.7 g dry weight [DW] m−2) followed by Lake Martin (polycephalous; 392.3 g DW m−2) and Orange Lake (polycephalous; 233.85 g DW m−2). Peak inflorescence biomass occurred in the winter for the Lake Martin and Orange Lake populations and in the summer for the Lake Columbus population. Inflorescence biomass in Lake Columbus had a positive relationship (r2 = 0.53) with warmer air temperatures. Emergent and submersed biomass generally had negative relationships with both photoperiod and temperature (r2 = 0.02 to 0.77) in all sites. Peak biomass was also negatively related to temperature and photoperiod. Results from this study indicate that there are differences in biomass allocation between the two growth forms of O. cubense and that growth can occur at temperatures below freezing. Low temperature tolerance may allow this species to expand its range farther north than previously suspected

The GBS code for tokamak scrape-off layer simulations

We describe a new version of GBS, a 3D global, flux-driven plasma turbulence code to simulate the turbulent dynamics in the tokamak scrape-off layer (SOL), superseding the code presented by Ricci et al. (2012) [14]. The present work is driven by the objective of studying SOL turbulent dynamics in medium size tokamaks and beyond with a high-fidelity physics model. We emphasize an intertwining framework of improved physics models and the computational improvements that allow them. The model extensions include neutral atom physics, finite ion temperature, the addition of a closed field line region, and a non-Boussinesq treatment of the polarization drift. GBS has been completely refactored with the introduction of a 3-D Cartesian communicator and a scalable parallel multigrid solver. We report dramatically enhanced parallel scalability, with the possibility of treating electromagnetic fluctuations very efficiently. The method of manufactured solutions as a verification process has been carried out for this new code version, demonstrating the correct implementation of the physical model. (C) 2016 Published by Elsevier Inc

Scrape Off Layer (SOL) transport and filamentary characteristics in high density tokamak regimes

A detailed cross-device investigation on the role of filamentary dynamics in high density regimes has been performed within the EUROfusion framework comparing ASDEX Upgrade (AUG) and TCV tokamaks. Both devices have run density ramp experiments at different levels of plasma current, keeping toroidal field or q95 constant in order to disentangle the role of parallel connection length and the current. During the scan at constant toroidal field, in both devices SOL profiles tend to develop a clear Scrape Off Layer (SOL) density shoulder at lower edge density whenever current is reduced. The different current behavior is substantially reconciled in terms of edge density normalized to Greenwald fraction. During the scan at constant q95 AUG exhibits a similar behaviour whereas in TCV no signature of upstream profile modification has been observed at lower level of currents. The latter behaviour has been ascribed to the lack of target density roll-over. The relation between upstream density profile modification and detachment condition has been investigated. For both devices the relation between blob-size and SOL density e-folding length is found independent of the plasma current, with a clear increase of blob-size with edge density normalized to Greenwald fraction observed. ASDEX Upgrade has also explored the filamentary behaviour in H-Mode. The experiments on AUG focused on the role of neutrals, performing discharges with and without the cryogenic pumps, highlighting how large neutral pressure not only in the divertor but at the midplane is needed in order to develop a H-Mode SOL profile shoulder in AUG