Hemodynamics in a Pediatric Ascending Aorta Using a Viscoelastic Pediatric Blood Model

Tailoring left ventricular assist device cannula implantation using coupled multi-scale multi-objective optimization

BACKGROUND

The frequent occurrence of thromboembolic cerebral events continues to limit the widespread implementation of Ventricular Assist Devices (VAD) despite continued advancements in VAD design and anti-coagulation treatments. Recent studies point to the optimal positioning of the outflow graft (OG) as a potential mitigator of post implantation thromboembolism.

OBJECTIVE

This study aims to examine the tailoring of the OG implantation orientation with the goal of minimizing the number of thrombi reaching the cerebral vessels by means of a formal shape optimization scheme incorporated into a multi-scale hemodynamics analysis.

METHODS

A 3-D patient-specific computational fluid dynamics model is loosely coupled in a two-way manner to a 0-D lumped parameter model of the peripheral circulation. A Lagrangian particle-tracking scheme models and tracks thrombi as non-interacting solid spheres. The loose coupling between CFD and LPM is integrated into a geometric shape optimization scheme which aims to optimize an objective function that targets a drop in cerebral embolization, and an overall reduction in particle residence times.

RESULTS

The results elucidate the importance of OG anastomosis orientation and placement particularly in the case that studied particle release from the OG, as a fivefold decrease in cerebral embolization was observed between the optimal and non-optimal implantations. Another case considered particle release from the ventricle and aortic root walls, in which optimal implantation was achieved with a shallow insertion angle. Particle release from all three origins was investigated in the third case, demonstrating that the optimal configurations were generally characterized by VAD flow directed along the central lumen of the aortic arch. Because optimal configurations depended on the anatomic origin of the thrombus, it is important to determine, in clinical studies, the most likely sites of thrombus formation in VAD patients.

Effect of Hematocrit and Elevated Beat Rate on the 12cc Penn State Pediatric Ventricular Assist Device

Congenital heart disease affects approximately 40,000 infants annually in the United States with 25% requiring invasive treatment. Due to limited number of donor hearts and treatment options available for children, pediatric ventricular assist devices (PVADs) are used as a bridge to transplant. The 12cc pneumatic Penn State PVAD is optimized to prevent platelet adhesion and thrombus formation at patient nominal conditions; however, children demonstrate variable blood hematocrit and elevated heart rates. Therefore, with pediatric patients exhibiting greater variability, particle image velocimetry is used to evaluate the PVAD with three non-Newtonian hematocrit blood analogs (20%, 40%, and 60%) and at two beat rates (75 and 120 bpm) to understand the device's performance. The flow fields demonstrate a strong inlet jet that transitions to a solid body rotation during diastole. During systole, the rotation dissipates and reorganizes into an outlet jet. This flow field is consistent across all hematocrits and beat rates but at a higher velocity magnitude during 120 bpm. There are also minor differences in flow field timing and surface washing due to hematocrit. Therefore, despite patient differences in hematocrit or required pumping output, thorough surface washing can be achieved in the PVAD by altering operating conditions, thus reducing platelet adhesion potential.

Next Generation Development of Hybrid Continuous Flow Pediatric Total Artificial Heart Technology: Design-Build-Test

To address the unmet clinical need for pediatric circulatory support, we are developing an operationally versatile, hybrid, continuous-flow, total artificial heart ("Dragon Heart"). This device integrates a magnetically levitated axial and centrifugal blood pump. Here, we utilized a validated axial flow pump, and we focused on the development of the centrifugal pump. A motor was integrated to drive the centrifugal pump, achieving 50% size reduction. The motor design was simulated by finite element analysis, and pump design improvement was attained by computational fluid dynamics. A prototype centrifugal pump was constructed from biocompatible 3D printed parts for the housing and machined metal parts for the drive system. Centrifugal prototype testing was conducted using water and then bovine blood. The fully combined device ( i.e. , axial pump nested inside of the centrifugal pump) was tested to ensure proper operation. We demonstrated the hydraulic performance of the two pumps operating in tandem, and we found that the centrifugal blood pump performance was not adversely impacted by the simultaneous operation of the axial blood pump. The current iteration of this design achieved a range of operation overlapping our target range. Future design iterations will further reduce size and incorporate complete and active magnetic levitation.

