Physics-Driven CFD Modeling of Complex Anatomical Cardiovascular Flows?A TCPC Case Study

Clinical situations for which 3D Printing is considered an appropriate representation or extension of data contained in a medical imaging examination: vascular conditions

BACKGROUND

Medical three-dimensional (3D) printing has demonstrated utility and value in anatomic models for vascular conditions. A writing group composed of the Radiological Society of North America (RSNA) Special Interest Group on 3D Printing (3DPSIG) provides appropriateness recommendations for vascular 3D printing indications.

METHODS

A structured literature search was conducted to identify all relevant articles using 3D printing technology associated with vascular indications. Each study was vetted by the authors and strength of evidence was assessed according to published appropriateness ratings.

RESULTS

Evidence-based recommendations for when 3D printing is appropriate are provided for the following areas: aneurysm, dissection, extremity vascular disease, other arterial diseases, acute venous thromboembolic disease, venous disorders, lymphedema, congenital vascular malformations, vascular trauma, vascular tumors, visceral vasculature for surgical planning, dialysis access, vascular research/development and modeling, and other vasculopathy. Recommendations are provided in accordance with strength of evidence of publications corresponding to each vascular condition combined with expert opinion from members of the 3DPSIG.

CONCLUSION

This consensus appropriateness ratings document, created by the members of the 3DPSIG, provides an updated reference for clinical standards of 3D printing for the care of patients with vascular conditions.

Optimal Renal Artery-Aorta Angulation Revealed by Flow Simulation

INTRODUCTION

In the circulatory system, the vessel branching angle may have hemodynamic consequences. We hypothesized that there is a hemodynamically optimal range for the renal artery's branching angle.

METHODS

Data on the posttransplant kinetics of estimated glomerular filtration rate (eGFR) were analyzed according to the donor and implant sides (right-to-right and left-to-right position; n = 46). The renal artery branching angle from the aorta of a randomly selected population was measured using an X-ray angiogram (n = 44). Computational fluid dynamics simulation was used to elucidate the hemodynamic effects of angulation.

RESULTS AND DISCUSSION

Renal transplant patients receiving a right donor kidney to the right side showed faster adaptation and higher eGFR values than those receiving a left donor kidney to the right side (eGFR: 65 ± 7 vs. 56 ± 6 mL/min/1.73 m2; p < 0.01). The average branching angle on the left side was 78° and that on the right side was 66°. Simulation results showed that the pressure, volume flow, and velocity were relatively constant between 58° and 88°, indicating that this range is optimal for the kidneys. The turbulent kinetic energy does not change significantly between 58° and 78°.

CONCLUSION

The results suggest that there is an optimal range for the renal artery's branching angle from the aorta where hemodynamic vulnerability caused by the degree of angulation is the lowest, which should be considered during kidney transplantations.

In-Silico and In-Vitro Analysis of the Novel Hybrid Comprehensive Stage II Operation for Single Ventricle Circulation

Single ventricle (SV) anomalies account for one-fourth of all congenital heart disease cases. The existing palliative treatment for this anomaly achieves a survival rate of only 50%. To reduce the trauma associated with surgical management, the hybrid comprehensive stage II (HCSII) operation was designed as an alternative for a select subset of SV patients with the adequate antegrade aortic flow. This study aims to provide better insight into the hemodynamics of HCSII patients utilizing a multiscale Computational Fluid Dynamics (CFD) model and a mock flow loop (MFL). Both 3D-0D loosely coupled CFD and MFL models have been tuned to match baseline hemodynamic parameters obtained from patient-specific catheterization data. The hemodynamic findings from clinical data closely match the in-vitro and in-silico measurements and show a strong correlation (r = 0.9). The geometrical modification applied to the models had little effect on the oxygen delivery. Similarly, the particle residence time study reveals that particles injected in the main pulmonary artery (MPA) have successfully ejected within one cardiac cycle, and no pathological flows were observed.

Numerical Investigation of the Fetal Left Heart Hemodynamics During Gestational Stages

Flow-driven hemodynamic forces on the cardiac tissues have critical importance, and have a significant role in the proper development of the heart. These mechanobiological mechanisms govern the cellular responses for the growth and remodeling of the heart, where the altered hemodynamic environment is believed to be a major factor that is leading to congenital heart defects (CHDs). In order to investigate the mechanobiological development of the normal and diseased hearts, identification of the blood flow patterns and wall shear stresses (WSS) on these tissues are required for an accurate hemodynamic assessment. In this study, we focus on the left heart hemodynamics of the human fetuses throughout the gestational stages. Computational fetal left heart models are created for the healthy fetuses using the ultrasound images at various gestational weeks. Realistic inflow boundary conditions are implemented in the models using the Doppler ultrasound measurements for resolving the specific blood flow waveforms in the mitral valve. Obtained results indicate that WSS and vorticity levels in the fetal left heart decrease with the development of the fetus. The maximum WSS around the mitral valve is determined around 36 Pa at the gestational week of 16. This maximum WSS decreases to 11 Pa at the gestational week of 27, indicating nearly three-times reduction in the peak shear stress. These findings reveal the highly dynamic nature of the left heart hemodynamics throughout the development of the human fetus and shed light into the relevance of hemodynamic environment and development of CHDs.

Experimental Validation of Enhanced Magnetic Resonance Imaging (EMRI) Using Particle Image Velocimetry (PIV)

Flow-sensitive four-dimensional Cardiovascular Magnetic Resonance Imaging (4D Flow CMR) has increasingly been utilised to characterise patients' blood flow, in association with patiens' state of health and disease, even though spatial and temporal resolutions still constitute a limit. Computational fluid dynamics (CFD) is a powerful tool that could expand these information and, if integrated with experimentally-obtained velocity fields, would enable to derive a large variety of the flow descriptors of interest. However, the accuracy of the flow parameters is highly influenced by the quality of the input data such as the anatomical model and boundary conditions typically derived from medical images including 4D Flow CMR. We previously proposed a novel approach in which 4D Flow CMR and CFD velocity fields are integrated to obtain an Enhanced 4D Flow CMR (EMRI), allowing to overcome the spatial-resolution limitation of 4D Flow CMR, and enable an accurate quantification of flow. In this paper, the proposed approach is validated in a U bend channel, an idealised model of the human aortic arch. The flow patterns were studied with 4D Flow CMR, CFD and EMRI, and compared with high resolution 2D PIV experiments obtained in pulsatile conditions. The main strengths and limitations of 4D Flow CMR and CFD were illustrated by exploiting the accuracy of PIV by comparing against PIV velocity fields. EMRI flow patterns showed a better qualitative and quantitative agreement with PIV results than the other techniques. EMRI enables to overcome the experimental limitations of MRI-based velocity measurements and the modelling simplifications of CFD, allowing an accurate prediction of complex flow patterns observed experimentally, while satisfying mass and momentum balance equations.

