Tokamak start-up is strongly dependent on the state of the initial plasma formed during plasma breakdown. To acquire a better understanding of the process and to estimate the influence of the impurity of beryllium on the ohmic heating tokamak start-up process, one-dimensional particle-in-cell coupled with a Monte Carlo collision method has been developed. The main aim is to investigate the plasma performance under various amounts of beryllium with different discharge parameters. Tokamak breakdown with the impurity of beryllium in the ohmic heating strategy has been simulated. The simulation results show that with the impurity of beryllium, the increase of plasma density is suppressed compared with the case without beryllium. The breakdown time is delayed by the impurity. Moreover, the successful breakdown has a much higher requirement on discharge parameters with a low electric field operational scenario, since in the low electric field discharge the influence of beryllium impurity is greater. As the plasma density increases, the effect of beryllium impurity on plasma becomes more critical. It indicates that impurities cannot be neglected in the high plasma density.

We demonstrate a fast adjoint-based method to optimise tokamak and stellarator equilibria against a pressure-driven instability known as the infinite-$n$ ideal ballooning mode. We present three finite-$\beta$ (the ratio of thermal to magnetic pressure) equilibria: one tokamak equilibrium and two stellarator equilibria that are unstable against the ballooning mode. Using the self-adjoint property of ideal magnetohydrodynamics, we construct a technique to rapidly calculate the change in the eigenvalue, a measure of ideal ballooning instability. Using the SIMSOPT optimisation framework, we then implement our fast adjoint gradient-based optimiser to minimise the eigenvalue and find stable equilibria for each of the three originally unstable equilibria.

The magnetohydrodynamic stability of plasmas with an anisotropic pressure component is analysed for a low magnetic field configuration of the large helical device. Magnetic equilibria are calculated by the anisotropic Neumann inverse moments equilibrium code, an extension of the three-dimensional variational moments equilibrium code. A modified version of the bi-Maxwellian is used to model the anisotropic particle velocity distribution. Magnetohydrodynamic stability calculations for the $n=1$ mode family are carried out by TERPSICHORE, which has been expanded by the Kruskal–Oberman energy principle. For on-axis particle deposition, the growth rate and plasma displacement show that the parallel dominant plasmas are significantly more stable than isotropic or perpendicular dominant plasmas. For off-axis particle deposition, the growth rate and the Mercier criterion in the peripheral region $\rho =0.9$, show that low field (LF) deposition perpendicular dominant plasmas are most unstable. For the most realistic off-axis deposition profile, it is found that parallel dominant plasmas are most stable for LF deposition, while perpendicular dominant plasmas are most stable for high field deposition. We conclude that, under low magnetic field conditions in the large helical device, tangential neutral beam injection heating has a stabilising influence on the plasma.

It is shown that the extensive literature on leaky modes in coronal magnetic flux tubes should be revised since those modes have no physical meaning.

This paper presents the calculation of the bounce-averaged drift of trapped particles in a near-axis framework for axisymmetric and quasisymmetric magnetic fields that possess up-down and stellarator symmetry, respectively. This analytic consideration provides important insight on the dependence of the bounce-averaged drift on the geometry and stability properties of the field. In particular, we show that although the maximum-$\mathcal {J}$ property is unattainable in quasisymmetric stellarators, one may approach it through increased plasma $\beta$ and triangular shaping, albeit going through a reduced precession scenario with potentially higher particle losses. The description of trapped particles allows us to calculate the available energy of trapped electrons analytically in two asymptotic regimes, providing insight into the behaviour of this measure of turbulence. It is shown that the available energy is intimately related to magnetohydrodynamics (MHD) stability, providing a potential synergy between this measure of gyrokinetic turbulence and MHD stability.