Preliminary in-silico analysis of vascular graft implantation configuration and surface modification

Vascular grafts are used to reconstruct congenital cardiac anomalies, redirect flow, and offer vascular access. Donor tissue, synthetic, or more recently tissue-engineered vascular grafts each carry limitations spanning compatibility, availability, durability and cost. Synthetic and tissue-engineered grafts offer the advantage of design optimization using in-silico or in-vitro modeling techniques. We focus on an in-silico parametric study to evaluate implantation configuration alternatives and surface finishing impact of a novel silicon-lined vascular graft. The model consists of a synthetic 3D-generic model of a graft connecting the internal carotid artery to the jugular vein. The flow is assumed unsteady, incompressible, and blood is modeled as a non-Newtonian fluid. A comparison of detached eddy turbulence and laminar modeling to determine the required accuracy needed found mild differences mainly dictated by the roughness level. The conduit walls are modeled as non-compliant and fixed. The shunt configurations considered, are straight and curved with varied surface roughness. Following a grid convergence study, two shunt configurations are analyzed to better understand flow distribution, peak shear locations, stagnation regions and eddy formation. The curved shunt was found to have lower peak and mean wall-shear stress, while resulting in lower flow power system and decreased power loss across the graft. The curved smooth surface shunt shows lower peak and mean wall-shear stress and lower power loss when compared to the straight shunt.

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.

Performance of a Hydrogel Coated Nitinol with Oligonucleotide-Modified Nanoparticles Within Turbulent Conditions of Blood-Contacting Devices

INTRODUCTION

Hydrogels offer a wide range of applications in the antithrombotic modification of biomedical devices. The functionalization of these hydrogels with potentially drug-laden nanoparticles in the context of deviceassociated turbulence is critically under-studied. Thus, the purpose of this study was to use a hydrogel-coating nitinol surface as a model to understand the functions of hydrogels and the capture of nanoparticles under clinically relevant flow conditions.

METHODS

Nitinol was coated by an oligonucleotide (ON) functionalized hydrogel. Nanoparticles were functionalized with complementary oligonucleotides (CONs). The capture of CONfunctionalized nanoparticles by the ON-functionalized hydrogel surfaces was studied under both static and dynamic attachment conditions. Fluorescent-labelling of nanoparticles was utilized to assess capture efficacy and resistance to removal by device-relevant flow conditions.

RESULTS

The specificity of the ON-CON bond was verified, exhibiting a dose-dependent attachment response. The hydrogel coating was resistant to stripping by flow, retaining >95% after exposure to one hour of turbulent flow. Attachment of nanoparticles to the hydrogel was higher in the static condition than under laminar flow (p < 0.01), but comparable to that of attachment under turbulent flow. Modified nitinol samples underwent one hour of flow treatment under both laminar and turbulent regimes and demonstrated decreased nanoparticle loss following static conjugation rather than turbulent conjugation (36.1% vs 53.8%, p < 0.05). There was no significant difference in nanoparticle functionalization by upstream injection between laminar and turbulent flow.

CONCLUSION

The results demonstrate promising potential of hydrogelfunctionalized nitinol for capturing nanoparticles using nucleic acid hybridization. The hydrogel structure and ONCON bond integrity both demonstrated a resistance to mechanical damage and loss of biomolecular functionalization by exposure to turbulence. Further investigation is warranted to highlight drug delivery and antithrombogenic modification applications of nanoparticle-functionalized hydrogels.

Computational Modeling of the Penn State Fontan Circulation Assist Device

