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)....

Celý popis

Uloženo v:
Podrobná bibliografie
Vydáno v:Journal of scientific computing Ročník 93; číslo 2; s. 42
Hlavní autoři: Cavalcante, T. M., Filho, R. J. M. Lira, Souza, A. C. R., Carvalho, D. K. E., Lyra, P. R. M.
Médium: Journal Article
Jazyk:angličtina
Vydáno: New York Springer US 01.11.2022
Springer Nature B.V
Témata:
ISSN:0885-7474, 1573-7691
On-line přístup:Získat plný text
Tagy: Přidat tag
Žádné tagy, Buďte první, kdo vytvoří štítek k tomuto záznamu!
Popis
Shrnutí: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.
Bibliografie: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