A fully implicit finite element method for bidomain models of cardiac electromechanics

2013-01-01
We propose a novel, monolithic, and unconditionally stable finite element algorithm for the bidomain-based approach to cardiac electromechanics. We introduce the transmembrane potential, the extracellular potential, and the displacement field as independent variables, and extend the common two-field bidomain formulation of electrophysiology to a three-field formulation of electromechanics. The intrinsic coupling arises from both excitation-induced contraction of cardiac cells and the deformation-induced generation of intra-cellular currents. The coupled reaction-diffusion equations of the electrical problem and the momentum balance of the mechanical problem are recast into their weak forms through a conventional isoparametric Galerkin approach. As a novel aspect, we propose a monolithic approach to solve the governing equations of excitation-contraction coupling in a fully coupled, implicit sense. We demonstrate the consistent linearization of the resulting set of non-linear residual equations. To assess the algorithmic performance, we illustrate characteristic features by means of representative three-dimensional initial-boundary value problems. The proposed algorithm may open new avenues to patient specific therapy design by circumventing stability and convergence issues inherent to conventional staggered solution schemes.
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING

Suggestions

A projection based variational multiscale method for a fluid–fluid interaction problem
Ağgül, Mustafa ; Eroğlu, Fatma Güler ; Kaya Merdan, Songül; Labovsky, Alexer E. (Elsevier BV, 2020-06-15)
The proposed method aims to approximate a solution of a fluid–fluid interaction problem in case of low viscosities. The nonlinear interface condition on the joint boundary allows for this problem to be viewed as a simplified version of the atmosphere–ocean coupling. Thus, the proposed method should be viewed as potentially applicable to air–sea coupled flows in turbulent regime. The method consists of two key ingredients. The geometric averaging approach is used for efficient and stable decoupling of the pr...
A modular regularized variational multiscale proper orthogonal decomposition for incompressible flows
Eroglu, Fatma G.; Kaya Merdan, Songül; Rebholz, Leo G. (Elsevier BV, 2017-10-01)
In this paper, we propose, analyze and test a post-processing implementation of a projection-based variational multiscale (VMS) method with proper orthogonal decomposition (POD) for the incompressible Navier-Stokes equations. The projection-based VMS stabilization is added as a separate post-processing step to the standard POD approximation, and since the stabilization step is completely decoupled, the method can easily be incorporated into existing codes, and stabilization parameters can be tuned independe...
The finite element method for MHD flow at high Hartmann numbers
Nesliturk, AI; Tezer, Münevver (Elsevier BV, 2005-01-01)
A stabilized finite element method using the residual-free bubble functions (RFB) is proposed for solving the governing equations of steady magnetohydrodynamic duct flow. A distinguished feature of the RFB method is the resolving capability of high gradients near the layer rep-ions without refining mesh. We show that the RFB method is stable by proving that the numerical method is coercive even not only at low values but also at moderate and high values of the Hartmann number. Numerical results confirming t...
A fully implicit finite element method for bidomain models of cardiac electrophysiology
Dal, Hüsnü; Göktepe, Serdar (Informa UK Limited, 2012-01-01)
This work introduces a novel, unconditionally stable and fully coupled finite element method for the bidomain system of equations of cardiac electrophysiology. The transmembrane potential phi(i) - phi(e) and the extracellular potential phi(e) are treated as independent variables. To this end, the respective reaction-diffusion equations are recast into weak forms via a conventional isoparametric Galerkin approach. The resultant nonlinear set of residual equations is consistently linearised. The method result...
A modulus gradient model for inhomogeneous materials with isotropic linear elastic constituents
Gülaşık, Hasan; Göktepe, Serdar; Gürses, Ercan (Elsevier BV, 2019-11-01)
The one-dimensional modulus gradient (E-grad) model proposed in Gülaşık et al. (2018) is extended to more general three-dimensional inhomogeneous materials with isotropic linear elastic constituents. In addition to the constitutive equations and balance relations, length scale dependent differential relations for the material parameters of isotropic linear elasticity are provided. The finite element formulation for axisymmetric problems is derived and a model problem of a soft cylindrical rod with a stiff s...
Citation Formats
H. Dal and S. Göktepe, “A fully implicit finite element method for bidomain models of cardiac electromechanics,” COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, pp. 323–336, 2013, Accessed: 00, 2020. [Online]. Available: https://hdl.handle.net/11511/35827.