Search NASA⌕ Search

SEARCH · Search NASA

Results for “peeling-ballooning mode”

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 OSTI2021

Predicting operational windows of ELMs suppression by resonant magnetic perturbations in the DIII-D and KSTAR tokamaks

A newly developed plasma response model, combining the nonlinear two-fluid MHD code TM1 and toroidal MHD code GPEC run in ideal mode, quantitatively predicts the narrow isolated q 95 windows (Δq 95 ~ 0.1) of ELM suppression by n = 1, 2 and 3 resonant magnetic perturbations (RMPs) in both DIII-D and KSTAR tokamaks across a wide range of plasma parameters. The key physics that unites both experimental observations and our simulations is the close alignment of essential resonant q-surfaces and the location of the top of the pedestal prior to an ELM. This alignment permits an applied RMP to produce field penetration due to the lower E × B rotation at the pedestal top rather than being screened. The model successfully predicts that narrow magnetic islands form when resonant field penetration occurs at the top of pedestal, and these islands are easily screened when q 95 moves off resonance, leading to very narrow windows of ELM suppression (typically Δq 95 ~ 0.1). Furthermore, the observed reduction in the pedestal height is also well captured by the calculated classical collisional transport across the island. We recover observed q 95 , β N and plasma shape dependence of ELM suppression due to the effect of magnetic islands on pedestal transport and Peeling-Ballooning- Mode (PBM) stability. Importantly, experiments do occasionally observe wide windows of ELM suppression (Δq 95 > 0.5). Our model reveals that at low pedestal-top density multiple islands open, leading to wide operational windows of ELM suppression consistent with experiment. The model indicates that wide q95 windows of ELM suppression can be achieved at substantially higher pedestal pressure with less confinement degradation in DIII-D by operating at higher toroidal mode number (n = 4) RMPs. As a result, this can have significant implications for the operation of the ITER ELM control coils for maintaining high confinement together with ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Isotope dependence of the type I ELMy H-mode pedestal in JET-ILW hydrogen and deuterium plasmas

The pedestal structure, edge transport and linear MHD stability have been analyzed in a series of JET with the ITER-like wall hydrogen (H) and deuterium (D) type I ELMy H-mode plasmas. The pedestal pressure is typically higher in D than in H at the same input power and gas rate, with the difference mainly due to lower density in H than in D (Maggi et al (JET Contributors) 2018 Plasma Phys. Control. Fusion 60 014045). A power balance analysis of the pedestal has shown that higher inter-ELM separatrix loss power is required in H than in D to maintain a similar pedestal top pressure. This is qualitatively consistent with a set of interpretative EDGE2D-EIRENE simulations for H and D plasmas, showing that higher edge particle and heat transport coefficients are needed in H than in D to match the experimental profiles. It has also been concluded that the difference in neutral penetration between H and D leads only to minor changes in the upstream density profiles and with trends opposite to experimental observations. This implies that neutral penetration has a minor role in setting the difference between H and D pedestals, but higher ELM and/or inter-ELM transport are likely to be the main players. The interpretative EDGE2D-EIRENE simulations, with simultaneous upstream and outer divertor target profile constraints, have indicated higher separatrix electron temperature in H than in D for a pair of discharges at low fueling gas rate and similar stored energy (which required higher input power in H than in D at the same gas rate). The isotope dependence of linear MHD pedestal stability has been found to be small, but if a higher separatrix temperature is considered in H than in D, this could lead to destabilization of peeling-ballooning modes and shrinking of the stability boundary, qualitatively consistent with the reduced pedestal confinement in H.

08 HYDROGEN↗
DOE OSTI2025

X-point effects on the ideal MHD modes in tokamaks in the description of dual-poloidal-region safety factor

