Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BOUT++”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

STORM: Scrape-off layer turbulence in tokamak fusion reactors

The scrape-off layer of a tokamak fusion reactor carries the plasma exhaust from the hot core plasma to the material surfaces of the reactor vessel. The heat loads imposed by the exhaust are a critical limit on the performance of fusion power plants. Turbulent transport of the plasma regulates the width of the scrape-off layer plasma and must be modelled to understand the intensity of these heat loads. STORM is a plasma turbulence code capable of simulating three dimensional turbulence across the full scrape-off layer of a tokamak fusion reactor, using a drift reduced, collisional fluid model. STORM uses mostly finite difference schemes, with a staggered grid in the direction parallel to the magnetic field. We describe the model, geometry and initialisation options used by STORM, as well as the numerical methods, which are implemented using the BOUT++ plasma simulation framework. BOUT++ has been enhanced alongside the development of STORM, providing better support for staggered grid methods. We summarise these enhancements, including a detailed explanation of the parallel derivative methods, which underwent a major update for version 4 of BOUT++.

BOUT++↗
DOE OSTI2024

Thermally unstable roosts influence winter torpor patterns in a threatened bat species

Abstract Many hibernating bats in thermally stable, subterranean roosts have experienced precipitous declines from white-nose syndrome (WNS). However, some WNS-affected species also use thermally unstable roosts during winter that may impact their torpor patterns and WNS susceptibility. From November to March 2017–19, we used temperature-sensitive transmitters to document winter torpor patterns of tricolored bats (Perimyotis subflavus) using thermally unstable roosts in the upper Coastal Plain of South Carolina. Daily mean roost temperature was 12.9 ± 4.9°C SD in bridges and 11.0 ± 4.6°C in accessible cavities with daily fluctuations of 4.8 ± 2°C in bridges and 4.0 ± 1.9°C in accessible cavities and maximum fluctuations of 13.8 and 10.5°C, respectively. Mean torpor bout duration was 2.7 ± 2.8 days and was negatively related to ambient temperature and positively related to precipitation. Bats maintained non-random arousal patterns focused near dusk and were active on 33.6% of tracked days. Fifty-one percent of arousals contained passive rewarming. Normothermic bout duration, general activity and activity away from the roost were positively related to ambient temperature, and activity away from the roost was negatively related to barometric pressure. Our results suggest ambient weather conditions influence winter torpor patterns of tricolored bats using thermally unstable roosts. Short torpor bout durations and potential nighttime foraging during winter by tricolored bats in thermally unstable roosts contrasts with behaviors of tricolored bats in thermally stable roosts. Therefore, tricolored bat using thermally unstable roosts may be less susceptible to WNS. More broadly, these results highlight the importance of understanding the effect of roost thermal stability on winter torpor patterns and the physiological flexibility of broadly distributed hibernating species.

Biodiversity & Conservation↗
DOE OSTI2022

Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D

The wide-pedestal quiescent high confinement mode discovered on DIII-D in recent years is a stationary and quiescent H-mode (QH-mode) with the pedestal width exceeding EPED prediction by at least 25%. Its characteristics, such as low rotation, high energy confinement and edge localized mode-free operation, make it an attractive operation mode for future reactors. Linear and nonlinear simulations using BOUT++ reduced two fluid MHD models and awere carried out to investigate the bursty broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling–ballooning (PB) mode γ PB < 0.04ω A ) located near the minimum in E r</:sub> well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave locates at smaller radius compared to E r</:sub> well propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitative agreement is likely due to the lack of trapped electron physics in current fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Characteristics of grassy ELMs and their impact on the divertor heat flux width

