Research Papers

Simulations Reveal Adverse Hemodynamics in Patients With Multiple Systemic to Pulmonary Shunts

[+] Author and Article Information
Mahdi Esmaily-Moghadam, Alison Marsden

Mechanical and Aerospace
Engineering Department,
University of California,
San Diego, CA 92093

Bari Murtuza, Tain-Yen Hsia

Cardiac Unit,
Great Ormond Street Hospital for Children and Institute of Child Health,
London WC1N 3JH, UK

1Part of the Modeling of Congenital Hearts Alliance (MOCHA). MOCHA Investigators: Edward Bove and Adam Dorfman (University of Michigan); Andrew Taylor, Alessandro Giardini, Sachin Khambadkone, Marc de Leval, Silvia Schievano, and T.-Y. Hsia (Institute of Child Health, UK); G. Hamilton Baker and Anthony Hlavacek (Medical University of South Carolina); Francesco Migliavacca, Giancarlo Pennati, and Gabriele Dubini (Politecnico di Milano, Italy); Richard Figliola and John McGregor (Clemson University); Alison Marsden (University of California, San Diego); and Irene Vignon-Clementel (National Institute of Research in Informatics and Automation, France).

Manuscript received July 17, 2013; final manuscript received December 2, 2014; published online January 29, 2015. Assoc. Editor: Francis Loth.

J Biomech Eng 137(3), 031001 (Mar 01, 2015) (12 pages) Paper No: BIO-13-1320; doi: 10.1115/1.4029429 History: Received July 17, 2013; Revised December 02, 2014; Online January 29, 2015

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.

Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.


