In hydrodynamic (HD) turbulence, an exact decomposition of the energy flux across scales has been derived that identifies the contributions associated with vortex stretching and strain self-amplification (Johnson, Phys. Rev. Lett., vol. 124, 2020 104501; J. Fluid Mech., vol. 922, 2021, A3) to the energy flux across scales. Here, we extend this methodology to general coupled advection–diffusion equations and, in particular, to homogeneous magnetohydrodynamic (MHD) turbulence. We show that several MHD subfluxes are related to each other by kinematic constraints akin to the Betchov relation in HD. Applied to data from direct numerical simulations, this decomposition allows for an identification of physical processes and for the quantification of their respective contributions to the energy cascade, as well as a quantitative assessment of their multi-scale nature through a further decomposition into single- and multi-scale terms. We find that vortex stretching is strongly depleted in MHD compared with HD, and the kinetic energy is transferred from large to small scales almost exclusively by the generation of regions of small-scale intense strain induced by the Lorentz force. In regions of large strain, current sheets are stretched by large-scale straining motion into regions of magnetic shear. This magnetic shear in turn drives extensional flows at smaller scales. Magnetic energy is transferred from large to small scales predominantly by the aforementioned current-sheet thinning in regions of high strain. The contributions from current-filament stretching – the analogue to vortex stretching – and from bending of magnetic field-lines into current filaments by vortical motion are both almost negligible, although the latter induces strong backscatter of magnetic energy. Consequences of these results for subgrid-scale turbulence modelling are discussed.

Physics of the plasma rotation driven by biasing in linear traps is analysed for two limiting cases. The first, relevant for traps with low effective viscosity, considers the line-tying effects to be responsible for the drive as well as for the dissipation of the angular momentum. Meanwhile, in long and thin traps with collisional plasma or developed turbulence, the radial transport of the angular momentum becomes its primary loss channel. The momentum flux goes into the scrape-off layer, which makes conditions there partially responsible for the achievable rotation limits.

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.

We present a detailed characterization of the structure and evolution of differentially rotating plasmas driven on the MAGPIE pulsed-power generator (1.4 MA peak current, 240 ns rise time). The experiments were designed to simulate physics relevant to accretion discs and jets on laboratory scales. A cylindrical aluminium wire array Z pinch enclosed by return posts with an overall azimuthal off-set angle was driven to produce ablation plasma flows that propagate inwards in a slightly off-radial trajectory, injecting mass, angular momentum and confining ram pressure to a rotating plasma column on the axis. However, the plasma is free to expand axially, forming a collimated, differentially rotating axial jet that propagates at ${\approx }100\,{\rm km}\,{\rm s}^{-1}$. The density profile of the jet corresponds to a dense shell surrounding a low-density core, which is consistent with the centrifugal barrier effect being sustained along the jet's propagation. We show analytically that, as the rotating plasma accretes mass, conservation of mass and momentum implies plasma radial growth scaling as $r \propto t^{1/3}$. As the characteristic moment of inertia increases, the rotation velocity is predicted to decrease and settle on a characteristic value ${\approx }20\,{\rm km}\,{\rm s}^{-1}$. We find that both predictions are in agreement with Thomson scattering and optical self-emission imaging measurements.

The paper presents experimental results from the SMOLA device, which was built in the Budker Institute of Nuclear Physics for the verification of the helical mirror confinement idea. This concept involves active control of axial losses from the confinement zone in an open magnetic trap through the use of multiple mirrors that move in the plasma frame of reference. The discussed experiments focused on determining the cumulative effect of a helical mirror system in combination with a short segment of a stronger magnetic field. Combination of these two methods of axial flow suppression results in higher efficiency compared with each method individually. Different combinations of the mirrors were tested. The most effective flow suppression was observed if the short mirror was placed between the confinement region and the helical mirror. In this configuration, an effective mirror ratio of $R_{{\rm eff}} = 32.6\pm 7.8$ was achieved, along with a more than three-fold increase in plasma density within the confinement region. The possibility of a cumulative effect of different types of magnetic mirrors offers a way to improve the confinement performance of the reactor-grade mirror confinement devices.

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.

