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.

Patient-Specific Atrial Hemodynamics of a Double Lumen Neonatal Cannula in Correct Caval Position.

Clinical success of pediatric veno-venous (VV) extracorporeal membrane oxygenation (ECMO) is associated with the double lumen cannula cardiovascular device design as well as its anatomic orientation in the atrium. The positions of cannula ports with respect to the vena cavae and the tricuspid valve are believed to play a significant role on device hemodynamics. Despite various improvements in ECMO catheters, especially for the neonatal and congenital heart patients, it is still challenging to select a catalogue size that would fit to most patients optimally. In effect, the local unfavorable blood flow characteristics of the cannula would translate to an overall loss of efficiency of the ECMO circuit. In this study, the complex flow regime of a neonatal double lumen cannula, positioned in a patient-specific right atrium, is presented for the first time in literature. A pulsatile computational fluid dynamics (CFD) solver that is validated for cardiovascular device flow regimes was used to perform the detailed flow, oxygenated blood transport, and site-specific blood damage analysis using an integrated cannula and right atrium model. A standard 13Fr double lumen cannula was scanned using micro-CT, reconstructed and simulated under physiologic flow conditions. User defined scalar transport equations allowed the quantification of the mixing and convection of oxygenated and deoxygenated blood as well as blood residence times and hemolysis build-up. Site-specific CFD analysis provided key insight into the hemodynamic challenges encountered in cannula design and the associated intra-atrial flow patterns. Due to neonatal flow conditions, an ultra high velocity infusion jet emanated from the infusion port and created a zone of major recirculation in the atrium. This flow regime influenced the delivery of the oxygenated blood to the tricuspid valve. Elevated velocities and complex gradients resulted in higher wall shear stresses (WSS) particularly at the infusion port having the highest value followed by the aspiration hole closest to the drainage port. Our results show that, in a cannula that is perfectly oriented in the atrium, almost 38% of the oxygenated blood is lost to the atrial circulation while only half of the blood from inferior vena cava (IVC) can reach to the tricuspid valve. As such, approximately 6% of venous blood from superior vena cava (SVC) can be delivered to tricuspid. High values of hemolysis index were observed with blood damage encountered around infusion hole (0.025%). These results warrant further improvements in the cannula design to achieve optimal performance of ECMO and better patient outcomes.

Continuum modeling of thrombus formation and growth under different shear rates.

Obstruction of blood flow due to thrombosis is a major cause of ischemic stroke, myocardial infarction, and in severe cases, mortality. In particular, in blood wetting medical devices, thrombosis is a common reason for failure. The prediction of thrombosis by understanding signaling pathways using computational models, lead to identify the risk of thrombus formation in blood-contacting devices and design improvements. In this study, a mathematical model of thrombus formation and growth is presented. A biochemical model of platelet activation and aggregation is developed to predict thrombus size and shape at the site of vascular injury. Computational fluid dynamics using the finite volume method is employed to compute the velocity and pressure fields which are influenced by the growing thrombi. The passive transport of platelets, agonists, the platelet activation kinetics, their adhesion to the growing thrombi and embolization of platelets are solved by a fully coupled set of convection-diffusion-reaction equations. The thrombogenic surface representing blood-contacting material or injured blood vessel was incorporated into the model as a surface flux boundary condition to initiate thrombus formation. The blood is considered as a Newtonian fluid, while the thrombus is treated as a porous medium. The results are compared with in vitro experiments of a microfluidic chamber at an initial inlet venous shear rate of 200s-1 using a pressure-inlet boundary condition. The thrombus development due to agonist concentrations and change in the shear rate as well as thromboembolism for this benchmark problem is successfully computed. Furthermore, to extend the current model to a physiologically relevant configuration, thrombus formation in a blood tube is simulated. Two different heterogeneous reaction rates for platelet aggregation are used to simulate thrombus growth under a constant inlet flow rate. The findings show that the thrombus shape can be substantially altered by the hemodynamic conditions, increase in the shear rate and due to the combined effects of shear induced platelet activation and the heterogeneous reaction rates. It is also concluded that the model is able to predict thrombus formation in different physiological and pathological hemodynamics.

A novel Fontan Y-graft for interrupted inferior vena cava and azygos continuation.

