Electromechanical feedback with reduced cellular connectivity alters electrical activity in an infarct injured left ventricle: a finite element model study

Computational cardiac physiology for new modelers: Origins, foundations, and future

Mathematical models of the cardiovascular system have come a long way since they were first introduced in the early 19th century. Driven by a rapid development of experimental techniques, numerical methods, and computer hardware, detailed models that describe physical scales from the molecular level up to organs and organ systems have been derived and used for physiological research. Mathematical and computational models can be seen as condensed and quantitative formulations of extensive physiological knowledge and are used for formulating and testing hypotheses, interpreting and directing experimental research, and have contributed substantially to our understanding of cardiovascular physiology. However, in spite of the strengths of mathematics to precisely describe complex relationships and the obvious need for the mathematical and computational models to be informed by experimental data, there still exist considerable barriers between experimental and computational physiological research. In this review, we present a historical overview of the development of mathematical and computational models in cardiovascular physiology, including the current state of the art. We further argue why a tighter integration is needed between experimental and computational scientists in physiology, and point out important obstacles and challenges that must be overcome in order to fully realize the synergy of experimental and computational physiological research.

On modeling the multiscale mechanobiology of soft tissues: Challenges and progress

Tissues grow and remodel in response to mechanical cues, extracellular and intracellular signals experienced through various biological events, from the developing embryo to disease and aging. The macroscale response of soft tissues is typically nonlinear, viscoelastic anisotropic, and often emerges from the hierarchical structure of tissues, primarily their biopolymer fiber networks at the microscale. The adaptation to mechanical cues is likewise a multiscale phenomenon. Cell mechanobiology, the ability of cells to transform mechanical inputs into chemical signaling inside the cell, and subsequent regulation of cellular behavior through intra- and inter-cellular signaling networks, is the key coupling at the microscale between the mechanical cues and the mechanical adaptation seen macroscopically. To fully understand mechanics of tissues in growth and remodeling as observed at the tissue level, multiscale models of tissue mechanobiology are essential. In this review, we summarize the state-of-the art modeling tools of soft tissues at both scales, the tissue level response, and the cell scale mechanobiology models. To help the interested reader become more familiar with these modeling frameworks, we also show representative examples. Our aim here is to bring together scientists from different disciplines and enable the future leap in multiscale modeling of tissue mechanobiology.

Computational modelling of mechano-electric feedback and its arrhythmogenic effects in human ventricular models

The current study aims to investigate the role of mechano-electric feedback (MEF) in healthy cardiac cycles and in cardiac arrhythmia using human ventricular models. The numerical formulation of stretch-activated channels (SACs) in terms of the fibre stretch of the myocardium is incorporated into the modified Hill model that describes the myocardium as an electro-visco-active material. Additionally, we propose models of SACs formulated in terms of the rate of stretch along fibre direction and the stretch along sheet direction. We analyze the effect of the three different models for SACs and different material properties on the regular cycles by using electrocardiogram and volume-time curves, and show that the each model of SACs has regionally different influences on the heart model. Moreover, we simulate 'commotio cordis' and 'precordial thump' and demonstrate that MEF plays a major role in the occurrence of fibrillation and defibrillation in the absence of the structural cardiac damage. Furthermore, we study the role of MEF in premature ventricular contraction when the blood pressure is disturbed.

Cardiac electro-mechanical activity in a deforming human cardiac tissue: modeling, existence-uniqueness, finite element computation and application to multiple ischemic disease