We report BOUT++ turbulence simulations are conducted for a 60 s steady-state long pulse high β p EAST grassy ELM discharge. BOUT++ linear simulations show that the unstable mode spectrum covers a range of toroidal mode numbers from low-n (n = 10–15) peeling–ballooning modes (P–B) to high-n (n = 40–80) drift-Alfvén instabilities. Nonlinear simulations show that the ELM crash is triggered by low-n peeling modes and fluctuation is generated at the peak pressure gradient position and radially spread outward into the scrape-off-layer, even though the drift-Alfvén instabilities dominate the linear growth phase. However, drift-Alfvén turbulence delays the onset of the grassy ELM and enhances the energy loss with the fluctuation extending to pedestal top region. Simulations further show that if the peeling drive is removed, the fluctuation amplitude drops by an order of magnitude and the ELM crashes disappear. The divertor heat flux width is ~2 times larger than the estimates based on the HD model and the Eich's ITPA multi-tokamak scaling (or empirical Eich scaling) due to the strong radial turbulence transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Modelling of nanometer scale dust grains in tokamak

Dust poses a serious threat to tokamak operation and safety. It is important to study the behaviour of dust grains under tokamak's discharge conditions, which depends heavily on their size and charge. Existing simulations mainly address issues on dust grains with radii larger than 1 μm, in which case, the drift effect due to electromagnetic fields can be safely ignored. For nanometer scale dust grains, however, the drift effect becomes significant and a new model based on guiding-centre system needs to be established. In this work, the NDS has been done under BOUT++ framework. The simulation contains two parts. Part one, NDS evaluates the charging and ablation processes of the dust grains. In the second part, the guiding-centre orbits of dust particles are tracked in tokamak plasmas, whose parameters are obtained from BOUT++, a highly desirable C++ code package for performing parallel plasma fluid simulations with an arbitrary number of equations in 3D curvilinear coordinates. The orbit of nanodust dynamics is described by guiding centre equations for simplicity, and these equations are numerically solved by conventional fourth-order Runge Kutta method. Simulations provide results such as trajectories and evolutions of dust particles with different sizes and velocities for different tokamak geometries. Results show tungsten dust grains with a radius of a few nanometers launched from outer midplane will oscillate before totally ablated in C-Mod. The oscillation in this case is driven by the ion drag force. Larger Nanodust with a radius of 100 nm, on the contrary, cannot be completely constrained by the electromagnetic field. In conclusion, the high plasma temperature and density in the seperatrix region causes severe dust ablation, resulting in total ablation within several ms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Turbulence spreading effects on the ELM size and SOL width

BOUT++ turbulence simulations were performed to investigate the impact of turbulence spreading on the edge localized mode (ELM) size and divertor heat flux width (λ q ) broadening in small ELM regimes. Here, this study is motivated by EAST experiments. BOUT++ linear simulations of a pedestal radial electric field (E r ) scan show that the dominant toroidal number mode (n) shifts from high-n to low-n, with a narrow mode spectrum, and the maximum linear growth rate increases as the pedestal E r well deepens. The nonlinear simulations show that as the net E × B pedestal flow increases, the pressure fluctuation level and its inward penetration beyond the top of the pedestal both increase. This leads to a transition from small ELMs to large ELMs. Both inward and outward turbulence spreading are sensitive to the scrape-off-layer (SOL) plasma profiles. The inward turbulence spreading increases for the steep SOL profiles, leading to increasing pedestal energy loss in the small ELM regime. The SOL width (λ q ) is significantly broadened progressing from the ELM-free to small ELM regime, due to the onset of strong radial turbulent transport. The extent of the SOL width (λ q ) broadening depends strongly on outward turbulence spreading. The fluctuation energy intensity flux Γ ε at the separatrix can be enhanced by increasing either pedestal E r flow shear or local SOL pressure gradient. The λ q is broadened as the fluctuation energy intensity flux Γ ε at the last close flux surface (LCFS) increases. Local SOL E × B flow shear will restrain outward turbulence spreading and the associated heat flux width broadening. Operating in H-mode with small ELMs has the potential to solve two critical problems: reducing the ELM size and broadening the SOL width.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Impact of pedestal density gradient and collisionality on ELM dynamics

