An inexact Newton algorithm for solving the tokamak edge plasma fluid equations on a multiply-connected domain (Q1346534)

Newton's method is combined with a preconditioned conjugate gradient-like algorithm and finite volume discretization to solve the steady-state two- dimensional tokamak edge plasma fluid equations. A numerical evaluation of the Jacobian is employed. Mesh sequencing, pseudo-transient continuation, and adaptive damping are used to increase the radius of convergence. The computations are performed on a multiply-connected curvilinear geometry in a fully coupled manner. Simulations of a high recycling divertor and a gaseous divertor on the DIII-D tokamak geometry are used to demonstrate algorithm performance.