In this study, the cardiac electro-mechanical model in a deforming domain is taken with the addition of mechanical feedback and stretch-activated channel current coupled with the ten Tusscher human ventricular cell level model that results in a coupled PDE-ODE system. The existence and uniqueness of such a coupled system in a deforming domain is proved. At first, the existence of a solution is proved in the deformed domain. The local existence of the solution is proved using the regularization and the Faedo-Galerkin technique. Then, the global existence is proved using the energy estimates in appropriate Banach spaces, Gronwall lemma, and the compactness procedure. The existence of the solution in an undeformed domain is proved using the lower semi-continuity of the norms. Uniqueness is proved using Young's inequality, Gronwall lemma, and the Cauchy-Schwartz inequality. For the application purpose, this model is applied to understand the electro-mechanical activity in ischemic cardiac tissue. It also takes care of the development of active tension, conductive, convective, and ionic feedback. The Second Piola-Kirchoff stress tensor arising in Lagrangian mapping between reference and moving frames is taken as a combination of active, passive, and volumetric components. We investigated the effect of varying strength of hyperkalemia and hypoxia, in the ischemic subregions of human cardiac tissue with local multiple ischemic subregions, on the electro-mechanical activity of healthy and ischemic zones. This system is solved numerically using the [Formula: see text] finite element method in space and the implicit-explicit Euler method in time. Discontinuities arising with the modeled multiple ischemic regions are treated to the desired order of accuracy by a simple regularization technique using the interpolating polynomials. We examined the cardiac electro-mechanical activity for several cases in multiple hyperkalemic and hypoxic human cardiac tissue. We concluded that local multiple ischemic subregions severely affect the cardiac electro-mechanical activity more, in terms of action potential (v) and mechanical parameters, intracellular calcium ion concentration [Formula: see text], active tension ([Formula: see text]), stretch ([Formula: see text]) and stretch rate ([Formula: see text]), of a healthy cell in its vicinity, compared to a single Hyperkalemic or Hypoxic subregion. The four moderate hypoxically generated ischemic subregions affect the waveform of the stretch along the fiber and the stretch rate more than a single severe ischemic subregion.

Numerical evaluation of cardiac mechanical markers as estimators of the electrical activation time

Recent advances in the development of noninvasive cardiac imaging technologies have made it possible to measure longitudinal and circumferential strains at a high spatial resolution also at intramural level. Local mechanical activation times derived from these strains can be used as noninvasive estimates of electrical activation, in order to determine, eg, the origin of premature ectopic beats during focal arrhythmias or the pathway of reentrant circuits. The aim of this work is to assess the reliability of mechanical activation time markers derived from longitudinal and circumferential strains, denoted by AT_ell and AT_ecc , respectively, by means of three-dimensional cardiac electromechanical simulations. These markers are compared against the electrical activation time (AT_v ), computed from the action potential waveform, and the reference mechanical activation markers derived from the active tension and fiber strain waveforms, denoted by AT_ta and AT_eff , respectively. Our numerical simulations are based on a strongly coupled electromechanical model, including bidomain representation of the cardiac tissue, mechanoelectric (ie, stretch-activated channels) and geometric feedbacks, transversely isotropic strain energy function for the description of passive mechanics and detailed membrane and excitation-contraction coupling models. The results have shown that, during endocardial and epicardial ectopic stimulations, all the mechanical markers considered are highly correlated with AT_v , exhibiting correlation coefficients larger than 0.8. However, during multiple endocardial stimulations, mimicking the ventricular sinus rhythm, the mechanical markers are less correlated with the electrical activation time, because of the more complex resulting excitation sequence. Moreover, the inspection of the endocardial and epicardial isochrones has shown that the AT_ell and AT_ecc mechanical activation sequences reproduce only some qualitative features of the electrical activation sequence, such as the areas of early and late activation, but in some cases, they might yield wrong excitation sources and significantly different isochrones patterns.

A finite element model of the cardiac ventricles with coupled circulation: Biventricular mesh generation with hexahedral elements, airbags and a functional mockup interface to the circulation

INTRODUCTION

Finite element (FE) mechanics models of the heart are becoming more sophisticated. However, there is lack of consensus about optimal element type and coupling of FE models to the circulation. We describe biventricular (left (LV) and right (RV) ventricles) FE mechanics model creation using hexahedral elements, airbags and a functional mockup interface (FMI) to lumped-parameter models of the circulation.