An Anterior Anastomosis for the Modified Fontan Connection: A Hemodynamic Analysis

This hemodynamic feasibility study examined total cavopulmonary connection (TCPC) designs connecting the extracardiac conduit to the anterior surface of pulmonary arteries (PAs) or superior vena cava (SVC) rather than to the inferior PA surface (traditional TCPC). The study involved twenty-five consecutive Fontan patients meeting inclusion criteria from a single institution. A virtual surgical platform mimicked the completed traditional TCPC and generated three anterior anastomosis designs: Anterior-PA, Middle-SVC, and SVC-Inn (Inn: innominate vein). Hemodynamic performance of anterior anastomosis designs was compared with the traditional TCPC regarding indexed power loss (iPL) and hepatic flow distribution (HFD). Compared to the traditional TCPC, the Anterior-PA design produces a similar iPL. The Middle-SVC design is also similar, though the iPL difference is positively correlated with the anastomosing height. The SVC-Inn design had significantly more iPL. The three anterior anastomosis designs did not have a significant difference in HFD (from traditional TCPC). Pulmonary flow distribution (PFD) has a stronger correlation with HFD from the anterior anastomosis designs than the traditional TCPC. This hemodynamic feasibility study examined anterior anastomosis, extracardiac TCPC designs that may offer surgeons clinical dexterity. The Anterior-PA design may be equivalent to the traditional TCPC. Fontan extracardiac conduit anastomosis just superior to the PAs (Middle-SVC) also preserves hemodynamic performance and avoids direct PA anastomosis. These designs could simplify surgical Fontan completion, and may particularly benefit patients requiring surgical dissection, having atypical PA orientation, or after PA stent angioplasty.

Assessment of biventricular hemodynamics and energy dynamics using lumen-tracking 4D flow MRI without contrast medium

BACKGROUND

Biventricular physiological interaction remains a challenging problem in cardiology. We developed a four-dimensional (4D) flow magnetic resonance imaging (MRI) scan and clinically available analysis protocol based on beat tracking of the cardiovascular lumen without contrast medium, which enabled measurement of the biventricular hemodynamics and energetic performance by calculating flow energy loss (EL) and kinetic energy (KE). The aim of this study was to observe the flow patterns and energy dynamics to reveal the physiology of the right and left ventricular systems.

METHODS

4D flow MRI studies were performed in 19 healthy volunteers including 11 male and 8 female. The right and left ventricular systems were segmented to visualize the flow patterns and to quantify the hemodynamics and energy dynamics.

RESULTS

A large vortex was observed in the left ventricle (LV), along the longitudinal axis, during end diastole and early systole. At early systole, the vortex appeared to facilitate smooth ejection with little EL. In contrast, in the right ventricle (RV), there were vortices near the free wall in both the short and long axes during the diastolic filling phase. Mean EL index during a single cardiac cycle in the right and left heart systems was 0.63 ± 0.16 (0.42-0.99) mW/m², and 1.02 ± 0.26 (0.58-1.58) mW/m², respectively. EL is inevitable loss caused by the vortex flow to facilitate smooth right and left ventricular function and left-sided EL tended to correlate positively with heart rate and right ventricular stroke volume. Kinetic energy at the aortic valve was influenced by LV end-diastolic volume/stroke volume. No gender difference was observed.

CONCLUSIONS

The RV appears to function as a regulator of the energy dynamics of the LV system.

Computerized virtual surgery based on computational fluid dynamics simulation for planning coronary revascularization with aortic root replacement in adult congenital heart disease: a case report

A 38-year-old woman presented with exertional dyspnea and chest compression. She had undergone repair of congenital supravalvular aortic stenosis at 8 years of age. Contrast-enhanced computed tomography showed re-stenosis in the ascending aorta, bilateral coronary arterial aneurysm, and a highly thickened left ventricular wall. Release of stenosis was necessary to avoid left ventricular functional deterioration; however, it could cause demand–supply mismatch in coronary flow due to substantial left ventricular hypertrophy. Sufficient statistical evidence was not available in this situation; therefore, computerized virtual surgery based on computational fluid dynamics (CFD) was performed to predict the postoperative hemodynamics. Consequently, root replacement with in situ Carrel patch coronary reconstruction was considered a better option than coronary artery graft bypass in the left-side coronary flow supply. The patient underwent root replacement with in situ Carrel patch coronary reconstruction as planned based on CFD without any complication and was discharged 15 days postoperatively.

Methodology for Hemodynamic Assessment of a Three-Dimensional Printed Patient-Specific Vascular Test Device

Assessing hemodynamics in vasculature is important for the development of cardiovascular diagnostic parameters and evaluation of medical devices. Benchtop experiments are a safe and comprehensive preclinical method for testing new diagnostic endpoints and devices within a controlled environment. Recent advances in three-dimensional (3D) printing have enhanced benchtop tests by allowing generation of patient-specific and pathophysiologic conditions. We used 3D printing, coupled with image processing and computer-aided design (CAD), to develop a patient-specific vascular test device from clinical data. The proximal pulmonary artery (PA) tree including the main, left, and right pulmonary arteries, with a stenosis within the left PA was selected as a representative anatomy for developing the vascular test device. Three test devices representing clinically relevant stenosis severities, 90%, 80%, and 70% area stenosis, were evaluated at different cardiac outputs (COs). A mock circulatory loop (MCL) generating pathophysiologic pulmonary pressure and flow was used to evaluate the hemodynamics within the devices. The dimensionless pressure drop–velocity ratio characteristic curves for the three stenosis severities were obtained. At a fixed CO, the dimensionless pressure drop increased nonlinearly with an increase in (a) the velocity ratio for a fixed stenosis severity and (b) the stenosis severity at a specific velocity ratio. The dimensionless pressure drop observed in vivo was similar (within 1%) to that measured in moderate area stenosis of 70% because both flows were viscous dominated. The hemodynamics of the 3D printed test device can be used for evaluating diagnostic endpoints and medical devices in a preclinical setting under realistic conditions.

