This paper reviews current understanding of key physics elements that control the H-mode pedestal structure, which exists at the boundary of magnetically confined plasmas. The structure of interest is the width, height and gradient of temperature, density and pressure profiles in the pedestal. Emphasis is placed on understanding obtained from combined experimental, theoretical and simulation work and on results observed on multiple machines. Pedestal profiles are determined by the self-consistent interaction of sources, transport and magnetohydrodynamic limits. The heat source is primarily from heat deposited in the core and flowing to the pedestal. This source is computed from modeling of experimental data and is generally well understood. Neutrals at the periphery of the plasma provide the dominant particle source in current machines. This source has a complex spatial structure, is very difficult to measure and is poorly understood. For typical H-mode operation, the achievable pedestal pressure is limited by repetitive, transient magnetohydrodynamic instabilities. First principles models of peeling–ballooning modes are generally able to explain the observed limits. In some regimes, instability occurs below the predicted limits and these remain unexplained. Several mechanisms have been identified as plausible sources of heat transport. These include neoclassical processes for ion heat transport and several turbulent processes, driven by the steep pedestal gradients, as sources of electron and ion heat transport. Reduced models have successfully predicted the pedestal or density at the pedestal top. Firming up understanding of heat and particle transport remains a primary challenge for developing more complete predictive pedestal models.
Purpose-led Publishing is a coalition of three not-for-profit publishers in the field of physical sciences: AIP Publishing, the American Physical Society and IOP Publishing.
Together, as publishers that will always put purpose above profit, we have defined a set of industry standards that underpin high-quality, ethical scholarly communications.
We are proudly declaring that science is our only shareholder.
ISSN: 1361-6587
Plasma Physics and Controlled Fusion is a monthly publication dedicated to the dissemination of original results on all aspects of plasma physics and associated science and technology.
Open all abstracts, in this tab
R J Groebner and S Saarelma 2023 Plasma Phys. Control. Fusion 65 073001
T D Arber et al 2015 Plasma Phys. Control. Fusion 57 113001
Particle-in-cell (PIC) methods have a long history in the study of laser-plasma interactions. Early electromagnetic codes used the Yee staggered grid for field variables combined with a leapfrog EM-field update and the Boris algorithm for particle pushing. The general properties of such schemes are well documented. Modern PIC codes tend to add to these high-order shape functions for particles, Poisson preserving field updates, collisions, ionisation, a hybrid scheme for solid density and high-field QED effects. In addition to these physics packages, the increase in computing power now allows simulations with real mass ratios, full 3D dynamics and multi-speckle interaction. This paper presents a review of the core algorithms used in current laser-plasma specific PIC codes. Also reported are estimates of self-heating rates, convergence of collisional routines and test of ionisation models which are not readily available elsewhere. Having reviewed the status of PIC algorithms we present a summary of recent applications of such codes in laser-plasma physics, concentrating on SRS, short-pulse laser-solid interactions, fast-electron transport, and QED effects.
A Pavone et al 2023 Plasma Phys. Control. Fusion 65 053001
This article reviews applications of Bayesian inference and machine learning (ML) in nuclear fusion research. Current and next-generation nuclear fusion experiments require analysis and modelling efforts that integrate different models consistently and exploit information found across heterogeneous data sources in an efficient manner. Model-based Bayesian inference provides a framework well suited for the interpretation of observed data given physics and probabilistic assumptions, also for very complex systems, thanks to its rigorous and straightforward treatment of uncertainties and modelling hypothesis. On the other hand, ML, in particular neural networks and deep learning models, are based on black-box statistical models and allow the handling of large volumes of data and computation very efficiently. For this reason, approaches which make use of ML and Bayesian inference separately and also in conjunction are of particular interest for today's experiments and are the main topic of this review. This article also presents an approach where physics-based Bayesian inference and black-box ML play along, mitigating each other's drawbacks: the former is made more efficient, the latter more interpretable.
Nathan Mackey et al 2024 Plasma Phys. Control. Fusion 66 055018
In curved magnetic geometries, field-aligned regions of enhanced plasma pressure and density, termed 'blobs,' move as coherent filaments across the magnetic field lines. Coherent blobs account for a significant fraction of transport at the edges of magnetic fusion experiments and arise in naturally-occurring space plasmas. This work examines the dynamics of blobs with a fully kinetic electromagnetic particle-in-cell code and with a drift-reduced fluid code. In low-beta regimes with moderate blob speeds, good agreement is found in the maximum blob velocity between the two simulation schemes and simple analytical estimates. The fully kinetic code demonstrates that blob speeds saturate near the initial sound speed, which is a regime outside the validity of the reduced fluid model.
D Hachmeister et al 2024 Plasma Phys. Control. Fusion 66 055016
In tokamaks, radial transport is ballooning, meaning it is enhanced at the low-field side (LFS). This work investigates the effect of the magnetic configuration on the high-field side (HFS) scrape-off layer. Our experiments involved L-mode and H-mode discharges at ASDEX Upgrade, in which we scanned the magnetic configuration from a lower to an upper single-null shape, thus varying the location of the secondary separatrix. We show that the secondary separatrix determines the width of the HFS scrape-off layer, meaning that the density is much lower in the region that is magnetically disconnected from the LFS scrape-off layer, outside the secondary separatrix. Furthermore, we observe that the large density often seen in the HFS divertor drastically decreases as the separation between the primary and secondary separatrices falls below a particular value. This value is different for L-mode and H-mode plasmas and closely matches the power decay length measured at the LFS midplane. We also show how the HFS scrape-off layer density is smaller in an upper single-null than in a lower single-null, when the ionic grad-B drift points down. This difference is likely caused by reversing the drifts in the active divertor when switching the active X-point from the bottom to the top. We further observe that the neutral density in the lower divertor also correlates with the plasma shape and the high-density region in the HFS scrape-off layer. During the shape scans analyzed here, the HFS divertor remained partially detached throughout, with transitory reattachment modulated by ELM activity in H-mode. This work provides novel experimental data that can be leveraged to further the modeling capabilities and understanding of scrape-off layer physics in highly shaped plasmas.
M Giacomin et al 2024 Plasma Phys. Control. Fusion 66 055010
In this work, we present first-of-their-kind nonlinear local gyrokinetic (GK) simulations of electromagnetic turbulence at mid-radius in the burning plasma phase of the conceptual high-β, reactor-scale, tight-aspect-ratio tokamak Spherical Tokamak for Energy Production (STEP). A prior linear analysis in Kennedy et al (2023 Nucl. Fusion63 126061) reveals the presence of unstable hybrid kinetic ballooning modes (KBMs), where inclusion of the compressional magnetic field fluctuation, , is crucial, and subdominant microtearing modes (MTMs) are found at binormal scales approaching the ion-Larmor radius. Local nonlinear GK simulations on the selected surface in the central core region suggest that hybrid KBMs can drive large turbulent transport, and that there is negligible turbulent transport from subdominant MTMs when hybrid KBMs are artificially suppressed (through the omission of ). Nonlinear simulations that include perpendicular equilibrium flow shear can saturate at lower fluxes that are more consistent with the available sources in STEP. This analysis suggests that hybrid KBMs could play an important role in setting the turbulent transport in STEP, and possible mechanisms to mitigate turbulent transport are discussed. Increasing the safety factor or the pressure gradient strongly reduces turbulent transport from hybrid KBMs in the cases considered here. Challenges of simulating electromagnetic turbulence in this high-β regime are highlighted. In particular the observation of radially extended turbulent structures in the absence of equilibrium flow shear motivates future advanced global GK simulations that include .
E Tonello et al 2024 Plasma Phys. Control. Fusion 66 065006
L-mode negative triangularity (NT) operation is a promising alternative to the positive triangularity (PT) H-mode as a high-confinement edge localised mode-free operational regime. In this work, two TCV Ohmic L-mode core density ramps with opposite triangularity are investigated using SOLPS-ITER modelling. This numerical study aims to investigate the power exhaust differences between NT and PT focusing, in particular, on the geometrical effect of triangularity. To disentangle the latter from differences related to cross-field transport, anomalous diffusivities for particle () and energy () transport are fixed to the same values in PT and NT. The simulation results clearly show dissimilar transport and accumulation of neutral particles in the scrape-off layer for the two configurations. This gives rise to different ionization sources in the edge and divertor regions and produces differences in the poloidal and cross-field fluxes, ultimately leading to different power and particle divertor fluxes in the two configurations. Simulations recover the experimental feature of a hotter and attached outer target ( ) in the NT scenario compared to the PT counterpart.
Yinlong Guo et al 2024 Plasma Phys. Control. Fusion 66 055012
The discrete and stochastic nature of the processes in the strong-field quantum electrodynamics (SF-QED) regime distinguishes them from classical ones. An important approach to identifying the SF-QED features is through the interaction of extremely intense lasers with plasma. Here, we investigate the seeded QED cascades driven by two counter-propagating laser pulses in the background of residual gases in a vacuum chamber via numerical simulations. We focus on the statistical distributions of positron yields from repeated simulations under various conditions. By increasing the gas density, the positron yields become more deterministic. Although the distribution stems from both the quantum stochastic effects and the fluctuations of the environment, the quantum stochastic effects can be identified via the width of the distribution and the exceptional yields, both of which are higher than the quantum-averaged results. The proposed method provides a statistical approach to identifying the quantum stochastic signatures in SFQED processes using high-power lasers and residual gases in the vacuum chamber.
Clemente Angioni 2021 Plasma Phys. Control. Fusion 63 073001
In this paper, the theory of collisional and turbulent transport of impurities in tokamak plasmas is reviewed. The results are presented with the aim of providing at the same time a historical reconstruction of the scientific progress and a complete description of the present theoretical knowledge, with a hopefully sufficiently complete reference to the works which have been published in the field in the last decades. After a general introduction on the physics challenges offered by the problem of impurity transport and their relevance for practical nuclear fusion energy, the theory of collisional transport is presented. Here a specific section is also dedicated to the transport parallel to the magnetic field lines. A complete review of the transport mechanisms produced by turbulence follows. The corresponding comparisons between theoretical predictions and experimental observations are also presented, highlighting the influence that the validation activities had in motivating further theoretical investigations. The paper is completed by a section on the direct interactions between collisional and turbulent transport and by a final specific review dedicated to the progress in the theory–based modelling activities. In the writing of this review paper, the main goal has been to combine readability with completeness and scientific rigour, providing a comprehensive list of references for deeper documentation on specific aspects.
D Darian et al 2019 Plasma Phys. Control. Fusion 61 085025
The collected current by spherical and cylindrical Langmuir probes immersed in an unmagnetized and collisionless non-Maxwellian plasma at rest are theoretically studied, and analytical expressions for the currents of attracted and repelled plasma particles are presented. We consider Kappa, Cairns and the generalized Kappa–Cairns distributions as possible models for the velocity field in the plasma. The current–voltage characteristics curves are displayed and discussed. Furthermore, comparisons with the collected currents in Maxwellian plasmas are given. The results of Particle-in-Cell (PIC) simulations of spherical and cylindrical probes in non-Maxwellian plasmas are also presented, and compared with the theoretical expressions. The results for the collected currents by the Langmuir probes obtained by PIC simulations are in good agreement with the corresponding analytical expressions.
Open all abstracts, in this tab
Frank J Wessel et al 2024 Plasma Phys. Control. Fusion 66 065028
A charge- and current-neutralized, intense-ion beam (IIB) will propagate undeflected in a magnetized plasma by the , or diamagnetic, drift for a distance large compared to the beam-ion gyro-radius. This propagation occurs because of collective plasma effects when the beam-energy density is sufficient to sustain the polarization-electric field, or diamagnetic-screening currents. Our principal interest is the drift, with . IIBs are characterized by parameters in the range: 0.1–2 MeV ion energy, 1–100 MA m−2 ion-current density, and 17–350 kW average power produced on repetitively-pulsed systems. Multi-MW average-beam power is possible with appropriate modifications to the power supply and ion-source/accelerator. IIBs can be focused geometrically, and/or magnetically, to spot sizes of the order of a few cm's and their angular divergence is typically, ∼15 mRadians, suitable for long-range propagation in a vacuum beam-line. IIBs lose energy and momentum in the same manner as particles injected by a neutral-beam injector (NBIs), hence could be useful for fusion applications, including: heating, current drive, fueling, profile modifications, etc and such applications have yet to be thoroughly studied. Summarized here are the physical principles accounting for IIB propagation; ion-source designs used to produce the IIB; and pulsed-power methods for energizing the IIB accelerator. This technology base informs scalable metrics for the ion-source/accelerator (excluding the HV power supply and interconnects) used in a conceptual injector that provides the same ion energy and injected power (1 MeV, 17 MW) as NBIs used on ITER. A comparison indicates that the IIB injector would be much smaller in size, lower cost, and have much greater efficiency, while also providing for real-time modulation of the beam energy and intensity and the use of small-diameter injection ports that could minimize fuel contamination and magnetic-field leakage between the IIB injector and tokamak.
K Habib et al 2024 Plasma Phys. Control. Fusion 66 065027
A magnetized nonthermal electron–positron-ion (e-p-i) plasma is considered to study the propagation properties of ion-acoustic solitary and shock waves in the presence of trapped positrons and electrons for the first time. The Schamel-κ (kappa) distribution function that describes plasma nonthermality and particle trapping is assumed to consider electrons and positrons. The diffusive effect of ion plasma fluid, which is responsible for shock dynamics, is taken into account. A nonlinear Schamel-Korteweg–de Vries-Burgers' (SKdVB) equation is derived by employing the reductive perturbation approach, and the solitary and shock wave solutions of the SKdVB equation have also been derived for different limiting cases. It is found that only positive potential nonlinear structures (for both solitary and shock waves) are formed in the proposed plasma system. The condition for stable solitons in the absence of dissipation is analyzed, and the nature of arbitrary amplitude solitary waves (obtained via the Sagdeev potential approach) is discussed. It is found through theoretical and numerical investigation that different plasma compositional parameters (such as the trapping effect of electrons (βe) and positrons (βp), the obliquity effect (θ), electron-to-ion number density ratio (µe), the magnetic field effect (via Ω) and the viscous effect (via η)) have a significant influence on the dynamics of ion-acoustic solitary and shock waves. The theoretical and numerical investigations in this study may be helpful in describing the nature of localized structures in different plasma contexts, e.g. space and astrophysical plasmas and experimental plasmas where electron–positron-ion plasmas exist.
T Fu et al 2024 Plasma Phys. Control. Fusion 66 065026
In a quasi-axisymmetric stellarator, a significant bootstrap current will result in the generation of low-order rational surfaces and three-dimensional (3D) magnetic islands. In this paper, the influence of plasma density profiles on the equilibrium magnetic islands for the Chinese first quasi-axisymmetric stellarator (CFQS) is investigated using the HINT code. It is found that the flattening of the core plasma density profile leads to a significant suppression of magnetic islands. When the peaking factor of plasma density is 1.19, complete suppression of magnetic islands occurs while maintaining excellent integrity of the magnetic surface even with the volume-averaged plasma beta <β> increase up to 2%. On the other hand, during the transition of a plasma density profile from flat to hollow, there is a reversal in the core bootstrap current, resulting in reduction of rotational transform values to pass through the rational surface. Hence, formation of magnetic islands in the core region. Therefore, effective inhibition of CFQS's magnetic islands can be achieved by appropriately controlling density profiles through methods like gas injection.
J Rueda-Rueda et al 2024 Plasma Phys. Control. Fusion 66 065025
In this paper we demonstrate how the inversion, in energy and major radius (E, R) coordinates, of imaging neutral particle analyser (INPA) measurements can be used to obtain the fast-ion distribution. The INPA is most sensitive to passing ions with energies in the range (20–150) keV and pitches near 0.5 in the core and 0.7 near the plasma edge. Inversion of synthetic signals, via 0th-order Tikhonov and Elastic Net regularization, were performed to demonstrate the capability of recovering the ground truth fast-ion 2D phase-space distribution resolved in major radius and energy, even in the presence of moderate noise levels (10%). Finally, we apply our method to measure the 2D phase-space distribution in an MHD quiescent plasma at ASDEX Upgrade and find good agreement with the slowing down fast-ion distribution predicted by TRANSP.
Grant Rutherford et al 2024 Plasma Phys. Control. Fusion 66 065024
High field side lower hybrid current drive (LHCD) is one potential candidate for efficient non-inductive current drive in tokamak power plants, and the first test of this technology will occur on the DIII-D tokamak during the 2024 campaign. Previous LFS launch experiments operated in the multi-pass regime and relied on scrape-off layer interactions to close the spectral gap. In the DIII-D experiment, single-pass damping is achievable via an upshift in the parallel refractive index caused by mode converting twice (slow fast slow). This mode conversion affects the ray trajectories and can lead to enhanced upshift depending on where mode conversion occurs. Compared to multi-pass absorption experiments, the optimization of launched and plasma parameters can be counter-intuitive: increased density may increase efficiency and smaller tend to damp closer to the separatrix. A hard x-ray camera installed to measure the bremsstrahlung (50–250 keV) radiation from LHCD-generated fast electrons is capable of verifying the trends reporting in this paper through comparison to the ray-tracing/Fokker–Planck codes GENRAY/CQL3D.
Open all abstracts, in this tab
R J Groebner and S Saarelma 2023 Plasma Phys. Control. Fusion 65 073001
This paper reviews current understanding of key physics elements that control the H-mode pedestal structure, which exists at the boundary of magnetically confined plasmas. The structure of interest is the width, height and gradient of temperature, density and pressure profiles in the pedestal. Emphasis is placed on understanding obtained from combined experimental, theoretical and simulation work and on results observed on multiple machines. Pedestal profiles are determined by the self-consistent interaction of sources, transport and magnetohydrodynamic limits. The heat source is primarily from heat deposited in the core and flowing to the pedestal. This source is computed from modeling of experimental data and is generally well understood. Neutrals at the periphery of the plasma provide the dominant particle source in current machines. This source has a complex spatial structure, is very difficult to measure and is poorly understood. For typical H-mode operation, the achievable pedestal pressure is limited by repetitive, transient magnetohydrodynamic instabilities. First principles models of peeling–ballooning modes are generally able to explain the observed limits. In some regimes, instability occurs below the predicted limits and these remain unexplained. Several mechanisms have been identified as plausible sources of heat transport. These include neoclassical processes for ion heat transport and several turbulent processes, driven by the steep pedestal gradients, as sources of electron and ion heat transport. Reduced models have successfully predicted the pedestal or density at the pedestal top. Firming up understanding of heat and particle transport remains a primary challenge for developing more complete predictive pedestal models.
A Pavone et al 2023 Plasma Phys. Control. Fusion 65 053001
This article reviews applications of Bayesian inference and machine learning (ML) in nuclear fusion research. Current and next-generation nuclear fusion experiments require analysis and modelling efforts that integrate different models consistently and exploit information found across heterogeneous data sources in an efficient manner. Model-based Bayesian inference provides a framework well suited for the interpretation of observed data given physics and probabilistic assumptions, also for very complex systems, thanks to its rigorous and straightforward treatment of uncertainties and modelling hypothesis. On the other hand, ML, in particular neural networks and deep learning models, are based on black-box statistical models and allow the handling of large volumes of data and computation very efficiently. For this reason, approaches which make use of ML and Bayesian inference separately and also in conjunction are of particular interest for today's experiments and are the main topic of this review. This article also presents an approach where physics-based Bayesian inference and black-box ML play along, mitigating each other's drawbacks: the former is made more efficient, the latter more interpretable.
Annick Pouquet 2023 Plasma Phys. Control. Fusion 65 033002
Nonlinear phenomena and turbulence are central to our understanding and modeling of the dynamics of fluids and plasmas, and yet they still resist analytical resolution in many instances. However, progress has been made recently, displaying a richness of phenomena, which was somewhat unexpected a few years back, such as double constant-flux cascades of the same invariant for both large and small scales, or the presence of non-Gaussian wings in large-scale fields, for fluids and plasmas. Here, I will concentrate on the direct measurement of the magnitude of dissipation and the evaluation of intermittency in a turbulent plasma using exact laws stemming from invariance principles and involving cross-correlation tensors with both the velocity and the magnetic fields. I will illustrate these points through scaling laws, together with data analysis from existing experiments, observations and numerical simulations. Finally, I will also briefly explore the possible implications for the validity and use of several modeling strategies.
J Citrin and P Mantica 2023 Plasma Phys. Control. Fusion 65 033001
In recent years tokamak experiments and modelling have increasingly indicated that the interaction between suprathermal (fast) ions and thermal plasma can lead to a reduction of turbulence and an improvement of confinement. The regimes in which this stabilization occurs are relevant to burning plasmas, and their understanding will inform reactor scenario optimization. This review summarizes observations, simulations, theoretical understanding, and open questions on this emerging topic.
S M Kaye et al 2021 Plasma Phys. Control. Fusion 63 123001
In this paper, we review the thermal plasma confinement and transport properties observed and predicted in low aspect ratio tokamaks, or spherical tokamaks (STs), which can depart significantly from those observed at higher aspect ratio. In particular, thermal energy confinement scalings show a strong, near linear dependence of energy confinement time on toroidal magnetic field, while the dependence on plasma current is more modest, the opposite of what is seen at higher aspect ratio. STs have revealed a very strong improvement in normalized confinement with decreasing collisionality, much stronger than at higher aspect ratio, which bodes well for an ST-based fusion pilot plant should this trend continue at an even lower collisionality than has already been accessed. These differences arise because of fundamental differences in transport in STs due to the more extreme toroidicity (i.e. reduced region of bad curvature), and to the relatively larger shearing rates, both of which can suppress electrostatic drift wave instabilities at both ion and electron gyroradius scales. In addition, electromagnetic effects are much stronger in STs because they operate at high βT. Gyrokinetic (GK) studies, coupled with low- and high-k turbulence measurements, have shed light on the underlying physics controlling transport. At lower βT, both ion- and electron-scale electrostatic drift turbulence may be responsible for transport. At higher βT, microtearing, kinetic ballooning, and hybrid trapped electron/kinetic ballooning modes increasingly play a role, and they have a much stronger impact in the core of ST plasmas than at higher aspect ratio. Flow shear affects the balance between ion- and electron-scale modes. Non-linear GK simulations find regimes where the electron heat flux decreases with decreasing collisionality, consistent with the experimental global normalized confinement scaling. The ST is unique in that the relatively low toroidal magnetic field allows for localized measurements of electron-scale turbulence, and this coupled with turbulence measurements at ion-scales has facilitated detailed comparisons with GK simulations. These data have provided compelling evidence for the presence of ion temperature gradient and electron temperature gradient turbulence in some plasmas, and direct experimental support for the impact of experimental actuators like rotation shear, density gradient and magnetic shear on turbulence and transport.
Open all abstracts, in this tab
Ju et al
Intense attosecond light sources are required for investigating ultrafast electron dynamics, such as in molecules and atoms, etc. However, so far their achievable photon energy remains at the keV level, thus limiting their range of applications. For instance, probing intranuclear dynamics requires photon energies in the gamma-ray regime. Here we propose a scheme for generating intense attosecond gamma-ray pulses by irradiating a thin solid-density foil with tailored intense vortex laser. In the interaction, the laser driven foil electrons are differentially accelerated within a quarter-wave field and become highly bunched. Three-dimensional particle-in-cell simulations show that they efficiently (at ∼ 30% conversion efficiency) radiate ultrashort (∼ 400as), high-flux (> 1010 photons per pulse), and ultrabrilliant (up to 1027 photons s−1 mm−2 mrad−2 per 0.1% BW at 1 MeV photon energy) gamma ray. With an X-ray driving laser, the scheme can even produce isolated bright sub-attosecond gamma-ray pulses that are useful for time-resolved nuclear spectroscopy and many other areas.
Kim et al
In this research, we develop a data-driven disruption predictor based on Bayesian deep probabilistic learning, capable of predicting disruptions and modeling uncertainty in KSTAR. Unlike conventional neural networks within a frequentist approach, Bayesian neural networks can quantify the uncertainty associated with their predictions, thereby enhancing the precision of disruption prediction by mitigating false alarm rates through uncertainty thresholding. Leveraging 0D plasma parameters from EFIT and diagnostic data, a Temporal Convolutional Network adept at handling multi-time scale data was utilized. The proposed framework demonstrates proficiency in predicting disruptions, substantiating its effectiveness through successful applications to KSTAR experimental data.
Barberis et al
This study explores the influence of sawtooth oscillations on the velocity space distribution of fast ions in tokamak plasma discharges. The relevant Fokker-Planck equation for fast ions is solved analytically. Two distinct effects arising from the temperature drop associated with a sawtooth crash and their impact on the distribution function of fast ions are considered. The first effect involves the modulation of the fusion alpha particle source on the timescale of the sawtooth period, linked to the drop in fusion yield resulting from the sawtooth temperature relaxations. The second effect is tied to the increase of the slowing-down time during the sawtooth ramp, causing particles born later in the sawtooth cycle to experience reduced slowing down compared to those born right after the crash, creating an accumulation-like mechanism at higher energies.
In regimes where the sawtooth period is shorter than the fast ion slowing-down time, the combined influence of these effects gives rise to fast ion distribution functions that transiently exhibit positive slopes in velocity space.
Scotti et al
The detachment cliff is a bifurcative transition to partial detachment recently discovered at the DIII-D tokamak. This work presents a database analysis of target parameters in L-mode and H-mode discharges to search for a detachment cliff at ASDEX Upgrade (AUG). Most of the transitions from attached to partially detached divertor conditions observed in H- and L-mode discharges in AUG show bifurcative-like characteristics that are consistent with the properties of the detachment cliff if the Bx∇B drift is directed towards the active X-point. In the operational space of power and density, the bifurcative transitions identified during an L-mode discharge occur at injected power and density higher than a threshold value Ptot > 0.7 MW and ne > 1.6 x 1019 m-3, respectively). Furthermore, the temperatures at which the transitions start are found to be insensitive to the injected impurity, to the injected power and to the value of the upstream density. Finally, the study of the evolution of the target parameters, of the intensity of the Dα line and of specific manometers and bolometer lines of sights shows that the physical process underlying the detachment cliff and the self-sustained divertor oscillations might be the same.
Song et al
The loss of ICRF-heated NBI ions in EAST tokamak was numerically investigated by ORBIT code simulations. The effects of collisions and ripples on particle losses taken into account, and the distributions of fast ions generated by different beams in combination with ICRF heating were calculated using the TRANSP code. Results showed that ICRF waves altered the orbital distributions of beam ions, causing an increase in trapped ions and fast ion losses. Additionally, for co-current injected beams, perpendicular injection resulted in higher fast ion losses in synergistic heating than tangential injection. The study also found that the synergistic effect of collisions and ripples enhanced fast ion losses, which were highly localized and generated a maximum heat load of 0.165 MW/m2 on the first wall. However, conducting synergy heating experiments at high plasma currents and low effective ion charge numbers can significantly reduce the loss of fast ions.
Open all abstracts, in this tab
J Rueda-Rueda et al 2024 Plasma Phys. Control. Fusion 66 065025
In this paper we demonstrate how the inversion, in energy and major radius (E, R) coordinates, of imaging neutral particle analyser (INPA) measurements can be used to obtain the fast-ion distribution. The INPA is most sensitive to passing ions with energies in the range (20–150) keV and pitches near 0.5 in the core and 0.7 near the plasma edge. Inversion of synthetic signals, via 0th-order Tikhonov and Elastic Net regularization, were performed to demonstrate the capability of recovering the ground truth fast-ion 2D phase-space distribution resolved in major radius and energy, even in the presence of moderate noise levels (10%). Finally, we apply our method to measure the 2D phase-space distribution in an MHD quiescent plasma at ASDEX Upgrade and find good agreement with the slowing down fast-ion distribution predicted by TRANSP.
Grant Rutherford et al 2024 Plasma Phys. Control. Fusion 66 065024
High field side lower hybrid current drive (LHCD) is one potential candidate for efficient non-inductive current drive in tokamak power plants, and the first test of this technology will occur on the DIII-D tokamak during the 2024 campaign. Previous LFS launch experiments operated in the multi-pass regime and relied on scrape-off layer interactions to close the spectral gap. In the DIII-D experiment, single-pass damping is achievable via an upshift in the parallel refractive index caused by mode converting twice (slow fast slow). This mode conversion affects the ray trajectories and can lead to enhanced upshift depending on where mode conversion occurs. Compared to multi-pass absorption experiments, the optimization of launched and plasma parameters can be counter-intuitive: increased density may increase efficiency and smaller tend to damp closer to the separatrix. A hard x-ray camera installed to measure the bremsstrahlung (50–250 keV) radiation from LHCD-generated fast electrons is capable of verifying the trends reporting in this paper through comparison to the ray-tracing/Fokker–Planck codes GENRAY/CQL3D.
S Munaretto et al 2024 Plasma Phys. Control. Fusion 66 065023
Achieving edge localized modes (ELMs) suppression in spherical tokamaks by applying resonant magnetic perturbations (RMPs) has proven challenging. The poloidal spectrum of the applied RMP is a key parameter that has an impact on the capability to mitigate and eventually suppress ELMs. In this work the resistive magnetohydrodynamic code MARS-F (Liu et al 2000 Phys. Plasmas7 3681) is used to evaluate the possibility of directly measuring the plasma response in MAST-U, and particularly its variation as function of the applied poloidal spectrum, in order to guide the experimental validation of the predicted best RMP configuration for ELM suppression. Toroidal mode number n = 2 RMP is considered to minimize the presence of sidebands, and to avoid the deleterious core coupling of n = 1. Singular Value Decomposition is used to highlight linearly independent structures in the simulated magnetic 3D fields and how those structures can be measured at the wall where the magnetic sensors are located. Alternative ways to measure the multimodal plasma response and how they can be used to infer the best RMP configuration to achieve ELM suppression are also presented, including the plasma displacement and the 3D footprints at the divertor plates.
Q Xia et al 2024 Plasma Phys. Control. Fusion 66 065022
Tokamak edge turbulence is crucial for the cross-field transport of particles and energy away from the separatrix. A better understanding of what affects the turbulence helps to control the heat flux to the divertor targets and the wall. One potentially important factor is the ion particle source in the divertor, as the neutral pathways and the ionisation source distributions are different depending on the divertor geometry, e.g. vertical- and horizontal-target configurations. Numerically, how to represent the sources and mimic the effects on the SOL in the simulations is still an open question. In this paper, we use a 3D turbulence code STORM, based on drift-reduced Braginskii equations, to study the effects of the divertor particle source distribution on turbulence in a simplified 3D slab geometry. The results show that it requires a large amount of divertor particle source to be peaked near the separatrix to alter the heat flux deposited on the target in attached conditions. This large non-uniform particle source can locally enhance the turbulence in the divertor volume, which redistributes the energy flux to the target and reduces the maximum amplitude. Meanwhile, the plasma profiles evaluated at the outboard midplane, such as the amplitudes and fluctuations of the density and temperature, are marginally changed. Another consequence of our results is that the prediction of the temperature difference between the outboard midplane and the target would be underestimated, if the calculation only considers the conductive heat flux and ignores this enhanced cross-field transport in the divertor.
J R Harrison et al 2024 Plasma Phys. Control. Fusion 66 065019
The integration of good core and edge/pedestal confinement with strong dissipation of heat and particles in the divertors is a significant challenge for the development of fusion energy. Alternative divertor configurations offer potential advantages by broadening the operational space where a device can operate with detached divertors and acceptable power exhaust. First results from MAST Upgrade are presented from high confinement mode experiments with outer divertors in the Super-X divertor configuration, showing that the outer divertors naturally detach when the Super-X is formed with no discernible impact on the plasma core and pedestal. These initial findings confirm predicted benefits of the Super-X configuration in terms of facilitating scenario integration.
S Guizzo et al 2024 Plasma Phys. Control. Fusion 66 065018
Negative triangularity (NT) tokamak configurations may be more susceptible to magneto-hydrodynamic instability, posing challenges for recent reactor designs centered around their favorable properties, such as improved confinement and operation free of edge-localized modes. In this work, we assess the vertical stability of plasmas with NT shaping and develop potential reactor solutions. When coupled with a conformal wall, NT equilibria are confirmed to be less vertically stable than equivalent positive triangularity (PT) configurations. Unlike PT, their vertical stability is degraded at higher poloidal beta. Furthermore, improvements in vertical stability at low aspect ratio do not translate to the NT geometry. NT equilibria are stabilized in PT vacuum vessels due to the increased proximity of the plasma and the wall on the outboard side, but this scenario is found to be undesirable due to reduced vertical gaps which give less spatial margin for control recovery. Instead, we demonstrate that informed positioning of passively conducting plates can lead to improved vertical stability in NT configurations on par with stability metrics expected in PT scenarios. An optimal setup for passive plates in highly elongated NT devices is presented, where plates on the outboard side of the device reduce vertical instability growth rates to 16% of their baseline value. For lower target elongations, integration of passive stabilizers with divertor concepts can lead to significant improvements in vertical stability. Plates on the inboard side of the device are also uniquely enabled in NT geometries, providing opportunity for spatial separation of vertical stability coils and passive stabilizers.
Qinglai Qiu et al 2024 Plasma Phys. Control. Fusion 66 065021
In 2021, EAST was equipped with a full-ring divertor coil to facilitate research on the fish tail divertor concept. Initially, it was observed that the coil current had a negligible ability to sweep the strike point. Conversely, when the amplitude and frequency of the alternating current were marginally increased, there was a significant interruption to plasma control. This perturbation was attributed to the poloidal control field's limited response rate to the coil's fluctuations. To address this issue, novel control methodologies were devised to ensure stable and effective sweeping of the strike point using the divertor coil. The devised methods are twofold: For high-frequency strike point control, a low-pass filter decoupling technique based on ISOFLUX control strategy enabled achieving a sweeping frequency of 100 Hz. This strategy allowed for consistent plasma management without compromising average stored energy or density regulation. Resulting from this proficient manipulation of the strike point, a reduction in the peak temperature of the divertor plate was observed. For low-frequency sweeping, a static multi-input multi-output decoupling approach was developed, facilitating concurrent sweeping of both the outer and inner strike points.
A P L Robinson 2024 Plasma Phys. Control. Fusion 66 065020
It is argued that fusion chain reactions in the D-D system is feasible with supra-thermal deuterons in the MeV regime, with new generations of deuterons being generated either via neutron–deuteron or proton–deuteron collisions. The propagation of supra-thermal deuterons in an infinite, hot, dense deuterium target was studied using a Monte Carlo method that includes multiple nuclear reactions, electron and ion stopping, along with neutron and proton knock-ons. Over a wide range of densities we observed significant, albeit sub-critical chain reactions in the multi-keV temperature regime. At very high densities (over 1000 gcm−3) and temperatures (over 40 keV) we observed chain reactions that reached criticality. These results suggest that there is a case to re-assess the potential of inertial confinement fusion based on deuterium-heavy targets.
Tommaso Barberis and Franco Porcelli 2024 Plasma Phys. Control. Fusion
This study explores the influence of sawtooth oscillations on the velocity space distribution of fast ions in tokamak plasma discharges. The relevant Fokker-Planck equation for fast ions is solved analytically. Two distinct effects arising from the temperature drop associated with a sawtooth crash and their impact on the distribution function of fast ions are considered. The first effect involves the modulation of the fusion alpha particle source on the timescale of the sawtooth period, linked to the drop in fusion yield resulting from the sawtooth temperature relaxations. The second effect is tied to the increase of the slowing-down time during the sawtooth ramp, causing particles born later in the sawtooth cycle to experience reduced slowing down compared to those born right after the crash, creating an accumulation-like mechanism at higher energies.
In regimes where the sawtooth period is shorter than the fast ion slowing-down time, the combined influence of these effects gives rise to fast ion distribution functions that transiently exhibit positive slopes in velocity space.
Luca Scotti et al 2024 Plasma Phys. Control. Fusion
The detachment cliff is a bifurcative transition to partial detachment recently discovered at the DIII-D tokamak. This work presents a database analysis of target parameters in L-mode and H-mode discharges to search for a detachment cliff at ASDEX Upgrade (AUG). Most of the transitions from attached to partially detached divertor conditions observed in H- and L-mode discharges in AUG show bifurcative-like characteristics that are consistent with the properties of the detachment cliff if the Bx∇B drift is directed towards the active X-point. In the operational space of power and density, the bifurcative transitions identified during an L-mode discharge occur at injected power and density higher than a threshold value Ptot > 0.7 MW and ne > 1.6 x 1019 m-3, respectively). Furthermore, the temperatures at which the transitions start are found to be insensitive to the injected impurity, to the injected power and to the value of the upstream density. Finally, the study of the evolution of the target parameters, of the intensity of the Dα line and of specific manometers and bolometer lines of sights shows that the physical process underlying the detachment cliff and the self-sustained divertor oscillations might be the same.