A fast, high-order solver for the Grad-Shafranov equation
arXiv:1210.2113 · doi:10.1016/j.jcp.2013.02.045
Abstract
We present a new fast solver to calculate fixed-boundary plasma equilibria in toroidally axisymmetric geometries. By combining conformal mapping with Fourier and integral equation methods on the unit disk, we show that high-order accuracy can be achieved for the solution of the equilibrium equation and its first and second derivatives. Smooth arbitrary plasma cross-sections as well as arbitrary pressure and poloidal current profiles are used as initial data for the solver. Equilibria with large Shafranov shifts can be computed without difficulty. Spectral convergence is demonstrated by comparing the numerical solution with a known exact analytic solution. A fusion-relevant example of an equilibrium with a pressure pedestal is also presented.
25 pages
Cited by in corpus (12)
- An energy-conserving and asymptotic-preserving charged-particle orbit implicit time integrator for arbitrary electromagnetic fields
- An adaptive fast multipole accelerated Poisson solver for complex geometries
- ECOM: a fast and accurate solver for toroidal axisymmetric MHD equilibria
- Neural network tokamak equilibria with incompressible flows
- A Hybridizable Discontinuous Galerkin solver for the Grad-Shafranov equation
- A mimetic spectral element solver for the Grad-Shafranov equation
- Adaptive Hybridizable Discontinuous Galerkin discretization of the Grad-Shafranov equation by extension from polygonal subdomains
- FEQIS: A free-boundary equilibrium solver for integrated modeling of tokamak plasmas
- ATEQ: Adaptive Toroidal Equilibrium code
- An adaptive discontinuous Petrov-Galerkin method for the Grad-Shafranov equation
- Perturbing an axisymmetric magnetic equilibrium to obtain a quasi-axisymmetric stellarator
- Accurate Derivative Evaluation for any Grad-Shafranov Solver