An efficient method for solving elliptic boundary element problems with application to the tokamak vacuum problem Original Research Article

A method for regularizing ill-posed Neumann Poisson-type problems based on applying operator transformations is presented. This method can be implemented in the context of the finite element method to compute the solution to inhomogeneous Neumann boundary conditions; it allows to overcome cases where the Neumann problem otherwise admits an infinite number of solutions. As a test application, we solve the Grad–Shafranov boundary problem in a toroidally symmetric geometry. Solving the regularized Neumann response problem is found to be several orders of magnitudes more efficient than solving the Dirichlet problem, which makes the approach competitive with the boundary element method without the need to derive a Green function. In the context of the boundary element method, the operator transformation technique can also be applied to obtain the response of the Grad–Shafranov equation from the toroidal Laplace image response matrix using a simple matrix transformation.