The dynamical behaviour of weakly nonlinear, low-frequency sound waves are investigated in a plasma composed of only positive and negative ions incorporating the effects of a weak external uniform magnetic field. In the plasma model the mass (temperature) of the positive ions is smaller (larger) than that of the negative ions. The dynamics of the nonlinear wave is shown to be governed by a novel nonlinear equation. The stationary plane wave (analytical and numerical) nonlinear analysis on the basis of experimental parameters reveals that the nonlinear wave does have quasi-periodic and chaotic solutions. The Poincarè return map analysis confirms these observed complex structures.

Numerical simulations of Alfvén modes driven by energetic particles are performed with the gyrokinetic (GK) global particle-in-cell code ORB5. A reversed shear equilibrium magnetic field is adopted. A simplified configuration with circular flux surfaces and large aspect ratio is considered. The nonlinear saturation of beta-induced Alfvén eigenmodes (BAE) is investigated. The roles of the wave–particle nonlinearity of the different species, i.e. thermal ions, electrons and energetic ions are described, in particular for their role in the saturation of the BAE and the generation of zonal flows. The nonlinear redistribution of the electron population is found to be important in increasing the BAE saturation level and the zonal flow amplitude.

The problem of plasma expansion into a vacuum is revisited with the addition of a finite boundary condition; an electrically insulated surface. As plasma expands towards a charge-accumulating surface, the leading electron cloud charges the surface negatively, which in turn repels electrons and attracts ions. This plasma–surface interaction is shown to result in a feedback process which accelerates the plasma expansion. In addition, we examine the decrease in (negative) surface potential and associated near-surface electron density. To investigate this plasma coupling with an electrically floating surface, we develop an analytic model including four neighbouring plasma regions: (i) undisturbed plasma, (ii) quasi-neutral self-similar expansion, (iii) ion front boundary layer and (iv) electron cloud. A key innovation in our approach is a self-contained analytic approximation of the ion front boundary layer, providing a spatially continuous electric field model for the early phase of bounded plasma expansion.

The excitation of low frequency dust acoustic (or dust density) waves in a dusty plasma can be driven by the flow of ions relative to dust. We consider the nonlinear development of the ion–dust streaming instability in a highly collisional plasma, where the ion and dust collision frequencies are a significant fraction of their corresponding plasma frequencies. This collisional parameter regime may be relevant to dusty plasma experiments under microgravity or ground-based conditions with high gas pressure. One-dimensional particle-in-cell simulations are presented, which take into account collisions of ions and dust with neutrals, and a background electric field that drives the ion flow. Ion flow speeds of the order of a few times thermal are considered. Waveforms of the dust density are found to have broad troughs and sharp crests in the nonlinear phase. The results are compared with the nonlinear development of the ion–dust streaming instability in a plasma with low collisionality.

A Hamiltonian six-field gyrofluid model is constructed, based on closure relations derived from the so-called ‘quasi-static’ gyrokinetic linear theory where the fields are assumed to propagate with a parallel phase velocity much smaller than the parallel particle thermal velocities. The main properties captured by this model, primarily aimed at exploring fundamental problems of interest for space plasmas such as the solar wind, are its ability to provide a reasonable agreement with kinetic theory for linear low-frequency modes, and at the same time to ensure a Hamiltonian structure in the absence of explicit dissipation. The model accounts for equilibrium temperature anisotropy, ion and electron finite Larmor radius corrections, electron inertia, magnetic fluctuations along the direction of a strong guide field and parallel Landau damping, introduced through a Landau-fluid modelling of the parallel heat transfers for both gyrocentre species. Remarkably, the quasi-static closure leads to exact and simple expressions for the nonlinear terms involving gyroaveraged electromagnetic fields and potentials. One of the consequences is that a rather natural identification of the Hamiltonian structure of the model becomes possible when Landau damping is neglected. A slight variant of the model consists of a four-field Hamiltonian reduction of the original six-field model, which is also used for the subsequent linear analysis. In the latter, the dispersion relations of kinetic Alfvén waves and the firehose instability are shown to be correctly reproduced, relatively far in the sub-ion range (depending on the plasma parameters), while the spectral range where the slow-wave dispersion relation and the field-swelling instabilities are precisely described is less extended. This loss of accuracy originates from the breaking of the condition of small phase velocity, relative to the parallel thermal velocity of the electrons (for kinetic Alfvén waves and firehose instability) or of the ions (in the case of the field-swelling instabilities).