Analysis of Inlet Velocity Profiles in Numerical Assessment of Fontan Hemodynamics

Computational fluid dynamic (CFD) simulations are widely utilized to assess Fontan hemodynamics that are related to long-term complications. No previous studies have systemically investigated the effects of using different inlet velocity profiles in Fontan simulations. This study implements real, patient-specific velocity profiles for numerical assessment of Fontan hemodynamics using CFD simulations. Four additional, artificial velocity profiles were used for comparison: (1) flat, (2) parabolic, (3) Womersley, and (4) parabolic with inlet extensions [to develop flow before entering the total cavopulmonary connection (TCPC)]. The differences arising from the five velocity profiles, as well as discrepancies between the real and each of the artificial velocity profiles, were quantified by examining clinically important metrics in TCPC hemodynamics: power loss (PL), viscous dissipation rate (VDR), hepatic flow distribution, and regions of low wall shear stress. Statistically significant differences were observed in PL and VDR between simulations using real and flat velocity profiles, but differences between those using real velocity profiles and the other three artificial profiles did not reach statistical significance. These conclusions suggest that the artificial velocity profiles (2)-(4) are acceptable surrogates for real velocity profiles in Fontan simulations, but parabolic profiles are recommended because of their low computational demands and prevalent applicability.

Evaluation of pulmonary artery stenosis in congenital heart disease patients using functional diagnostic parameters: An in vitro study

Congenital pulmonary artery (PA) stenosis is often associated with abnormal PA hemodynamics including increased pressure drop (Δp) and reduced asymmetric flow (Q), which may result in right ventricular dysfunction. We propose functional diagnostic parameters, pressure drop coefficient (CDP), energy loss (E_loss), and normalized energy loss (E¯_loss) to characterize pulmonary hemodynamics, and evaluate their efficacy in delineating stenosis severity using in vitro experiments. Subject-specific test sections including the main PA (MPA) bifurcating into left and right PAs (LPA, RPA) with a discrete LPA stenosis were manufactured from cross-sectional imaging and 3D printing. Three clinically-relevant stenosis severities, 90% area stenosis (AS), 80% AS, and 70% AS, were evaluated at different cardiac outputs (COs). A benchtop flow loop simulating pulmonary hemodynamics was used to measure Q and Δp within the test sections. The experimental Δp-Q characteristics along with clinical data were used to obtain pathophysiologic conditions and compute the diagnostic parameters. The pathophysiologic Q_LPA decreased as the stenosis severity increased at a fixed CO. CDP_LPA, E_loss,LPA (absolute), and E¯_loss,LPA (absolute) increased with an increase in LPA stenosis severity at a fixed CO. Importantly, CDP_LPA and E¯_loss,LPA had reduced variability with CO, and distinct values for each LPA stenosis severity. Under variable CO, a) CDP_LPA values were 14.5-21.0 (70% AS), 60.7- 2.2 (80% AS), ≥ 261.6 (90% AS), and b) E¯_loss,LPA values (in mJ per Q_LPA) were -501.9 to -1023.8 (70% AS), -1247.6 to -1773.0 (80% AS), -1934.5 (90% AS). Hence, CDP_LPA and E¯_loss,LPA are expected to assess the true functional severity of PA stenosis.

The effect of inlet and outlet boundary conditions in image-based CFD modeling of aortic flow

Background

Computational modeling of cardiovascular flow is a growing and useful field, but such simulations usually require the researcher to guess the flow’s inlet and outlet conditions since they are difficult and expensive to measure. It is critical to determine the amount of uncertainty introduced by these assumptions in order to evaluate the degree to which cardiovascular flow simulations are accurate. Our work begins to address this question by examining the sensitivity of flow to several different assumed velocity inlet and outlet conditions in a patient-specific aorta model.

Methods

We examined the differences between plug flow, parabolic flow, linear shear flows, skewed cubic flow profiles, and Womersley flow at the inlet. Only the shape of the inlet velocity profile was varied—all other parameters were identical among these simulations. Secondary flow in the form of a counter-rotating pair of vortices was also added to parabolic axial flow to study its effect on the solution. In addition, we examined the differences between two-element Windkessel, three element Windkessel and the outflow boundary conditions. In these simulations, only the outlet boundary condition was varied.

Results

The results show axial and in-plane velocities are considerably different close to the inlet for the cases with different inlet velocity profile shapes. However, the solutions are qualitatively similar beyond 1.75D, where D is the inlet diameter. This trend is also observed in other quantities such as pressure and wall shear stress. Normalized root-mean-square deviation, a measure of axial velocity magnitude differences between the different cases, generally decreases along the streamwise coordinate. The linear shear inlet velocity boundary condition and plug velocity boundary condition solution exhibit the highest time-averaged wall shear stress, approximately \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$8\%$$\end{document}8% higher than the parabolic inlet velocity boundary condition. Upstream of 1D from the inlet, adding secondary flow has a significant impact on temporal wall shear stress distributions. This is especially observable during diastole, when integrated wall shear stress magnitude varies about \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$26\%$$\end{document}26% between simulations with and without secondary flow. The results from the outlet boundary condition study show the Windkessel models differ from the outflow boundary condition by as much as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$18\%$$\end{document}18% in terms of time-averaged wall shear stress. Furthermore, normalized root-mean-square deviation of axial velocity magnitude, a measure of deviation between Windkessel and the outflow boundary condition, increases along the streamwise coordinate indicating larger variations near outlets.

Conclusion

It was found that the selection of inlet velocity conditions significantly affects only the flow region close to the inlet of the aorta. Beyond two diameters distal to the inlet, differences in flow solution are small. Although additional studies must be performed to verify this result, the data suggest that it is important to use patient-specific inlet conditions primarily if the researcher is concerned with the details of the flow very close to the inlet. Similarly, the selection of outlet conditions significantly affects the flow in the vicinity of the outlets. Upstream of five diameters proximal to the outlet, deviations between the outlet boundary conditions examined are insignificant. Although the inlet and outlet conditions only affect the flow significantly in their respective neighborhoods, our study indicates that outlet conditions influence a larger percentage of the solution domain.

Mechanical axial flow blood pump to support cavopulmonary circulation

