Virtual surgeries in patients with congenital heart disease: a multi-scale modelling test case

Passive performance evaluation and validation of a viscous impeller pump for subpulmonary fontan circulatory support

Patients with single ventricle defects undergoing the Fontan procedure eventually face Fontan failure. Long-term cavopulmonary assist devices using rotary pump technologies are currently being developed as a subpulmonary power source to prevent and treat Fontan failure. Low hydraulic resistance is a critical safety requirement in the event of pump failure (0 RPM) as a modest 2 mmHg cavopulmonary pressure drop can compromise patient hemodynamics. The goal of this study is therefore to assess the passive performance of a viscous impeller pump (VIP) we are developing for Fontan patients, and validate flow simulations against in-vitro data. Two different blade heights (1.09 mm vs 1.62 mm) and a blank housing model were tested using a mock circulatory loop (MCL) with cardiac output ranging from 3 to 11 L/min. Three-dimensional flow simulations were performed and compared against MCL data. In-silico and MCL results demonstrated a pressure drop of < 2 mmHg at a cardiac output of 7 L/min for both blade heights. There was good agreement between simulation and MCL results for pressure loss (mean difference - 0.23 mmHg 95% CI [0.24-0.71]). Compared to the blank housing model, low wall shear stress area and oscillatory shear index on the pump surface were low, and mean washout times were within 2 s. This study demonstrated the low resistance characteristic of current VIP designs in the failed condition that results in clinically acceptable minimal pressure loss without increased washout time as compared to a blank housing model under normal cardiac output in Fontan patients.

Inverse problems in blood flow modeling: A review

Mathematical and computational modeling of the cardiovascular system is increasingly providing non-invasive alternatives to traditional invasive clinical procedures. Moreover, it has the potential for generating additional diagnostic markers. In blood flow computations, the personalization of spatially distributed (i.e., 3D) models is a key step which relies on the formulation and numerical solution of inverse problems using clinical data, typically medical images for measuring both anatomy and function of the vasculature. In the last years, the development and application of inverse methods has rapidly expanded most likely due to the increased availability of data in clinical centers and the growing interest of modelers and clinicians in collaborating. Therefore, this work aims to provide a wide and comparative overview of literature within the last decade. We review the current state of the art of inverse problems in blood flows, focusing on studies considering fully dimensional fluid and fluid-solid models. The relevant physical models and hemodynamic measurement techniques are introduced, followed by a survey of mathematical data assimilation approaches used to solve different kinds of inverse problems, namely state and parameter estimation. An exhaustive discussion of the literature of the last decade is presented, structured by types of problems, models and available data.

Virtual Transcatheter Interventions for Peripheral Pulmonary Artery Stenosis in Williams and Alagille Syndromes

Background

Despite favorable outcomes of surgical pulmonary artery (PA) reconstruction, isolated proximal stenting of the central PAs is common clinical practice for patients with peripheral PA stenosis in association with Williams and Alagille syndromes. Given the technical challenges of PA reconstruction and the morbidities associated with transcatheter interventions, the hemodynamic consequences of all treatment strategies must be rigorously assessed. Our study aims to model, assess, and predict hemodynamic outcomes of transcatheter interventions in these patients.

Methods and Results

Isolated proximal and “extensive” interventions (stenting and/or balloon angioplasty of proximal and lobar vessels) were performed in silico on 6 patient‐specific PA models. Autoregulatory adaptation of the cardiac output and downstream arterial resistance was modeled in response to intervention‐induced hemodynamic perturbations. Postintervention computational fluid dynamics predictions were validated in 2 stented patients and quantitatively assessed in 4 surgical patients. Our computational methods accurately predicted postinterventional PA pressures, the primary indicators of success for treatment of peripheral PA stenosis. Proximal and extensive treatment achieved median reductions of 14% and 40% in main PA systolic pressure, 27% and 56% in pulmonary vascular resistance, and 10% and 45% in right ventricular stroke work, respectively.

Conclusions

In patients with Williams and Alagille syndromes, extensive transcatheter intervention is required to sufficiently reduce PA pressures and right ventricular stroke work. Transcatheter therapy was shown to be ineffective for long‐segment stenosis and pales hemodynamically in comparison with published outcomes of surgical reconstruction. Regardless of the chosen strategy, a virtual treatment planning platform could identify lesions most critical for optimizing right ventricular afterload.

Individual variability in animal-specific hemodynamic compensation following myocardial infarction

Ventricular enlargement and heart failure are common in patients who survive a myocardial infarction (MI). There is striking variability in the degree of post-infarction ventricular remodeling, however, and no one factor or set of factors have been identified that predicts heart failure risk well. Sympathetic activation directly and indirectly modulates hypertrophic stimuli by altering both neurohormonal milieu and ventricular loading. In a recent study, we developed a method to identify the balance of reflex compensatory mechanisms employed by individual animals following MI based on measured hemodynamics. Here, we conducted prospective studies of acute myocardial infarction in rats to test the degree of variability in reflex compensation as well as whether responses to pharmacologic agents targeted at those reflex mechanisms could be anticipated in individual animals. We found that individual animals use very different mixtures of reflex compensation in response to experimental coronary ligation. Some of these mechanisms were related - animals that compensated strongly with venoconstriction tended to exhibit a decrease in the contractility of the surviving myocardium and those that increased contractility tended to exhibit venodilation. Furthermore, some compensatory mechanisms - such as venoconstriction - increased the extent of predicted ventricular enlargement. Unfortunately, initial reflex responses to infarction were a poor predictor of subsequent responses to pharmacologic agents, suggesting that customizing pharmacologic therapy to individuals based on an initial response will be challenging.

Computational Fluid Dynamics Support for Fontan Planning in Minutes, Not Hours: The Next Step in Clinical Pre-Interventional Simulations

