The centrifugal mirror confinement scheme incorporates supersonic rotation of a plasma into a magnetic mirror device. This concept has been shown experimentally to drastically decrease parallel losses and increase plasma stability as compared with prior axisymmetric mirrors. MCTrans++ is a dimensionless (0-D) scoping tool which rapidly models experimental operating points in the Centrifugal Mirror Fusion Experiment (CMFX) at the University of Maryland. In the low-collisionality regime, parallel losses can be modelled analytically. A confining potential is set up that is partially ambipolar and partially centrifugal. Due to the stabilizing effects of flow shear, the perpendicular losses can be modelled as classical. Radiation losses such as bremsstrahlung and cyclotron emission are taken into account. A neutrals model is included, and, in some circumstances, charge-exchange losses are found to exceed all other loss mechanisms. We use the SUNDIALS ARKODE library to solve the underlying equations of this model; the resulting software is suitable for scanning large parameter spaces, and can also be used to model time-dependent phenomena such as a capacitive discharge. MCTrans++ has been used to verify results from prior centrifugal mirrors, create an experimental plan for CMFX and find configurations for future reactor-scale fusion devices.

BOUT++ turbulence simulations were performed to investigate the impact of turbulence spreading on the edge localized mode (ELM) size and divertor heat flux width $({\lambda _q})$ broadening in small ELM regimes. This study is motivated by EAST experiments. BOUT++ linear simulations of a pedestal radial electric field (Er) scan show that the dominant toroidal number mode (n) shifts from high-n to low-n, with a narrow mode spectrum, and the maximum linear growth rate increases as the pedestal Er well deepens. The nonlinear simulations show that as the net E × B pedestal flow increases, the pressure fluctuation level and its inward penetration beyond the top of the pedestal both increase. This leads to a transition from small ELMs to large ELMs. Both inward and outward turbulence spreading are sensitive to the scrape-off-layer (SOL) plasma profiles. The inward turbulence spreading increases for the steep SOL profiles, leading to increasing pedestal energy loss in the small ELM regime. The SOL width $({\lambda _q})$ is significantly broadened progressing from the ELM-free to small ELM regime, due to the onset of strong radial turbulent transport. The extent of the SOL width $({\lambda _q})$ broadening depends strongly on outward turbulence spreading. The fluctuation energy intensity flux ${\varGamma _\varepsilon }$ at the separatrix can be enhanced by increasing either pedestal Er flow shear or local SOL pressure gradient. The ${\lambda _q}$ is broadened as the fluctuation energy intensity flux ${\varGamma _\varepsilon }$ at the last close flux surface (LCFS) increases. Local SOL E × B flow shear will restrain outward turbulence spreading and the associated heat flux width broadening. Operating in H-mode with small ELMs has the potential to solve two critical problems: reducing the ELM size and broadening the SOL width.

This paper is based on a report at the 2nd International Fusion Plasma Conference & 13th International Conference on Open Magnetic Systems for Plasma Confinement (iFPC & OS 2023), August 21–25, 2023, Busan, Korea and provides a brief overview of the status of work at the Budker Institute on the study of hot plasma confinement in open-type magnetic traps with a linear axisymmetric configuration. The main attention is paid to key problems: magnetohydrodynamics (MHD) stability in regimes with extremely high relative pressure, longitudinal electronic thermal conductivity, stability with respect to the development of kinetic modes and transverse transport. This paper provides an overview of the methods we are developing to address these problems, the experimental and theoretical results achieved and plans for future development. The last section of the article provides brief information about the preliminary design of the gas-dynamic multiple-mirror trap device, the development of which has been completed.