OBJECTIVES: To evaluate the hemodynamicdynamic advantage of a new Fontan surgical template that is intended for complex single-ventricle patients with interrupted inferior vena cava-azygos and hemi-azygos continuation. The new technique has emerged from a comprehensive pre-surgical simulation campaign conducted to facilitate a balanced hepatic flow and somatic Fontan pathway growth after Kawashima procedure. METHODS: For 9 patients, aged 2 to18 years, majority having poor preoperative oxygen saturation, a pre-surgical computational fluid dynamics customization is conducted. Both the traditional Fontan pathways and the proposed novel Y-graft templates are considered. Numerical model was validated against in vivo phase-contrast magnetic resonance imaging data and in vitro experiments. RESULTS: The proposed template is selected and executed for 6 out of the 9 patients based on its predicted superior hemodynamic performance. Pre-surgical simulations performed for this cohort indicated that flow from the hepatic veins (HEP) do not reach to the desired lung. The novel Y-graft template, customized via a right- or left-sided displacement of the total cavopulmonary connection anastomosis location resulted a drastic increase in HEP flow to the desired lung. Orientation of HEP to azygos direct shunt is found to be important as it can alter the flow pattern from 38% in the caudally located direct shunt to 3% in the cranial configuration with significantly reversed flow. The postoperative measurements prove that oxygen saturation increased significantly (P-value = 0.00009) to normal levels in 1 year follow-up. CONCLUSIONS: The new Y-graft template, if customized for the individual patient, is a viable alternative to the traditional surgical pathways. This template addresses the competing hemodynamic design factors of low physiological venous pressure, high postoperative oxygen saturation, low energy loss and balanced hepatic growth factor distribution possibly assuring adequate lung development. DATE AND NUMBER OF IRB APPROVAL: 25 October 2019, 280011928-604.01.01.

Hemodynamic performance limits of the neonatal Double-Lumen cannula.

Venovenous extracorporeal membrane oxygenation (VV-ECMO) is the preferred surgical intervention for patients suffering from severe cardiorespiratory failure, also encountered in SARS-Cov-2 management. The key component of VV-ECMO is the double-lumen cannula (DLC) that enables single-site access. The biofluid dynamics of this compact device is particularly challenging for neonatal patients due to high Reynolds numbers, tricuspid valve location and right-atrium hemodynamics. In this paper we present detailed findings of our comparative analysis of the right-atrial hemodynamics and salient design features of the 13Fr Avalon Elite DLC (as the clinically preferred neonatal cannula) with the alternate Origen DLC design, using experimentally validated computational fluid dynamics. Highly accurate 3D-reconstructions of both devices were obtained through an integrated optical coherence tomography and micro-CT imaging approach. Both cannula configurations displayed complex flow structures inside the atrium, superimposed over predominant recirculation regimes. We found that the Avalon DLC performed significantly better than the Origen alternative, by capturing 80% and 94% of venous blood from the inferior and superior vena cavae, respectively and infusing the oxygenated blood with an efficiency of more than 85%. The micro-scale geometric design features of the Avalon DLC that are associated with superior hemodynamics were investigated through 14 parametric cannula configurations. These simulations showed that the strategic placement of drainage holes, the smooth infusion blood stream diverter and efficient distribution of the venous blood capturing area between the vena cavae are associated with robust blood flow performance. Nevertheless, our parametric results indicate that there is still room for further device optimization beyond the performance measurements for both Avalon and Origen DLC in this study. In particular, the performance envelope of malpositioned cannula and off-design conditions require additional blood flow simulations for analysis.

Patient-Specific Hemodynamics of New Coronary Artery Bypass Configurations.

PURPOSE: This study aims to quantify the patient-specific hemodynamics of complex conduit routing configurations of coronary artery bypass grafting (CABG) operation which are specifically suitable for off-pump surgeries. Coronary perfusion efficacy and local hemodynamics of multiple left internal mammary artery (LIMA) with sequential and end-to-side anastomosis are investigated. Using a full anatomical model comprised of aortic arch and coronary artery branches the optimum perfusion configuration in multi-vessel coronary artery stenosis is desired. METHODOLOGY: Two clinically relevant CABG configurations are created using a virtual surgical planning tool where for each configuration set, the stenosis level, anastomosis distance and angle were varied. A non-Newtonian computational fluid dynamics solver in OpenFOAM incorporated with resistance boundary conditions representing the coronary perfusion physiology was developed. The numerical accuracy is verified and results agreed well with a validated commercial cardiovascular flow solver and experiments. For segmental performance analysis, new coronary perfusion indices to quantify deviation from the healthy scenario were introduced. RESULTS: The first simulation configuration set;-a CABG targeting two stenos sites on the left anterior descending artery (LAD), the LIMA graft was capable of 31 mL/min blood supply for all the parametric cases and uphold the healthy LAD perfusion in agreement with the clinical experience. In the second end-to-side anastomosed graft configuration set;-the radial artery graft anastomosed to LIMA, a maximum of 64 mL/min flow rate in LIMA was observed. However, except LAD, the obtuse marginal (OM) and second marginal artery (m2) suffered poor perfusion. In the first set, average wall shear stress (WSS) were in the range of 4 to 35 dyns/cm2 for in LAD. Nevertheless, for second configuration sets the WSS values were higher as the LIMA could not supply enough blood to OM and m2. CONCLUSION: The virtual surgical configurations have the potential to improve the quality of operation by providing quantitative surgical insight. The degree of stenosis is a critical factor in terms of coronary perfusion and WSS. The sequential anastomosis can be done safely if the anastomosis angle is less than 90 degrees regardless of degree of stenosis. The smaller proposed perfusion index value, O(0.04 - 0) × 102, enable us to quantify the post-op hemodynamic performance by comparing with the ideal healthy physiological flow.