Computational fluid dynamics (CFD) modeling may aid in planning of invasive interventions in Fontan patients. Clinical application of current CFD techniques is however limited by complexity and long computation times. Therefore, we validated a “lean” CFD method to magnetic resonance imaging (MRI) and an “established” CFD method, ultimately aiming to reduce complexity to enable predictive CFD during ongoing interventions. Fifteen Fontan patients underwent MRI for CFD modeling. The differences between lean and established approach, in hepatic and total flow percentage to the left pulmonary artery (%LPA), power loss and relative wall shear stress area were 1.5 ± 4.0%, -0.17 ± 1.1%, -0.055 ± 0.092 mW and 1.1 ± 1.4%. Compared with MRI, the lean and established method showed a bias in %LPA of -1.9 ± 3.4% and -1.8 ± 3.1%. Computation time was for the lean and established approach 3.0 ± 2.0 min and 7.0 ± 3.4 h, respectively. We conclude that the proposed lean method provides fast and reliable results for future CFD support during interventions.

Predictive Models for Pulmonary Artery Size in Fontan Patients

We developed models of pulmonary artery (PA) size in Fontan patients as a function of age and body surface area (BSA) using linear regression and breakpoint analyses based on data from 43 Fontan patients divided into two groups: the extracardiac conduit (ECC) group (n = 24) and the non-ECC group (n = 19). Model predictions were compared against those of a non-Fontan control group (n = 18) and published literature. We observed strong positive correlations of the mean PA diameter with BSA (r = 0.9, p < 0.05) and age (r = 0.88, p < 0.05) in the ECC group. The absolute percentage differences between our BSA and age model predictions against published literature were less than 16% and 20%, respectively. Predicted PA size for Fontan patients was consistently smaller than the control group. These models may serve as useful references for clinicians and be utilized to construct 3D anatomic models that correspond to patient body size or age.

Engineering Perspective on Cardiovascular Simulations of Fontan Hemodynamics: Where Do We Stand with a Look Towards Clinical Application

BACKGROUND

Cardiovascular simulations for patients with single ventricles undergoing the Fontan procedure can assess patient-specific hemodynamics, explore surgical advances, and develop personalized strategies for surgery and patient care. These simulations have not yet been broadly accepted as a routine clinical tool owing to a number of limitations. Numerous approaches have been explored to seek innovative solutions for improving methodologies and eliminating these limitations.

PURPOSE

This article first reviews the current state of cardiovascular simulations of Fontan hemodynamics. Then, it will discuss the technical progress of Fontan simulations with the emphasis of its clinical impact, noting that substantial improvements have been made in the considerations of patient-specific anatomy, flow, and blood rheology. The article concludes with insights into potential future directions involving clinical validation, uncertainty quantification, and computational efficiency. The advancements in these aspects could promote the clinical usage of Fontan simulations, facilitating its integration into routine clinical practice.

Is Doppler Echocardiography Adequate for Surgical Planning of Single Ventricle Patients?

PURPOSE

Surgical planning has shown great potential for optimizing outcomes for patients affected by single ventricle (SV) malformations. Phase-contrast magnetic resonance imaging (PC-MRI) is the routine technique used for flow acquisition in the surgical planning paradigm. However, PC-MRI may suffer from possible artifacts in certain cases; furthermore, this technology may not be readily available for patients in low and lower-middle-income countries. Therefore, this study aims to investigate the effectiveness of using Doppler echocardiography (echo-Doppler) for flow acquisitions of SV surgical planning.

METHODS

This study included eight patients whose blood flow data was acquired by both PC-MRI and echo-Doppler. A virtual surgery platform was used to generate two surgical options for each patient: (1) a traditional Fontan conduit and (2) a Y-graft. Computational fluid dynamics (CFD) simulations were conducted using the two flow acquisitions to assess clinically relevant hemodynamic metrics: indexed power loss (iPL) and hepatic flow distribution (HFD).

RESULTS

Differences exist in flow data acquired by PC-MRI and echo-Doppler, but no statistical significance was obtained. Flow fields, therefore, exhibit discrepancies between simulations using flow acquisitions by PC-MRI and echo-Doppler. In virtual surgery, the two surgical options were ranked based on these metrics. No difference was observed in the ranking of surgical options between using different flow acquisitions.

CONCLUSION

Doppler echocardiography is an adequate alternative approach to acquire flow data for SV surgical planning. This finding encourages broader usage of SV surgical planning with echo-Doppler when MRI may present artifacts or is not available, especially in low and lower-middle-income countries.

Application of physics-based flow models in cardiovascular medicine: Current practices and challenges

Personalized physics-based flow models are becoming increasingly important in cardiovascular medicine. They are a powerful complement to traditional methods of clinical decision-making and offer a wealth of physiological information beyond conventional anatomic viewing using medical imaging data. These models have been used to identify key hemodynamic biomarkers, such as pressure gradient and wall shear stress, which are associated with determining the functional severity of cardiovascular diseases. Importantly, simulation-driven diagnostics can help researchers understand the complex interplay between geometric and fluid dynamic parameters, which can ultimately improve patient outcomes and treatment planning. The possibility to compute and predict diagnostic variables and hemodynamics biomarkers can therefore play a pivotal role in reducing adverse treatment outcomes and accelerate development of novel strategies for cardiovascular disease management.

The impact of compliance on Stage 2 uni-ventricular heart circulation: An experimental assessment of the Bidirectional Glenn

The Bidirectional Glenn (BDG) or cavopulmonary connection is typically undertaken to volume unload the single ventricle in an effort to preserve ventricular and atrioventricular valve function. The geometry of this surgical palliation has been shown to influence the fluid energy loss as well as the distribution of flow that enters through the superior vena cava. In-vitro and in-silico studies to date have been performed on rigid wall models, while this investigation looks at the impact of flexible thin walled models versus rigid walls. Rigid and compliant models of two patient-specific Glenn geometries were fabricated and tested under various flow conditions, within a biosimulator capable of replicating patient specific flow conditions. It was found that the compliant models exhibit greater levels of energy loss compared to the rigid models. Along with these findings greater levels of turbulence was found in both compliant models compared to their rigid counterparts under ultrasound examinations. This shows that vessel compliance has a significant impact on the hemodynamics within hypoplastic left heart syndrome.

Computational fluid dynamics: a primer for congenital heart disease clinicians

Computational fluid dynamics has become an important tool for studying blood flow dynamics. As an in-silico collection of methods, computational fluid dynamics is noninvasive and provides numerical values for the most important parameters of blood flow, such as velocity and pressure that are crucial in hemodynamic studies. In this primer, we briefly explain the basic theory and workflow of the two most commonly applied computational fluid dynamics techniques used in the congenital heart disease literature: the finite element method and the finite volume method. We define important terminology and include specific examples of how using these methods can answer important clinical questions in congenital cardiac surgery planning and perioperative patient management.