Laser-driven ion acceleration is well established using solid targets mainly in the target normal sheath acceleration regime. To follow the increasing repetition rate available on high-intensity lasers, the use of high-density gas targets has been explored in the past decade. When interacting with targets reaching densities close to the critical one, the laser pulse can trigger different acceleration mechanisms such as Collisionless Shock Acceleration (CSA) or hole boring. Particle-in-cell simulations using ideal target profiles show that CSA can accelerate a collimated, narrow energy spread and few hundreds of megaelectronvolts ion beam on the laser axis. Nevertheless, in real experiments, the laser will not only interact with an overcritical, thin plasma slab with sharp density gradients, but also with lower density regions surrounding the core of the gas jet, extending to several hundreds of micrometres. The interaction of the laser with these lower density wings will lead to nonlinear effects that will reduce the available energy to drive the shock in the high-density region of the target. Optically tailoring this target could mitigate that issue. Recent experiments conducted on different laser facilities aimed at testing several tailoring configurations. We first tested a scheme with a copropagating picosecond prepulse to create a lower density plasma channel to facilitate the propagation of the main pulse, while the second one was a transverse tailoring driven by nanosecond laser pulses to generate blast waves and form a high-density plasma slab. The main results will be presented here and the methods compared.

Recently, stochastic motion of 3.5-MeV alpha particles with orbits that vary between locally trapped and locally passing states (transitioning particles) in a Helias reactor was observed numerically. This validated the theoretical predictions that (i) the stochastic diffusion represents a mechanism of considerable delayed loss of fast ions in quasi-isodynamic stellarators and (ii) it is possible to prevent the escape of particles to the wall by closing the separatrix between the locally trapped and passing states. It was concluded that, in principle, the separatrix could be made closed, resulting in a reduction of fast ion losses, by compensating for the effect of the helical component of the magnetic field $(1, 1)$ with an enhanced ‘anti-helical’ harmonic $(1, -1)$; the enlargement of this harmonic was proposed previously for other reasons. This possibility is explored in this work. Equations of previous relevant works were generalised to include the $(1, -1)$ harmonic. Calculations were carried out for several magnitudes of the ratio of anti-helical to helical magnetic field harmonics. Positive effects were found already at the smallest anti-helical harmonic considered: when the latter ratio is 0.25, transitioning particles with the smallest and intermediate pitch parameters are confined and, moreover, their fraction decreased. When the ratio is 0.85, almost all transitioning particles are confined and their fraction is minimal; well-confined localised orbits dominate at the smallest pitch parameters.

It is known that a nonlinear Schrödinger equation describes the self-modulation of a large amplitude circularly polarized wave in relativistic electron–positron plasmas in the weakly and strongly magnetized limits. Here, we show that such an equation can be written as a modified second Painlevé equation, producing accelerated propagating wave solutions for those nonlinear plasmas. This solution even allows the plasma wave to reverse its direction of propagation. The acceleration parameter depends on the plasma magnetization. This accelerating solution is different to the usual soliton solution propagating at constant speed.

The paper presents an investigation of the plasma fluctuation in the SMOLA helical mirror, which is suspected to be responsible for anomalous scattering. The helical mirror confinement is effective when the ion mean free path is equal to the helix pitch length. This condition can be satisfied in hot collisionless plasma only by anomalous scattering. The wave, which scatters the passing ions, is considered to receive energy from the trapped ions. The oscillations of the electric field in the helically symmetric plasma were observed in experiment. The oscillations have both regular highly correlated and chaotic components. The dependency of the regular component frequency on the Alfvén velocity is linear for $V_{\rm A} < 2.8 \times 10^6\ \text {m}\ \text {s}^{-1}$ and constant for higher values. It is shown experimentally that the condition for the wave to be in phase resonance with the trapped ions is satisfied in a specific region of the plasma column for the highly correlated component. The amplitude of the chaotic component (up to $3\ \text {V}\ \text {cm}^{-1}$) is higher than the estimated electric field required for the ion scattering.

Omnigenity is a desirable property of toroidal magnetic fields that ensures confinement of trapped particles. Confining charged particles is a basic requirement for any fusion power plant design, but it can be difficult to satisfy with the non-axisymmetric magnetic fields used by the stellarator approach. Every ideal magnetohydrodynamic equilibrium previously found to approximate omnigenity has been either axisymmetric, quasi-symmetric or has poloidally closed contours of magnetic field strength $B$. However, general omnigenous equilibria are a much larger design space than these subsets. A new model is presented and employed in the DESC stellarator optimization suite to represent and discover the full parameter space of omnigenous equilibria. Although exact omnigenity aside from quasi-symmetry is impossible, these results reveal that excellent particle confinement can be achieved in practice. Examples far from quasi-symmetry with poloidally, helically and toroidally closed $B$ contours are attained with DESC and shown to have low neoclassical collisional transport and fast particle losses.

