Numerical analysis of an asymptotic-preserving scheme for anisotropic elliptic equations
For computational scientists solving stiff anisotropic problems, this work provides rigorous convergence guarantees for AP-schemes, though the approach is incremental.
The paper proves ε-independent convergence results for asymptotic-preserving schemes applied to highly anisotropic elliptic equations, ensuring robustness across a wide range of stiffness parameters.
The main purpose of the present paper is to study from a numerical analysis point of view some robust methods designed to cope with stiff (highly anisotropic) elliptic problems. The so-called asymptotic-preserving schemes studied in this paper are very efficient in dealing with a wide range of $\varepsilon$-values, where $0 < \varepsilon \ll 1$ is the stiffness parameter, responsible for the high anisotropy of the problem. In particular, these schemes are even able to capture the macroscopic properties of the system, as $\varepsilon$ tends towards zero, while the discretization parameters remain fixed. The objective of this work shall be to prove some $\varepsilon$-independent convergence results for these numerical schemes and put hence some more rigor in the construction of such AP-methods.