Preserving the accuracy of numerical methods discretizing anisotropic elliptic problems

In this paper we study the loss of precision of numerical methods discretizing anisotropic problems and propose alternative approaches free from this drawback. The deterioration of the accuracy is observed when the coordinates and the mesh are unrelated to the anisotropy direction. While this issue is commonly addressed by increasing the scheme approximation order, we demonstrate that, though the gains are evident, the precision of these numerical methods remain far from optimal and limited to moderate anisotropy strengths. This is analysed and explained by an amplification of the approximation error related to the anisotropy strength. We propose an approach consisting in the introduction of an auxiliary variable aimed at removing the amplification of the discretization error. By this means the precision of the numerical approximation is demonstrated to be independent of the anisotropy strength.

[1]  Cory D. Hauck,et al.  An asymptotic-preserving semi-Lagrangian algorithm for the time-dependent anisotropic heat transport equation , 2014, J. Comput. Phys..

[2]  Barry Koren,et al.  Finite-difference schemes for anisotropic diffusion , 2014, J. Comput. Phys..

[3]  Min Tang,et al.  Uniformly Convergent Scheme for Strongly Anisotropic Diffusion Equations with Closed Field Lines , 2018, SIAM J. Sci. Comput..

[4]  M. Ottaviani,et al.  A flux-coordinate independent field-aligned approach to plasma turbulence simulations , 2013, Comput. Phys. Commun..

[5]  Fabrice Deluzet,et al.  On the numerical resolution of anisotropic equations with high order differential operators arising in plasma physics , 2019, J. Comput. Phys..

[6]  Fabrice Deluzet,et al.  A Two Field Iterated Asymptotic-Preserving Method for Highly Anisotropic Elliptic Equations , 2019, Multiscale Model. Simul..

[7]  Fabrice Deluzet,et al.  Iterative Solvers for Elliptic Problems with Arbitrary Anisotropy Strengths , 2018, Multiscale Model. Simul..

[8]  Jacek Narski,et al.  Asymptotic Preserving scheme for strongly anisotropic parabolic equations for arbitrary anisotropy direction , 2013, Comput. Phys. Commun..

[9]  Ahmed Ratnani,et al.  Anisotropic Diffusion in Toroidal geometries , 2016 .

[10]  Claudia Negulescu,et al.  An asymptotic-preserving method for highly anisotropic elliptic equations based on a Micro-Macro decomposition , 2011, J. Comput. Phys..

[11]  Phillip Colella,et al.  A Fourth-Order Accurate Finite-Volume Method with Structured Adaptive Mesh Refinement for Solving the Advection-Diffusion Equation , 2012, SIAM J. Sci. Comput..

[12]  Barry Koren,et al.  Finite-volume scheme for anisotropic diffusion , 2016, J. Comput. Phys..

[13]  Guillaume Latu,et al.  Comparison of Numerical Solvers for Anisotropic Diffusion Equations Arising in Plasma Physics , 2015, J. Sci. Comput..

[14]  Fabrice Deluzet,et al.  An Asymptotic Preserving Scheme for Strongly Anisotropic Elliptic Problems , 2009, Multiscale Model. Simul..

[15]  Sibylle Günter,et al.  Finite element and higher order difference formulations for modelling heat transport in magnetised plasmas , 2007, J. Comput. Phys..

[16]  Sibylle Günter,et al.  Modelling of heat transport in magnetised plasmas using non-aligned coordinates , 2005 .

[17]  Min Tang,et al.  An asymptotic preserving method for strongly anisotropic diffusion equations based on field line integration , 2017, J. Comput. Phys..

[18]  I. Babuska,et al.  Locking effects in the finite element approximation of elasticity problems , 1992 .