The SPARC tokamak is a critical next step towards commercial fusion energy. SPARC is designed as a high-field ($B_0 = 12.2$ T), compact ($R_0 = 1.85$ m, $a = 0.57$ m), superconducting, D-T tokamak with the goal of producing fusion gain $Q>2$ from a magnetically confined fusion plasma for the first time. Currently under design, SPARC will continue the high-field path of the Alcator series of tokamaks, utilizing new magnets based on rare earth barium copper oxide high-temperature superconductors to achieve high performance in a compact device. The goal of $Q>2$ is achievable with conservative physics assumptions ($H_{98,y2} = 0.7$) and, with the nominal assumption of $H_{98,y2} = 1$, SPARC is projected to attain $Q \approx 11$ and $P_{\textrm {fusion}} \approx 140$ MW. SPARC will therefore constitute a unique platform for burning plasma physics research with high density ($\langle n_{e} \rangle \approx 3 \times 10^{20}\ \textrm {m}^{-3}$), high temperature ($\langle T_e \rangle \approx 7$ keV) and high power density ($P_{\textrm {fusion}}/V_{\textrm {plasma}} \approx 7\ \textrm {MW}\,\textrm {m}^{-3}$) relevant to fusion power plants. SPARC's place in the path to commercial fusion energy, its parameters and the current status of SPARC design work are presented. This work also describes the basis for global performance projections and summarizes some of the physics analysis that is presented in greater detail in the companion articles of this collection.

We present an energy-conserving discontinuous Galerkin scheme for the full-$f$ electromagnetic gyrokinetic system in the long-wavelength limit. We use the symplectic formulation and solve directly for $\unicode[STIX]{x2202}A_{\Vert }/\unicode[STIX]{x2202}t$, the inductive component of the parallel electric field, using a generalized Ohm’s law derived directly from the gyrokinetic equation. Linear benchmarks are performed to verify the implementation and show that the scheme avoids the Ampère cancellation problem. We perform a nonlinear electromagnetic simulation in a helical open-field-line system as a rough model of the tokamak scrape-off layer using parameters from the National Spherical Torus Experiment (NSTX). This is the first published nonlinear electromagnetic gyrokinetic simulation on open field lines. Comparisons are made to a corresponding electrostatic simulation.

A gyrokinetic model is presented that can properly describe large and small amplitude electromagnetic fluctuations occurring on scale lengths ranging from the electron Larmor radius to the equilibrium perpendicular pressure gradient scale length, and the arbitrarily large deviations from thermal equilibrium that are present in the plasma periphery of tokamak devices. The formulation of the gyrokinetic model is based on a second-order accurate description of the single charged particle dynamics, derived from Lie perturbation theory, where the fast particle gyromotion is decoupled from the slow drifts assuming that the ratio of the ion sound Larmor radius to the perpendicular equilibrium pressure scale length is small. The collective behaviour of the plasma is obtained by a gyrokinetic Boltzmann equation that describes the evolution of the gyroaveraged distribution function. The collisional effects are included by a nonlinear gyrokinetic Dougherty collision operator. The gyrokinetic model is then developed into a set of coupled fluid equations referred to as the gyrokinetic moment hierarchy. To obtain this hierarchy, the gyroaveraged distribution function is expanded onto a Hermite–Laguerre velocity-space polynomial basis. Then, the gyrokinetic equation is projected onto the same basis yielding the spatial and temporal evolution of the Hermite–Laguerre expansion coefficients. A closed set of fluid equations for the lowest-order coefficients is presented. The Hermite–Laguerre projection is performed accurately at arbitrary perpendicular wavenumber values. Finally, the self-consistent evolution of the electromagnetic fields is described by a set of gyrokinetic Maxwell equations derived from a variational principle where the velocity integrals are explicitly evaluated.

