A fast, accurate, and easy to implement Kapur-Rokhlin quadrature scheme for singular integrals in axisymmetric geometries
arXiv:2211.00493 · doi:10.1017/S002237782300020X
Abstract
Many applications in magnetic confinement fusion require the efficient calculation of surface integrals with singular integrands. The singularity subtraction approaches typically used to handle such singularities are complicated to implement and low order accurate. In contrast, we demonstrate that the Kapur-Rokhlin quadrature scheme is well-suited for the logarithmically singular integrals encountered for a toroidally axisymmetric confinement system, is easy to implement, and is high order accurate. As an illustration, we show how to apply this quadrature scheme for the efficient and accurate calculation of the normal component of the magnetic field due to the plasma current on the plasma boundary, via the virtual casing principle.
arXiv comment21 pages, 4 figures
References in corpus (7)
- Quadrature by Expansion: A New Method for the Evaluation of Layer Potentials
- An improved current potential method for fast computation of stellarator coil shapes
- Efficient magnetic fields for supporting toroidal plasmas
- ECOM: a fast and accurate solver for toroidal axisymmetric MHD equilibria
- Taylor States in Stellarators: A Fast High-order Boundary Integral Solver
- Tokamak elongation: how much is too much? II Numerical results
- Efficient high-order singular quadrature schemes in magnetic fusion