This paper establishes a mean-field equation set and an energy theorem to provide a theoretical basis in view of the development of self-consistent, physics-based turbulent transport models for mean-field transport codes. A rigorous averaging procedure identifies the exact form of the perpendicular turbulent fluxes which are modelled by ad hoc diffusive terms in mean-field transport codes, next to other closure terms which are not commonly considered. Earlier work suggested that the turbulent $E\times B$ particle and heat fluxes, which are thus identified to be important closure terms, can be modelled to reasonable accuracy using the kinetic energy in the $E\times B$ velocity fluctuations ($k_{E}$). The related enstrophy led to further modelling improvements in an initial study, although further analysis is required. To support this modelling approach, transport equations are derived analytically for both quantities. In particular, an energy theorem is established in which the various source and sink terms of $k_{E}$ are shown to couple to mean-field and turbulent parallel kinetic energy, kinetic energy in the other perpendicular velocity components, the thermal energy and the magnetic energy. This provides expressions for the interchange, drift-wave and Reynolds stress terms amongst others. Note that most terms in these energy equations are in turn closure terms. It is suggested to evaluate these terms using reference data from detailed turbulence code simulations in future work.

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.

Satellite data analysis of a compressed gyro-scale current sheet prior to magnetic reconnection in the magnetotail shows that electrostatic lower hybrid waves localized to the region of a transverse ambipolar electric field at the centre of the current sheet are driven by $\boldsymbol{E} \times \boldsymbol{B}$ velocity shear and result from compression. The presence and location of shear-driven waves around the centre of the current sheet, where the magnetic field reverses and the density gradient is minimal, is consistent with our model. This is notable because the free energy source is the curvature of the electron $\boldsymbol{E} \times \boldsymbol{B}$ flow and not the density gradient. Laboratory experiments and particle-in-cell (PIC) simulations have shown that shear-driven lower hybrid fluctuations are capable of producing anomalous cross-field transport (viscosity) and resistivity, which can trigger magnetic reconnection. We estimate the terms in the generalized Ohm's Law directly from MMS data as the spacecraft cross a gyro-scale current sheet. Our analysis shows that the wave effects (resistivity, diffusion and viscosity) and pressure anisotropy effects are comparable. We also find that the quasi-static electric field gradient is correlated with a non-gyrotropic electron distribution function, which is consistent with our model. Furthermore, theoretical arguments suggest agyrotropy is an indicator of the possibility for magnetic reconnection to occur.

The properties of the edge ion-scale turbulence are studied using the beam emission spectroscopy (BES) diagnostic on MAST. Evidence of the formation of large-scale high-amplitude coherent structures, filamentary density blobs and holes, 2–4 cm inside the plasma separatrix is presented. Measurements of radial velocity and skewness of the density fluctuations indicate that density holes propagate radially inwards, with the skewness profile peaking at 7–10 cm inside the separatrix. Poloidal velocities of the density fluctuations measured using cross-correlation time delay estimation (CCTDE) are found to exhibit an intermittent behaviour. Zonal-flow analysis reveals the presence of poloidally symmetric coherent oscillations – low-frequency (LF) zonal flows and geodesic acoustic modes (GAM). Shearing rates of the observed zonal flows are found to be comparable to the turbulence decorrelation rate. The observed bursts in density-fluctuation power are followed by quiescent periods with a transient increase in the power of sheared flows. Three-wave interactions between broadband turbulence and a GAM are illustrated using the autobispectral technique. It is shown that the zonal flows and the density-fluctuation field are nonlinearly coupled and LF zonal flows mediate the energy transfer from high- to low-frequency density fluctuations.