In recent years, several optimal dynamos have been discovered. They minimize the magnetic energy dissipation or, equivalently, maximize the growth rate at a fixed magnetic Reynolds number. In the optimal dynamo of Willis (Phys. Rev. Lett., vol. 109, 2012, 251101), we find mean-field dynamo action for planar averages. One component of the magnetic field grows exponentially while the other decays in an oscillatory fashion near onset. This behaviour is different from that of an $\unicode[STIX]{x1D6FC}^{2}$ dynamo, where the two non-vanishing components of the planar averages are coupled and have the same growth rate. For the Willis dynamo, we find that the mean field is excited by a negative turbulent magnetic diffusivity, which has a non-uniform spatial profile near onset. The temporal oscillations in the decaying component are caused by the corresponding component of the diffusivity tensor being complex when the mean field is decaying and, in this way, time dependent. The growing mean field can be modelled by a negative magnetic diffusivity combined with a positive magnetic hyperdiffusivity. In two other classes of optimal dynamos of Chen et al. (J. Fluid Mech., vol. 783, 2015, pp. 23–45), we find, to some extent, similar mean-field dynamo actions. When the magnetic boundary conditions are mixed, the two components of the planar averaged field grow at different rates when the dynamo is 15 % supercritical. When the mean magnetic field satisfies homogeneous boundary conditions (where the magnetic field is tangential to the boundary), mean-field dynamo action is found for one-dimensional averages, but not for planar averages. Despite having different spatial profiles, both dynamos show negative turbulent magnetic diffusivities. Our finding suggests that negative turbulent magnetic diffusivities may support a broader class of dynamos than previously thought, including these three optimal dynamos.

SPARC is designed to be a high-field, medium-size tokamak aimed at achieving net energy gain with ion cyclotron range-of-frequencies (ICRF) as its primary auxiliary heating mechanism. Empirical predictions with conservative physics indicate that SPARC baseline plasmas would reach $Q\approx 11$, which is well above its mission objective of $Q>2$. To build confidence that SPARC will be successful, physics-based integrated modelling has also been performed. The TRANSP code coupled with the theory-based trapped gyro-Landau fluid (TGLF) turbulence model and EPED predictions for pedestal stability find that $Q\approx 9$ is attainable in standard H-mode operation and confirms $Q > 2$ operation is feasible even with adverse assumptions. In this analysis, ion cyclotron waves are simulated with the full wave TORIC code and alpha heating is modelled with the Monte–Carlo fast ion NUBEAM module. Detailed analysis of expected turbulence regimes with linear and nonlinear CGYRO simulations is also presented, demonstrating that profile predictions with the TGLF reduced model are in reasonable agreement.

The beam–plasma instability can be addressed as a reduced model in several contexts of plasma physics, from space to fusion plasma. In this paper, we review and refine some nonlinear features of this model. Specifically, by analysing the dependence of the nonlinear velocity spread as a function of the linear growth rate, we discuss the effective size of the resonance in view of its role in the spectral overlap at saturation. The relevance of this characterization relies on the necessity of a quantitative determination of the overlap degree to discriminate among different transport regimes of the self-consistent dynamics. The analysis is enriched with a study of the phase-space dynamics by means of the Lagrangian coherent structure technique, in order to define the transport barriers of the system describing the relevant features of the overlap process. Finally, we discuss relevant features related to the mode saturation levels.

We discuss a family of Vlasov–Maxwell equilibrium distribution functions for current sheet equilibria that are intermediate cases between the Harris sheet and the force-free (or modified) Harris sheet. These equilibrium distribution functions have potential applications to space and astrophysical plasmas. The existence of these distribution functions had been briefly discussed by Harrison & Neukirch (Phys. Rev. Lett., vol. 102, (2009a), 135003), but here it is shown that their approach runs into problems in the limit where the guide field goes to zero. The nature of this problem will be discussed and an alternative approach will be suggested that avoids the problem. This is achieved by considering a slight variation of the magnetic field profile, which allows a smooth transition between the Harris and force-free Harris sheet cases.

