A mathematical model of circulation in the presence of the bidirectional cavopulmonary anastomosis in children with a univentricular heart

A multi-dimensional CFD framework for fast patient-specific fractional flow reserve prediction

Fractional flow reserve (FFR) is considered as the gold standard for diagnosing coronary myocardial ischemia. Existing 3D computational fluid dynamics (CFD) methods attempt to predict FFR noninvasively using coronary computed tomography angiography (CTA). However, the accuracy and efficiency of the 3D CFD methods in coronary arteries are considerably limited. In this work, we introduce a multi-dimensional CFD framework that improves the accuracy of FFR prediction by estimating 0D patient-specific boundary conditions, and increases the efficiency by generating 3D initial conditions. The multi-dimensional CFD models contain the 3D vascular model for coronary simulation, the 1D vascular model for iterative optimization, and the 0D vascular model for boundary conditions expression. To improve the accuracy, we utilize clinical parameters to derive 0D patient-specific boundary conditions with an optimization algorithm. To improve the efficiency, we evaluate the convergence state using the 1D vascular model and obtain the convergence parameters to generate appropriate 3D initial conditions. The 0D patient-specific boundary conditions and the 3D initial conditions are used to predict FFR (FFRC). We conducted a retrospective study involving 40 patients (61 diseased vessels) with invasive FFR and their corresponding CTA images. The results demonstrate that the FFRC and the invasive FFR have a strong linear correlation (r = 0.80, p < 0.001) and high consistency (mean difference: 0.014 ±0.071). After applying the cut-off value of FFR (0.8), the accuracy, sensitivity, specificity, positive predictive value, and negative predictive value of FFRC were 88.5%, 93.3%, 83.9%, 84.8%, and 92.9%, respectively. Compared with the conventional zero initial conditions method, our method improves prediction efficiency by 71.3% per case. Therefore, our multi-dimensional CFD framework is capable of improving the accuracy and efficiency of FFR prediction significantly.

Recasting Current Knowledge of Human Fetal Circulation: The Importance of Computational Models

Computational hemodynamic simulations are becoming increasingly important for cardiovascular research and clinical practice, yet incorporating numerical simulations of human fetal circulation is relatively underutilized and underdeveloped. The fetus possesses unique vascular shunts to appropriately distribute oxygen and nutrients acquired from the placenta, adding complexity and adaptability to blood flow patterns within the fetal vascular network. Perturbations to fetal circulation compromise fetal growth and trigger the abnormal cardiovascular remodeling that underlies congenital heart defects. Computational modeling can be used to elucidate complex blood flow patterns in the fetal circulatory system for normal versus abnormal development. We present an overview of fetal cardiovascular physiology and its evolution from being investigated with invasive experiments and primitive imaging techniques to advanced imaging (4D MRI and ultrasound) and computational modeling. We introduce the theoretical backgrounds of both lumped-parameter networks and three-dimensional computational fluid dynamic simulations of the cardiovascular system. We subsequently summarize existing modeling studies of human fetal circulation along with their limitations and challenges. Finally, we highlight opportunities for improved fetal circulation models.

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.

Computational fluid dynamics investigation of the novel hybrid comprehensive stage II operation

Background

The hybrid comprehensive stage 2 (HCS2) procedure is a novel palliative operation applicable to a select subset of single ventricle patients with adequate native antegrade aortic flow to the upper body. Flow to the descending aorta, through the pulmonary outlet and ductal arch, is influenced by a stented intrapulmonary baffle connecting the branch pulmonary arteries. We used computational fluid dynamics (CFD) to elucidate the hemodynamic characteristics of this reconstruction.

Methods

We used multiscale CFD analysis of a synthetic, patient-derived HCS2 anatomic configuration with unsteady laminar flow conditions and a non-Newtonian blood model to quantify the resultant hemodynamics. The 3-dimensional CFD model was coupled to a 0-dimensional lumped parameter model of the peripheral circulation to determine the required boundary conditions.

Results

For the specific anatomy studied, the intrapulmonary baffle did not obstruct flow from the pulmonary trunk to ductal arch as long as the distance between the anterior pulmonary artery wall and baffle wall exceeded ∼7 mm. Vortex shedding off of the baffle wall did not develop, because of the short distance to the ductal arch. The stented baffle experienced significantly uneven "inward" loading from the systemic side. Pulmonary outlet flow separation distal to the baffle produced a low-speed recirculation region.

Conclusions

Hemodynamic patterns in this complex anatomy are generally favorable. Low flow recirculation could be mitigated by preoperative shape optimization. Calculated inward stresses on the pulmonary baffle can be used in the future to study baffle stent deformation, which is expected to be small.

Predictive Modeling of Secondary Pulmonary Hypertension in Left Ventricular Diastolic Dysfunction

Diastolic dysfunction is a common pathology occurring in about one third of patients affected by heart failure. This condition may not be associated with a marked decrease in cardiac output or systemic pressure and therefore is more difficult to diagnose than its systolic counterpart. Compromised relaxation or increased stiffness of the left ventricle induces an increase in the upstream pulmonary pressures, and is classified as secondary or group II pulmonary hypertension (2018 Nice classification). This may result in an increase in the right ventricular afterload leading to right ventricular failure. Elevated pulmonary pressures are therefore an important clinical indicator of diastolic heart failure (sometimes referred to as heart failure with preserved ejection fraction, HFpEF), showing significant correlation with associated mortality. However, accurate measurements of this quantity are typically obtained through invasive catheterization and after the onset of symptoms. In this study, we use the hemodynamic consistency of a differential-algebraic circulation model to predict pulmonary pressures in adult patients from other, possibly non-invasive, clinical data. We investigate several aspects of the problem, including the ability of model outputs to represent a sufficiently wide pathologic spectrum, the identifiability of the model's parameters, and the accuracy of the predicted pulmonary pressures. We also find that a classifier using the assimilated model parameters as features is free from the problem of missing data and is able to detect pulmonary hypertension with sufficiently high accuracy. For a cohort of 82 patients suffering from various degrees of heart failure severity, we show that systolic, diastolic, and wedge pulmonary pressures can be estimated on average within 8, 6, and 6 mmHg, respectively. We also show that, in general, increased data availability leads to improved predictions.