We are developing a collapsible, percutaneously inserted, axial flow blood pump to support the cavopulmonary circulation in infants with a failing single ventricle physiology. An initial design of the impeller for this axial flow blood pump was performed using computational fluid dynamics analysis, including pressure-flow characteristics, scalar stress estimations, blood damage indices, and fluid force predictions. A plastic prototype was constructed for hydraulic performance testing, and these experimental results were compared with the numerical predictions. The numerical predictions and experimental findings of the pump performance demonstrated a pressure generation of 2–16 mm Hg for 50–750 ml/min over 5,500–7,500 RPM with deviation found at lower rotational speeds. The axial fluid forces remained below 0.1 N, and the radial fluid forces were determined to be virtually zero due to the centered impeller case. The scalar stress levels remained below 250 Pa for all operating conditions. Blood damage analysis yielded a mean residence time of the released particles, which was found to be less than 0.4 seconds for both flow rates that were examined, and a maximum residence time was determined to be less than 0.8 seconds. We are in the process of designing a cage with hydrodynamically shaped filament blades to act as a diffuser and optimizing the impeller blade shape to reduce the flow vorticity at the pump outlet. This blood pump will improve the clinical treatment of patients with failing Fontan physiology and provide a unique catheter-based therapeutic approach as a bridge to recovery or transplantation.

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

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

Study of the transdentinal diffusion of bioactive molecules

In this work the mass transfer characteristics in a µ-tube that simulates a simplified dentinal tubule geometry are numerically investigated. The aim is to assess the key features that affect transdentinal diffusion of substances and consequently to define the necessary quantitative and qualitative issues related to a specific bioactive agent before its potential application in clinical practice. CFD simulations were performed in an S-shaped tapered micro-tube, while the code was validated using the non-intrusive optical measuring technique Laser Induced Fluorescence (LIF). As the phenomenon is one-dimensional, diffusion dominated and strongly dependent on the molecular size, the time needed for the concentration of released molecules to attain a required value can be controlled by their initial concentration. Thus, we propose a model, which is successfully verified by experimental data using a dentinal disc and which given the type of applied molecules and their critical pulpal concentration is able to estimate the initial concentration to be imposed.

Synchronization of the Flow and Pressure Waves Obtained With Non-Simultaneous Multipoint Measurements

The use of measured data as boundary conditions renders hemodynamic simulations more patient-specific. However, synchronized acquisition of data at multiple locations is often difficult in clinical practice. This study proposes a method for resynchronizing measured data for use as boundary conditions for flow simulations using frequency analyses, and discusses the optimal cut-off frequency for differentiating cardiac and respiratory variation in hemodynamic data during resynchronization. To demonstrate the utility of the method, a Fontan circulation, which is the final palliative result with single-ventricle physiology, was used. The results suggest that it is optimal to set a cut-off frequency that gives a local minimum in the power spectrum that is slightly lower than the peak frequency of the heartbeat. Additionally, the total energy loss depended on the cut-off frequency, although the overall flow patterns appeared to be similar. The method is applicable to cardiovascular systems other than the Fontan circulation, where hemodynamic data with multifactorial fluctuations are required at various locations but simultaneous measurements are not possible.

A systematic review of image segmentation methodology, used in the additive manufacture of patient-specific 3D printed models of the cardiovascular system

Background

Shortcomings in existing methods of image segmentation preclude the widespread adoption of patient-specific 3D printing as a routine decision-making tool in the care of those with congenital heart disease. We sought to determine the range of cardiovascular segmentation methods and how long each of these methods takes.

Methods

A systematic review of literature was undertaken. Medical imaging modality, segmentation methods, segmentation time, segmentation descriptive quality (SDQ) and segmentation software were recorded.

Results

Totally 136 studies met the inclusion criteria (1 clinical trial; 80 journal articles; 55 conference, technical and case reports). The most frequently used image segmentation methods were brightness thresholding, region growing and manual editing, as supported by the most popular piece of proprietary software: Mimics (Materialise NV, Leuven, Belgium, 1992–2015). The use of bespoke software developed by individual authors was not uncommon. SDQ indicated that reporting of image segmentation methods was generally poor with only one in three accounts providing sufficient detail for their procedure to be reproduced.

Conclusions and implication of key findings

Predominantly anecdotal and case reporting precluded rigorous assessment of risk of bias and strength of evidence. This review finds a reliance on manual and semi-automated segmentation methods which demand a high level of expertise and a significant time commitment on the part of the operator. In light of the findings, we have made recommendations regarding reporting of 3D printing studies. We anticipate that these findings will encourage the development of advanced image segmentation methods.

Effects of intraluminal thrombus on patient-specific abdominal aortic aneurysm hemodynamics via stereoscopic particle image velocity and computational fluid dynamics modeling

The pathology of the human abdominal aortic aneurysm (AAA) and its relationship to the later complication of intraluminal thrombus (ILT) formation remains unclear. The hemodynamics in the diseased abdominal aorta are hypothesized to be a key contributor to the formation and growth of ILT. The objective of this investigation is to establish a reliable 3D flow visualization method with corresponding validation tests with high confidence in order to provide insight into the basic hemodynamic features for a better understanding of hemodynamics in AAA pathology and seek potential treatment for AAA diseases. A stereoscopic particle image velocity (PIV) experiment was conducted using transparent patient-specific experimental AAA models (with and without ILT) at three axial planes. Results show that before ILT formation, a 3D vortex was generated in the AAA phantom. This geometry-related vortex was not observed after the formation of ILT, indicating its possible role in the subsequent appearance of ILT in this patient. It may indicate that a longer residence time of recirculated blood flow in the aortic lumen due to this vortex caused sufficient shear-induced platelet activation to develop ILT and maintain uniform flow conditions. Additionally, two computational fluid dynamics (CFD) modeling codes (Fluent and an in-house cardiovascular CFD code) were compared with the two-dimensional, three-component velocity stereoscopic PIV data. Results showed that correlation coefficients of the out-of-plane velocity data between PIV and both CFD methods are greater than 0.85, demonstrating good quantitative agreement. The stereoscopic PIV study can be utilized as test case templates for ongoing efforts in cardiovascular CFD solver development. Likewise, it is envisaged that the patient-specific data may provide a benchmark for further studying hemodynamics of actual AAA, ILT, and their convolution effects under physiological conditions for clinical applications.

