Asymptotic-preserving methods for an anisotropic model of electrical potential in a tokamak

A 2D nonlinear model for the electrical potential in the edge plasma in a tokamak generates a stiff problem due to the low resistivity in the direction parallel to the magnetic field lines. An asymptotic-preserving method based on a micro-macro decomposition is studied in order to have a well-posed problem, even when the parallel resistivity goes to $0$. Numerical tests with a finite difference scheme show a bounded condition number for the linearised discrete problem solved at each time step, which confirms the theoretical analysis on the continuous problem.Comment: 8 page

Numerical resolution of an anisotropic non-linear diffusion problem

International audienceThis paper is devoted to the numerical resolution of an anisotropic non-linear diffusion problem involving a small parameter ε, defined as the anisotropy strength reciprocal. In this work, the anisotropy is carried by a variable vector function b. The equation being supplemented with Neumann boundary conditions, the limit ε → 0 is demonstrated to be a singular perturbation of the original diffusion equation. To address efficiently this problem, an Asymptotic-Preserving scheme is derived. This numerical method does not require the use of coordinates adapted to the anisotropy direction and exhibits an accuracy as well as a computational cost independent of the anisotropy strength

Degenerate anisotropic elliptic problems and magnetized plasma simulations

This paper is devoted to the numerical approximation of a degenerate anisotropic elliptic problem. The numerical method is designed for arbitrary space-dependent anisotropy directions and does not require any specially adapted coordinate system. It is also designed to be equally accurate in the strongly and the mildly anisotropic cases. The method is applied to the Euler-Lorentz system, in the drift-fluid limit. This system provides a model for magnetized plasmas

Asymptotic-Preserving methods and multiscale models for plasma physics

The purpose of the present paper is to provide an overview of As ymptotic- Preserving methods for multiscale plasma simulations by ad dressing three sin- gular perturbation problems. First, the quasi-neutral lim it of fluid and kinetic models is investigated in the framework of non magnetized as well as magne- tized plasmas. Second, the drift limit for fluid description s of thermal plasmas under large magnetic fields is addressed. Finally efficient nu merical resolutions of anisotropic elliptic or diffusion equations arising in ma gnetized plasma simu- lation are reviewed

Asymptotic-Preserving scheme for a bi-fluid Euler-Lorentz model

The present work is devoted to the simulation of a strongly magnetized plasma considered as a mixture of an ion fluid and an electron fluid. For the sake of simplicity, we assume that the model is isothermal and described by Euler equations coupled with a term representing the Lorentz force. Moreover we assume that both Euler systems are coupled through a quasi-neutrality constraint. The numerical method which is described in the present document is based on an Asymptotic-Preserving semi-discretization in time of a variant of this two-fluid Euler-Lorentz model with a small perturbation of the quasi-neutrality constraint. Firstly, we present the two-fluid model and the motivations for introducing a small perturbation into the quasi-neutrality equation, then we describe the time semi-discretization of the perturbed model and a fully-discrete finite volume scheme based on it. Finally, we present some numerical results which have been obtained with this method