Skip to Main content Skip to Navigation
New interface
Journal articles

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

Michaël Leguèbe 1, 2 Clair Poignard 2, 1 Lisl Weynans 3, 2 
1 MONC - Modélisation Mathématique pour l'Oncologie
IMB - Institut de Mathématiques de Bordeaux, Institut Bergonié [Bordeaux], Inria Bordeaux - Sud-Ouest
3 MEMPHIS - Modeling Enablers for Multi-PHysics and InteractionS
IMB - Institut de Mathématiques de Bordeaux, Inria Bordeaux - Sud-Ouest
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.
Document type :
Journal articles
Complete list of metadata
Contributor : Lisl Weynans Connect in order to contact the contributor
Submitted on : Sunday, May 31, 2015 - 11:47:47 PM
Last modification on : Thursday, January 20, 2022 - 5:31:38 PM

Links full text




Michaël Leguèbe, Clair Poignard, Lisl Weynans. A second-order Cartesian method for the simulation of electropermeabilization cell models. Journal of Computational Physics, 2015, 292, pp.26. ⟨10.1016/⟩. ⟨hal-01158377⟩



Record views