Implementation and Calibration of a Deep Neural Network to Predict Parameters of Left Ventricular Systolic Function Based on Pulmonary and Systemic Arterial Pressure Signals

The evaluation of cardiac contractility by the assessment of the ventricular systolic elastance function is clinically challenging and cannot be easily obtained at the bedside. In this work, we present a framework characterizing left ventricular systolic function from clinically readily available data, including systemic and pulmonary arterial pressure signals. We implemented and calibrated a deep neural network (DNN) consisting of a multi-layer perceptron with 4 fully connected hidden layers and with 16 neurons per layer, which was trained with data obtained from a lumped model of the cardiovascular system modeling different levels of cardiac function. The lumped model included a function of circulatory autoregulation from carotid baroreceptors in pulsatile conditions. Inputs for the DNN were systemic and pulmonary arterial pressure curves. Outputs from the DNN were parameters of the lumped model characterizing left ventricular systolic function, especially end-systolic elastance. The DNN adequately performed and accurately recovered the relevant hemodynamic parameters with a mean relative error of less than 2%. Therefore, our framework can easily provide complex physiological parameters of cardiac contractility, which could lead to the development of invaluable tools for the clinical evaluation of patients with severe cardiac dysfunction.

A mathematical model of human heart including the effects of heart contractility varying with heart rate changes

Incorporating the intrinsic variability of heart contractility varying with heart rate into the mathematical model of human heart would be useful for addressing the dynamical behaviors of human cardiovascular system, but models with such features were rarely reported. This study focused on the development and evaluation of a mathematical model of the whole heart, including the effects of heart contractility varying with heart rate changes. This model was developed based on a paradigm and model presented by Ottesen and Densielsen, which was used to model ventricular contraction. A piece-wise function together with expressions for time-related parameters were constructed for modeling atrial contraction. Atrial and ventricular parts of the whole heart model were evaluated by comparing with models from literature, and then the whole heart model were assessed through coupling with a simple model of the systemic circulation system and the pulmonary circulation system. The results indicated that both atrial and ventricular parts of the whole heart model could reasonably reflect their contractility varying with heart rate changes, and the whole heart model could exhibit major features of human heart. Results of the parameters variation studies revealed the correlations between the parameters in the whole heart model and performances (including the maximum pressure and the stroke volume) of every chamber. These results would be useful for helping users to adjust parameters in special applications.

Lumped parameter model for hemodynamic simulation of congenital heart diseases

The recent development of computer technology has made it possible to simulate the hemodynamics of congenital heart diseases on a desktop computer. However, multi-scale modeling of the cardiovascular system based on computed tomographic and magnetic resonance images still requires long simulation times. The lumped parameter model is potentially beneficial for real-time bedside simulation of congenital heart diseases. In this review, we introduce the basics of the lumped parameter model (time-varying elastance chamber model combined with modified Windkessel vasculature model) and illustrate its usage in hemodynamic simulation of congenital heart diseases using examples such as hypoplastic left heart syndrome and Fontan circulation. We also discuss the advantages of the lumped parameter model and the problems for clinical use.

Numerical Simulation of a Realistic Total Cavo-pulmonary Connection: Effect of Unbalanced Pulmonary Resistances on Hydrodynamic Performance

Total cavo pulmonary connection (TCPC) is one of the surgical techniques adopted to compensate the failure of the right heart in pediatric patients. The main goal of this procedure is the realization of a configuration for the caval veins and for the pulmonary arteries that can guarantee as low as possible pressure losses and appropriate lung perfusion. Starting from this point of view, a realistic TCPC with extracardiac conduit (TECPC) is investigated by means of Computational Fluid Dynamics (CFD) to evaluate the pressure loss under different pressure conditions, simulating different vessel resistances, on the pulmonary arteries. A total flow of 3 L/min, with a distribution between the inferior vena cava (IVC) and the superior vena cava (SVC) equal to 6/4, was investigated; three different boundary conditions for the pressure were imposed, resulting in three simulations in steady-state conditions, to the right pulmonary artery (RPA) and to the left pulmonary artery (LPA), simulating a balanced (deltaP(LPA-RPA) = 0 mmHg) and two unbalanced pulmonary resistances to blood flow (a pressure difference deltaP(LPA-RPA) = +/- 2 mmHg, respectively). The geometry for the TECPC was realized according to MRI derived physiological values for the vessels and for the configuration adopted for the anastomosis (the extra-cardiac conduit was inclined 22 degrees towards the left pulmonary artery with respect to the IVC axis). The computed power losses agree with previous in vitro Particle Image Velocimetry investigations. The results show that a higher resistance on the LPA causes the greater pressure loss for the TECPC under study, while the minimum pressure loss can be achieved balancing the pulmonary resistances, subsequently obtaining a balanced flow repartition towards the lungs.

Physical and Computer Modelling of Blood Flow in a Systemic-to-pulmonary Shunt

The aim of this work was the application of computer and physical in vitro simulation methods for estimating surgery procedure hemodynamics. The modified Blalock-Taussig (mB-T) palliative surgical procedure is performed to increase the pulmonary blood flow in children with congenital heart defects. Such a systemic-to-pulmonary shunt yields substantial modification in the blood flow within the large blood vessels. The objective of the present study was to investigate basic characteristics of the flow, flow pattern and pressure-flow efficiency, before and after opening of the mB-T graft.

Methods

The model was based on the vessel geometry obtained from the Visible Human Project and included the arch of aorta, the three arteries branching from the arch, the pulmonary trunck, and the left and right pulmonary arteries. The graft was added between the left subclavian artery and the left pulmonary artery. The glass model of the vessels was produced and investigated in a physical model of the cardiovascular system with an artificial ventricular device as the blood pump. Flow rate and hydrostatic pressure were measured at the inlet to and outlets from the glass model and in a few points within the system. Laser flow visualization was also performed. Computer simulations were done using the boundary conditions from the physical model.

Results

