Characterization of neonatal aortic cannula jet flow regimes for improved cardiopulmonary bypass

Impact of a Multidisciplinary Research Team Approach to Prevent Avoidable Mistakes for Neonatal CPB Population

Significant progress has been made in neonatal cardiopulmonary bypass (CPB) over the years. At Penn State Health Children's Hospital, we have established a multidisciplinary research team that brings clinicians, engineers, scientists, research nurses, neuromonitoring technicians, perfusionists, and students from various departments to help reduce adverse outcomes following CPB in neonates. With the help of this team, we evaluate each CPB component in simulated conditions identical to those used in clinical practice. The objective of this review is to demonstrate the results of these translational projects and present critical mistakes to avoid for neonatal CPB patients.

Manual Carotid Compression is a Viable Alternative for Reduction of Cerebral Microemboli

BACKGROUND

Stroke is a devastating complication of cardiovascular surgeries, and the risk is particularly high for those requiring cardiopulmonary bypass (CPB). Embolic particles generated during the unclamping of the aortic cross-clamp may enter the cerebral circulation, lodging in small vessels. External manual compression of the carotid arteries is a non-invasive technique that has been proposed for cerebral protection during CPB procedures but is not widely deployed.

METHODS

The aim of this study is to assess the potential for cerebral emboli reduction with carotid compression using an in vitro model. Experiments were performed with a glass aortic arch model in a mock cardiovascular circuit. Small fluorescent particles were released into the circulation with and without carotid compression, and the particles visualized in the aortic midplane. The number of particles in the aorta and arch branch vessels were counted from the images before, during and following the release of carotid compression for durations of 10, 15 and 20 s. A gamma variate function was fit to the data to describe the bolus dynamics.

RESULTS

Carotid compression for 10 s reduces the number of embolic articles entering the carotid arteries by over 75%. A compression duration of 15-20 s does not result in greater particle reduction than one of 10 s.

CONCLUSION

Brief compression of the common carotid arteries during cardiovascular interventions has the potential to dramatically reduce the number of cerebral emboli and should be investigated further.

Flow Dynamics and Mixing in Extracorporeal Support: A Study of the Return Cannula

Cannulation strategies in medical treatment such as in extracorporeal life support along with the associated cannula position, orientation and design, affects the mixing and the mechanical shear stress appearing in the flow field. This in turn influences platelet activation state and blood cell destruction. In this study, a co-flowing confined jet similar to a return cannula flow configuration found in extracorporeal membrane oxygenation was investigated experimentally. Cannula diameters, flow rate ratios between the jet and the co-flow and cannula position were studied using Particle Image Velocimetry and Planar Laser Induced Fluorescence. The jet was turbulent for all but two cases, in which a transitional regime was observed. The mixing, governed by flow entrainment, shear layer induced vortices and a backflow along the vessel wall, was found to require 9–12 cannula diameters to reach a fully homogeneous mixture. This can be compared to the 22–30 cannula diameters needed to obtain a fully developed flow. Although not significantly affecting mixing characteristics, cannula position altered the development of the flow structures, and hence the shear stress characteristics.

Aortic cannula orientation and flow impacts embolic trajectories: computational cardiopulmonary bypass

INTRODUCTION

Emboli events are associated with the aortic cannula insertion and final position in the ascending aorta. However, the impact of subtle changes in aortic cannula movement and flow influencing embolic transport throughout the aortic arch is not well understood. The present study evaluated the aortic cannula's outflow and orientation effect on emboli entering the aortic branch arteries.

METHODS

A simplified aortic computational model was anteriorly cannulated in the distal ascending aorta with a 21-French straight aortic cannula, and two orientations were analysed by injecting gaseous and solid emboli at pump flows 2, 3 and 5 L/minute. The first aortic cannula orientation (forward flow cannula) was directed towards the lesser curvature. The second aortic cannula orientation (rear flow cannula) was tilted slightly backwards by 15°, providing flow in the retrograde direction.

RESULTS

