The heating effect of electromagnetic waves in ion cyclotron range of frequencies (ICRFs) in magnetic confinement fusion device is different in different plasma conditions. In order to evaluate the ICRF heating effect in different plasma conditions, we conducted a series of experiments and corresponding TRANSP simulations on the EAST tokamak. Both simulation and experimental results show that the effect of ICRF heating is poor at low core electron density. The decrease in electron density changes the left-handed electric field near the resonant layer, resulting in a significant decrease in the power absorbed by the hydrogen fundamental resonance. However, quite a few experiments must be performed in plasma conditions with low electron density. It is necessary to study how to make ICRF heating best in low electron density plasma. Through a series of simulation scans of the parallel refractive index (n//) of the ICRF antenna, it is concluded that the change of the ICRF antenna n// will lead to the change of the left-handed electric field, which will change the fundamental absorption of ICRF power by the hydrogen minority ions. Fully considering the coupling of ion cyclotron wave at the tokamak boundary and the absorption in the plasma core, optimizing the ICRF antenna structure and selecting appropriate parameters such as parallel refractive index, minority ion concentration, resonance layer position, plasma current and core electron temperature can ensure better heating effect in the ICRF heating experiments in the future EAST upgrade. These results have important implications for the enhancement of the auxiliary heating effect of EAST and other tokamaks.

This paper presents a Hammir tandem mirror confinement performance analysis based on Realta Fusion’s first-of-a-kind model for axisymmetric magnetic mirror fusion performance. This model uses an integrated end plug simulation model including, heating, equilibrium and transport combined with a new formulation of the plasma operation contours (POPCONs) technique for the tandem mirror central cell. Using this model in concert with machine learning optimization techniques, it is shown that an end plug utilizing high temperature superconducting magnets and modern neutral beams enables a classical tandem mirror pilot plant producing a fusion gain Q > 5. The approach here represents an important advance in tandem mirror design. The high-fidelity end plug model enables calculations of heating and transport in the highly non-Maxwellian end plug to be made more accurately. The detailed end plug modelling performed in this work has highlighted the importance of classical radial transport and neutral beam absorption efficiency on end plug viability. The central cell POPCON technique allows consideration of a wide range of parameters in the relatively simple near-Maxwellian central cell, facilitating the selection of more optimal central cell plasmas. These advances make it possible to find more conservative classical tandem mirror fusion pilot plant operating points with lower temperatures, neutral beam energies and end plug performance requirements than designs in the literature. Despite being more conservative, it is shown that these operating points have sufficient confinement performance to serve as the basis of a viable fusion pilot plant provided that they can be stabilized against magnetohydrodynamic and trapped particle modes.

Plasma-terminating disruptions represent a critical outstanding issue for reactor-relevant tokamaks. ITER will use shattered pellet injection (SPI) as its disruption mitigation system to reduce heat loads, vessel forces and to suppress the formation of runaway electrons. In this paper we demonstrate that reduced kinetic modelling of SPI is capable of capturing the major experimental trends in ASDEX Upgrade SPI experiments, such as dependence of the radiated energy fraction on neon content, or the current quench dynamics. Simulations are also consistent with the experimental observation of no runaway electron generation with neon and mixed deuterium–neon pellet composition. We also show that statistical variations in the fragmentation process only have a notable impact on the disruption dynamics at intermediate neon doping, as was observed in experiments.

Electron-only reconnection (E-REC) is a process recently observed in the Earth’s magnetosheath, where magnetic reconnection occurs at electron kinetic scales, and ions do not couple to the reconnection process. Electron-only reconnection is likely to have a significant impact on the energy conversion and dissipation of turbulence cascades at kinetic scales in some settings. This paper investigates E-REC under different intensities of strong guide fields (the ratio between the guide field and the in-plane asymptotic field strength is 5, 10 and 20, respectively) via two-dimensional fully kinetic particle-in-cell simulations, focusing on electron heating. The simulations are initialized with a force-free current sheet equilibrium under various intensities of strong guide fields. Similarly to previous experimental studies, electron temperature anisotropy along separatrices is observed, which is found to be mainly caused by the variations of parallel temperature. Both regions of anisotropy and parallel temperature increase/decrease along separatrices become thinner with increasing guide fields. Besides, we find a transition from a quadrupolar to a hexapolar (six-polar) to an octopolar (eight-polar) structure in temperature anisotropy and parallel temperature as the guide field intensifies. Non-Maxwellian electron velocity distribution functions (EVDFs) at different locations in the three simulations are observed. Our results show that parallel electron velocity varies notably with different guide field intensities and finite parallel electron heat flux density is observed. The three simulations exhibit features of the Chew–Goldberger–Low theory, with the level of consistency increasing as the guide field strength increases. This explains the electron parallel temperature variations and the shape of the EVDFs observed along the separatrices. This work may provide insights into the understanding of electron heating and parallel heat flux density in E-REC observed in the turbulent magnetosheath.

The shear Alfvén wave (SAW) continuum plays a critical role in the stability of energetic particle-driven Alfvén eigenmodes. We develop a theoretical framework to analyze the SAW continuum in three-dimensional (3-D) quasisymmetric magnetic fields, focusing on its implications for stellarator design. By employing a near-axis model and degenerate perturbation theory, the continuum equation is solved, highlighting unique features in 3-D configurations, such as the interactions between spectral gaps. Numerical examples validate the theory, demonstrating the impact of flux-surface shaping and quasisymmetric field properties on continuum structure. The results provide insights into optimizing stellarator configurations to minimize resonance-driven losses of energetic particles. This work establishes a basis for incorporating Alfvénic stability considerations into the stellarator design process, demonstrated through optimization of a quasihelical configuration to avoid high-frequency spectral gaps.

Eddy currents play a significant role in the evolution of tokamak plasmas and must therefore be correctly taken into account in time-dependent simulations. In this paper, a computational method for solving the evolution of tokamak plasma considering eddy currents utilising VMEC (Hirshman & Whitson, Phys. Fluids, vol. 26, 1983, pp. 3553–3568), a commonly used static magnetohydrodynamic equilibrium solver, is proposed. This method is convenient since it does not modify the equilibrium solver internally and achieves convergence calculation through external processing. By allowing the components of the magnetic field to be treated separately, this method provides convergence for cases with displacements in arbitrary directions, which has been difficult to achieve with the previous methods.

Magnetic reconnection, a fundamental plasma process, is pivotal in understanding energy conversion and particle acceleration in astrophysical systems. While extensively studied in two-dimensional (2-D) configurations, the dynamics of reconnection in three-dimensional (3-D) systems remains under-explored. In this work, we extend the classical tearing mode instability to three dimensions by introducing a modulation along the otherwise uniform direction in a 2-D equilibrium, given by $g(y)$, mimicking a flux-tube-like configuration. We perform linear stability analysis (both analytically and numerically) and direct numerical simulations to investigate the effects of three-dimensionality. Remarkably, we find that a tearing-like instability arises in three dimensions as well, even without the presence of guide fields. Further, our findings reveal that the 3-D tearing instability exhibits reduced growth rates compared with two dimensions by a factor of $\int g(y)^{1/2} {\rm d}y\,/\int {\rm d}y$, with the dispersion relation maintaining similar scaling characteristics. We show that the modulation introduces spatially varying resistive layer properties, which influence the reconnection dynamics.

For small-shear helical-axis stellarators, linear ideal-magnetohydrodynamic (MHD) stability calculations and full-torus, nonlinear, electromagnetic gyrokinetic (GK) simulations (the latter with this unprecedented combination of objectives in stellarator GKs) in their linear phase are shown to yield well agreeing spatio-temporal structures of unstable, globally extended perturbations. Likewise, good agreement is found for their dependence on the plasma pressure and the vacuum-field magnetic well in plasma equilibria with identical gradient lengths of the temperature and density profiles. In the nonlinear phase, these perturbations with MHD signatures entail deformations of the magnetic surfaces, growing magnetic islands which rotate in the electron diamagnetic direction and, eventually, lead to ergodisation of a larger part of the magnetic surfaces.

The kinetic stability of collisionless, sloshing beam-ion ($45^\circ$ pitch angle) plasma is studied in a three-dimensional (3-D) simple magnetic mirror, mimicking the Wisconsin high-temperature superconductor axisymmetric mirror experiment. The collisional Fokker–Planck code CQL3D-m provides a slowing-down beam-ion distribution to initialize the kinetic-ion/fluid-electron code Hybrid-VPIC, which then simulates free plasma decay without external heating or fuelling. Over $1$–$10\;\mathrm{\unicode{x03BC} s}$, drift-cyclotron loss-cone (DCLC) modes grow and saturate in amplitude. The DCLC scatters ions to a marginally stable distribution with gas-dynamic rather than classical-mirror confinement. Sloshing ions can trap cool (low-energy) ions in an electrostatic potential well to stabilize DCLC, but DCLC itself does not scatter sloshing beam-ions into the said well. Instead, cool ions must come from external sources such as charge-exchange collisions with a low-density neutral population. Manually adding cool $\mathord {\sim } 1\;\mathrm{keV}$ ions improves beam-ion confinement several-fold in Hybrid-VPIC simulations, which qualitatively corroborates prior measurements from real mirror devices with sloshing ions.

One of the critical challenges in future high-current tokamaks is the avoidance of runaway electrons during disruptions. Here, we investigate disruptions mitigated with combined deuterium and noble gas injection in SPARC. We use multi-objective Bayesian optimisation of the densities of the injected material, taking into account limits on the maximum runaway current, the transported fraction of the heat loss and the current quench time. The simulations are conducted using the numerical framework Dream (disruption runaway electron analysis model). We show that during deuterium operation, runaway generation can be avoided with material injection, even when we account for runaway electron generation from deuterium–deuterium induced Compton scattering. However, when including the latter, the region in the injected-material-density space corresponding to successful mitigation is reduced. During deuterium–tritium operation, acceptable levels of runaway current and transported heat losses are only obtainable at the highest levels of achievable injected deuterium densities. Furthermore, disruption mitigation is found to be more favourable when combining deuterium with neon, compared with deuterium combined with helium or argon.

The propagation and absorption of the slow waves in the plasma of the Joint European Torus (JET) tokamak have been investigated by ray tracing. The study aims to obtain a qualitative notion of the penetration into the plasma and absorption of the slow wave excited by the A2 ITER-like antenna. The slow waves are radiated by antennas in the ion cyclotron resonance frequency inverted minority heating or mode conversion heating regimes. It has been discovered that the rays propagate in the toroidal direction over a significant distance, up to 6 × 103 cm, from the antenna. Spreading in the peripheral plasma, mainly between the separatrix and the wall, they slowly shift in the poloidal direction and can reach the divertor region. The change in equilibrium of the JET tokamak has a strong influence on both the propagation and absorption of slow waves. Absorption of the slow waves is caused by ion–electron collisions and Landau damping. In the minority heating regimes, the slow waves are strongly damped in the cyclotron resonance of minority ions even at very low minority density.

Electron energisation by magnetic reconnection has historically been studied in the Lagrangian guiding-centre framework. Insights from such studies include that Fermi acceleration in magnetic islands can accelerate electrons to high energies. An alternative Eulerian fluid formulation of electron energisation was recently used to study electron energisation during magnetic reconnection in the absence of magnetic islands. Here, we use particle-in-cell simulations to compare the Eulerian and Lagrangian models of electron energisation in a set-up where reconnection leads to magnetic island formation. We find the largest energisation at the edges of magnetic islands. There, energisation related to the diamagnetic drift dominates in the Eulerian model, while the Fermi related term dominates in the Lagrangian model. The models predict significantly different energisation rates locally. A better agreement is found after integrating over the simulation domain. We show that strong magnetic curvature can break the magnetic moment conservation assumed by the Lagrangian model, leading to erroneous results. The Eulerian fluid model is a complete fluid description and accurately models bulk energisation. However, local measurements of its constituent energisation terms need not reflect locations where plasma is heated or accelerated. The Lagrangian guiding centre model can accurately describe the energisation of particles, but it cannot describe the evolution of the fluid energy. We conclude that while both models can be valid, they describe two fundamentally different quantities, and care should be taken when choosing which model to use.

Optimal-mode theory (Landreman et al. 2015 J. Plasma Phys. 81, 905810501) can be used to derive upper bounds on growth rates of local gyrokinetic instabilities (Helander & Plunk 2021 Phys. Rev. Lett. 127, 155001). These bounds follow from thermodynamic principles (specifically on the Helmholtz free energy) (Helander & Plunk Phys. Rev. Lett. 127, 2021, p. 155001), and thus apply to any instability and geometry, independently of many plasma parameters. In this work, we compare these upper bounds with the growth rates of linear gyrokinetic eigenmodes. Experimentally relevant scenarios of density-gradient- and ion-temperature-gradient-driven instabilities are considered. The difference between the upper bounds and the numerically computed growth rates is always positive, as it must be, but depends strongly on the instability in question and on the geometry of the magnetic field. The nature of this difference can be analysed by examining the contributions of optimal modes to gyrokinetic eigenmodes. This approach exploits the completeness and orthogonality properties of optimal modes.

The bootstrap current in stellarators can be presented as a sum of a collisionless value given by the Shaing–Callen asymptotic formula and an off-set current, which non-trivially depends on plasma collisionality and radial electric field. Using NEO-2 modeling, analytical estimates and semi-analytical studies with the help of a propagator method, it is shown that the off-set current in the $1/\nu$ regime does not converge with decreasing collisionality $\nu _\ast$ but rather shows oscillations over $\log \nu _\ast$ with an amplitude of the order of the bootstrap current in an equivalent tokamak. The convergence to the Shaing–Callen limit appears in regimes with significant orbit precession, in particular, due to a finite radial electric field, where the off-set current decreases as $\nu _\ast ^{3/5}$. The off-set current strongly increases in case of nearly aligned magnetic field maxima on the field line where it diverges as $\nu _\ast ^{-1/2}$ in the $1/\nu$ regime and saturates due to the precession at a level exceeding the equivalent tokamak value by ${v_E^\ast }^{-1/2}$, where $v_E^\ast$ is the perpendicular Mach number. The latter off-set, however, can be minimized by further aligning the local magnetic field maxima and by fulfilling an extra integral condition of “equivalent ripples” for the magnetic field. A criterion for the accuracy of this alignment and of ripple equivalence is derived. In addition, the possibility of the bootstrap effect at the magnetic axis caused by the above off-set is also discussed.

When analysing stellarator configurations, it is common to perform an asymptotic expansion about the magnetic axis. This so-called near-axis expansion is convenient for the same reason asymptotic expansions often are, namely, it reduces the dimension of the problem. This leads to convenient and quickly computed expressions of physical quantities, such as quasisymmetry and stability criteria, which can be used to gain further insight. However, it has been repeatedly found that the expansion diverges at high orders in the distance from axis, limiting the physics the expansion can describe. In this paper, we show that the near-axis expansion diverges in vacuum due to ill-posedness and that it can be regularised to improve its convergence. Then, using realistic stellarator coil sets, we demonstrate numerical convergence of the vacuum magnetic field and flux surfaces to the true values as the order increases. We numerically find that the regularisation improves the solutions of the near-axis expansion under perturbation, and we demonstrate that the radius of convergence of the vacuum near-axis expansion is correlated with the distance from the axis to the coils.

Even if the magnetic field in a stellarator is integrable, phase-space integrability for energetic particle guiding-center trajectories is not guaranteed. Both trapped and passing particle trajectories can experience convective losses, caused by wide phase-space island formation, and diffusive losses, caused by phase-space island overlap. By locating trajectories that are closed in the angle coordinate but not necessarily closed in the radial coordinate, we can quantify the magnitude of the perturbation that results in island formation. We characterize island width and island overlap in quasihelical (QH) and quasiaxisymmetric (QA) equilibria with finite plasma pressure $\beta$ for both trapped and passing energetic particles. For trapped particles in QH, low-shear toroidal precession frequency profiles near zero result in wide island formation. While QA transit frequencies do not cross through the zero resonance, we observe that island overlap is more likely since higher shear results in the crossing of more low-order resonances.

In this work, we describe the use of a 1D-2V quasi-neutral hybrid electrostatic PIC with Monte-Carlo Coulomb collisions and non-uniform magnetic field to model the parallel transport and confinement in an axisymmetric tandem mirror device. End-plugs, based on simple-mirrors, are positioned at each end of the device and fueled with neutral beams (25 and 100 keV) to produce a sloshing ion population and increase the density of the end-plugs relative to the central cell. Results show the formation of a potential difference barrier between the central cell and the end-plugs. This potential confines a large fraction of the low energy thermal ions in the central cell which would otherwise be lost in a simple mirror, demonstrating the advantage of the beam-driven tandem mirror configuration relative to simple mirrors. In addition, we explore the effect of end-plug electron temperature on the confinement time of the device and compare it with theoretical estimates. Finally, we discuss the limitations of the code in its present form and describe the next logical steps to improve its predictive capability such as a fully nonlinear Fokker–Planck collision operator, multiply nested flux surface solutions and modeling the exhaust region up to the wall.

Progress in understanding multi-scale collisionless plasma phenomena requires employing tools which balance computational efficiency and physics fidelity. Collisionless fluid models are able to resolve spatio-temporal scales that are unfeasible with fully kinetic models. However, constructing such models requires truncating the infinite hierarchy of moment equations and supplying an appropriate closure to approximate the unresolved physics. Data-driven methods have recently begun to see increased application to this end, enabling a systematic approach to constructing closures. Here, we use sparse regression to search for heat flux closures for one-dimensional electrostatic plasma phenomena. We examine OSIRIS particle-in-cell simulation data of Landau-damped Langmuir waves and two-stream instabilities. Sparse regression consistently identifies six terms as physically relevant, together regularly accounting for more than 95 % of the variation in the heat flux. We further quantify the relative importance of these terms under various circumstances and examine their dependence on parameters such as thermal speed and growth/damping rate. The results are discussed in the context of previously known collisionless closures and linear collisionless theory.

We investigate the one-dimensional non-relativistic Weibel instability through the capture of anisotropic pressure tensor dynamics using an implicit 10-moment fluid model that employs the electromagnetic Darwin approximation. The results obtained from the 10-moment model are compared with an implicit particle-in-cell simulation. The linear growth rates obtained from the numerical simulations are in good agreement with the theoretical fluid and kinetic dispersion relations. The fluid dispersion relations are derived using Maxwell’s equations and the Darwin approximation. We also show that the magnetohydrodynamic approximation can be used to model the Weibel instability if one accounts for an anisotropic pressure tensor and unsteady terms in the generalised Ohm’s law. In addition, we develop a preliminary theory for the saturation magnetic field strength of the Weibel instability, showing good agreement with the numerical results.

We derive a set of simplified equations that can be used for numerical studies of reduced magnetohydrodynamic turbulence within a small patch of the radially expanding solar wind. We allow the box to be either stationary in the Sun’s frame or to be moving at an arbitrary velocity along the background magnetic-field lines, which we take to be approximately radial. We focus in particular on the case in which the box moves at the same speed as outward-propagating Alfvén waves. To aid in the design and optimization of future numerical simulations, we express the equations in terms of scalar potentials and Clebsch coordinates. The equations we derive will be particularly useful for conducting high-resolution numerical simulations of reflection-driven magnetohydrodynamic turbulence in the solar wind, and may also be useful for studying turbulence within other astrophysical outflows.