Kernel and image of the Biot-Savart operator and their applications in stellarator designs

We consider the Biot-Savart operator acting on $W^{-\frac{1}{2},2}$ regular, div-free, surface currents $j$ \begin{gather} \nonumber \operatorname{BS}(j)(x)=\frac{1}{4\pi}\int_{\Sigma}j(y)\times \frac{x-y}{|x-y|^3}d\sigma(y)\text{, }x\in \Omega \end{gather} where $\Sigma$ is a connected surface to which $j$ is tangent and where $\Omega$ is the finite domain bounded by $\Sigma$. We answer two questions regarding this operator. i) We provide an algorithm which converges (theoretically) exponentially fast to an element of the kernel of the Biot-Savart operator, as well as characterise the elements of the kernel of the Biot-Savart operator in terms of certain solutions to exterior boundary value problems. This allows one to explicitly exploit the non-uniqueness of the coil reconstruction process in the context of stellarator designs. ii) We provide a simple, concise characterisation of the image of the Biot-Savart operator. This allows to define a 2-step current reconstruction procedure to obtain surface currents which approximate to arbitrary precision a prescribed target magnetic field within the plasma region of a stellarator device. The first step does not require computing integrals involving singular integral kernels of the form $\frac{x-y}{|x-y|^3}$ but may have a potentially slow convergence rate, while the second step requires the computation of integrals involving singular integral kernels but in turn has (theoretically) an exponential convergence rate. This approximation procedure always leads to approximating surface currents which are as poloidal as possible.