Efficient high-order singular quadrature schemes in magnetic fusion

by   Dhairya Malhotra, et al.

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 in terms of integral equation expressions. Based on Biot-Savart-like formulae, these integrals contain singular integrands. The regularization method commonly used to address the computation of various singular surface integrals along general toroidal surfaces is low-order accurate, and therefore requires a dense computational mesh in order to obtain sufficient accuracy. In this work, we present a fast, high-order quadrature scheme for the efficient computation of these integrals. Several numerical examples are provided demonstrating the computational efficiency and the high-order accurate convergence. A corresponding code for use in the community has been publicly released.


page 7

page 10

page 11

page 12


Computing weakly singular and near-singular integrals in high-order boundary elements

We present algorithms for computing weakly singular and near-singular in...

A High-Order-Accurate 3D Surface Integral Equation Solver for Uniaxial Anisotropic Media

This paper introduces a high-order accurate surface integral equation me...

High-accuracy mesh-free quadrature for trimmed parametric surfaces and volumes

This work presents a high-accuracy, mesh-free, generalized Stokes theore...

Improved Algorithm for Reconstructing Singular Connection in Multi-Block CFD Applications

An improved algorithm is proposed for the reconstruction of singular con...

High order low-bit Sigma-Delta quantization for fusion frames

We construct high order low-bit Sigma-Delta (ΣΔ) quantizers for the vect...

General-purpose kernel regularization of boundary integral equations via density interpolation

This paper presents a general high-order kernel regularization technique...