A lattice Boltzmann approach to mathematical modeling of myocardial perfusion.

A mathematical model of myocardial perfusion based on the lattice Boltzmann method (LBM) is proposed and its applicability is investigated in both healthy and diseased cases. The myocardium is conceptualized as a porous material in which the transport and mass transfer of a contrast agent in blood flow is studied. The results of myocardial perfusion obtained using LBM in 1D and 2D are confronted with previously reported results in the literature and the results obtained using the mixed-hybrid finite element method. Since LBM is not suitable for simulating flow in heterogeneous porous media, a simplified and computationally efficient 1D-analog approach to 2D diseased case is proposed and its applicability discussed.

Patient-specific closed-loop model of the fontan circulation: Calibration and validation.

The Fontan circulation, designed for managing patients with a single functional ventricle, presents challenges in long-term outcomes. Computational methods offer potential solutions, yet their application in cardiology practice remains largely unexplored. Our aim was to assess the ability of a patient-specific, closed-loop, reduced-order blood flow model to simulate pulsatile blood flow in the Fontan circulation. Using one-dimensional models, we simulated the aorta, superior and inferior venae cavae, and right and left pulmonary arteries, while lumping heart chambers and remaining vessels into zero-dimensional models. The model was calibrated with patient-specific haemodynamic data from combined cardiac catheterisation and magnetic resonance exams, using a novel physics-based stepwise methodology involving simpler open-loop models. Testing on a 10-year-old, anesthetised patient, demonstrated the model's capability to replicate pulsatile pressure and flow in the larger vessels and ventricular pressure. Average relative errors in mean pressure and flow were 2.9 % and 3.6 %, with average relative point-to-point errors (RPPE) in pressure and flow at 5.2 % and 16.0 %. Comparing simulation results to measurements, mean aortic pressure and flow values were 50.7 vs. 50.4 mmHg and 41.6 vs. 41.9 ml/s, respectively, while ventricular pressure values were 28.7 vs. 27.4 mmHg. The model accurately described time-varying ventricular volume with a RPPE of 2.9 %, with mean, minimum, and maximum ventricular volume values for simulation results vs. measurements at 59.2 vs. 58.2 ml, 38.0 vs. 37.6 ml, and 76.0 vs. 74.4 ml, respectively. It provided physiologically realistic predictions of haemodynamic changes from pulmonary vasodilation and atrial fenestration opening. The new model and calibration methodology are freely available, offering a platform to virtually investigate the Fontan circulation's response to clinical interventions and explore potential mechanisms of Fontan failure. Future efforts will concentrate on broadening the model's applicability to a wider range of patient populations and clinical scenarios, as well as testing its effectiveness.

Altered Aortic Hemodynamics and Relative Pressure in Patients with Dilated Cardiomyopathy.

Ventricular-vascular interaction is central in the adaptation to cardiovascular disease. However, cardiomyopathy patients are predominantly monitored using cardiac biomarkers. The aim of this study is therefore to explore aortic function in dilated cardiomyopathy (DCM). Fourteen idiopathic DCM patients and 16 controls underwent cardiac magnetic resonance imaging, with aortic relative pressure derived using physics-based image processing and a virtual cohort utilized to assess the impact of cardiovascular properties on aortic behaviour. Subjects with reduced left ventricular systolic function had significantly reduced aortic relative pressure, increased aortic stiffness, and significantly delayed time-to-pressure peak duration. From the virtual cohort, aortic stiffness and aortic volumetric size were identified as key determinants of aortic relative pressure. As such, this study shows how advanced flow imaging and aortic hemodynamic evaluation could provide novel insights into the manifestation of DCM, with signs of both altered aortic structure and function derived in DCM using our proposed imaging protocol.

Prediction of Ventricular Mechanics After Pulmonary Valve Replacement in Tetralogy of Fallot by Biomechanical Modeling: A Step Towards Precision Healthcare.