BOUT++ turbulence simulations are conducted to capture the underlying physics of small ELM characteristics achieved by increasing separatrix density via controlling strike points from vertical to horizontal divertor plates for three EAST discharges. BOUT++ linear simulations show that the most unstable modes change from high-n ideal ballooning modes to intermediate-n peeling–ballooning modes and eventually to peeling–ballooning stable plasmas in the pedestal. Nonlinear simulations show that the fluctuation is saturated at a high level for the lowest separatrix density. The ELM size decreases with increasing separatrix density, until the fraction of this energy lost during the ELM crash becomes less than 1% of the pedestal stored energy, leading to small ELMs. Simulations indicate that small ELMs can be triggered either by the marginally peeling–ballooning instability near the peak pressure gradient position inside the pedestal or by a local instability in the pedestal foot with a larger separatrix density gradient. The pedestal collisionality scan for type-I ELMs with steep pedestal density gradient shows that both linear growth rate and ELM size decrease with increasing collisionality. On the contrary, the pedestal collisionality and pedestal density width scan with a weak pedestal density gradient indicate small ELMs can either be triggered by a high-n ballooning mode or by a low-n peeling mode in a low collisionality region 0.04–0.1. We report the simulations indicate the weaker the linear unstable modes near marginal stability with small linear growth rate, the lower nonlinearly saturated fluctuation intensity and the smaller turbulence spreading from the linear unstable zone to stable zone in the nonlinear saturation phase, leading to small ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Edge localized mode characteristics and divertor heat flux during stationary and transient phase for CFETR hybrid scenario

The study of edge localized mode (ELM) behavior, stationary heat flux and transient heat flux during the ELM crash phase is performed for a China Fusion Engineering Test Reactor (CFETR) 1 GW hybrid mode operation scenario ($R$ = 7.2 m, ${B_{\text{T}}}$ = 6.5 T, ${I_{\text{p}}}$ = 13.78 MA). Modeling and simulation start with a scenario obtained by multi-code integrated modeling on the one modeling framework for integrated tasks framework. Linear stability and nonlinear simulations of ELM dynamics are carried out using the BOUT++ six-field reduced magnetohydrodynamic module, which show a much smaller ELM energy loss (Δ ELM ~ 0.13%) compared to that of a Type-I ELM. Parametric analysis of the weak linear growth rate and small ELM energy loss characteristics shed light on physics corresponding to a grassy ELM regime for CFETR 1 GW hybrid scenario. The transient heat flux on the divertor target during this small ELM phase is investigated using BOUT++. Here we found that upstream radial transport in the scrape-off-layer (SOL) induced by small ELMs is weak, which keeps it in the drift-dominated region. However, the heat flux width is still broadened to ${\lambda _{\text{q}}}$ = 4.64 mm by the increase of separatrix temperature during ELM nonlinear evolution. The impact of transient peak heat load and ELM energy fluence on tungsten melting and net erosion rate of divertor target is evaluated for the first-time using physics-based transport that connects the pedestal with the SOL. Energy fluence caused by a single ELM pulse is below the tungsten melting limit, while tungsten erosion would exceed the material requirements. We conclude that external mitigation methods, such as divertor detachment and advanced divertor geometry are likely needed for the steady state operation of CFETR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental and theoretical study of weakly coherent mode in I-mode edge plasmas in the EAST tokamak

The I-mode is a promising operation mode for fusion in the future, featuring high-temperature and low-density confinement, but the reason why the temperature and density are decoupled remains an important aspect to be explored. The experimental results from the experimental advanced superconducting tokamak (EAST) showed that the weakly coherent mode (WCM) is directly related to sustaining the I-mode and that the peak amplitude of the WCM is proportional to the temperature in the pedestal. Here, simulating the experimental data from EAST with the six-field model of BOUT++, we find a density perturbation close to the frequency of the WCM observed in the experiment. By testing all the physical terms in this model, we find that the density perturbation and particle transport are directly related to the drift Alfvén wave (DAW) mode. Additionally we use the SymPIC program (Xiao et al Plasma Sci. Technol. 20 110501; Phys. Plasmas 22 112504; Plasma Sci. Technol. 23 055102) to simulate the same experimental data and find that the frequency range of the WCM is close to both experimental and BOUT++ results. Therefore, the WCM of the I-mode can be considered to be driven by the DAW, which helps improve the transport of the I-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Understanding the formation of a low-pressure pedestal in the presence of a strong internal transport barrier in DIII-D high β p plasmas