METHODS

Cardiac MRI (CMR) was performed in two healthy volunteers and a single patient with ischemic heart disease (IHD). CMR images were segmented and surfaced, meshing with hexahedral elements was performed with a "thin butterfly with septum" topology. LV and RV inflow and outflow airbags were coupled to lumped-parameter circulation models with an FMI interface. Pulmonary constriction (PAC) and vena cava occlusion (VCO) were simulated and end-systolic pressure-volume relations (ESPVR) were calculated.

RESULTS

Mesh construction was prompt with representative contouring and mesh adjustment requiring 32 and 26 min Respectively. The numbers of elements ranged from 4104 to 5184 with a representative Jacobian of 1.0026 ± 0.4531. Agreement between CMR-based surfaces and mesh was excellent with root-mean-squared error of 0.589 ± 0.321 mm. The LV ESPVR slope was 3.37 ± 0.09 in volunteers but 2.74 in the IHD patient. The effect of PAC and VCO on LV ESPVR was consistent with ventricular interaction (p = 0.0286).

CONCLUSION

Successful co-simulation using a biventricular FE mechanics model with hexahedral elements, airbags and an FMI interface to lumped-parameter model of the circulation was demonstrated. Future studies will include comparison of element type and study of cardiovascular pathologies and device therapies.

A clinical-in silico study on the effectiveness of multipoint bicathodic and cathodic-anodal pacing in cardiac resynchronization therapy

Up to one-third of patients undergoing cardiac resynchronization therapy (CRT) are nonresponders. Multipoint bicathodic and cathodic-anodal left ventricle (LV) stimulations could overcome this clinical challenge, but their effectiveness remains controversial. Here we evaluate the performance of such stimulations through both in vivo and in silico experiments, the latter based on computer electromechanical modeling. Seven patients, all candidates for CRT, received a quadripolar LV lead. Four stimulations were tested: right ventricular (RVS); conventional single point biventricular (S-BS); multipoint biventricular bicathodic (CC-BS) and multipoint biventricular cathodic-anodal (CA-BS). The following parameters were processed: QRS duration; maximal time derivative of arterial pressure (dPdt_max); systolic arterial pressure (P_sys); and stroke volume (SV). Echocardiographic data of each patient were then obtained to create an LV geometric model. Numerical simulations were based on a strongly coupled Bidomain electromechanical coupling model. Considering the in vivo parameters, when comparing S-BS to RVS, there was no significant decrease in SV (from 45 ± 11 to 44 ± 20 ml) and 6% and 4% increases of dPdt_max and P_sys, respectively. Focusing on in silico parameters, with respect to RVS, S-BS exhibited a significant increase of SV, dPdt_max and P_sys. Neither the in vivo nor in silico results showed any significant hemodynamic and electrical difference among S-BS, CC-BS and CA-BS configurations. These results show that CC-BS and CA-BS yield a comparable CRT performance, but they do not always yield improvement in terms of hemodynamic parameters with respect to S-BS. The computational results confirmed the in vivo observations, thus providing theoretical support to the clinical experiments.

Reperfused vs. nonreperfused myocardial infarction: when to use which model

There is a lack of understanding in the cardiac remodeling field regarding the use of nonreperfused myocardial infarction (MI) and reperfused MI in animal models of MI. This Perspectives summarizes the consensus of the authors regarding how to select the optimum model for your experiments and is a part of ongoing efforts to establish rigor and reproducibility in cardiac physiology research.

Biomechanics of infarcted left ventricle: a review of modelling

