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

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 results in a symmetric set of equations, which reduces the computational time significantly compared to the conventional solution algorithms. The proposed method is inherently modular and can be combined with phenomenological or ionic models across the cell membrane. The efficiency of the method and the comparison of its computational cost with respect to the simplified monodomain models are demonstrated through representative numerical examples.
COMPUTER METHODS IN BIOMECHANICS AND BIOMEDICAL ENGINEERING

Suggestions

A fully implicit finite element method for bidomain models of cardiac electromechanics
Dal, Hüsnü; Göktepe, Serdar (Elsevier BV, 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 gen...
An orthotropic viscoelastic material model for passive myocardium: theory and algorithmic treatment
CANSIZ, F. Baris Can; Dal, Hüsnü; KALISKE, Michael (Informa UK Limited, 2015-08-18)
This contribution presents a novel constitutive model in order to simulate an orthotropic rate-dependent behaviour of the passive myocardium at finite strains. The motivation for the consideration of orthotropic viscous effects in a constitutive level lies in the disagreement between theoretical predictions and experimentally observed results. In view of experimental observations, the material is deemed as nearly incompressible, hyperelastic, orthotropic and viscous. The viscoelastic response is formulated ...
A Kalman filter-based approach to reduce the effects of geometric errors and the measurement noise in the inverse ECG problem
Aydin, Umit; Serinağaoğlu Doğrusöz, Yeşim (Springer Science and Business Media LLC, 2011-09-01)
In this article, we aimed to reduce the effects of geometric errors and measurement noise on the inverse problem of Electrocardiography (ECG) solutions. We used the Kalman filter to solve the inverse problem in terms of epicardial potential distributions. The geometric errors were introduced into the problem via wrong determination of the size and location of the heart in simulations. An error model, which is called the enhanced error model (EEM), was modified to be used in inverse problem of ECG to compens...
An algorithm for constructing various kinds of nanojunctions using zig-zag and armchair nanotubes
Tasci, Emre; Erkoç, Şakir (American Scientific Publishers, 2007-04-01)
A method for generating various forms of junctions involving armchair and zig-zag nanotubes, firstly introduced by Zsoldos et al., is developed to cover all types of armchair and zig-zag nanotubes in a systematical way. This method can also be used to produce nanogears and toothed canals. The method is explained and flowcharts are included to aid in programming into a code.
A unifying grid approach for solving potential flows applicable to structured and unstructured grid configurations
Cete, A. Ruhsen; Yuekselen, M. Adil; Kaynak, Uenver (Elsevier BV, 2008-01-01)
In this study, an efficient numerical method is proposed for unifying the structured and unstructured grid approaches for solving the potential flows. The new method, named as the "alternating cell directions implicit - ACDI", solves for the structured and unstructured grid configurations equally well. The new method in effect applies a line implicit method similar to the Line Gauss Seidel scheme for complex unstructured grids including mixed type quadrilateral and triangle cells. To this end, designated al...
Citation Formats
H. Dal and S. Göktepe, “A fully implicit finite element method for bidomain models of cardiac electrophysiology,” COMPUTER METHODS IN BIOMECHANICS AND BIOMEDICAL ENGINEERING, pp. 645–656, 2012, Accessed: 00, 2020. [Online]. Available: https://hdl.handle.net/11511/38345.