The flux coordinates with dual-region safety factor (q) in the poloidal direction are developed in this work. The X-point effects on the ideal MHD modes in tokamaks are then analyzed using this coordinate system. Since the X-point effects mainly affect the edge region, the modes localized at the tokamak edge are particularly examined. Two types of modes are studied. The first is related to the conventional peeling or peeling-ballooning modes. The mode existence aligned with the local magnetic field in the poloidally core region as observed experimentally is confirmed. The X points are shown to contribute to a stabilizing effect for the conventionally treated modes with the surface-averaged q and with the tokamak edge portion truncated. The other is the axisymmetric modes localized in the vicinity of X points, which can affect the cross-field-line transport near the X points. The existence of axisymmetric modes points to the possibility of applying a toroidally axisymmetric resonant magnetic perturbation (RMP) in the X-point area for mitigating the edge localized modes, which can be an alternative to the current RMP design. The dual q description also has important implications for the existing non-axisymmetric RMP concept. It helps to understand why the RMP suppression of edge localized modes is difficult to achieve in the double-null tokamak configurations and points to the possibility of further improving the current RMP concept by considering the alignment to the local q.

Fourier analysis↗
DOE OSTI2020

Grassy ELM regime at low pedestal collisionality in high-power tokamak plasma

Spontaneous mitigation of giant ELMs and appearance of grassy ELMs have been observed repeatedly at low pedestal collisionality (${\nu }_{\mathrm{e},\mathrm{p}\mathrm{e}\mathrm{d}}^{{\ast}}\sim 0.15$) in the high-power (P inj > 13 MW) hybrid scenario in DIII-D. Higher β p and higher q 95 appear to be beneficial to achieving the grassy ELM regime. The grassy ELM H-mode plasma shows high energy confinement performance (H 98y2 up to over 1.5) under the conditions of high neutral beam torque and high core rotation. The pedestal width appears to exceed the EPED1.0 model prediction by more than 50%. Pedestal stability analysis performed with the ELITE code indicates that the stability against low-n kink/peeling modes is improved with increased auxiliary heating power and the operational point in the grassy ELM regime is located near the ballooning boundary. The pedestal stability characteristics during the grassy ELM crashes have been investigated in comparison with the giant ELM crashes based on plasma profiles experimentally measured with high resolution and accuracy. It has been found that the underlying mechanism for the observed small-amplitude ELM crashes is mainly the expansion of the ballooning stability boundary induced by an initial radially localized collapse in the pedestal, which helps to stop the growth of instabilities and further collapse of the pedestal. The effect of electron density pedestal on mitigating edge localized instabilities has been analyzed by numerical simulation, suggesting that the electron density pedestal characterized by high n e,sep /n e,ped and low density gradient helps to stabilize peeling-ballooning modes because of a low pressure pedestal gradient and to lower the ballooning boundary mainly because of a low ion diamagnetic frequency in the pedestal region, thus triggering ballooning instabilities and producing the intrinsic grassy ELMs. Numerical simulation of the Chinese fusion engineering test reactor (CFETR) with the SOLPS code indicates that the separatrix density might be insensitive to the electron diffusivity in the pedestal region and increase with the power flowing from the core region to the edge region. Furthermore, pedestal stability analysis suggests that the flat density pedestal with high separatrix density obtained in the high-power plasma in CFETR would make the operational point close to the ballooning boundary, which is considered to help destabilize ballooning instabilities and facilitate the access to the grassy ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Importance of resistivity on edge-localized mode onset in spherical tokamaks

In this work, we explore the impact of realistic plasma resistivity on the linear stability of peeling-ballooning (PB) modes in tokamak discharges with low-aspect ratio. For this study we consider discharges that are subject to edge-localized modes (ELMs) in the National Spherical Torus Experiment (NSTX). Employing the state of the art extended-magnetohydrodynamic (MHD) code M3D-C1 it is demonstrated that non-ideal effects can significantly affect PB stability thresholds in NSTX discharges. In particular, robust resistive PB modes are found to exist well before the ideal PB stability threshold is met. These novel results can explain why ideal-MHD theory often does not accurately describe ELM onset in spherical torus configurations, and also present a valuable basis for the development of a predictive model for ELMs in low-aspect ratio tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Overview of the TCV tokamak program: Scientific progress and facility upgrades

The TCV tokamak is augmenting its unique historical capabilities (strong shaping, strong electron heating) with ion heating, additional electron heating compatible with high densities, and variable divertor geometry, in a multifaceted upgrade program designed to broaden its operational range without sacrificing its fundamental flexibility. The TCV program is rooted in a three-pronged approach aimed at ITER support, explorations towards DEMO, and fundamental research. A 1 MW, tangential neutral beam injector (NBI) was recently installed and promptly extended the TCV parameter range, with record ion temperatures and toroidal rotation velocities and measurable neutral-beam current drive. ITER-relevant scenario development has received particular attention, with strategies aimed at maximizing performance through optimized discharge trajectories to avoid MHD instabilities, such as peeling-ballooning and neoclassical tearing modes. Experiments on exhaust physics have focused particularly on detachment, a necessary step to a DEMO reactor, in a comprehensive set of conventional and advanced divertor concepts. The specific theoretical prediction of an enhanced radiation region between the two X-points in the low-field-side snowflake-minus configuration was experimentally confirmed. Fundamental investigations of the power decay length in the scrape-off layer (SOL) are progressing rapidly, again in widely varying configurations and in both D and He plasmas; in particular, the double decay length in L-mode limited plasmas was found to be replaced by a single length at high SOL resistivity. Experiments on disruption mitigation by massive gas injection and electron-cyclotron resonance heating (ECRH) have begun in earnest, in parallel with studies of runaway electron generation and control, in both stable and disruptive conditions; a quiescent runaway beam carrying the entire electrical current appears to develop in some cases. Developments in plasma control have benefited from progress in individual controller design and have evolved steadily towards controller integration, mostly within an environment supervised by a tokamak profile control simulator. TCV has demonstrated effective wall conditioning with ECRH in He in support of the preparations for JT-60SA operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A study of resistive peeling–ballooning modes across spherical tokamaks

We investigate how non-ideal-magnetohydrodynamics (MHD) effects, in particular plasma resistivity, impact the peeling–ballooning stability thresholds in spherical tokamaks. This analysis follows the discovery of resistive kink-peeling modes in ELMing National Spherical Torus Experiment (NSTX) discharges. In the present study we extend this modeling to ELMing pulses in the Mega Ampere Spherical Tokamak (MAST) and MAST—Upgrade (MAST-U), where we find a clear resistive scaling for peeling–ballooning modes. While in NSTX ideal-MHD predicts stability for ELMing discharges, in MAST-U we find that the plasma is slightly unstable to peeling–ballooning modes, but is fully stabilized once diamagnetic effects are considered in terms of a growth rate normalization. A resistive power law scaling is calculated for these modes on MAST-U, which lies in between that of tearing modes and resistive interchange modes. A comparison between M3D-C1 and NIMROD shows reasonable agreement for this scaling. Resistivity destabilizes the modes and the peeling–ballooning unstable domain is considerably expanded in both, MAST and MAST-U. In addition to the MAST/-U pulses we also analyze resistive PB stability in a NSTX-similarity discharge on DIII-D. While having a different aspect ratio from NSTX, this discharge uses NSTX-like shaping parameters, toroidal field and plasma current. By considering these discharges alongside NSTX cases, we identify conditions influencing the onset of resistive peeling–ballooning modes. Furthermore, our findings indicate that magnetic shear in the pedestal region is closely linked to the emergence of resistive edge modes.

MAST-U↗
DOE OSTI2024

Observation of a new pedestal stability regime in MAST Upgrade H-mode plasmas

Abstract The first pedestal stability and structure analysis on the new MAST Upgrade (MAST-U) spherical tokamak H-mode plasmas is presented. Our results indicate that MAST-U pedestals are close to the low toroidal mode number ( n ) peeling branch of the peeling-ballooning instability, in contrast with MAST H-mode pedestals which were deeply in the high- n ballooning branch. This offers the possibility of reaching the ELM-free quiescent H-mode (Burrell et al 2005 Plasma Phys. Control. Fusion 47 B37–B52) or high-performance super H-mode (Snyder et al 2015 Nucl. Fusion 55 083026; Snyder et al 2019 Nucl. Fusion 59 086017) regimes. In addition, the coupling between the peeling and ballooning branches is weak in MAST-U, suggesting that a path to very high pedestal pressure gradient at high density may exist with sufficient heating power. A possible explanation for the differences between MAST and MAST-U pedestal stability is given in terms of plasma shaping parameters, in particular squareness and elongation, as well as the pedestal top temperature and collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Overview of the ASDEX Upgrade results