In Vitro Validation of Patient-Specific Hemodynamic Simulations in Coronary Aneurysms Caused by Kawasaki Disease

To perform experimental validation of computational fluid dynamics (CFD) applied to patient specific coronary aneurysm anatomy of Kawasaki disease. We quantified hemodynamics in a patient-specific coronary artery aneurysm physical phantom under physiologic rest and exercise flow conditions. Using phase contrast MRI (PCMRI), we acquired 3-component flow velocity at two slice locations in the aneurysms. We then performed numerical simulations with the same geometry and inflow conditions, and performed qualitative and quantitative comparisons of velocities between experimental measurements and simulation results. We observed excellent qualitative agreement in flow pattern features. The quantitative spatially and temporally varying differences in velocity between PCMRI and CFD were proportional to the flow velocity. As a result, the percent discrepancy between simulation and experiment was relatively constant regardless of flow velocity variations. Through 1D and 2D quantitative comparisons, we found a 5-17% difference between measured and simulated velocities. Additional analysis assessed wall shear stress differences between deformable and rigid wall simulations. This study demonstrated that CFD produced good qualitative and quantitative predictions of velocities in a realistic coronary aneurysm anatomy under physiological flow conditions. The results provide insights on factors that may influence the level of agreement, and a set of in vitro experimental data that can be used by others to compare against CFD simulation results. The findings of this study increase confidence in the use of CFD for investigating hemodynamics in the specialized anatomy of coronary aneurysms. This provides a basis for future hemodynamics studies in patient-specific models of Kawasaki disease.

INFLUENCE OF EXERCISE ON THREE TYPES OF FONTAN MODIFICATIONS: COMPARISON BY NUMERICAL SIMULATIONS

Low energy efficiency of Fontan-type connections in exercise conditions may lead to limited post-operative exercise capacity of patients. The objective of this study was to compare the energy efficiency of different types of Fontan modifications in exercise conditions. In this study, the flow rate of the inferior vena cava (IVC) was increased artificially to simulate exercise conditions. The control volume power loss and energy efficiency were investigated in three types of Fontan modifications using the method of computational fluid dynamics (CFD). In the intracardiac tunnel Fontan, the power loss was 7.9 times of that calculated in resting state when the flow rate of IVC increased to three times of that measured at the right pulmonary artery flow ratio of 50%. The energy efficiency of it decreased prominently from 91.9% to 79.9%. On the contrary, in the extracardiac Fontan with direct cavopulmonary anastomosis (direct Fontan), the magnitude of power loss was only 2.9 times of that in resting state and the energy efficiency changed only from 91.8% to 89.1%. The results of extracardiac conduit Fontan fell in between. It showed that direct Fontan had higher energy efficiency than intracardiac tunnel or extracardiac conduit Fontan in exercise conditions.

Computational fluid dynamics simulation to evaluate aortic coarctation gradient with contrast-enhanced CT

Coarctation of aorta (CoA) is a narrowing of the aorta leading to a pressure gradient (ΔP) across the coarctation, increased afterload and reduced peripheral perfusion pressures. Indication to invasive treatment is based on values of maximal (systolic) trans-coarctation ΔP. A computational fluid dynamic (CFD) approach is herein presented for the non-invasive haemodynamic assessment of ΔP across CoA. Patient-specific CFD simulations were created from contrast-enhanced computed tomography (CT) and appropriate flow boundary conditions. Computed ΔP was validated with invasive intravascular trans-CoA pressure measurements. Haemodynamic indices, including pressure loss coefficient (PLc), time-averaged wall shear stress (TAWSS) and oscillatory shear index (OSI), were also quantified. CFD-estimated ΔP values were comparable to the invasive ones. Moreover, the aorta proximal to CoA was exposed to altered TAWSS and OSI suggesting hypertension. PLc was found as a further geometric marker of CoA severity. Finally, CFD-estimated ΔP confirmed a significant reduction after percutaneous balloon dilatation and stenting of the CoA in one patient (e.g. from ΔP∼52 mmHg to ΔP∼3 mmHg). The validation of the ΔP computations with catheterisation measurements suggests that CFD simulation, based on CT-derived anatomical data, is a useful tool to readily quantify CoA severity.

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

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

Numerical blood flow simulation in surgical corrections: what do we need for an accurate analysis?

BACKGROUND

Computational fluid dynamics has been increasingly used in congenital heart surgery to simulate pathophysiological blood flow, investigate surgical options, or design medical devices. Several commercial and research computational or numerical codes have been developed. They present different approaches to numerically solve the blood flow equations, raising the question whether these numerical codes are equally reliable to achieve accurate simulation results. Accordingly, we sought to examine the influence of numerical code selection in several complex congenital cardiac operations.

MATERIAL AND METHODS

The main steps of blood flow simulations are detailed (geometrical mesh, boundary conditions, and solver numerical methods) for congenital cardiac operations of increasing complexity. The first case tests different numerical solutions against an analytical, or exact, solution. In the second case, the three-dimensional domain is a patient-specific superior cavopulmonary anastomosis. As an analytical solution does not exist in such a complex geometry, different numerical solutions are compared. Finally, a realistic case of a systemic-to-pulmonary shunt is presented with both geometrically and physiologically challenging conditions. For all, solutions from a commercially available code and an open-source research code are compared.

RESULTS

In the first case, as the mesh or solver numerical method is refined, the simulation results for both codes converged to the analytical solution. In the second example, velocity differences between the two codes are greater when the resolution of the mesh were lower and less refined. The third case with realistic anatomy reveals that the pulsatile complex flow is very similar for both codes.

CONCLUSIONS

The precise setup of the numerical cases has more influence on the results than the choice of numerical codes. The need for detailed construction of the numerical model that requires high computational cost depends on the precision needed to answer the biomedical question at hand and should be assessed for each problem on a combination of clinically relevant patient-specific geometry and physiological conditions.

Irrigation of human prepared root canal--ex vivo based computational fluid dynamics analysis

AIM

To analyze the influence of the needle type, insertion depth, and irrigant flow rate on irrigant flow pattern, flow velocity, and apical pressure by ex-vivo based endodontic irrigation computational fluid dynamics (CFD) analysis.

METHODS