Forward flow cannula produced a primary arch flow, whereas rear flow cannula produced a secondary arch flow resulting in four times longer emboli arch resident times than forward flow cannula. The rear flow cannula had the highest percentage of gaseous emboli entering the brachiocephalic artery of 8%, 12% and 36% (at 2, 3 and 5 L/minute, respectively). Rear flow cannula provided a positive aortic branch arterial flow at all pump flows, whereas at forward flow cannula, the brachiocephalic artery experienced retrograde flows of -1.0% (3 L/minute) and -4.0% (5 L/minute), with the left common carotid -0.23% (5 L/minute). No significant number of solid emboli entered the aortic branch arteries.

CONCLUSION

This numerical study illustrated distinct trajectory behaviours between gaseous and solid emboli where slight changes in aortic cannula orientation influenced idealised emboli direction with higher pump flows magnifying the effects.

Hemodynamics of neonatal double lumen cannula malposition

Objective:

Malposition of dual lumen cannula is a frequent and challenging complication in neonates and plays a significant role in shaping the in vitro device hemodynamics. This study aims to analyze the effect of the dual lumen cannula malposition on right-atrial hemodynamics in neonatal patients using an experimentally validated computational fluid dynamics model.

Methods:

A computer model was developed for clinically approved dual lumen cannula (13Fr Origen Biomedical, Austin, Texas, USA) oriented inside the atrium of a 3-kg neonate with normal venous return. Atrial hemodynamics and dual lumen cannula malposition were systematically simulated for two rotations (antero-atrial and atrio-septal) and four translations (two intravascular movements along inferior vena cava and two dislodged configurations in the atrium). A multi-domain compartmentalized mesh was prepared to allow the site-specific evaluation of important hemodynamic parameters. Transport of each blood stream, blood damage levels, and recirculation times are quantified and compared to dual lumen cannula in proper position.

Results:

High recirculation levels (39 ± 4%) in malpositioned cases resulted in poor oxygen saturation where maximum recirculation of up to 42% was observed. Apparently, Origen dual lumen cannula showed poor inferior vena cava blood–capturing efficiency (48 ± 8%) but high superior vena cava blood–capturing efficiency (86 ± 10%). Dual lumen cannula malposition resulted in corresponding changes in residence time (1.7 ± 0.5 seconds through the tricuspid). No significant differences in blood damage were observed among the simulated cases compared to normal orientation. Compared to the correct dual lumen cannula position, both rotational and translational displacements of the dual lumen cannula resulted in significant hemodynamic differences.

Conclusion:

Rotational or translational movement of dual lumen cannula is the determining factor for atrial hemodynamics, venous capturing efficiency, blood residence time, and oxygenated blood delivery. Results obtained through computational fluid dynamics methodology can provide valuable foresight in assessing the performance of the dual lumen cannula in patient-specific configurations.

Hemodynamic effects of perfusion level of peripheral ECMO on cardiovascular system

Background

Peripheral ECMO is an effective cardiopulmonary support in clinical. The perfusion level could directly influence the performances and complications. However, there are few studies on the effects of the perfusion level on hemodynamics of peripheral ECMO.

Methods

The geometric model of cardiovascular system with peripheral ECMO was established. The blood assist index was used to classify the perfusion level of the ECMO. The flow pattern from the aorta to the femoral artery and their branches, blood flow rate from aorta to brain and limbs, flow interface, harmonic index of blood flow, wall shear stress and oscillatory shear index were chosen to evaluate the hemodynamic effects of peripheral ECMO.

Results

The results demonstrated that the flow rate of aorta outlets increased and perfusion condition had been improved. And the average flow to the upper limbs and brain has a positive correlation with BAI (r = 0.037, p < 0.05), while there is a negative correlation with lower limbs (r = − 0.054, p < 0.05). The HI has negative correlation with BAI (p < 0.05, r < 0). The blood interface is further from the heart with the BAI decrease. And the average WSS has negative correlation with BAI (p < 0.05, r = − 0.983) at the bifurcation of femoral aorta and has positive correlation with BAI (p < 0.05, r = 0.99) at the inner aorta. The OSI under different BAI is higher (reaching 0.4) at the inner wall of the aortic arch, the descending aorta and the femoral access.

Conclusions

The pathogenesis of peripheral ECMO with different perfusion levels varies; its further research will be thorough and extensive.

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

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

Mytilus galloprovincialis as a smart micro-pump

