Overtaking collisions of large-amplitude solitons are investigated via fluid simulations for a plasma consisting of cold ions and Boltzmann-distributed electrons. To achieve this, a new fluid simulation code is presented. In addition, a novel approach to construct soliton initial conditions is developed. Using these ideas, initial conditions are combined that allows the simulation of overtaking collisions. It is shown that, in the small-amplitude regime, simulation results agree well with the two-soliton solution obtained from reductive perturbation theory. Interestingly, in the large amplitude regime, both the slow and fast solitons re-emerge after the collision with no significant change, showing that the collisions remain elastic. A comparison between reductive perturbation analysis and the simulations show that the only significant effect of higher order nonlinearities on overtaking collisions is a reduction in the magnitude of the phase shifts of both solitons.

Electron-cyclotron resonance heating (ECRH) is the main heating mechanism in the Wendelstein 7-X (W7-X) stellarator. Although second-harmonic ECRH (X2) has been used routinely for plasma startup, startup at third harmonic (X3) is known to be much more difficult. In this work, we investigate the energy gain of particles during nonlinear wave–particle interaction for conditions relevant to second- and third-harmonic startups in W7-X. We take into account both the beam and the ambient magnetic field inhomogeneities. The latter is shown to significantly increase the mean energy gain resulting from a single wave–particle resonant interaction. In W7-X-like conditions, the improvement in maximum gained energy is up to 4 times the analogous uniform magnetic field case. However, this improvement is not enough to ensure X3 startup. The optimal magnetic field inhomogeneity length scale for average energy gain and start up in W7-X-like conditions is found to be in the range of $1$ to $3\ {\rm km}^{-1}$. A possibility of using multiple beams with neighbouring resonances is also considered. A considerable enhancement of the energy gain is demonstrated.

The particle-in-cell (PIC) method is successfully used to study magnetized plasmas. However, this requires large computational costs and limits simulations to short physical run times and often to set-ups of less than three spatial dimensions. Traditionally, this is circumvented either via hybrid-PIC methods (adopting massless electrons) or via magneto-hydrodynamic-PIC methods (modelling the background plasma as a single charge-neutral magneto-hydrodynamical fluid). Because both methods preclude modelling important plasma-kinetic effects, we introduce a new fluid-PIC code that couples a fully explicit and charge-conserving multi-fluid solver to the PIC code SHARP through a current-coupling scheme and solve the full set of Maxwell's equations. This avoids simplifications typically adopted for Ohm's law and enables us to fully resolve the electron temporal and spatial scales while retaining the versatility of initializing any number of ion, electron or neutral species with arbitrary velocity distributions. The fluid solver includes closures emulating Landau damping so that we can account for this important kinetic process in our fluid species. Our fluid-PIC code is second-order accurate in space and time. The code is successfully validated against several test problems, including the stability and accuracy of shocks and the dispersion relation and damping rates of waves in unmagnetized and magnetized plasmas. It also matches growth rates and saturation levels of the gyro-scale and intermediate-scale instabilities driven by drifting charged particles in magnetized thermal background plasmas in comparison with linear theory and PIC simulations. This new fluid-SHARP code is specially designed for studying high-energy cosmic rays interacting with thermal plasmas over macroscopic time scales.

The gyrokinetic theory of the residual flow, in the electrostatic limit, is revisited, with optimized stellarators in mind. We consider general initial conditions for the problem, and identify cases that lead to a non-zonal residual electrostatic potential, i.e. one having a significant component that varies within a flux surface. We investigate the behaviour of the ‘intermediate residual’ in stellarators, a measure of the flow that remains after geodesic acoustic modes have damped away, but before the action of the slower damping that is caused by unconfined particle orbits. The case of a quasi-isodynamic stellarator is identified as having a particularly large such residual, owing to the small orbit width achieved by optimization.

We perform experiments in the Large Plasma Device at the University of California, Los Angeles, studying how different end-electrode biasing schemes modify the radial potential profile in the machine. We impose biasing profiles of different polarities and gradient signs on a set of five concentric electrodes placed 12 m downstream from the plasma source. We find that imposing concave-down profiles (negative potential radial gradient) on the electrodes creates radial potential profiles halfway up the plasma column that are comparable to those imposed on the electrodes and a few electron temperature in height, regardless of the biasing polarity. On the other hand, imposing concave-up profiles (positive potential radial gradient) leads to non-monotonic radial potential profiles. This observation can be explained by the current drawn through the electrodes and the parallel plasma resistivity, highlighting their important role in controlling the rotation of plasma. Concave-down plasma potential profiles, obtained by drawing electrons on the axis, are predicted to drive azimuthal drift velocities that can approach significant fractions of the ion sound speed in the central region of the plasma column.