To address the increasing number of failing Fontan patients, Penn State University and the Penn State Hershey Medical Center are developing a centrifugal blood pump for long-term mechanical support. Computational fluid dynamics (CFD) modeling of the Penn State Fontan Circulatory Assist Device (FCAD) was performed to understand hemodynamics within the pump and its potential for hemolysis and thrombosis. CFD velocity and pressure results were first validated against experimental data and found to be within the standard deviations of the velocities and within 5% of the pressures. Further simulations performed with a human blood model found that most of the fluid domain was subjected to low shear stress (<50 Pa), with areas of highest stress around the rotor blade tips that increased with pump flow rate and rotor speed (138-178 Pa). However, the stresses compared well to previous CFD studies of commercial blood pumps and remained mostly below common thresholds of hemolysis and platelet activation. Additionally, few regions of low shear rate were observed within the FCAD, signifying minimal potential for platelet adhesion. These results further emphasize the FCAD's potential that has been observed previously in experimental and animal studies.

Development of the Centrifugal Blood Pump for a Hybrid Continuous Flow Pediatric Total Artificial Heart: Model, Make, Measure

Clinically-available blood pumps and total artificial hearts for pediatric patients continue to lag well behind those developed for adults. We are developing a hybrid, continuous-flow, magnetically levitated, pediatric total artificial heart (TAH). The hybrid TAH design integrates both an axial and centrifugal blood pump within a single, compact housing. The centrifugal pump rotates around the separate axial pump domain, and both impellers rotate around a common central axis. Here, we concentrate our development effort on the centrifugal blood pump by performing computational fluid dynamics (CFD) analysis of the blood flow through the pump. We also conducted transient CFD analyses (quasi-steady and transient rotational sliding interfaces) to assess the pump's dynamic performance conditions. Through modeling, we estimated the pressure generation, scalar stress levels, and fluid forces exerted on the magnetically levitated impellers. To further the development of the centrifugal pump, we also built magnetically-supported prototypes and tested these in an in vitro hydraulic flow loop and via 4-h blood bag hemolytic studies (n = 6) using bovine blood. The magnetically levitated centrifugal prototype delivered 0–6.75 L/min at 0–182 mmHg for 2,750–4,250 RPM. Computations predicted lower pressure-flow performance results than measured by testing; axial and radial fluid forces were found to be <3 N, and mechanical power usage was predicted to be <5 Watts. Blood damage indices (power law weighted exposure time and scalar stress) were <2%. All data trends followed expectations for the centrifugal pump design. Six peaks in the pressure rise were observed in the quasi-steady and transient simulations, correlating to the blade passage frequency of the 6-bladed impeller. The average N.I.H value (n = 6) was determined to be 0.09 ± 0.02 g/100 L, which is higher than desired and must be addressed through design improvement. These data serve as a strong foundation to build upon in the next development phase, whereby we will integrate the axial flow pump component.

Parametric investigation of an injection-jet self-powered Fontan circulation

Approximately \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1/2500$$\end{document}1/2500 babies are born with only one functioning ventricle and the Fontan is the third and, ideally final staged palliative operation for these patients. This altered circulation is prone to failure with survival rates below \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$50\%$$\end{document}50% into adulthood. Chronically elevated inferior vena cava (IVC) pressure is implicated as one cause of the mortality and morbidity in this population. An injection jet shunt (IJS) drawing blood-flow directly from the aortic arch to significantly lower IVC pressure is proposed. A computer-generated 3D model of a 2–4 year old patient with a fenestrated Fontan and a cardiac output of 2.3 L/min was generated. The detailed 3D pulsatile hemodynamics are resolved in a zero-dimensional lumped parameter network tightly-coupled to a 3D computational fluid dynamics model accounting for non-Newtonian blood rheology and resolving turbulence using large eddy simulation. IVC pressure and systemic oxygen saturation were tracked for various IJS-assisted Fontan configurations, altering design parameters such as shunt and fenestration diameters and locations. A baseline “failing” Fontan with a 4 mm fenestration was tuned to have an elevated IVC pressure (+ 17.8 mmHg). Enlargement of the fenestration to 8 mm resulted in a 3 mmHg IVC pressure drop but an unacceptable reduction in systemic oxygen saturation below 80%. Addition of an IJS with a 2 mm nozzle and minor volume load to the ventricle improved the IVC pressure drop to 3.2 mmHg while increasing systemic oxygen saturation above 80%. The salutary effects of the IJS to effectively lower IVC pressure while retaining acceptable levels of oxygen saturation are successfully demonstrated.