Personalized Perioperative Multi-scale, Multi-physics Heart Simulation of Double Outlet Right Ventricle

For treatment of complex congenital heart disease, computer simulation using a three-dimensional heart model may help to improve outcomes by enabling detailed preoperative evaluations. However, no highly integrated model that accurately reproduces a patient's pathophysiology, which is required for this simulation has been reported. We modelled a case of complex congenital heart disease, double outlet right ventricle with ventricular septal defect and atrial septal defect. From preoperative computed tomography images, finite element meshes of the heart and torso were created, and cell model of cardiac electrophysiology and sarcomere dynamics was implemented. The parameter values of the heart model were adjusted to reproduce the patient's electrocardiogram and haemodynamics recorded preoperatively. Two options of in silico surgery were performed using this heart model, and the resulting changes in performance were examined. Preoperative and postoperative simulations showed good agreement with clinical records including haemodynamics and measured oxyhaemoglobin saturations. The use of a detailed sarcomere model also enabled comparison of energetic efficiency between the two surgical options. A novel in silico model of congenital heart disease that integrates molecular models of cardiac function successfully reproduces the observed pathophysiology. The simulation of postoperative state by in silico surgeries can help guide clinical decision-making.

Novel surgical strategy for complicated pulmonary stenosis using haemodynamic analysis based on a virtual operation with numerical flow analysis

OBJECTIVES

A novel surgical strategy using haemodynamic analyses based on virtual operations with computational simulations has been introduced for complicated pulmonary stenosis. We evaluated the efficacy of this strategy.

METHODS

Six patients were enrolled. Before surgery, the optimal pulmonary arteries were constructed based on computational fluid dynamics using 3-dimensional computed tomography. Energy loss (EL, mW) and wall shear stress (WSS, Pa) were calculated. We compared the shapes of preoperative and optimal pulmonary arteries to determine the surgical strategy, including the incision line and the shape of the patch (virtual surgery). EL and WSS were compared between virtual and actual surgeries using flow analysis.

RESULTS

In both the virtual and actual surgeries, postoperative EL tended to be lower than the preoperative EL, although there were no significant differences (P = 0.12 and P = 0.17, respectively). The mean WSS in the virtual surgery was significantly reduced from 112 ± 130 Pa to 25 ± 24 Pa (P = 0.028). After the actual surgery, the mean WSS was also significantly reduced to 30 ± 23 Pa (P = 0.047). There were no significant differences in the values for EL and WSS before and after surgery or between virtual and actual surgery (P = 0.94 and P = 0.85, respectively).

CONCLUSIONS

Pulmonary artery plasty, using computational fluid dynamics based on virtual surgery, is an efficient surgical strategy. This novel strategy can easily and successfully provide an optimal pulmonary artery plasty equivalent to that using the conventional approach, which is based on the surgeon's personal experience and judgement.

Cognitive Imaging

This article, an extended version of ICCI*CC-2017 paper, co-authored by biomedical engineers specializing in brain blood circulation modeling and by experts in meaning-based NLP. This article suggests a cognitive computing technology for medical imaging analysis that removes image artifacts resulting in visual deviations from reality, such as discontinuous blood vessels or two vessels shown merged when they are not. It is implemented by supplying the pertinent knowledge that humans have to the computer and letting it initiate the corrective post-processing. The existing OST resource is centered on the ontology that is made to accommodate the domain with a minor adjustment effort; however, any ontology can be used, as demonstrated in this article. The examples from the ontology demonstrate the disparities between what the image shows and what the human knows. The computer detects them autonomously and can initiate the appropriate post-processing. If and when this cognitive imaging prevails, the post-processed images may replace the current ones as legitimate artifact-free MRIs.

Virtual medicine: Utilization of the advanced cardiac imaging patient avatar for procedural planning and facilitation

Advances in imaging technology have led to a paradigm shift from planning of cardiovascular procedures and surgeries requiring the actual patient in a "brick and mortar" hospital to utilization of the digitalized patient in the virtual hospital. Cardiovascular computed tomographic angiography (CCTA) and cardiovascular magnetic resonance (CMR) digitalized 3-D patient representation of individual patient anatomy and physiology serves as an avatar allowing for virtual delineation of the most optimal approaches to cardiovascular procedures and surgeries prior to actual hospitalization. Pre-hospitalization reconstruction and analysis of anatomy and pathophysiology previously only accessible during the actual procedure could potentially limit the intrinsic risks related to time in the operating room, cardiac procedural laboratory and overall hospital environment. Although applications are specific to areas of cardiovascular specialty focus, there are unifying themes related to the utilization of technologies. The virtual patient avatar computer can also be used for procedural planning, computational modeling of anatomy, simulation of predicted therapeutic result, printing of 3-D models, and augmentation of real time procedural performance. Examples of the above techniques are at various stages of development for application to the spectrum of cardiovascular disease processes, including percutaneous, surgical and hybrid minimally invasive interventions. A multidisciplinary approach within medicine and engineering is necessary for creation of robust algorithms for maximal utilization of the virtual patient avatar in the digital medical center. Utilization of the virtual advanced cardiac imaging patient avatar will play an important role in the virtual health care system. Although there has been a rapid proliferation of early data, advanced imaging applications require further assessment and validation of accuracy, reproducibility, standardization, safety, efficacy, quality, cost effectiveness, and overall value to medical care.

Inverse problems in reduced order models of cardiovascular haemodynamics: aspects of data assimilation and heart rate variability

Inverse problems in cardiovascular modelling have become increasingly important to assess each patient individually. These problems entail estimation of patient-specific model parameters from uncertain measurements acquired in the clinic. In recent years, the method of data assimilation, especially the unscented Kalman filter, has gained popularity to address computational efficiency and uncertainty consideration in such problems. This work highlights and presents solutions to several challenges of this method pertinent to models of cardiovascular haemodynamics. These include methods to (i) avoid ill-conditioning of the covariance matrix, (ii) handle a variety of measurement types, (iii) include a variety of prior knowledge in the method, and (iv) incorporate measurements acquired at different heart rates, a common situation in the clinic where the patient state differs according to the clinical situation. Results are presented for two patient-specific cases of congenital heart disease. To illustrate and validate data assimilation with measurements at different heart rates, the results are presented on a synthetic dataset and on a patient-specific case with heart valve regurgitation. It is shown that the new method significantly improves the agreement between model predictions and measurements. The developed methods can be readily applied to other pathophysiologies and extended to dynamical systems which exhibit different responses under different sets of known parameters or different sets of inputs (such as forcing/excitation frequencies).

