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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Upright exercise or supine lower body negative pressure exercise maintains exercise responses after bed rest

Adaptation to bed rest or space flight is accompanied by an impaired ability to exercise in an upright position. We hypothesized that a daily, 30-min bout of intense, interval exercise in upright posture or supine against lower body negative pressure (LBNP) would maintain upright exercise heart rate and respiratory responses after bed rest. Twenty-four men (31 +/- 3 yr) underwent 5 d of 6 degree head-down tilt: eight performed no exercise (CON), eight performed upright treadmill exercise (UPex), and eight performed supine treadmill exercise against LBNP at -51.3 +/- 0.4 mm Hg (LBNPex). Submaximal treadmill exercise responses (56, 74, and 85% of VO2peak) were measured pre- and post-bed rest. In CON, submaximal heart rate, respiratory exchange ratio, and ventilation were significantly greater (P < or = 0.05) after bed rest. In UPex and LBNPex, submaximal exercise responses were similar pre- and post-bed rest. Our results indicate that a daily 30-min bout of intense, interval upright exercise training or supine exercise training against LBNP is sufficient to maintain upright exercise responses after 5 d of bed rest. These results may have important implications for the development of exercise countermeasures during space flight.

NASA Center ARC↗
NASA NTRSDate not supplied

Light as a Reactive Countermeasure to Sleep Inertia: Translating Laboratory Findings to the Field

Sleep inertia describes the brief period of impaired alertness, mood, and cognitive performance experienced after waking. Under laboratory settings, light exposure upon waking during a habitual sleep period has been shown to improve sleep inertia symptoms. We investigated whether a field-deployable light-emitting device would help to mitigate sleep inertia in a real-world setting. Thirty-six participants (18 female; 26.6 years ± 6.1) completed an at-home, within-subject, randomized crossover study. Participants followed their habitual sleep-wake schedule for five nights before an adaptation and experimental night. Forty-five minutes after bedtime on the experimental night, participants received a phone call and were instructed to wear light-emitting glasses with the light either on (light condition) or off (control). A 5-minute psychomotor vigilance task (PVT), the Karolinska Sleepiness Scale (KSS), visual analog scales of mood (VASmood), and a 3-minute descending subtraction task (DST) were performed starting at +2, +12, +22, and +32 minutes after the call. Participants then went back to sleep and were called 45 minutes after lights out for the opposite condition. A series of mixed-effects models were performed with fixed effects of test bout, condition, test bout × condition, and a random effect of participant. Covariates included pre-sleep baseline scores, randomization order, sex, and sleep history. Participants rated themselves as more alert and energetic in the light condition compared to the control condition (VASalert-sleepy p = .01; VASlethargic-energetic p = .001). There was no effect of condition for DST outcomes, but there was a significant improvement in DST total responses in the light condition in a subset of participants waking from N3 (p = .03). There was a significant effect of condition for PVT outcomes, with faster responses (p < .001) and fewer lapses (p < .001) in the control condition. Our results under naturalistic at-home settings suggest that, similar to the in-laboratory study findings, the light intervention improved subjective alertness and mood, while working memory improved after waking from N3. Future studies of light interventions should include measures of visual acuity and comfort to assess the full feasibility and efficacy of interventions in real-world environments.

sleep inertia↗
NASA NTRSDate not supplied

Light as a Reactive Countermeasure to Sleep Inertia: Translating Laboratory Findings to the Field