Mathematical modelling in biomechanics of infarcted left ventricle (LV) serves as an indispensable tool for remodelling mechanism exploration, LV biomechanical property estimation and therapy assessment after myocardial infarction (MI). However, a review of mathematical modelling after MI has not been seen in the literature so far. In the paper, a systematic review of mathematical models in biomechanics of infarcted LV was established. The models include comprehensive cardiovascular system model, essential LV pressure-volume and stress-stretch models, constitutive laws for passive myocardium and scars, tension models for active myocardium, collagen fibre orientation optimization models, fibroblast and collagen fibre growth/degradation models and integrated growth-electro-mechanical model after MI. The primary idea, unique characteristics and key equations of each model were identified and extracted. Discussions on the models were provided and followed research issues on them were addressed. Considerable improvements in the cardiovascular system model, LV aneurysm model, coupled agent-based models and integrated electro-mechanical-growth LV model are encouraged. Substantial attention should be paid to new constitutive laws with respect to stress-stretch curve and strain energy function for infarcted passive myocardium, collagen fibre orientation optimization in scar, cardiac rupture and tissue damage and viscoelastic effect post-MI in the future.

Multisite pacing and myocardial scars: a computational study

Cardiac resynchronization therapy (CRT) is a frequently effective treatment modality for dyssynchronous heart failure, however, 30% of patients do not respond, usually due to suboptimal activation of the left ventricle (LV). Multisite pacing (MSP) may increase the response rate, but its effect in the presence of myocardial scars is not fully understood. We use a computational model to study the outcome of MSP in an LV with scars in two different locations and of two different sizes. The LV was stimulated from anterior, posterior and lateral locations individually and in pairs, while a septal stimulation site represented right ventricular (RV) pacing. Intraventricular pressures were measured, and outcomes evaluated in terms of maximum LV pressure gradient (dP/dt_max)- change compared to isolated RV pacing. The best result obtained using various LV pacing locations included a combination of sites remote from scars and the septum. The highest dP/dt_max increase was achieved, regardless of scar size, using MSP with one pacing site located on the LV free wall opposite to the scar and one site opposite to the septum. These in silico modelling results suggest that making placement of pacing electrodes dependent on location of scarring, may alter acute haemodynamics and that such modelling may contribute to future CRT optimization.

Electromechanical Model to Predict Cardiac Resynchronization Therapy

Cardiac resynchronization therapy (CRT) can substantially improve dyssynchronous heart failure and reduce mortality. However, one-third of the CRT patients derive no measurable benefit from CRT, due to suboptimal placement of the left ventricular (LV) lead. We introduce a pipeline for improved CRT-therapy by creating an electromechanical model using patient-specific geometric parameters allowing individualization of therapy. The model successfully mimics expected changes when variables for tension, stiffness, and conduction are entered. Changing LV pacing site had a notable effect on maximum pressure gradient (dP/dt_max) in the presence of cardiac scarring, causing non-uniform excitation propagation through the LV. Tailoring CRT to the individual requires simulations with patient-specific biventricular meshes including cardiac geometry and conductivity properties.

High Spatial Resolution Multi-Organ Finite Element Modeling of Ventricular-Arterial Coupling

While it has long been recognized that bi-directional interaction between the heart and the vasculature plays a critical role in the proper functioning of the cardiovascular system, a comprehensive study of this interaction has largely been hampered by a lack of modeling framework capable of simultaneously accommodating high-resolution models of the heart and vasculature. Here, we address this issue and present a computational modeling framework that couples finite element (FE) models of the left ventricle (LV) and aorta to elucidate ventricular-arterial coupling in the systemic circulation. We show in a baseline simulation that the framework predictions of (1) LV pressure-volume loop, (2) aorta pressure-diameter relationship, (3) pressure-waveforms of the aorta, LV, and left atrium (LA) over the cardiac cycle are consistent with the physiological measurements found in healthy human. To develop insights of ventricular-arterial interactions, the framework was then used to simulate how alterations in the geometrical or, material parameter(s) of the aorta affect the LV and vice versa. We show that changing the geometry and microstructure of the aorta model in the framework led to changes in the functional behaviors of both LV and aorta that are consistent with experimental observations. On the other hand, changing contractility and passive stiffness of the LV model in the framework also produced changes in both the LV and aorta functional behaviors that are consistent with physiology principles.

