Developing transmission line equations of oxygen transport for predicting oxygen distribution in the arterial system

Blood flow and emboli transport patterns during venoarterial extracorporeal membrane oxygenation: A computational fluid dynamics study

PROBLEM

Despite advances in Venoarterial Extracorporeal Membrane Oxygenation (VA-ECMO), a significant mortality rate persists due to complications. The non-physiological blood flow dynamics of VA-ECMO may lead to neurological complications and organ ischemia. Continuous retrograde high-flow oxygenated blood enters through a return cannula placed in the femoral artery which opposes the pulsatile deoxygenated blood ejected by the left ventricle (LV), which impacts upper body oxygenation and subsequent hyperoxemia. The complications underscore the critical need to comprehend the impact of VA-ECMO support level and return cannula size, as mortality remains a significant concern.

AIM

The aim of this study is to predict and provide insights into the complications associated with VA-ECMO using computational fluid dynamics (CFD) simulations. These complications will be assessed by characterising blood flow and emboli transport patterns through a comprehensive analysis of the influence of VA-ECMO support levels and arterial return cannula sizes.

METHODS

Patient-specific 3D aortic and major branch models, derived from a male patient's CT scan during VA-ECMO undergoing respiratory dysfunction, were analyzed using CFD. The investigation employed species transport and discrete particle tracking models to study ECMO blood (oxygenated) mixing with LV blood (deoxygenated) and to trace emboli transport patterns from potential sources (circuit, LV, and aorta wall). Two cannula sizes (15 Fr and 19 Fr) were tested alongside varying ECMO pump flow rates (50%, 70%, and 90% of the total cardiac output).

RESULTS

Cannula size did not significantly affect oxygen transport. At 90% VA-ECMO support, all arteries distal to the aortic arch achieved 100% oxygen saturation. As support level decreased, oxygen transport to the upper body also decreased to a minimum saturation of 73%. Emboli transport varied substantially between emboli origin and VAECMO support level, with the highest risk of cerebral emboli coming from the LV with a 15 Fr cannula at 90% support.

CONCLUSION

Arterial return cannula sizing minimally impacted blood oxygen distribution; however, it did influence the distribution of emboli released from the circuit and aortic wall. Notably, it was the support level alone that significantly affected the mixing zone of VA-ECMO and cardiac blood, subsequently influencing the risk of embolization of the cardiogenic source and oxygenation levels across various arterial branches.

From Experiments to Simulation: Shear-Induced Responses of Red Blood Cells to Different Oxygen Saturation Levels

Red blood cells (RBC) carry and deliver oxygen (O₂) to peripheral tissues through different microcirculatory regions where they are exposed to various levels of shear stress (SS). O₂ affinity of hemoglobin (Hb) decreases as the blood enters the microcirculation. This phenomenon determines Hb interactions with RBC membrane proteins that can further regulate the structure of cytoskeleton and affect the mechanical properties of cells. The goal of this study is to evaluate shear-induced RBC deformability and simulate RBC dynamics in blood flow under oxygenated and deoxygenated conditions. Venous blood samples from healthy donors were oxygenated with ambient air or deoxygenated with 100% nitrogen gas for 10 min and immediately applied into an ektacytometer (LORRCA). RBC deformability was measured before and after the application of continuous 5 Pa SS for 300 s by LORRCA and recorded as elongation index (EI) values. A computational model was generated for the simulation of blood flow in a real carotid artery section. EI distribution throughout the artery and its relationships with velocity, pressure, wall SS and viscosity were determined by computational tools. RBC deformability significantly increased in deoxygenation compared to oxygenated state both before and after 5 Pa SS implementation (p < 0.0001). However, EI values after continuous SS were not significant at higher SS levels (>5.15 Pa) in deoxygenated condition. Simulation results revealed that the velocity gradient dominates the generation of SS and the shear thinning effect of blood has a minor effect on it. Distribution of EI was calculated during oxygenation/deoxygenation which is 5-10 times higher around the vessel wall compared to the center of the lumen for sections of the pulsatile flow profile. The extent of RBC deformability increases as RBCs approach to the vessel wall in a real 3D artery model and this increment is higher for deoxygenated condition compared to the oxygenated state. Hypoxia significantly increases shear-induced RBC deformability. RBCs could regulate their own mechanical properties in blood flow by increasing their deformability in hypoxic conditions. Computational tools can be applied for defining hypoxia-mediated RBC deformability changes to monitor blood flow in hypoxic tissues.