The opening of the mB-T graft changed flow distribution in all branches (including inflow). A complex flow pattern with large eddies and channelling of the flow in the vicinity of the graft and within it was observed in flow visualization and in computer simulations. Because of that complexity the local measurements of hydrostatic pressure at the vessel wall could not predict the average flow rate. The reversed flow in the graft was observed during the systole.

Conclusions

The complex flow pattern developed in the physical model of the mB-T graft. The channelling of the flow and the formation of large eddies may yield high shear stress and modify blood properties. The rigid wall model can describe only some flow characteristics observed in vivo. Computer simulation is a very fast and accurate method which permits earlier qualification of cardiac surgeons on how to change cardiac vascular blood flow after operations.

Computer and Physical Modeling of Blood Circulation Pump Support for a New Field of Application in Palliative Surgery

Objectives

One of the most popular palliative procedures performed to increase pulmonary blood flow in children with congenital heart defects is a shunt operation (Blalock-Taussig graft or Glenn procedure), which creates the new blood channel to the pulmonary artery. The main problem with this kind of surgery is poor shunt effectiveness and the lack of possibility to regulate the flow. The aim of this work is to use advanced computer simulation methods to study the effectiveness of a new idea to introduce a small axial blood pump into a Blalock-Taussig (B-T) or Glenn shunt in order to control the blood flow and prevent any increase in the graft stenosis.

Methods

Physical and computer 3-D simulation based on a finite element mesh (FEM) model was applied. Studies for optimization of the shunt and hybrid shunt with pump were performed for different stages of the disease.

Results and Conclusion

The graft with the axial pump creates good conditions for the vascular system and pulmonary artery blood flow as well as regulating blood pressure under variable conditions caused by palliative procedures. Its use permits the afterload of the left heart ventricle to be decreased. A palliative procedure is only a temporary solution. When a child grows, while the graft size is fixed, the blood flow through this graft may be not sufficient under changing hemodynamic conditions. The use of an axial pump for regulating the blood flow volume, during palliative procedures, allows to obtain the optimal flow conditions in pulmonary artery and safely wait on the final cardiac surgery correction later. However, the use of a pump mounted inside the graft increased hemodynamic resistance, which caused the flow to decrease up to 70% in the graft when the axial pump was not working.

Patient-specific parameter estimation in single-ventricle lumped circulation models under uncertainty

Computational models of cardiovascular physiology can inform clinical decision-making, providing a physically consistent framework to assess vascular pressures and flow distributions, and aiding in treatment planning. In particular, lumped parameter network (LPN) models that make an analogy to electrical circuits offer a fast and surprisingly realistic method to reproduce the circulatory physiology. The complexity of LPN models can vary significantly to account, for example, for cardiac and valve function, respiration, autoregulation, and time-dependent hemodynamics. More complex models provide insight into detailed physiological mechanisms, but their utility is maximized if one can quickly identify patient specific parameters. The clinical utility of LPN models with many parameters will be greatly enhanced by automated parameter identification, particularly if parameter tuning can match non-invasively obtained clinical data. We present a framework for automated tuning of 0D lumped model parameters to match clinical data. We demonstrate the utility of this framework through application to single ventricle pediatric patients with Norwood physiology. Through a combination of local identifiability, Bayesian estimation and maximum a posteriori simplex optimization, we show the ability to automatically determine physiologically consistent point estimates of the parameters and to quantify uncertainty induced by errors and assumptions in the collected clinical data. We show that multi-level estimation, that is, updating the parameter prior information through sub-model analysis, can lead to a significant reduction in the parameter marginal posterior variance. We first consider virtual patient conditions, with clinical targets generated through model solutions, and second application to a cohort of four single-ventricle patients with Norwood physiology. Copyright © 2016 John Wiley & Sons, Ltd.

Mathematical Analysis and Physical Profile of Blalock-Taussig Shunt and Sano Modification Procedure in Hypoplastic Left Heart Syndrome: Review of the Literature and Implications for the Anesthesiologist

The first stage of surgical treatment for hypoplastic left heart syndrome (HLHS) includes the creation of artificial systemic-to-pulmonary connections to provide pulmonary blood flow. The modified Blalock-Taussig (mBT) shunt has been the technique of choice for this procedure; however, a right ventricle-pulmonary artery (RV-PA) shunt has been introduced into clinical practice with encouraging but still conflicting outcomes when compared with the mBT shunt. The aim of this study is to explore mathematical modeling as a tool for describing physical profiles that could assist the surgical team in predicting complications related to stenosis and malfunction of grafts in an attempt to find correlations with clinical outcomes from clinical studies that compared both surgical techniques and to assist the anesthesiologist in making decisions to manage patients with this complex cardiac anatomy. Mathematical modeling to display the physical characteristics of the chosen surgical shunt is a valuable tool to predict flow patterns, shear stress, and rate distribution as well as energetic performance at the graft level and relative to ventricular efficiency. Such predictions will enable the surgical team to refine the technique so that hemodynamic complications be anticipated and prevented, and are also important for perioperative management by the anesthesia team.

Partial hepatectomy hemodynamics changes: Experimental data explained by closed-loop lumped modeling

The liver function may be degraded after partial liver ablation surgery. Adverse liver hemodynamics have been shown to be associated to liver failure. The link between these hemodynamics changes and ablation size is however poorly understood. This article proposes to explain with a closed-loop lumped model the hemodynamics changes observed during twelve surgeries in pigs. The portal venous tree is modeled with a pressure-dependent variable resistor. The variables measured, before liver ablation, are used to tune the model parameters. Then, the liver partial ablation is simulated with the model and the simulated pressures and flows are compared with post-operative measurements. Fluid infusion and blood losses occur during the surgery. The closed-loop model presented accounts for these blood volume changes. Moreover, the impact of blood volume changes and the liver lobe mass estimations on the simulated variables is studied. The typical increase of portal pressure, increase of liver pressure loss, slight decrease of portal flow and major decrease in arterial flow are quantitatively captured by the model for a 75% hepatectomy. It appears that the 75% decrease in hepatic arterial flow can be explained by the resistance increase induced by the surgery, and that no hepatic arterial buffer response (HABR) mechanism is needed to account for this change. The different post-operative states, observed in experiments, are reproduced with the proposed model. Thus, an explanation for inter-subjects post-operative variability is proposed. The presented framework can easily be adapted to other species circulations and to different pathologies for clinical hepatic applications.

