Commun. Comput. Phys., 16 (2014), pp. 841-891.


A Nominally Second-Order Cell-Centered Finite Volume Scheme for Simulating Three-Dimensional Anisotropic Diffusion Equations on Unstructured Grids

Pascal Jacq 1, Pierre-Henri Maire 1*, Remi Abgrall 2

1 CEA/CESTA, 15 Avenue des Sablieres CS 60001, 33116 Le Barp cedex, France.
2 Institut fur Mathematik, Universitat Zurich, CH-8057 Zurich, Switzerland.

Received 31 May 2013; Accepted (in revised version) 17 March 2014
Available online 24 July 2014
doi:10.4208/cicp.310513.170314a

Abstract

We present a finite volume based cell-centered method for solving diffusion equations on three-dimensional unstructured grids with general tensor conduction. Our main motivation concerns the numerical simulation of the coupling between fluid flows and heat transfers. The corresponding numerical scheme is characterized by cell-centered unknowns and a local stencil. Namely, the scheme results in a global sparse diffusion matrix, which couples only the cell-centered unknowns. The space discretization relies on the partition of polyhedral cells into sub-cells and on the partition of cell faces into sub-faces. It is characterized by the introduction of sub-face normal fluxes and sub-face temperatures, which are auxiliary unknowns. A sub-cell-based variational formulation of the constitutive Fourier law allows to construct an explicit approximation of the sub-face normal heat fluxes in terms of the cell-centered temperature and the adjacent sub-face temperatures. The elimination of the sub-face temperatures with respect to the cell-centered temperatures is achieved locally at each node by solving a small and sparse linear system. This system is obtained by enforcing the continuity condition of the normal heat flux across each sub-cell interface impinging at the node under consideration. The parallel implementation of the numerical algorithm and its efficiency are described and analyzed. The accuracy and the robustness of the proposed finite volume method are assessed by means of various numerical test cases.

AMS subject classifications: 65M08, 65F10, 68W10, 76R50

Notice: Undefined variable: pac in /var/www/html/readabs.php on line 165
Key words: Finite volume methods, unstructured grids, anisotropic diffusion, parallel computing.

*Corresponding author.
Email: pascal.jacq@inria.fr (P. Jacq), maire@celia.u-bordeaux1.fr (P.-H. Maire), remi.abgrall@math.uzh.ch (R. Abgrall)
 

The Global Science Journal