Blalock, A., and Taussig, H. B., 1945, “The Surgical Treatment of Malformations of the Heart,” J. Am. Med. Assoc., 128(3), pp. 189–202. [CrossRef]
Tamisier, D., Vouhe, P., Vernant, F., Leca, F., Massot, C., and Neveux, J., 1990, “Modified Blalock–Taussig Shunts: Results in Infants Less Than 3 Months of Age,” Ann. Thorac. Surg., 49(5), pp. 797–801. [CrossRef] [PubMed]
Glenn, W. W., 1958, “Circulatory Bypass of the Right Side of the Heart: Shunt Between Superior Vena Cava and Distal Right Pulmonary Artery-Report of Clinical Application,” N. Engl. J. Med., 259(3), pp. 117–120. [CrossRef] [PubMed]
Fontan, F., and Baudet, E., 1971, “Surgical Repair of Tricuspid Atresia,” Thorax, 26(3), pp. 240–248. [CrossRef] [PubMed]
Norwood, W., Jacobs, M., and Murphy, J., 1992, “Fontan Procedure for Hypoplastic Left Heart Syndrome,” Ann. Thorac. Surg., 54(6), pp. 1025–1030. [CrossRef] [PubMed]
Shachar, G. B., Fuhrman, B. P., Wang, Y., Lucas, R., and Lock, J. E., 1982, “Rest and Exercise Hemodynamics After the Fontan Procedure,” Circulation, 65(6), pp. 1043–1048. [CrossRef] [PubMed]
Gaynor, J. W., Mahle, W. T., Cohen, M. I., Ittenbach, R. F., DeCampli, W. M., Steven, J. M., Nicolson, S. C., and Spray, T. L., 2002, “Risk Factors for Mortality After the Norwood Procedure,” Eur. J. Cardiothorac. Surg., 22(1), pp. 82–89. [CrossRef] [PubMed]
Bartram, U., Grünenfelder, J., and Praagh, R. V., 1997, “Causes of Death After the Modified Norwood Procedure: A Study of 122 Postmortem Cases,” Ann. Thorac. Surg., 64(6), pp. 1795–1802. [CrossRef] [PubMed]
Corno, A., Mazzera, E., Marino, B., Parisi, F., and Marcelletti, C., 1986, “Simultaneous Patency of Ductus Arteriosus and Surgical Shunt in Pulmonary Atresia With Intact Ventricular Septum: A Cause of Acute Myocardial Failure?,” Scand. Cardiovasc. J., 20(2), pp. 123–127. [CrossRef]
Zahorec, M., Hrubsova, Z., Skrak, P., Poruban, R., Nosal, M., and Kovacikova, L., 2011, “A Comparison of Blalock–Taussig Shunts With and Without Closure of the Ductus Arteriosus in Neonates With Pulmonary Atresia,” Ann. Thorac. Surg., 92(2), pp. 653–658. [CrossRef] [PubMed]
Migliavacca, F., Pennati, G., Dubini, G., Fumero, R., Pietrabissa, R., Urcelay, G., Bove, E., Hsia, T., and de Leval, M., 2001, “Modeling of the Norwood Circulation: Effects of Shunt Size, Vascular Resistances, and Heart Rate,” Am. J. Physiol.: Heart Circ. Physiol., 280(5), pp. H2076–H2086. [PubMed]
Pennati, G., Migliavacca, F., Dubini, G., and Bove, E., 2010, “Modeling of Systemic-to-Pulmonary Shunts in Newborns With a Univentricular Circulation: State of the Art and Future Directions,” Prog. Pediatr. Cardiol., 30(1–2), pp. 23–29. [CrossRef]
Esmaily-Moghadam, M., Migliavacca, F., Vignon-Clementel, I., Hsia, T., and Marsden, A., 2012, “Optimization of Shunt Placement for the Norwood Surgery Using Multi-Domain Modeling,” ASME J. Biomech. Eng., 134(5), p. 051002. [CrossRef]
Esmaily-Moghadam, M., Hsia, T.-Y., and Marsden, A., 2014, “The Assisted Bidirectional Glenn: A Novel Surgical Approach for First Stage Single Ventricle Heart Palliation,” J. Thorac. Cardiovasc. Surg. (in press). [CrossRef]
Bove, E., Migliavacca, F., de Leval, M., Balossino, R., Pennati, G., Lloyd, T., Khambadkone, S., Hsia, T., and Dubini, G., 2008, “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,” J. Thorac. Cardiovasc. Surg., 136(2), pp. 312–320.e2. [CrossRef] [PubMed]
Laganà, K., Balossino, R., Migliavacca, F., Pennati, G., Bove, E., de Leval, M., and Dubini, G., 2005, “Multiscale Modeling of the Cardiovascular System: Application to the Study of Pulmonary and Coronary Perfusions in the Univentricular Circulation,” J. Biomech., 38(5), pp. 1129–1141. [CrossRef] [PubMed]
Pennati, G., Migliavacca, F., Gervaso, F., and Dubini, G., 2004, “Assessment by Computational and In Vitro Studies of the Blood Flow Rate Through Modified Blalock–Taussig Shunts,” Cardiol. Young, 14(Suppl. 3), pp. 24–29. [CrossRef] [PubMed]
Pennati, G., Fiore, G. B., Migliavacca, F., Laganà, K., Fumero, R., and Dubini, G., 2001, “In Vitro Steady-Flow Analysis of Systemic-to-Pulmonary Shunt Haemodynamics,” J. Biomech., 34(1), pp. 23–30. [CrossRef] [PubMed]
Esmaily-Moghadam, M., Vignon-Clementel, I., Figliola, R., and Marsden, A., 2013, “A Modular Numerical Method for Implicit 0D/3D Coupling in Cardiovascular Finite Element Simulations,” J. Comput. Phys., 224, pp. 63–79. [CrossRef]
Bluestein, D., Niu, L., Schoephoerster, R., and Dewanjee, M., 1997, “Fluid Mechanics of Arterial Stenosis: Relationship to the Development of Mural Thrombus,” Ann. Biomed. Eng., 25(2), pp. 344–356. [CrossRef] [PubMed]
Esmaily-Moghadam, M., Hsia, T.-Y., and Marsden, A., 2013, “A Non-Discrete Method for Computation of Residence Time in Fluid Mechanics Simulations,” Phys. Fluids, 25(11), p. 110802. [CrossRef]
Kunov, M., Steinman, D., and Ethier, C., 1996, “Particle Volumetric Residence Time Calculations in Arterial Geometries,” ASME J. Biomech. Eng., 118(2), pp. 158–164. [CrossRef]
Suh, G., Les, A., Tenforde, A., Shadden, S., Spilker, R., Yeung, J., Cheng, C., Herfkens, R., Dalman, R., and Taylor, C., 2011, “Quantification of Particle Residence Time in Abdominal Aortic Aneurysms Using Magnetic Resonance Imaging and Computational Fluid Dynamics,” Ann. Biomed. Eng., 39(2), pp. 864–883. [CrossRef] [PubMed]
Schmidt, J., Delp, S., Sherman, M., Taylor, C., Pande, V., and Altman, R., 2008, “The Simbios National Center: Systems Biology in Motion,” Proc. IEEE Inst. Electr. Electron Eng., 96(8), pp. 1266–1280. [CrossRef] [PubMed]
Valant, A. Z., Žiberna, L., Papaharilaou, Y., Anayiotos, A., and Georgiou, G., 2011, “The Influence of Temperature on Rheological Properties of Blood Mixtures With Different Volume Expanders-Implications in Numerical Arterial Hemodynamics Simulations,” Rheol. Acta, 50(4), pp. 389–402. [CrossRef]
Jansen, K., Whiting, C., and Hulbert, G., 2000, “A Generalized-[Alpha] Method for Integrating the Filtered Navier–Stokes Equations With a Stabilized Finite Element Method,” Comput. Methods Appl. Mech. Eng., 190(3–4), pp. 305–319. [CrossRef]
Franca, L., and Frey, S., 1992, “Stabilized Finite Element Methods: II. The Incompressible Navier–Stokes Equations,” Comput. Methods Appl. Mech. Eng., 99(2–3), pp. 209–233. [CrossRef]
Brooks, A., and Hughes, T., 1982, “Streamline Upwind/Petrov–Galerkin Formulations for Convection Dominated Flows With Particular Emphasis on the Incompressible Navier–Stokes Equations,” Comput. Methods Appl. Mech. Eng., 32(1–3), pp. 199–259. [CrossRef]
Esmaily-Moghadam, M., Bazilevs, Y., and Marsden, A., 2014, “A Bi-Partitioned Iterative Algorithm for Solving Linear Systems Arising From Incompressible Flow Problems,” Comput. Methods Appl. Mech. Eng. (in press). [CrossRef]
Esmaily-Moghadam, M., Bazilevs, Y., and Marsden, A. L., 2013, “A New Preconditioning Technique for Implicitly Coupled Multidomain Simulations With Applications to Hemodynamics,” Comput. Mech., 52(5), pp. 1141–1152. [CrossRef]
Esmaily-Moghadam, M., Bazilevs, Y., and Marsden, A., 2014, “Impact of Data Distribution on the Parallel Performance of Iterative Linear Solvers With Emphasis on CFD of Incompressible Flows,” Comput. Mech. (published online). [CrossRef]
Esmaily-Moghadam, M., Bazilevs, Y., Hsia, T., Vignon-Clementel, I., and Marsden, A., 2011, “A Comparison of Outlet Boundary Treatments for Prevention of Backflow Divergence With Relevance to Blood Flow Simulations,” Comput. Mech., 48(3), pp. 277–291. [CrossRef]
Bazilevs, Y., Gohean, J., Hughes, T., Moser, R., and Zhang, Y., 2009, “Patient-Specific Isogeometric Fluid-Structure Interaction Analysis of Thoracic Aortic Blood Flow due to Implantation of the Jarvik 2000 Left Ventricular Assist Device,” Comput. Methods Appl. Mech. Eng., 198(45–46), pp. 3534–3550. [CrossRef]
Long, C., Esmaily-Moghadam, M., Marsden, A., and Bazilevs, Y., 2014, “Computation of Residence Time in the Simulation of Pulsatile Ventricular Assist Devices,” Comput. Mech., 54(4), pp. 911–919. [CrossRef]
Figueroa, C., 2006, “A Coupled-Momentum Method to Model Blood Flow and Vessel Deformation in Human Arteries: Applications in Disease Research and Simulation-Based Medical Planning,” Ph.D. thesis, Stanford University.
He, X., and Ku, D., 1996, “Pulsatile Flow in the Human Left Coronary Artery Bifurcation: Average Conditions,” ASME J. Biomech. Eng., 118(1), pp. 74–82. [CrossRef]
Taylor, C., Hughes, T., and Zarins, C., 1999, “Effect of Exercise on Hemodynamic Conditions in the Abdominal Aorta,” J. Vasc. Surg., 29(6), pp. 1077–1089. [CrossRef] [PubMed]
Barnea, O., Santamore, W., Rossi, A., Salloum, E., Chien, S., and Austin, E., 1998, “Estimation of Oxygen Delivery in Newborns With a Univentricular Circulation,” Circulation, 98(14), pp. 1407–1413. [CrossRef] [PubMed]
Holzapfel, G. A., 2000, Nonlinear Solid Mechanics: A Continuum Approach for Engineering, John Wiley & Sons, Ltd., West Sussex, UK.
Tezduyar, T. E., Sathe, S., Schwaab, M., and Conklin, B. S., 2008, “Arterial Fluid Mechanics Modeling With the Stabilized Space–Time Fluid–Structure Interaction Technique,” Int. J. Numer. Methods Fluids, 57(5), pp. 601–629. [CrossRef]
Long, C., Hsu, M., Bazilevs, Y., and Feinstein, J., 2012, “A. Marsden, Fluid-Structure Interaction Simulations of the Fontan Procedure Using Variable Wall Properties,” Int. J. Numer. Methods Biomed. Eng., 28(5), pp. 513–527. [CrossRef]
Coogan, J. S., Humphrey, J. D., and Figueroa, C. A., 2013, “Computational Simulations of Hemodynamic Changes Within Thoracic, Coronary, and Cerebral Arteries Following Early Wall Remodeling in Response to Distal Aortic Coarctation,” Biomech. Model. Mechanobiol., 12(1), pp. 79–93. [CrossRef] [PubMed]
Bazilevs, Y., Hsu, M., Benson, D., Sankaran, S., and Marsden, A., 2009, “Computational Fluid-Structure Interaction: Methods and Application to a Total Cavopulmonary Connection,” Comput. Mech., 45(1), pp. 77–89. [CrossRef]
Bazilevs, Y., Calo, V., Hughes, T., and Zhang, Y., 2008, “Isogeometric Fluid-Structure Interaction: Theory, Algorithms, and Computations,” Comput. Mech., 43(1), pp. 3–37. [CrossRef]
Tezduyar, T. E., Sathe, S., Keedy, R., and Stein, K., 2006, “Space-Time Finite Element Techniques for Computation of Fluid-Structure Interactions,” Comput. Methods Appl. Mech. Eng., 195(17), pp. 2002–2027. [CrossRef]
Tezduyar, T., Aliabadi, S., Behr, M., Johnson, A., and Mittal, S., 1993, “Parallel Finite-Element Computation of 3D Flows,” Computer, 26(10), pp. 27–36. [CrossRef]
Stein, K., Tezduyar, T., and Benney, R., 2003, “Mesh Moving Techniques for Fluid-Structure Interactions With Large Displacements,” ASME J. Appl. Mech., 70(1), pp. 58–63. [CrossRef]
Rouleau, L., Farcas, M., Tardif, J., Mongrain, R., and Leask, R., 2010, “Endothelial Cell Morphologic Response to Asymmetric Stenosis Hemodynamics: Effects of Spatial Wall Shear Stress Gradients,” ASME J. Biomech. Eng., 132(8), p. 081013. [CrossRef]
Xu, Z., Chen, N., Kamocka, M., Rosen, E., and Alber, M., 2008, “A Multiscale Model of Thrombus Development,” J. R. Soc., Interface, 5(24), pp. 705–722. [CrossRef]
Meng, H., Wang, Z., Hoi, Y., Gao, L., Metaxa, E., Swartz, D., and Kolega, J., 2007, “Complex Hemodynamics at the Apex of an Arterial Bifurcation Induces Vascular Remodeling Resembling Cerebral Aneurysm Initiation,” Stroke, 38(6), pp. 1924–1931. [CrossRef] [PubMed]
Turitto, V., and Hall, C., 1998, “Mechanical Factors Affecting Hemostasis and Thrombosis,” Thromb. Res., 92 (6, Supplement 2), pp. S25–S31. [CrossRef] [PubMed]
Holme, P., Orvim, U., Hamers, M., Solum, N., Brosstad, F., Barstad, R., and Sakariassen, K., 1997, “Shear-Induced Platelet Activation and Platelet Microparticle Formation at Blood Flow Conditions as in Arteries With a Severe Stenosis,” Arterioscler., Thromb., Vasc. Biol., 17(4), pp. 646–653. [CrossRef]
Yin, W., Shanmugavelayudam, S., and Rubenstein, D., 2011, “The Effect of Physiologically Relevant Dynamic Shear Stress on Platelet and Endothelial Cell Activation,” Thromb. Res., 127(3), pp. 235–241. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 1