Human upper canine root canal was prepared using rotary files. Contrast fluid was introduced in the root canal and scanned by computed tomography (CT) providing a three-dimensional object that was exported to the computer-assisted design (CAD) software. Two probe points were established in the apical portion of the root canal model for flow velocity and pressure measurement. Three different CAD models of 27G irrigation needles (closed-end side-vented, notched open-end, and bevel open-end) were created and placed at 25, 50, 75, and 95% of the working length (WL). Flow rates of 0.05, 0.1, 0.2, 0.3, and 0.4 mL/s were simulated. A total of 60 irrigation simulations were performed by CFD fluid flow solver.

RESULTS

Closed-end side-vented needle required insertion depth closer to WL, regarding efficient irrigant replacement, compared to open-end irrigation needle types, which besides increased velocity produced increased irrigant apical pressure. For all irrigation needle types and needle insertion depths, the increase of flow rate was followed by an increased irrigant apical pressure.

CONCLUSIONS

The human root canal shape obtained by CT is applicable in the CFD analysis of endodontic irrigation. All the analyzed values -irrigant flow pattern, velocity, and pressure - were influenced by irrigation needle type, as well as needle insertion depth and irrigant flow rate.

Fontan conversion templates: patient-specific hemodynamic performance of the lateral tunnel versus the intraatrial conduit with fenestration

Intraatrial-conduit Fontan is considered a modification of both extracardiac and lateral-tunnel Fontan. In this study, the patient-specific hemodynamic performance of intraatrial-conduit and lateral-tunnel Fontan with fenestration, considered as conversion templates, was investigated based on the authors' patient cohort. Pulsatile computational fluid dynamics simulations were performed using patient-specific models of intraatrial-conduit and lateral-tunnel Fontan patients. Real-time "simultaneous" inferior and superior vena cava, pulmonary artery, and fenestration flow waveforms were acquired from ultrasound. Multiple hemodynamic performance indices were investigated, with particular focus on evaluation of the pulsatile flow performance. Power loss inside the lateral-tunnel Fontan appeared to be significantly higher than with the intraatrial-conduit Fontan for patient-specific cardiac output and normalized connection size. Inclusion of the 4-mm fenestration at a 0.24 L/min mean flow resulted in a lower cavopulmonary pressure gradient and less time-averaged power loss for both Fontan connections. Flow structures within the intraatrial conduit were notability more uniform than within the lateral tunnel. Hepatic flow majorly favored the left lung in both surgical connections: conversion from lateral-tunnel to intraatrial-conduit Fontan resulted in better hemodynamics with less power loss, a lower pressure gradient, and fewer stagnant flow zones along the conduit. This patient-specific computational case study demonstrated superior hemodynamics of intraatrial-conduit Fontan over those of lateral-tunnel Fontan with or without fenestration and improved performance after conversion of the lateral tunnel to the intraatrial conduit. The geometry-specific effect of the nonuniform hepatic flow distribution may motivate new rationales for the surgical design.

Computational fluid dynamics models and congenital heart diseases

Mathematical modeling is a powerful tool to investigate hemodynamics of the circulatory system. With improving imaging techniques and detailed clinical investigations, it is now possible to construct patient-specific models of reconstructive surgeries for the treatment of congenital heart diseases. These models can help clinicians to better understand the hemodynamic behavior of different surgical options for a treated patient. This review outlines recent advances in mathematical modeling in congenital heart diseases, the discoveries and limitations these models present, and future directions that are on the horizon.

An integrated approach to patient-specific predictive modeling for single ventricle heart palliation

In patients with congenital heart disease and a single ventricle (SV), ventricular support of the circulation is inadequate, and staged palliative surgery (usually 3 stages) is needed for treatment. In the various palliative surgical stages individual differences in the circulation are important and patient-specific surgical planning is ideal. In this study, an integrated approach between clinicians and engineers has been developed, based on patient-specific multi-scale models, and is here applied to predict stage 2 surgical outcomes. This approach involves four distinct steps: (1) collection of pre-operative clinical data from a patient presenting for SV palliation, (2) construction of the pre-operative model, (3) creation of feasible virtual surgical options which couple a three-dimensional model of the surgical anatomy with a lumped parameter model (LPM) of the remainder of the circulation and (4) performance of post-operative simulations to aid clinical decision making. The pre-operative model is described, agreeing well with clinical flow tracings and mean pressures. Two surgical options (bi-directional Glenn and hemi-Fontan operations) are virtually performed and coupled to the pre-operative LPM, with the hemodynamics of both options reported. Results are validated against postoperative clinical data. Ultimately, this work represents the first patient-specific predictive modeling of stage 2 palliation using virtual surgery and closed-loop multi-scale modeling.

Numerical and experimental investigation of pulsatile hemodynamics in the total cavopulmonary connection

Computational fluid dynamics (CFD) tools have been extensively applied to study the hemodynamics in the total cavopulmonary connection (TCPC) in patients with only a single functioning ventricle. Without the contraction of a sub-pulmonary ventricle, pulsatility of flow through this connection is low and variable across patients, which is usually neglected in most numerical modeling studies. Recent studies suggest that such pulsatility can be non-negligible and can be important in hemodynamic predictions. The goal of this work is to compare the results of an in-house numerical methodology for simulating pulsatile TCPC flow with experimental results. Digital particle image velocimetry (DPIV) was acquired on TCPC in vitro models to evaluate the capability of the CFD tool in predicting pulsatile TCPC flow fields. In vitro hemodynamic measurements were used to compare the numerical prediction of power loss across the connection. The results demonstrated the complexity of the pulsatile TCPC flow fields and the validity of the numerical approach in simulating pulsatile TCPC flow dynamics in both idealized and complex patient specific models.

Mechanotransduction in embryonic vascular development

A plethora of biochemical signals provides spatial and temporal cues that carefully orchestrate the complex process of vertebrate embryonic development. The embryonic vasculature develops not only in the context of these biochemical cues, but also in the context of the biomechanical forces imparted by blood flow. In the mature vasculature, different blood flow regimes induce distinct genetic programs, and significant progress has been made toward understanding how these forces are perceived by endothelial cells and transduced into biochemical signals. However, it cannot be assumed that paradigms that govern the mature vasculature are pertinent to the developing embryonic vasculature. The embryonic vasculature can respond to the mechanical forces of blood flow, and these responses are critical in vascular remodeling, certain aspects of sprouting angiogenesis, and maintenance of arterial-venous identity. Here, we review data regarding mechanistic aspects of endothelial cell mechanotransduction, with a focus on the response to shear stress, and elaborate upon the multifarious effects of shear stress on the embryonic vasculature. In addition, we discuss emerging predictive vascular growth models and highlight the prospect of combining signaling pathway information with computational modeling. We assert that correlation of precise measurements of hemodynamic parameters with effects on endothelial cell gene expression and cell behavior is required for fully understanding how blood flow-induced loading governs normal vascular development and shapes congenital cardiovascular abnormalities.