Hydrodynamic performance of the marine mussel, Mytilus galloprovincialis, is studied with time-resolved particle image velocimetry. We evaluated inhalant flow, exhalant jet flow, suction performance and flow control capabilities of the mussels quantitatively. Inhalant flow structures of mussels are measured at the coronal plane for the first time in literature. Nutrient fluid is convected into the mussel by three-dimensional sink flow. Inhalant velocity reaches its highest magnitude inside the mussel mantle while it is accelerating outward from the mussels. We calculated pressure gradient at the coronal plane. As inhalant flow approaches the mussel shell tip, suction force generated by the inhalant flow increases and becomes significant at the shell tip. Likewise, exhalant jet flow regimes were studied for 17 mussels. Mussels can control their exhalant jet flow structure from a single potential core region to double potential core region or vice versa. Peak exhalant jet velocity generated by the mussels changes between 2.77 cm s-1 and 11.1 cm s-1 as a function of mussel cavity volume. Measurements of hydrodynamic dissipation at the sagittal plane revealed no interaction between the inhalant and exhalant jet flow, indicating energy-efficient synchronized pumping mechanism. This efficient pumping mechanism is associated with the flow-turning angle between inhalant and exhalant jet flows, ∼90° (s.d. 12°).

Endoleak Assessment Using Computational Fluid Dynamics and Image Processing Methods in Stented Abdominal Aortic Aneurysm Models

Endovascular aortic aneurysm repair (EVAR) is a predominant surgical procedure to reduce the risk of aneurysm rupture in abdominal aortic aneurysm (AAA) patients. Endoleak formation, which eventually requires additional surgical reoperation, is a major EVAR complication. Understanding the etiology and evolution of endoleak from the hemodynamic perspective is crucial to advancing the current posttreatments for AAA patients who underwent EVAR. Therefore, a comprehensive flow assessment was performed to investigate the relationship between endoleak and its surrounding pathological flow fields through computational fluid dynamics and image processing. Six patient-specific models were reconstructed, and the associated hemodynamics in these models was quantified three-dimensionally to calculate wall stress. To provide a high degree of clinical relevance, the mechanical stress distribution calculated from the models was compared with the endoleak positions identified from the computed tomography images of patients through a series of imaging processing methods. An endoleak possibly forms in a location with high local wall stress. An improved stent graft (SG) structure is conceived accordingly by increasing the mechanical strength of the SG at peak wall stress locations. The presented analytical paradigm, as well as numerical analysis using patient-specific models, may be extended to other common human cardiovascular surgeries.

cMRI-BED: A novel informatics framework for cardiac MRI biomarker extraction and discovery applied to pediatric cardiomyopathy classification

Background

Pediatric cardiomyopathies are a rare, yet heterogeneous group of pathologies of the myocardium that are routinely examined clinically using Cardiovascular Magnetic Resonance Imaging (cMRI). This gold standard powerful non-invasive tool yields high resolution temporal images that characterize myocardial tissue. The complexities associated with the annotation of images and extraction of markers, necessitate the development of efficient workflows to acquire, manage and transform this data into actionable knowledge for patient care to reduce mortality and morbidity.

Methods

We develop and test a novel informatics framework called cMRI-BED for biomarker extraction and discovery from such complex pediatric cMRI data that includes the use of a suite of tools for image processing, marker extraction and predictive modeling. We applied our workflow to obtain and analyze a dataset of 83 de-identified cases and controls containing cMRI-derived biomarkers for classifying positive versus negative findings of cardiomyopathy in children. Bayesian rule learning (BRL) methods were applied to derive understandable models in the form of propositional rules with posterior probabilities pertaining to their validity. Popular machine learning methods in the WEKA data mining toolkit were applied using default parameters to assess cross-validation performance of this dataset using accuracy and percentage area under ROC curve (AUC) measures.

Results

The best 10-fold cross validation predictive performance obtained on this cMRI-derived biomarker dataset was 80.72% accuracy and 79.6% AUC by a BRL decision tree model, which is promising from this type of rare data. Moreover, we were able to verify that mycocardial delayed enhancement (MDE) status, which is known to be an important qualitative factor in the classification of cardiomyopathies, is picked up by our rule models as an important variable for prediction.

Conclusions

Preliminary results show the feasibility of our framework for processing such data while also yielding actionable predictive classification rules that can augment knowledge conveyed in cardiac radiology outcome reports. Interactions between MDE status and other cMRI parameters that are depicted in our rules warrant further investigation and validation. Predictive rules learned from cMRI data to classify positive and negative findings of cardiomyopathy can enhance scientific understanding of the underlying interactions among imaging-derived parameters.

