Computational modeling of passive myocardium

Göktepe, Serdar
Wong, Jonathan
Kuhl, Ellen
This work deals with the computational modeling of passive myocardial tissue within the framework ofmixed, non-linear finite element methods. We consider a recently proposed, convex, anisotropic hyperelastic model that accounts for the locally orthotropic micro-structure of cardiac muscle. A coordinate-free representation of anisotropy is incorporated through physically relevant invariants of the Cauchy-Green deformation tensors and structural tensors of the corresponding material symmetry group. This model, which has originally been designed for exactly incompressible deformations, is extended towards entirely three-dimensional inhomogeneous deformations by additively decoupling the strain energy function into volumetric and isochoric parts along with the multiplicative split of the deformation gradient. This decoupled constitutive structure is then embedded in a mixed finite element formulation through a three-field Hu-Washizu functional whose simultaneous variation with respect to the independent pressure, dilatation, and placement fields results in the associated Euler-Lagrange equations, thereby minimizing the potential energy. This weak form is then consistently linearized for uniform-pressure elements within the framework of an implicit finite element method. To demonstrate the performance of the proposed approach, we present a three-dimensional finite element analysis of a generic biventricular heart model, subjected to physiological ventricular pressure. The parameters employed in the numerical analysis are identified by solving an optimization problem based on six simple shear experiments on explanted cardiac tissue. Copyright (C) 2010 John Wiley & Sons, Ltd.


Computational modeling of electrocardiograms: A finite element approach toward cardiac excitation
Kotikanyadanam, Mohan; Göktepe, Serdar; Kuhl, Ellen (Wiley, 2010-05-01)
The objective of this work is the computational simulation of a patient-specific electrocardiogram (EKG) using a novel, robust, efficient, and modular finite element-based simulation tool for cardiac electrophysiology. We apply a two-variable approach in terms of a fast action potential and a slow recovery variable, whereby the latter phenomenologically summarizes the concentration of ionic currents. The underlying algorithm is based on a staggered solution scheme in which the action potential is introduced...
Computational modeling of chemo-electro-mechanical coupling: A novel implicit monolithic finite element approach
Wong, J.; Göktepe, Serdar; Kuhl, E. (Wiley, 2013-10-01)
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 chemic...
The differential quadrature solution of nonlinear reaction-diffusion and wave equations using several time-integration schemes
Meral, Gulnihal; Tezer, Münevver (Wiley, 2011-04-01)
Three different time-integration schemes, namely the finite difference method (FDM) with a relaxation parameter, the least-squares method (LSM) and the finite element method (FEM), are applied to the differential quadrature (DQM) solution of one-dimensional nonlinear reaction-diffusion and wave equations. In the solution procedure, the space derivatives are discretized using DQM, which may also be used without the need of boundary conditions. The aim of the paper is to find computationally more efficient ti...
Modeling of dislocation-grain boundary interactions in a strain gradient crystal plasticity framework
ÖZDEMİR, İZZET; Yalçınkaya, Tuncay (Springer Science and Business Media LLC, 2014-08-01)
This paper focuses on the continuum scale modeling of dislocation-grain boundary interactions and enriches a particular strain gradient crystal plasticity formulation (convex counter-part of Yal double dagger inkaya et al., J Mech Phys Solids 59:1-17, 2011; Int J Solids Struct 49:2625-2636, 2012) by incorporating explicitly the effect of grain boundaries on the plastic slip evolution. Within the framework of continuum thermodynamics, a consistent extension of the model is presented and a potential type non-...
Asymptotic behavior of linear impulsive integro-differential equations
Akhmet, Marat; YILMAZ, Oğuz (Elsevier BV, 2008-08-01)
Asymptotic equilibria of linear integro-differential equations and asymptotic relations between solutions of linear homogeneous impulsive differential equations and those of linear integro-differential equations are established. A new Gronwall-Bellman type lemma for integro-differential inequalities is proved. An example is given to demonstrate the validity of one of the results.
Citation Formats
S. Göktepe, J. Wong, and E. Kuhl, “Computational modeling of passive myocardium,” INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, pp. 1–12, 2011, Accessed: 00, 2020. [Online]. Available: