Search NASA⌕ Search

SEARCH · Search NASA

Results for “Electron Fluid Model”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Modeling the effect of MHD activity on runaway electron generation during SPARC disruptions

Magnetohydrodynamic (MHD) instabilities and runaway electrons (REs) interact in several ways, making it important to self-consistently model these interactions for accurate predictions of RE generation and the design of mitigation strategies, such as massive gas injection (MGI). Using M3D-C1 – an extended MHD code with a RE fluid model – we investigate the effects of 3-D nonlinear MHD activity, material injection, and 2-D axisymmetric vertical displacement events (VDEs) on RE evolution during disruptions on SPARC – a high-field, high-current tokamak designed to achieve a fusion gain Q > 1. Several cases, comprising different combinations of neon (Ne) and deuterium (D 2 ) injection, are considered. Here, our results demonstrate key effects that arise from the self-consistent RE + MHD coupling, such as an initial increase in RE generation due to MHD instability growth, decreased saturation energies of the m/n = 1/1 mode driving sawteeth-like activity, RE losses in stochastic magnetic fields, and subsequent RE confinement and plateau formation due to re-healing of flux surfaces. Large RE plateaus (>5 MA) are obtained with Ne-only injection (2-5 × 10 21 atoms), while combined D 2 + Ne injection (2 × 10 21 Ne atoms; 1.8 × 10 22 D 2 molecules) produces a lower RE current (<2 MA). With D 2 + Ne injection, a post thermal quench "cold" VDE terminates the RE beam, preventing a steady plateau. These simulations couple REs, 3-D MHD instabilities, MGI, and axisymmetric VDEs for the first time in SPARC disruption simulations and represent a crucial step in understanding RE generation and mitigation in high-current devices like SPARC.

