Intrinsic toroidal rotation driven by turbulent and neoclassical processes in tokamak plasmas from global gyrokinetic simulations
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Recently, it was shown that strong electron thermionic emission from material walls could result in the formation of an “inverse sheath”, which prevents the flow of cold ions to the wall. Such regimes look very favorably from the point of view of plasma-material interactions at the edge of magnetic fusion devices, where the problem of the erosion of plasma facing components under ion irradiation is one of the key issues for developing of future magnetic fusion reactors. However, it is not clear whether such regimes are compatible with edge plasma parameters and heat removal requirements in fusion reactors.To address the issue of practicality of the “inverse sheath” regime for edge tokamak plasma conditions, we perform a set of numerical simulations with 2D edge plasma transport code UEDGE for a DIII-D-like geometry and magnetic configuration.To describe both “standard” and “inverse sheath” conditions within the framework of the UEDGE code (which does not consider the sheath region per se), at the material surfaces, we apply effective boundary conditions which emulate both “standard” and“inverse sheath” regimes. We demonstrate that for the same input parameters, spatial distributions of edge plasma parameters corresponding to detached divertor and“inverse sheath” regimes are similar, with only a few minor differences. We discuss the compatibility of “inverse sheath” regimes with core plasma parameters.
The International Tokamak Physics Activity (ITPA) has developed a comprehensive overview of the results of its coordinated R&D activities in fusion physics and diagnostics implemented during the past two decades with the aim of developing an improved physics basis for the operation of tokamak burning plasma experiments. Here, an introduction is presented to key concepts in the physics of burning plasmas in tokamaks, together with an overview of the scope of the principal areas of physics R&D pursued by the ITPA Topical Physics Groups in which the critical issues, major areas of progress and most significant recent results are highlighted. This paper introduces the following collection of articles, which encompass detailed presentations of the progress achieved by the Topical Groups in preparing the physics basis for the operation of tokamak burning plasma experiments expected to come into operation in the 2030s.
Integration of divertor detachment with a high-performance (β N ~ 3, β p > 2, H 98 ~ 1.5) core plasma has been demonstrated in DIII-D high-β p (poloidal beta) plasmas associated with a sustained core internal transport barrier (ITB) and an H-mode edge transport barrier (ETB). Such good core-edge integration has been achieved for both neon and nitrogen seeding, for both favorable and unfavorable B-field directions, independently from the impurity puffing locations, though these variations play important roles on divertor characteristics. Compared to the standard H-mode plasmas, the high-β p plasma exhibits a much wider window of detachment compatible with high confinement core. Fully detached divertor plasmas with low plasma temperature (Te < 5 eV), low particle flux, and low heat flux across the entire divertor target plate were obtained by using nitrogen seeding. This detached high-β p plasma is compatible with a newly developed detachment control system which can help optimize the nitrogen gas flow rate. Several features, i.e., the high edge safety factor in the high-β p scenario, impurity injection, closed divertor and reduced heating power requirement due to the high confinement, facilitate the achievement of full divertor detachment at lower density. Instead of degrading global performance, the divertor detachment facilitates the access to an even stronger ITB at large radius with a relatively weak ETB through self-organized synergy between ITB and ETB, leading to sustained high confinement. The strengthening of the large-radius ITB compensates for the ETB degradation associated with divertor detachment. In addition, a weak ETB naturally has smaller edge localized modes (ELMs). In particular, with neon injection, a long-period no-ELM H-mode phase has been achieved simultaneously with high-performance core and partially detached divertor plasmas. Furthermore, these results demonstrate the possibility of integrating excellent core plasma performance with an effective divertor solution, an essential step toward steady-state operation of reactor-grade plasmas.
Not Available
This paper describes the extensive progress that has been made in the understanding of tokamak pedestal physics since the 2007 publication of ‘Progress in the ITER Physics Basis’ (Ikeda 2007 Nucl. Fusion 47 E01–S500). It serves as Chapter 3 of the 2025 Nuclear Fusion Special Issue titled ‘On the Path to Tokamak Burning Plasma Operation’ (Campbell et al 2025 Nucl. Fusion ). This review was compiled by the pedestal and edge physics (PEP) community affiliated with the International Tokamak Physics Activity organization. It attempts to collect in one place citations to the majority of published literature on the pedestal physics topics that will be most important for the operation of a future power producing burning plasma tokamak. These include citations to publications describing the physics of the pedestal plasmas in many operating tokamaks worldwide and the pedestal physics projections for several near-term future devices including ITER. Descriptions of experimental results, interpretive modeling and predictive extrapolations are integrated together and comprehensive references are provided. This review is organized around four primary technical sections, viz.: pedestal structure, edge localized mode (ELM) characteristics, ELM control and regimes without large ELMs. Key results from many of the references are described briefly and set into the tokamak burning plasma power plant context. In addition, different perspectives on pedestal physics topics that are currently under debate within the community are also described, to provide guidance on needs for future research. Finally, attempts are made to describe conclusions from all of this progress consistent with discussions by the pedestal physics community at this time. The goal of this review is to provide a useful reference document for pedestal physics researchers going forward toward operation of a burning tokamak fusion plasma.
In this study, a new simulation model for tokamak boundary plasma, SOLT3D, is implemented in the BOUT++ framework. The simulation model includes a set of dynamic equations describing collisional boundary plasma and neutral gas in the tokamak scrape-off layer and divertor region. The model is verified against standard linear plasma instabilities and available nonlinear results. For L-mode like conditions, SOLT3D reproduces characteristics of boundary plasma turbulent fluctuations that are similar to published experimental data, in terms of the amplitude and spatial dependence of the fluctuations. It also reproduces realistic plasma fluxes on material surfaces and expected Bohm-like effective radial transport. Plasma fluctuations inferred from the simulations imply inevitably a significant level of intrinsic 'noise' for axisymmetric tokamak plasma transport modeling, introducing errors on the order of unity. In particular, the toroidally averaged atomic rates below 5–10 eV are strongly modified by turbulent plasma fluctuations, which should significantly affect the standard axisymmetric modeling of the tokamak edge plasma and divertor.
Steering tokamak plasmas is commonly viewed as a way to avoid disruptions and runaway electrons. Plasma steering sounds as safe as driving to work but will be shown to more closely resemble driving at high speed through a dense fog on an icy road. The long time required to terminate an ITER discharge compared to time over which dangers can be foreseen is analogous to driving in a dense fog. Here, the difficulty of regaining plasma control if it is lost resembles driving on an icy road. Disruptions and runaways are associated with three issues—a solution to one tends to complicate the solution to the other two: loss of plasma position control, excessive heat deposition, and wall melting due to runaway electrons. All three risks must be addressed for ITER to achieve its mission and essentially eliminated before tokamak power plants can be deployed.
Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.
Abstract Establishment of an efficient central solenoid (CS) free tokamak plasma start-up method may lead to an economical fusion reactor. CS-free start-up using lower hybrid (LH) waves has been studied on the TST-2 spherical tokamak. Plasma current of about a quarter of CS-driven discharges has been obtained fully non-inductively using the outer-midplane and top LH launchers. Recently, an outer-off-midplane LH launcher was developed to achieve higher plasma current by optimizing for core absorption and minimal fast electron losses. Using the (outer-)off-midplane launcher, fully non-inductive plasma current start-up up to about 8 kA was achieved. Coupled ray-tracing and Fokker–Planck simulation was performed on equilibria reconstructed with an extended MHD model. It was found that the experimentally observed plasma current was in reasonable agreement with the numerical simulation. The simulation predicted appreciable orbit losses for the off-midplane launcher driven discharge at the present parameters, which was consistent with the experimentally observed x-ray radiation characteristics. The simulation showed that the current density was saturated for the present off-midplane launcher discharges and higher density and higher LH power was necessary to achieve higher plasma current.
We have constructed a nonhydrodynamic coronal model for calculating radiation from impurity atoms in a heated plasma. Some recent developments in the calculation of dielectronic recombination rate coefficients and collisional excitation rate coefficients are included. The model is applied to oxygen impurity radiation during the first few milliseconds of a TFR Tokamak plasma discharge, and good agreement with experimental results is obtained. Estimates of total line and continuum radiation from the oxygen impurity are given. It is shown that impurity radiation represents a considerable energy loss.
The paper presents computer simulation studies of burst injection of carbon and tungsten dust particles in DIII-D-like edge plasmas. The injection causes a large transient influx of the low- and high-Z impurities associated with the dust ablation in the plasmas. The dust transport and the effects of the ablated impurities on the edge plasma dynamics in a modern mid-size tokamak geometry are investigated for low- and high-power plasma discharge conditions. The core plasma contamination with dust-ablated impurities and the factors affecting it are evaluated.
During ICRF antenna operation, complex interactions between turbulent density profiles, nonlinear RF sheaths, and RF-induced convective transport are observed to alter plasma density in the tokamak edge [D’Ippolito et al., Nucl. Fusion 38, 1543 (1998)]. In this work, we explore the physics of such interactions via numerical modeling, using a nonlinear EM/plasma/sheath code (VSim) and profiles obtained from a fluid plasma turbulence code (Hermes) in a 3D slab domain containing biased side-wall limiters. RF-rectified sheath formation on antenna and limiter surfaces is observed as electromagnetic waves launched by the antenna are refracted through the turbulent density profile. On transport timescales, such sheath potentials have been shown to influence both the mean species density and its RMS fluctuation spectrum [Smithe et al., these proceedings]. On the faster RF timescales, we demonstrate that the converse is also true – regions of high plasma density near material surfaces give rise to the highest sheath potential amplitudes. When density is turbulent and spatially nonuniform, localized regions of high sheath potential (hotspots) may develop where high-density filaments intersect material surfaces. Such hotspots are of particular concern as sources of impurity sputtering, and we explore their behavior in response to changes both to the local plasma density and to antenna operating parameters and structure. Related results exploring the role of Faraday shields and/or enclosing structures in suppressing high sheath potentials for other devices (e.g. SPARC) will also be shown.
Abstract Future tokamak devices that aim to create conditions relevant to power plant operations must consider strategies for mitigating damage to plasma facing components in the divertor. One of the goals of MAST-U tokamak operations is to inform these considerations by researching advanced divertor configurations that aid stable plasma detachment. Machine design, scenario planning and detachment control would all greatly benefit from tools that enable rapid calculation of scenario-relevant quantities given some input parameters. This paper presents a method for generating large, simulated scrape-off layer data sets, which was applied to generate a data set of steady-state Hermes-3 simulations of the MAST-U tokamak. A machine learning model was constructed using a Bayesian approach to hyperparameter optimisation to predict diagnosable output quantities given control-relevant input features. The resulting best-performing model, which is based on a feedforward neural network, achieves high accuracy when predicting electron temperature at the divertor target and carbon impurity radiation front position and runs in around 1 ms in inference mode. Techniques for interpreting the predictions made by the model were applied, and a high-resolution parameter scan of upstream conditions was performed to demonstrate the utility of rapidly generating accurate predictions using the emulator. This work represents a step forward in the design of machine learning-driven emulators of tokamak exhaust simulation codes in operational modes relevant to divertor detachment control and plasma scenario design.
Abstract not provided.