Clinical indicators of heart function are often limited in their ability to accurately evaluate the current mechanical state of the myocardium. Biomechanical modeling has been shown to be a promising tool in addition to clinical indicators. By providing a patient-specific measure of myocardial active stress (contractility), biomechanical modeling can enhance the precision of the description of patient's pathophysiology at any given point in time. In this work we aim to explore the ability of biomechanical modeling to predict the response of ventricular mechanics to the progressively decreasing afterload in repaired tetralogy of Fallot (rTOF) patients undergoing pulmonary valve replacement (PVR) for significant residual right ventricular outflow tract obstruction (RVOTO). We used 19 patient-specific models of patients with rTOF prior to pulmonary valve replacement (PVR), denoted as PSMpre, and patient-specific models of the same patients created post-PVR (PSMpost)-both created in our previous published work. Using the PSMpre and assuming cessation of the pulmonary regurgitation and a progressive decrease of RVOT resistance, we built relationships between the contractility and RVOT resistance post-PVR. The predictive value of such in silico obtained relationships were tested against the PSMpost, i.e. the models created from the actual post-PVR datasets. Our results show a linear 1-dimensional relationship between the in silico predicted contractility post-PVR and the RVOT resistance. The predicted contractility was close to the contractility in the PSMpost model with a mean (± SD) difference of 6.5 (± 3.0)%. The relationships between the contractility predicted by in silico PVR vs. RVOT resistance have a potential to inform clinicians about hypothetical mechanical response of the ventricle based on the degree of pre-operative RVOTO.

Right-left ventricular shape variations in tetralogy of Fallot: associations with pulmonary regurgitation.

BACKGROUND: Relationships between right ventricular (RV) and left ventricular (LV) shape and function may be useful in determining optimal timing for pulmonary valve replacement in patients with repaired tetralogy of Fallot (rTOF). However, these are multivariate and difficult to quantify. We aimed to quantify variations in biventricular shape associated with pulmonary regurgitant volume (PRV) in rTOF using a biventricular atlas. METHODS: In this cross-sectional retrospective study, a biventricular shape model was customized to cardiovascular magnetic resonance (CMR) images from 88 rTOF patients (median age 16, inter-quartile range 11.8-24.3 years). Morphometric scores quantifying biventricular shape at end-diastole and end-systole were computed using principal component analysis. Multivariate linear regression was used to quantify biventricular shape associations with PRV, corrected for age, sex, height, and weight. Regional associations were confirmed by univariate correlations with distances and angles computed from the models, as well as global systolic strains computed from changes in arc length from end-diastole to end-systole. RESULTS: PRV was significantly associated with 5 biventricular morphometric scores, independent of covariates, and accounted for 12.3% of total shape variation (p < 0.05). Increasing PRV was associated with RV dilation and basal bulging, in conjunction with decreased LV septal-lateral dimension (LV flattening) and systolic septal motion towards the RV (all p < 0.05). Increased global RV radial, longitudinal, circumferential and LV radial systolic strains were significantly associated with increased PRV (all p < 0.05). CONCLUSION: A biventricular atlas of rTOF patients quantified multivariate relationships between left-right ventricular morphometry and wall motion with pulmonary regurgitation. Regional RV dilation, LV reduction, LV septal-lateral flattening and increased RV strain were all associated with increased pulmonary regurgitant volume. Morphometric scores provide simple metrics linking mechanisms for structural and functional alteration with important clinical indices.

Biomechanical Modeling to Inform Pulmonary Valve Replacement in Tetralogy of Fallot Patients After Complete Repair.

