Narrow Results

This paper is concerned with the mathematical analysis of a coupled elliptic–parabolic system modeling the interaction between the propagation of electric potential and subsequent deformation of the cardiac tissue. The problem consists in a reaction–diffusion system governing the dynamics of ionic quantities, intra- and extra-cellular potentials, and the linearized elasticity equations are adopted to describe the motion of an incompressible material. The coupling between muscle contraction, biochemical reactions and electric activity is introduced with a so-called active strain decomposition framework, where the material gradient of deformation is split into an active (electrophysiology-dependent) part and an elastic (passive) one. Under the assumption of linearized elastic behavior and a truncation of the updated nonlinear diffusivities, we prove existence of weak solutions to the underlying coupled reaction–diffusion system and uniqueness of regular solutions. The proof of existence is based on a combination of parabolic regularization, the Faedo–Galerkin method, and the monotonicity-compactness method of Lions. A finite element formulation is also introduced, for which we establish existence of discrete solutions and show convergence to a weak solution of the original problem. We close with a numerical example illustrating the convergence of the method and some features of the model.

We follow a formal homogenization approach to investigate the effects of mechanical deformations in electrophysiology models relying on a bidomain description of ionic motion at the microscopic level. To that purpose, we extend these microscopic equations to take into account the mechanical deformations, and proceed by recasting the problem in the framework of classical two-scale homogenization in periodic media, and identifying the equations satisfied by the first coefficients in the formal expansions. The homogenized equations reveal some interesting effects related to the microstructure — and associated with a specific cell problem to be solved to obtain the macroscopic conductivity tensors — in which mechanical deformations play a nontrivial role, i.e. they do not simply lead to a standard bidomain problem posed in the deformed configuration. We then present detailed numerical illustrations of the homogenized model with coupled cardiac electrical–mechanical simulations — all the way to ECG simulations — albeit without taking into account the abundantly-investigated effect of mechanical deformations in ionic models, in order to focus here on other effects. And in fact our numerical results indicate that these other effects are numerically of a comparable order, and therefore cannot be disregarded.