A numerical study of the homo-interactions between two falling droplets and between two rising bubbles in a strongly coupled dusty plasma medium is presented in this article. The strongly coupled dusty plasma is considered as a viscoelastic fluid using the generalized hydrodynamic fluid model formalism. Two factors that affect homo-interactions are taken into account: the initial spacing and the coupling strength of the medium. Three different spacings between two droplets are simulated: widely, medium and closely. In each case, the coupling strength has been given as mild–strong and strong. It is shown that the overall dynamic is governed by the competition between the acceleration of two droplets/bubbles due to gravity and the interaction due to the closeness of the droplets/bubbles. Particularly in viscoelastic fluids, apart from the initial separation, shear waves originating from rotating vortices are responsible for the closeness of two droplets or bubbles. Several two-dimensional simulations have been carried out. This work is a continuation of the work done in Parts 1 (Dharodi & Das, J. Plasma Phys., vol. 87, issue 2, 2021, 905870216) and 2 (Dharodi, J. Plasma Phys., vol. 87, issue 4, 2021, 905870402).

Hot electrons produced by parametric instabilities can negatively affect the thermonuclear gain of a direct-drive inertial confinement fusion capsule. A Monte Carlo transport is coupled to the hydrodynamics MULTI-IFE code in order to study the hot electron transport in a high-gain target implosion. The thermonuclear energy produced by the directly driven implosion of a reference target drops from 27 MJ without hot electrons to 13.5 kJ when the hot electrons are taken into account. A clear relationship is established between the level of hot electrons produced and the degradation of the thermonuclear yield produced by the target. It is shown that adding a relatively short laser spike (shock ignition) to the main pulse restores the thermonuclear energy released by the implosion. Different shock-ignition windows are presented depending on the duration of the spike. It appears that the longer the duration of the spike, the lower the power required for shock ignition. Finally, a 300 ps-duration spike requires only 50 TW to restore a thermonuclear energy output of 20 MJ.

When simulating runaway electron dynamics in tokamak disruptions, fluid models with lower numerical cost are often preferred to more accurate kinetic models. The aim of this work is to compare fluid and kinetic simulations of a large variety of different disruption scenarios in ITER. We consider both non-activated and activated scenarios; for the latter, we derive and implement kinetic sources for the Compton scattering and tritium beta decay runaway electron generation mechanisms in our simulation tool Dream (Hoppe et al., Comput. Phys. Commun., vol. 268, 2021, 108098). To achieve a diverse set of disruption scenarios, Bayesian optimization is used to explore a range of massive material injection densities for deuterium and neon. The cost function is designed to distinguish between successful and unsuccessful disruption mitigation based on the runaway current, current quench time and transported fraction of the heat loss. In the non-activated scenarios, we find that fluid and kinetic disruption simulations can have significantly different runaway electron dynamics, due to an overestimation of the runaway seed by the fluid model. The primary cause of this is that the fluid hot-tail generation model neglects superthermal electron transport losses during the thermal quench. In the activated scenarios, the fluid and kinetic models give similar predictions, which can be explained by the significant influence of the activated sources on the runaway dynamics and the seed.

Microinstabilities drive turbulent fluctuations in inhomogeneous, magnetised plasmas. In the context of magnetic confinement fusion devices, this leads to an enhanced transport of particles, momentum and energy, thereby degrading confinement. In this work, we describe an application of the adjoint method to efficiently determine variations of gyrokinetic linear growth rates on a general set of external parameters in the local $\delta f$-gyrokinetic model. We then offer numerical verification of this approach. When coupled with gradient-based techniques, this methodology can facilitate the optimisation process for the microstability of the confined plasmas across a high-dimensional parameter space. We present a numerical demonstration wherein the ion-temperature-gradient instability growth rate in a tokamak plasma is minimised with respect to flux surface shaping parameters. The adjoint method approach demonstrates a significant computational speed-up compared with a finite-difference gradient calculation.