BACKGROUND: A biomechanical model of the heart can be used to incorporate multiple data sources (electrocardiography, imaging, invasive hemodynamics). The purpose of this study was to use this approach in a cohort of patients with tetralogy of Fallot after complete repair (rTOF) to assess comparative influences of residual right ventricular outflow tract obstruction (RVOTO) and pulmonary regurgitation on ventricular health. METHODS: Twenty patients with rTOF who underwent percutaneous pulmonary valve replacement (PVR) and cardiovascular magnetic resonance imaging were included in this retrospective study. Biomechanical models specific to individual patient and physiology (before and after PVR) were created and used to estimate the RV myocardial contractility. The ability of models to capture post-PVR changes of right ventricular (RV) end-diastolic volume (EDV) and effective flow in the pulmonary artery (Qeff) was also compared with expected values. RESULTS: RV contractility before PVR (mean 66 ± 16 kPa, mean ± standard deviation) was increased in patients with rTOF compared with normal RV (38-48 kPa) (P < 0.05). The contractility decreased significantly in all patients after PVR (P < 0.05). Patients with predominantly RVOTO demonstrated greater reduction in contractility (median decrease 35%) after PVR than those with predominant pulmonary regurgitation (median decrease 11%). The model simulated post-PVR decreased EDV for the majority and suggested an increase of Qeff-both in line with published data. CONCLUSIONS: This study used a biomechanical model to synthesize multiple clinical inputs and give an insight into RV health. Individualized modeling allows us to predict the RV response to PVR. Initial data suggest that residual RVOTO imposes greater ventricular work than isolated pulmonary regurgitation.

Left Ventricular Torsion Obtained Using Equilibrated Warping in Patients with Repaired Tetralogy of Fallot.

Patients after surgical repair of Tetralogy of Fallot (rTOF) may suffer a decrease in left ventricular (LV) function. The aim of our study is to evaluate a novel method of assessing LV torsion in patients with rTOF, as an early indicator of systolic LV dysfunction. Motion tracking based on image registration regularized by the equilibrium gap principle, known as equilibrated warping, was employed to assess LV torsion. Seventy-six cases of rTOF and ten normal controls were included. The group of controls was assessed for reproducibility using both equilibrated warping and standard clinical tissue tracking software (CVI42, version 5.10.1, Calgary, Canada). Patients were dichotomized into two groups: normal vs. loss of torsion. Torsion by equilibrated warping was successfully obtained in 68 of 76 (89%) patients and 9 of 10 (90%) controls. For equilibrated warping, the intra- and interobserver coefficients of variation were 0.095 and 0.117, respectively, compared to 0.260 and 0.831 for tissue tracking by standard clinical software. The intra- and inter-observer intraclass correlation coefficients for equilibrated warping were 0.862 and 0.831, respectively, compared to 0.992 and 0.648 for tissue tracking. Loss of torsion was noted in 32 of the 68 (47%) patients with rTOF. There was no difference in LV or RV volumes or ejection fraction between these groups. The assessment of LV torsion by equilibrated warping is feasible and shows good reliability. Loss of torsion is common in patients with rTOF and its robust assessment might contribute into uncovering heart failure in an earlier stage.

HARP-I: A Harmonic Phase Interpolation Method for the Estimation of Motion From Tagged MR Images.

We proposed a novel method called HARP-I, which enhances the estimation of motion from tagged Magnetic Resonance Imaging (MRI). The harmonic phase of the images is unwrapped and treated as noisy measurements of reference coordinates on a deformed domain, obtaining motion with high accuracy using Radial Basis Functions interpolations. Results were compared against Shortest Path HARP Refinement (SP-HR) and Sine-wave Modeling (SinMod), two harmonic image-based techniques for motion estimation from tagged images. HARP-I showed a favorable similarity with both methods under noise-free conditions, whereas a more robust performance was found in the presence of noise. Cardiac strain was better estimated using HARP-I at almost any motion level, giving strain maps with less artifacts. Additionally, HARP-I showed better temporal consistency as a new method was developed to fix phase jumps between frames. In conclusion, HARP-I showed to be a robust method for the estimation of motion and strain under ideal and non-ideal conditions.

Monitoring of cardiovascular physiology augmented by a patient-specific biomechanical model during general anesthesia. A proof of concept study.