Pulsatile venous waveform quality affects the conduit performance in functional and "failing" Fontan circulations

OBJECTIVE

To investigate the effect of pulsatility of venous flow waveform in the inferior and superior caval vessels on the performance of functional and "failing" Fontan patients based on two primary performance measures - the conduit power loss and the distribution of inferior caval flow (hepatic factors) to the lungs.

METHODS

Doppler angiography flows were acquired from two typical extra-cardiac conduit "failing" Fontan patients, aged 13 and 25 years, with ventricle dysfunction. Using computational fluid dynamics, haemodynamic efficiencies of "failing", functional, and in vitro-generated mechanically assisted venous flow waveforms were evaluated inside an idealised total cavopulmonary connection with a caval offset. To investigate the effect of venous pulsatility alone, cardiac output was normalised to 3 litres per minute in all cases. To quantify the pulsatile behaviour of venous flows, two new performance indices were suggested.

RESULTS

Variations in the pulsatile content of venous waveforms altered the conduit efficiency notably. High-frequency and low-amplitude oscillations lowered the pulsatile component of the power losses in "failing" Fontan flow waveforms. Owing to the offset geometry, hepatic flow distribution depended strongly on the ratio of time-dependent caval flows and the pulsatility content rather than mixing at the junction. "Failing" Fontan flow waveforms exhibited less balanced hepatic flow distribution to lungs.

CONCLUSIONS

The haemodynamic efficiency of single-ventricle circulation depends strongly on the pulsatility of venous flow waveforms. The proposed performance indices can be calculated easily in the clinical setting in efforts to better quantify the energy efficiency of Fontan venous waveforms in pulsatile settings.

Virtual Interventions for Image-based Blood Flow Computation

Image-based blood flow computation provides great promise for evaluation of vascular devices and assessment of surgical procedures. However, many previous studies employ idealized arterial and device models or only patient-specific models from the image data after device deployment, since the tools for model construction are unavailable or limited and tedious to use. Moreover, in contrast to retrospective studies from existing data, there is a pressing need for prospective analysis with the goal of surgical planning. Therefore, it is necessary to construct models with deployed devices in a fast, virtual and interactive fashion. The goal of this paper is to develop new geometric methods to deploy stents or stent grafts virtually to patient-specific geometric models constructed from a 3D segmentation of medical images. A triangular surface representing the vessel lumen boundary is extracted from the segmentation. The diseased portion is either clipped and replaced by the surface of a deployed device or rerouted in the case of a bypass graft. For diseased arteries close to bifurcations, bifurcated device models are generated. A method to map a 2D strut pattern on the surface of a device is also presented. We demonstrate three applications of our methods in personalized surgical planning for aortic aneurysms, aortic coarctation, and coronary artery stenosis using blood flow computation. Our approach enables prospective model construction and may help to expand the throughput required by routine clinical uses in the future.

An evaluation of dynamic outlet boundary conditions in a 1D fluid dynamics model

When modeling the cardiovascular system, the use of boundary conditions that closely represent the interaction between the region of interest and the surrounding vessels and organs will result in more accurate predictions. An often overlooked feature of outlet boundary conditions is the dynamics associated with regulation of the distribution of pressure and flow. This study implements a dynamic impedance outlet boundary condition in a one-dimensional fluid dynamics model using the pulmonary vasculature and respiration (feedback mechanism) as an example of a dynamic system. The dynamic boundary condition was successfully implemented and the pressure and flow were predicted for an entire respiration cycle. The cardiac cycles at maximal expiration and inspiration were predicted with a root mean square error of 0.61 and 0.59 mm Hg, respectively.

Visualization of flow structures in Fontan patients using 3-dimensional phase contrast magnetic resonance imaging

OBJECTIVE

Our objective was to analyze 3-dimensional (3D) blood flow patterns within the total cavopulmonary connection (TCPC) using in vivo phase contrast magnetic resonance imaging (PC MRI).

METHODS

Sixteen single-ventricle patients were prospectively recruited at 2 leading pediatric institutions for PC MRI evaluation of their Fontan pathway. Patients were divided into 2 groups. Group 1 comprised 8 patients with an extracardiac (EC) TCPC, and group 2 comprised 8 patients with a lateral tunnel (LT) TCPC. A coronal stack of 5 to 10 contiguous PC MRI slices with 3D velocity encoding (5-9 ms resolution) was acquired and a volumetric flow field was reconstructed.

RESULTS

Analysis revealed large vortices in LT TCPCs and helical flow structures in EC TCPCs. On average, there was no difference between LT and EC TCPCs in the proportion of inferior vena cava flow going to the left pulmonary artery (43% ± 7% vs 46% ± 5%; P = .34). However, for EC TCPCs, the presence of a caval offset was a primary determinant of inferior vena caval flow distribution to the pulmonary arteries with a significant bias to the offset side.

CONCLUSIONS

3D flow structures within LT and EC TCPCs were reconstructed and analyzed for the first time using PC MRI. TCPC flow patterns were shown to be different, not only on the basis of LT or EC considerations, but with significant influence from the superior vena cava connection as well. This work adds to the ongoing body of research demonstrating the impact of TCPC geometry on the overall hemodynamic profile.

Computer-Aided Patient-Specific Coronary Artery Graft Design Improvements Using CFD Coupled Shape Optimizer

