A Finite Volume Method for the Approximation of Diffusion Operators on Distorted Meshes

doi:10.1006/jcph.2000.6466

Journal of Computational Physics

Volume 160, Issue 2, 20 May 2000, Pages 481-499

https://doi.org/10.1006/jcph.2000.6466 Get rights and content

Abstract

A new finite volume method is presented for discretizing general linear or nonlinear elliptic second-order partial-differential equations with mixed boundary conditions. The advantage of this method is that arbitrary distorted meshes can be used without the numerical results being altered. The resulting algorithm has more unknowns than standard methods like finite difference or finite element methods. However, the matrices that need to be inverted are positive definite, so the most powerful linear solvers can be applied. The method has been tested on a few elliptic and parabolic equations, either linear, as in the case of the standard heat diffusion equation, or nonlinear, as in the case of the radiation diffusion equation and the resistive diffusion equation with Hall term.

References (17)

I. Faille
A control volume method to solve an elliptic equation on a two-dimensional irregular meshing
Comput. Meth. Appl. Mech. Eng.
(1992)
J.E. Morel et al.
A cell centered Lagrangian-mesh diffusion differencing scheme
J. Comput. Phys.
(1992)
F. Hermeline
Une méthode de volumes finis pour les équations elliptiques du second ordre
C. R. Acad. Sci. Paris
(1998)
D.S. Kershaw
Differencing of the diffusion equation in Lagrangian hydrodynamic codes
J. Comput. Phys.
(1981)
M. Shashkov et al.
Support-operator finite-difference algorithms for general elliptic problems
J. Comput. Phys.
(1995)
P.G. Ciarlet
The Finite Element Method for Elliptic Problems
(1978)
D.J. Ewing et al.
Rules governing the numbers of nodes and elements in a finite element mesh
Int. J. Num. Meth. Eng.
(1970)
P.L. George et al.
Triangulation de Delaunay et maillage
(1997)

There are more references available in the full text version of this article.

Cited by (186)

