Monte Carlo Simulation Code for Solving Radial Fluid Equations in Toroidal Plasmas

Abstract

We develop a new Monte Carlo simulation method to calculate steady-state solutions of fluid equations for edge plasmas. To confirm the computational principle of the new method, benchmark tests including nonlinear problems in one dimensional (i.e., radial) coordinate space are attempted in the first trial; the code based on the method is called DIPS-1D. We confirm that DIPS-1D is useful for solving a Dirichlet problem and that the solution given by the present method provides sufficient numerical accuracy.