We model the electrical behavior of several biological cells under external stimuli by extending and computationally improving the multiple traces formulation introduced in Henríquez et al. [Numer. Math. 136 (2016) 101-145]. Therein, the electric potential and current for a single cell are retrieved through the coupling of boundary integral operators and non-linear ordinary differential systems of equations. Yet, the low-order discretization scheme presented becomes impractical when accounting for interactions among multiple cells. In this note, we consider multi-cellular systems and show existence and uniqueness of the resulting non-linear evolution problem in finite time. Our main tools are analytic semigroup theory along with mapping properties of boundary integral operators in Sobolev spaces. Thanks to the smoothness of cellular shapes, solutions are highly regular at a given time. Hence, spectral spatial discretization can be employed, thereby largely reducing the number of unknowns. Time-space coupling is achieved via a semi-implicit time-stepping scheme shown to be stable and second order convergent. Numerical results in two dimensions validate our claims and match observed biological behavior for the Hodgkin-Huxley dynamical model.
|Number of pages||44|
|Journal||ESAIM: Mathematical Modelling and Numerical Analysis|
|State||Published - 1 Mar 2018|
- Biological cells
- Exponential convergence
- Multiple traces formulation
- Semi-implicit time stepping