In silico parametric analysis of femoro-jugular venovenous ECMO and return cannula dynamics

Abstract

Background: Increasingly, computational fluid dynamics (CFD) is helping explore the impact of variables like: cannula design/size/position/flow rate and patient physiology on venovenous (VV) extracorporeal membrane oxygenation (ECMO). Here we use a CFD model to determine what role cardiac output (CO) plays and to analyse return cannula dynamics. Methods: Using a patient-averaged model of the right atrium and venae cava, we virtually inserted a 19Fr return cannula and a 25Fr drainage cannula. Running large eddy simulations, we assessed cardiac output at: 3.5–6.5 L/min and ECMO flow rate at: 2–6 L/min. We analysed recirculation fraction (Rf), time-averaged wall shear stress (TAWSS), pressure, velocity, and turbulent kinetic energy (TKE) and extracorporeal flow fraction (EFF = ECMO flow rate/CO). Results: Increased ECMO flow rate and decreased CO (high EFF) led to increased Rf (R = 0.98, log fit). Negative pressures developed in the venae cavae at low CO and high ECMO flow (high CR). Mean return cannula TAWSS was >10 Pa for all ECMO flow rates, with majority of the flow exiting the tip (94.0–95.8 %). Conclusions: Our results underpin the strong impact of CO on VV ECMO. A simple metric like EFF, once supported by clinical data, might help predict Rf for a patient at a given ECMO flow rate. The return cannula imparts high shear stresses on the blood, largely a result of the internal diameter.