Tunable Blood Shunt for Neonates With Complex Congenital Heart Defects

Despite advancements in procedures and patient care, mortality rates for neonatal recipients of the Norwood procedure, a palliation for single ventricle congenital malformations, remain high due to the use of a fixed-diameter blood shunt. In this study, a new geometrically tunable blood shunt was investigated to address limitations of the current treatment paradigm (e.g., Modified Blalock-Taussig Shunt) by allowing for controlled modulation of blood flow through the shunt to accommodate physiological changes due to the patient’s growth. First, mathematical and computational cardiovascular models were established to investigate the hemodynamic requirements of growing neonatal patients with shunts and to inform design criteria for shunt diameter changes. Then, two stages of prototyping were performed to design, build and test responsive hydrogel systems that facilitate tuning of the shunt diameter by adjusting the hydrogel’s degree of crosslinking. We examined two mechanisms to drive crosslinking: infusion of chemical crosslinking agents and near-UV photoinitiation. The growth model showed that 15–18% increases in shunt diameter were required to accommodate growing patients’ increasing blood flow; similarly, the computational models demonstrated that blood flow magnitudes were in agreement with previous reports. These target levels of diameter increases were achieved experimentally with model hydrogel systems. We also verified that the photocrosslinkable hydrogel, composed of methacrylated dextran, was contact-nonhemolytic. These results demonstrate proof-of-concept feasibility and reflect the first steps in the development of this novel blood shunt. A tunable shunt design offers a new methodology to rebalance blood flow in this vulnerable patient population during growth and development.

In-silico analysis of outflow graft implantation orientation and cerebral thromboembolism incidence for full LVAD support

We investigate tailoring cannula implantation angles of left ventricle assist devices (LVAD) to reduce cerebral embolism risk for full LVAD support. We resolve pulsatile hemodynamics with a multi-scale computational fluid dynamics model coupled to a Lagrangian scheme tracking 2-5 mm particles for three cannula implantations. Blood is modeled as non-Newtonian. Cerebral flow distribution is altered depending on anastomosis angle and comparison of means embolization rates between steady and unsteady flow models show that unsteady modeling is more accurate even in the full LVAD support case. Intermediate angle implantation yields lowest cerebral embolization incidence of 11% vs 29% for normal and 36% for shallow implantation.

The effects of non-Newtonian blood modeling and pulsatility on hemodynamics in the food and drug administration's benchmark nozzle model

BACKGROUND

Computational fluid dynamics (CFD) is an important tool for predicting cardiovascular device performance. The FDA developed a benchmark nozzle model in which experimental and CFD data were compared, however, the studies were limited by steady flows and Newtonian models.

OBJECTIVE

Newtonian and non-Newtonian blood models will be compared under steady and pulsatile flows to evaluate their influence on hemodynamics in the FDA nozzle.

METHODS

CFD simulations were validated against the FDA data for steady flow with a Newtonian model. Further simulations were performed using Newtonian and non-Newtonian models under both steady and pulsatile flows.

RESULTS

CFD results were within the experimental standard deviations at nearly all locations and Reynolds numbers. The model differences were most evident at Re = 500, in the recirculation regions, and during diastole. The non-Newtonian model predicted blunter upstream velocity profiles, higher velocities in the throat, and differences in the recirculation flow patterns. The non-Newtonian model also predicted a greater pressure drop at Re = 500 with minimal differences observed at higher Reynolds numbers.

CONCLUSIONS

An improved modeling framework and validation procedure were used to further investigate hemodynamics in geometries relevant to cardiovascular devices and found that accounting for blood's non-Newtonian and pulsatile behavior can lead to large differences in predictions in hemodynamic parameters.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Computational Fluid Dynamics Study of Cerebral Thromboembolism Risk in Ventricular Assist Device Patients: Effects of Pulsatility and Thrombus Origin

