Generalising the Elsässer variables, we introduce the $Q$-variables. These are more flexible than the Elsässer variables, because they also allow us to track waves with phase speeds different than the Alfvén speed. We rewrite the magnetohydrodynamics (MHD) equations with these $Q$-variables. We consider also the linearised version of the resulting MHD equations in a uniform plasma, and recover the classical Alfvén waves, but also separate the fast and slow magnetosonic waves into upward- and downward-propagating waves. Moreover, we show that the $Q$-variables may also track the upward- and downward-propagating surface Alfvén waves in a non-uniform plasma, displaying the power of our generalisation. In the end, we lay the mathematical framework for driving solar wind models with a multitude of wave drivers.

We investigate the effects of different guard geometries on the currents to needle-type Langmuir probes. The results are based on particle-in-cell numerical simulations. We show that if the guard length is less than 6–8 Debye lengths there can be a significant effect on the currents to the probe. A guard radius should not be larger than the Debye length, otherwise it can also significantly affect the currents. However, since guard radii are often close to the probe radius, the second condition is usually satisfied.

This study investigates the effect of a transverse magnetic field on high-voltage pulsed discharge in helium at a pressure of 30 Torr. A simple two-dimensional fluid model that describes the high-voltage pulsed discharge in helium in a transverse weak magnetic field (B = 0.4 T) is presented, which uses an empirical relation to account for the magnetic field. The results of using the empirical relation for the effective field agree well with the experimental results. The dynamics of discharge development in the presence of the magnetic field is also investigated. The magnetic field does not significantly affect the gas-discharge development dynamics in helium at a pressure of 30 Torr.

Radio frequency quadrupole coolers (RFQCs) are very suitable to cool ion beams with moderate energy spread, typically ions of exotic nuclear species (like $^{132}$Sn$^{1+}$) as in the Selective Production of Exotic Species project at the Laboratori Nazionali di Legnaro, whose ion source supplies 40 keV ions. Beam dynamics includes ion–gas collisions (with a balance of cooling and diffusion effects), acceleration and deceleration and radiofrequency confinement, which can be supplemented by static magnetic field effects. Insertion of a prototype RFQC in the solenoid of the Eltrap machine is also discussed here, with innovations in beam extraction and in modelling, now also based on stochastic equations. Practical consideration on gas pumping and voltage distribution are also included. Typical limits of RFQC are discussed, with special attention to the extracted beam root mean square emittance, which is shown to strongly depend not only on cooler parameters, but also on extraction optics.

The theory of the terahertz (THz) waves emission from a near-critical plasma slab under the action of the focused p-polarized laser pulse is developed. The spectral, angular and energy characteristics of the THz signal are studied as functions of the focusing degree and the incidence angle of laser radiation, as well as the density and thickness of the plasma slab. It is shown that the extremely strong increase in the energy of the THz signal (up to millijoule level) and conversion rate (up to 10 %) occurs at the almost normal incidence of the ultra-short, tightly focused p-polarized laser pulse on the thin plasma slab with the near-critical density and rare electron collisions.

The heating rate of plasma electrons induced by external fields or other processes can be used as an experimental tool to measure fundamental plasma properties such as electrical conductivity or electron–ion collision rates. We have developed a technique that can measure electron heating rates in ultracold neutral plasmas (UNPs) with $\sim 10\,\%$ precision while simultaneously referencing the measurement to a calibrated amount of heating. This technique uses a sequence of applied electric fields in four sections: to control the ratio of electrons to ions in the UNP; to provide a time for the application of fields that cause electron heating and subsequent thermalization of the electrons after the application of those fields; to extract electrons from the UNP using a method sensitive to electron temperature that allows the measurement of electron heating; and to extract the remaining electrons to measure the total electron (and therefore ion) number. The primary signal used to measure the heating rate is the measurement of the number of electrons that escape in the third section of the experiment as a larger number of escaping electrons indicates a larger amount of heating. We illustrate the use of this technique by measuring electron heating caused by high-frequency radiofrequency (RF) fields. In addition to the main technique, several subtechniques to calibrate the electron temperature, electron density, amount of heating and applied RF field amplitude were developed as well.

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.

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.

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.

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.

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.

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.

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.