The formation and the dynamics of coherent magnetic field structures in the context of laser plasma interaction has attracted considerable attention. In the literature the formation of these structures has, however, mostly been reported in the wake of a laser pulse propagating in an underdense plasma medium (Bulanov et al., Phys. Rev. Lett., vol. 76, 1996, pp. 3562–3565; Nakamura & Mima Phys. Rev. Lett., vol. 100, 2008, 205006; Bulanov et al., Plasma Phys. Rep., vol. 31, no. 5, 2005, pp. 369–381; Naumova et al., Phys. Plasmas, vol. 8, no. 9, 2001, pp. 4149–4155; Nakamura et al., Phys. Rev. Lett., vol. 105, no. 13, 2010, 135002). The study here focuses on the formation of coherent structures by an intense laser pulse when it interacts with an overdense plasma medium. The laser in this case gets reflected and partially dumps its energy to the lighter electrons species. Particle-in-cell simulation studies have been carried out in two dimensions to show that the energetic electrons (generated at the critical layer and having relativistic energies), together with the background plasma electrons often self-organize to form distinct electron current vortices. These electron vortices have associated magnetic fields with monopolar or dipolar symmetries depending on the rotation profile of the electron current. The formation, stability and dynamics of these structures in the context of overdense plasma is of special importance as they provide a possibility of energy transport into those regions of plasma which are inaccessible to lasers. For such applications, higher energy content (involvement of relativistic electrons in their formation) of these structures is desirable. It is shown that their salient propagation characteristics even at relativistic energies follow the rules of electron magnetohydrodynamics (EMHD) (Isichenko & Marnachev, Sov. Phys. JETP, vol. 66, 1987, p. 702; Biskamp et al., Phys. Rev. Lett., vol. 76, 1996, p. 1264) (Generalized - EMHD Yadav et al., Phys. Plasmas, vol. 15, no. 6, 2008, 062308; Yadav et al., Phys. Plasmas, vol. 16, no. 4, 2009, 040701) for homogeneous (inhomogeneous) plasma medium, respectively.

This paper presents the results of studies of plasma flow parameters on the PF-1000U facility. A distinctive feature of this facility is the ability to create profiled initial gas distributions using gas puffs. In the experiments described, a combined system for filling the vacuum chamber with a working gas was used, in which an additional injection of various gases (deuterium, helium, neon and their mixtures) into the axial region of the chamber prefilled with deuterium was performed using a pulse valve. Thus, both the pinching processes and, accordingly, the generation of axial plasma flows and the conditions of their propagation in the background gas of the facility chamber were affected. Regimes with the generation of compact stable plasma formations propagating over long distances were found. The results obtained can be used in laboratory modelling of astrophysical jets from young stellar objects.

In order to inform core performance projections and divertor design, the baseline SPARC tokamak plasma discharge is evaluated for its expected H-mode access, pedestal pressure and edge-localized mode (ELM) characteristics. A clear window for H-mode access is predicted for full field DT plasmas, with the available 25 MW of design auxiliary power. Additional alpha heating is likely needed for H-mode sustainment. Pressure pedestal predictions in the developed H-mode are surveyed using the EPED model. The projected SPARC pedestal would be limited dominantly by peeling modes and may achieve pressures in excess of 0.3 MPa at a density of approximately 3 × 1020 m−3. High pedestal pressure is partially enabled by strong equilibrium shaping, which has been increased as part of recent design iterations. Edge-localized modes (ELMs) with >1 MJ of energy are projected, and approaches for reducing the ELM size, and thus the peak energy fluence to divertor surfaces, are under consideration. The high pedestal predicted for SPARC provides ample margin to satisfy its high fusion gain (Q) mission, so that even if ELM mitigation techniques result in a 2× reduction of the pedestal pressure, Q > 2 is still predicted.