This study aims to (i) demonstrate the efficacy of a new surgical planning framework for complex cardiovascular reconstructions, (ii) develop a computational fluid dynamics (CFD) coupled multi-dimensional shape optimization method to aid patient-specific coronary artery by-pass graft (CABG) design and, (iii) compare the hemodynamic efficiency of the sequential CABG, i.e., raising a daughter parallel branch from the parent CABG in patient-specific 3D settings. Hemodynamic efficiency of patient-specific complete revascularization scenarios for right coronary artery (RCA), left anterior descending artery (LAD), and left circumflex artery (LCX) bypasses were investigated in comparison to the stenosis condition. Multivariate 2D constraint optimization was applied on the left internal mammary artery (LIMA) graft, which was parameterized based on actual surgical settings extracted from 2D CT slices. The objective function was set to minimize the local variation of wall shear stress (WSS) and other hemodynamic indices (energy dissipation, flow deviation angle, average WSS, and vorticity) that correlate with performance of the graft and risk of re-stenosis at the anastomosis zone. Once the optimized 2D graft shape was obtained, it was translated to 3D using an in-house "sketch-based" interactive anatomical editing tool. The final graft design was evaluated using an experimentally validated second-order non-Newtonian CFD solver incorporating resistance based outlet boundary conditions. 3D patient-specific simulations for the healthy coronary anatomy produced realistic coronary flows. All revascularization techniques restored coronary perfusions to the healthy baseline. Multi-scale evaluation of the optimized LIMA graft enabled significant wall shear stress gradient (WSSG) relief (~34%). In comparison to original LIMA graft, sequential graft also lowered the WSSG by 15% proximal to LAD and diagonal bifurcation. The proposed sketch-based surgical planning paradigm evaluated the selected coronary bypass surgery procedures based on acute hemodynamic readjustments of aorta-CA flow. This methodology may provide a rational to aid surgical decision making in time-critical, patient-specific CA bypass operations before in vivo execution.

Numerical, hydraulic, and hemolytic evaluation of an intravascular axial flow blood pump to mechanically support Fontan patients

Currently available mechanical circulatory support systems are limited for adolescent and adult patients with a Fontan physiology. To address this growing need, we are developing a collapsible, percutaneously-inserted, axial flow blood pump to support the cavopulmonary circulation in Fontan patients. During the first phase of development, the design and experimental evaluation of an axial flow blood pump was performed. We completed numerical modeling of the pump using computational fluid dynamics analysis, hydraulic testing of a plastic pump prototype, and blood bag experiments (n=7) to measure the levels of hemolysis produced by the pump. Statistical analyses using regression were performed. The prototype with a 4-bladed impeller generated a pressure rise of 2-30 mmHg with a flow rate of 0.5-4 L/min for 3000-6000 RPM. A comparison of the experimental performance data to the numerical predictions demonstrated an excellent agreement with a maximum deviation being less than 6%. A linear increase in the plasma-free hemoglobin (pfHb) levels during the 6-h experiments was found, as desired. The maximum pfHb level was measured to be 21 mg/dL, and the average normalized index of hemolysis was determined to be 0.0097 g/100 L for all experiments. The hydraulic performance of the prototype and level of hemolysis are indicative of significant progress in the design of this blood pump. These results support the continued development of this intravascular pump as a bridge-to-transplant, bridge-to-recovery, bridge-to-hemodynamic stability, or bridge-to-surgical reconstruction for Fontan patients.

Cavopulmonary assist for the univentricular Fontan circulation: von Kármán viscous impeller pump

OBJECTIVE

In a univentricular Fontan circulation, modest augmentation of existing cavopulmonary pressure head (2-5 mm Hg) would reduce systemic venous pressure, increase ventricular filling, and thus substantially improve circulatory status. An ideal means of providing mechanical cavopulmonary support does not exist. We hypothesized that a viscous impeller pump, based on the von Kármán viscous pump principle, is optimal for this role.

METHODS

A 3-dimensional computational model of the total cavopulmonary connection was created. The impeller was represented as a smooth 2-sided conical actuator disk with rotation in the vena caval axis. Flow was modeled under 3 conditions: (1) passive flow with no disc; (2) passive flow with a nonrotating disk, and (3) induced flow with disc rotation (0-5K rpm). Flow patterns and hydraulic performance were examined for each case. Hydraulic performance for a vaned impeller was assessed by measuring pressure increase and induced flow over 0 to 7K rpm in a laboratory mock loop.

RESULTS

A nonrotating actuator disc stabilized cavopulmonary flow, reducing power loss by 88%. Disk rotation (from baseline dynamic flow of 4.4 L/min) resulted in a pressure increase of 0.03 mm Hg. A further increase in pressure of 5 to 20 mm Hg and 0 to 5 L/min flow was obtained with a vaned impeller at 0 to 7K rpm in a laboratory mock loop.

CONCLUSIONS

A single viscous impeller pump stabilizes and augments cavopulmonary flow in 4 directions, in the desired pressure range, without venous pathway obstruction. A viscous impeller pump applies to the existing staged protocol as a temporary bridge-to-recovery or -transplant in established univentricular Fontan circulations and may enable compressed palliation of single ventricle without the need for intermediary surgical staging or use of a systemic-to-pulmonary arterial shunt.

Optimization of inflow waveform phase-difference for minimized total cavopulmonary power loss

The Fontan operation is a palliative surgical procedure performed on children, born with congenital heart defects that have yielded only a single functioning ventricle. The total cavo-pulmonary connection (TCPC) is a common variant of the Fontan procedure, where the superior vena cava (SVC) and inferior vena cava (IVC) are routed directly into the pulmonary arteries (PA). Due to the limited pumping energy available, optimized hemodynamics, in turn, minimized power loss, inside the TCPC pathway is required for the best optimal surgical outcomes. To complement ongoing efforts to optimize the anatomical geometric design of the surgical Fontan templates, here, we focused on the characterization of power loss changes due to the temporal variations in between SVC and IVC flow waveforms. An experimentally validated pulsatile computational fluid dynamics solver is used to quantify the effect of phase-shift between SVC and IVC inflow waveforms and amplitudes on internal energy dissipation. The unsteady hemodynamics of two standard idealized TCPC geometries are presented, incorporating patient-specific real-time PC-MRI flow waveforms of "functional" Fontan patients. The effects of respiration and pulsatility on the internal energy dissipation of the TCPC pathway are analyzed. Optimization of phase-shift between caval flows is shown to lead to lower energy dissipation up to 30% in these idealized models. For physiological patient-specific caval waveforms, the power loss is reduced significantly (up to 11%) by the optimization of all three major harmonics at the same mean pathway flow (3 L/min). Thus, the hemodynamic efficiency of single ventricle circuits is influenced strongly by the caval flow waveform quality, which is regulated through respiratory dependent physiological pathways. The proposed patient-specific waveform optimization protocol may potentially inspire new therapeutic applications to aid postoperative hemodynamics and improve the well being of the Fontan patients.