Simulations reveal adverse hemodynamics in patients with multiple systemic to pulmonary shunts.

For newborns diagnosed with pulmonary atresia or severe pulmonary stenosis leading to insufficient pulmonary blood flow, cyanosis can be mitigated with placement of a modified Blalock-Taussig shunt (MBTS) between the innominate and pulmonary arteries. In some clinical scenarios, patients receive two systemic-to-pulmonary connections, either by leaving the patent ductus arteriosus (PDA) open or by adding an additional central shunt (CS) in conjunction with the MBTS. This practice has been motivated by the thinking that an additional source of pulmonary blood flow could beneficially increase pulmonary flow and provide the security of an alternate pathway in case of thrombosis. However, there have been clinical reports of premature shunt occlusion when more than one shunt is employed, leading to speculation that multiple shunts may in fact lead to unfavorable hemodynamics and increased mortality. In this study, we hypothesize that multiple shunts may lead to undesirable flow competition, resulting in increased residence time (RT) and elevated risk of thrombosis, as well as pulmonary overcirculation. Computational fluid dynamics-based multiscale simulations were performed to compare a range of shunt configurations and systematically quantify flow competition, pulmonary circulation, and other clinically relevant parameters. In total, 23 cases were evaluated by systematically changing the PDA/CS diameter, pulmonary vascular resistance (PVR), and MBTS position and compared by quantifying oxygen delivery (OD) to the systemic and coronary beds, wall shear stress (WSS), oscillatory shear index (OSI), WSS gradient (WSSG), and RT in the pulmonary artery (PA), and MBTS. Results showed that smaller PDA/CS diameters can lead to flow conditions consistent with increased thrombus formation due to flow competition in the PA, and larger PDA/CS diameters can lead to insufficient OD due to pulmonary hyperfusion. In the worst case scenario, it was found that multiple shunts can lead to a 160% increase in RT and a 10% decrease in OD. Based on the simulation results presented in this study, clinical outcomes for patients receiving multiple shunts should be critically investigated, as this practice appears to provide no benefit in terms of OD and may actually increase thrombotic risk.

[1]  B. Marino,et al.  Simultaneous patency of ductus arteriosus and surgical shunt in pulmonary atresia with intact ventricular septum. A cause of acute myocardial failure? , 1986, Scandinavian journal of thoracic and cardiovascular surgery.

[2]  Jeanette P. Schmidt,et al.  The Simbios National Center: Systems Biology in Motion , 2008, Proceedings of the IEEE.

[3]  A. L. Marsden,et al.  Computation of residence time in the simulation of pulsatile ventricular assist devices , 2014 .

[4]  K S Sakariassen,et al.  Shear-induced platelet activation and platelet microparticle formation at blood flow conditions as in arteries with a severe stenosis. , 1997, Arteriosclerosis, thrombosis, and vascular biology.

[5]  Marek Behr,et al.  Parallel finite-element computation of 3D flows , 1993, Computer.

[6]  Yuri Bazilevs,et al.  A bi-partitioned iterative algorithm for solving linear systems arising from incompressible flow problems , 2015 .

[7]  Danny Bluestein,et al.  Fluid mechanics of arterial stenosis: Relationship to the development of mural thrombus , 1997, Annals of Biomedical Engineering.

[8]  D A Steinman,et al.  Particle volumetric residence time calculations in arterial geometries. , 1996, Journal of biomechanical engineering.

[9]  T. Hughes,et al.  Effect of exercise on hemodynamic conditions in the abdominal aorta. , 1999, Journal of vascular surgery.

[10]  T. Tezduyar,et al.  Mesh Moving Techniques for Fluid-Structure Interactions With Large Displacements , 2003 .

[11]  Giancarlo Pennati,et al.  assessment by computational and in vitro studies of the blood flow rate through modified blalock-taussig shunts , 2004, Cardiology in the Young.

[12]  M. Jacobs,et al.  Fontan procedure for hypoplastic left heart syndrome. , 1992, The Annals of thoracic surgery.

[13]  Georgios C. Georgiou,et al.  The influence of temperature on rheological properties of blood mixtures with different volume expanders—implications in numerical arterial hemodynamics simulations , 2011 .

[14]  A. Marsden,et al.  The assisted bidirectional Glenn: a novel surgical approach for first-stage single-ventricle heart palliation. , 2015, The Journal of thoracic and cardiovascular surgery.

[15]  Leonie Rouleau,et al.  Endothelial cell morphologic response to asymmetric stenosis hemodynamics: effects of spatial wall shear stress gradients. , 2010, Journal of biomechanical engineering.

[16]  O Barnea,et al.  Estimation of oxygen delivery in newborns with a univentricular circulation. , 1998, Circulation.

[17]  A. L. Marsden,et al.  Impact of data distribution on the parallel performance of iterative linear solvers with emphasis on CFD of incompressible flows , 2015 .

[18]  S. Nicolson,et al.  Risk factors for mortality after the Norwood procedure. , 2002, European journal of cardio-thoracic surgery : official journal of the European Association for Cardio-thoracic Surgery.

[19]  V T Turitto,et al.  Mechanical factors affecting hemostasis and thrombosis. , 1998, Thrombosis research.

[20]  H. Taussig,et al.  The Surgical Treatment of Malformations of the Heart: In Which There Is Pulmonary Stenosis or Pulmonary Atresia , 1945 .

