The total cavopulmonary connection resistance: a significant impact on single ventricle hemodynamics at rest and exercise

Effects of cardiac growth on electrical dyssynchrony in the single ventricle patient

Single ventricle patients, including those with hypoplastic left heart syndrome (HLHS), typically undergo three palliative heart surgeries culminating in the Fontan procedure. HLHS is associated with high rates of morbidity and mortality, and many patients develop arrhythmias, electrical dyssynchrony, and eventually ventricular failure. However, the correlation between ventricular enlargement and electrical dysfunction in HLHS physiology remains poorly understood. Here we characterize the relationship between growth and electrophysiology in HLHS using computational modeling. We integrate a personalized finite element model, a volumetric growth model, and a personalized electrophysiology model to perform controlled in silico experiments. We show that right ventricle enlargement negatively affects QRS duration and interventricular dyssynchrony. Conversely, left ventricle enlargement can partially compensate for this dyssynchrony. These findings have potential implications on our understanding of the origins of electrical dyssynchrony and, ultimately, the treatment of HLHS patients.

Predicting the effects of surgically determined parameters on exercise tolerance in Fontan patients

The Fontan physiology is a surgically created circulation for patients with a single functioning ventricle. Patients with this circulation tend to have lower exercise tolerance compared to those with a normal circulation. Important computational and experimental work has been done to investigate this reduction in exercise tolerance. However, there are few systematic modeling studies that focus on the effect of several surgically determined parameters within the same framework. We propose a mathematical model to describe the Fontan circulation under exercise. We then formulate a heuristic based on clinical data from Fontan patients to estimate exercise tolerance. The model is used to investigate the effect of three important surgically determined parameters on exercise tolerance: the systemic arterial compliance, the systemic-venous to pulmonary-venous fenestration, and the resistance of the total cavopulmonary connection.

Evaluation of the total hydrodynamic energy loss using 4D flow MRI in a case with Fontan failure

Fontan Failure (FF) is a common problem for single-ventricle patients as they reach adulthood. Although several mechanisms may cause FF, an optimized blood flow stream through the surgical conduits is essential to avoid excessive energy loss (EL). Recent clinical studies showed EL is related to the quality of life, exercise capacity, and hepatic function since the single-ventricle feeds pulmonary and systemic circulation serially. 4D flow MRI effectively estimates EL in Fontan circulation and allows clinicians to compare the effectiveness of the treatment strategy concerning pre-intervention. Here, we present 26-year-old women with FF who had normal cardiac catheterization findings and were treated according to high EL definitions that are measured through 4D flow MRI.

In vitro hemodynamic performance of a blood pump for self-powered venous assist in univentricular hearts

Univentricular heart anomalies represent a group of severe congenital heart defects necessitating early surgical intervention in infancy. The Fontan procedure, the final stage of single-ventricle palliation, establishes a serial connection between systemic and pulmonary circulation by channeling venous return to the lungs. The absence of the subpulmonary ventricle in this peculiar circulation progressively eventuates in failure, primarily due to chronic elevation in inferior vena cava (IVC) pressure. This study experimentally validates the effectiveness of an intracorporeally-powered venous ejector pump (VEP) in reducing IVC pressure in Fontan patients. The VEP exploits a fraction of aortic flow to create a jet-venturi effect for the IVC, negating the external power requirement and driveline infections. An invitro Fontan mock-up circulation loop is developed and the impact of VEP design parameters and physiological conditions is assessed using both idealized and patient-specific total cavopulmonary connection (TCPC) phantoms. The VEP performance in reducing IVC pressure exhibited an inverse relationship with the cardiac output and extra-cardiac conduit (ECC) size and a proportional relationship with the transpulmonary pressure gradient (TPG) and mean arterial pressure (MAP). The ideal VEP with fail-safe features provided an IVC pressure drop of 1.82 ± 0.49, 2.45 ± 0.54, and 3.12 ± 0.43 mm Hg for TPG values of 6, 8, and 10 mm Hg, respectively, averaged over all ECC sizes and cardiac outputs. Furthermore, the arterial oxygen saturation was consistently maintained above 85% during full-assist mode. These results emphasize the potential utility of the VEP to mitigate elevated venous pressure in Fontan patients.

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.

Fontan Heart: Insight Into the Physiological Role of the Right Heart

BACKGROUND

Cardiac output (CO) is almost normal in children born without a functional right ventricle (RV), and a Fontan repair, so why is RV dysfunction such a clinical problem? We tested the hypotheses that increased pulmonary vascular resistance (PVR) is the dominant factor and volume expansion by any means is of limited benefit.

METHODS

We removed the RV from a previously used MATLAB model and altered vascular volume, venous compliance (Cv), PVR, and measures of left ventricular (LV) systolic and diastolic function. CO and regional vascular pressures were the primary outcome measures.

RESULTS

RV removal decreased CO by 25%, and raised mean systemic filling pressure (MSFP). A 10 mL/kg increase in stressed volume only moderately increased CO with or without the RV. Decreasing systemic Cv increased CO but also markedly increased pulmonary venous pressure. With no RV, increasing PVR had the greatest effect on CO. Increasing LV function had little benefit.

CONCLUSIONS

Model data indicate that increasing PVR dominates the decrease in CO in Fontan physiology. Increasing stressed volume by any means only moderately increased CO and increasing LV function had little effect. Decreasing systemic Cv unexpectedly markedly increased pulmonary venous pressures even with the RV intact.

Beyond CFD: Emerging methodologies for predictive simulation in cardiovascular health and disease

