Patient specific fluid–structure ventricular modelling for integrated cardiac care

Acute vasodilator response testing in the adult Fontan circulation using non-invasive 4D Flow MRI: a proof-of-principle study

BACKGROUND

Pulmonary vasodilator therapy in Fontan patients can improve exercise tolerance. We aimed to assess the potential for non-invasive testing of acute vasodilator response using four-dimensional (D) flow MRI during oxygen inhalation.

MATERIALS AND METHODS

Six patients with well-functioning Fontan circulations were prospectively recruited and underwent cardiac MRI. Ventricular anatomical imaging and 4D Flow MRI were acquired at baseline and during inhalation of oxygen. Data were compared with six age-matched healthy volunteers with 4D Flow MRI scans acquired at baseline.

RESULTS

All six patients tolerated the MRI scan well. The dominant ventricle had a left ventricular morphology in all cases. On 4D Flow MRI assessment, two patients (Patients 2 and 6) showed improved cardiac filling with improved preload during oxygen administration, increased mitral inflow, increased maximum E-wave kinetic energy, and decreased systolic peak kinetic energy. Patient 1 showed improved preload only. Patient 5 showed no change, and patient 3 had equivocal results. Patient 4, however, showed a decrease in preload and cardiac filling/function with oxygen.

DISCUSSION

Using oxygen as a pulmonary vasodilator to assess increased pulmonary venous return as a marker for positive acute vasodilator response would provide pre-treatment assessment in a more physiological state - the awake patient. This proof-of-concept study showed that it is well tolerated and has shown changes in some stable patients with a Fontan circulation.

A partition of unity approach to fluid mechanics and fluid-structure interaction

For problems involving large deformations of thin structures, simulating fluid-structure interaction (FSI) remains a computationally expensive endeavour which continues to drive interest in the development of novel approaches. Overlapping domain techniques have been introduced as a way to combine the fluid-solid mesh conformity, seen in moving-mesh methods, without the need for mesh smoothing or re-meshing, which is a core characteristic of fixed mesh approaches. In this work, we introduce a novel overlapping domain method based on a partition of unity approach. Unified function spaces are defined as a weighted sum of fields given on two overlapping meshes. The method is shown to achieve optimal convergence rates and to be stable for steady-state Stokes, Navier-Stokes, and ALE Navier-Stokes problems. Finally, we present results for FSI in the case of 2D flow past an elastic beam simulation. These initial results point to the potential applicability of the method to a wide range of FSI applications, enabling boundary layer refinement and large deformations without the need for re-meshing or user-defined stabilization.

Sensitivity Analysis of a Left Ventricle Model in the Context of Intraventricular Dyssynchrony

The objective of the current study was to propose a sensitivity analysis of a 3D left ventricle model in order to assess the influence of parameters on myocardial mechanical dispersion. A finite element model of LV electro-mechanical activity was proposed and a screening method was used to evaluate the sensitivity of model parameters on the standard deviation of time to peak strain. Results highlight the importance of propagation parameters associated with septal and lateral segments activation. Simulated curves were compared to myocardial strains, obtained from echocardiography of one healthy subject and one patient diagnosed with intraventricular dyssynchrony and coronary artery disease. Results show a close match between simulation and clinical strains and illustrate the model ability to reproduce myocardial strains in the context of intraventricular dyssynchrony.

Fully coupled fluid-electro-mechanical model of the human heart for supercomputers

In this work, we present a fully coupled fluid-electro-mechanical model of a 50th percentile human heart. The model is implemented on Alya, the BSC multi-physics parallel code, capable of running efficiently in supercomputers. Blood in the cardiac cavities is modeled by the incompressible Navier-Stokes equations and an arbitrary Lagrangian-Eulerian (ALE) scheme. Electrophysiology is modeled with a monodomain scheme and the O'Hara-Rudy cell model. Solid mechanics is modeled with a total Lagrangian formulation for discrete strains using the Holzapfel-Ogden cardiac tissue material model. The three problems are simultaneously and bidirectionally coupled through an electromechanical feedback and a fluid-structure interaction scheme. In this paper, we present the scheme in detail and propose it as a computational cardiac workbench.

In vivo estimation of passive biomechanical properties of human myocardium