After a 26-month vent ASDEX Upgrade (AUG) went back in operation with a newly designed upper W-divertor suitable for alternative divertor configurations (featuring in-vessel coils and cryo-pump). Parameter scans and an extensive set of measurements were obtained and their interpretation is ongoing. Prompted by the ITER wall change, dedicated experiments on non-boronized plasma startup were contrasted to that employing asymmetric and more symmetric boronizations. The asymmetric boronization proved to be as beneficial as the more symmetric one, which is in contrast to previous model calculations assuming perfect sticking of boron (measurements suggest sticking ≈ 0.3). In the startup phase also the impurity influxes at the outboard limiters were investigated contrasting the unboronized case featuring cold edges (low- Z radiation) to the boronized case, in which the lifetime of the boron layers could be estimated. Pedestal stability investigations revealed that the quasi-continuous exhaust (QCE) regime is obtained when ballooning modes are active in the vicinity of the separatrix and the global peeling-ballooning stability is high enough. Thus, at high enough shaping and high gas flux both can be achieved and QCE is a consequence. The closely related enhanced D-alpha (EDA) mode is not clearly distinguishable from QCE, e.g. the quasi-coherent mode characteristic for EDA also shows up in QCE. In QCE the impurity transport is behaving benign as could be measured for Ne with a novel analysis method making use of a comprehensive set of CXRS measurements. For high radiative fractions the regime of the X-point radiator (XPR) is accessible at AUG and the understanding of its access conditions and behaviour is further developed. Due to the localized radiative cooling at the X-point the XPR can be well diagnosed and thus controlled. For negative triangularity shapes, further experiments at increased shaping resulted in strongly heated L-mode plasmas avoiding ELMs. Two integrated modelling approaches towards ITER suggested that core W-accumulation will be no issue for ITER and that the fusion yield in ITER may be Q = 12 (i.e. ITPA20-IL scaling is too pessimistic). Further, investigations of the ITER ramp-down in AUG provide insights into maintaining position control. Various aspects of shattered pellet injection were investigated in AUG and one of the results show that with increasing Ne fraction the radiation during the current quench increases and the current decay becomes faster.

ASDEX Upgrade↗
DOE OSTI2024

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The effect of triangularity on predicted pedestals for the CFETR

In this work, the dominant unstable magnetohydrodynamic mode and its stable region in the pedestal of the China Fusion Engineering Test Reactor are studied by numerical modeling with the peeling-ballooning theory over a wide range of triangularities (δ) and collisionalities (υ*). A new accessible stable region is found at δ<−0.1 for large βN,ped and υ*. This new stable region expands toward lower βN,ped and υ* with decreasing δ and is totally covered by the peeling unstable region with δ=−0.5 because of the increasing trapped particle fraction (ft,eff). The sensitivity of this new stable region to the kinetic ballooning mode constraint and elongation is studied. For negative and low δ<0.1, the boundary of the first stable region is determined from the ballooning mode. For δ>0.1, the peeling mode becomes dominant as the boundary approaches low s with low υ*, while the ballooning mode is still dominant at the boundary with high υ*. When δ increases beyond 0.46, the first stable region is expanded, and access to the second stable region of the ballooning mode opens up. The entire newly connected region of the first and second stable regions keeps expanding with further increases in δ until δ=0.6. Beyond this value, the ballooning mode becomes more unstable with increasing δ, while the peeling mode is approximately unchanged with increasing δ in this range. The change in the dominant mode and the stable region with increasing δ can be explained by the trade-off between the stabilization effect from the deeper poloidal magnetic well and destabilization due to the enlarged drive term.

Physics↗
DOE OSTI2020

Identification of a network of nonlinear interactions as a mechanism triggering the onset of edge localized modes

Edge localized modes (ELMs) remain a critical issue for the lifetime of in-vessel components of fusion reactors such as ITER. We study ELMs triggered when the pedestal parameters are below the peeling-ballooning limit. These ELMs occur while the pedestal recovers after an ELM crash. We identify a mechanism by which a non-linear perturbation involving many active triads—a network of non-linear interactions – leads to the triggering of these ELMs. Each active triad leads to the non-linear transfer of energy between three waves, such that the network of non-linear interactions lead to a transfer of energy between many waves. The triggering of this non-linear perturbation systematically coincides with a pedestal perturbation induced by the neutral beam injection beams and localized near the magnetic surface of safety factor q = 5. Among the non-linear interactions detected during the ELM by a bi-coherence analysis, the most active triad involves a slow frequency that is qualitatively consistent with a geodesic acoustic mode located at the same location, i.e., near the q = 5 surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Progress from ASDEX Upgrade experiments in preparing the physics basis of ITER operation and DEMO scenario development