New imaging tools in cardiovascular medicine: computational fluid dynamics and 4D flow MRI

Blood flow imaging is a novel technology in cardiovascular medicine and surgery. Today, two types of blood flow imaging tools are available: measurement-based flow visualization including 4D flow MRI (or 3D cine phase-contrast magnetic resonance imaging), or echocardiography flow visualization software, and computer flow simulation modeling based on computational fluid dynamics (CFD). MRI and echocardiography flow visualization provide measured blood flow but have limitations in temporal and spatial resolution, whereas CFD flow calculates the flow according to assumptions instead of flow measurement, and it has sufficiently fine resolution up to the computer memory limit, and it enables even virtual surgery when combined with computer graphics. Blood flow imaging provides profound insight into the pathophysiology of cardiovascular diseases, because it quantifies and visualizes mechanical stress on the vessel walls or heart ventricle. Wall shear stress (WSS) is a stress on the endothelial wall caused by the near wall blood flow, and it is thought to be a predictor of atherosclerosis progression in coronary or aortic diseases. Flow energy loss (EL) is the loss of blood flow energy caused by viscous friction of turbulent diseased flow, and it is expected to be a predictor of ventricular workload on various heart diseases including heart valve disease, cardiomyopathy, and congenital heart diseases. Blood flow imaging can provide useful information for developing predictive medicine in cardiovascular diseases, and may lead to breakthroughs in cardiovascular surgery, especially in the decision-making process.

On the choice of outlet boundary conditions for patient-specific analysis of aortic flow using computational fluid dynamics

Boundary conditions (BCs) are an essential part in computational fluid dynamics (CFD) simulations of blood flow in large arteries. Although several studies have investigated the influence of BCs on predicted flow patterns and hemodynamic wall parameters in various arterial models, there is a lack of comprehensive assessment of outlet BCs for patient-specific analysis of aortic flow. In this study, five different sets of outlet BCs were tested and compared using a subject-specific model of a normal aorta. Phase-contrast magnetic resonance imaging (PC-MRI) was performed on the same subject and velocity profiles extracted from the in vivo measurements were used as the inlet boundary condition. Computational results obtained with different outlet BCs were assessed in terms of their agreement with the PC-MRI velocity data and key hemodynamic parameters, such as pressure and flow waveforms and wall shear stress related indices. Our results showed that the best overall performance was achieved by using a well-tuned three-element Windkessel model at all model outlets, which not only gave a good agreement with in vivo flow data, but also produced physiological pressure waveforms and values. On the other hand, opening outlet BCs with zero pressure at multiple outlets failed to reproduce any physiologically relevant flow and pressure features.

Role of Computational Modelling in Planning and Executing Interventional Procedures for Congenital Heart Disease

Increasingly, computational modelling and numerical simulations are used to help plan complex surgical and interventional cardiovascular procedures in children and young adults with congenital heart disease. From its origins more than 30 years ago, surgical planning with analysis of flow hemodynamics and energy loss/efficiency has helped design and implement many modifications to existing techniques. On the basis of patient-specific medical imaging, surgical planning allows accurate model production that can then be manipulated in a virtual surgical environment, with the proposed solutions finally tested with advanced computational fluid dynamics to evaluate the results. Applications include a broad range of congenital heart disease, including patients with single-ventricle anatomy undergoing staged palliation, those with arch obstruction, with double outlet right ventricle, or with tetralogy of Fallot. In the present work, we focus on clinical applications of this exciting field. We describe the framework for these techniques, including brief descriptions of the engineering principles applied and the interaction between "benchtop" data with medical decision-making. We highlight some early insights learned from pioneers over the past few decades, including refinements in Fontan baffle geometries and configurations. Finally, we offer a glimpse into exciting advances that are presently being explored, including use of modelling for transcatheter interventions. In this era of personalized medicine, computational modelling and surgical planning allows patient-specific tailoring of interventions to optimize clinical outcomes.

Tetralogy of Fallot Surgical Repair: Shunt Configurations, Ductus Arteriosus and the Circle of Willis

In this study, hemodynamic performance of three novel shunt configurations that are considered for the surgical repair of tetralogy of Fallot (TOF) disease are investigated in detail. Clinical experience suggests that the shunt location, connecting angle, and its diameter can influence the post-operative physiology and the neurodevelopment of the neonatal patient. An experimentally validated second order computational fluid dynamics (CFD) solver and a parametric neonatal diseased great artery model that incorporates the ductus arteriosus (DA) and the full patient-specific circle of Willis (CoW) are employed. Standard truncated resistance CFD boundary conditions are compared with the full cerebral arterial system, which resulted 21, -13, and 37% difference in flow rate at the brachiocephalic, left carotid, and subclavian arteries, respectively. Flow splits at the aortic arch and cerebral arteries are calculated and found to change with shunt configuration significantly for TOF disease. The central direct shunt (direct shunt) has pulmonary flow 5% higher than central oblique shunt (oblique shunt) and 23% higher than modified Blalock Taussig shunt (RPA shunt) while the DA is closed. Maximum wall shear stress (WSS) in the direct shunt configuration is 9 and 60% higher than that of the oblique and RPA shunts, respectively. Patent DA, significantly eliminated the pulmonary flow control function of the shunt repair. These results suggests that, due to the higher flow rates at the pulmonary arteries, the direct shunt, rather than the central oblique, or right pulmonary artery shunts could be preferred by the surgeon. This extended model introduced new hemodynamic performance indices for the cerebral circulation that can correlate with the post-operative neurodevelopment quality of the patient.

Effect of stenosis eccentricity on the functionality of coronary bifurcation lesions-a numerical study