As a promising scenario for fusion reactors, the high poloidal-beta (β P ) scenario is characterized by a strong large radius internal transport barrier (ITB), which significantly enhances the overall confinement quality and the bootstrap current fraction for fully non-inductive operation. It is frequently observed that in the presence of a strong ITB, the pedestal height is lower and is accompanied by small edge localized modes (ELMs), which further improves the compatibility of a high performance core with an edge solution. A mechanism for the formation of the low pedestal is proposed in this paper. It is found that the strong ITB creates an off-axis bootstrap current to clamp the local safety factor q, and thus the magnetic shear in the outer core/pedestal region is increased. Gyrokinetic simulations with the CGYRO code show that the higher magnetic shear brings the experimental profiles into the range where the growth rate of drift-wave instabilities and thus transport is higher, and therefore a lower pedestal gradient is expected. Here, the combination of low pedestal and high magnetic shear further enhances the turbulent transport across the whole pedestal, consistent with power balance analysis. Such a positive feedback mechanism ultimately results in a lower pressure pedestal as observed in experiments. Under such a low pedestal, linear simulations with BOUT++ predict the growth rates of peeling–ballooning modes to be lower across the whole toroidal mode number spectra, and the nonlinear BOUT++ simulation exhibits lower saturated fluctuation intensity as well, consistent with the experimentally observed lower ELM size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Turbulence simulations of drift-Alfvén wave effects on edge-localized modes and divertor target heat flux width

A series of BOUT++ turbulence simulations are performed for two Pre-Fusion Power Operation (PFPO-1,2) phases of the ITER Research plan proclaimed in 2019. Linear simulation results show that PFPO-1 is unstable to Peeling-ballooning modes (PB), while the PFPO-2 is unstable to the coupling of PB and Drift-Alfvén wave (DAW) instabilities. The linear results are qualitatively consistent with the dispersion relation of DAW. To study the influence of the DAW on the edge-localized mode (ELM) crash, the PFPO-2 is adopted in this section. Different from the grassy ELM in [Li et al 2022 Nucl. Fusion 62, 096030], nonlinear simulations show that the ELM size of PFPO-2 is almost one-third of the grassy ELM, representing a distinct small ELM. However, simulations then show that if the PB instability is removed, the fluctuation amplitude drops by an order of magnitude and the ELM crash disappears, which is in accord with the theory in [Xu et al 2010 Physical Review Letters, Vol. 105, 175005] and the results in [Li et al 2022 Nucl. Fusion 62, 096030], confirming that the PB instability is a necessary condition for an ELM crash. Furthermore, removing the DAW drive also suppresses ELM crashes, implying that PB instability is necessary but insufficient for PFPO-2 ELM and that DAW could amplify PB-driven turbulence. In addition, by integrating results of both PFPO phases by BOUT++ turbulence analysis with heat flux width (λ q ) and electron thermal diffusivity (χ e ), these simulations establish a predictive framework for heat flux width and its dominant mechanisms, thereby providing a comprehensive scheme for the prediction of the heat flux width. Moreover, simulations indicate that DAW driving can increase the transport coefficient by enhancing the turbulent transport, leading to a broadened heat flux width once the transport coefficient exceeds

Drift-Alfvén wave↗
DOE OSTI2025

Q2 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q2: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis: $\circ$ Run CESOL with BOUT++/Hermes-3 and immersed boundary condition to directly map to wall: • Run BOUT++/Hermes-3 through the IPS workflow to find radial particle and energy diffusivities to match either the Eich or the physics-based scaling of the SOL heat flux width, and • Expand source of first wall heat flux to include charged particles, neutrals, and radiation from the core+edge. 2. Generate medium fidelity parametrized CAD: $\circ$ Develop the TRACER tool to read an existing CAD, regenerate the geometry based on vertex location and connectivity information, define vertex translation and parameters needed for scaling the CAD, and $\circ$ Utilize the FreeGS code to determine CAT PF coil placement, including minimizing the number of coils, coil current, and electromechanical stresses. 3. Utilize plasma loading for engineering analysis: $\circ$ Couple the plasma loading to input for OpenFOAM and demonstrate initial test of thermal analysis of CAT first wall loading with typical DCLL blanket component cooling boundary conditions. 4. Demonstrate nuclear analysis: $\circ$ Apply initial analysis of tritium transport in DCLL blanket by evaluating spatially resolved tritium generation rates, tritium diffusion and convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Experimental and simulation study of target biasing effects on plasma transport in linear plasma device MPS-LD