During general anesthesia (GA), direct analysis of arterial pressure or aortic flow waveforms may be inconclusive in complex situations. Patient-specific biomechanical models, based on data obtained during GA and capable to perform fast simulations of cardiac cycles, have the potential to augment hemodynamic monitoring. Such models allow to simulate Pressure-Volume (PV) loops and estimate functional indicators of cardiovascular (CV) system, e.g. ventricular-arterial coupling (Vva), cardiac efficiency (CE) or myocardial contractility, evolving throughout GA. In this prospective observational study, we created patient-specific biomechanical models of heart and vasculature of a reduced geometric complexity for n = 45 patients undergoing GA, while using transthoracic echocardiography and aortic pressure and flow signals acquired in the beginning of GA (baseline condition). If intraoperative hypotension (IOH) appeared, diluted norepinephrine (NOR) was administered and the model readjusted according to the measured aortic pressure and flow signals. Such patients were a posteriori assigned into a so-called hypotensive group. The accuracy of simulated mean aortic pressure (MAP) and stroke volume (SV) at baseline were in accordance with the guidelines for the validation of new devices or reference measurement methods in all patients. After NOR administration in the hypotensive group, the percentage of concordance with 10% exclusion zone between measurement and simulation was >95% for both MAP and SV. The modeling results showed a decreased Vva (0.64±0.37 vs 0.88±0.43; p = 0.039) and an increased CE (0.8±0.1 vs 0.73±0.11; p = 0.042) in hypotensive vs normotensive patients. Furthermore, Vva increased by 92±101%, CE decreased by 13±11% (p < 0.001 for both) and contractility increased by 14±11% (p = 0.002) in the hypotensive group post-NOR administration. In this work we demonstrated the application of fast-running patient-specific biophysical models to estimate PV loops and functional indicators of CV system using clinical data available during GA. The work paves the way for model-augmented hemodynamic monitoring at operating theatres or intensive care units to enhance the information on patient-specific physiology.

Investigation of phase-contrast magnetic resonance imaging underestimation of turbulent flow through the aortic valve phantom: experimental and computational study using lattice Boltzmann method.

OBJECTIVE: The accuracy of phase-contrast magnetic resonance imaging (PC-MRI) measurement is investigated using a computational fluid dynamics (CFD) model with the objective to determine the magnitude of the flow underestimation due to turbulence behind a narrowed valve in a phantom experiment. MATERIALS AND METHODS: An acrylic stationary flow phantom is used with three insertable plates mimicking aortic valvular stenoses of varying degrees. Positive and negative horizontal fluxes are measured at equidistant slices using standard PC-MRI sequences by 1.5T and 3T systems. The CFD model is based on the 3D lattice Boltzmann method (LBM). The experimental and simulated data are compared using the Bland-Altman-derived limits of agreement. Based on the LBM results, the turbulence is quantified and confronted with the level of flow underestimation. RESULTS: LBM gives comparable results to PC-MRI for valves up to moderate stenosis on both field strengths. The flow magnitude through a severely stenotic valve was underestimated due to signal void in the regions of turbulent flow behind the valve, consistently with the level of quantified turbulence intensity. DISCUSSION: Flow measured by PC-MRI is affected by noise and turbulence. LBM can simulate turbulent flow efficiently and accurately, it has therefore the potential to improve clinical interpretation of PC-MRI.

Dobutamine stress testing in patients with Fontan circulation augmented by biomechanical modeling.

Understanding (patho)physiological phenomena and mechanisms of failure in patients with Fontan circulation-a surgically established circulation for patients born with a functionally single ventricle-remains challenging due to the complex hemodynamics and high inter-patient variations in anatomy and function. In this work, we present a biomechanical model of the heart and circulation to augment the diagnostic evaluation of Fontan patients with early-stage heart failure. The proposed framework employs a reduced-order model of heart coupled with a simplified circulation including venous return, creating a closed-loop system. We deploy this framework to augment the information from data obtained during combined cardiac catheterization and magnetic resonance exams (XMR), performed at rest and during dobutamine stress in 9 children with Fontan circulation and 2 biventricular controls. We demonstrate that our modeling framework enables patient-specific investigation of myocardial stiffness, contractility at rest, contractile reserve during stress and changes in vascular resistance. Hereby, the model allows to identify key factors underlying the pathophysiological response to stress in these patients. In addition, the rapid personalization of the model to patient data and fast simulation of cardiac cycles make our framework directly applicable in a clinical workflow. We conclude that the proposed modeling framework is a valuable addition to the current clinical diagnostic XMR exam that helps to explain patient-specific stress hemodynamics and can identify potential mechanisms of failure in patients with Fontan circulation.