The objective of the present paper is to investigate the constancy of the topological invariant, denoted the non-barotropic generalized cross-helicity in the case of non-ideal magnetohydrodynamics (MHD). Existing work considers only ideal barotropic MHD and ideal non-barotropic MHD. Here, we consider dissipative processes in the form of thermal conduction, finite electrical conductivity and viscosity and the effect of these processes on the cross-helicity conservation. An analytical approach has been adopted to obtain the mathematical expressions for the time derivative of the cross-helicity. Obtained results show that the generalized cross-helicity is not conserved in the non-ideal MHD limit and indicate which processes affect the helicity and which do not. Furthermore, we indicate the configurations in which this topological constant is conserved despite the dissipative processes. Some examples and applications are also given.

Charge exchange recombination spectroscopy (CXRS) measures the radial electric field in the pedestal by measuring the impurity density, temperature and flow. Combined outboard and inboard CXRS measurements allow poloidal variations that arise due to the poloidal variation of the magnetic field to be determined. At present, impurity neoclassical pedestal models avoid the complications of treating finite poloidal gyroradius effects by assuming the impurity charge number is large compared with the main ion charge number. These models are extended slightly by retaining the impurity radial pressure gradient to demonstrate that no substantial effect occurs due to impurity diamagnetic effects. More importantly, the neoclassical model is significantly extended to obtain a more comprehensive treatment of the main ions in the plateau and banana regimes. A parallel impurity momentum equation is derived that is consistent with previous results in the banana regime and reduces to the proper large aspect ratio form required in the plateau regime. The implications for interpreting the CXRS measurements are discussed by writing all results in terms of the gradient drive and poloidal flow.

Pellet injection is used for fuelling and controlling discharges in tokamaks, and it is foreseen in ITER. During pellet injection, a movement of the ablated material towards the low-field side (or outward major radius direction) occurs because of the inhomogeneity of the magnetic field. Due to the complexity of the theoretical models, computer codes developed to simulate the cross-field drift are computationally expensive. Here, we present a one-dimensional semi-analytical model for the radial displacement of ablated material after pellet injection, taking into account both the Alfvén and ohmic currents which shortcircuit the charge separation creating the drift. The model is suitable for rapid calculation of the radial drift displacement, and can be useful for e.g. modelling of disruption mitigation via pellet injection.

The paper presents experimental results from the SMOLA device on the testing of the helical mirror confinement hypothesis. Helical mirror confinement is the technique of an active control of axial plasma losses from a confinement zone by multiple magnetic mirrors that move along the axis in the reference frame of the plasma that experiences $\boldsymbol{E} \times \boldsymbol{B}$ rotation due to an applied radial electric field. Theory predicts that a helical mirror will provide an axial force that modifies the plasma flow and, simultaneously, density pinching to the axis. The force direction depends on the plasma rotation direction. Experimental data on the axial plasma losses at different direction of the magnetic mirror movement are presented. If the trapped ions move in the direction opposite to the direction of the axial losses, then the particle flux reduces in the broad range of the plasma density. The confinement improves with the increase of the fraction of the trapped particles (effective mirror ratio was up to $R_{{\rm eff}}=5.8\pm 1.4$). If the trapped ions move in the same direction as the axial losses, then the flux depends on density. At intermediate densities, the integral flux through the transport section rises compared to the plasma flowing through the straight magnetic field. The effective mirror ratio is lower and does not significantly depend on the fraction of the trapped particles (effective mirror ratio at intermediate density was $R_{{\rm eff}}=3.3\pm 0.8$).

Using highly resolved direct numerical simulations we examine the statistical properties of the local energy flux rate $\varPi _\ell (x)$ towards small scales for three isotropic turbulent magnetohydrodynamic flows, which differ in strength and structure of the magnetic field. We analyse the cascade process in both kinetic and magnetic energy, disentangling the different flux contributions to the overall energy dynamics. The results show that the probability distribution of the local energy flux develops long tails related to extreme events, similar to the hydrodynamic case. The different terms of the energy flux display different properties and show sensitivity to the type of the flow examined. We further examine the joint probability density function between the local energy flux and the gradients of the involved fields. The results point out a correlation with the magnetic field gradients, showing, however, a dispersion much stronger than what is observed in hydrodynamic flows. Finally, it is also shown that the local energy flux shows some dependence on the local amplitude of the magnetic field. The present results have implications for subgrid-scale models, which we discuss.