Quasilinear treatments are widely used for tokamaks to evaluate radio frequency (rf) heating and current drive. Even though the core of a tokamak plasma is weakly collisional, the solution of the linearized kinetic equation is evaluated using unperturbed collisionless trajectories while often treating successive poloidal circuits of the passing (and trapped) particles as uncorrelated or nearly so. In addition, the most important effect of tokamak geometry, the mirror force, is usually mistreated or ignored when obtaining the solution. These concerning aspects of rf treatments are clarified by considering lower hybrid heating and current drive to illustrate that the electrons in resonance with the applied rf are enclosed by narrow collisional boundary layers, and that tokamak geometry makes it necessary to retain poloidal variation when solving a weakly collisional linearized kinetic equation. Other aspects such as collisional boundary layers at the trapped–passing boundary, cyclotron resonances, and the limitations of quasilinear theory are also considered. The new insights lead to a fundamentally different formulation and interpretation of the solution of the linearized Fokker–Planck equation used for rf quasilinear theory in a tokamak, while retaining many of the features that have contributed to its successful application to rf heating and current drive.

Large-amplitude current-driven plasma instabilities, which can transition to the Buneman instability, were observed in one-dimensional simulations to generate high-energy back-streaming ions. We investigate the saturation of multi-dimensional plasma instabilities and its effects on energetic ion formation. Such ions directly impact spacecraft thruster lifetimes and are associated with magnetic reconnection and cosmic ray inception. An Eulerian Vlasov–Poisson solver employing the grid-based direct kinetic method is used to study the growth and saturation of 2D2V collisionless, electrostatic current-driven instabilities spanning two dimensions each in the configuration (D) and velocity (V) spaces supporting ion and electron phase-space transport. Four stages characterise the electric potential evolution in such instabilities: linear modal growth, harmonic growth, accelerated growth via quasi-linear mechanisms alongside nonlinear fill-in and saturated turbulence. Its transition and isotropisation process bears considerable similarities to the development of hydrodynamic turbulence. While a tendency to isotropy is observed in the plasma waves, followed by electron and then ion phase spaces after several ion-acoustic periods, the formation of energetic back-streaming ions is more limited in the 2D2V than in the 1D1V simulations. Plasma waves formed by two-dimensional electrostatic kinetic instabilities can propagate in the direction perpendicular to the net electron drift. Thus, large-amplitude multi-dimensional waves generate high-energy transverse-streaming ions and eventually limit energetic backward-streaming ions along the longitudinal direction. The multi-dimensional study sheds light on interactions between longitudinal and transverse electrostatic plasma instabilities, as well as fundamental characteristics of the inception and sustenance of unmagnetised plasma turbulence.

Filament-based coil optimizations are performed for several quasi-helical stellarator configurations, beginning with the one from Landreman & Paul (Phys. Rev. Lett., vol. 128, 2022, 035001), demonstrating that precise quasi-helical symmetry can be achieved with realistic coils. Several constraints are placed on the shape and spacing of the coils, such as low curvature and sufficient plasma–coil distance for neutron shielding. The coils resulting from this optimization have a maximum curvature 0.8 times that of the coils of the Helically Symmetric eXperiment (HSX) and a mean squared curvature 0.4 times that of the HSX coils when scaled to the same plasma minor radius. When scaled up to reactor size and magnetic field strength, no fast particle losses were found in the free-boundary configuration when simulating 5000 alpha particles launched at $3.5\,\mathrm {MeV}$ on the flux surface with a normalized toroidal flux of $s=0.5$. An analysis of the tolerance of the coils to manufacturing errors is performed using a Gaussian process model, and the coils are found to maintain low particle losses for smooth, large-scale errors up to amplitudes of approximately $0.15\,\mathrm {m}$. Another coil optimization is performed for the Landreman–Paul configuration with the additional constraint that the coils are purely planar. Visual inspection of the Poincaré plot of the resulting magnetic field-lines reveal that the planar modular coils alone do a poor job of reproducing the target equilibrium. Additional non-planar coil optimizations are performed for the quasi-helical configuration with $5\,\%$ volume-averaged plasma beta from Landreman et al. (Phys. Plasma, vol. 29, issue 8, 2022, 082501), and a similar configuration also optimized to satisfy the Mercier criterion. The finite beta configurations had larger fast-particle losses, with the free-boundary Mercier-optimized configuration performing the worst, losing approximately $5.5\,\%$ of alpha particles launched at $s=0.5$.