Linear plasma devices (LPDs) are important experimental platforms for investigating plasma–material interactions (PMI). In PMI experiments, it has been found that applying a target bias not only effectively modifies the incident ion energy, but also induces significant changes in the electron density and electron temperature, whereby the evolution of these plasma parameters is primarily governed by plasma transport processes. However, at present, the physical process and mechanism underlying such bias-induced variations remain unclear. In this work, biasing experiments under argon plasma discharge conditions were first carried out on the MPS-LD device. For the corresponding experiments, an electric potential model was newly developed based on the BOUT++ LPD module, enabling self-consistent simulations of plasma transport under biased conditions. Numerical simulations were then performed to reproduce the experimental results and to validate the accuracy of the proposed model. Finally, by combining experimental measurements with numerical simulations, a bias-voltage scan was performed to investigate how the electron density and electron temperature vary with the bias voltage (U bias ). The results show that applying negative bias decreases the target electron density (n e,T ) while increasing the target electron temperature (T e,T ). In contrast, positive bias increases both n e,T and T e,T ; however, at high positive bias, n e,T first reaches a maximum and subsequently decreases with further increases in U bias . The underlying physical mechanisms are analyzed using particle flux, momentum, and energy conservation. It indicates that the applied bias regulates the parallel electric field, thereby changing ion and electron velocities, and consequently affecting the electron density. At high positive bias, the ion velocity is further influenced by ion viscosity, leading to the reversal in n e,T . Meanwhile, the enhanced parallel electric field drives stronger currents, significantly increasing ion–electron frictional work and converting the input bias power into electron energy, which raises the electron temperature. In conclusion, these results contribute to a deeper understanding of the effects and mechanisms of biasing on plasma transport in the MPS-LD device.

BOUT++ simulation↗
DOE OSTI2025

Investigation of divertor detachment induced through neon seeding and density ramp on HL-3