Physics-based computational models of the cardiovascular system are increasingly used to simulate hemodynamics, tissue mechanics, and physiology in evolving healthy and diseased states. While predictive models using computational fluid dynamics (CFD) originated primarily for use in surgical planning, their application now extends well beyond this purpose. In this review, we describe an increasingly wide range of modeling applications aimed at uncovering fundamental mechanisms of disease progression and development, performing model-guided design, and generating testable hypotheses to drive targeted experiments. Increasingly, models are incorporating multiple physical processes spanning a wide range of time and length scales in the heart and vasculature. With these expanded capabilities, clinical adoption of patient-specific modeling in congenital and acquired cardiovascular disease is also increasing, impacting clinical care and treatment decisions in complex congenital heart disease, coronary artery disease, vascular surgery, pulmonary artery disease, and medical device design. In support of these efforts, we discuss recent advances in modeling methodology, which are most impactful when driven by clinical needs. We describe pivotal recent developments in image processing, fluid-structure interaction, modeling under uncertainty, and reduced order modeling to enable simulations in clinically relevant timeframes. In all these areas, we argue that traditional CFD alone is insufficient to tackle increasingly complex clinical and biological problems across scales and systems. Rather, CFD should be coupled with appropriate multiscale biological, physical, and physiological models needed to produce comprehensive, impactful models of mechanobiological systems and complex clinical scenarios. With this perspective, we finally outline open problems and future challenges in the field.

Haemodynamic performance of 16-20-mm extracardiac Goretex conduits in adolescent Fontan patients at rest and during simulated exercise

OBJECTIVES

To date, it is not known if 16-20-mm extracardiac conduits are outgrown during somatic growth from childhood to adolescence. This study aims to determine total cavopulmonary connection (TCPC) haemodynamics in adolescent Fontan patients at rest and during simulated exercise and to assess the relationship between conduit size and haemodynamics.

METHODS

Patient-specific, magnetic resonance imaging-based computational fluid dynamic models of the TCPC were performed in 51 extracardiac Fontan patients with 16-20-mm conduits. Power loss, pressure gradient and normalized resistance were quantified in rest and during simulated exercise. The cross-sectional area (CSA) (mean and minimum) of the vessels of the TCPC was determined and normalized for flow rate (mm2/l/min). Peak (predicted) oxygen uptake was assessed.

RESULTS

The median age was 16.2 years (Q1-Q3 14.0-18.2). The normalized mean conduit CSA was 35-73% smaller compared to the inferior and superior vena cava, hepatic veins and left/right pulmonary artery (all P < 0.001). The median TCPC pressure gradient was 0.7 mmHg (Q1-Q3 0.5-0.8) and 2.0 (Q1-Q3 1.4-2.6) during rest and simulated exercise, respectively. A moderate-strong inverse non-linear relationship was present between normalized mean conduit CSA and TCPC haemodynamics in rest and exercise. TCPC pressure gradients of ≥1.0 at rest and ≥3.0 mmHg during simulated exercise were observed in patients with a conduit CSA ≤ 45 mm2/l/min and favourable haemodynamics (<1 mmHg during both rest and exercise) in conduits ≥125 mm2/l/min. Normalized TCPC resistance correlated with (predicted) peak oxygen uptake.

CONCLUSIONS

Extracardiac conduits of 16-20 mm have become relatively undersized in most adolescent Fontan patients leading to suboptimal haemodynamics.

4D flow cardiovascular magnetic resonance derived energetics in the Fontan circulation correlate with exercise capacity and CMR-derived liver fibrosis/congestion

AIM

This study explores the relationship between in vivo 4D flow cardiovascular magnetic resonance (CMR) derived blood flow energetics in the total cavopulmonary connection (TCPC), exercise capacity and CMR-derived liver fibrosis/congestion.

BACKGROUND

The Fontan circulation, in which both caval veins are directly connected with the pulmonary arteries (i.e. the TCPC) is the palliative approach for single ventricle patients. Blood flow efficiency in the TCPC has been associated with exercise capacity and liver fibrosis using computational fluid dynamic modelling. 4D flow CMR allows for assessment of in vivo blood flow energetics, including kinetic energy (KE) and viscous energy loss rate (EL).

METHODS

Fontan patients were prospectively evaluated between 2018 and 2021 using a comprehensive cardiovascular and liver CMR protocol, including 4D flow imaging of the TCPC. Peak oxygen consumption (VO2) was determined using cardiopulmonary exercise testing (CPET). Iron-corrected whole liver T1 (cT1) mapping was performed as a marker of liver fibrosis/congestion. KE and EL in the TCPC were computed from 4D flow CMR and normalized for inflow. Furthermore, blood flow energetics were compared between standardized segments of the TCPC.

RESULTS

Sixty-two Fontan patients were included (53% male, 17.3 ± 5.1 years). Maximal effort CPET was obtained in 50 patients (peak VO2 27.1 ± 6.2 ml/kg/min, 56 ± 12% of predicted). Both KE and EL in the entire TCPC (n = 28) were significantly correlated with cT1 (r = 0.50, p = 0.006 and r = 0.39, p = 0.04, respectively), peak VO2 (r = - 0.61, p = 0.003 and r = - 0.54, p = 0.009, respectively) and % predicted peak VO2 (r = - 0.44, p = 0.04 and r = - 0.46, p = 0.03, respectively). Segmental analysis indicated that the most adverse flow energetics were found in the Fontan tunnel and left pulmonary artery.

CONCLUSIONS

Adverse 4D flow CMR derived KE and EL in the TCPC correlate with decreased exercise capacity and increased levels of liver fibrosis/congestion. 4D flow CMR is promising as a non-invasive screening tool for identification of patients with adverse TCPC flow efficiency.