In this paper, the modelling and manufacturing of a two-stage supersonic gas jet nozzle enabling the formation of adaptive plasma concentration profiles for injection and acceleration of electrons using few-cycle laser beams are presented. The stages are modelled using the rhoSimpleFoam algorithm of the OpenFOAM computational fluid dynamics software. The first 200–300 ${\rm \mu}$m diameter nozzle stage is dedicated to 1 % N2 + He gas jet formation and electron injection. By varying the pressure between the first and second stages of the injectors, the electron injection location could be adjusted, and the maximum acceleration distance could be ensured. By changing the concentration of the nitrogen in the gas mixture, the charge of the accelerated electrons could be controlled. The second nozzle stage is designed for acceleration in fully ionised He or hydrogen gas and forms the optimal plasma concentration for bubble formation depending on the laser pulse energy, duration and focused beam diameter. In order to reduce the diameter of the plasma profile formed by the first nozzle and the concentration drop gap between the two nozzles, a one-side straight section was introduced in the first nozzle. The shock wave reflected from the straight section of the wall propagates parallel to the shock wave of the intersecting supersonic jets and ensures a minimal gap between the jets. The second-stage longitudinal plasma concentration profile could have an increasing gas density gradient to compensate for dephasing between the electron bunch and the plasma wave due to wave shortening with increasing plasma concentration.

A modulational instability of nonlinearly interacting electron whistlers and magnetosonic perturbations is studied in the present paper. For typical parameters, there is no modulational instability. However, modulational instability appears in special cases. For example, when the whistler wavenumber is small enough, there is modulational instability. Its growth rate decreases as the angle between the external magnetic field and the perturbed wave's direction increases, while it increases as the whistler wavenumber increases. It is also found that there is no modulational instability when the whistler wavenumber is larger than a critical value ($k_0 > 0.05$), in which the perturbed wave frequency increases as the angle between the external magnetic field and the perturbed wave's direction increases when the angle between the external magnetic field and the perturbed wave's direction is large enough. Whereas, the perturbed wave frequency first increases as the whistler wavenumber increases, reaches a peak value and then decreases as whistler wavenumber increases.

A high-speed solenoid-type pulsed valve is applied to a water-fuelled magnetron sputtering thruster, which can lead a compact and neutralizer-free electric propulsion device by generating a thrust due to the ejection of the sputtered target material, in order to minimize the size of the gas feeding system. The pulsed valve is opened for approximately 5 msec and the gaseous water is introduced into the source via a gas line connected to the source. A dc high voltage is firstly applied between the cathode and the anode; then the discharge is sustained for the period in which a gas pressure enough for the discharge is maintained and a change in the discharge mode is observed during the discharge pulse. To investigate the impact of the mode change on the thrust generation, a temporally resolved measurement of the force is equivalently performed by installing a pulsed voltage power supply and by changing the pulse width of the discharge voltage. The results show the temporal reduction of the thrust-to-power ratio during a pulsed discharge due to the nonlinear correlation between ion energy and sputtering yield. It is suggested that the use of a pulsed valve would be useful for downsizing of the electric propulsion system, while an optimization to properly control the gas flow rate is needed for further development.

The frequency spectra of the Trivelpiece–Gould modes of a waveguide partially filled with non-neutral plasma are determined numerically by solving the dispersion equation. The modes having azimuthal number $m = 1$ are considered. The results are presented for the entire acceptable range of electron densities, magnetic field strengths, for different values of the charge neutralization coefficient. The Cherenkov resonance condition of an ion with a diocotron mode having a finite value of the longitudinal wave vector was studied. The characteristics of resonant low-frequency electron–ion instability caused by relative azimuth motion of electrons and ions in crossed fields and by the anisotropy of the distribution function of ions are discussed. Ions are created by ionization of residual gas in the plasma volume. Due to the anisotropy, instability occurs not only in the vicinity of the resonance, but also outside it. For typical values of plasma parameters in experiments, estimations of the frequency growth rate are given. A conclusion is drawn that this instability can be the cause of the low-frequency oscillations observed in linear devices with non-neutral plasma produced in an electron beam channel.