An overview of recent results obtained at the tokamak ASDEX Upgrade (AUG) is given. A work flow for predictive profile modelling of AUG discharges was established which is able to reproduce experimental H-mode plasma profiles based on engineering parameters only. In the plasma center, theoretical predictions on plasma current redistribution by a dynamo effect were confirmed experimentally. For core transport, the stabilizing effect of fast ion distributions on turbulent transport is shown to be important to explain the core isotope effect and improves the description of hollow low-Z impurity profiles. The L–H power threshold of hydrogen plasmas is not affected by small helium admixtures and it increases continuously from the deuterium to the hydrogen level when the hydrogen concentration is raised from 0 to 100%. One focus of recent campaigns was the search for a fusion relevant integrated plasma scenario without large edge localised modes (ELMs). Results from six different ELM-free confinement regimes are compared with respect to reactor relevance: ELM suppression by magnetic perturbation coils could be attributed to toroidally asymmetric turbulent fluctuations in the vicinity of the separatrix. Stable improved confinement mode plasma phases with a detached inner divertor were obtained using a feedback control of the plasma β. The enhanced D α H-mode regime was extended to higher heating power by feedback controlled radiative cooling with argon. The quasi-coherent exhaust regime was developed into an integrated scenario at high heating power and energy confinement, with a detached divertor and without large ELMs. Small ELMs close to the separatrix lead to peeling-ballooning stability and quasi continuous power exhaust. Helium beam density fluctuation measurements confirm that transport close to the separatrix is important to achieve the different ELM-free regimes. Based on separatrix plasma parameters and interchange-drift-Alfvén turbulence, an analytic model was derived that reproduces the experimentally found important operational boundaries of the density limit and between L- and H-mode confinement. Feedback control for the X-point radiator (XPR) position was established as an important element for divertor detachment control. Stable and detached ELM-free phases with H-mode confinement quality were obtained when the XPR was moved 10 cm above the X-point. Investigations of the plasma in the future flexible snow-flake divertor of AUG by means of first SOLPS-ITER simulations with drifts activated predict beneficial detachment properties and the activation of an additional strike point by the drifts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Impact of shape on pedestal characteristics in the wide pedestal quiescent H-mode in the DIII-D tokamak

Experiments on the DIII-D tokamak have advanced the operational limits of wide pedestal quiescent H-mode (WPQH) plasmas towards increased ITER relevance by simultaneously demonstrating well-matched plasma shape and net zero injected torque. Wide pedestal QHmodes are a compelling candidate regime for a future power producing device because they maintain a stationary pedestal without ELMs via additional edge transport. The pedestal is wider than what would be predicted from kinetic ballooning mode physics due to enhanced edge transport generated by broadband turbulence, a limit cycle oscillation, or some combination thereof. Here, compared to the double null shape, the lower single null shape is observed to have a lower density, narrower pedestal width, larger density fluctuations over a broad range of wavenumber, and operates closer to the peeling-ballooning instability boundary calculated from the simple pedestal scaling, Δ ψΝ , which is still observed to be wider than the EPED prediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

LLNL FESP Theory Highlights: March 2024

Predictions of heat load widths λ q based on particle orbits alone are very pessimistic. This paper shows that pedestal peeling-ballooning (P-B) MHD turbulence broadens the stable SOL by the transport, or spreading, of fluctuation energy from the pedestal. λ q is seen to increase with Γ ε , the fluctuation energy density flux. We elucidate the fundamental physics of the spreading process. Γ ε increases with pressure fluctuation correlation length. P-B turbulence is seen to be especially effective at spreading, on account of its large effective mixing length. Spreading is shown to be a multiscale process, which is enhanced by the synergy of large and small-scale modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The physics basis to integrate an MHD stable, high-power hybrid scenario to a cool divertor for steady-state reactor operation