Sleep inertia describes the brief period of impaired alertness, mood, and cognitive performance experienced after waking. Under laboratory settings, light exposure upon waking during a habitual sleep period has been shown to improve sleep inertia symptoms. We investigated whether a field-deployable light-emitting device would help to mitigate sleep inertia in a real-world setting.\ Thirty-six participants (18 female; 26.6 years ± 6.1) completed an at-home, within-subject, randomized crossover study. Participants followed their habitual sleep-wake schedule for five nights before an adaptation and experimental night. Forty-five minutes after bedtime on the experimental night, participants received a phone call and were instructed to wear light-emitting glasses with the light either on (light condition) or off (control). A 5-minute psychomotor vigilance task (PVT), the Karolinska Sleepiness Scale (KSS), visual analog scales of mood (VASmood), and a 3-minute descending subtraction task (DST) were performed starting at +2, +12, +22, and +32 minutes after the call. Participants then went back to sleep and were called 45 minutes after lights out for the opposite condition. A series of mixed-effects models were performed with fixed effects of test bout, condition, test bout × condition, and a random effect of participant. Covariates included pre-sleep baseline scores, randomization order, sex, and sleep history. Participants rated themselves as more alert and energetic in the light condition compared to the control condition (VASalert-sleepy p = .01; VASlethargic-energetic p = .001). There was no effect of condition for DST outcomes, but there was a significant improvement in DST total responses in the light condition in a subset of participants waking from N3 (p = .03). There was a significant effect of condition for PVT outcomes, with faster responses (p < .001) and fewer lapses (p < .001) in the control condition. Our results under naturalistic at-home settings suggest that, similar to the in-laboratory study findings, the light intervention improved subjective alertness and mood, while working memory improved after waking from N3. Future studies of light interventions should include measures of visual acuity and comfort to assess the full feasibility and efficacy of interventions in real-world environments.

sleep inertia↗
NASA NTRSDate not supplied

Light as a Reactive Countermeasure to Sleep Inertia: Translating Laboratory Findings to the Field

Sleep inertia describes the brief period of impaired alertness, mood, and cognitive performance experienced after waking. Under laboratory settings, light exposure upon waking during a habitual sleep period has been shown to improve sleep inertia symptoms. We investigated whether a field-deployable light-emitting device would help to mitigate sleep inertia in a real-world setting. Thirty-six participants (18 female; 26.6 years ± 6.1) completed an at-home, within-subject, randomized crossover study. Participants followed their habitual sleep-wake schedule for five nights before an adaptation and experimental night. Forty-five minutes after bedtime on the experimental night, participants received a phone call and were instructed to wear light-emitting glasses with the light either on (light condition) or off (control). A 5-minute psychomotor vigilance task (PVT), the Karolinska Sleepiness Scale (KSS), visual analog scales of mood (VASmood), and a 3-minute descending subtraction task (DST) were performed starting at +2, +12, +22, and +32 minutes after the call. Participants then went back to sleep and were called 45 minutes after lights out for the opposite condition. A series of mixed-effects models were performed with fixed effects of test bout, condition, test bout × condition, and a random effect of participant. Covariates included pre-sleep baseline scores, randomization order, sex, and sleep history. Participants rated themselves as more alert and energetic in the light condition compared to the control condition (VASalert-sleepy p = .01; VASlethargic-energetic p = .001). There was no effect of condition for DST outcomes, but there was a significant improvement in DST total responses in the light condition in a subset of participants waking from N3 (p = .03). There was a significant effect of condition for PVT outcomes, with faster responses (p < .001) and fewer lapses (p < .001) in the control condition. Our results under naturalistic at-home settings suggest that, similar to the in-laboratory study findings, the light intervention improved subjective alertness and mood, while working memory improved after waking from N3. Future studies of light interventions should include measures of visual acuity and comfort to assess the full feasibility and efficacy of interventions in real-world environments.

sleep inertia↗
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↗
DOE OSTI2020

Impact of equilibrium radial electric field on energy loss process after pedestal collapse

An impact of the equilibrium radial electric field on energy loss processes after pedestal collapse is numerically investigated using the BOUT++ framework. Using linear stability analysis, the resistive ballooning mode is shown to be stabilized by the radial shear of the equilibrium radial electric field. On the other hand, the energy loss level after the pedestal collapse increases if the equilibrium radial electric field is taken into account. The spatio-temporal and phase diagram analyses show that the equilibrium radial electric field partially cancels the fluctuation-driven toroidally axisymmetric radial electric field and weakens the E × B shearing rate after pedestal collapse, weakening the turbulence suppression by vortex shearing. Finally, the equilibrium radial electric field therefore increases turbulence intensity in nonlinear cyclic oscillations among pressure gradient, E × B shearing rate, and turbulence intensity, which gives rise to subsequent bursts of turbulent transport and increases the energy loss level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Pedestal origin and extrapolation of high-density small edge-localised-modes peak parallel energy fluence in ITER and SPARC