This study investigates the hypothesis that by surgically manipulating the outflow graft (OG) implantation during ventricle assist device placement, it may be possible to reduce the risk of cerebral embolism. We investigate this hypothesis using a computational approach on a patient-specific basis under fully pulsatile hemodynamics with a multiscale computational fluid dynamics model incorporating a coupled Eulerian-Lagrangian scheme that effectively tracks emboli in the fluid domain. Blood is modeled as a non-Newtonian fluid based on the hematocrit level. Preliminary flow analysis shows that depending on the anastomosis angle the left ventricular assist device (LVAD) can enhance the flow to the cerebral circulation by nearly 31%. Z-test results suggest that unsteady-flow modeling ought to be an integral part of any cardiovascular simulation with residual ventricular function. Assuming unsteady-flow conditions, a shallow LVAD outflow graft anastomosis angle is the most optimal if thrombi are released from the aortic-root reducing cerebral embolization incidence to 15.5% and from the ventricle to 17%, while a more pronounced anastomosis angle becomes advantageous when particles originate from the LVAD with an embolization rate of 16.9%. Overall, computations suggest that a pronounced LVAD anastomosis angle is the better implementation. Unsteady modeling is shown to be necessary for the presence of significant antegrade aortic-root flow which induces cyclical flow patterns due to residual pulsatility. On the other hand, depending on thrombus origin and ventricular assist devices (VAD) anastomosis angle there is a strong tradeoff in embolization rates.

Role of the left coronary artery geometry configuration in atherosusceptibility: CFD simulations considering sPTT model for blood

The achievement of clinically viable methodologies to simulate the hemodynamics in patient-specific coronary arteries is still a major challenge. Therefore, the novelty of this work is attained by the introduction of the viscoelastic property of blood in the numerical simulations, to study the role of the left coronary artery (LCA) geometry configuration in the atherosusceptibility. Apparently healthy patients were used and four different methodologies were tested. The methodology giving the most accurate results at the same time of having the lowest computational time is the one considering the viscoelastic property of blood and computational fluid dynamics. A Pearson correlation analysis was used to highlight relationships between geometric configuration and hemodynamic descriptors based on the simulated wall shear stress (WSS). The left main stem (LMS) has the greatest atherosusceptibility followed by the left anterior descending artery (LAD) since the relative residence time (RRT) average values are 3.81 and 3.70 Pa^-1, respectively. The geometric parameters with relevant contribution to directional flow change are the cross-sectional areas, especially the one of LMS segment (A_LMS), and the curvature of LMS segment. For LMS and LAD segments, when A_LMS increases, blood flow disturbance (r = 0.81 in LMS and r = 0.74 in LAD) and atherosusceptibility (r = 0.84 in LMS and r = 0.85 in LAD) increases. When the LMS curvature decreases, the WSS magnitude (r = 0.80 in LMS and r = 0.83 in LAD) decreases, and disturbance (r=-0.80 in LMS and r=-0.91 in LAD) and atherosusceptibility (r=-0.74 in LMS and r=-0.74 in LAD) increases.

Computational fluid dynamics: a primer for congenital heart disease clinicians

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

Non-Newtonian Effects on Patient-Specific Modeling of Fontan Hemodynamics

The Fontan procedure is a common palliative surgery for congenital single ventricle patients. In silico and in vitro patient-specific modeling approaches are widely utilized to investigate potential improvements of Fontan hemodynamics that are related to long-term complications. However, there is a lack of consensus regarding the use of non-Newtonian rheology, warranting a systematic investigation. This study conducted in silico patient-specific modeling for twelve Fontan patients, using a Newtonian and a non-Newtonian model for each patient. Differences were quantified by examining clinically relevant metrics: indexed power loss (iPL), indexed viscous dissipation rate (iVDR), hepatic flow distribution (HFD), and regions of low wall shear stress (A_WSS). Four sets of "non-Newtonian importance factors" were calculated to explore their effectiveness in identifying the non-Newtonian effect. No statistical differences were observed in iPL, iVDR, and HFD between the two models at the population-level, but large inter-patient variations exist. Significant differences were detected regarding A_WSS, and its correlations with non-Newtonian importance factors were discussed. Additionally, simulations using the non-Newtonian model were computationally faster than those using the Newtonian model. These findings distinguish good importance factors for identifying non-Newtonian rheology and encourage the use of a non-Newtonian model to assess Fontan hemodynamics.