Inverse estimation of cardiac activation times via gradient-based optimization

Computational modeling may provide a quantitative framework for integrating multiscale data to gain insight into mechanisms of heart disease, identify and test pharmacological and electrical therapy and interventions, and support clinical decisions. Patient-specific computational cardiac models can help guide such procedures, and cardiac inverse modeling is a promising alternative to adequately personalize these models. Indeed, full cardiac inverse modeling is currently becoming computationally feasible; however, fundamental work to assess the feasibility of emerging techniques is still needed. In this study, we use a partial differential equation-constrained optimal control approach to numerically investigate the identifiability of an initial activation sequence from synthetic (partial) observations of the extracellular potential using the bidomain approximation and 2D representations of cardiac tissue. Our results demonstrate that activation times and duration of several stimuli can be recovered even with high levels of noise, that it is sufficient to sample the observations at the electrocardiogram-relevant sampling frequency of 1 kHz, and that spatial resolutions that are coarser than the standard in electrophysiological simulations can be used. The optimization of activation times is still effective when synthetic data are generated with a different cell membrane kinetics model than optimized for. The findings thus indicate that the presented approach has potential for finding activation sequences from clinical data modalities, as an extension to existing cardiac imaging approaches.

Effects of mechanical feedback on the stability of cardiac scroll waves: A bidomain electro-mechanical simulation study

In this work, we investigate the influence of cardiac tissue deformation on re-entrant wave dynamics. We have developed a 3D strongly coupled electro-mechanical Bidomain model posed on an ideal monoventricular geometry, including fiber direction anisotropy and stretch-activated currents (SACs). The cardiac mechanical deformation influences the bioelectrical activity with two main mechanical feedback: (a) the geometric feedback (GEF) due to the presence of the deformation gradient in the diffusion coefficients and in a convective term depending on the deformation rate and (b) the mechano-electric feedback (MEF) due to SACs. Here, we investigate the relative contribution of these two factors with respect to scroll wave stability. We extend the previous works [Keldermann et al., Am. J. Physiol. Heart Circ. Physiol. 299, H134-H143 (2010) and Hu et al., PLoS One 8(4), e60287 (2013)] that were based on the Monodomain model and a simple non-selective linear SAC, while here we consider the full Bidomain model and both selective and non-selective components of SACs. Our simulation results show that the stability of cardiac scroll waves is influenced by MEF, which in case of low reversal potential of non-selective SACs might be responsible for the onset of ventricular fibrillation; GEF increases the scroll wave meandering but does not determine the scroll wave stability.

The role of infarct transmural extent in infarct extension: A computational study

Infarct extension, a process involving progressive extension of the infarct zone (IZ) into the normally perfused border zone (BZ), leads to continuous degradation of the myocardial function and adverse remodelling. Despite carrying a high risk of mortality, detailed understanding of the mechanisms leading to BZ hypoxia and infarct extension remains unexplored. In the present study, we developed a 3D truncated ellipsoidal left ventricular model incorporating realistic electromechanical properties and fibre orientation to examine the mechanical interaction among the remote, infarct and BZs in the presence of varying infarct transmural extent (TME). Localized highly abnormal systolic fibre stress was observed at the BZ, owing to the simultaneous presence of moderately increased stiffness and fibre strain at this region, caused by the mechanical tethering effect imposed by the overstretched IZ. Our simulations also demonstrated the greatest tethering effect and stress in BZ regions with fibre direction tangential to the BZ-remote zone boundary. This can be explained by the lower stiffness in the cross-fibre direction, which gave rise to a greater stretching of the IZ in this direction. The average fibre strain of the IZ, as well as the maximum stress in the sub-endocardial layer, increased steeply from 10% to 50% infarct TME, and slower thereafter. Based on our stress-strain loop analysis, we found impairment in the myocardial energy efficiency and elevated energy expenditure with increasing infarct TME, which we believe to place the BZ at further risk of hypoxia. Copyright © 2016 John Wiley & Sons, Ltd.

