A Multipoint Flux Approximation with a Diamond Stencil and a Non-Linear Defect Correction Strategy for the Numerical Solution of Steady State Diffusion Problems in Heterogeneous and Anisotropic Media Satisfying the Discrete Maximum Principle

In the present paper, we solve the steady state diffusion equation in 3D domains by means of a cell-centered finite volume method that uses a Multipoint Flux Approximation with a Diamond Stencil and a Non-Linear defect correction strategy (MPFA-DNL) to guarantee the Discrete Maximum Principle (DMP)....

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Journal of scientific computing Jg. 93; H. 2; S. 42
Hauptverfasser: Cavalcante, T. M., Filho, R. J. M. Lira, Souza, A. C. R., Carvalho, D. K. E., Lyra, P. R. M.
Format: Journal Article
Sprache:Englisch
Veröffentlicht: New York Springer US 01.11.2022
Springer Nature B.V
Schlagworte:
Algorithms
Anisotropic media
Anisotropy
Approximation
Computational Mathematics and Numerical Analysis
Defects
Diamonds
Diffusion rate
Finite element analysis
Finite volume method
Mathematical analysis
Mathematical and Computational Engineering
Mathematical and Computational Physics
Mathematics
Mathematics and Statistics
Maximum principle
Partial differential equations
Robustness (mathematics)
Steady state
Tensors
Theoretical
Heterogeneous and anisotropic media
3D diffusion problems
MPFA-DNL
Discrete maximum principle (DMP)
Unstructured tetrahedral meshes
ISSN:0885-7474, 1573-7691
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:In the present paper, we solve the steady state diffusion equation in 3D domains by means of a cell-centered finite volume method that uses a Multipoint Flux Approximation with a Diamond Stencil and a Non-Linear defect correction strategy (MPFA-DNL) to guarantee the Discrete Maximum Principle (DMP). Our formulation is based in the fact that the flux of MPFA methods can be split into two different parts: a Two Point Flux Approximation (TPFA) component and the Cross-Diffusion Terms (CDT). In the linear MPFA-D method, this split is particularly simple since it lies at the core of the original method construction. In this context, we introduce a non-linear defect correction, aiming to mitigate, whenever necessary, the contributions from the CDT, avoiding, this way, spurious oscillations and DMP violations. Our new MPFA-DNL scheme is locally conservative and capable of dealing with arbitrary anisotropic diffusion tensors and unstructured meshes, without harming the second order convergence rates of the original MPFA-D. To appraise the accuracy and robustness of our formulation, we solve some benchmark problems found in literature. In this paper, we restrict ourselves to tetrahedral meshes, even though, in principle, there is no restriction to extend the method to other polyhedral control volumes.
Bibliographie:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
ISSN:0885-7474
1573-7691
DOI:10.1007/s10915-022-01978-6