Impairments in Pulmonary Function in Fontan Patients: Their Causes and Consequences

Patients with a Fontan circulation lack a sub-pulmonary ventricle with pulmonary blood flow passively redirected to the lungs. In the Fontan circulation, ventilation has a significant influence on pulmonary blood flow and cardiac output both at rest and with exercise. Children and adults with a Fontan circulation have abnormalities in lung function. In particular, restrictive ventilatory patterns, as measured by spirometry, and impaired gas transfer, as measured by the diffusing capacity of carbon monoxide, have been frequently observed. These abnormalities in lung function are associated with reduced exercise capacity and quality of life. Moderate to severe impairment in lung volumes is independently associated with reduced survival in adults with congenital heart disease. Skeletal and inspiratory muscle weakness has also been reported in patients with a Fontan circulation, with the prospect of improving respiratory muscle function through exercise training programs. In this review, we will present data on cardiopulmonary interactions in the Fontan circulation, the prevalence and severity of impaired lung function, and respiratory muscle function in this population. We will discuss potential causes for and consequence of respiratory impairments, and their impact on exercise capacity and longer-term Fontan outcome. We aim to shed light on possible strategies to reduce morbidity by improving respiratory function in this growing population of patients.

Segmental assessment of blood flow efficiency in the total cavopulmonary connection using four-dimensional flow magnetic resonance imaging: vortical flow is associated with increased viscous energy loss rate

Aims

To study flow-related energetics in multiple anatomical segments of the total cavopulmonary connection (TCPC) in Fontan patients from four-dimensional (4D) flow magnetic resonance imaging (MRI), and to study the relationship between adverse flow patterns and segment-specific energetics.

Methods and results

Twenty-six extracardiac Fontan patients underwent 4D flow MRI of the TCPC. A segmentation of the TCPC was automatically divided into five anatomical segments [conduit, superior vena cava (SVC), right/left pulmonary artery (LPA), and the Fontan confluence]. The presence of vortical flow in the pulmonary arteries or Fontan confluence was qualitatively scored. Kinetic energy (KE), viscous energy loss rate, and vorticity were calculated from the 4D flow MRI velocity field and normalized for segment length and/or inflow. Energetics were compared between segments and the relationship between vortical flow and segment cross-sectional area (CSA) with segment-specific energetics was determined. Vortical flow in the LPA (n = 6) and Fontan confluence (n = 12) were associated with significantly higher vorticity (P = 0.001 and P = 0.015, respectively) and viscous energy loss rate (P = 0.046 and P = 0.04, respectively) compared to patients without vortical flow. The LPA and conduit segments showed the highest KE and viscous energy loss rate, while most favourable energetics were observed in the SVC. Conduit CSA inversely correlated with KE (r = −0.614, P = 0.019) and viscous energy loss rate (r = −0.652, P = 0.011).

Conclusions

Vortical flow in the Fontan confluence and LPA associated with significantly increased viscous energy loss rate. Four-dimensional flow MRI-derived energetics may be used as a screening tool for direct, MRI-based assessment of flow efficiency in the TCPC.

Total vascular resistance increases during volume-unloading in asymptomatic single ventricle patients

OBJECTIVE

While the surgical stages of single ventricle (SV) palliation serve to separate pulmonary venous and systemic venous return, and to volume-unload the SV, staged palliation also results in transition from parallel to series circulation, increasing total vascular resistance. How this transition affects pressure loading of the SV is as yet unreported.

METHODS

We performed a retrospective chart review of Stage I, II, and III cardiac catheterization (CC) and echocardiographic data from 2001-2017 in all SV pts, with focus on systemic, pulmonary, and total vascular resistance (SVR, PVR, TVR respectively). Longitudinal analyses were performed with log-transformed variables. Effects of SVR-lowering medications were analyzed using Wilcoxon rank-sum testing.

RESULTS

There were 372 total patients who underwent CC at a Stage I (median age of 4.4 months, n=310), Stage II (median age 2.7 years, n = 244), and Stage III (median age 7.3 years, n = 113). Total volume loading decreases with progression to Stage III (P< 0.001). While PVR gradually increases from Stage II to Stage III, and SVR increases from Stage I to Stage III, TVR dramatically increases with progress towards series circulation. TVR was not affected by use of systemic vasodilator therapy. TVR, PVR, SVR, and CI did not correlate with indices of SV function at Stage III.

CONCLUSIONS

TVR steadily increases with an increasing contribution from SVR over progressive stages. TVR was not affected by systemic vasodilator agents. TVR did not correlate with echo-based indices of SV function. Further studies are needed to see if modulating TVR can improve exercise tolerance and outcomes.

Ventricular myocardial response to exercise in patients with Fontan circulation

BACKGROUND

Exercise stress echocardiography has been used to assess myocardial reserve in various heart diseases. This study examined the ventricular myocardial response to exercise in Fontan patients using exercise stress echocardiography.

METHODS

Twenty-five Fontan patients and 19 control subjects underwent semi-supine bicycle exercise stress echocardiography in this prospective, single-center, cross-sectional study. Pulsed-wave Doppler tissue imaging peak systolic (s') and diastolic (e') velocities, longitudinal strain and systolic strain rate, and early diastolic strain rate data at rest and at peak exercise were obtained for the systemic ventricle. The myocardial reserve of functional parameters was calculated as the difference between peak exercise and rest.

RESULTS