Joint influence of transmural heterogeneities and wall deformation on cardiac bioelectrical activity: A simulation study

The aim of this work is to investigate, by means of numerical simulations, the influence of myocardial deformation due to muscle contraction and relaxation on the cardiac repolarization process in presence of transmural intrinsic action potential duration (APD) heterogeneities. The three-dimensional electromechanical model considered consists of the following four coupled components: the quasi-static transversely isotropic finite elasticity equations for the deformation of the cardiac tissue; the active tension model for the intracellular calcium dynamics and cross-bridge binding; the anisotropic Bidomain model for the electrical current flow through the deforming cardiac tissue; the membrane model of ventricular myocytes, including stretch-activated channels. The numerical simulations are based on our finite element parallel solver, which employs Multilevel Additive Schwarz preconditioners for the solution of the discretized Bidomain equations and Newton-Krylov methods for the solution of the discretized non-linear finite elasticity equations. Our findings show that: (i) the presence of intrinsic transmural cellular APD heterogeneities is not fully masked by electrotonic current flow or by the presence of the mechanical deformation; (ii) despite the presence of transmural APD heterogeneities, the recovery process follows the activation sequence and there is no significant transmural repolarization gradient; (iii) with or without transmural APD heterogeneities, epicardial electrograms always display the same wave shape and discordance between the polarity of QRS complex and T-wave; (iv) the main effects of the mechanical deformation are an increase of the dispersion of repolarization time and APD, when computed over the total cardiac domain and over the endo- and epicardial surfaces, while there is a slight decrease along the transmural direction.

Multiphysics and multiscale modelling, data-model fusion and integration of organ physiology in the clinic: ventricular cardiac mechanics

With heart and cardiovascular diseases continually challenging healthcare systems worldwide, translating basic research on cardiac (patho)physiology into clinical care is essential. Exacerbating this already extensive challenge is the complexity of the heart, relying on its hierarchical structure and function to maintain cardiovascular flow. Computational modelling has been proposed and actively pursued as a tool for accelerating research and translation. Allowing exploration of the relationships between physics, multiscale mechanisms and function, computational modelling provides a platform for improving our understanding of the heart. Further integration of experimental and clinical data through data assimilation and parameter estimation techniques is bringing computational models closer to use in routine clinical practice. This article reviews developments in computational cardiac modelling and how their integration with medical imaging data is providing new pathways for translational cardiac modelling.

Biomechanics of Cardiac Function

The heart pumps blood to maintain circulation and ensure the delivery of oxygenated blood to all the organs of the body. Mechanics play a critical role in governing and regulating heart function under both normal and pathological conditions. Biological processes and mechanical stress are coupled together in regulating myocyte function and extracellular matrix structure thus controlling heart function. Here, we offer a brief introduction to the biomechanics of left ventricular function and then summarize recent progress in the study of the effects of mechanical stress on ventricular wall remodeling and cardiac function as well as the effects of wall mechanical properties on cardiac function in normal and dysfunctional hearts. Various mechanical models to determine wall stress and cardiac function in normal and diseased hearts with both systolic and diastolic dysfunction are discussed. The results of these studies have enhanced our understanding of the biomechanical mechanism in the development and remodeling of normal and dysfunctional hearts. Biomechanics provide a tool to understand the mechanism of left ventricular remodeling in diastolic and systolic dysfunction and guidance in designing and developing new treatments.

An integrated electromechanical-growth heart model for simulating cardiac therapies

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.

Increased cell membrane capacitance is the dominant mechanism of stretch-dependent conduction slowing in the rabbit heart: a computational study

