No Arabic abstract
It is demonstrated that finite-pressure, approximately quasi-axisymmetric stellarator equilibria can be directly constructed (without numerical optimization) via perturbations of given axisymmetric equilibria. The size of such perturbations is measured in two ways, via the fractional external rotation and, alternatively, via the relative magnetic field strength, i.e. the average size of the perturbed magnetic field, divided by the unperturbed field strength. It is found that significant fractional external rotational transform can be generated by quasi-axisymmetric perturbations, with a similar value of the relative field strength, despite the fact that the former scales more weakly with the perturbation size. High mode number perturbations are identified as a candidate for generating such transform with local current distributions. Implications for the development of a general non-perturbative solver for optimal stellarator equilibria is discussed.
An asymptotic expansion is performed to obtain quasi-axisymmetric magnetic configurations that are weakly non-axisymmetric. A large space of solutions is identified, which satisfy the condition of quasi-axisymmetry on a single magnetic flux surface, while (non-axisymmetric) globally quasi-axisymmetric solutions are shown to not exist, agreeing with the conclusions of previous theoretical work. The solutions found are shown to be geometrically constrained at low aspect ratio or high toroidal period number. Solutions satisfying the more general condition of omnigeneity (generalized quasi-axisymmetry) are also shown to exist, and it is found that quasi-axisymmetric deformations can be superposed with an omnigenous solution, while preserving the property of omnigeneity, effectively extending the space of good configurations. A numerical solution of the first order quasi-axisymmetry problem is demonstrated and compared with solutions found with a widely used MHD equilibrium solver, independently verifying that quasi-axisymmetry is satisfied at the appropriate order. It is thereby demonstrated that approximately quasi-axisymmetric solutions can be directly constructed, i.e. without using numerical search algorithms.
A novel, compact, quasi-axisymmetric configuration is presented which exhibits low fast-particle losses and is stable to ideal MHD instabilities. The design has fast-particle loss rates below 8% for flux surfaces within the half-radius, and is shown to have an MHD-stability limit of a normalised pressure of $langlebetarangle=3%$ where $langlebetarangle$ is volume averaged. The flux surfaces at various plasma betas and currents as calculated using the SPEC equilibrium code are presented. Neoclassical transport coefficients are shown to be similar to an equivalent tokamak, with a distinct banana regime at half-radius. An initial coil design study is presented to assess the feasibility of this configuration as a fusion-relevant experiment.
Abel Inversion is currently used in laser-plasma studies in order to estimate the electronic density $n_e$ from the phase-shift map $delta phi$ obtained via interferometry. The main limitation of the Abel method is due to the assumption of axial symmetry of the electronic density, which is often hardly fulfilled. In this paper we present an improvement to the Abel inversion technique in which the axial symmetry condition is relaxed by means of a truncated Legendre Polinomial expansion in the azimutal angle. With the help of simulated interferograms, we will show that the generalized Abel inversion generates accurate densities maps when applied to non axisymmetric density sources.
A non-relativistic multi-fluid plasma axisymmetric equilibrium model was developed recently to account for the presence of an energetic electron fluid in addition to thermal electron and ion fluids. The equilibrium formulation of a multi-fluid plasma with relativistic energetic electrons is developed and reported in this paper. Relativistic effects in a fluid model approximation can appear in two ways: due to a large macroscopic fluid velocity comparable to the speed of light and large particles microscopic random motion which becomes significant if the temperature becomes comparable to or larger than the electron rest mass-energy. It is found that the axial component of relativistic generalized angular momentum can be used to describe relativistic axisymmetric equilibrium. The formulation is applied to a four-fluid plasma composed of a relativistic energetic electron fluid, a thermal electron fluid, and fluids of two thermal ion species (e.g. proton and boron ions). The four-fluid density expression which is consistent with the electrostatic potential is obtained and applied in the computation. An example equilibrium approximating a four-fluid plasma recently observed in a solenoid-free ECRH sustained spherical torus plasma is calculated and presented. A second equilibrium that extends the energetic electron temperature of the first example to 679keV is calculated revealing significant relativistic effects.
Ray tracing codes are useful to study the electromagnetic wave propagation and absorption in the geometrical optics approximation. In magnetized fusion plasma community, most ray tracing codes assume the plasma density and temperature be functions of the magnetic flux and study waves only inside the last closed flux surface, which are sufficient for the present day tokamak. However, they are difficult to be used for configurations with open magnetic field line plasmas, such as mirror machine and field-reversed-configuration (FRC). We develop a ray tracing code in cylindrical coordinates $(r,phi,z)$ to support arbitrary axisymmetric configurations with both closed and open field lines plasmas. For wave propagation, the cold plasma dispersion relation is usually sufficient, and we require the magnetic field ${bf B}(r,z)$ and species densities $n_{s0}(r,z)$ profiles as input. For wave absorption, we require a further temperature $T_{s0}(r,z)$ profile to solve a hot kinetic plasma dispersion relation. In difference to other ray tracing codes which calculate the imaginary part of wave vector ${bf k}_{perp,i}$ for wave absorption, we calculate the imaginary part of wave frequency $omega_i$, which is shown to be equivalent with the former technique under weak damping approximation. The code can use either numerical or analytical equilibrium. Examples and benchmarks with electron cyclotron wave, lower hybrid wave and ion cyclotron wave for tokamak, spherical tokamak (ST), FRC and mirror machine are shown.