Available energy (Æ), which quantifies the maximum amount of thermal energy that may be liberated and converted into instabilities and turbulence, has shown to be a useful metric for predicting saturated energy fluxes in trapped-electron-mode-driven turbulence. Here, we calculate and investigate the Æ in the analytical tokamak equilibria introduced by Miller et al. (Phys. Plasmas, vol. 5, issue, 4, 1998, pp. 973–978). The Æ of trapped electrons reproduces various trends also observed in experiments; negative shear, increasing Shafranov shift, vertical elongation and negative triangularity can all be stabilising, as indicated by a reduction in Æ, although it is strongly dependent on the chosen equilibrium. Comparing Æ with saturated energy flux estimates from the TGLF (trapped gyro-Landau fluid) model, we find fairly good correspondence, showcasing that Æ can be useful to predict trends. We go on to investigate Æ and find that negative triangularity is especially beneficial in vertically elongated configurations with positive shear or low gradients. Furthermore, we extract a gradient-threshold-like quantity from Æ and find that it behaves similarly to gyrokinetic gradient thresholds: it tends to increase linearly with magnetic shear, and negative triangularity leads to an especially high threshold. We next optimise the device geometry for minimal Æ and find that the optimum is strongly dependent on equilibrium parameters, for example, magnetic shear or pressure gradient. Investigating the competing effects of increasing the density gradient, the pressure gradient, and decreasing the shear, we find regimes that have steep gradients yet low Æ, and that such a regime is inaccessible in negative-triangularity tokamaks.

This paper documents the design and construction of MUSE, the world's first permanent magnet (PM) stellarator and the first quasi-axisymmetric experiment. The purpose of MUSE is to develop and assess a new way of building optimised stellarators that uses simple planar coils PMs. Our PM optimisation algorithm consists of initialising a geometry to pack dipoles densely, running the FAMUS code to minimise surface field error subject to PM constraints and applying discrete jumps to reach a physically realisable solution. FAMUS treats the PM system as a set of ideal point dipoles. From there we construct finite-volume magnet towers to be housed in 3D-printed PM holders. We describe the design of the PM holders, which were validated by laser metrology. We analyse the effects of finite permeability, sensitivity to perturbations and magnetostatic forces. An exact analytic formula for the magnetic field from a finite-volume PM tower is presented to compute PM–PM forces and stress on the PM holder. Stellarator construction is complete and experiments are underway.

In this paper, we review recent and ongoing work by our group on numerical simulations of relativistic jets. Relativistic outflows in astrophysics are related to dilute, high energy plasmas, with physical conditions out of the reach of current laboratory capabilities. Simulations are thus imperative for the study of these objects. We present a number of such scenarios that have been studied by our group at the Universitat de València. In particular, we have focused on the evolution of extragalactic outflows through galactic and intergalactic environments, deceleration by interaction with stars or clouds or the propagation of jets in X-ray binaries and interaction with stellar winds from massive companions. All also share their role as particle acceleration sites and production of non-thermal radiation throughout the electromagnetic spectrum. Therefore, our work is not only aimed at understanding the impact of outflows on their environments and thus their role in galaxy and cluster evolution, but also the nature and capabilities of these sites as generators of high- and very-high-energy radiation and cosmic rays.

Methods for reducing the radius, temperature and space charge of a non-neutral plasma are usually reported for conditions which approximate an ideal Penning Malmberg trap. Here, we show that (i) similar methods are still effective under surprisingly adverse circumstances: we perform strong drive regime (SDR) compression and SDREVC in a strong magnetic mirror field using only 3 out of 4 rotating wall petals. In addition, we demonstrate (ii) an alternative to SDREVC, using e-kick instead of evaporative cooling (EVC) and (iii) an upper limit for how much plasma can be cooled to $T<20\ \mathrm {K}$ using EVC. This limit depends on the space charge, not on the number of particles or the plasma density.

We study the stability of plasma wake wave and the properties of density-downramp injection in an electron-driven plasma accelerator. In this accelerator type, a short high-current electron bunch (generated by a conventional accelerator or a laser–wakefield acceleration stage) drives a strongly nonlinear plasma wake wave (blowout), and accelerated electrons are injected into it using a sharp density transition which leads to the elongation of the wake. The accelerating structure remains highly stable until the moment some electrons of the driver reach almost zero energy, which corresponds to the best interaction length for optimal driver-to-plasma energy transfer efficiency. For a particular driver, this efficiency can be optimised by choosing appropriate plasma density. Studying the dependence of the current of the injected bunch on driver and plasma parameters, we show that it does not depend on the density downramp length as long as the condition for trapping is satisfied. Most importantly, we find that the current of the injected bunch primarily depends on just one parameter which combines both the properties of the driver (its current and duration) and the plasma density.