Based on the chemical reaction model proposed by Park, the ‘blackout’ of a reentry vehicle is studied in this paper. The temperature, pressure and electron density distribution characteristics around the reentry vehicle were simulated at various flight speeds and altitudes by USim. Subsequently, the scattering matrix method was used to study the transmission characteristics of terahertz waves in ‘blackout’. The simulation results show that the temperature around the aircraft is mainly affected by speed, the pressure is mainly affected by the altitude and electron density is affected by both of these factors. The calculation results show that the transmission characteristics of terahertz waves in plasma are mainly affected by electron density, while the effects of temperature and pressure cannot be ignored either.

The dispersion relation of a surface wave generated by a drifting plasma in an infinite duct surrounded by vacuum is derived non-relativistically by means of the Vlasov equation. The kinematic boundary condition imposed on the distribution function, the specular reflection conditions on the four sides of duct, can be satisfied by placing an infinite number of fictitious surface charge sheets spaced by the duct widths. The surface wave mode is specifically the transverse magnetic mode, often called the surface polariton, which propagates with phasor $\exp ({{\rm i}k_zz-{\rm i}\omega t})$. The method of placing appropriate fictitious surface charge sheets enables one to treat the surface waves in semi-infinite, slab and duct plasmas simultaneously on an equal footing, kinetically. The streaming effect manifests itself through the Doppler-shifted frequency and a correction-like term ${u^2}/{c^2}$, where u is the streaming velocity and c is the speed of light.

The Weibel instability is investigated theoretically and numerically under three scenarios: counterstreaming electron beams in background plasma, an electron–positron beam and an electron–proton beam in background plasma. These models occur widely in laboratory and astrophysical environments. The Weibel instability growth rates are determined numerically from the corresponding cold-fluid dispersion relations, which are confirmed with two-dimensional particle-in-cell simulations. The maximum growth rates for the counterstreaming beams in background plasma are an order of magnitude smaller than the maximum growth rates for the beams cases in the same range of density ratios and beam energies. The maximum growth rate for the electron–positron beam case is shown to be at most a factor $\sqrt {2}$ greater than the electron–proton beam case with similar dispersion behaviours. A non-monotonic relation is found between the maximum Weibel instability growth rates and the electron–positron beam energy, suggesting that increasing beam energies does not entail an increase in the Weibel instability growth rate.

Magnetic reconnection is explored on the Terrestrial Reconnection Experiment (TREX) for asymmetric inflow conditions and in a configuration where the absolute rate of reconnection is set by an external drive. Magnetic pileup enhances the upstream magnetic field of the high-density inflow, leading to an increased upstream Alfvén speed and helping to lower the normalized reconnection rate to values expected from theoretical consideration. In addition, a shock interface between the far upstream supersonic plasma inflow and the region of magnetic flux pileup is observed, important to the overall force balance of the system, thereby demonstrating the role of shock formation for configurations including a supersonically driven inflow. Despite the specialized geometry where a strong reconnection drive is applied from only one side of the reconnection layer, previous numerical and theoretical results remain robust and are shown to accurately predict the normalized rate of reconnection for the range of system sizes considered. This experimental rate of reconnection is dependent on system size, reaching values as high as 0.8 at the smallest normalized system size applied.

The evolution of the radial profile of the rotation of a weakly magnetized plasma column is investigated experimentally in a radio-frequency argon plasma at low pressure when a strong electron current is emitted by large emissive cathodes. Current injection from large emissive cathodes over a background plasma column (with a plasma density of a few $10^{18}\ \textrm {m}^{-3}$) is characterized. Radial scans of the ion velocity show that a continuous control of the rotation profile may be obtained using two spatial configurations for the locations of the emissive cathodes (either in the centre or at the edge of the plasma column). The rotation profile results from the electric drift velocity, damped by the drag exerted on the ions. We demonstrate that the evolution of the rotation profile with the injected current is then controlled by the modification of the plasma potential profile in the presence of strongly emissive cathodes.