Models of cardiac mechanics are increasingly used to investigate cardiac physiology. These models are characterized by a high level of complexity, including the particular anisotropic material properties of biological tissue and the actively contracting material. A large number of independent simulation codes have been developed, but a consistent way of verifying the accuracy and replicability of simulations is lacking. To aid in the verification of current and future cardiac mechanics solvers, this study provides three benchmark problems for cardiac mechanics. These benchmark problems test the ability to accurately simulate pressure-type forces that depend on the deformed objects geometry, anisotropic and spatially varying material properties similar to those seen in the left ventricle and active contractile forces. The benchmark was solved by 11 different groups to generate consensus solutions, with typical differences in higher-resolution solutions at approximately 0.5%, and consistent results between linear, quadratic and cubic finite elements as well as different approaches to simulating incompressible materials. Online tools and solutions are made available to allow these tests to be effectively used in verification of future cardiac mechanics software.

No abstract

No abstract

The mechanisms underlying the ST segment shifts associated with subendocardial ischemia remain unclear. The aim of this paper is to shed further light on the subject through numerical simulations of these shifts. A realistic three-dimensional model of the ventricles, including fiber rotation and anisotropy, is embedded in a nonhomogeneous torso model. A simplification of the bidomain model is used to calculate only the ST segment shift, assuming known values of the transmembrane potential during the plateau and rest phases. A similar simulation is performed in two dimensions. The simulation results suggest that subendocardial ischemia can be located by ST segment shift on the epicardial and torso surfaces. It is shown that ST elevation is associated with the transmural ischemic boundary, while ST depression is associated with the lateral ischemic boundaries.

The bidomain model, coupled with accurate models of cell membrane kinetics, is generally believed to provide a reasonable basis for numerical simulations of cardiac electrophysiology. Because of changes occurring in very short time intervals and over small spatial domains, discretized versions of these models must be solved on fine computational grids, and small time-steps must be applied. This leads to huge computational challenges that have been addressed by several authors. One popular way of reducing the CPU demands is to approximate the bidomain model by the monodomain model, and thus reducing a two by two set of partial differential equations to one scalar partial differential equation; both of which are coupled to a set of ordinary differential equations modeling the cell membrane kinetics. A reduction in CPU time of two orders of magnitude has been reported. It is the purpose of the present paper to provide arguments that such a reduction is not present when order-optimal numerical methods are applied. Theoretical considerations and numerical experiments indicate that the reduction factor of the CPU requirements from bidomain to monodomain computations, using order-optimal methods, typically is about 10 for simple cell models and less than two for more complex cell models.

No abstract

An emerging class of models has been developed in recent years to predict cardiac growth and remodeling (G&R). We recently developed a cardiac G&R constitutive model that predicts remodeling in response to elevated hemodynamics loading, and a subsequent reversal of the remodeling process when the loading is reduced. Here, we describe the integration of this G&R model to an existing strongly coupled electromechanical model of the heart. A separation of timescale between growth deformation and elastic deformation was invoked in this integrated electromechanical-growth heart model. To test our model, we applied the G&R scheme to simulate the effects of myocardial infarction in a realistic left ventricular (LV) geometry using the finite element method. We also simulate the effects of a novel therapy that is based on alteration of the infarct mechanical properties. We show that our proposed model is able to predict key features that are consistent with experiments. Specifically, we show that the presence of a non-contractile infarct leads to a dilation of the left ventricle that results in a rightward shift of the pressure volume loop. Our model also predicts that G&R is attenuated by a reduction in LV dilation when the infarct stiffness is increased.

Two mathematical models are part of the foundation of Computational neurophysiology; (a) the Cable equation is used to compute the membrane potential of neurons, and, (b) volume-conductor theory describes the extracellular potential around neurons. In the standard procedure for computing extracellular potentials, the transmembrane currents are computed by means of (a) and the extracellular potentials are computed using an explicit sum over analytical point-current source solutions as prescribed by volume conductor theory. Both models are extremely useful as they allow huge simplifications of the computational efforts involved in computing extracellular potentials. However, there are more accurate, though computationally very expensive, models available where the potentials inside and outside the neurons are computed simultaneously in a self-consistent scheme. In the present work we explore the accuracy of the classical models (a) and (b) by comparing them to these more accurate schemes. The main assumption of (a) is that the ephaptic current can be ignored in the derivation of the Cable equation. We find, however, for our examples with stylized neurons, that the ephaptic current is comparable in magnitude to other currents involved in the computations, suggesting that it may be significant—at least in parts of the simulation. The magnitude of the error introduced in the membrane potential is several millivolts, and this error also translates into errors in the predicted extracellular potentials. While the error becomes negligible if we assume the extracellular conductivity to be very large, this assumption is, unfortunately, not easy to justify a priori for all situations of interest.