The paper presents experimental results from the SMOLA device that is the first facility with a helical mirror section of the magnetic field. This device is built in the Budker Institute of Nuclear Physics for the verification of the helical mirror confinement idea that is the recently introduced technique of the active control of axial losses from a confinement zone. Theory predicts that with rotating plasma, a helical mirror will provide suppression of the axial plasma flow and, simultaneously, density pinching to the axis. Experiments demonstrated that plasma density at the exit from the transport section is suppressed with activation of the helical field, the effect is significant and highly reproducible. The most pronounced effect is observed on the plasma periphery, where the mirror ratio is the highest. The integral suppression ratio reaches 2–2.5 in the discussed experiments. Experimental results are compared with simplified theoretical estimates. The integral suppression ratio matches the simple theoretical estimates even if the transversal diffusion is neglected.

Impurity temperature screening is a favourable neoclassical phenomenon involving an outward radial flux of impurity ions from the core of fusion devices. Quasisymmetric magnetic fields lead to intrinsically ambipolar neoclassical fluxes that give rise to temperature screening for low enough $\unicode[STIX]{x1D702}^{-1}\equiv d\ln n/d\ln T$. In contrast, neoclassical fluxes in generic stellarators will depend on the radial electric field, which is predicted to be inward for ion-root plasmas, potentially leading to impurity accumulation. Here, we examine the impurity particle flux in a number of approximately quasisymmetric stellarator configurations and parameter regimes while varying the amount of symmetry breaking in the magnetic field. For the majority of this work, neoclassical fluxes have been obtained using the SFINCS drift-kinetic equation solver with electrostatic potential $\unicode[STIX]{x1D6F7}=\unicode[STIX]{x1D6F7}(r)$, where $r$ is a flux-surface label. Results indicate that achieving temperature screening is possible, but unlikely, at low-collisionality reactor-relevant conditions in the core. Thus, the small departures from symmetry in nominally quasisymmetric stellarators are large enough to significantly alter the neoclassical impurity transport. A further look at the magnitude of these fluxes when compared to a gyro-Bohm turbulence estimate suggests that neoclassical fluxes are small in configurations optimized for quasisymmetry when compared to turbulent fluxes.

We present a Monte Carlo collisional scheme that models single Compton scattering between leptons and photons in particle-in-cell codes. The numerical implementation of Compton scattering can deal with macro-particles of different weights and conserves momentum and energy in each collision. Our scheme is validated through two benchmarks for which exact analytical solutions exist: the inverse Compton spectra produced by an electron scattering with an isotropic photon gas and the photon–electron gas equilibrium described by the Kompaneets equation. It provides new opportunities for numerical investigation of plasma phenomena where a significant population of high-energy photons is present in the system.

In support of the ICRF experiments planned on the Wendelstein 7-X (W7-X) stellarator, i.e. fast ion generation, wall conditioning, target plasma production and heating, a first experimental study on plasma production has been made in the Uragan-2M (U-2M) stellarator using W7-X-like two-strap antenna. In all the experiments, antenna monopole phasing was used. The W7-X-like antenna operation with launched radiofrequency power of ~100 kW have been performed in helium (p = (4–14) × 10−2 Pa) with the vacuum vessel walls pre-loaded with hydrogen. Production of plasma with a density higher than 1012 cm−3 was observed near the first harmonic of the hydrogen cyclotron frequency. Operation at first hydrogen harmonic is feasible in W7-X future ICRF experiments.

Recently, the validity range of the approximations commonly used in neoclassical calculation has been reconsidered. One of the primary motivations behind this trend is observation of an impurity hole in LHD (Large Helical Device), i.e. the formation of an extremely hollow density profile of an impurity ion species, such as carbon $\text{C}^{6+}$, in the plasma core region where a negative radial electric field ($E_{r}$) is expected to exist. Recent studies have shown that the variation of electrostatic potential on the flux surface, $\unicode[STIX]{x1D6F7}_{1}$, has significant impact on neoclassical impurity transport. Nevertheless, the effect of $\unicode[STIX]{x1D6F7}_{1}$ has been studied with radially local codes and the necessity of global calculation has been suggested. Thus, we have extended a global neoclassical code, FORTEC-3D, to simulate impurity transport in an impurity hole plasma including $\unicode[STIX]{x1D6F7}_{1}$ globally. Independently of the $\unicode[STIX]{x1D6F7}_{1}$ effect, an electron root of the ambipolar condition for the impurity hole plasma has been found by global simulation. Hence, we have considered two different cases, each with a positive (global) and a negative (local) solution of the ambipolar condition, respectively. Our result provides another support that $\unicode[STIX]{x1D6F7}_{1}$ has non-negligible impact on impurity transport. However, for the ion-root case, the radial $\text{C}^{6+}$ flux is driven further inwardly by $\unicode[STIX]{x1D6F7}_{1}$. For the electron-root case, on the other hand, the radial particle $\text{C}^{6+}$ flux is outwardly enhanced by $\unicode[STIX]{x1D6F7}_{1}$. These results indicate that how $\unicode[STIX]{x1D6F7}_{1}$ affects the radial particle transport crucially depends on the profile of the ambipolar-$E_{r}$, which is found to be susceptible to $\unicode[STIX]{x1D6F7}_{1}$ itself and the global effects.