The excitation conditions of the magnetorotational instability (MRI) are studied for axially unbounded Taylor–Couette (TC) flows of various gap widths between the cylinders. The cylinders are considered as made from both perfect-conducting or insulating material and the conducting fluid with a finite but small magnetic Prandtl number rotates with a quasi-Keplerian velocity profile. The solutions are optimized with respect to the wavenumber and the Reynolds number of the rotation of the inner cylinder. For the axisymmetric modes, we find the critical Lundquist number of the applied axial magnetic field: the lower, the wider the gap between the cylinders. A similar result is obtained for the induced cell structure: the wider the gap, the more spherical the cells are. The marginal rotation rate of the inner cylinder – for a fixed size of the outer cylinder – always possesses a minimum for not too wide and not too narrow gap widths. For perfect-conducting walls the minimum lies at $r_{{\rm in}}\simeq 0.4$, where $r_{{\rm in}}$ is the ratio of the radii of the two rotating cylinders. The lowest magnetic field amplitudes to excite the instability are required for TC flows between perfect-conducting cylinders with gaps corresponding to $r_{{\rm in}}\simeq ~0.2$. For even wider and also for very thin gaps the needed magnetic fields and rotation frequencies are shown to become rather huge. Also the non-axisymmetric modes with $|m|=1$ have been considered. Their excitation generally requires stronger magnetic fields and higher magnetic Reynolds numbers in comparison with those for the axisymmetric modes. If TC experiments with too slow rotation for the applied magnetic fields yield unstable modes of any azimuthal symmetry, such as the currently reported Princeton experiment (Wang et al., Phys. Rev. Lett., vol. 129, 115001), then also other players, including axial boundary effects, than the MRI-typical linear combination of current-free fields and differential rotation should be in the game.

Between a plasma and a solid target lies a positively charged sheath of several Debye lengths $\lambda _{D}$ width, typically much smaller than the characteristic length scale $L$ of the main plasma. This scale separation implies that the asymptotic limit $\epsilon = \lambda _{D} / L \rightarrow 0$ is useful to solve for the plasma-sheath system. In this limit, the Bohm criterion must be satisfied at the sheath entrance. A new derivation of the kinetic criterion, admitting a general ion velocity distribution, is presented. It is proven that, for $\epsilon \rightarrow 0$, the distribution of the velocity component normal to the target, $v_x$, and its first derivative must vanish for $|v_x| \rightarrow 0$ at the sheath entrance. These two conditions can be subsumed into a third integral one after it is integrated by parts twice. A subsequent interchange of the limits $\epsilon \rightarrow 0$ and $|v_x| \rightarrow 0$ is invalid, leading to a divergence which underlies the misconception that the criterion gives undue importance to slow ions.

We analysed the excitation of a surface magnetoplasmon wave by the mode conversion of a p-polarized laser beam over a rippled semiconductor (n-type)-free space interface. The pump surface magnetoplasmon wave exerts a ponderomotive force on the free electrons in the semiconductor, imparting a linear oscillatory velocity at the laser modulation frequency to them. This linear oscillatory velocity couples with the modulated electron density to produce a current density, which develops a resonant surface magnetoplasmon wave in the terahertz region. The amplitude of the terahertz surface magnetoplasmon wave can be tuneable with an external magnetic field and the semiconductor's temperature.

Recent developments in the design of magnetic confinement fusion devices have allowed the construction of exceptionally optimized stellarator configurations. The near-axis expansion in particular has been proven to enable the construction of magnetic configurations with good confinement properties while taking only a fraction of the usual computation time to generate optimized magnetic equilibria. However, not much is known about the overall features of fast-particle orbits computed in such analytical, yet simplified, equilibria when compared with those originating from accurate equilibrium solutions. This work aims to assess and demonstrate the potential of the near-axis expansion to provide accurate information on particle orbits and to compute loss fractions in moderate to high aspect ratios. The configurations used here are all scaled to fusion-relevant parameters and approximate quasi-symmetry to various degrees. This allows us to understand how deviations from quasi-symmetry affect particle orbits and what are their effects on the estimation of the loss fraction. Guiding-centre trajectories of fusion-born alpha particles are traced using gyronimo and SIMPLE codes under the NEAT framework, showing good numerical agreement. Discrepancies between near-axis and magnetohydrodynamic fields have minor effects on passing particles but significant effects on trapped particles, especially in quasi-helically symmetric magnetic fields. Effective expressions were found for estimating orbit widths and passing–trapped separatrix in quasi-symmetric near-axis fields. Loss fractions agree in the prompt losses regime but diverge afterwards.

Laser-driven neutron sources are routinely produced by the interaction of laser-accelerated protons with a converter. They present complementary characteristics to those of conventional accelerator-based neutron sources (e.g. short pulse durations, enabling novel applications like radiography). We present here results from an experiment aimed at performing a global characterization of the neutrons produced using the Titan laser at the Jupiter Laser Facility (Livermore, USA), where protons were accelerated from 23 $\mathrm {\mu }$m thick plastic targets and directed onto a LiF converter to produce neutrons. For this purpose, several diagnostics were used to measure these neutron emissions, such as CR-39, activation foils, time-of-flight detectors and direct measurement of $^7{\rm Be}$ residual activity in the LiF converters. The use of these different, independently operating diagnostics enables comparison of the various measurements performed to provide a robust characterization. These measurements led to a neutron yield of $2.0\times 10^{9}$ neutrons per shot with a modest angular dependence, close to that simulated.

A fusion neutron source (FNS) based on the gas-dynamic trap (GDT, Budker Institute, Novosibirsk) is considered for confinement of two-species plasma heated by neutral beam injection in a regime where the fast ion distribution function is far from Maxwellian. Kinetic instabilities are expected to develop in this regime, and in this paper we investigate the ion-cyclotron instability evolving in moderate densities of pure hydrogen and mixed deuterium–hydrogen target plasmas. The properties of the studied unstable mode, such as its azimuthal wavenumbers, propagation direction and its being affected by changes in the bulk plasma density and composition, allow us to identify it as the drift cyclotron loss cone (DCLC) instability. This mode scatters fast ions and thereby leads to drops in diamagnetic flux signals and increases longitudinal energy and particle losses, with the average energy of the lost ions estimated to be far above the temperature of warm Maxwellian ions. Our interpretation is that the unstable wave grows due to interaction with the fast ions located near the loss cone in the velocity space and scatters them. Applying the method of suppressing the DCLC instability by filling the loss cone with warm plasma, we have determined the values of plasma density and deuterium percentage that allow us to suppress the DCLC instability in the GDT. These findings justify using mixed bulk plasmas in fusion neutron source operation.

A systematic theory of the asymptotic expansion of the magnetohydrostatics (MHS) equilibrium in the distance from the magnetic axis is developed to include arbitrary smooth currents near the magnetic axis. Compared with the vacuum and the force-free system, an additional magnetic differential equation must be solved to obtain the pressure-driven currents. It is shown that there exist variables in which the rest of the MHS system closely mimics the vacuum system. Thus, a unified treatment of MHS fields is possible. The mathematical structure of the near-axis expansions to arbitrary order is examined carefully to show that the double-periodicity of physical quantities in a toroidal domain can be satisfied order by order. The essential role played by the leading-order Birkhoff–Gustavson normal form in solving the magnetic differential equations is highlighted. Several explicit examples of vacuum, force-free and MHS equilibrium in different geometries are presented.

Electron heating and acceleration in collisionless shocks is a long-standing problem. Rapid isotropization of heated electrons cannot be explained solely by the cross-shock potential. In addition, the macroscopic cross-shock potential prevents efficient reflection and injection into the diffusive acceleration regime. Recent observations have shown that small-scale electric fields are present in the shock front, together with the large-scale cross-shock potential. These small-scale fields have been found also in the upstream and downstream regions. Electron heating in shocks is produced by the combined action of the large- and small-scale fields. The large-scale potential determines the energy transferred to the electrons. The small-scale electrostatic fields scatter electrons. Here we study the scattering of electrons on the typical waveforms, namely solitary bipolar spikes and wavepackets. The main effect is the generation of backstreaming electrons with large pitch angles. It is found that wavepackets are more efficient in electron reflection in the interaction of electrons both with a single spike and with multiple spikes.