Computation of Equilibria of a Flowing Two-fluid in Two Dimensions

摘要

Two-fluid equilibria encompass more physics than the Grad-Shafranov (GS) equation for static equilibria, or even the flowing MHD model. However the two-fluid system is more complicated, and, worse yet, is a singular perturbation problem. The latter difficulty is overcome using the "nearby-fluids" ordering. A "1.5D" solution method has been used to interpret results from the TCS experiment. These results, summarized here, exhibit trends indicative of the improved stability and transport observed experimentally. An algorithm for solving 2D equilibria has been developed based on a relaxation method. The magnetic flux function, governed by the extended GS equation, is updated by successive-overrelaxation, while the toroidal field and flow components and the density are updated using a Newton-Raphson-like method.