Identification of in vivo passive biomechanical properties of healthy human myocardium from regular clinical data is essential for subject-specific modelling of left ventricle (LV). In this work, myocardium was defined by Holzapfel-Ogden constitutive law. Therefore, the objectives of the study were (a) to estimate the ranges of the constitutive parameters for healthy human myocardium using non-invasive routine clinical data, and (b) to investigate the effect of geometry, LV end-diastolic pressure (EDP) and fibre orientations on estimated values. In order to avoid invasive measurements and additional scans, LV cavity volume, measured from routine MRI, and empirical pressure-normalised-volume relation (Klotz-curve) were used as clinical data. Finite element modelling, response surface method and genetic algorithm were used to inversely estimate the constitutive parameters. Due to the ill-posed nature of the inverse optimisation problem, the myocardial properties was extracted by identifying the ranges of the parameters, instead of finding unique values. Additional sensitivity studies were carried out to identify the effect of LV EDP, fibre orientation and geometry on estimated parameters. Although uniqueness of the solution cannot be achieved, the normal ranges of the parameters produced similar mechanical responses within the physiological ranges. These information could be used in future computational studies for designing heart failure treatments. Graphical abstract.

In vivo estimation of passive biomechanical properties of human myocardium

Identification of in vivo passive biomechanical properties of healthy human myocardium from regular clinical data is essential for subject-specific modelling of left ventricle (LV). In this work, myocardium was defined by Holzapfel-Ogden constitutive law. Therefore, the objectives of the study were (a) to estimate the ranges of the constitutive parameters for healthy human myocardium using non-invasive routine clinical data, and (b) to investigate the effect of geometry, LV end-diastolic pressure (EDP) and fibre orientations on estimated values. In order to avoid invasive measurements and additional scans, LV cavity volume, measured from routine MRI, and empirical pressure-normalised-volume relation (Klotz-curve) were used as clinical data. Finite element modelling, response surface method and genetic algorithm were used to inversely estimate the constitutive parameters. Due to the ill-posed nature of the inverse optimisation problem, the myocardial properties was extracted by identifying the ranges of the parameters, instead of finding unique values. Additional sensitivity studies were carried out to identify the effect of LV EDP, fibre orientation and geometry on estimated parameters. Although uniqueness of the solution cannot be achieved, the normal ranges of the parameters produced similar mechanical responses within the physiological ranges. These information could be used in future computational studies for designing heart failure treatments.

Patient-specific modeling of left ventricular electromechanics as a driver for haemodynamic analysis

Aims

Models of blood flow in the left ventricle (LV) and aorta are an important tool for analysing the interplay between LV deformation and flow patterns. Typically, image-based kinematic models describing endocardial motion are used as an input to blood flow simulations. While such models are suitable for analysing the hemodynamic status quo, they are limited in predicting the response to interventions that alter afterload conditions. Mechano-fluidic models using biophysically detailed electromechanical (EM) models have the potential to overcome this limitation, but are more costly to build and compute. We report our recent advancements in developing an automated workflow for the creation of such CFD ready kinematic models to serve as drivers of blood flow simulations.

Methods and results

EM models of the LV and aortic root were created for four pediatric patients treated for either aortic coarctation or aortic valve disease. Using MRI, ECG and invasive pressure recordings, anatomy as well as electrophysiological, mechanical and circulatory model components were personalized.

Results

The implemented modeling pipeline was highly automated and allowed model construction and execution of simulations of a patient’s heartbeat within 1 day. All models reproduced clinical data with acceptable accuracy.

Conclusion

Using the developed modeling workflow, the use of EM LV models as driver of fluid flow simulations is becoming feasible. While EM models are costly to construct, they constitute an important and nontrivial step towards fully coupled electro-mechano-fluidic (EMF) models and show promise as a tool for predicting the response to interventions which affect afterload conditions.

Two-dimensional intraventricular flow pattern visualization using the image-based computational fluid dynamics

The image-based computational fluid dynamics (IB-CFD) technique, as the combination of medical images and the CFD method, is utilized in this research to analyze the left ventricle (LV) hemodynamics. The research primarily aims to propose a semi-automated technique utilizing some freely available and commercial software packages in order to simulate the LV hemodynamics using the IB-CFD technique. In this research, moreover, two different physiological time-resolved 2D models of a patient-specific LV with two different types of aortic and mitral valves, including the orifice-type valves and integrated with rigid leaflets, are adopted to visualize the process of developing intraventricular vortex formation and propagation. The blood flow pattern over the whole cardiac cycle of two models is also compared to investigate the effect of utilizing different valve types in the process of the intraventricular vortex formation. Numerical findings indicate that the model with integrated valves can predict more complex intraventricular flow that can match better the physiological flow pattern in comparison to the orifice-type model.

Heart blood flow simulation: a perspective review