Numerical investigation of the haemodynamics in the human fetal umbilical vein/ductus venosus based on the experimental data

Abortion of the fetus due to a disease, in an early stage of pregnancy, has been dramatically increased in the last decades. There is a still lack of knowledge on the various types of diseases which lead fetus to a vulnerable circumstance. The transport of oxygenated blood from the placenta to the human fetus has been an important clinical feature in Doppler velocimetry studies, especially the ductus venosus (DV). The DV connects intra-abdominal portion of the umbilical vein and the inferior vena cava (IVC) at the inlet of the right atrium and is, therefore, important when examining the fetus state of health. An abnormal flow in the DV can indicate a fetal disease such as, chromosomal abnormalities, cardiac defect, hypoxaemia and intrauterine growth restriction (IUGR). The blood flow in the fetal circulation has not been investigated much in detail. The blood flow in the fetal circulation provides necessary information for physician to make a suitable decision on abortion or alternative medical practice before or even after birth. The present study performed a comparative study to quantify the blood velocity in DV by a combination approach based on 3D computational simulation and Doppler measurement. The results showed that the velocity value in DV is significant and can be considered as an indicator of any kind of disease in fetal. The nodal displacement of the model was also analysed. It shows that DV tolerates a higher level of displacement compared with the other regions of the model, whereas the nodal pressure shows different results as the lowest values are located in DV.

Hemodynamic Effects of Ventricular Assist Device Implantation on Norwood, Glenn, and Fontan Circulation: A Simulation Study

The growing population of failing single-ventricle (SV) patients might benefit from ventricular assist device (VAD) support as a bridge to heart transplantation. However, the documented experience is limited to isolated case reports. Considering the complex and different physiopathology of Norwood, Glenn, and Fontan patients and the lack of established experience, the aim of this work is to realize and test a lumped parameter model of the cardiovascular system able to simulate SV hemodynamics and VAD implantation effects to support clinical decision. Hemodynamic and echocardiographic data of 30 SV patients (10 Norwood, 10 Glenn, and 10 Fontan) were retrospectively collected and used to simulate patients' baseline. Then, the effects of VAD implantation were simulated. Simulation results suggest that the implantation of VAD: (i) increases the cardiac output and the mean arterial systemic pressure in all the three palliation conditions (Norwood 77.2 and 19.7%, Glenn 38.6 and 32.2%, and Fontan 17.2 and 14.2%); (ii) decreases the SV external work (Norwood 55%, Glenn 35.6%, and Fontan 41%); (iii) decreases the pressure pulsatility index (Norwood 65.2%, Glenn 81.3%, and Fontan 64.8%); (iv) increases the pulmonary arterial pressure in particular in the Norwood circulation (Norwood 39.7%, Glenn 12.1% and Fontan 3%); and (v) decreases the atrial pressure (Norwood 2%, Glenn 10.6%, and Fontan 8.6%). Finally, the VAD work is lower in the Norwood circulation (30.4 mL·mm Hg) in comparison with Fontan (40.3 mL·mm Hg) and to Glenn (64.5 mL·mm Hg) circulations. The use of VAD in SV physiology could be helpful to bridge patients to heart transplantations by increasing the CO and unloading the SV with a decrement of the atrial pressure and the SV external work. The regulation of the pulmonary flow is challenging because the Pap is increased by the presence of VAD. The hemodynamic changes are different in the different SV palliation step. The use of numerical models could be helpful to support patient and VAD selection to optimize the clinical outcome.

Exercise capacity in the Bidirectional Glenn physiology: Coupling cardiac index, ventricular function and oxygen extraction ratio

In Bi-directional Glenn (BDG) physiology, the superior systemic circulation and pulmonary circulation are in series. Consequently, only blood from the superior vena cava is oxygenated in the lungs. Oxygenated blood then travels to the ventricle where it is mixed with blood returning from the lower body. Therefore, incremental changes in oxygen extraction ratio (OER) could compromise exercise tolerance. In this study, the effect of exercise on the hemodynamic and ventricular performance of BDG physiology was investigated using clinical patient data as inputs for a lumped parameter model coupled with oxygenation equations. Changes in cardiac index, Qp/Qs, systemic pressure, oxygen extraction ratio and ventricular/vascular coupling ratio were calculated for three different exercise levels. The patient cohort (n=29) was sub-grouped by age and pulmonary vascular resistance (PVR) at rest. It was observed that the changes in exercise tolerance are significant in both comparisons, but most significant when sub-grouped by PVR at rest. Results showed that patients over 2 years old with high PVR are above or close to the upper tolerable limit of OER (0.32) at baseline. Patients with high PVR at rest had very poor exercise tolerance while patients with low PVR at rest could tolerate low exercise conditions. In general, ventricular function of SV patients is too poor to increase CI and fulfill exercise requirements. The presented mathematical model provides a framework to estimate the hemodynamic performance of BDG patients at different exercise levels according to patient specific data.

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.

Hemodynamic performance of the Fontan circulation compared with a normal biventricular circulation: a computational model study

The physiological limitations of the Fontan circulation have been extensively addressed in the literature. Many studies emphasized the importance of pulmonary vascular resistance in determining cardiac output (CO) but gave little attention to other cardiovascular properties that may play considerable roles as well. The present study was aimed to systemically investigate the effects of various cardiovascular properties on clinically relevant hemodynamic variables (e.g., CO and central venous pressure). To this aim, a computational modeling method was employed. The constructed models provided a useful tool for quantifying the hemodynamic effects of any cardiovascular property of interest by varying the corresponding model parameters in model-based simulations. Herein, the Fontan circulation was studied compared with a normal biventricular circulation so as to highlight the unique characteristics of the Fontan circulation. Based on a series of numerical experiments, it was found that 1) pulmonary vascular resistance, ventricular diastolic function, and systemic vascular compliance play a major role, while heart rate, ventricular contractility, and systemic vascular resistance play a secondary role in the regulation of CO in the Fontan circulation; 2) CO is nonlinearly related to any single cardiovascular property, with their relationship being simultaneously influenced by other cardiovascular properties; and 3) the stability of central venous pressure is significantly reduced in the Fontan circulation. The findings suggest that the hemodynamic performance of the Fontan circulation is codetermined by various cardiovascular properties and hence a full understanding of patient-specific cardiovascular conditions is necessary to optimize the treatment of Fontan patients.