Interventional cardiologists still rely heavily on angiography for the evaluation of coronary lesion severity, despite its poor correlation with the presence of ischemia. In order to improve the accuracy of the current diagnostic procedures, an understanding of the relative influence of geometric characteristics on the induction of ischemia is required. This idea is especially important for coronary bifurcation lesions (CBLs), whose treatment is complex and is associated with high rates of peri- and post-procedural clinical events. Overall, it is unclear which geometric and morphological parameters of CBLs influence the onset of ischemia. More specifically, the effect of stenosis eccentricity is unknown. Computational fluid dynamic simulations, under a geometric multiscale framework, were executed for seven CBL configurations within the left main coronary artery bifurcation. Both concentric and eccentric stenosis profiles of mild to severe constriction were considered. By using a geometric multiscale framework, the fractional flow reserve, which is the gold-standard clinical diagnostic index, could be calculated and was compared between the eccentric and concentric profiles for each case. The results suggested that for configurations where the supplying vessel is stenosed, eccentricity could have a notable effect on and therefore be an important factor that influences configuration functionality.

Adaptive outflow boundary conditions improve post-operative predictions after repair of peripheral pulmonary artery stenosis

Peripheral pulmonary artery stenosis (PPS) is a congenital abnormality resulting in pulmonary blood flow disparity and right ventricular hypertension. Despite recent advance in catheter-based interventions, surgical reconstruction is still preferred to treat complex PPS. However optimal surgical strategies remain unclear. It would be of great benefit to be able to predict post-operative hemodynamics to assist with surgical planning toward optimizing outcomes. While image-based computational fluid dynamics has been used in cardiovascular surgical planning, most studies have focused on the impact of local geometric changes on hemodynamic performance. Previous experimental studies suggest morphological changes in the pulmonary arteries not only alter local hemodynamics but also lead to distal pulmonary adaptation. In this proof of concept study, a constant shear stress hypothesis and structured pulmonary trees are used to derive adaptive outflow boundary conditions for post-operative simulations. Patient-specific simulations showed the adaptive outflow boundary conditions by the constant shear stress model to provide better predictions of pulmonary flow distribution than the conventional strategy of maintaining outflow boundary conditions. On average, the relative difference, when compared to the gold standard clinical test, in blood flow distribution to the right lung is reduced from 20 to 4 %. This suggests adaptive outflow boundary conditions should be incorporated into post-operative modeling in patients with complex PPS.

Data assimilation and modelling of patient-specific single-ventricle physiology with and without valve regurgitation

A closed-loop lumped parameter model of blood circulation is considered for single-ventricle shunt physiology. Its parameters are estimated by an inverse problem based on patient-specific haemodynamics measurements. As opposed to a black-box approach, maximizing the number of parameters that are related to physically measurable quantities motivates the present model. Heart chambers are described by a single-fibre mechanics model, and valve function is modelled with smooth opening and closure. A model for valve prolapse leading to valve regurgitation is proposed. The method of data assimilation, in particular the unscented Kalman filter, is used to estimate the model parameters from time-varying clinical measurements. This method takes into account both the uncertainty in prior knowledge related to the parameters and the uncertainty associated with the clinical measurements. Two patient-specific cases - one without regurgitation and one with atrioventricular valve regurgitation - are presented. Pulmonary and systemic circulation parameters are successfully estimated, without assumptions on their relationships. Parameters governing the behaviour of heart chambers and valves are either fixed based on biomechanics, or estimated. Results of the inverse problem are validated qualitatively through clinical measurements or clinical estimates that were not included in the parameter estimation procedure. The model and the estimation method are shown to successfully capture patient-specific clinical observations, even with regurgitation, such as the double peaked nature of valvular flows and anomalies in electrocardiogram readings. Lastly, biomechanical implications of the results are discussed.

Control of respiration-driven retrograde flow in the subdiaphragmatic venous return of the Fontan circulation

Respiration influences the subdiaphragmatic venous return in the total cavopulmonary connection (TCPC) of the Fontan circulation whereby both the inferior vena cava (IVC) and hepatic vein flows can experience retrograde motion. Controlling retrograde flows could improve patient outcomes. Using a patient-specific model within a Fontan mock circulatory system with respiration, we inserted a valve into the IVC to examine its effects on local hemodynamics while varying retrograde volumes by changing vascular impedances. A bovine valved conduit reduced IVC retrograde flow to within 3% of antegrade flow in all cases. The valve closed only under conditions supporting retrograde flow and its effects on local hemodynamics increased with larger retrograde volume. Liver and TCPC pressures improved only when the valve leaflets were closed whereas cycle-averaged pressures improved only slightly (<1 mm Hg). Increased pulmonary vascular resistance raised mean circulation pressures, but the valve functioned and cardiac output improved and stabilized. Power loss across the TCPC improved by 12%-15% (p < 0.05) with a valve. The effectiveness of valve therapy is dependent on patient vascular impedance.

Hemodynamic Impact of Superior Vena Cava Placement in the Y-Graft Fontan Connection

BACKGROUND

A Fontan Y-shaped graft using a commercially available aortoiliac graft has been used to connect the inferior vena cava (IVC) to the pulmonary arteries. This modification of the Fontan procedure seeks to improve hepatic flow distribution (HFD) to the lungs. However, patient-specific anatomical restrictions might limit the space available for graft placement. Altering the superior vena cava (SVC) positioning is hypothesized to provide more space for an optimal connection, avoiding caval flow collision. Computational modeling tools were used to retrospectively study the effect of SVC placement on Y-graft hemodynamics.

METHODS

Patient-specific anatomies (N = 10 patients) and vessel flows were reconstructed from retrospective cardiac magnetic resonance (CMR) images after Fontan Y-graft completion. Alternative geometries were created using a virtual surgery environment, altering the SVC position and the offset in relation to the Y-graft branches. Geometric characterization and computational fluid dynamics simulations were performed. Hemodynamic factors (power loss and HFD) were computed.

RESULTS

Patients with a higher IVC return showed less sensitivity to SVC positioning. Patients with low IVC flow showed varied HFD results, depending on SVC location. Balanced HFD values (50% to each lung) were obtained when the SVC lay completely between the Y-graft branches. The effect on power loss was patient specific.

CONCLUSIONS

SVC positioning with respect to the Y-graft affects HFD, especially in patients with lower IVC flow. Careful positioning of the SVC at the time of a bidirectional Glenn (BDG) procedure based on patient-specific anatomy can optimize the hemodynamics of the eventual Fontan completion.