Inter- and intra-observer reliability were both high for exercise stress echocardiography measurements. Compared with controls, Fontan patients had significantly lower s', e', longitudinal systolic strain and strain rate, and early diastolic longitudinal strain rate at rest and at peak exercise as well as reduced myocardial reserve.

CONCLUSIONS

Fontan patients have markedly reduced myocardial reserve during exercise. The use of exercise stress echocardiography assessment may improve the clinical management of Fontan patients.

Evaluating the Longevity of the Fontan Pathway

Children born with single ventricle physiology who undergo Fontan palliation face a diverse set of long-term complications. However, patient follow-up has in large part been limited to single institutional experiences without uniform application of diagnostic modalities to screen for relevant outcomes. Additionally, the use of different graft materials and variable surgical technique as part of the Fontan procedure has further complicated the evaluation of single ventricle patients. The purpose of this review is to define the changes in the Fontan pathway specific to the graft material used and its relationship to patient outcomes. As a means of introduction, we briefly review the historical evolution of the Fontan procedure with a focus on the intent behind design changes and incorporation of different biomaterials. We further delineate changes to the Fontan pathway which include the development of stenosis, differential growth, thrombosis, and calcification. Ultimately, the recognition of the changes noted within the Fontan pathway need to be assessed relative to their impact on patient hemodynamics, functional capacity, and Fontan-associated comorbidities.

Fluid-Structure Interaction Simulation of an Intra-Atrial Fontan Connection

Simple Summary

A fluid-structure interaction (FSI) simulation of an intra-atrial Fontan connection was performed. Power loss and pressure drop results fluctuated less during the FSI simulation than during the simulation run with rigid walls, but there were no observable differences in time-averaged pressure drop, connection power loss or hepatic flow distribution. These results suggested that employing a rigid wall is a reasonable assumption when evaluating time-averaged hemodynamic quantities of the Fontan connection under resting breath-held flow conditions.

Abstract

Total cavopulmonary connection (TCPC) hemodynamics has been hypothesized to be associated with long-term complications in single ventricle heart defect patients. Rigid wall assumption has been commonly used when evaluating TCPC hemodynamics using computational fluid dynamics (CFD) simulation. Previous study has evaluated impact of wall compliance on extra-cardiac TCPC hemodynamics using fluid-structure interaction (FSI) simulation. However, the impact of ignoring wall compliance on the presumably more compliant intra-atrial TCPC hemodynamics is not fully understood. To narrow this knowledge gap, this study aims to investigate impact of wall compliance on an intra-atrial TCPC hemodynamics. A patient-specific model of an intra-atrial TCPC is simulated with an FSI model. Patient-specific 3D TCPC anatomies were reconstructed from transverse cardiovascular magnetic resonance images. Patient-specific vessel flow rate from phase-contrast magnetic resonance imaging (MRI) at the Fontan pathway and the superior vena cava under resting condition were prescribed at the inlets. From the FSI simulation, the degree of wall deformation was compared with in vivo wall deformation from phase-contrast MRI data as validation of the FSI model. Then, TCPC flow structure, power loss and hepatic flow distribution (HFD) were compared between rigid wall and FSI simulation. There were differences in instantaneous pressure drop, power loss and HFD between rigid wall and FSI simulations, but no difference in the time-averaged quantities. The findings of this study support the use of a rigid wall assumption on evaluation of time-averaged intra-atrial TCPC hemodynamic metric under resting breath-held condition.

A novel method for hemodynamic analysis of penile erection

Measurement of blood flow velocity through the cavernosal arteries via penile color Doppler ultrasound (PDUS) is the most common objective method for the assessment of erectile function. However, in some clinical cases, this method needs to be augmented via the invasive intracavernosal pressure (ICP) measurement, which is arguably a more direct index for erectile function. The aim of this study is to develop a lumped parameter model (LPM) of the penile circulation mechanism integrated to a pulsatile, patient-specific, bi-ventricular circulation system to estimate ICP values non-invasively. PDUS data obtained from four random patients with erectile dysfunction are used to develop patient-specific LPMs. Cardiac output is estimated from the body surface area. Systemic pressure is obtained by a sphygmomanometer. Through the appropriate parameter set determined by optimization, patient-specific ICP values are predicted with only using PDUS data and validated by pre- and post-papaverine injection cavernosometry measurements. The developed model predicts the ICP with an average error value of 3 mmHg for both phases. Penile size change during erection is predicted with a ~15% error, according to the clinical size measurements. The developed mathematical model has the potential to be used as an effective non-invasive tool in erectile function evaluation, expanding the existing clinical decision parameters significantly.

Target Flow-Pressure Operating Range for Designing a Failing Fontan Cavopulmonary Support Device

Fontan operation as the current standard of care for the palliation of single ventricle defects results in significant late complications. Using a mechanical circulatory device for the right circulation to serve the function of the missing subpulmonary ventricle could potentially stabilize the failing Fontan circulation. This study aims to elucidate the hydraulic operating regions that should be targeted for designing cavopulmonary blood pumps. By integrating numerical analysis and available clinical information, the interaction of the cavopulmonary support via the IVC and full assist configurations with a wide range of simulated adult failing scenarios was investigated; with IVC and full assist corresponding to the inferior venous return or the entire venous return, respectively, being routed through the device. We identified the desired hydraulic operating regions for a cavopulmonary assist device by clustering all head pressures and corresponding pump flows that result in hemodynamic improvement for each simulated failing Fontan physiology. Results show that IVC support can produce beneficial hemodynamics in only a small fraction of failing Fontan scenarios. Cavopulmonary assist device could increase cardiac index by 35% and decrease the inferior vena cava pressure by 45% depending on the patient's pre-support hemodynamic state and surgical configuration of the cavopulmonary assist device (IVC or full support). The desired flow-pressure operating regions we identified can serve as the performance criteria for designing cavopulmonary assist devices as well as evaluating off-label use of commercially available left-side blood pumps for failing Fontan cavopulmonary support.

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.

