An improved convection scheme applied to recombining divertor plasma flows (Q1286941)

A high-order finite volume convection scheme, in conjunction with a monotonicity preserving flux-limiter, is applied to a combined tokamak edge plasma/Navier-Stokes neutral transport model. This is a highly nonlinear system of convection-diffusion reaction equations which describe the partially ionized boundary layer plasma of a tokamak fusion reactor. The solutions of interest contain a sharp ionization front. We use the improved convective discretization within the context of matrix-free Newton-Krylov algorithm. It is demonstrated that a matrix-free Newton-Krylov implementation, where the preconditioner is derived using first-order upwind convective differencing, provides savings in both memory requirements and CPU time. \(\copyright\) Academic Press.