Computational modeling of the Food and Drug Administration's benchmark centrifugal blood pump

In order to simulate hemodynamics within centrifugal blood pumps and to predict pump hemolysis, CFD simulations must be thoroughly validated against experimental data. They must also account for and accurately model the specific working fluid in the pump, whether that is a blood-analog solution to match an experimental PIV study or animal blood in a hemolysis experiment. Therefore, the Food and Drug Administration (FDA) benchmark centrifugal blood pump and its database of experimental PIV and hemolysis data were used to thoroughly validate CFD simulations of the same blood pump. A Newtonian blood model was first used to compare to the PIV data with a blood analog fluid while hemolysis data were compared using a power-law hemolysis model fit to porcine blood data. A viscoelastic blood model was then incorporated into the CFD solver to investigate the importance of modeling blood's viscoelasticity in centrifugal pumps. The established computational framework, including a dynamic rotating mesh, animal blood-specific fluid properties and hemolysis modeling, and a k-ω SST turbulence model, was shown to more accurately predict pump pressure heads, velocity fields, and hemolysis compared to previously published CFD studies of the FDA centrifugal pump. The CFD simulations were able to match the FDA pressure and hemolysis data for multiple pump operating conditions, with the CFD results being within the standard deviations of the experimental results. While CFD radial velocity profiles between the impeller blades also compared well to the PIV velocity results, more work is still needed to address the large variability among both experimental and computational predictions of velocity in the diffuser outlet jet. Small differences were observed between the Newtonian and viscoelastic blood models in pressure head and hemolysis at the higher flow rate cases (FDA Conditions 4 and 5) but were more significant at lower flow rate and pump impeller speeds (FDA Condition 1). These results suggest that the importance of accounting for blood's viscoelasticity may be dependent on the specific blood pump operating conditions. This detailed computational framework with improved modeling techniques and an extensive validation procedure will be used in future CFD studies of centrifugal blood pumps to aid in device design and predictions of their biological responses.

Accounting for residence-time in blood rheology models: do we really need non-Newtonian blood flow modelling in large arteries?

Patient-specific computational fluid dynamics (CFD) is a promising tool that provides highly resolved haemodynamics information. The choice of blood rheology is an assumption in CFD models that has been subject to extensive debate. Blood is known to exhibit shear-thinning behaviour, and non-Newtonian modelling has been recommended for aneurysmal flows. Current non-Newtonian models ignore rouleaux formation, which is the key player in blood's shear-thinning behaviour. Experimental data suggest that red blood cell aggregation and rouleaux formation require notable red blood cell residence-time (RT) in a low shear rate regime. This study proposes a novel hybrid Newtonian and non-Newtonian rheology model where the shear-thinning behaviour is activated in high RT regions based on experimental data. Image-based abdominal aortic and cerebral aneurysm models are considered and highly resolved CFD simulations are performed using a minimally dissipative solver. Lagrangian particle tracking is used to define a backward particle RT measure and detect stagnant regions with increased rouleaux formation likelihood. Our novel RT-based non-Newtonian model shows a significant reduction in shear-thinning effects and provides haemodynamic results qualitatively identical and quantitatively close to the Newtonian model. Our results have important implications in patient-specific CFD modelling and suggest that non-Newtonian models should be revisited in large artery flows.

Patient-specific multiscale computational fluid dynamics assessment of embolization rates in the hybrid Norwood: effects of size and placement of the reverse Blalock-Taussig shunt

The hybrid Norwood operation is performed to treat hypoplastic left heart syndrome. Distal arch obstruction may compromise flow to the brain. In a variant of this procedure, a synthetic graft (reverse Blalock-Taussig shunt) is placed between the pulmonary trunk and innominate artery to improve upper torso blood flow. Thrombi originating in the graft may embolize to the brain. In this study, we used computational fluid dynamics and particle tracking to investigate the patterns of particle embolization as a function of the anatomic position of the reverse Blalock-Taussig shunt. The degree of distal arch obstruction and position of particle origin influence embolization probabilities to the cerebral arteries. Cerebral embolization probabilities can be reduced by as much as 20% by optimizing graft position, for a given arch geometry, degree of distal arch obstruction, and particle origin. There is a tradeoff, however, between cerebral pulmonary and coronary embolization probabilities.