Cavopulmonary assist: Long-term reversal of the Fontan paradox

OBJECTIVE

Fontan circulatory inefficiency can be addressed by replacing the missing subpulmonary power source to reverse the Fontan paradox. An implantable cavopulmonary assist device is described that will simultaneously reduce systemic venous pressure and increase pulmonary arterial pressure, improving preload and cardiac output, in a univentricular Fontan circulation on a long-term basis.

METHODS

A rotary blood pump that was based on the von Karman viscous pump was designed for implantation into the total cavopulmonary connection (TCPC). It will impart modest pressure energy to augment Fontan flow without risk of obstruction. In the event of rotational failure, it is designed to default to a passive flow diverter. Pressure-flow performance was characterized in vitro in a Fontan mock circulatory loop with blood analog.

RESULTS

The pump performed through the fully specified operating range, augmenting flow in all 4 directions of the TCPC. Pressure rise of 6 to 8 mm Hg was readily achieved, ranging to 14 mm Hg at highest speed (5600 rpm). Performance was consistent across a wide range of cardiac outputs. In stalled condition (0 rpm), there was no discernible pressure loss across the TCPC.

CONCLUSIONS

A blood pump technology is described that can reverse the Fontan paradox and may permit a surgical strategy of long-term biventricular maintenance of a univentricular Fontan circulation. The technology is intended for Fontan failure in which right-sided circulatory inefficiencies predominate and ventricular systolic function is preserved. It may also apply before clinical Fontan failure as health maintenance to preempt the progression of Fontan disease.

PROPOSED OVAL ST. JUDE MEDICAL VALVE: EFFECT OF HEART RATE

In this study, the hemodynamic performance of the conventional St. Jude Medical (SJM) valve and our proposed design known as the oval SJM valve are studied and compared. These studies are based on a wide range of physiological heart rates, i.e., 70–130[Formula: see text]bpm, in the opening phase. We designed and developed a precise computational platform to assess the hemodynamics of bileaflet mechanical heart valves for laminar and turbulent regimes. Also, as one of the fundamental changes applied to the conventional SJM vales, the housing is considered oval similar to oval shape of annulus. Results clearly indicate hemodynamic improvements in the proposed design over the SJM valve. The improvements are characterized by lower shear stress and wall shear stress distributions around the valve and leaflets, and lower valve pressure drop compared to that of the conventional SJM model. The proposed design shows potential and merits additional development.

Effect of Location on Tracheal Intubation Safety in Cardiac Disease-Are Cardiac ICUs Safer?

OBJECTIVES

Evaluate differences in tracheal intubation-associated events and process variances (i.e., multiple intubation attempts and oxygen desaturation) between pediatric cardiac ICUs and noncardiac PICUs in children with underlying cardiac disease.

DESIGN

Retrospective cohort study using a multicenter tracheal intubation quality improvement database (National Emergency Airway Registry for Children).

SETTING

Thirty-six PICUs (five cardiac ICUs, 31 noncardiac ICUs) from July 2012 to March 2016.

PATIENTS

Children with medical or surgical cardiac disease who underwent intubation in an ICU.

INTERVENTIONS

None.

MEASUREMENTS AND MAIN RESULTS

Our primary outcome was the rate of any adverse tracheal intubation-associated event. Secondary outcomes were severe tracheal intubation-associated events, multiple tracheal intubation attempt rates, and oxygen desaturation. There were 1,502 tracheal intubations in children with underlying cardiac disease (751 in cardiac ICUs, 751 in noncardiac ICUs) reported. Cardiac ICUs and noncardiac ICUs had similar proportions of patients with surgical cardiac disease. Patients undergoing intubation in cardiac ICUs were younger (median age, 1 mo [interquartile range, 0-6 mo]) compared with noncardiac ICUs (median 3 mo [interquartile range, 1-11 mo]; p < 0.001). Tracheal intubation-associated event rates were not different between cardiac ICUs and noncardiac ICUs (16% vs 19%; adjusted odds ratio, 0.74; 95% CI, 0.54-1.02; p = 0.069). However, in a sensitivity analysis comparing cardiac ICUs with mixed ICUs (i.e., ICUs caring for children with either general pediatric or cardiac diseases), cardiac ICUs had decreased odds of adverse events (adjusted odds ratio, 0.71; 95% CI, 0.52-0.97; p = 0.033). Rates of severe tracheal intubation-associated events and multiple attempts were similar. Desaturations occurred more often during intubation in cardiac ICUs (adjusted odds ratio, 1.61; 95% CI, 1.04-1.15; p = 0.002).

CONCLUSIONS

In children with underlying cardiac disease, rates of adverse tracheal intubation-associated events were not lower in cardiac ICUs as compared to noncardiac ICUs, even after adjusting for differences in patient characteristics and care models.

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.

Model-Based Comparison of the Normal and Fontan Circulatory Systems-Part III

BACKGROUND

For patients with the Fontan circulatory arrangement, angiotensin-converting enzyme inhibition, guanylate cyclase activation, phosphodiesterase 5 inhibition, and endothelin receptor antagonism have so far resulted in little or no improvement in [Formula: see text] or peak cardiac index (CI), suggesting that our understanding of the factors that most impact the exercise hemodynamics is incomplete.