Modeling biology spanning different scales: an open challenge

It is coming nowadays more clear that in order to obtain a unified description of the different mechanisms governing the behavior and causality relations among the various parts of a living system, the development of comprehensive computational and mathematical models at different space and time scales is required. This is one of the most formidable challenges of modern biology characterized by the availability of huge amount of high throughput measurements. In this paper we draw attention to the importance of multiscale modeling in the framework of studies of biological systems in general and of the immune system in particular.

Effect of high altitude exposure on the hemodynamics of the bidirectional Glenn physiology: modeling incremented pulmonary vascular resistance and heart rate

The considerable blood mixing in the bidirectional Glenn (BDG) physiology further limits the capacity of the single working ventricle to pump enough oxygenated blood to the circulatory system. This condition is exacerbated under severe conditions such as physical activity or high altitude. In this study, the effect of high altitude exposure on hemodynamics and ventricular function of the BDG physiology is investigated. For this purpose, a mathematical approach based on a lumped parameter model was developed to model the BDG circulation. Catheterization data from 39 BDG patients at stabilized oxygen conditions was used to determine baseline flows and pressures for the model. The effect of high altitude exposure was modeled by increasing the pulmonary vascular resistance (PVR) and heart rate (HR) in increments up to 80% and 40%, respectively. The resulting differences in vascular flows, pressures and ventricular function parameters were analyzed. By simultaneously increasing PVR and HR, significant changes (p <0.05) were observed in cardiac index (11% increase at an 80% PVR and 40% HR increase) and pulmonary flow (26% decrease at an 80% PVR and 40% HR increase). Significant increase in mean systemic pressure (9%) was observed at 80% PVR (40% HR) increase. The results show that the poor ventricular function fails to overcome the increased preload and implied low oxygenation in BDG patients at higher altitudes, especially for those with high baseline PVRs. The presented mathematical model provides a framework to estimate the hemodynamic performance of BDG patients at different PVR increments.

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.

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

The objective of this work is to perform a virtual planning of surgical repairs in patients with congenital heart diseases--to test the predictive capability of a closed-loop multi-scale model. As a first step, we reproduced the pre-operative state of a specific patient with a univentricular circulation and a bidirectional cavopulmonary anastomosis (BCPA), starting from the patient's clinical data. Namely, by adopting a closed-loop multi-scale approach, the boundary conditions at the inlet and outlet sections of the three-dimensional model were automatically calculated by a lumped parameter network. Successively, we simulated three alternative surgical designs of the total cavopulmonary connection (TCPC). In particular, a T-junction of the venae cavae to the pulmonary arteries (T-TCPC), a design with an offset between the venae cavae (O-TCPC) and a Y-graft design (Y-TCPC) were compared. A multi-scale closed-loop model consisting of a lumped parameter network representing the whole circulation and a patient-specific three-dimensional finite volume model of the BCPA with detailed pulmonary anatomy was built. The three TCPC alternatives were investigated in terms of energetics and haemodynamics. Effects of exercise were also investigated. Results showed that the pre-operative caval flows should not be used as boundary conditions in post-operative simulations owing to changes in the flow waveforms post-operatively. The multi-scale approach is a possible solution to overcome this incongruence. Power losses of the Y-TCPC were lower than all other TCPC models both at rest and under exercise conditions and it distributed the inferior vena cava flow evenly to both lungs. Further work is needed to correlate results from these simulations with clinical outcomes.

Use of mathematical modeling to compare and predict hemodynamic effects between hybrid and surgical Norwood palliations for hypoplastic left heart syndrome

BACKGROUND

Combining bilateral pulmonary artery banding with arterial duct stenting, the hybrid approach achieves stage 1 palliation for hypoplastic left heart syndrome with different flow characteristics than those after the surgical Norwood procedures. Accordingly, we used computational modeling to assess some of these differences, including influence on systemic and cerebral oxygen deliveries.

METHODS AND RESULTS

A 3-dimensional computational model of hybrid palliation was developed by the finite volume method, along with models of the Norwood operation with a modified Blalock-Tausig or right ventricle-to-pulmonary artery shunt. Hybrid circulation was modeled with a 7-mm ductal stent and bilateral pulmonary artery banding to a 2-mm diameter. A 3.5-mm conduit was used in the Blalock-Tausig shunt model, whereas a 5-mm conduit was used in the right ventricle-to-pulmonary artery shunt model. Coupled to all the models was an identical hydraulic network that described the entire circulatory system based on pre-stage 2 hemodynamics. This clinically validated multiscale approach predicts flow dynamics, as well as global cardiac output, mixed venous oxygen saturation, and systemic and cerebral oxygen delivery. Compared with either of the Norwood models, the hybrid palliation had higher pulmonary-to-systemic flow ratio and lower cardiac output. Total systemic oxygen delivery was markedly reduced in the hybrid palliation (Blalock-Tausig shunt 591, right ventricle-to-pulmonary artery shunt 640, and hybrid 475 mL · min(-1) · m(-2)). Cerebral oxygen delivery was similarly lower in the hybrid palliation.

CONCLUSIONS

These computational results suggest that the hybrid approach may provide inferior systemic and cerebral oxygen deliveries compared with either of the 2 surgical Norwood procedures before stage 2 palliation.

Review of zero-D and 1-D models of blood flow in the cardiovascular system

Background

Zero-dimensional (lumped parameter) and one dimensional models, based on simplified representations of the components of the cardiovascular system, can contribute strongly to our understanding of circulatory physiology. Zero-D models provide a concise way to evaluate the haemodynamic interactions among the cardiovascular organs, whilst one-D (distributed parameter) models add the facility to represent efficiently the effects of pulse wave transmission in the arterial network at greatly reduced computational expense compared to higher dimensional computational fluid dynamics studies. There is extensive literature on both types of models.