A new self-consistent 1D scrape-off layer model has been recently developed in BOUT++ framework, named SD1D, which includes equations for various particle species (e.g. main plasma, neutrals and impurities) and couples open databases like ADAS and AMJUEL. It is able to quickly and effectively simulate divertor detachment experiments. In this work, a typical detachment experiment (shot #6270) on HL-3 with neon seeding is simulated using the SD1D code. It is found that the target electron temperature and the target ion saturation current in the simulations are consistent with experimental results measured by Langmuir probes on the target plate. The variation of D α radiation intensity in the divertor is qualitatively similar to the measured D α signal. Following the experimental validations, different upstream densities are set in the simulations to study the impurity distribution under different plasma density conditions. It is found that increasing upstream density can be helpful for the control of the neon radiation front (closer to the target). In this work we also compare two detachment regimes in simulations. Based on the same initial experimental parameters (shot #6270) on HL-3, a scan of upstream density and a scan of neon seeding rate are carried out respectively. It is found that the role of atomic and molecular processes is different in the two detachment regimes. The current density roll-over is ascribed to a drop in the divertor ion source, and the variation of D α radiation intensity via different excitation channels is associated with the relevant collisional reaction sources.

BOUT++↗
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 OSTI2020

Impact of plasma density/collisionality on divertor heat flux width

Both ASDEX-Upgrade (AUG) data and the generalized HD (GHD) model showed that the scrape-off width broadens as the density/collisionality increases [1, 2]. A series of BOUT++ transport simulations are performed to study the physics of the scaling characteristics of the divertor heat flux width vs density/collisionality via a plasma density scan with either fixed pressure profile or fixed temperature profile inside separatrix. Additionally, the simulations show that even in the drift dominated regime, the divertor heat flux width can be broadened due to the transition of the SOL residence time from the parallel particle flow time to the enhanced parallel conduction time as the collisionality/density increases as posited in the GHD model. In addition, the heat flux width is found to be proportional to the square root of ion mass for low collisionality while it has a weakly dependence on ion mass for high collisionality. Furthermore, our simulations show that as the density increases, the radial electric field (E r ) well shallows, which potentially weakens E r × B flow shear stabilization of turbulence at high density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effects of radial transport on divertor power and particle flux widths under different operational regimes in EAST

Here, a study of the effects of radial transport on the outer divertor particle and power flux widths (λ js and λ q ) is carried out using plasmas under different operational regimes in the Experimental Advanced Superconducting Tokamak (EAST). In the EAST experiments, the λ js values measured during the grassy intra- and inter-ELM phases are similar, and the averaged value of λ js during grassy ELMy discharges is found to be smaller than that of the intra-ELM phase of type-I ELMy discharges but larger than that of the inter-ELM phase of type-I ELMy discharges. Simulations of scans of the radial particle and heat transport coefficients, D and χ e , performed using the BOUT++ transport code show that both λq and λjs increase with D and χe, especially when D or χ e is larger than a threshold value, indicating that background turbulence starts to matter in the determination of λ q and λ js when D or χ e surpasses the threshold value. A comparison between the simulation of the D scan and the experimental results shows that the different values of λ js obtained under different plasma operational regimes are probably due to the different intensities of background turbulence, which is beyond the scope of Goldston's heuristic drift-based model as reported by Goldston et al (2012 Nucl. Fusion 52 013009), since their model was derived by assuming that turbulent transport does not exceed a certain level. The ratio of λ js to λ q is greatly affected by the radial transport; however, simulations carried out using Scrape-off Layer Plasma Simulation (SOLPS) show that divertor geometry and plasma density have big influences on λ js , and thus could also affect the ratio of λ js to λ q .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Edge-localized-mode simulation in CFETR steady-state scenario

We report the EPED1 model and self-consistent core-pedestal coupling in integrated modeling are used to design the pedestal structure of the China Fusion Engineering Testing Reactor (CFETR) steady-state scenario. The key parameters, such as β p and q 95 , are based on the grassy edge-localized-mode (ELM) experimental database. In this work, we use the BOUT++ six-field two-fluid code to simulate the onset of the ELM in the CFETR steady-state scenario. The ELM size is around 0.2% in nonlinear simulations, which is in the experimental range of the grassy ELM discharges, 0.1%–1% observed in multiple tokamak devices. Linear and nonlinear simulations show that the dominant high-n ballooning modes peak around n = 40. Compared to type-I ELM crashing dynamics, grassy ELM crashing has a smaller initial crash and is then followed by three phases of turbulence spreading, which are dominated by multi-modes, a high-n mode of n = 45 and low-n mode of n = 5, respectively. In contras to type-I ELM, the perturbation of the high-n mode has a narrow width around ψ = 0.95, and magnetic island formation and reconnection occur only beyond ψ = 0.95, leading to a small initial crash. Mode–mode interaction in the multi-mode coexistence stage stops the growth of individual modes and reduces the transport of particles and heat, and these are the two reasons why the ELM size is small. In–out asymmetry of transient heat flux with a ratio of E out /E in = 3.5 is found during grassy ELM crash. The rise and delay times of the heat flux match the calculation from the free-streaming model. To evaluate the erosion of the divertor target, the energy fluence at the outer divertor target is calculated, which is 0.029 MJ m -2 , 5.5 times smaller than the tungsten melting limit 0.16 MJ m -2 . The calculated energy fluency still follows the experimental scaling law from type-I ELM experiments. The fluctuation eddies in the toroidal direction show a filament structure at the outer mid-plane. Parallel heat flux patterns with a toroidal mode number n = 10 are found at the outer divertor with an amplitude of 680 MW m -2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