[21]  F. Fontan,et al.  Surgical repair of tricuspid atresia , 1971, Thorax.

[22]  R. Poruban,et al.  A comparison of Blalock-Taussig shunts with and without closure of the ductus arteriosus in neonates with pulmonary atresia. , 2011, The Annals of thoracic surgery.

[23]  Wei Yin,et al.  The effect of physiologically relevant dynamic shear stress on platelet and endothelial cell activation. , 2011, Thrombosis research.

[24]  Alison L. Marsden,et al.  A modular numerical method for implicit 0D/3D coupling in cardiovascular finite element simulations , 2013, J. Comput. Phys..

[25]  C. A. Figueroa,et al.  Computational simulations of hemodynamic changes within thoracic, coronary, and cerebral arteries following early wall remodeling in response to distal aortic coarctation , 2012, Biomechanics and Modeling in Mechanobiology.

[26]  P. Vouhé,et al.  Modified Blalock-Taussig shunts: results in infants less than 3 months of age. , 1990, The Annals of thoracic surgery.

[27]  U Bartram,et al.  Causes of death after the modified Norwood procedure: a study of 122 postmortem cases. , 1997, The Annals of thoracic surgery.

[28]  Christopher P. Cheng,et al.  Quantification of Particle Residence Time in Abdominal Aortic Aneurysms Using Magnetic Resonance Imaging and Computational Fluid Dynamics , 2011, Annals of Biomedical Engineering.

[29]  Alison L. Marsden,et al.  Optimization of shunt placement for the Norwood surgery using multi-domain modeling. , 2012, Journal of biomechanical engineering.

[30]  Zhijie Wang,et al.  Complex Hemodynamics at the Apex of an Arterial Bifurcation Induces Vascular Remodeling Resembling Cerebral Aneurysm Initiation , 2007, Stroke.

[31]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[32]  Tayfun E. Tezduyar,et al.  Space-time finite element techniques for computation of fluid-structure interactions , 2005 .

[33]  Tain-Yen Hsia,et al.  A non-discrete method for computation of residence time in fluid mechanics simulations. , 2013, Physics of fluids.

[34]  D. Ku,et al.  Pulsatile flow in the human left coronary artery bifurcation: average conditions. , 1996, Journal of biomechanical engineering.

[35]  Gerhard A. Holzapfel,et al.  Nonlinear Solid Mechanics: A Continuum Approach for Engineering Science , 2000 .

[36]  Giancarlo Pennati,et al.  Use of mathematic modeling to compare and predict hemodynamic effects of the modified Blalock-Taussig and right ventricle-pulmonary artery shunts for hypoplastic left heart syndrome. , 2008, The Journal of thoracic and cardiovascular surgery.

[37]  Zhiliang Xu,et al.  A multiscale model of thrombus development , 2008, Journal of The Royal Society Interface.

[38]  Giancarlo Pennati,et al.  Modeling of systemic-to-pulmonary shunts in newborns with a univentricular circulation: State of the art and future directions , 2010 .

[39]  G. Hulbert,et al.  A generalized-α method for integrating the filtered Navier–Stokes equations with a stabilized finite element method , 2000 .

[40]  F. Migliavacca,et al.  Multiscale modeling of the cardiovascular system: application to the study of pulmonary and coronary perfusions in the univentricular circulation. , 2005, Journal of biomechanics.

[41]  F. Migliavacca,et al.  Modeling of the Norwood circulation: effects of shunt size, vascular resistances, and heart rate. , 2001, American journal of physiology. Heart and circulatory physiology.

[42]  W W GLENN,et al.  FURTHER EXPERIMENTS ON LONG TERM SURVIVORS AFTER CIRCULATORY BYPASS OF THE RIGHT SIDE OF THE HEART. , 1958, Surgery, gynecology & obstetrics.

[43]  A. Marsden,et al.  A comparison of outlet boundary treatments for prevention of backflow divergence with relevance to blood flow simulations , 2011 .

[44]  Yuri Bazilevs,et al.  A new preconditioning technique for implicitly coupled multidomain simulations with applications to hemodynamics , 2013 .

[45]  G B Fiore,et al.  In vitro steady-flow analysis of systemic-to-pulmonary shunt haemodynamics. , 2001, Journal of biomechanics.

[46]  Yuri Bazilevs,et al.  Computational fluid–structure interaction: methods and application to a total cavopulmonary connection , 2009 .

[47]  T. Tezduyar,et al.  Arterial fluid mechanics modeling with the stabilized space–time fluid–structure interaction technique , 2008 .

[48]  Thomas J. R. Hughes,et al.  Patient-specific isogeometric fluid–structure interaction analysis of thoracic aortic blood flow due to implantation of the Jarvik 2000 left ventricular assist device , 2009 .

[49]  James E. Lock,et al.  Rest and exercise hemodynamics after the Fontan procedure. , 1982 .

[50]  T. Hughes,et al.  Isogeometric fluid-structure interaction: theory, algorithms, and computations , 2008 .

[51]  L. Franca,et al.  Stabilized finite element methods. II: The incompressible Navier-Stokes equations , 1992 .

[52]  C. C. Long,et al.  Fluid–structure interaction simulations of the Fontan procedure using variable wall properties , 2012, International journal for numerical methods in biomedical engineering.