Method and Results

The purpose of this review article is to summarise published 0D and 1D models of the cardiovascular system, to explore their limitations and range of application, and to provide an indication of the physiological phenomena that can be included in these representations. The review on 0D models collects together in one place a description of the range of models that have been used to describe the various characteristics of cardiovascular response, together with the factors that influence it. Such models generally feature the major components of the system, such as the heart, the heart valves and the vasculature. The models are categorised in terms of the features of the system that they are able to represent, their complexity and range of application: representations of effects including pressure-dependent vessel properties, interaction between the heart chambers, neuro-regulation and auto-regulation are explored. The examination on 1D models covers various methods for the assembly, discretisation and solution of the governing equations, in conjunction with a report of the definition and treatment of boundary conditions. Increasingly, 0D and 1D models are used in multi-scale models, in which their primary role is to provide boundary conditions for sophisticate, and often patient-specific, 2D and 3D models, and this application is also addressed. As an example of 0D cardiovascular modelling, a small selection of simple models have been represented in the CellML mark-up language and uploaded to the CellML model repository http://models.cellml.org/. They are freely available to the research and education communities.

Conclusion

Each published cardiovascular model has merit for particular applications. This review categorises 0D and 1D models, highlights their advantages and disadvantages, and thus provides guidance on the selection of models to assist various cardiovascular modelling studies. It also identifies directions for further development, as well as current challenges in the wider use of these models including service to represent boundary conditions for local 3D models and translation to clinical application.

Boundary conditions of patient-specific fluid dynamics modelling of cavopulmonary connections: possible adaptation of pulmonary resistances results in a critical issue for a virtual surgical planning

Cavopulmonary connections are surgical procedures used to treat a variety of complex congenital cardiac defects. Virtual pre-operative planning based on in silico patient-specific modelling might become a powerful tool in the surgical decision-making process. For this purpose, three-dimensional models can be easily developed from medical imaging data to investigate individual haemodynamics. However, the definition of patient-specific boundary conditions is still a crucial issue. The present study describes an approach to evaluate the vascular impedance of the right and left lungs on the basis of pre-operative clinical data and numerical simulations. Computational fluid dynamics techniques are applied to a patient with a bidirectional cavopulmonary anastomosis, who later underwent a total cavopulmonary connection (TCPC). Multi-scale models describing the surgical region and the lungs are adopted, while the flow rates measured in the venae cavae are used at the model inlets. Pre-operative and post-operative conditions are investigated; namely, TCPC haemodynamics, which are predicted using patient-specific pre-operative boundary conditions, indicates that the pre-operative balanced lung resistances are not compatible with the TCPC measured flows, suggesting that the pulmonary vascular impedances changed individually after the surgery. These modifications might be the consequence of adaptation to the altered pulmonary blood flows.

Pulmonary regurgitation: the effects of varying pulmonary artery compliance, and of increased resistance proximal or distal to the compliance

BACKGROUND

Pulmonary regurgitation is common after repair of tetralogy of Fallot, predisposing to right ventricular dilatation and potentially fatal arrhythmias. Magnetic resonance studies of such patients led us to hypothesize that the amount of regurgitation, in the absence of an effective valve, depends on pulmonary arterial compliance and on the location of resistance relative to the compliance.

METHODS AND RESULTS

Using a pre-existing mathematical model representing the cardiovascular system, removal of the virtual pulmonary valve gave a triphasic pulmonary artery flow curve similar in shape to those recorded in patients with free regurgitation, with a regurgitant fraction of 30%. There was no reversal of flow at pulmonary capillary level, the regurgitant volume originating entirely from the compliance of the virtual pulmonary arteries and arterioles. Doubling their compliance increased regurgitation to 35%, whereas halving it decreased regurgitation to 23%. Doubling the total pulmonary vascular resistance by increasing arteriolar resistance increased regurgitation to 46%, whereas doubling it by simulating pulmonary annular stenosis proximal to the compliance limited regurgitation to 10%, but at the cost of a 32 mmHg peak systolic pressure drop.

CONCLUSIONS

The model supported our hypotheses, indicating the relevance to pulmonary regurgitation of previously overlooked variables. The virtual pulmonary regurgitation was exacerbated by pulmonary artery compliance and by elevated resistance distal to it, but was limited by more proximal resistance. These relationships merit careful clinical investigation as they would have implications for the initial management, subsequent investigation and decisions on re-intervention in patients with pulmonary regurgitation.

Comparison of Particle Image Velocimetry and Phase Contrast MRI in a Patient-Specific Extracardiac Total Cavopulmonary Connection

Particle image velocimetry (PIV) and phase contrast magnetic resonance imaging (PC-MRI) have not been compared in complex biofluid environments. Such analysis is particularly useful to investigate flow structures in the correction of single ventricle congenital heart defects, where fluid dynamic efficiency is essential. A stereolithographic replica of an extracardiac total cavopulmonary connection (TCPC) is studied using PIV and PC-MRI in a steady flow loop. Volumetric two-component PIV is compared to volumetric three-component PC-MRI at various flow conditions. Similar flow structures are observed in both PIV and PC-MRI, where smooth flow dominates the extracardiac TCPC, and superior vena cava flow is preferential to the right pulmonary artery, while inferior vena cava flow is preferential to the left pulmonary artery. Where three-component velocity is available in PC-MRI studies, some helical flow in the extracardiac TCPC is observed. Vessel cross sections provide an effective means of validation for both experiments, and velocity magnitudes are of the same order. The results highlight similarities to validate flow in a complex patient-specific extracardiac TCPC. Additional information obtained by velocity in three components further describes the complexity of the flow in anatomic structures.

Modeling the Fontan circulation: where we are and where we need to go

