High-order monotone finite-volume schemes for 1D elliptic problems

Xavier Blanc (1) , Francois Hermeline (2, 3) , Emmanuel Labourasse (2, 3) , Julie Patela (1, 2)

1 LJLL (UMR_7598) - Laboratoire Jacques-Louis Lions
2 CEA/DAM - CEA/DAM [Arpajon]
3 LIHPC - Laboratoire en Informatique Haute Performance pour le Calcul et la simulation

When solving numerically an elliptic problem, it is important in most applications that the scheme used preserves the positivity of the solution. When using finite volume schemes on deformed mesh, the question has been solved rather recently. Such schemes are usually (at most) second order convergent, and nonlinear. On the other hand, many high-order schemes have been proposed, that do not ensure positivity of the solution. In this paper we propose a very high-order monotone (that is, positivity preserving) numerical method for elliptic problems in 1D. We prove that this method converges to an arbitrary order and is indeed monotone. We also show how to handle discontinuous sources or diffusion coefficients, while keeping the order of convergence. We assess the new scheme, on several test problems, with arbitrary (regular, distorted, random) meshes.

Mots clés: High order; Monotonicity; Finite volume method; Elliptic problems

Domaines: Equations aux dérivées partielles [math.AP]; Analyse numérique [math.NA]

Xavier Blanc, Francois Hermeline, Emmanuel Labourasse, Julie Patela. High-order monotone finite-volume schemes for 1D elliptic problems. 2021. ⟨cea-03421015v2⟩