Datta, Rishabh [Massachusetts Inst. of Technology ↗
DOE OSTI2020

Revealing 3D Morphological and Chemical Evolution Mechanisms of Metals in Molten Salt by Multimodal Microscopy

Growing interest in molten salts as effective high-temperature heat-transfer fluids for sustainable energy systems drives a critical need to fundamentally understand the interactions between metals and molten salts. This work utilizes the multimodal microscopy methods of synchrotron X-ray nanotomography and electron microscopy to investigate the 3D morphological and chemical evolution of two-model systems, pure nickel metal and Ni-20Cr binary alloy, in a representative molten salt (KCl-MgCl2 50–50 mol %, 800 °C). In both systems, unexpected shell-like structures formed because of the presence of more noble tungsten, suggesting a potential route of using Ni–W alloys for enhanced molten-salt corrosion resistance. The binary alloy Ni-20Cr developed a bicontinuous porous structure, reassembling functional porous metals manufactured by dealloying. This work elucidates better mechanistic understanding of corrosion in molten salts, which can contribute to the design of more reliable alloys for molten salt applications including next-generation nuclear and solar power plants and opens the possibility of using molten salts to fabricate functional porous materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Estimates of global recycling coefficients for LTX- β discharges

We report the first observation of global recycling coefficient R near 0.5 in the Lithium Tokamak eXperiment-β (LTX-β), significantly below the minimum R previously reported in other devices. In a series of experiments with varied Li wall conditioning, estimates of the recycling coefficient have been made using a Lyman-α array and DEGAS2 modeling. A progressive reduction in Lyman-α emission with increased lithium and an increase in edge electron temperature are observed. It is also observed that with increasing Li coating thickness, the effective particle confinement time τp* is reduced and approaches TRANSP calculated energy confinement time (τE), with τp* near τE,TRANSP for the lowest recycling coefficients. Edge temperatures approaching core plasma temperatures, first reported in LTX, can now be directly connected to estimates of the recycling coefficient and qualitatively agree with previous UEDGE simulations. The particle flux to the limiting surfaces appears to be significantly reduced in comparison with fluid scrape-off layer (SOL) models, indicating that a large fraction of the SOL ions are mirror trapped. SOL collisionality drops more than an order of magnitude below the banana regime boundary, indicating the importance of kinetic effects. Full-f 1x2v gyrokinetic simulations of SOL field lines with the GKEYLL code indicate that the fraction of ions trapped along field lines increases as collisionality drops, as a result of increased lithium evaporation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Electron Beam Irradiation for Water Treatment of Per- and Polyfluoroalkyl Substances (PFAS)

Per- and polyfluoroalkyl substances (PFAS) are widely used but are now considered a water contamination risk. Fermi National Accelerator Laboratory’s IARC group has demonstrated that passing water through electron beam radiation can destroy PFAS. This project aims to develop a process to scale the system for bulk water treatment. The key is efficient radiation usage, ensuring that all water receives only the minimum dose. Software was developed for this purpose, consisting of a computational fluid dynamics model (CFD) in COMSOL, which calculates particle trajectories through the radiation area. A MATLAB script integrates the radiation dose of these particles, and statistical analysis is performed to evaluate the radiation utilization efficiency. These models are validated with a flow test where colored dye is injected and optically tracked. Radiation dose is measured by testing under an e-beam to measure the degradation of a PFAS analog.

Mueller, Scott [Northern Illinois U.]↗
DOE OSTI2026

Beam-water Interaction System for PFAS Treatment

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants that resist conventional water treatment due to the strength of their carbon fluorine bonds. Electron-beam irradiation has shown potential for PFAS degradation; however, its effectiveness is limited by shallow energy penetration, restricting treatment to the upper region of the water surface. This project addresses this limitation through the design of a compact mixing chamber and recirculation system to enhance surface renewal and improve treatment efficiency. The system, constrained to approximately 10 cm 20 cm 20 cm, was developed using computer-aided design, computational fluid dynamics (CFD), and MATLAB-based hydraulic analysis to evaluate flow behavior, head loss, and pump compatibility. A perforated-baffle geometry and controlled flow range of 1 3 GPM were implemented to promote uniform mixing while maintaining a stable free surface for electron-beam exposure. Modeling results indicate improved flow distribution and reduced dead zones, supporting more consistent interaction between contaminated water and the beam. The final design enhances hydrodynamic performance within system constraints, contributing to scalable electron-beam water treatment technologies for improved PFAS degradation.

Ramos, Gerardo [Northern Illinois U.]↗
DOE OSTI2021

Electric-Drive Vehicle Power Electronics Thermal Management: Current Status, Challenges, and Future Directions

Effective thermal management of traction-drive power electronics is critical to the advancement of electric-drive vehicles and is necessary for increasing power density and improving reliability. Replacing traditional silicon devices with more efficient, higher temperature, higher voltage, and higher frequency wide-bandgap (WBG) devices will enable increased power density but will result in higher device heat fluxes. Compact packaging of high-temperature WBG devices near low-temperature-rated components creates thermal management challenges that need to be addressed for future power-dense systems. This paper summarizes the thermal performance of on-road automotive power electronics thermal management systems and provides thermal performance and pumping-power metrics for select vehicles. Thermal analyses reveal that the package/conduction resistance dominates the total thermal resistance (for existing automotive systems). We model advanced packaging concepts and compare the results with existing packaging designs to quantify their thermal performance enhancements. Double-side-cooled configurations that do not use thermal interface materials are package concepts predicted to provide a low junction-to-fluid thermal resistance (compared to current packages). Dielectric-fluid-cooled concepts enable a redesign of the package to reduce the package resistance, can be implemented in single- and two-phase cooling approaches, and allow for cooling of passive components (e.g., capacitors) and bus bars.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Final Report for Center for Tokamak Transient Simulations at USU

Providing plasma fluid codes like NIMROD with continuum drift kinetic (CDK) physics that is quantitatively valid and computationally feasible throughout the spatial domain is difficult. Work at Utah State University (USU), in collaboration with the Center for Tokamak Transient Simulations (CTTS), focused on applying CDK closures in disruption-related calculations. Three examples where kinetic physics is paramount are (1) the electron stress tensor closure in Ohms law for accurately describing neoclassical tearing mode (NTM) evolution, (2) runaway electron (RE) density (nRE) and current (jRE) moments in NIMROD’s extended MHD model for self-consistent evolution of RE populations during disruptions and, (3) energetic ion effects on a myriad of MHD instabilities. While NTM simulations and continuum and PIC approaches to energetic ions in NIMROD have been a major goals of USU’s closure work for several years, the development of self-consistent CDK RE capability in NIMROD was started and extended considerably during the CTTS effort. Some goals of CDK RE in NIMROD are to explore the effects of the 2D relativistic phase space in 4D simulations and compare with NIMROD’s fluid RE model. Four publications and two PhD theses came out of the USU CTTS effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Formation of turing patterns in strongly magnetized electric discharges

Abstract Pattern formation and self-organization in many biological and non-biological systems can be explained through Turing’s activator-inhibitor model. Here we show how this model can be employed to describe the formation of filamentary structures in a low-pressure electric discharge exposed to a strong magnetic field. Theoretical investigation reveals that the fluid equations describing a magnetized plasma can be rearranged to take the mathematical form of Turing’s activator-inhibitor model. Numerical simulations based on the equations derived from this approach could reproduce the various patterns observed in the experiments. Also, it is shown that a density imbalance between electrons and ions exists in the bulk of the magnetized plasma that generates an electric field structure transverse to the applied magnetic field. This electric field is responsible for the stability of the filamentary patterns in the magnetized plasma over time scales much longer than the characteristic time scales of the electric discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Progress and innovations in the TCV tokamak research programme

Research on the Tokamak à Configuration Variable addresses a wide range of key questions relevant to ITER and future fusion power plants. Over the past two years, highly productive experimental campaigns have led to major advances across several areas: the ITER baseline scenario; pedestal properties in low-collisionality, peeling-limited conditions; and the development of high- β N , non-inductive regimes. Alternative high-confinement scenarios have likewise received significant attention, with remarkable progress in quasi-continuous exhaust operation, X-point radiator plasmas, and negative triangularity configurations. Substantial achievements were also made in the mitigation or benign termination of runaway electron beams, in elucidating fast-ion loss mechanisms, and in improving exhaust behaviour in both conventional and alternative divertor geometries. These experimental results have been strongly supported by advances in modelling and their direct application to the experiment, ranging from gyrokinetic simulations of core and pedestal turbulence to fluid-based studies of scrape-off layer and divertor physics in diverse geometries. Plasma control has taken on an increasingly important role, with model-based and data-driven approaches now closely intertwined with physics studies. This article provides a overview of these recent activities, together with a brief outlook on forthcoming upgrades and next steps.

EPFL↗
DOE OSTI2024

Do not forget the electrons: Extending moderately-sized nuclear networks for multidimensional hydrodynamic codes

Context.Nuclear networks are widely used coupled with hydrodynamical simulations of explosive scenarios to account for the change of nuclear species and energy generation rate due to nuclear reactions. In this way, there is a feedback mechanism between the hydrodynamical state and the nuclear processes. Unfortunately, the timescale of nuclear reactions is orders of magnitude smaller than the dynamical timescale that drives hydrodynamical simulations. Therefore, these nuclear networks are usually very small, reduced in most cases to a dozen elements, especially when simulations are carried out in more than one dimension. Aims.We present here an extended nuclear network, with 90 species, designed for being coupled with hydrodynamic simulations, which includes neutrons, protons, electrons, positrons, and the corresponding neutrino and anti-neutrino emission. This network is also coupled with temperature, making it extremely robust and, together with its size, unique of its kind. The inclusion of electron captures on free protons makes the network very appropriate for multidimensional studies of Type Ia supernova explosions, especially when the exploding object is a massive white dwarf. Methods.We perform several tests that are relevant to simulate explosive scenarios, such as Type Ia supernovae and core-collapse supernovae. We compare the results of the 90 nuclei network with a standardα-chain network with 14 elements to evaluate the differences in the energy generation rate. We also evaluate the relevance of including the electrons in the network in terms of generated yields and how it affects the pressure of a degenerate fluid such as that of white dwarfs. The results obtained with the 90-nuclei network have been verified with a much larger 2000-nuclei network built from REACLIB (WinNet), in terms of nuclear energy generation rate, pressure, and produced yields. Results.The results obtained with the proposed medium-sized network compare fairly well, to a few percent, with those computed withWinNetin scenarios reproducing the gross physical conditions of current Type Ia supernova explosion models. In those cases where the carbon and oxygen fuel ignites at high density, the high-temperature plateau typical of the nuclear statistical equilibrium regime is well defined and stable, allowing large integration time steps. We show that the inclusion of electron captures on free protons substantially improves the estimation of the electron fraction of the mixture. Therefore, the pressure is better determined than in networks where electron captures are excluded, which will ultimately lead to more reliable hydrodynamic models. Explosive combustion of helium at low density, occurring near the surface layer of a white dwarf, is also better described with the proposed network, which gives nuclear energy generation rates much closer toWinNetthan typical reduced alpha networks. Conclusions.A nuclear network withN= 90 species, including electrons, aimed at multidimensional calculations of supernova explosions is described and verified. The proposed network is suitable for the study of Type Ia supernova explosions because it provides better values of pressure and electron abundance than other existing networks with smaller or even a similar size but without including electron capture processes.

Astronomy & Astrophysics↗
DOE OSTI2023

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗
DOE OSTI2021

Impact of collisionality on turbulence in the edge of tokamak plasma using 3D global simulations

Collisionality is one of the key parameters in determining turbulent transport in the plasma edge, regulating phenomena such as 'shoulder formation', separation of scale lengths in the scrape-off layer (SOL), turbulence damping and zonal flow dynamics. Understanding its role is therefore of primary importance for future reactors like ITER. Getting reliable predictions and a better characterization of plasma flow properties when varying collisionality remains, however, a critical challenge for the simulations. This paper focuses on the impact of varying collisionality in a non-isothermal three-dimensional fluid model of the plasma edge. A high field side limited configuration encompassing open and closed magnetic field lines with parameters typical of a medium-sized tokamak is considered. The present model can consistently account for the variations of collisionality and its impact on both the parallel resistivity η ∥ and the ion and electron parallel thermal conductivities χ ∥e,i . Details on mean flow and turbulence properties are given. Changing collisionality leads to significant changes in the flow properties both on the mean and fluctuating quantities. In particular, lowering collisionality decreases the size of coherent structures, the fluctuation levels of turbulence, and steepens the density and temperature equilibrium profiles around the separatrix leading to a global reduction of the turbulent transport. The SOL width is observed to increase with collisionality, eventually resulting in the disappearance of the scale lengths separation between near and far SOL, consistently with previous experimental observations. At low collisionality, where the presence of narrow feature is well-established, a contribution of heat conduction increases up to compete with heat convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Density wavenumber spectrum measurements, synthetic diagnostic development, and tests of quasilinear turbulence modeling in the core of electron-heated DIII-D H-mode plasmas

Abstract Measurements of the turbulent density wavenumber spectrum, δ n ˆ e ( k ⊥ ) , using the Doppler Back-Scattering (DBS) diagnostic are reported from DIII-D H-mode plasmas with electron cyclotron heating as the only auxiliary heating method. These electron-heated plasmas have low collisionality, ν e ∗ < 1 , T e / T i > 1 , and zero injected torque—a regime expected to be relevant for future fusion devices. We probe density fluctuations in the core ( ρ ≈ 0.7) over a broad wavenumber range, 0.5 ⩽ k ⊥ ⩽ 16 cm −1 ( 0.1 ⩽ k ⊥ ρ s ⩽ 5 ), to characterize plasma instabilities and compare with theoretical predictions. We present a novel synthetic DBS diagnostic to relate the back-scattered power spectrum, P s ( k ⊥ ) —which is directly measured by DBS—to the underlying electron density fluctuation spectrum, δ n ˆ e ( k ⊥ ) . The synthetic DBS P s ( k ⊥ ) spectrum is calculated by combining the SCOTTY beam-tracing code with a model δ n ˆ e ( k ⊥ ) predicted either analytically or numerically. In this work we use the quasi-linear code Trapped Gyro-Landau Fluid (TGLF) to approximate the δ n ˆ e ( k ⊥ ) spectrum. We find that TGLF, using the experimental profiles, is capable of closely reproducing the DBS measurements. Both the DBS measurements and the TGLF-DBS synthetic diagnostic show a wavenumber spectrum with variable decay. The measurements show weak decay ( k −0.6 ) for k < 3.5 cm −1 , with k −2.6 at intermediate- k ( 3.5 ⩽ k ⩽ 8.5 cm −1 ), and rapid decay ( k −9.4 ) for k > 8.5 cm −1 . Scans of physics parameters using TGLF suggest that the normalized ∇ T e scale-length, R / L T e , is an important factor for distinguishing microturbulence regimes in these plasmas. A combination of DBS observations and TGLF simulations indicate that fluctuations remain peaked at ITG-scales (low k ) while R / L T e -driven TEM/ETG-type modes (intermediate/high k ) are marginally sub-dominant.

synthetic diagnostics↗
DOE OSTI2022

Quantifying Equiaxed vs Epitaxial Solidification in Laser Melting of CMSX-4 Single Crystal Superalloy

The competition between epitaxial vs. equiaxed solidification has been investigated in CMSX-4 single crystal superalloy during laser melting as practiced in additive manufacturing. Single-track laser scans were performed on a powder-free surface of directionally solidified CMSX-4 alloy with several combinations of laser power and scanning velocity. Electron backscattered diffraction (EBSD) mapping facilitated identification of new orientations, i.e., “stray grains” that nucleated within the fusion zone along with their area fraction and spatial distribution. Using high-fidelity computational fluid dynamics simulations, both the temperature and fluid velocity fields within the melt pool were estimated. This information was combined with a nucleation model to determine locations where nucleation has the highest probability to occur in melt pools. In conformance with general experience in metals additive manufacturing, the as-solidified microstructure of the laser-melted tracks is dominated by epitaxial grain growth; nevertheless, stray grains were evident in elongated melt pools. It was found that, though a higher laser scanning velocity and lower power are generally helpful in the reduction of stray grains, the combination of a stable keyhole and minimal fluid velocity further mitigates stray grains in laser single tracks.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Validation of the GFS model for gyrokinetic stability of NSTX pedestal data

This study presents a large database validation of the gyro fluid system (GFS) model for linear gyrokinetic stability for high-mode (H-mode) edge transport barrier conditions in the national spherical torus experiment (NSTX) tokamak. The database of linear stability calculations with the CGYRO gyrokinetic code was produced using plasma profile measurements from NSTX discharges to identify kinetic ballooning modes (KBM), trapped electron modes (TEM), and micro-tearing modes (MTM) that limit the pressure profile gradient in the H-mode barrier. A novel Bayesian optimization approach determines optimal resolution parameters for GFS specifically for spherical tokamak pedestal conditions. Our results demonstrate that GFS, with optimized resolution, can achieve accurate linear stability analysis in NSTX pedestal conditions for reduced resolution compared to CGYRO. GFS can accurately find the KBM, TEM, and MTM instability branches. Parametric analysis reveals that GFS accuracy in this extreme pedestal parameter range is degraded for low magnetic shear and near the separatrix conditions. These findings establish GFS as a fast linear eigenmode solver for spherical tokamak pedestal gyrokinetic stability and demonstrate a systematic methodology for determining the optimum resolution settings.

Yang, Minglei [Oak Ridge National Laboratory (ORNL↗
DOE OSTI2024

Runaway electron plateau current profile reconstruction from synchrotron imaging and Ar-II line polarization angle measurements in DIII-D

Abstract Current profile reconstructions are obtained for high current ( I p ≃ 550 kA) post-disruption runaway electron (RE) plateau plasmas in DIII-D. Two novel methods of measuring the RE current profile in high-current RE plateaus are introduced and compared: localization of the q = 2 rational surface using visible synchrotron emission (SE) imaging and the measurement of the polarization angle of line-integrated Ar-II line emission. The two methods are found to be consistent with each other within the data uncertainties. Different simulations of the RE current profile are compared with the measurements: the toroidal fluid RE model is found to best fit the data, within the measurement uncertainties. In addition to introducing two novel methods to measure the RE current profile and validating present simulation capabilities, this work demonstrates that instabilities can grow at q = 2 and q = 1 surfaces without necessarily causing a RE final loss instability. Numerical simulations are also presented to elucidate the role of these instabilities on synchrotron emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