The Fontan procedure and its subsequent modifications over the past 30 years can be described as a class of surgical procedures for patients born with complex congenital heart disease exhibiting a single-ventricle physiology. The long-term outcome for children currently undergoing a Fontan procedure remains worrisome because of multiple late morbidities observed. Despite significant modeling efforts spanning three decades, improvements to the Fontan procedure have occurred without comprehensive validation from these modeling studies. Careful examination shows that modeling studies to date offer only a "glimpse through a keyhole" into understanding and modeling a representative range of the variations in anatomy and physiology that exist in Fontan patients. Suggestions for future investigations are provided.

Simulations of congenital septal defect closure and reactivity testing in patient-specific models of the pediatric pulmonary vasculature: A 3D numerical study with fluid-structure interaction

Clinical imaging methods are highly effective in the diagnosis of vascular pathologies, but they do not currently provide enough detail to shed light on the cause or progression of such diseases, and would be hard pressed to foresee the outcome of surgical interventions. Greater detail of and prediction capabilities for vascular hemodynamics and arterial mechanics are obtained here through the coupling of clinical imaging methods with computational techniques. Three-dimensional, patient-specific geometric reconstructions of the pediatric proximal pulmonary vasculature were obtained from x-ray angiogram images and meshed for use with commercial computational software. Two such models from hypertensive patients, one with multiple septal defects, the other who underwent vascular reactivity testing, were each completed with two sets of suitable fluid and structural initial and boundary conditions and used to obtain detailed transient simulations of artery wall motion and hemodynamics in both clinically measured and predicted configurations. The simulation of septal defect closure, in which input flow and proximal vascular stiffness were decreased, exhibited substantial decreases in proximal velocity, wall shear stress (WSS), and pressure in the post-op state. The simulation of vascular reactivity, in which distal vascular resistance and proximal vascular stiffness were decreased, displayed negligible changes in velocity and WSS but a significant drop in proximal pressure in the reactive state. This new patient-specific technique provides much greater detail regarding the function of the pulmonary circuit than can be obtained with current medical imaging methods alone, and holds promise for enabling surgical planning.

Effects of atrial contraction, atrioventricular interaction and heart valve dynamics on human cardiovascular system response

Various simulation models of different complexity have been proposed to model the dynamic response of the human cardiovascular system. In a related paper we proposed an improved numerical model to study the dynamic response of the cardiovascular system, and the pressures, volumes and flow-rates in the four chambers of the heart, which included the effects of atrial contraction, atrioventricular interaction, and heart valve dynamics. This paper investigates the effects of each one of these aspects of the model on the overall dynamic system response. The dynamic response is studied under different situations, with and without including the effect of various features of the model, and these situations are studied and compared among themselves and to detailed aspects of expected healthy-system response. As an important contribution with potential clinical applications, this paper examines the corresponding effects of atrioventricular interaction, and heart valve opening and closing dynamics to the general system dynamic response. This isolation of physical cause-effect relationships is difficult to study with purely experimental methods. The simulation results agree well with results in the open literature. Comparison shows that introduction of these new features greatly improves the simulation accuracy of the effects of a, v and c waves, and in predicting regurgitant valve flow, the dichrotic notch, and E/A velocity ratio.

A concentrated parameter model for the human cardiovascular system including heart valve dynamics and atrioventricular interaction

Numerical modeling of the human cardiovascular system has always been an active research direction since the 19th century. In the past, various simulation models of different complexities were proposed for different research purposes. In this paper, an improved numerical model to study the dynamic function of the human circulation system is proposed. In the development of the mathematical model, the heart chambers are described with a variable elastance model. The systemic and pulmonary loops are described based on the resistance-compliance-inertia concept by considering local effects of flow friction, elasticity of blood vessels and inertia of blood in different segments of the blood vessels. As an advancement from previous models, heart valve dynamics and atrioventricular interaction, including atrial contraction and motion of the annulus fibrosus, are specifically modeled. With these improvements the developed model can predict several important features that were missing in previous numerical models, including regurgitant flow on heart valve closure, the value of E/A velocity ratio in mitral flow, the motion of the annulus fibrosus (called the KG diaphragm pumping action), etc. These features have important clinical meaning and their changes are often related to cardiovascular diseases. Successful simulation of these features enhances the accuracy of simulations of cardiovascular dynamics, and helps in clinical studies of cardiac function.

Multiscale modeling of the cardiovascular system: application to the study of pulmonary and coronary perfusions in the univentricular circulation

The objective of this study is to compare the coronary and pulmonary blood flow dynamics resulting from two configurations of systemic-to-pulmonary artery shunts currently utilized during the Norwood procedure: the central (CS) and modified Blalock Taussig (MBTS) shunts. A lumped parameter model of the neonatal cardiovascular circulation and detailed 3-D models of the shunt based on the finite volume method were constructed. Shunt sizes of 3, 3.5 and 4 mm were considered. A multiscale approach was adopted to prescribe appropriate and realistic boundary conditions for the 3-D models of the Norwood circulation. Results showed that the average shunt flow rate is higher for the CS option than for the MBTS and that pulmonary flow increases with shunt size for both options. Cardiac output is higher for the CS option for all shunt sizes. Flow distribution between the left and the right pulmonary arteries is not completely balanced, although for the CS option the discrepancy is low (50-51% of the pulmonary flow to the right lung) while for the MBTS it is more pronounced with larger shunt sizes (51-54% to the left lung). The CS option favors perfusion to the right lung while the MBTS favors the left. In the CS option, a smaller percentage of aortic flow is distributed to the coronary circulation, while that percentage rises for the MBTS. These findings may have important implications for coronary blood flow and ventricular function.

Adaptation to mechanical load determines shape and properties of heart and circulation: the CircAdapt model