The constructed idealized models: (a) distal MBTS and a 4.0 mm PDA, (b) proximal MBTS and a 3.0 mm PDA, and (c) proximal MBTS and a 4.0 mm CS. The region between CS/PDA and MBTS inside the PA (red) and the region in the MBTS (blue) are used to compare the results of studied cases.

Grahic Jump Location
Fig. 2

The LPN, which is coupled to the MBTS anatomy, contains blocks for the upper body arteries (UBA), upper body capillary bed (UBB), upper body veins (UBV), pulmonary artery bed (PAB), pulmonary vein bed (PVB), lower body arteries (LBA), lower body capillary bed (LBB), lower body veins (LBV), two coronary arteries (CA1 and CA2), coronary capillary bed (CB), coronary veins (CV), left atrium (LA), right atrium (RA), and single ventricle (SV). The ascending aorta (AA), descending aorta (AoD), brachiocephalic artery (BA), right common carotid artery (RCCA), left common carotid artery (LCCA), left subclavian artery (LSA), left PA (LPA), right PA (RPA), and right coronary artery (RCA) are shown in the 3D model. Note the left coronary artery, which is omitted here to make the schematic less crowded, is connected to an LPN block identical to that of the RCA.

Grahic Jump Location
Fig. 3

Wall thickness is calculated by solving the Laplace equation (left figure). Then, the wall of the fluid domain mesh is extruded in the normal direction to generate the solid domain mesh, i.e., vasculature wall (middle and right figures).