Computational modeling and engineering in pediatric and congenital heart disease

PURPOSE OF REVIEW

Recent methodological advances in computational simulations are enabling increasingly realistic simulations of hemodynamics and physiology, driving increased clinical utility. We review recent developments in the use of computational simulations in pediatric and congenital heart disease, describe the clinical impact in modeling in single-ventricle patients, and provide an overview of emerging areas.

RECENT FINDINGS

Multiscale modeling combining patient-specific hemodynamics with reduced order (i.e., mathematically and computationally simplified) circulatory models has become the de-facto standard for modeling local hemodynamics and 'global' circulatory physiology. We review recent advances that have enabled faster solutions, discuss new methods (e.g., fluid structure interaction and uncertainty quantification), which lend realism both computationally and clinically to results, highlight novel computationally derived surgical methods for single-ventricle patients, and discuss areas in which modeling has begun to exert its influence including Kawasaki disease, fetal circulation, tetralogy of Fallot (and pulmonary tree), and circulatory support.

SUMMARY

Computational modeling is emerging as a crucial tool for clinical decision-making and evaluation of novel surgical methods and interventions in pediatric cardiology and beyond. Continued development of modeling methods, with an eye towards clinical needs, will enable clinical adoption in a wide range of pediatric and congenital heart diseases.

The effect of resolution on viscous dissipation measured with 4D flow MRI in patients with Fontan circulation: Evaluation using computational fluid dynamics

Viscous dissipation inside Fontan circulation, a parameter associated with the exercise intolerance of Fontan patients, can be derived from computational fluid dynamics (CFD) or 4D flow MRI velocities. However, the impact of spatial resolution and measurement noise on the estimation of viscous dissipation is unclear. Our aim was to evaluate the influence of these parameters on viscous dissipation calculation. Six Fontan patients underwent whole heart 4D flow MRI. Subject-specific CFD simulations were performed. The CFD velocities were down-sampled to isotropic spatial resolutions of 0.5mm, 1mm, 2mm and to MRI resolution. Viscous dissipation was compared between (1) high resolution CFD velocities, (2) CFD velocities down-sampled to MRI resolution, (3) down-sampled CFD velocities with MRI mimicked noise levels, and (4) in-vivo 4D flow MRI velocities. Relative viscous dissipation between subjects was also calculated. 4D flow MRI velocities (15.6 ± 3.8 cm/s) were higher, although not significantly different than CFD velocities (13.8 ± 4.7 cm/s, p=0.16), down-sampled CFD velocities (12.3 ± 4.4 cm/s, p=0.06) and the down-sampled CFD velocities with noise (13.2 ± 4.2 cm/s, p=0.06). CFD-based viscous dissipation (0.81 ± 0.55 mW) was significantly higher than those based on down-sampled CFD (0.25 ± 0.19 mW, p=0.03), down-sampled CFD with noise (0.49 ± 0.26 mW, p=0.03) and 4D flow MRI (0.56 ± 0.28 mW, p=0.06). Nevertheless, relative viscous dissipation between different subjects was maintained irrespective of resolution and noise, suggesting that comparison of viscous dissipation between patients is still possible.

Patient-Specific Surgical Planning, Where Do We Stand? The Example of the Fontan Procedure

The Fontan surgery for single ventricle heart defects is a typical example of a clinical intervention in which patient-specific computational modeling can improve patient outcome: with the functional heterogeneity of the presenting patients, which precludes generic solutions, and the clear influence of the surgically-created Fontan connection on hemodynamics, it is acknowledged that individualized computational optimization of the post-operative hemodynamics can be of clinical value. A large body of literature has thus emerged seeking to provide clinically relevant answers and innovative solutions, with an increasing emphasis on patient-specific approaches. In this review we discuss the benefits and challenges of patient-specific simulations for the Fontan surgery, reviewing state of the art solutions and avenues for future development. We first discuss the clinical impact of patient-specific simulations, notably how they have contributed to our understanding of the link between Fontan hemodynamics and patient outcome. This is followed by a survey of methodologies for capturing patient-specific hemodynamics, with an emphasis on the challenges of defining patient-specific boundary conditions and their extension for prediction of post-operative outcome. We conclude with insights into potential future directions, noting that one of the most pressing issues might be the validation of the predictive capabilities of the developed framework.

Measurement and modeling of coronary blood flow

Ischemic heart disease that comprises both coronary artery disease and microvascular disease is the single greatest cause of death globally. In this context, enhancing our understanding of the interaction of coronary structure and function is not only fundamental for advancing basic physiology but also crucial for identifying new targets for treating these diseases. A central challenge for understanding coronary blood flow is that coronary structure and function exhibit different behaviors across a range of spatial and temporal scales. While experimental studies have sought to understand this feature by isolating specific mechanisms, in tandem, computational modeling is increasingly also providing a unique framework to integrate mechanistic behaviors across different scales. In addition, clinical methods for assessing coronary disease severity are continuously being informed and updated by findings in basic physiology. Coupling these technologies, computational modeling of the coronary circulation is emerging as a bridge between the experimental and clinical domains, providing a framework to integrate imaging and measurements from multiple sources with mathematical descriptions of governing physical laws. State-of-the-art computational modeling is being used to combine mechanistic models with data to provide new insight into coronary physiology, optimization of medical technologies, and new applications to guide clinical practice.

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.

Relationship of single ventricle filling and preload to total cavopulmonary connection hemodynamics

BACKGROUND

Single ventricle lesions are associated with gradual attrition after surgical palliation with the total cavopulmonary connection (TCPC). Ventricular dysfunction is frequently noted, particularly impaired diastolic performance. This study seeks to relate TCPC hemodynamic energy losses to single ventricle volumes and filling characteristics.

METHODS

Cardiac magnetic resonance (CMR) data were retrospectively analyzed for 30 single ventricle patients at an average age of 12.7 ± 4.8 years. Cine ventricular short-axis scans were semiautomatically segmented for all cardiac phases. Ventricular volumes, ejection fraction, peak filling rate, peak ejection rate, and time to peak filling were calculated. Corresponding patient-specific TCPC geometry was acquired from a stack of transverse CMR images; relevant flow rates were segmented from through-plane phase contrast CMR data at TCPC inlets and outlets. The TCPC indexed power loss was calculated from computational fluid dynamics simulations using a validated custom solver. Time-averaged flow conditions and rigid vessel walls were assumed in all cases. Pearson correlations were used to detect relationships between variables, with p less than 0.05 considered significant.