Abstract Coupling a high-performance core to a low heat flux divertor is a crucial step for ITER and a Fusion Pilot Plant or DEMO. Experiments in DIII-D recently expanded the steady-state hybrid scenario to high density and divertor impurity injection to study the feasibility of a radiating mantle solution. This work presents the physics basis for trade-offs between density, current drive and stability to tearing modes (TMs) at high β. EC power is crucial to tailor the plasma profiles into a passively stable state, and to eject impurities from the core. Off-axis EC depositions decrease the heating efficiency, but calculated electron heat transport coefficients show that this effect is partially mitigated by improved confinement inside the EC deposition. Additionally, the reduction in pressure is recovered by increasing the density. This favourable scaling of confinement with density was observed in high power plasmas for years, and this work provides a comprehensive explanation. ELITE predictions indicate that a path in peeling-ballooning stability opens up for certain conditions of density, power, q 95 and shaping, allowing the edge pressure to continue increasing without encountering a limit. In the core, calculated anomalous fast-ion diffusion coefficients are consistent with density fluctuation measurements in the toroidicity-induced Alfvén eigenmode range, showing that smaller fast-ion losses contribute to the enhanced confinement at high density. The edge integration study shows that divertor heat loads can be reduced with Ne and Ar injection, but this eventually triggers a cascade of n = 1, 2, 3 core TMs. We can now show that impurity radiation in the core is small and it is not the cause for the drop in confinement at high Ar and Ne injection rates. The overlap between the core TMs is consistent with the loss of pressure as estimated by the Belt model for the coupled rational surfaces. Optimization of these trade-offs has achieved plasmas with sustained H 98y2 = 1.7, f GW = 0.7 and ∼85% mantle radiation. The scenario and its variations at higher density and on- vs off-axis EC heating has been studied as a candidate for an integrated solution for several reactor designs, such as ITER, ARC, and the ARIES-ACT1 case, showing promising results in terms of fusion power and gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

How fluctuation intensity flux drives SOL expansion

Abstract Predictions of heat load widths λ q based on particle orbits alone are very pessimistic. This paper shows that pedestal peeling-ballooning (P-B) magnetohydrodynamic (MHD) turbulence broadens the stable scrape-off layer (SOL) by the transport, or spreading, of fluctuation energy from the pedestal. λ q is seen to increase with Γ ε , the fluctuation energy density flux. We elucidate the fundamental physics of the spreading process. Γ ε increases with pressure fluctuation correlation length. P-B turbulence is seen to be especially effective at spreading, on account of its large effective mixing length. Spreading is shown to be a multiscale process, which is enhanced by the synergy of large and small-scale modes. Pressure fluctuation skewness correlates well with the spreading flux–with the zero crossing of skewness and Γ ε spatially coincident–suggesting the role of coherent fluctuation structures and the presence of intermittency in λ q broadening. λ q ∼ B p − 1 scaling persists for the broadened SOL. We show that the spreading flux increases for increasing pedestal pressure gradient ∇ P 0 and for decreasing pedestal collisionality υ ped ∗ . This trend is due to the dominance of peeling modes for large ∇ P 0 and low υ ped ∗ . Ultimately, we see that a state of weak MHD turbulence, as for small ELMs, is very attractive for heat load management. Our findings have transformative implications for future fusion reactor designs and call for experimental investigations to validate the observed trends.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Development of a neural network model for peeling–ballooning stability analysis in the KSTAR tokamak pedestals

The neural network model, MISHKA-NN is developed to mitigate the computational burden associated with the linear ideal magnetohydrodynamic (MHD) stability analysis of the pedestal based on the peeling–ballooning (P–B) model. By utilizing both 1D plasma profiles (current density, pressure gradient, and safety factor) and 0D parameters (plasma geometry, total current, and toroidal mode number), the model predicts linear growth rate of edge-localized ideal MHD instability in a given equilibrium state. By enabling the prediction of each instability within a second, the model reduces the time required for plotting a pedestal P–B stability diagram (j - α diagram) from approximately 100 CPU hours to a few CPU minutes. Notably, even with the utilization of parametric pressure and current profiles and plasma boundary shapes for the training dataset, the model shows a satisfactory level of performance in benchmarking the j - α diagram for the reconstructed equilibrium from a KSTAR tokamak experiment. We anticipate the model to serve as a versatile alternative to 2D linear MHD stability codes, alleviating numerical costs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