The effects of the electron energy distribution function (EEDF) on the structure of a dusty plasma sheath are investigated. Here, it is assumed that the electrons obey a Druyvesteyn-type distribution with a parameter $x$ controlling the shape of the EEDF. The Druyvesteyn-like distribution tends to a Maxwellian distribution as $x$ varies from 2 to 1. Using the orbital motion limited theory, the incident electron current on the dust is evaluated for a given $x$. The results of numerical simulations are compared with those of a Maxwellian distribution. It was found that the sheath dynamics depends strongly on the magnitude of $x$. The sheath thickness increases monotonically with increasing $x$. However, the absolute dust charge decreases and, as a result, the accelerating ion drag force is weakened and thus the dust number density is enhanced. For a plasma with a Druyvesteyn-like distribution, the Bohm speed is a function of $x$ and increases with increasing $x$.

The incompressibility constraint for fluid flow was imposed by Lagrange in the so-called Lagrangian variable description using his method of multipliers in the Lagrangian (variational) formulation. An alternative is the imposition of incompressibility in the Eulerian variable description by a generalization of Dirac’s constraint method using noncanonical Poisson brackets. Here it is shown how to impose the incompressibility constraint using Dirac’s method in terms of both the canonical Poisson brackets in the Lagrangian variable description and the noncanonical Poisson brackets in the Eulerian description, allowing for the advection of density. Both cases give the dynamics of infinite-dimensional geodesic flow on the group of volume preserving diffeomorphisms and explicit expressions for this dynamics in terms of the constraints and original variables is given. Because Lagrangian and Eulerian conservation laws are not identical, comparison of the various methods is made.

Massive material injection has been proposed as a way to mitigate the formation of a beam of relativistic runaway electrons that may result from a disruption in tokamak plasmas. In this paper we analyse runaway generation observed in eleven ASDEX Upgrade discharges where disruption was triggered using massive gas injection. We present numerical simulations in scenarios characteristic of on-axis plasma conditions, constrained by experimental observations, using a description of the runaway dynamics with a self-consistent electric field and temperature evolution in two-dimensional momentum space and zero-dimensional real space. We describe the evolution of the electron distribution function during the disruption, and show that the runaway seed generation is dominated by hot-tail in all of the simulated discharges. We reproduce the observed dependence of the current dissipation rate on the amount of injected argon during the runaway plateau phase. Our simulations also indicate that above a threshold amount of injected argon, the current density after the current quench depends strongly on the argon densities. This trend is not observed in the experiments, which suggests that effects not captured by zero-dimensional kinetic modelling – such as runaway seed transport – are also important.

The head-on collision between two dust-acoustic solitary waves in a non-magnetized, collisionless and strongly coupled dust plasma has been studied. The application scope of the analytical solution of the head-on collision is given in the present paper by using the particle-in-cell simulation method. It is noted that the analytical results are valid if the amplitudes of both of the solitary waves are small enough. The effects of the coupling parameters on both the head-on collision and the waveform are also studied in the present paper.

Owing to its high magnetic field, high power, and compact size, the SPARC experiment will operate with divertor conditions at or above those expected in reactor-class tokamaks. Power exhaust at this scale remains one of the key challenges for practical fusion energy. Based on empirical scalings, the peak unmitigated divertor parallel heat flux is projected to be greater than 10 GW m−2. This is nearly an order of magnitude higher than has been demonstrated to date. Furthermore, the divertor parallel Edge-Localized Mode (ELM) energy fluence projections (~11–34 MJ m−2) are comparable with those for ITER. However, the relatively short pulse length (~25 s pulse, with a ~10 s flat top) provides the opportunity to consider mitigation schemes unsuited to long-pulse devices including ITER and reactors. The baseline scenario for SPARC employs a ~1 Hz strike point sweep to spread the heat flux over a large divertor target surface area to keep tile surface temperatures within tolerable levels without the use of active divertor cooling systems. In addition, SPARC operation presents a unique opportunity to study divertor heat exhaust mitigation at reactor-level plasma densities and power fluxes. Not only will SPARC test the limits of current experimental scalings and serve for benchmarking theoretical models in reactor regimes, it is also being designed to enable the assessment of long-legged and X-point target advanced divertor magnetic configurations. Experimental results from SPARC will be crucial to reducing risk for a fusion pilot plant divertor design.