With circulatory pathology, patient-specific simulation of hemodynamics is required to minimize invasiveness for diagnosis, treatment planning, and followup. We investigated the advantages of a smart combination of often already known hemodynamic principles. The CircAdapt model was designed to simulate beat-to-beat dynamics of the four-chamber heart with systemic and pulmonary circulation while incorporating a realistic relation between pressure-volume load and tissue mechanics and adaptation of tissues to mechanical load. Adaptation was modeled by rules, where a locally sensed signal results in a local action of the tissue. The applied rules were as follows: For blood vessel walls, 1) flow shear stress dilates the wall and 2) tensile stress thickens the wall; for myocardial tissue, 3) strain dilates the wall material, 4) larger maximum sarcomere length increases contractility, and 5) contractility increases wall mass. The circulation was composed of active and passive compliances and inertias. A realistic circulation developed by self-structuring through adaptation provided mean levels of systemic pressure and flow. Ability to simulate a wide variety of patient-specific circumstances was demonstrated by application of the same adaptation rules to the conditions of fetal circulation followed by a switch to the newborn circulation around birth. It was concluded that a few adaptation rules, directed to normalize mechanical load of the tissue, were sufficient to develop and maintain a realistic circulation automatically. Adaptation rules appear to be the key to reduce dramatically the number of input parameters for simulating circulation dynamics. The model may be used to simulate circulation pathology and to predict effects of treatment.

Theoretical analysis of rest and exercise hemodynamics in patients with total cavopulmonary connection

The objective of this study was to determine the impact of a total cavopulmonary connection on the main hemodynamic quantities, both at rest and during exercise, when compared with normal biventricular circulation. The analysis was performed by means of a mathematical model of the cardiovascular system. The model incorporates the main parameters of systemic and pulmonary circulation, the pulsating heart, and the action of arterial and cardiopulmonary baroreflex mechanisms. Furthermore, the effect of changes in intrathoracic pressure on venous return is also incorporated. Finally, the response to moderate dynamic exercise is simulated, including the effect of a central command, local metabolic vasodilation, and the "muscle pump" mechanism. Simulations of resting conditions indicate that the action of baroreflex regulatory mechanisms alone can only partially compensate for the absence of the right heart. Cardiac output and mean systemic arterial pressure at rest show a large decrease compared with the normal subject. More acceptable hemodynamic quantity values are obtained by combining the action of regulatory mechanisms with a chronic change in parameters affecting mean filling pressure. With such changes assumed, simulations of the response to moderate exercise show that univentricular circulation exhibits a poor capacity to increase cardiac output and to sustain aerobic metabolism, especially when the oxygen consumption rate is increased above 1.2-1.3 l/min. The model ascribes the poor response to exercise in these patients to the incapacity to sustain venous return caused by the high resistance to venous return and/or to exhaustion of volume compensation reserve.

Modeling of the Norwood circulation: effects of shunt size, vascular resistances, and heart rate

Hypoplastic left heart syndrome is the most common lethal cardiac malformation of the newborn. Its treatment, apart from heart transplantation, is the Norwood operation. The initial procedure for this staged repair consists of reconstructing a circulation where a single outlet from the heart provides systemic perfusion and an interpositioning shunt contributes blood flow to the lungs. To better understand this unique physiology, a computational model of the Norwood circulation was constructed on the basis of compartmental analysis. Influences of shunt diameter, systemic and pulmonary vascular resistance, and heart rate on the cardiovascular dynamics and oxygenation were studied. Simulations showed that 1) larger shunts diverted an increased proportion of cardiac output to the lungs, away from systemic perfusion, resulting in poorer O2 delivery, 2) systemic vascular resistance exerted more effect on hemodynamics than pulmonary vascular resistance, 3) systemic arterial oxygenation was minimally influenced by heart rate changes, 4) there was a better correlation between venous O2 saturation and O2 delivery than between arterial O2 saturation and O2 delivery, and 5) a pulmonary-to-systemic blood flow ratio of 1 resulted in optimal O2 delivery in all physiological states and shunt sizes.

Analysis by mathematical model of haemodynamic data in the failing Fontan circulation

Several late complications jeopardize the clinical performance of recipients of the Fontan operation. The underlying causes have been referred to disturbed flow dynamics in the cavopulmonary connections. Presumably, the large pressure drops occurring in the inferior and superior connections play a pivotal role in the pressure level of the entire circulation, especially in the venous. To address this issue, we retrospectively reviewed catheterization data of six patients with failing Fontan circulation and compared them with those of six patients with functioning Fontan circulation. The impact on the systemic and pulmonary pressure of the increase in the cavopulmonary connection resistances was studied through a steady-state mathematical model of the univentricular closed-loop circulation. In the patients with failing Fontan, pressure in the venae cavae was found to be significantly higher, especially at the inferior cava (19.3 +/- 2.2 versus 12.5 +/- 2.3 mmHg) with the pressure drop at the inferior cavopulmonary connection significantly increased (4.7 +/- 3.1 versus 0.33 +/- 0.82 mmHg). The proposed mathematical model permits us to clearly relate the pressure increase in the venae cavae to an increased resistance in the cavopulmonary connections. Therefore, the present analysis confirms that, to avoid possible congestion of venous circulation, the definitive palliation of univentricular heart should not cause pressure drops at the cavopulmonary connections.

Use of mathematical model to predict hemodynamics in cavopulmonary anastomosis with persistent forward flow

BACKGROUND

The bidirectional cavopulmonary anastomosis with additional pulmonary blood flow is used as a staged procedure or a definitive palliation of univentricular hearts. In this paper the flow competition occurring between the caval and the pulmonary flows is investigated. The hemodynamics in the superior vena cava and the blood flow distribution into the lungs, as well as the systemic arterial oxygen availability, are correlated with the severity of the right ventricle outflow tract obstruction and the pulmonary arteriolar resistance.

MATERIALS AND METHODS

Computer models of the pre- and postoperative hemodynamics of univentricular hearts were developed. The effects of increasing severity of the right ventricle outflow tract obstruction, with a pulmonary arteriolar resistance ranging from 0.8 to 7.9 nonindexed Woods units, were simulated.

RESULTS

The study indicates that the presence of an additional pulmonary blood flow from the native pulmonary artery may be beneficial. Since an excessive additional blood flow may cause central venous hypertension, its optimal value should be chosen according to the value of pulmonary arteriolar resistance. The model was utilized to simulate four clinical cases.

CONCLUSIONS

The simulations show that the model can predict the postoperative hemodynamics and could therefore be usefully applied to predict quantitatively the effect of the native pulmonary blood flow following bidirectional cavopulmonary anastomosis.