Show/Hide Menu
Hide/Show Apps
anonymousUser
Logout
Türkçe
Türkçe
Search
Search
Login
Login
OpenMETU
OpenMETU
About
About
Open Science Policy
Open Science Policy
Frequently Asked Questions
Frequently Asked Questions
Communities & Collections
Communities & Collections
Computational modeling of electrochemical coupling: A novel finite element approach towards ionic models for cardiac electrophysiology
Date
2011-01-01
Author
Wong, Jonathan
Göktepe, Serdar
Kuhl, Ellen
Metadata
Show full item record
This work is licensed under a
Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License
.
Item Usage Stats
8
views
0
downloads
We propose a novel, efficient finite element solution technique to simulate the electrochemical response of excitable cardiac tissue. We apply a global-local split in which the membrane potential of the electrical problem is introduced globally as a nodal degree of freedom, while the state variables of the chemical problem are treated locally as internal variables on the integration point level. This particular discretization is efficient and highly modular since different cardiac cell models can be incorporated in a straightforward way through only minor local modifications on the constitutive level. Here, we derive the underlying algorithmic framework for a recently proposed ionic model for human ventricular cardiomyocytes, and demonstrate its integration into an existing nonlinear finite element infrastructure. To ensure unconditional algorithmic stability, we apply an implicit backward Euler scheme to discretize the evolution equations for both the electrical potential and the chemical state variables in time. To increase robustness and guarantee optimal quadratic convergence, we suggest an incremental iterative Newton-Raphson scheme and illustrate the consistent linearization of the weak form of the excitation problem. This particular solution strategy allows us to apply an adaptive time stepping scheme, which automatically generates small time steps during the rapid upstroke, and large time steps during the plateau, the repolarization, and the resting phases. We demonstrate that solving an entire cardiac cycle for a real patient-specific geometry characterized through a transmembrane potential, four ion concentrations, thirteen gating variables, and fifteen ionic currents requires computation times of less than ten minutes on a standard desktop computer.
Subject Keywords
Nonlinear reaction-diffusion systems
,
Finite element method
,
Adaptive time stepping
,
Electrophysiology
,
Electrochemistry
,
Computational biophysics
URI
https://hdl.handle.net/11511/48457
Journal
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING
DOI
https://doi.org/10.1016/j.cma.2011.07.003
Collections
Department of Civil Engineering, Article