Grahic Jump Location
Fig. 4

Deformation of the vessel walls at the peak systole

Grahic Jump Location
Fig. 5

Systemic (top plot) and coronary oxygen deliveries (bottom plot) for proximal (solid) and distal (dashed) MBTS and normal PVR (black) and high PVR (red) versus PDA/CS diameter. Note the single MBTS corresponds to DPDA = 0. The results of three simulations with CS are shown with red circles, in which DPDA denotes CS diameter. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Grahic Jump Location
Fig. 6

Time averaged velocity field in normal PVR (top row) and high PVR (bottom row), and a single proximal MBTS (left column), a proximal MBTS and a 2 mm PDA (middle column), and a proximal MBTS and a 4 mm PDA (right column)

Grahic Jump Location
Fig. 7

(a) Pressure contours for normal PVR (top row) and high PVR (bottom row), and a single proximal MBTS (left column), a proximal MBTS and a 2 mm PDA (middle column) cases, and a proximal MBTS and a 4 mm PDA (right column). (b) Averaged pressure in the PA outlets (top plot) and in the AA (bottom plot). See Fig. 5 caption for more details.

Grahic Jump Location
Fig. 8

(a) RT contours for normal PVR (top row) and high PVR (bottom row), and a single proximal MBTS (left column), a proximal MBTS and a 2 mm PDA (middle column), and a proximal MBTS and a 4 mm PDA (right column). (b) Spatially averaged RT in the PA segment (top plot) and in the MBTS (bottom plot). See Fig. 5 caption for more details.

Grahic Jump Location
Fig. 9

(a) WSS contours for normal PVR (top row) and high PVR (bottom row), and a single proximal MBTS (left column), a proximal MBTS and a 2 mm PDA (middle column), and a proximal MBTS and a 4 mm PDA (right column). (b) Spatially averaged WSS in the PA segment (top plot) and in the MBTS (bottom plot). See Fig. 5 caption for more details.

Grahic Jump Location
Fig. 10

(a) WSSG contours for normal PVR (top row) and high PVR (bottom row), and a single proximal MBTS (left column), a proximal MBTS and a 2 mm PDA (middle column), and a proximal MBTS and a 4 mm PDA (right column). (b) Spatially averaged WSSG in the PA segment (top plot) and in the MBTS (bottom plot). See Fig. 5 caption for more details.

Grahic Jump Location
Fig. 11

(a) OSI contours for normal PVR (top row) and high PVR (bottom row), and a single proximal MBTS (left column), a proximal MBTS and a 2 mm PDA (middle column), and a proximal MBTS and a 4 mm PDA (right column). (b) Spatially averaged OSI in the PA segment (top plot) and in the MBTS (bottom plot). See Fig. 5 caption for more details.




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In