Asynchronous Pumping of a Pulsatile Ventricular Assist Device in a Pediatric Anastomosis Model

BACKGROUND

Both pulsatile and continuous flow ventricular assist devices are being developed for pediatric congenital heart defect patients. Pulsatile devices are often operated asynchronously with the heart in either an "automatic" or a fixed beat rate mode. However, most studies have only investigated synchronized ejection.

METHODS

A previously validated viscoelastic blood solver is used to investigate the parameters of pulsatility, power loss, and graft failure in a pediatric aortic anastomosis model.

RESULTS

Pulsatility was highest with synchronized flow and lowest at a 90° phase shift. Power loss decreased at 90° and 180° phase shifts but increased at a 270° phase shift. Similar regions of potential intimal hyperplasia and graft failure were seen in all cases but with phase-shifted ejection leading to higher wall shear stress on the anastomotic floor and oscillatory shear index on the anastomotic toe.

CONCLUSION

The ranges of pulsatility and hemodynamics that can result clinically using asynchronous pulsatile devices were investigated in a pediatric anastomosis model. These results, along with the different postoperative benefits of pump modulation, can be used to design an optimal weaning protocol.

Modeling Outcomes: Modified Aortic Arch Advancement for Neonatal Hypoplastic Arch

OBJECTIVE

Numerous surgical approaches regarding aortic arch advancement for neonatal arch hypoplasia have been described. These repairs can be classified into two categories: those that incorporate a patch and those that do not. The decision between repairs remains largely experiential, rather than empirical, because of the limited number of reported outcomes. We report early outcomes from neonates undergoing modified aortic arch advancement with an anterior patch and our experience using computational fluid dynamic modeling to better understand the hemodynamic consequences associated with this repair.

METHODS

A retrospective review of neonates undergoing aortic arch advancement with anterior patch in 2014 at a single institution was performed. Anatomical, perioperative, and follow-up data were collected. Three-dimensional cardiac magnetic resonance images were used to generate computational fluid dynamic models of the modified anterior patch and direct end-to-side repairs. Cardiac waveform inputs were simulated and hemodynamic analyzed.

RESULTS

Ten neonates underwent modified aortic arch advancement. No hemodynamically significant gradients were observed at a median follow-up of 0.77 (0.30-1.2) years. Asymmetrical flow was observed in the end-to-side repair, whereas more concentric laminar flow was observed throughout the modified model. Spatial variations in velocities immediately distal to the anastomosis were greater in the end-to-side model (0.35 vs 0.17 m/s, P < 0.001). Time-averaged variations in wall shear stress during systole were greater in the end-to-side model at the same location (3.44 vs 1.98 dynes/cm, P < 0.001).

CONCLUSIONS

Early outcomes after the use of an anterior patch for neonatal hypoplastic aortic arch repair show favorable hemodynamic outcomes.

Comparison of Blood Viscoelasticity in Pediatric and Adult Cardiac Patients

Evidence is accumulating that blood flow patterns in the cardiovascular system and in cardiovascular devices do, in some instances, depend on blood viscoelasticity. Thus, to better understand the challenges to providing circulatory support and surgical therapies for pediatric and adult patients, viscous and elastic components of complex blood viscoelasticity of 31 pediatric patients were compared to those of 29 adult patients with a Vilastic-3 rheometer. A random effects model with categorical age covariates found statistically significant differences between pediatric and adult patients for log viscosity (p = 0.005). Log strain (p < 0.0001) and hematocrit (p < 0.0001) effects were also significant, as were the hematocrit-by-log-strain (p = 0.0006) and age-by-log strain (p = 0.001) interactions. The hematocrit-by-age interaction was not significant. For log elasticity, age differences were insignificant (p = 0.39). The model for log elasticity had significant log strain (p < 0.0001), log strain squared (p < 0.0001) and hematocrit (p < 0.0001) effects, as well as hematocrit-by-log-strain and hematocrit-by-log-strain-squared interactions (p = 0.014). A model for log viscosity with continuous age was also fit to the data, which can be used to refine cardiovascular device design and operation to the age of the patient. We conclude that there are distinct differences between pediatric and adult blood viscosity, as well as substantial variation within the pediatric population, that may impact the performance of devices and procedures.

