Monte Carlo Solver for Partly Calculating a Solution to the Poisson Equation in Three-Dimensional Curvilinear Coordinates

Abstract

We develop a new simulation code for solving the Poisson equation, based on Monte Carlo methods. When static resonant magnetic perturbations (RMPs) are used in tokamak plasma to mitigate or suppress edge-localized modes, the RMPs generate an electric field in the ergodized edge region. The electrostatic potential should be calculated only in the edge region to reduce the computational cost of solving the Poisson equation in the complicated three-dimensional magnetic structure, which is assumed to be fixed in time. In this study, we propose a basic idea for evaluating an electrostatic potential given by the Poisson equation in only a part of the domain in curvilinear coordinates. This Poisson solver allows for the boundary condition to be set not only inside the selected region in which the potential is evaluated, but also outside the selected region. Several benchmarks for the developed code are also presented.