Related papers: A fast, high-order solver for the Grad-Shafranov e…
A fast algorithm (linear in the degrees of freedom) for the solution of linear variable-coefficient rational-order fractional integral and differential equations is described. The approach is related to the ultraspherical method for…
We present a general construction for exact analytic Taylor states in axisymmetric toroidal geometries. In this construction, the Taylor equilibria are fully determined by specifying the aspect ratio, elongation, and triangularity of the…
A high order, deterministic direct numerical method is proposed for the nonrelativistic $2D_{\bf x} \times 3D_{\bf v}$ Vlasov-Maxwell system, coupled with Fokker-Planck-Landau type operators. Such a system is devoted to the modelling of…
Exact solutions to the static equilibrium magnetohydrodynamic equations are presented and discussed for both axially and helically reduced systems. For both symmetries, physical restrictions on the solutions are discussed and it is seen…
The theory of plasma physics offers a number of nontrivial examples of partial differential equations, which can be successfully treated with symmetry methods. We propose the Grad-Shafranov equation which may illustrate the reciprocal…
Self-gravitational force calculation for infinitesimally thin disks is important for studies on the evolution of galactic and protoplanetary disks. Although high-order methods have been developed for hydrodynamic and magneto-hydrodynamic…
The Grad-Shafranov equation is solved using spectral elements for tokamak equilibrium with toroidal rotation. The Grad-Shafranov solver builds upon and extends the NIMEQ code [Howell and Sovinec, Comput. Phys. Commun. 185 (2014) 1415]…
The tangential layers are characterized by a bulk plasma velocity and a magnetic field that are perpendicular to the gradient direction. They have been extensively described in the frame of the Magneto-Hydro-Dynamic (MHD) theory. But the…
We consider the linear stability problem for a 3D cylindrically symmetric equilibrium of the relativistic Vlasov-Maxwell system that describes a collisionless plasma. For an equilibrium whose distribution function decreases monotonically…
Starting from the equations of collisionless linear gyrokinetics for magnetised plasmas with an imposed inhomogeneous magnetic field, we present the first known analytical, closed-form solution for the resulting velocity-space integrals in…
In this paper, we present an interior penalty discontinuous Galerkin finite element scheme for solving diffusion problems with strong anisotropy arising in magnetized plasmas for fusion applications. We demonstrate the accuracy produced by…
We present an efficient numerical method, inspired by transformation optics, for solving the Poisson equation in complex and arbitrarily shaped geometries. The approach operates by mapping the physical domain to a uniform computational…
This work deals with the numerical approximation of plasmas which are confined by the effect of a fast oscillating magnetic field (see \cite{Bostan2012}) in the Vlasov model. The presence of this magnetic field induces oscillations (in…
Several problems in magnetically confined fusion, such as the computation of exterior vacuum fields or the decomposition of the total magnetic field into separate contributions from the plasma and the external sources, are best formulated…
We construct analytic solutions to the generalized Grad-Shafranov equation, which incorporates both toroidal and poloidal flows. This is achieved by adopting a general linearizing ansatz for the free-function terms of the equation and…
Standard reduced models often fail to adequately describe the complex geometric response of tokamak plasmas to strong toroidal rotation. In this work, we present VEQ-R, a computationally efficient spectral solver designed to calculate…
We consider the periodic problem for two-fluid non-isentropic Euler-Maxwell systems in plasmas. By means of suitable choices of symmetrizers and an induction argument on the order of the time-space derivatives of solutions in energy…
We study strategies for increasing the precision in the blurring models by maintaining a complexity in the related numerical linear algebra procedures (matrix-vector product, linear system solution, computation of eigenvalues etc.) of the…
Fusion plasma and space plasma are typical non-equilibrium and nonlinear systems, with the interactions between different species well described by the Vlasov-Fokker-Planck (VFP) equations. The transport of mass, momentum, energy, and…
In this article we develop a high order accurate method to solve the incompressible boundary layer equations in a provably stable manner.~We first derive continuous energy estimates,~and then proceed to the discrete setting.~We formulate…