No Arabic abstract
Recent tokamak experiments employing off-axis, non-inductive current drive have found that a large central current hole can be produced. The current density is measured to be approximately zero in this region, though in principle there was sufficient current drive power for the central current density to have gone significantly negative. Recent papers have used a large aspect-ratio expansion to show that normal MHD equilibria (with axisymmetric nested flux surfaces, non-singular fields, and monotonic peaked pressure profiles) can not exist with negative central current. We extend that proof here to arbitrary aspect ratio, using a variant of the virial theorem to derive a relatively simple integral constraint on the equilibrium. However, this constraint does not, by itself, exclude equilibria with non-nested flux surfaces, or equilibria with singular fields and/or hollow pressure profiles that may be spontaneously generated.
The effect of small deviations from a Maxwellian equilibrium on turbulent momentum transport in tokamak plasmas is considered. These non-Maxwellian features, arising from diamagnetic effects, introduce a strong dependence of the radial flux of co-current toroidal angular momentum on collisionality: As the plasma goes from nearly collisionless to weakly collisional, the flux reverses direction from radially inward to outward. This indicates a collisionality-dependent transition from peaked to hollow rotation profiles, consistent with experimental observations of intrinsic rotation.
The Stepped Pressure Equilibrium Code (SPEC) [Hudson et al., Phys. Plasmas 19, 112502 (2012)] has been successful in the construction of equilibria in 3D configurations that contain a mixture of flux surfaces, islands and chaotic magnetic field lines. In this model, the plasma is sliced into sub-volumes separated by ideal interfaces, and in each volume the magnetic field is a Beltrami field. In the cases where the system is far from possessing a continuous symmetry, such as in stellarators, the existence of solutions to a stepped-pressure equilibrium with given constraints, such as a multi-region relaxed MHD minimum energy state, is not guaranteed but is often taken for granted. Using SPEC, we have studied two different scenarios in which a solution fails to exist in a slab with analytic boundary perturbations. We found that with a large boundary perturbation, a certain interface becomes fractal, corresponding to the break up of a Kolmogorov-Arnold-Moser (KAM) surface. Moreover, an interface can only support a maximum pressure jump while a solution of the magnetic field consistent with the force balance condition can be found. An interface closer to break-up can support a smaller pressure jump. We discovered that the pressure jump can push the interface closer to being non-smooth through force balance, thus significantly decreasing the maximum pressure it can support. Our work shows that a convergence study must be performed on a SPEC equilibrium with interfaces close to break-up. These results may also provide insights into the choice of interfaces and have applications in finding out the maximum pressure a machine can support.
Neutral beam injection or ion cyclotron resonance heating induces pressure anisotropy. The axisymmetric plasma equilibrium code HELENA has been upgraded to include anisotropy and toroidal flow. With both analytical and numerical methods, we have studied the determinant factors in anisotropic equilibria and their impact on flux surfaces, magnetic axis shift, the displacement of pressures and density contours from flux surface. With $p_parallel/p_perp approx 1.5$, $p_perp$ can vary 20% on $s=0.5$ flux surface, in a MAST like equilibrium. We have also re-evaluated the widely applied approximation to anisotropy in which $p^*=(p_parallel + p_perp)/2$, the average of parallel and perpendicular pressure, is taken as the approximate isotropic pressure. We find the reconstructions of the same MAST discharge with $p_parallel/p_perp approx 1.25$, using isotropic and anisotropic model respectively, to have a 3% difference in toroidal field but a 66% difference in poloidal current.
A new force balance model for the EFIT magnetohydrodynamic equilibrium technique for tokamaks is presented which includes the full toroidal flow and anisotropy changes to the Grad-Shafranov equation. The free functions are poloidal flux functions and all non-linear contributions to the toroidal current density are treated iteratively. The parallel heat flow approximation chosen for the model is that parallel temperature is a flux function and that both parallel and perpendicular pressures may be described using parallel and perpendicular temperatures. This choice for the fluid thermodynamics has been shown elsewhere to be the same as a guiding centre kinetic solution of the same problem under the same assumptions. The model reduces identically to the static and isotropic Grad-Shafranov equation in the appropriate limit as different flux functions are set to zero. An analytical solution based on a modified Soloviev solution for non-zero toroidal flow and anisotropy is also presented. The force balance model has been demonstrated in the code EFIT TENSOR, a branch of the existing code EFIT++. Benchmark results for EFIT TENSOR are presented and the more complicated force balance model is found to converge to force balance similarly to the usual EFIT model and with comparable speed.
The NASA Magnetospheric Multiscale mission has made in situ diffusion region and kinetic-scale resolution measurements of asymmetric magnetic reconnection for the first time, in the Earths magnetopause. The principal theoretical tool currently used to model collisionless asymmetric reconnection is particle-in-cell simulations. Many particle-in-cell simulations of asymmetric collisionless reconnection start from an asymmetric Harris-type magnetic field but with distribution functions that are not exact equilibrium solutions of the Vlasov equation. We present new and exact equilibrium solutions of the Vlasov-Maxwell system that are self-consistent with one-dimensional asymmetric current sheets, with an asymmetric Harris-type magnetic field profile, plus a constant nonzero guide field. The distribution functions can be represented as a combination of four shifted Maxwellian distribution functions. This equilibrium describes a magnetic field configuration with more freedom than the previously known exact solution and has different bulk flow properties.