Related papers: Constructing stellarators with quasisymmetry to hi…
In this paper, we present a new static and time-dependent MagnetoHydroDynamic (MHD) equilibrium code, TokaMaker, for axisymmetric configurations of magnetized plasmas, based on the well-known Grad-Shafranov equation. This code utilizes…
A new approach to high pressure magnetically-confined plasmas is necessary to design efficient fusion devices. This paper presents an equilibrium combining two solutions of the Grad-Shafranov equation, which describes the…
Numerical computation of the ideal Magnetohydrodynamic (MHD) equilibrium magnetic field is at the base of stellarator optimisation and provides the starting point for solving more sophisticated Partial Differential Equations (PDEs) like…
In this work, we utilize new coil objectives for stellarator optimization with autodifferentiation, including pointwise and net coil-coil forces and torques. We use these methods to perform the first large-scale optimization of planar…
In this research, to solve the large indefinite least squares problem, we firstly transform its normal equation into a sparse block three-by-three linear systems, then use GMRES method with an accelerated preconditioner to solve it. The…
To address the magnetization dynamics in ferromagnetic materials described by the Landau-Lifshitz-Gilbert equation under large damping parameters, a third-order accurate numerical scheme is developed by building upon a second-order method…
We demonstrate a fast adjoint-based method to optimize 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…
The problem of finding an optimal curve for the target magnetic axis of a stellarator is addressed. Euler-Lagrange equations are derived for finite length three-dimensional curves that extremise their bending energy while yielding fixed…
When analyzing stellarator configurations, it is common to perform an asymptotic expansion about the magnetic axis. This so-called near-axis expansion is convenient for the same reason asymptotic expansions often are, namely, it reduces the…
In this paper, we apply the method of approximate transformation groups proposed by Baikov, Gaziziv and Ibragimov, to compute the first-order approximate symmetry for the Gardner equations with the small parameters. We compute the optimal…
We present a new numerical model for solving the Chew-Goldberger-Low system of equations describing a bi-Maxwellian plasma in a magnetic field. Heliospheric and geospace environments are often observed to be in an anisotropic state with…
We present a high-order spatial discretization of a continuum gyrokinetic Vlasov model in axisymmetric tokamak edge plasma geometries. Such models describe the phase space advection of plasma species distribution functions in the absence of…
The design of divertor targets and baffles for optimal heat and particle exhaust from magnetically confined fusion plasmas requires a combination of fast, low-fidelity models (such as EMC3-lite [1]) for scoping studies and high-fidelity…
The quasilinear theory describes the resonant interaction between particles and waves with two coupled equations: one for the evolution of the particle probability density function(\textit{pdf}), the other for the wave spectral energy…
The development of stellarators that use permanent magnet arrays to shape their confining magnetic fields has been a topic of recent interest, but the requirements for how such magnets must be shaped, manufactured, and assembled remain to…
A technical note is given regarding our previous laboratory plasma-astrophysical studies [C.-S. Jao et al., High Energy Density Physics 32, 31-43 (2019) and Y. Chen et al., Nucl. Instrum. Methods Phys. Res., Sect. A 903, 119 (2018)]. In…
The close vicinity of neutron stars remains poorly constrained by observations. Although plenty of data are available for the peculiar class of pulsars we are still unable to deduce the underlying plasma distribution in their magnetosphere.…
By making use of the spatial shape of paired photons, parametric down-conversion allows the generation of two-photon entanglement in a multidimensional Hilbert space. How much entanglement can be generated in this way? In principle, the…
We present numerical solutions for stationary and axisymmetric equilibriums of compact stars associated with extremely strong magnetic fields. The interior of the compact stars is assumed to satisfy ideal magnetohydrodynamic (MHD)…
Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for…