A second-order Cartesian method for the simulation of electropermeabilization cell models

Michaël Leguèbe 1, 2 Clair Poignard 2, 1 L. Weynans 3, 2
1 MONC - Modélisation Mathématique pour l'Oncologie
IMB - Institut de Mathématiques de Bordeaux, Inria Bordeaux - Sud-Ouest, Institut Bergonié - CRLCC Bordeaux
3 MEMPHIS - Modeling Enablers for Multi-PHysics and InteractionS
Inria Bordeaux - Sud-Ouest, IMB - Institut de Mathématiques de Bordeaux
Abstract : In this paper, we present a new finite differences method to simulate electropermeabilization models, like the model of Neu and Krassowska or the recent model of Kavian et al. These models are based on the evolution of the electric potential in a cell embedded in a conducting medium. The main feature lies in the transmission of the voltage potential across the cell membrane: the jump of the potential is proportional to the normal flux thanks to the well-known Kirchoff law. An adapted scheme is thus necessary to accurately simulate the voltage potential in the whole cell, notably at the membrane separating the cell from the outer medium. We present a second-order finite differences scheme in the spirit of the method introduced by Cisternino and Weynans for elliptic problems with immersed interfaces. This is a Cartesian grid method based on the accurate discretization of the fluxes at the interface, through the use of additional interface unknowns. The main novelty of our present work lies in the fact that the jump of the potential is proportional to the flux, and therefore is not explicitly known. The original use of interface unknowns makes it possible to discretize the transmission conditions with enough accuracy to obtain a second-order spatial convergence. We prove the second-order spatial convergence in the stationary linear one-dimensional case, and the first-order temporal convergence for the dynamical non-linear model in one dimension. We then perform numerical experiments in two dimensions that corroborate these results.
Type de document :
Article dans une revue
Journal of Computational Physics, Elsevier, 2015, 292, pp.26. 〈10.1016/j.jcp.2015.03.028〉
Liste complète des métadonnées

https://hal.inria.fr/hal-01158377
Contributeur : Lisl Weynans <>
Soumis le : dimanche 31 mai 2015 - 23:47:47
Dernière modification le : mardi 1 mai 2018 - 18:50:02

Identifiants

Collections

Citation

Michaël Leguèbe, Clair Poignard, L. Weynans. A second-order Cartesian method for the simulation of electropermeabilization cell models. Journal of Computational Physics, Elsevier, 2015, 292, pp.26. 〈10.1016/j.jcp.2015.03.028〉. 〈hal-01158377〉

Partager

Métriques

Consultations de la notice

246