We present the elements required to construct two devices used in an undergraduate plasma physics laboratory. The materials and construction costs of the sources, the vacuum systems and probe drives and electrical circuits are presented in detail in the text and the first appendix. We also provide the software for probe motion and data acquisition as well as the electrical schematics for key components. Experiments which have been performed are listed and two (resonance cones and whistler waves) are described in greater detail. The machines are flexible and original research is possible.

Utilizing the similarity between the spinor representation of the Dirac and the Maxwell equations that has been recognized since the early days of relativistic quantum mechanics, a quantum lattice algorithm (QLA) representation of unitary collision-stream operators of Maxwell's equations is derived for both homogeneous and inhomogeneous media. A second-order accurate 4-spinor scheme is developed and tested successfully for two-dimensional (2-D) propagation of a Gaussian pulse in a uniform medium whereas for normal (1-D) incidence of an electromagnetic Gaussian wave packet onto a dielectric interface requires 8-component spinors because of the coupling between the two electromagnetic polarizations. In particular, the well-known phase change, field amplitudes and profile widths are recovered by the QLA asymptotic profiles without the imposition of electromagnetic boundary conditions at the interface. The QLA simulations yield the time-dependent electromagnetic fields as the wave packet enters and straddles the dielectric boundary. QLA involves unitary interleaved non-commuting collision and streaming operators that can be coded onto a quantum computer: the non-commutation being the very reason why one perturbatively recovers the Maxwell equations.

Statistical properties of the scrape-off layer plasma fluctuations are studied in ohmically heated plasmas in the Alcator C-Mod tokamak. For the first time, plasma fluctuations as well as parameters that describe the fluctuations are compared across measurements from a mirror Langmuir probe (MLP) and from gas-puff imaging (GPI) that sample the same plasma discharge. This comparison is complemented by an analysis of line emission time-series data, synthesized from the MLP electron density and temperature measurements. The fluctuations observed by the MLP and GPI typically display relative fluctuation amplitudes of order unity together with positively skewed and flattened probability density functions. Such data time series are well described by an established stochastic framework that models the data as a superposition of uncorrelated, two-sided exponential pulses. The most important parameter of the process is the intermittency parameter, $\gamma = {\tau _{d}} / {\tau _{w}}$, where ${\tau _{d}}$ denotes the duration time of a single pulse and ${\tau _{w}}$ gives the average waiting time between consecutive pulses. Here we show, using a new deconvolution method, that these parameters can be consistently estimated from different statistics of the data. We also show that the statistical properties of the data sampled by the MLP and GPI diagnostic are very similar. Finally, a synthetic GPI signal using only plasma parameters sampled by the MLP shows qualitatively different fluctuation statistics from the measured GPI signal.

The multiphase simulations are conducted with the kinetic-magnetohydrodynamics hybrid code MEGA to investigate the spatial and the velocity distributions of lost fast ions due to the Alfvén eigenmode (AE) bursts in the Large Helical Device plasmas. It is found that fast ions are lost along the divertor region with helical symmetry both before and during the AE burst except for the promptly lost particles. On the other hand, several peaks are present in the spatial distribution of lost fast ions along the divertor region. These peaks along the divertor region can be attributed to the deviation of the fast-ion orbits from the magnetic surfaces due to the grad-B and the curvature drifts. For comparison with the velocity distribution of lost fast ions measured by the fast-ion loss detector (FILD), the ‘numerical FILD’ which solves the Newton–Lorentz equation is constructed in the MEGA code. The velocity distribution of lost fast ions detected by the numerical FILD during AE burst is in good qualitative agreement with the experimental FILD measurements. During the AE burst, fast ions with high energy (100–180 keV) are detected by the numerical FILD, while co-going fast ions lost to the divertor region are the particles with energy lower than 50 keV.

Here, we propose a thermodynamic model for dusty plasma, where the dust is confined in a small volume within a large plasma background by external fields. In this model, the parameters of dust, e.g. Helmholtz energy, pressure, entropy and enthalpy, etc. can be calculated for given dust density and temperature. The model is solved analytically in the mean field (gaseous) limit and various processes associated with the gaseous phase of dust, e.g. adiabatic/isothermal/constant internal energy expansion/compression, specific heat, free expansion within the plasma background, and the dispersion of novel acoustic waves are studied. Some predictions of the model, e.g. electrostatic pressure of the dust and the isothermal equation of state, have been earlier verified in experiments and numerical simulations. The model is compared with an earlier thermodynamic model of dusty plasma proposed by Hamaguchi and Farouki.