METHODS

To facilitate comparisons with clinical reports of the exercise performance of preadolescent Fontan patients, we rescaled our previously reported computational models of a two-year-old normal child and similarly aged Fontan patient, extended our Fontan model to capture the nonlinear relationship between flow and resistance quantified from previous computational fluid dynamic analyses of the total cavopulmonary connection (TCPC), and added respiration as well as skeletal muscle contraction.

RESULTS

(1) Without respiration, the computational model for both the normal and the Fontan cannot attain the values for CI at peak exercise reported in the clinical literature, (2) because flow through the TCPC is much greater during inspiration than during expiration, the effect on the CI of the dynamic (flow-related) TCPC resistance is much more dramatic during exercise than it is in breath-hold mode at rest, and (3) coupling breathing with skeletal muscle contraction leads to the highest augmentation of cardiac output, that is, the skeletal muscle pump is most effective when the intrathoracic pressure is at a minimum-at peak inspiration.

CONCLUSIONS

Novel insights emerge when a Fontan model incorporating dynamic TCPC resistance, full respiration, and skeletal muscle contraction can be compared to the model of the normal.

Effect of Fontan geometry on exercise haemodynamics and its potential implications

OBJECTIVE

Exercise intolerance afflicts Fontan patients with total cavopulmonary connections (TCPCs) causing a reduction in quality of life. Optimising TCPC design is hypothesised to have a beneficial effect on exercise capacity. This study investigates relationships between TCPC geometries and exercise haemodynamics and performance.

METHODS

This study included 47 patients who completed metabolic exercise stress test with cardiac magnetic resonance (CMR). Phase-contrast CMR images were acquired immediately following supine lower limb exercise. Both anatomies and exercise vessel flow rates at ventilatory anaerobic threshold (VAT) were extracted. The vascular modelling toolkits were used to analyse TCPC geometries. Computational simulations were performed to quantify TCPC indexed power loss (iPL) at VAT.

RESULTS

A highly significant inverse correlation was found between the TCPC diameter index, which factors in the narrowing of TCPC vessels, with iPL at VAT (r=-0.723, p<0.001) but positive correlations with exercise performance variables, including minute oxygen consumption (VO₂) at VAT (r=0.373, p=0.01), VO₂ at peak exercise (r=0.485, p=0.001) and work at VAT/weight (r=0.368, p=0.01). iPL at VAT was negatively correlated with VO₂ at VAT (r=-0.337, p=0.02), VO₂ at peak exercise (r=-0.394, p=0.007) and work at VAT/weight (r=-0.208, p=0.17).

CONCLUSIONS

Eliminating vessel narrowing in TCPCs and reducing elevated iPL at VAT could enhance exercise tolerance for patients with TCPCs. These findings could help plan surgical or catheter-based strategies to improve patients' exercise capacity.

Hemodynamics of patient-specific aorta-pulmonary shunt configurations

Optimal hemodynamics in aorta-pulmonary shunt reconstruction is essential for improved post-operative recovery of the newborn congenital heart disease patient. However, prior to in vivo execution, the prediction of post-operative hemodynamics is extremely challenging due to the interplay of multiple confounding physiological factors. It is hypothesized that the post-operative performance of the surgical shunt can be predicted through computational blood flow simulations that consider patient size, shunt configuration, cardiac output and the complex three-dimensional disease anatomy. Utilizing only the routine patient-specific pre-surgery clinical data sets, we demonstrated an intelligent decision-making process for a real patient having pulmonary artery atresia and ventricular septal defect. For this patient, a total of 12 customized candidate shunt configurations are contemplated and reconstructed virtually using a sketch-based computer-aided anatomical editing tool. Candidate shunt configurations are evaluated based on the parameters that are computed from the flow simulations, which include 3D flow complexity, outlet flow splits, shunt patency, coronary perfusion and energy loss. Our results showed that the modified Blalock-Taussig (mBT) shunt has 12% higher right pulmonary artery (RPA) and 40% lower left pulmonary artery (LPA) flow compared to the central shunt configuration. Also, the RPA flow regime is distinct from the LPA, creating an uneven flow split at the pulmonary arteries. For all three shunt sizes, right mBT innominate and central configurations cause higher pulmonary artery (PA) flow and lower coronary artery pressure than right and left mBT subclavian configurations. While there is a trade-off between energy loss, flow split and coronary artery pressure, overall, the mBT shunts provide sufficient PA perfusion with higher coronary artery pressures and could be preferred for similar patients having PA overflow risk. Central shunts would be preferred otherwise particularly for cases with very low PA overflow risk.

Ventricular-Arterial Function and Coupling in the Adult Fontan Circulation

Background

In adult Fontan patients, ventricular or arterial dysfunction may impact homeostasis of the Fontan circulation and predispose to heart failure. We sought to characterize ventricular‐arterial (VA) properties in adult Fontan patients.

Methods and Results

Adult Fontan patients (n=170), including those with right (SRV, n=57) and left (SLV, n=113) dominant ventricular morphology, were compared to age, sex, and body size matched controls (n=170). Arterial function, load‐insensitive measures of contractility, VA coupling, diastolic function, and ventricular efficiency were assessed. Compared to controls, Fontan patients had similar arterial (Ea), but lower end‐systolic ventricular (Ees), elastance, preload recruitable stroke work and peak power index, impaired VA coupling, eccentric remodeling, reduced ventricular efficiency and increased diastolic stiffness (P<0.05 for all). Ventricular efficiency declined steeply with higher heart rate in Fontan, but not control, patients. Among Fontan patients (n=123) and controls (n=162) with preserved cardiac index (CI; ≥2.5 L/min per m²), Fontan patients had worse contractility than controls, but CI was preserved owing to relative tachycardia, lower afterload, and eccentric remodeling. However, 25% of Fontan patients had reduced CI and were distinguished from those with preserved CI by less‐eccentric remodeling and worse diastolic function, rather than more‐impaired contractility.