Volume loading of the cardiac ventricles is known to slow electrical conduction in the rabbit heart, but the mechanisms remain unclear. Previous experimental and modeling studies have investigated some of these mechanisms, including stretch-activated membrane currents, reduced gap junctional conductance, and altered cell membrane capacitance. In order to quantify the relative contributions of these mechanisms, we combined a monomain model of rabbit ventricular electrophysiology with a hyperelastic model of passive ventricular mechanics. First, a simplified geometric model with prescribed homogeneous deformation was used to fit model parameters and characterize individual MEF mechanisms, and showed good qualitative agreement with experimentally measured strain-CV relations. A 3D model of the rabbit left and right ventricles was then compared with experimental measurements from optical electrical mapping studies in the isolated rabbit heart. The model was inflated to an end-diastolic pressure of 30 mmHg, resulting in epicardial strains comparable to those measured in the anterior left ventricular free wall. While the effects of stretch activated channels did alter epicardial conduction velocity, an increase in cellular capacitance was required to explain previously reported experimental results. The new results suggest that for large strains, various mechanisms can combine and produce a biphasic relationship between strain and conduction velocity. However, at the moderate strains generated by high end-diastolic pressure, a stretch-induced increase in myocyte membrane capacitance is the dominant driver of conduction slowing during ventricular volume loading.

Applications of computational modeling in cardiac surgery

Although computational modeling is common in many areas of science and engineering, only recently have advances in experimental techniques and medical imaging allowed this tool to be applied in cardiac surgery. Despite its infancy in cardiac surgery, computational modeling has been useful in calculating the effects of clinical devices and surgical procedures. In this review, we present several examples that demonstrate the capabilities of computational cardiac modeling in cardiac surgery. Specifically, we demonstrate its ability to simulate surgery, predict myofiber stress and pump function, and quantify changes to regional myocardial material properties. In addition, issues that would need to be resolved in order for computational modeling to play a greater role in cardiac surgery are discussed.

An H∞ strategy for strain estimation in ultrasound elastography using biomechanical modeling constraint

The purpose of ultrasound elastography is to identify lesions by reconstructing the hardness characteristics of tissue reconstructed from ultrasound data. Conventional quasi-static ultrasound elastography is easily applied to obtain axial strain components along the compression direction, with the results inverted to represent the distribution of tissue hardness under the assumption of constant internal stresses. However, previous works of quasi-static ultrasound elastography have found it difficult to obtain the lateral and shear strain components, due to the poor lateral resolution of conventional ultrasound probes. The physical nature of the strain field is a continuous vector field, which should be fully described by the axial, lateral, and shear strain components, and the clinical value of lateral and shear strain components of deformed tissue is gradually being recognized by both engineers and clinicians. Therefore, a biomechanical-model-constrained filtering framework is proposed here for recovering a full displacement field at a high spatial resolution from the noisy ultrasound data. In our implementation, after the biomechanical model constraint is integrated into the state-space equation, both the axial and lateral displacement components can be recovered at a high spatial resolution from the noisy displacement measurements using a robust H∞ filter, which only requires knowledge of the worst-case noise levels in the measurements. All of the strain components can then be calculated by applying a gradient operator to the recovered displacement field. Numerical experiments on synthetic data demonstrated the robustness and effectiveness of our approach, and experiments on phantom data and in-vivo clinical data also produced satisfying results.

Free Tools and Strategies for the Generation of 3D Finite Element Meshes: Modeling of the Cardiac Structures

The Finite Element Method is a well-known technique, being extensively applied in different areas. Studies using the Finite Element Method (FEM) are targeted to improve cardiac ablation procedures. For such simulations, the finite element meshes should consider the size and histological features of the target structures. However, it is possible to verify that some methods or tools used to generate meshes of human body structures are still limited, due to nondetailed models, nontrivial preprocessing, or mainly limitation in the use condition. In this paper, alternatives are demonstrated to solid modeling and automatic generation of highly refined tetrahedral meshes, with quality compatible with other studies focused on mesh generation. The innovations presented here are strategies to integrate Open Source Software (OSS). The chosen techniques and strategies are presented and discussed, considering cardiac structures as a first application context.