Experimental analysis and simulations with the BOUT++ code show that small edge-localised modes (ELMs) in reactor-relevant high-density regimes originate in a region close to the separatrix and only marginally perturb the pedestal structure. The measured divertor peak parallel energy fluence (ε ∥,peak ) for a database of small ELM scenarios in DIII-D and ASDEX Upgrade can be reproduced, within 40 % accuracy on average, if an ad hoc modification of the Eich peak parallel ELM energy fluence model is applied to account for the small ELM pedestal birth location. This allows for first-order extrapolation of small-ELM divertor ε ∥,peak to ITER and SPARC, resulting in values that satisfy the nominal melting threshold of tungsten monoblocks of 12 MJ m −2 . The findings reported in this study, both via modelling and direct measurements, constitute a step forward in assessing small ELMs in high edge-collisionality scenarios as a viable plasma regime for the operation of next-generation fusion machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The Remarkable 2024 North Atlantic Mid‐Season Hurricane Lull

Abstract The 2024 North Atlantic (hereafter Atlantic) hurricane season started quickly, with the earliest Category 5 on record (Beryl) and three hurricanes forming through 14 August. Following Ernesto's dissipation on 20 August, the Atlantic hurricane season became extremely quiet during the climatological peak of hurricane season, with only one Category 2 hurricane (Francine) and one tropical storm through 23 September. Several environmental factors likely contributed to this unexpected, prolonged lull. During mid‐to‐late August, subseasonal conditions were broadly favorable for Atlantic hurricanes, but a northward shift in African easterly wave emergence latitude yielded fewer tropical cyclone seed disturbances that also traversed unfavorably cool ocean water. During early‐to‐mid September, subseasonal variability driven by the Madden‐Julian oscillation was less conducive to hurricane activity, with several bouts of increased vertical wind shear across the central Atlantic. Throughout most of the lull, the tropical Atlantic was anomalously dry and subsident, suppressing hurricane formation chances. Plain Language Summary The 2024 North Atlantic (hereafter Atlantic) hurricane season was anticipated to be extremely active, with seasonal forecasts calling for a potentially historic season. Following a busy start including the earliest Category 5 Atlantic hurricane on record (Beryl), the season became very quiet between 20 August–23 September, with only one Category 2 hurricane (Francine) and one tropical storm forming. Following the lull, the Atlantic became extremely busy again, with seven hurricanes forming from 25 September through the end of the season–the most on record for this period. This paper focuses on the surprising lull and attributes it to several different factors, including a late‐August northward shift in African easterly wave tracks (organized thunderstorm complexes spanning thousands of kilometers that often serve as Atlantic hurricane seeds). This northward shift brought these systems out over colder water and into a drier airmass, suppressing their hurricane formation chances. During the first 3 weeks of September, a large‐scale tropical phenomenon known as the Madden‐Julian oscillation, resulted in increased sinking motion, drier mid‐levels and increased vertical wind shear–all of which typically reduce Atlantic hurricane activity. Overall, the lull demonstrates how seed disturbances and other environmental factors can strongly modulate aggregate seasonal tropical cyclone activity. Key Points Following an active start, the 2024 Atlantic hurricane was quiet during the climatological peak with only 1 tropical storm and 1 hurricane The lull was not anticipated by seasonal forecasting agencies, which called for an extremely active 2024 Atlantic hurricane season Factors reducing storm count included a northward shift in African easterly waves, broad‐scale subsidence, and reduced mid‐level moisture

Klotzbach, P J↗
DOE OSTI2022

Modeling study of divertor particle and heat flux asymmetries for EAST H-mode discharges

The BOUT++ transport code is run to study the effects of plasma drifts on the divertor out-in asymmetries (DOIAs) of particle and heat fluxes and their decay widths for EAST lower single null H-mode discharges. The diamagnetic drift seems to have no effects on the DOIAs of total particle and heat fluxes due to its divergence-free nature. However, it could significantly increase the DOIAs of peak particle and heat fluxes and the flux decay widths. The E × B drift is found to induce a large plasma flow to the divertor region, enhancing the DOIAs of both total and peak particle and heat fluxes and the flux decay widths. Both the radial and poloidal components of the E × B drift are necessary in increasing the DOIAs, however, the radial E × B drift seems to play a more important role. The effects on the DOIAs caused by both diamagnetic and E × B drifts are reversed with the reverse of toroidal magnetic field. The heat flux decay width λ q and spreading width S q are important physical and engineering parameters for the divertors and could be obtained by fitting the heat flux profiles at divertor targets. The λ q at the outer target from the simulation case with all drifts could well match with the multi-machine scaling proposed by Eich and the DOIA of λ q is in reasonable agreement with the scaling proposed by Goldston.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulations of heat fluxes in an ELMy H-mode discharge on HL-2A

