par Vanden Abeele, David ;Degrez, Gérard
Référence AIAA journal, 38, 2, page (234-242)
Publication Publié, 2000-02
Article révisé par les pairs
Résumé : A new approach to the numerical modeling of high-pressure inductive plasmas is presented. The governing magnetohydrodynamic equations are discretized in a second-order accurate finite volume manner. A pressure-stabilized flowfield solver is introduced as an alternative to the staggered-mesh solvers used in traditional algorithms. It is argued that the widely used integral boundary formulation for the electric field is computationally expensive and cannot be incorporated into an efficient iterative solution procedure. A far-field formulation of the electric field is adopted instead, such that state-of-the-art iterative methods can be applied to speed up the calculation. The discretized equations are solved through a damped Picard method and an approximate and a full Newton method. Efficient linear algebra methods are used to solve the linear systems arising from the iterative methods. An appropriate linearization of the (strongly positive) joule heating source term is found to be important for the convergence at the linear level. The new model is tested on a 3-species argon and an 11-species air inductive plasma computation. The proposed damped Picard iterative method is found to be very robust during the initial iterations, but does not converge well thereafter. The approximate Newton method converges substantially better. Finally, the full Newton method yields rapid quadratic convergence rates at the expense of an increase in storage. Taking into account both speed of convergence and memory use, the approximate Newton method is considered to be optimal.