A mimetic interpolation-free cell-centered finite volume scheme for the 2D and 3D heterogeneous anisotropic diffusion equations
2024, Journal of Computational and Applied Mathematics
We propose a mimetic interpolation-free cell-centered finite volume scheme for diffusion equations on arbitrary two-dimensional (2D) and three dimensional (3D) unstructured, possibly non star-shaped or nonplanar meshes. This scheme borrows ideas from the well known mimetic finite difference methods. This interpolation-free means that we do not need to introduce auxiliary unknowns in the process of constructing the scheme, and the final scheme has only cell-centered unknowns. After carefully handling continuity of solutions and normal fluxes, we construct piecewise linear approximations on each cell. This makes our scheme applicable to arbitrary discontinuities and arbitrary mesh topologies. Interpolation-free makes the program easy to implement and the 2D and 3D programs have a unified framework. It is worth mentioning that the traditional star-shaped mesh assumption has been abolished, and the new scheme has more relaxed restrictions on meshes. Some representative and even challenging numerical examples verify that the scheme is robust and almost has the optimal convergence order on some benchmark polygonal and polyhedral meshes.
A finite volume method to solve the Poisson equation with jump conditions and surface charges: Application to electroporation
2024, Journal of Computational Physics
Efficient numerical schemes for solving the Poisson equation with jump conditions are of great interest for a variety of problems, including the modeling of electroporation phenomena and filamentary discharges. In this paper, we propose a modification to a finite volume scheme, namely the discrete dual finite volume method, in order to account for jump conditions with surface charges, i.e. with a source term. Our numerical tests demonstrate second-order convergence even with highly distorted meshes. We then apply the proposed method to model electroporation phenomena in biological cells by proposing a model that considers the thickness of the cell membrane as a separate domain, which differs from the literature. We show the advantages of the proposed method in this context through numerical experiments.
Biquadratic element discrete duality finite volume method for solving elliptic equations on quadrilateral mesh
2024, Journal of Computational Physics
Nowadays numerical methods for handling partial differential equations is so vast that it is almost taken for granted. In this paper, we propose a discrete duality finite volume scheme with biquadratic elements and its application in solving two-dimensional elliptic equation on quadrilateral meshes. To test the robustness and efficiency of this method, we investigate several examples from physical and engineering fields with different parameters and coefficients. Numerical experiments, encompassing linear and nonlinear elliptic equations with constant or variable coefficients, reveal that the new method achieves optimal convergence in the continuous $H^{1}$ -norm, with a 2nd-order convergence, and in the $L^{2}$ -norm, with a 3rd-order convergence. Moreover, comparison with the existing biquadratic element finite volume scheme shows that the above method has an advantage of accuracy in the discrete $L^{\infty}$ -norm on several distorted meshes.
Composite finite volume schemes for the diffusion equation on unstructured meshes
2024, Computers and Mathematics with Applications
We present a finite volume scheme for the anisotropic diffusion equation. The scheme is based on a reformulation of the diffusion equation as an advection equation. We prove that it is first order consistent and stable under a parabolic CFL condition. We propose a second order extension with similar properties. We also propose a third order extension. Numerical tests are provided, confirming the expected properties of the scheme.
A nonlinear scheme preserving maximum principle for heterogeneous anisotropic diffusion equation
2024, Journal of Computational and Applied Mathematics
We present a new nonlinear method in order to preserve the maximum principle for the diffusion problem with heterogeneous anisotropic coefficient. It is well-known that for diffusion problems with general discontinuous coefficients, the existing cell-centered finite volume schemes on general meshes cannot unconditionally satisfy the discrete maximum principle (DMP), i.e., there are always certain severe restrictions on the mesh-cell geometry or locations of discontinuity in order to preserve DMP. To deal with the issue we propose a nonlinear method of modifying discrete flux of linear second-order schemes into a new conservative flux. Then the resulting scheme keeps conservation and accuracy, and has the structure preserving maximum principle without imposing those severe restrictions. Specifically, we will start with the diamond scheme whose flux is linear and has second-order accuracy, but has no the structure preserving maximum principle. Then the nonlinear correction method is applied to get a scheme satisfying the discrete maximum principle unconditionally. We give some theoretical analysis including the coercivity property and a prior estimation. Numerical results show that the accuracy of our new scheme is superior to the existing scheme preserving maximum principle, and in some cases the accuracy can even be improved by one order of magnitude.
A dynamic multiscale model of cerebral blood flow and autoregulation in the microvasculature
2023, Applied Mathematical Modelling
Models of the micro-circulatory blood flow in the brain can play a key role in understanding the variety of cerebrovascular diseases that occur in the microvasculature. These conditions are often linked to structural modifications in the vessel network, alterations in the blood flow patterns, as well as impairment in the autoregulatory response, all of which are pathological changes that the model should be able to address if it were to have any clinical value. Furthermore, if the model results were to be validated against clinical MRI data, the model simulations need to be computationally feasible when used on networks on the scale of an MRI voxel. This requires some form of an upscaling approach that bypasses the need for an explicit architectural representation of the whole network while maintaining the relevant anatomical connections. To this end, we developed a hybrid multiscale model of blood flow and autoregulation that traces the dynamic changes in blood flow, volume, and pressure in the cortical microvasculature, where the discrete topology of the penetrating vessels is preserved, and these are then appropriately coupled to the homogenised capillary bed by a spatially distributing support function in the terminal endings. In contrast to the other multiscale models, the model developed here accounts for the dynamic physiological phenomena of the blood flow and the autoregulation processes in the microvessels. We show how the adaptive meshing scheme for the capillary bed developed in this study can be employed to ensure a scale-invariant coupling formulation and numerically accurate simulations, all without compromising the computational feasibility of the model. A statistically accurate cortical network on the scale of an MRI voxel is generated, and the model parameter values are calibrated using a Monte Carlo Filtering analysis to ensure that the model results are physiologically informed. The model is found to be able to capture the steep pressure gradients that have been reported to occur at coupling interfaces. Furthermore, in response to an upstream pressure drop, the network is found to be able to recover cerebral blood flow while exhibiting the characteristic autoregulatory behaviour in terms of changes in vessel calibre and the biphasic flow response. Overall, the model developed here offers a high-quality characterisation of dynamic flow and autoregulation in the microvasculature at improved computational efficiency and lays the ground for whole-brain dynamic simulations.

View all citing articles on Scopus

View full text

Regular ArticleA Finite Volume Method for the Approximation of Diffusion Operators on Distorted Meshes

Abstract

Comput. Meth. Appl. Mech. Eng.

J. Comput. Phys.

C. R. Acad. Sci. Paris

J. Comput. Phys.

J. Comput. Phys.

The Finite Element Method for Elliptic Problems

Rules governing the numbers of nodes and elements in a finite element mesh

Int. J. Num. Meth. Eng.

Triangulation de Delaunay et maillage

Regular Article
A Finite Volume Method for the Approximation of Diffusion Operators on Distorted Meshes