Conclusions

Adult Fontan patients have contractile and diastolic dysfunction with normal afterload, impaired VA coupling, and reduced ventricular efficiency with heightened sensitivity to heart rate. Maintenance of CI is dependent on lower afterload, eccentric remodeling, and relative preservation of diastolic function. These data contribute to our understanding of circulatory physiology in adult Fontan patients.

Effect of heart rate on the hemodynamics of bileaflet mechanical heart valves’ prostheses (St. Jude Medical) in the aortic position and in the opening phase: A computational study

To date, to the best of the authors’ knowledge, in almost all of the studies performed around the hemodynamics of bileaflet mechanical heart valves, a heart rate of 70–72 beats/min has been considered. In fact, the heart rate of ~72 beats/min does not represent the entire normal physiological conditions under which the aortic or prosthetic valves function. The heart rates of 120 or 50 beats/min may lead to hemodynamic complications, such as plaque formation and/or thromboembolism in patients. In this study, the hemodynamic performance of the bileaflet mechanical heart valves in a wide range of normal and physiological heart rates, that is, 60–150 beats/min, was studied in the opening phase. The model considered in this study was a St. Jude Medical bileaflet mechanical heart valve with the inner diameter of 27 mm in the aortic position. The hemodynamics of the native valve and the St. Jude Medical valve were studied in a variety of heart rates in the opening phase and the results were carefully compared. The results indicate that peak values of the velocity profile downstream of the valve increase as heart rate increases, as well as the location of the maximum velocity changes with heart rate in the St. Jude Medical valve model. Also, the maximum values of shear stress and wall shear stresses downstream of the valve are proportional to heart rate in both models. Interestingly, the maximum shear stress and wall shear stress values in both models are in the same range when heart rate is <90 beats/min; however, these values significantly increase in the St. Jude Medical valve model when heart rate is >90 beats/min (up to ~40% growth compared to that of the native valve). The findings of this study may be of importance in the hemodynamic performance of bileaflet mechanical heart valves. They may also play an important role in design improvement of conventional prosthetic heart valves and the design of the next generation of prosthetic valves, such as percutaneous valves.

Haemodynamic impact of stent implantation for lateral tunnel Fontan stenosis: a patient-specific computational assessment

BACKGROUND

The physiological importance of the lateral tunnel stenosis in the Fontan pathway for children with single ventricle physiology can be difficult to determine. The impact of the stenosis and stent implantation on total cavopulmonary connection resistance has not been characterized, and there are no clear guidelines for intervention. Methods and results A computational framework for haemodynamic assessment of stent implantation in patients with lateral tunnel stenosis was developed. Cardiac magnetic resonances images were reconstructed to obtain total cavopulmonary connection anatomies before stent implantation. Stents with 2-mm diameter increments were virtually implanted in each patient to understand the impact of stent diameter. Numerical simulations were performed in all geometries with patient-specific flow rates. Exercise conditions were simulated by doubling and tripling the lateral tunnel flow rate. The resulting total cavopulmonary connection vascular resistances were computed. A total of six patients (age: 14.4 ± 3.1 years) with lateral tunnel stenosis were included for preliminary analysis. The mean baseline resistance was 1.54 ± 1.08 WU · m(2) and dependent on the stenosis diameter. It was further exacerbated during exercise. It was observed that utilising a stent with a larger diameter lowered the resistance, but the resistance reduction diminished at larger diameters.

CONCLUSIONS

Using a computational framework to assess the severity of lateral tunnel stenosis and the haemodynamic impact of stent implantation, it was observed that stenosis in the lateral tunnel pathway was associated with higher total cavopulmonary connection resistance than unobstructed pathways, which was exacerbated during exercise. Stent implantation could reduce the resistance, but the improvement was specific to the minimum diameter.

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.

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.

Pulmonary vasodilator therapy and early postoperative outcome after modified Fontan operation

Although mortality is low after the modified Fontan procedure, there is a significant percentage of patients with prolonged postoperative recovery. The objective of this study is to evaluate the usefulness of postoperative administration of oral sildenafil and inhaled nitric oxide on early postoperative outcome. A prospective interventional and comparison study with a historical cohort was conducted. Between January, 2010 and March, 2013, 16 patients received oral sildenafil during immediate modified Fontan postoperative period. Inhaled nitric oxide was also administered if the patient was kept intubated 12 hours after surgery. Early postoperative outcome was compared with a historical cohort of 32 patients on whom the modified Fontan procedure was performed between March, 2000 and December, 2009. Postoperative administration of sildenafil and nitric oxide had no influence on early postoperative outcome after the modified Fontan procedure in terms of duration of pleural effusions, mechanical ventilation time, length of stay in the ICU, and length of hospital stay.

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.

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.

Model-Based Comparison of the Normal and Fontan Circulatory Systems—Part II

Background:

In the absence of an accessible chronic animal model of the Fontan circulation, computational modeling can provide insights into this unique circulatory arrangement, especially how differently it behaves from the normal mammalian circulation. Many groups have focused on refining a single element of the entire Fontan circulation—the total cavopulmonary connection (TCPC). Yet, only modest improvements in transplant-free survival have resulted. From an engineering perspective, optimizing the performance of a complex, multiparameter system requires an understanding of how the performance is affected by the full set of system parameters.