The importance of the pericardium for cardiac biomechanics: from physiology to computational modeling.

The human heart is enclosed in the pericardial cavity. The pericardium consists of a layered thin sac and is separated from the myocardium by a thin film of fluid. It provides a fixture in space and frictionless sliding of the myocardium. The influence of the pericardium is essential for predictive mechanical simulations of the heart. However, there is no consensus on physiologically correct and computationally tractable pericardial boundary conditions. Here, we propose to model the pericardial influence as a parallel spring and dashpot acting in normal direction to the epicardium. Using a four-chamber geometry, we compare a model with pericardial boundary conditions to a model with fixated apex. The influence of pericardial stiffness is demonstrated in a parametric study. Comparing simulation results to measurements from cine magnetic resonance imaging reveals that adding pericardial boundary conditions yields a better approximation with respect to atrioventricular plane displacement, atrial filling, and overall spatial approximation error. We demonstrate that this simple model of pericardial-myocardial interaction can correctly predict the pumping mechanisms of the heart as previously assessed in clinical studies. Utilizing a pericardial model not only can provide much more realistic cardiac mechanics simulations but also allows new insights into pericardial-myocardial interaction which cannot be assessed in clinical measurements yet.

An adaptive hybridizable discontinuous Galerkin approach for cardiac electrophysiology.

Cardiac electrophysiology simulations are numerically challenging because of the propagation of a steep electrochemical wave front and thus require discretizations with small mesh sizes to obtain accurate results. In this work, we present an approach based on the hybridizable discontinuous Galerkin method (HDG), which allows an efficient implementation of high-order discretizations into a computational framework. In particular, using the advantage of the discontinuous function space, we present an efficient p-adaptive strategy for accurately tracking the wave front. The HDG allows to reduce the overall degrees of freedom in the final linear system to those only on the element interfaces. Additionally, we propose a rule for a suitable integration accuracy for the ionic current term depending on the polynomial order and the cell model to handle high-order polynomials. Our results show that for the same number of degrees of freedom, coarse high-order elements provide more accurate results than fine low-order elements. Introducing p-adaptivity further reduces computational costs while maintaining accuracy by restricting the use of high-order elements to resolve the wave front. For a patient-specific simulation of a cardiac cycle, p-adaptivity reduces the average number of degrees of freedom by 95% compared to the nonadaptive model. In addition to reducing computational costs, using coarse meshes with our p-adaptive high-order HDG method also simplifies practical aspects of mesh generation and postprocessing.

Exploring kinetic energy as a new marker of cardiac function in the single ventricle circulation.

Ventricular volumetric ejection fraction (VV EF) is often normal in patients with single ventricle circulations despite them experiencing symptoms related to circulatory failure. We sought to determine if kinetic energy (KE) could be a better marker of ventricular performance. KE was prospectively quantified using four-dimensional flow MRI in 41 patients with a single ventricle circulation (aged 0.5-28 yr) and compared with 43 healthy volunteers (aged 1.5-62 yr) and 14 patients with left ventricular (LV) dysfunction (aged 28-79 yr). Intraventricular end-diastolic blood was tracked through systole and divided into ejected and residual blood components. Two ejection fraction (EF) metrics were devised based on the KE of the ejected component over the total of both the ejected and residual components using 1) instantaneous peak KE to assess KE EF or 2) summating individual peak particle energy (PE) to assess PE EF. KE EF and PE EF had a smaller range than VV EF in healthy subjects (97.9 ± 0.8 vs. 97.3 ± 0.8 vs. 60.1 ± 5.2%). LV dysfunction caused a fall in KE EF ( P = 0.01) and PE EF ( P = 0.0001). VV EF in healthy LVs and single ventricle hearts was equivalent; however, KE EF and PE EF were lower ( P < 0.001) with a wider range indicating a spectrum of severity. Those reporting the greatest symptomatic impairment (New York Heart Association II) had lower PE EF than asymptomatic subjects ( P = 0.0067). KE metrics are markers of healthy cardiac function. PE EF may be useful in grading dysfunction. NEW & NOTEWORTHY Kinetic energy (KE) represents the useful work of the heart in ejecting blood. This article details the utilization of KE indexes to assess cardiac function in health and a variety of pathophysiological conditions. KE ejection fraction and particle energy ejection fraction (PE EF) showed a narrow range in health and a lower wider range in disease representing a spectrum of severity. PE EF was altered by functional status potentially offering the opportunity to grade dysfunction.

