Nora Wieczorek i Masdeu: Numerical model of cadiac electromechanics

Numerical model of cadiac electromechanics

In this project, we develop a finite element method (FME) formulation that solves the electrophysiological problem on a three dimensional piece of tissue. This problem is modeled by an electromechanical model that includes the activation of a tension that depends on the cell's transmembrane potential and induces the contraction of the tissue. After applying an implicit time discretization, we end with a nonlinear system that depends on the position and potential at each node that we solve using the Newton-Raphson’s method at each time iteration. Using this resolution methodology, we present a full implicit scheme. We also implement a faster, less accurate way of solving the coupled system with a staggered scheme, first computing the change of potential, and then actualizing the position of every node. Then, we simulate the electrophysiological model to observe the effect of the grid affects on the results. Finally, using the staggered algorithm, we simulate the propagation of plane wave and the subsequent tissue contraction.