Growth and hemodynamics after early embryonic aortic arch occlusion

The majority of severe clinically significant forms of congenital heart disease (CHD) are associated with great artery lesions, including hypoplastic, double, right or interrupted aortic arch morphologies. While fetal and neonatal interventions are advancing, their potential ability to restore cardiac function, optimal timing, location, and intensity required for intervention remain largely unknown. Here, we combine computational fluid dynamics (CFD) simulations with in vivo experiments to test how individual pharyngeal arch artery hemodynamics alter as a result of local interventions obstructing individual arch artery flow. Simulated isolated occlusions within each pharyngeal arch artery were created with image-derived three-dimensional (3D) reconstructions of normal chick pharyngeal arch anatomy at Hamburger-Hamilton (HH) developmental stages HH18 and HH24. Acute flow redistributions were then computed using in vivo measured subject-specific aortic sinus inflow velocity profiles. A kinematic vascular growth-rendering algorithm was then developed and implemented to test the role of changing local wall shear stress patterns in downstream 3D morphogenesis of arch arteries. CFD simulations predicted that altered pressure gradients and flow redistributions were most sensitive to occlusion of the IVth arches. To evaluate these simulations experimentally, a novel in vivo experimental model of pharyngeal arch occlusion was developed and implemented using two-photon microscopy-guided femtosecond laser-based photodisruption surgery. The right IVth arch was occluded at HH18, and resulting diameter changes were followed for up to 24 h. Pharyngeal arch diameter responses to acute hemodynamic changes were predicted qualitatively but poorly quantitatively. Chronic growth and adaptation to hemodynamic changes, however, were predicted in a subset of arches. Our findings suggest that this complex biodynamic process is governed through more complex forms of mechanobiological vascular growth rules. Other factors in addition to wall shear stress or more complex WSS rules are likely important in the long-term arterial growth and patterning. Combination in silico/experimental platforms are essential for accelerating our understanding and prediction of consequences from embryonic/fetal cardiovascular occlusions and lay the foundation for noninvasive methods to guide CHD diagnosis and fetal intervention.

Novel Straight Type Aortic Cannula with Spiral Flow Inducing Design

Aortic cannula is one of major factors leading to adverse events such as thrombosis and atherosclerosis development during open heart surgery. This is due to oxygenated blood outflow with high velocity jet from heart lung machine when exiting the cannula tip into ascending aorta. It was discovered, and validated by several researchers that blood flow out of the left ventricle into the aorta is spiral in nature. In this study, a novel design in which internal profile of the cannula was made to induce spiral flow were tested by way of numerical simulation, and compared against existing commercial cannula. Three designs were tested, which differed in number of groove employed. Among the cannula model designs, cannula design with 4 grooves yielded the lowest value of maximum wall shear stress at testing tube with 3.778 Pa and highest value of area weighted helicity density at 40 mm from cannula tips with 11.829 m/s2. Overall, spiral cannula models were showed highly potential in inducing spiral flow, and also the effect on blood hemolysis is acceptable.

Hydrodynamic evaluation of aortic cardiopulmonary bypass cannulae using particle image velocimetry

The high velocity jet from aortic arterial cannulae used during cardiopulmonary bypass potentially causes a "sandblasting" injury to the aorta, increasing the possibility of embolisation of atheromatous plaque. We investigated a range of commonly available dispersion and non-dispersion cannulae, using particle image velocimetry. The maximum velocity of the exit jet was assessed 20 and 40 mm from the cannula tip at flow rates of 3 and 5 L/min. The dispersion cannulae had lower maximum velocities compared to the non-dispersion cannulae. Dispersion cannulae had fan-shaped exit profiles and maximum velocities ranged from 0.63 to 1.52 m/s when measured at 20 mm and 5 L/min. Non-dispersion cannulae had maximum velocities ranging from 1.52 to 3.06 m/s at 20 mm and 5 L/min, with corresponding narrow velocity profiles. This study highlights the importance of understanding the hydrodynamic performance of these cannulae as it may help in selecting the most appropriate cannula to minimize the risk of thromboembolic events or aortic injury.

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.