Boundary Integral Formulation of the Cell-by-Cell Model of Cardiac Electrophysiology

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: de Souza, Giacomo Rosilho, Krause, Rolf, Pezzuto, Simone
Natura: Preprint
Pubblicazione: 2023
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866916302581399552
author de Souza, Giacomo Rosilho
Krause, Rolf
Pezzuto, Simone
author_facet de Souza, Giacomo Rosilho
Krause, Rolf
Pezzuto, Simone
contents We propose a boundary element method for the accurate solution of the cell-by-cell bidomain model of electrophysiology. The cell-by-cell model, also called Extracellular-Membrane-Intracellular (EMI) model, is a system of reaction-diffusion equations describing the evolution of the electric potential within each domain: intra- and extra-cellular space and the cellular membrane. The system is parabolic but degenerate because the time derivative is only in the membrane domain. In this work, we adopt a boundary-integral formulation for removing the degeneracy in the system and recast it to a parabolic equation on the membrane. The formulation is also numerically advantageous since the number of degrees of freedom is sensibly reduced compared to the original model. Specifically, we prove that the boundary-element discretization of the EMI model is equivalent to a system of ordinary differential equations, and we consider a time discretization based on the multirate explicit stabilized Runge-Kutta method. We numerically show that our scheme convergences exponentially in space for the single-cell case. We finally provide several numerical experiments of biological interest.
format Preprint
id arxiv_https___arxiv_org_abs_2302_05281
institution arXiv
publishDate 2023
record_format arxiv
spellingShingle Boundary Integral Formulation of the Cell-by-Cell Model of Cardiac Electrophysiology
de Souza, Giacomo Rosilho
Krause, Rolf
Pezzuto, Simone
Numerical Analysis
Tissues and Organs
We propose a boundary element method for the accurate solution of the cell-by-cell bidomain model of electrophysiology. The cell-by-cell model, also called Extracellular-Membrane-Intracellular (EMI) model, is a system of reaction-diffusion equations describing the evolution of the electric potential within each domain: intra- and extra-cellular space and the cellular membrane. The system is parabolic but degenerate because the time derivative is only in the membrane domain. In this work, we adopt a boundary-integral formulation for removing the degeneracy in the system and recast it to a parabolic equation on the membrane. The formulation is also numerically advantageous since the number of degrees of freedom is sensibly reduced compared to the original model. Specifically, we prove that the boundary-element discretization of the EMI model is equivalent to a system of ordinary differential equations, and we consider a time discretization based on the multirate explicit stabilized Runge-Kutta method. We numerically show that our scheme convergences exponentially in space for the single-cell case. We finally provide several numerical experiments of biological interest.
title Boundary Integral Formulation of the Cell-by-Cell Model of Cardiac Electrophysiology
topic Numerical Analysis
Tissues and Organs
url https://arxiv.org/abs/2302.05281