RESULTS

Ventricular end-diastolic (R = -0.48) and stroke volumes (R = -0.37) had significant negative correlations with the natural logarithm of a flow-independent measure of power loss. This power loss measure also had a significant positive relationship to time to peak filling rate (normalized to cycle time; R = 0.67).

CONCLUSIONS

Flow-independent TCPC power loss is inversely related with ventricular end-diastolic and stroke volumes. Elevated power losses may contribute to impaired diastolic filling and limited preload reserve in single ventricle patients.

Hemodynamic effects of left pulmonary artery stenosis after superior cavopulmonary connection: a patient-specific multiscale modeling study

OBJECTIVE

Currently, no quantitative guidelines have been established for treatment of left pulmonary artery (LPA) stenosis. This study aims to quantify the effects of LPA stenosis on postoperative hemodynamics for single-ventricle patients undergoing stage II superior cavopulmonary connection (SCPC) surgery, using a multiscale computational approach.

METHODS

Image data from 6 patients were segmented to produce 3-dimensional models of the pulmonary arteries before stage II surgery. Pressure and flow measurements were used to tune a 0-dimensional model of the entire circulation. Postoperative geometries were generated through stage II virtual surgery; varying degrees of LPA stenosis were applied using mesh morphing and hemodynamics assessed through coupled 0-3-dimensional simulations. To relate metrics of stenosis to clinical classifications, pediatric cardiologists and surgeons ranked the degrees of stenosis in the models. The effects of LPA stenosis were assessed based on left-to-right pulmonary artery flow split ratios, mean pressure drop across the stenosis, cardiac pressure-volume loops, and other clinically relevant parameters.

RESULTS

Stenosis of >65% of the vessel diameter was required to produce a right pulmonary artery:LPA flow split <30%, and/or a mean pressure drop of >3.0 mm Hg, defined as clinically significant changes.

CONCLUSIONS

The effects of <65% stenosis on SCPC hemodynamics and physiology were minor and may not justify the increased complexity of adding LPA arterioplasty to the SCPC operation. However, in the longer term, pulmonary augmentation may affect outcomes of the Fontan completion surgery, as pulmonary artery distortion is a risk factor that may influence stage III physiology.

A methodological paradigm for patient-specific multi-scale CFD simulations: from clinical measurements to parameter estimates for individual analysis

A new framework for estimation of lumped (for instance, Windkessel) model parameters from uncertain clinical measurements is presented. The ultimate aim is to perform patient-specific haemodynamic analysis. This framework is based on sensitivity analysis tools and the sequential estimation approach of the unscented Kalman filter. Sensitivity analysis and parameter estimation are performed in lumped parameter models, which act as reduced order surrogates of the 3D domain for haemodynamic analysis. While the goal of sensitivity analysis is to assess potential identifiability problems, the unscented Kalman filter estimation leads to parameter estimates based on clinical measurements and modelling assumptions. An application of such analysis and parameter estimation methodology is demonstrated for synthetic and real data. Equality constraints on various physiological parameters are enforced. Since the accuracy of the Windkessel parameter estimates depends on the lumped parameter representativeness, the latter is iteratively improved by running few 3D simulations while simultaneously improving the former. Such a method is applied on a patient-specific aortic coarctation case. Less than 3% and 9% errors between the clinically measured quantities and 3D simulation results for rest and stress are obtained, respectively. Knowledge on how these Windkessel parameters change from rest to stress can thus be learned by such an approach. Lastly, it is demonstrated that the proposed approach is capable of dealing with a wide variety of measurements and cases where the pressure and flow clinical measurements are not taken simultaneously.

Surgical planning of the total cavopulmonary connection: robustness analysis

In surgical planning of the Fontan connection for single ventricle physiologies, there can be differences between the proposed and implemented options. Here, we developed a surgical planning framework that help determine the best performing option and ensures that the results will be comparable if there are slight geometrical variations. Eight patients with different underlying anatomies were evaluated in this study; surgical variations were created for each connection by changing either angle, offset or baffle diameter. Computational fluid dynamics were performed and the energy efficiency (indexed power loss-iPL) and hepatic flow distribution (HFD) computed. Differences with the original connection were evaluated: iPL was not considerably affected by the changes in geometry. For HFD, the single superior vena cava (SVC) connections presented less variability compared to the other anatomies. The Y-graft connection was the most robust overall, while the extra-cardiac connections showed dependency to offset. Bilateral SVC and interrupted inferior vena cava with azygous continuation showed high variability in HFD. We have developed a framework to assess the robustness of a surgical option for the TCPC; this will be useful to assess the most complex cases where pre-surgery planning could be most beneficial to ensure an efficient and robust hemodynamic performance.

Patient-specific assessment of cardiovascular function by combination of clinical data and computational model with applications to patients undergoing Fontan operation

The assessment of cardiovascular function is becoming increasingly important for the care of patients with single-ventricle defects. However, most measurement methods available in the clinical setting cannot provide a separate measure of cardiac function and loading conditions. In the present study, a numerical method has been proposed to compensate for the limitations of clinical measurements. The main idea was to estimate the parameters of a cardiovascular model by fitting model simulations to patient-specific clinical data via parameter optimization. Several strategies have been taken to establish a well-posed parameter optimization problem, including clinical data-matched model development, parameter selection based on an extensive sensitivity analysis, and proper choice of parameter optimization algorithm. The numerical experiments confirmed the ability of the proposed parameter optimization method to uniquely determine the model parameters given an arbitrary set of clinical data. The method was further tested in four patients undergoing the Fontan operation. Obtained results revealed a prevalence of ventricular abnormalities in the patient cohort and at the same time demonstrated the presence of marked inter-patient differences and preoperative to postoperative changes in cardiovascular function. Because the method allows a quick assessment and makes use of clinical data available in clinical practice, its clinical application is promising.

Computational fluid dynamics in Fontan patients to evaluate power loss during simulated exercise

OBJECTIVE