Cardiovascular disease (CVD), the leading cause of death today, incorporates a wide range of cardiovascular system malfunctions that affect heart functionality. It is believed that the hemodynamic loads exerted on the cardiovascular system, the left ventricle (LV) in particular, are the leading cause of CVD initiation and propagation. Moreover, it is believed that the diagnosis and prognosis of CVD at an early stage could reduce its high mortality and morbidity rate. Therefore, a set of robust clinical cardiovascular assessment tools has been introduced to compute the cardiovascular hemodynamics in order to provide useful insights to physicians to recognize indicators leading to CVD and also to aid the diagnosis of CVD. Recently, a combination of computational fluid dynamics (CFD) and different medical imaging tools, image-based CFD (IB-CFD), has been widely employed for cardiovascular functional assessment by providing reliable hemodynamic parameters. Even though the capability of CFD to provide reliable flow dynamics in general fluid mechanics problems has been widely demonstrated for many years, up to now, the clinical implications of the IB-CFD patient-specific LVs have not been applicable due to its limitations and complications. In this paper, we review investigations conducted to numerically simulate patient-specific human LV over the past 15 years using IB-CFD methods. Firstly, we divide different studies according to the different LV types (physiological and different pathological conditions) that have been chosen to reconstruct the geometry, and then discuss their contributions, methodologies, limitations, and findings. In this regard, we have studied CFD simulations of intraventricular flows and related cardiology insights, for (i) Physiological patient-specific LV models, (ii) Pathological heart patient-specific models, including myocardial infarction, dilated cardiomyopathy, hypertrophic cardiomyopathy and hypoplastic left heart syndrome. Finally, we discuss the current stage of the IB-CFD LV simulations in order to mimic realistic hemodynamics of patient-specific LVs. We can conclude that heart flow simulation is on the right track for developing into a useful clinical tool for heart function assessment, by (i) incorporating most of heart structures' (such as heart valves) operations, and (ii) providing useful diagnostic indices based hemodynamic parameters, for routine adoption in clinical usage.

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.

A multiscale model for the study of cardiac biomechanics in single-ventricle surgeries: a clinical case

Complex congenital heart disease characterized by the underdevelopment of one ventricular chamber (single ventricle (SV) circulation) is normally treated with a three-stage surgical repair. This study aims at developing a multiscale computational framework able to couple a patient-specific three-dimensional finite-element model of the SV to a patient-specific lumped parameter (LP) model of the whole circulation, in a closed-loop fashion. A sequential approach was carried out: (i) cardiocirculatory parameters were estimated by using a fully LP model; (ii) ventricular material parameters and unloaded geometry were identified by means of the stand-alone, three-dimensional model of the SV; and (iii) the three-dimensional model of SV was coupled to the LP model of the circulation, thus closing the loop and creating a multiscale model. Once the patient-specific multiscale model was set using pre-operative clinical data, the virtual surgery was performed, and the post-operative conditions were simulated. This approach allows the analysis of local information on ventricular function as well as global parameters of the cardiovascular system. This methodology is generally applicable to patients suffering from SV disease for surgical planning at different stages of treatment. As an example, a clinical case from stage 1 to stage 2 is considered here.

Using 4D Cardiovascular Magnetic Resonance Imaging to Validate Computational Fluid Dynamics: A Case Study

Computational fluid dynamics (CFD) can have a complementary predictive role alongside the exquisite visualization capabilities of 4D cardiovascular magnetic resonance (CMR) imaging. In order to exploit these capabilities (e.g., for decision-making), it is necessary to validate computational models against real world data. In this study, we sought to acquire 4D CMR flow data in a controllable, experimental setup and use these data to validate a corresponding computational model. We applied this paradigm to a case of congenital heart disease, namely, transposition of the great arteries (TGA) repaired with arterial switch operation. For this purpose, a mock circulatory loop compatible with the CMR environment was constructed and two detailed aortic 3D models (i.e., one TGA case and one normal aortic anatomy) were tested under realistic hemodynamic conditions, acquiring 4D CMR flow. The same 3D domains were used for multi-scale CFD simulations, whereby the remainder of the mock circulatory system was appropriately summarized with a lumped parameter network. Boundary conditions of the simulations mirrored those measured in vitro. Results showed a very good quantitative agreement between experimental and computational models in terms of pressure (overall maximum % error = 4.4% aortic pressure in the control anatomy) and flow distribution data (overall maximum % error = 3.6% at the subclavian artery outlet of the TGA model). Very good qualitative agreement could also be appreciated in terms of streamlines, throughout the cardiac cycle. Additionally, velocity vectors in the ascending aorta revealed less symmetrical flow in the TGA model, which also exhibited higher wall shear stress in the anterior ascending aorta.

Implicit Partitioned Cardiovascular Fluid-Structure Interaction of the Heart Cycle Using Non-newtonian Fluid Properties and Orthotropic Material Behavior