This paper is devoted to experimental studies of plasma confinement with high relative pressure ($\beta$) in the Gas Dynamic Trap (BINP, Novosibirsk). In previous high-$\beta$ confinement studies a maximum local $\beta = 0.6$ was achieved in the fast-ion turning point, contributed to by a beam-driven population of fast ions with an anisotropic distribution function. In this study the axial magnetic field profile was modified to bring the turning points closer to one another, which effectively increased the energy density of plasma and pushed the $\beta$ value higher. Experiments were performed for two non-standard magnetic configurations, where the axial fast-ion confinement region length was reduced by 1.5 and 2 times compared with the standard configuration. The average values of $\langle \beta _{\perp } \rangle$ over the plasma central cross-section were found to be 0.1 and 0.18, respectively, for the two configurations, with the latter value significantly exceeding the $\langle \beta _{\perp } \rangle =0.08$ of the standard configuration, in which the previous record was set. Moreover, halving the fast ion confinement region almost doubled the D–D fusion proton flux from the trap centre compared with the standard configuration. The electron temperature in both new magnetic configurations was only slightly smaller than in the standard configuration. In addition, an effect of Alfvén ion–cyclotron instability (AICI) development on the pressure in the turning points is discussed. Presumably, with some decrease in magnetic field an evolving AICI does not result in considerable pressure axial redistribution, so the pressure maximum is in the turning points’ vicinity despite the instability.

In this paper, we investigate the kinetic stability of classical, collisional plasma – that is, plasma in which the mean-free-path $\lambda$ of constituent particles is short compared with the length scale $L$ over which fields and bulk motions in the plasma vary macroscopically, and the collision time is short compared with the evolution time. Fluid equations are typically used to describe such plasmas, since their distribution functions are close to being Maxwellian. The small deviations from the Maxwellian distribution are calculated via the Chapman–Enskog (CE) expansion in $\lambda /L \ll 1$, and determine macroscopic momentum and heat fluxes in the plasma. Such a calculation is only valid if the underlying CE distribution function is stable at collisionless length scales and/or time scales. We find that at sufficiently high plasma $\beta$, the CE distribution function can be subject to numerous microinstabilities across a wide range of scales. For a particular form of the CE distribution function arising in strongly magnetised plasma (viz. plasma in which the Larmor periods of particles are much smaller than collision times), we provide a detailed analytic characterisation of all significant microinstabilities, including peak growth rates and their associated wavenumbers. Of specific note is the discovery of several new microinstabilities, including one at sub-electron-Larmor scales (the ‘whisper instability’) whose growth rate in certain parameter regimes is large compared with other instabilities. Our approach enables us to construct the kinetic stability maps of classical, two-species collisional plasma in terms of $\lambda$, the electron inertial scale $d_e$ and the plasma $\beta$. This work is of general consequence in emphasising the fact that high-$\beta$ collisional plasmas can be kinetically unstable; for strongly magnetised CE plasmas, the condition for instability is $\beta \gtrsim L/\lambda$. In this situation, the determination of transport coefficients via the standard CE approach is not valid.

We investigate the linear and nonlinear evolution of the current-driven ion-acoustic instability in a collisionless plasma via two-dimensional (2-D) Vlasov–Poisson numerical simulations. We initialise the system in a stable state and gradually drive it towards instability with an imposed, weak external electric field, thus avoiding physically unrealisable super-critical initial conditions. A comprehensive analysis of the nonlinear evolution of ion-acoustic turbulence (IAT) is presented, including the detailed characteristics of the evolution of the particles’ distribution functions, (2-D) wave spectrum and the resulting anomalous resistivity. Our findings reveal the dominance of 2-D quasi-linear effects around saturation, with nonlinear effects, such as particle trapping and nonlinear frequency shifts, becoming pronounced during the later stages of the system's nonlinear evolution. Remarkably, the Kadomtsev–Petviashvili (KP) spectrum is observed immediately after the saturation of the instability. Another crucial and noteworthy result is that no steady saturated nonlinear state is ever reached: strong ion heating suppresses the instability, which implies that the anomalous resistivity associated with IAT is transient and short-lived, challenging earlier theoretical results. Towards the conclusion of the simulation, electron-acoustic waves are triggered by the formation of a double layer and strong modifications to the particle distribution induced by IAT.