A parallel second order Cartesian method for elliptic interface problems
WEYNANS, Lisl
Institut de Mathématiques de Bordeaux [IMB]
Modélisation et calculs pour l'électrophysiologie cardiaque [CARMEN]
Institut de Mathématiques de Bordeaux [IMB]
Modélisation et calculs pour l'électrophysiologie cardiaque [CARMEN]
WEYNANS, Lisl
Institut de Mathématiques de Bordeaux [IMB]
Modélisation et calculs pour l'électrophysiologie cardiaque [CARMEN]
< Reduce
Institut de Mathématiques de Bordeaux [IMB]
Modélisation et calculs pour l'électrophysiologie cardiaque [CARMEN]
Language
en
Article de revue
This item was published in
Communications in Computational Physics. 2012-06-12, vol. 12, p. 1562-1587
Global Science Press
Abstract
Nous présentons dans ce rapport une méthode cartésienne parallèle pour résoudre des problèmes elliptiques avec des interfaces complexes. Cette méthode est basée sur un schéma aux différences finies et est d'ordre deux dans ...Read more >
Nous présentons dans ce rapport une méthode cartésienne parallèle pour résoudre des problèmes elliptiques avec des interfaces complexes. Cette méthode est basée sur un schéma aux différences finies et est d'ordre deux dans tout le domaine. Nous utilisons un schéma standart à cinq points pour la discrétisation de l'opérateur elliptique sur tous les points de grille, couplé à une discrétisation des conditions de transmission à travers l'interface. L'originalité de la méthode consiste en l'utilisation d'inconnues additionnelles situées sur l'interface, là où sont imposées les conditions de transmission. Nous décrivons d'abord la méthode elle-même et les détails de sa parallélisation réalisée avec la librairie PETSc. Puis nous présentons des validations numériques en deux dimensions, accompagnées de comparaisons avec d'autres méthodes du même type, et une étude numérique de sa version parallèle.Read less <
English Abstract
We present in this paper a parallel Cartesian method to solve elliptic problems with complex immersed interfaces. This method is based on a finite difference scheme and second order accurate in the whole domain. We use a ...Read more >
We present in this paper a parallel Cartesian method to solve elliptic problems with complex immersed interfaces. This method is based on a finite difference scheme and second order accurate in the whole domain. We use a standard five point scheme for the discretization of the elliptic operator on all grid points, coupled with the discretization of transmission conditions across the interface. The originality of the method lies in the use of additional unknowns located on the interface, where the transmission conditions are imposed. We firstly describe the method itself and the details of its parallelization performed with the PETSc library. Then we present numerical validations in two dimensions, assorted with comparisons to other related methods, and a numerical study of the parallelized method.Read less <
English Keywords
Elliptic interface problem
Second order scheme
Interface unknowns
Cartesian method
Origin
Hal imported