Although image-based methods like MRI are well-developed, numerical simulation can help to understand human heart function. This function results from a complex interplay of biochemistry, structural mechanics, and blood flow. The complexity of the entire system often causes one of the three parts to be neglected, which limits the truth to reality of the reduced model. This paper focuses on the interaction of myocardial stress distribution and ventricular blood flow during diastole and systole in comparison to a simulation of the same patient-specific geometry with a given wall movement (Spiegel, Strömungsmechanischer Beitrag zur Planung von Herzoperationen, 2009). The orthotropic constitutive law proposed by Holzapfel et al. (Philos. Trans. R. Soc. Lond. Ser. A, 367:3445-3475, 2009) was implemented in a finite element package to model the passive behavior of the myocardium. Then, this law was modified for contraction. Via the ALE method, the structural model was coupled to a flow model which incorporates blood rheology and the circulatory system (Oertel, Prandtl-Essentials of Fluid Mechanics, 3rd edn, Springer Science + Business Media, 2010; Oertel et al., Modelling the Human Cardiac Fluid Mechanics, 3rd edn, Universitätsverlag Karlsruhe, 2009). Comparison reveals a good quantitative and qualitative agreement with respect to fluid flow. The motion of the myocardium is consistent with physiological observations. The calculated stresses and the distribution are within the physiological range and appear to be reasonable. The coupled model presented contains many features essential to cardiac function. It is possible to calculate wall stresses as well as the characteristic ventricular fluid flow. Based on the simulations we derive two characteristics to assess the health state quantitatively including solid and fluid mechanical aspects.

Quasi-static image-based immersed boundary-finite element model of left ventricle under diastolic loading

Finite stress and strain analyses of the heart provide insight into the biomechanics of myocardial function and dysfunction. Herein, we describe progress toward dynamic patient-specific models of the left ventricle using an immersed boundary (IB) method with a finite element (FE) structural mechanics model. We use a structure-based hyperelastic strain-energy function to describe the passive mechanics of the ventricular myocardium, a realistic anatomical geometry reconstructed from clinical magnetic resonance images of a healthy human heart, and a rule-based fiber architecture. Numerical predictions of this IB/FE model are compared with results obtained by a commercial FE solver. We demonstrate that the IB/FE model yields results that are in good agreement with those of the conventional FE model under diastolic loading conditions, and the predictions of the LV model using either numerical method are shown to be consistent with previous computational and experimental data. These results are among the first to analyze the stress and strain predictions of IB models of ventricular mechanics, and they serve both to verify the IB/FE simulation framework and to validate the IB/FE model. Moreover, this work represents an important step toward using such models for fully dynamic fluid-structure interaction simulations of the heart. © 2014 The Authors. International Journal for Numerical Methods in Engineering published by John Wiley & Sons, Ltd.

Congenital left ventricular outpouchings: a systematic review of 839 cases and introduction of a novel classification after two centuries

BACKGROUND

Congenital left ventricular outpouchings (LVOs) are reported under five overlapping and poorly defined terms including left ventricular accessory chamber, left ventricular aneurysm (LVA), left ventricular diverticulum (LVD), double-chambered LV, and accessory left ventricle. Diagnostic criteria are frequently mixed and not mutually exclusive. They convey no information regarding treatment strategy and prognosis.

OBJECTIVES

The aim of this systematic review is to provide a clear and inclusive classification, with therapeutic and prognostic implications, for congenital LVOs.

DATA SOURCES

We performed three separate sets of search on three subjects including "congenital left ventricular outpouchings," "important and simply measurable markers of left ventricular function," and "relationship of mechanics of intraventricular blood flow and optimal vortex formation in left ventricle and elliptical geometry of LV."

STUDY ELIGIBILITY CRITERIA

We enrolled case series, review articles, and case reports with literature review. All types of acquired LVO's were excluded.

STUDY APPRAISAL AND SYNTHESIS METHODS

We studied the abstracts of all searched articles. We focused on diagnostic criteria and patients' outcome. To examine the validity and reliability of the novel classification, fifteen previous studies were revisited using the novel classification.

RESULTS

A total of 20 papers from 11 countries fulfilled our inclusion criteria. The age of patients ranged from prenatal age to geriatric age range. Diagnostic criteria were clearly stated only for two of the above five terms (i.e., congenital LVA and congenital LVD). Cases with mixed diagnostic criteria were frequent.Elliptical geometry of left ventricle was found to have significant impact on effective blood flow mechanics in LV. A simple inclusive classification for congenital LVOs, with therapeutic and prognostic implications, was introduced.

CONCLUSION

The cornerstone of this classification is elliptical LV geometry. Large-type IIc LVO have dismal prognosis, if left untreated. LVO type I and small LVO type IIa have the best prognosis.