Computational modeling of chemo-electro-mechanical coupling: A novel implicit monolithic finite element approach

Wong J., Göktepe S., Kuhl E.

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, cilt.29, ss.1104-1133, 2013 (SCI-Expanded) identifier identifier identifier

  • Yayın Türü: Makale / Tam Makale
  • Cilt numarası: 29
  • Basım Tarihi: 2013
  • Doi Numarası: 10.1002/cnm.2565
  • Dergi Adı: INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING
  • Derginin Tarandığı İndeksler: Science Citation Index Expanded (SCI-EXPANDED), Scopus
  • Sayfa Sayıları: ss.1104-1133
  • Anahtar Kelimeler: multifield, multiscale, finite element method, electrochemistry, electromechanics, HUMAN VENTRICULAR TISSUE, CARDIAC ELECTROPHYSIOLOGY, ELECTROMECHANICAL MODEL, ACTIVE CONTRACTION, HEART-FAILURE, BEATING HEART, MUSCLE, SIMULATION, EXCITATION, CELLS
  • Orta Doğu Teknik Üniversitesi Adresli: Evet

Özet

Computational modeling of the human heart allows us to predict how chemical, electrical, and mechanical fields interact throughout a cardiac cycle. Pharmacological treatment of cardiac disease has advanced significantly over the past decades, yet it remains unclear how the local biochemistry of an individual heart cell translates into global cardiac function. Here, we propose a novel, unified strategy to simulate excitable biological systems across three biological scales. To discretize the governing chemical, electrical, and mechanical equations in space, we propose a monolithic finite element scheme. We apply a highly efficient and inherently modular global-local split, in which the deformation and the transmembrane potential are introduced globally as nodal degrees of freedom, whereas the chemical state variables are treated locally as internal variables. To ensure unconditional algorithmic stability, we apply an implicit backward Euler finite difference scheme to discretize the resulting system in time. To increase algorithmic robustness and guarantee optimal quadratic convergence, we suggest an incremental iterative Newton-Raphson scheme. The proposed algorithm allows us to simulate the interaction of chemical, electrical, and mechanical fields during a representative cardiac cycle on a patient-specific geometry, robust and stable, with calculation times on the order of 4days on a standard desktop computer. Copyright (c) 2013 John Wiley & Sons, Ltd.