Methods:

We evaluated the hemodynamic impact of nine physiological perturbations in the two-year-old (yo) patient with hypoplastic left heart syndrome having a Fontan rearrangement (using our previously described lumped-parameter multicompartment model of both pulmonary and systemic circulations). In cases where comparison is appropriate, we evaluated the hemodynamic impact of analogous pathophysiologies in the normal two-year-olds. We operated the model in open-loop mode in order to expose the magnitude of the impact of uncompensated physiological perturbations.

Results:

Without the benefit of compensatory mechanisms, a valvar regurgitant fraction of 50% is sufficient to drop the cardiac index (CI) to 2.0 L/min/m2 or less. Aortopulmonary collateral flow of 0.6 L/min (1.1 L/min/m2) or 0.5 L/min (0.9 L/min/m2), sufficient to raise the ratio of pulmonary flow to systemic flow (Qp/Qs) to no higher than 1.2 or 1.5 (fenestration present or absent, respectively), is the maximum which could be tolerated (CI = 2.0 L/min/m2) without the help of compensatory mechanisms. Ventricular end-diastolic elastance (stiffness) changes have dramatic effects on CI in a Fontan circulatory arrangement.

Conclusions:

Several components of the Fontan circulation other than the TCPC actually have equal, or greater, impact on CI under certain conditions.

Computational Analysis of the Mechanical Behaviors of Hemiarch and Total Arch Replacements

Surgery for aortic dissections or aneurysms can be extended into the aortic arch by hemiarch replacement (HAR) or total arch replacement (TAR). Although cardiovascular surgeons have been performing HAR and TAR for decades, the mechanical properties of HAR and TAR are not well understood. This study investigates the mechanical behaviors and stress distributions in HAR and TAR using a hybrid fluid-structure interaction analysis that combines computational fluid dynamics and structural static analysis. Geometrical information on the aortas of 11 subjects was extracted from contrast-enhanced computed tomography (CT) scan data. The CT images were imported into medical image processing software to reconstruct 3D models of the aortas. A 3D finite element model was employed to simulate aortas that receive HAR or TAR. The deformation of the great vessels and the stress distributions at both the vessels and the aortic grafts were calculated. The numerical results revealed that the aortas following TAR exhibited a lower level of stress than those following HAR. Higher stresses may cause arterial wall injury and increase the risk of rupture. Finite element analysis of the aortas and the aortic grafts provides useful information that helps physicians better understand the potential problems that may arise after various surgical procedures.

Transition from fetal to neonatal circulation: Modeling the effect of umbilical cord clamping

Hemodynamics of the fetal to neonatal transition are orchestrated through complex physiological changes and results in cardiovascular adaptation to the adult biventricular circulation. Clinical practice during this critical period can influence vital organ physiology for normal newborns, premature babies and congenital heart defect patients. Particularly, the timing of the cord clamping procedure, immediate (ICC) vs. delayed cord clamping (DCC), is hypothesized to be an important factor for the transitory fetal hemodynamics. The clinical need for a quantitative understanding of this physiology motivated the development of a lumped parameter model (LPM) of the fetal cardio-respiratory system covering the late-gestation to neonatal period. The LPM was validated with in vivo clinical data and then used to predict the effects of cord clamping procedures on hemodynamics and vital gases. Clinical time-dependent resistance functions to simulate the vascular changes were introduced. For DCC, placental transfusion (31.3 ml) increased neonatal blood volume by 11.7%. This increased blood volume is reflected in an increase in preload pressures by ~20% compared to ICC, which in turn increased the cardiac output (CO) by 20% (COICC=993 ml/min; CODCC=1197 ml/min). Our model accurately predicted dynamic flow patterns in vivo. DCC was shown to maintain oxygenation if the onset of pulmonary respiration was delayed or impaired. On the other hand, a significant 25% decrease in oxygen saturations was observed when applying ICC under the same physiological conditions. We conclude that DCC has a significant impact on newborn hemodynamics, mainly because of the improved blood volume and the sustained placental respiration.

A simulation protocol for exercise physiology in Fontan patients using a closed loop lumped-parameter model

BACKGROUND

Reduced exercise capacity is nearly universal among Fontan patients, though its etiology is not yet fully understood. While previous computational studies have attempted to model Fontan exercise, they did not fully account for global physiologic mechanisms nor directly compare results against clinical and physiologic data.

METHODS

In this study, we developed a protocol to simulate Fontan lower-body exercise using a closed-loop lumped-parameter model describing the entire circulation. We analyzed clinical exercise data from a cohort of Fontan patients, incorporated previous clinical findings from literature, quantified a comprehensive list of physiological changes during exercise, translated them into a computational model of the Fontan circulation, and designed a general protocol to model Fontan exercise behavior. Using inputs of patient weight, height, and if available, patient-specific reference heart rate (HR) and oxygen consumption, this protocol enables the derivation of a full set of parameters necessary to model a typical Fontan patient of a given body-size over a range of physiologic exercise levels.

RESULTS

In light of previous literature data and clinical knowledge, the model successfully produced realistic trends in physiological parameters with exercise level. Applying this method retrospectively to a set of clinical Fontan exercise data, direct comparison between simulation results and clinical data demonstrated that the model successfully reproduced the average exercise response of a cohort of typical Fontan patients.

CONCLUSION

This work is intended to offer a foundation for future advances in modeling Fontan exercise, highlight the needs in clinical data collection, and provide clinicians with quantitative reference exercise physiologies for Fontan patients.

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.

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.