Development of a platelet adhesion transport equation for a computational thrombosis model

Thrombosis is a significant issue for cardiovascular device development and use. While thrombosis models are available, very few are device-related and none have been thoroughly validated experimentally. Here, we introduce a surface adherent platelet transport equation into a continuum model to account for the biomaterial interface/blood interaction. Using a rotating disc system and polyurethane-urea material, we characterize steady and pulsatile flow fields using laser Doppler velocimetry. In vitro measurements of platelet adhesion are used in combination with the LDV data to provide further experimental validation. The rotating disc system is computationally studied using the device-induced thrombosis model with the surface platelet adherent transport equation. The results indicate that the flow field is in excellent agreement to the experimental LDV data and that the platelet adhesion simulations are in good agreement with the in vitro platelet data. These results provide good evidence that this transport equation can be used to express the relationship between blood and a biomaterial if the correct platelet adhesion characteristics are known for the biomaterial. Further validation is necessary with other materials.

Computational modelling for congenital heart disease: how far are we from clinical translation?

Computational models of congenital heart disease (CHD) have become increasingly sophisticated over the last 20 years. They can provide an insight into complex flow phenomena, allow for testing devices into patient-specific anatomies (pre-CHD or post-CHD repair) and generate predictive data. This has been applied to different CHD scenarios, including patients with single ventricle, tetralogy of Fallot, aortic coarctation and transposition of the great arteries. Patient-specific simulations have been shown to be informative for preprocedural planning in complex cases, allowing for virtual stent deployment. Novel techniques such as statistical shape modelling can further aid in the morphological assessment of CHD, risk stratification of patients and possible identification of new ‘shape biomarkers’. Cardiovascular statistical shape models can provide valuable insights into phenomena such as ventricular growth in tetralogy of Fallot, or morphological aortic arch differences in repaired coarctation. In a constant move towards more realistic simulations, models can also account for multiscale phenomena (eg, thrombus formation) and importantly include measures of uncertainty (ie, CIs around simulation results). While their potential to aid understanding of CHD, surgical/procedural decision-making and personalisation of treatments is undeniable, important elements are still lacking prior to clinical translation of computational models in the field of CHD, that is, large validation studies, cost-effectiveness evaluation and establishing possible improvements in patient outcomes.

Continuous and Pulsatile Pediatric Ventricular Assist Device Hemodynamics with a Viscoelastic Blood Model

To investigate the effects of pulsatile and continuous pediatric ventricular assist (PVAD) flow and pediatric blood viscoelasticity on hemodynamics in a pediatric aortic graft model. Hemodynamic parameters of pulsatility, along with velocity and wall shear stress (WSS), are analyzed and compared between Newtonian and viscoelastic blood models at a range of physiological pediatric hematocrits using computational fluid dynamics. Both pulsatile and continuous PVAD flow lead to a decrease in pulsatility (surplus hemodynamic energy, ergs/cm(3)) compared to healthy aortic flow but with continuous PVAD pulsatility up to 2.4 times lower than pulsatile PVAD pulsatility at each aortic outlet. Significant differences are also seen between the two flow modes in velocity and WSS. The higher velocity jet during systole with pulsatile flow leads to higher WSSs at the anastomotic toe and at the aortic branch bifurcations. The lower velocity but continuous flow jet leads to a much different flow field and higher WSSs into diastole. Under a range of physiological pediatric hematocrit (20-60%), both velocity and WSS can vary significantly with the higher hematocrit blood model generally leading to higher peak WSSs but also lower WSSs in regions of flow separation. The large decrease in pulsatility seen from continuous PVAD flow could lead to complications in pediatric vascular development while the high WSSs during peak systole from pulsatile PVAD flow could lead to blood damage. Both flow modes lead to similar regions prone to intimal hyperplasia resulting from low time-averaged WSS and high oscillatory shear index.