Non-invasive Model-Based Assessment of Passive Left-Ventricular Myocardial Stiffness in Healthy Subjects and in Patients with Non-ischemic Dilated Cardiomyopathy.

Patient-specific modelling has emerged as a tool for studying heart function, demonstrating the potential to provide non-invasive estimates of tissue passive stiffness. However, reliable use of model-derived stiffness requires sufficient model accuracy and unique estimation of model parameters. In this paper we present personalised models of cardiac mechanics, focusing on improving model accuracy, while ensuring unique parametrisation. The influence of principal model uncertainties on accuracy and parameter identifiability was systematically assessed in a group of patients with dilated cardiomyopathy ([Formula: see text]) and healthy volunteers ([Formula: see text]). For all cases, we examined three circumferentially symmetric fibre distributions and two epicardial boundary conditions. Our results demonstrated the ability of data-derived boundary conditions to improve model accuracy and highlighted the influence of the assumed fibre distribution on both model fidelity and stiffness estimates. The model personalisation pipeline-based strictly on non-invasive data-produced unique parameter estimates and satisfactory model errors for all cases, supporting the selected model assumptions. The thorough analysis performed enabled the comparison of passive parameters between volunteers and dilated cardiomyopathy patients, illustrating elevated stiffness in diseased hearts.

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.

Age-related changes in intraventricular kinetic energy: a physiological or pathological adaptation?

Aging has important deleterious effects on the cardiovascular system. We sought to compare intraventricular kinetic energy (KE) in healthy subjects of varying ages with subjects with ventricular dysfunction to understand if changes in energetic momentum may predispose individuals to heart failure. Four-dimensional flow MRI was acquired in 35 healthy subjects (age: 1-67 yr) and 10 patients with left ventricular (LV) dysfunction (age: 28-79 yr). Healthy subjects were divided into age quartiles (1st quartile: <16 yr, 2nd quartile: 17-32 yr, 3rd quartile: 33-48 yr, and 4th quartile: 49-64 yr). KE was measured in the LV throughout the cardiac cycle and indexed to ventricular volume. In healthy subjects, two large peaks corresponding to systole and early diastole occurred during the cardiac cycle. A third smaller peak was seen during late diastole in eight adults. Systolic KE (P = 0.182) and ejection fraction (P = 0.921) were preserved through all age groups. Older adults showed a lower early peak diastolic KE compared with children (P < 0.0001) and young adults (P = 0.025). Subjects with LV dysfunction had reduced ejection fraction (P < 0.001) and compared with older healthy adults exhibited a similar early peak diastolic KE (P = 0.142) but with the addition of an elevated KE in diastasis (P = 0.029). In healthy individuals, peak diastolic KE progressively decreases with age, whereas systolic peaks remain constant. Peak diastolic KE in the oldest subjects is comparable to those with LV dysfunction. Unique age-related changes in ventricular diastolic energetics might be physiological or herald subclinical pathology.

Estimation of passive and active properties in the human heart using 3D tagged MRI.

Advances in medical imaging and image processing are paving the way for personalised cardiac biomechanical modelling. Models provide the capacity to relate kinematics to dynamics and-through patient-specific modelling-derived material parameters to underlying cardiac muscle pathologies. However, for clinical utility to be achieved, model-based analyses mandate robust model selection and parameterisation. In this paper, we introduce a patient-specific biomechanical model for the left ventricle aiming to balance model fidelity with parameter identifiability. Using non-invasive data and common clinical surrogates, we illustrate unique identifiability of passive and active parameters over the full cardiac cycle. Identifiability and accuracy of the estimates in the presence of controlled noise are verified with a number of in silico datasets. Unique parametrisation is then obtained for three datasets acquired in vivo. The model predictions show good agreement with the data extracted from the images providing a pipeline for personalised biomechanical analysis.