A 2D/3D Discrete Duality Finite Volume Scheme. Application to ECG simulation

Yves Coudiere (1) , Charles Pierre (2) , Olivier Rousseau (3) , Rodolphe Turpault (1)

1 LMJL - Laboratoire de Mathématiques Jean Leray
2 LMAP - Laboratoire de Mathématiques et de leurs Applications [Pau]
3 Department of Mathematics and Statistics [Ottawa]

This paper presents a 2D/3D discrete duality finite volume method for solving heterogeneous and anisotropic elliptic equations on very general unstructured meshes. The scheme is based on the definition of discrete divergence and gradient operators that fulfill a duality property mimicking the Green formula. As a consequence, the discrete problem is proved to be well-posed, symmetric and positive-definite. Standard numerical tests are performed in 2D and 3D and the results are discussed and compared with P1 finite elements ones. At last, the method is used for the resolution of a problem arising in biomathematics: the electrocardiogram simulation on a 2D mesh obtained from segmented medical images.

Keywords: {Finite volumes; Anisotropic heterogeneous diffusion; electrocardiology; 3D discrete duality method

Domains: Numerical Analysis [math.NA]

Yves Coudiere, Charles Pierre, Olivier Rousseau, Rodolphe Turpault. A 2D/3D Discrete Duality Finite Volume Scheme. Application to ECG simulation. International Journal on Finite Volumes, 2009, 6 (1), pp.1-24. ⟨hal-00328251⟩