In order to study the distribution and evolution of the transient heat flux on HL-2A during edge-localized-mode (ELM) bursts, the BOUT++ electromagnetic six-field two-fluid model is used to simulate the pedestal collapse under the lower single-null divertor geometry. The equilibrium profiles of HL-2A ELMy H-mode discharge No. 24 953 are adopted as the initial condition in the original case. In this instance, linear analysis shows that the resistive ballooning mode (RBM) and drift-Alfven wave are unstable to this equilibrium, and RBM is the dominant instability. The evolutions of the radial heat fluxes at the outer mid-plane and heat fluxes to the inner and outer targets during the ELM event are presented. Six more equilibria are constructed based on the original case to find out the influence of the pedestal profiles on the peak electron heat flux. The results indicate that the heat flux increases with temperature and/or density, and the theoretical analysis and simulation results consistently show that the heat flux q ∥e is proportional to n e0,SEP $T^{\frac32}_{e0,SEP}$.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Two-stage crash process in resistive drift ballooning mode driven ELM crash

We report a two-stage crash process in edge localized mode (ELM) driven by resistive drift-ballooning modes (RDBMs) numerically simulated in a full annular torus domain with a scale-separated four-field reduced MHD (RMHD) model using the BOUT++ framework. In the early nonlinear phase, the small first crash is triggered by linearly unstable RDBMs, and m/n=2/1 magnetic islands are nonlinearly excited by nonlinear coupling of RDBMs as well as their higher harmonics. Here, m is the poloidal mode number, n is the toroidal mode number, the q = 2 rational surface exists near the pressure gradient peak, and q is the safety factor. Simultaneously, middle-n RDBM turbulence develops but is poloidally localized around X-points of the magnetic islands, leading to the small energy loss. The second large crash occurs in the late nonlinear phase. Higher harmonics of m/n=2/1 magnetic islands well develop around the q = 2 surface via nonlinear coupling and make the magnetic field stochastic by magnetic island overlapping. Turbulence heat transport develops at X-points of higher harmonics of m/n=2/1 magnetic islands, resulting in the turbulence spreading in the poloidal direction. The large second crash is triggered when the turbulence covers the whole poloidal region so that the magnetic island generation and magnetic field stochastization before the large crash can be interpreted as ELM precursors. It is concluded that the ELM trigger is attributed to the turbulent spreading in the poloidal direction in synchronization with the magnetic field stochastization and the crash is driven by E × B convection rather than the conventional Rechester–Rosenbluth anomalous electron heat transport.

Seto, H. (ORCID:0000000259278643)↗
DOE OSTI2024

Theoretical and global simulation analysis of collisional microtearing modes

Microtearing modes (MTMs) are suggested as a candidate for anomalous thermal transport in tokamak H-mode discharges. This study investigates MTMs in tokamak plasmas, employing simulations in the BOUT++ framework. It simplifies and linearizes the governing equations in detailed linear simulations. The study meticulously evaluates various conductivity models under diverse plasma conditions and collision regimes. The research thoroughly assesses different conductivity models across a range of plasma conditions and collision regimes. A unified dispersion relation that includes both MTM and Drift-Alfvén Wave (DAW) instabilities is derived, showing that DAW and MTM instabilities occur at varying distances from the rational surface. Specifically, MTMs become unstable near the rational surface but stabilize farther away, while drift-Alfvén instability appears farther from the rational surface. The study also re-derives MTM dispersion relations using Ohm's law and the vorticity equation, providing a thorough analysis of electromagnetic and electrostatic interactions in tokamaks. Global simulations demonstrate an inverse correlation between MTM growth rates and collisionality, and a direct correlation with temperature gradients. The nonalignment of the rational surface with the peak ω*e stabilizes the MTMs. Nonlinear simulations highlight electron temperature relaxation as the primary saturation mechanism for MTMs, with magnetic flutter identified as the dominant mode of electron thermal transport.

Fan, K. (ORCID:0000000227518809)↗