Exercise intolerance is common in total cavopulmonary connection (TCPC) patients. It has been suggested that power loss (Ploss) inside the TCPC plays a role in reduced exercise performance. Our objective is to establish the role of Ploss inside the TCPC during increased flow, simulating exercise in a patient-specific way.

METHODS

Cardiac MRI (CMR) was used to obtain flow rates from the caval veins during rest and increased flow, simulating exercise with dobutamine. A 3D reconstruction of the TCPC was created using CMR data. Computational fluid dynamics (CFD) simulations were performed to calculate Ploss inside the TCPC structure for rest and stress conditions. To reflect the flow distribution during exercise, a condition where inferior caval vein (IVC) flow was increased twofold compared with rest was added. 29 TCPC patients (15 intra-atrial lateral tunnel (ILT) and 14 extracardiac conduit (ECC)) were included.

RESULTS

Mean Ploss at rest was 1.36 ± 0.94 (ILT) and 3.20 ± 1.26 (ECC) mW/m(2) (p<0.001), 2.84 ±1.95 (ILT) and 8.41 ± 3.77 (ECC) mW/m(2) (p<0.001) during dobutamine and 5.21 ± 3.50 (ILT) and 15.28 ± 8.30 (ECC) mW/m(2) (p=0.001) with twofold IVC flow. The correlation between cardiac index and Ploss was exponential (ILT: R(2)=0.811, p<0.001; ECC: R(2)=0.690, p<0.001).

CONCLUSIONS

Ploss inside the TCPC structure is limited but increases with simulated exercise. This relates to the anatomy of TCPC and the surgical technique used. In all flow conditions, ILT patients have lower Ploss than ECC patients. We did not find a relationship between Ploss and exercise capacity.

Transient hemodynamic changes upon changing a BCPA into a TCPC in staged Fontan operation: a computational model study

The clinical benefits of the Fontan operation in treating single-ventricle defects have been well documented. However, perioperative mortality or morbidity remains a critical problem. The purpose of the present study was to identify the cardiovascular factors that dominate the transient hemodynamic changes upon the change of a bidirectional cavopulmonary (Glenn) anastomosis (BCPA) into a total cavopulmonary connection (TCPC). For this purpose, two computational models were constructed to represent, respectively, a single-ventricle circulation with a BCPA and that with a TCPC. A series of model-based simulations were carried out to quantify the perioperative hemodynamic changes under various cardiovascular conditions. Obtained results indicated that the presence of a low pulmonary vascular resistance and/or a low lower-body vascular resistance is beneficial to the increase in transpulmonary flow upon the BCPA to TCPC change. Moreover, it was found that ventricular diastolic dysfunction and mitral valve regurgitation, despite being well-known risk factors for poor postoperative outcomes, do not cause a considerable perioperative reduction in transpulmonary flow. The findings may help physicians to assess the perioperative risk of the TCPC surgery based on preoperative measurement of cardiovascular function.

Mock circulatory system of the Fontan circulation to study respiration effects on venous flow behavior

We describe an in vitro model of the Fontan circulation with respiration to study subdiaphragmatic venous flow behavior. The venous and arterial connections of a total cavopulmonary connection (TCPC) test section were coupled with a physical lumped parameter (LP) model of the circulation. Intrathoracic and subdiaphragmatic pressure changes associated with normal breathing were applied. This system was tuned for two patients (5 years, 0.67 m2; 10 years, 1.2 m2) to physiological values. System function was verified by comparison to the analytical model on which it was based and by consistency with published clinical measurements. Overall, subdiaphragmatic venous flow was influenced by respiration. Flow within the arteries and veins increased during inspiration but decreased during expiration, with retrograde flow in the inferior venous territories. System pressures and flows showed close agreement with the analytical LP model (p < 0.05). The ratio of the flow rates occurring during inspiration to expiration were within the clinical range of values reported elsewhere. The approach used to set up and control the model was effective and provided reasonable comparisons with clinical data.

Modeling single ventricle physiology: review of engineering tools to study first stage palliation of hypoplastic left heart syndrome

First stage palliation of hypoplastic left heart syndrome, i.e., the Norwood operation, results in a complex physiological arrangement, involving different shunting options (modified Blalock-Taussig, RV-PA conduit, central shunt from the ascending aorta) and enlargement of the hypoplastic ascending aorta. Engineering techniques, both computational and experimental, can aid in the understanding of the Norwood physiology and their correct implementation can potentially lead to refinement of the decision-making process, by means of patient-specific simulations. This paper presents some of the available tools that can corroborate clinical evidence by providing detailed insight into the fluid dynamics of the Norwood circulation as well as alternative surgical scenarios (i.e., virtual surgery). Patient-specific anatomies can be manufactured by means of rapid prototyping and such models can be inserted in experimental set-ups (mock circulatory loops) that can provide a valuable source of validation data as well as hydrodynamic information. Such models can be tuned to respond to differing the patient physiologies. Experimental set-ups can also be compatible with visualization techniques, like particle image velocimetry and cardiovascular magnetic resonance, further adding to the knowledge of the local fluid dynamics. Multi-scale computational models include detailed three-dimensional (3D) anatomical information coupled to a lumped parameter network representing the remainder of the circulation. These models output both overall hemodynamic parameters while also enabling to investigate the local fluid dynamics of the aortic arch or the shunt. As an alternative, pure lumped parameter models can also be employed to model Stage 1 palliation, taking advantage of a much lower computational cost, albeit missing the 3D anatomical component. Finally, analytical techniques, such as wave intensity analysis, can be employed to study the Norwood physiology, providing a mechanistic perspective on the ventriculo-arterial coupling for this specific surgical scenario.

Lagrangian analysis of hemodynamics data from FSI simulation

We present the computation of lagrangian-based flow characterization measures for time-dependent, deformable-wall, finite-element blood flow simulations. Applicability of the algorithm is demonstrated in a fluid-structure interaction simulation of blood flow through a total cavopulmonary connection (Fontan procedure), and results are compared with a rigid-vessel simulation. Specifically, we report on several important lagrangian-based measures including flow distributions, finite-time Lyapunov exponent fields, particle residence time, and exposure time calculations. Overall, strong similarity in lagrangian measures of the flow between deformable and rigid-vessel models was observed.