Pulmonary Hemodynamics Simulations Before Stage 2 Single Ventricle Surgery: Patient-Specific Parameter Identification and Clinical Data Assessment

Abstract Single ventricle heart defects involve pathologies in which the heart has only one functional pumping chamber. In these conditions, treatment consists of three staged procedures. At stage 1 pulmonary flow is provided through an artificial shunt from the systemic circulation. Representative hemodynamics models able to explore different virtual surgical options can be built based on pre-operative imaging and patient data. In this context, the specification of boundary conditions is necessary to compute pressure and flow in the entire domain. However, these boundary conditions are rarely the measured variables. Moreover, to take into account the rest of the circulation outside of the three-dimensional modeled domain, a number of reduced order models exist. A simplified method is presented to iteratively, but automatically, tune reduced model parameters from hemodynamic data clinically measured before stage 2 surgery. Patient-specific local hemodynamics around the distal systemic-to-pulmonary shunt anastomosis and the connected pulmonary arteries are also analyzed. Multi-scale models of pre-stage 2 single ventricle patients are developed, including a 3D model of shunt-pulmonary connection and a number of pulmonary arteries. For each pulmonary outlet a total downstream resistance is identified, consistent with measured flow split and pressures. Target pressures such as minimum, maximum or average over one or both lungs are considered, depending on the clinical measurement. When possible, both steady and pulsatile identifications are performed. The methodology is demonstrated with six patient-specific models: the clinical target data are well-matched, except for one case where clinical data were subsequently found inconsistent. Inhomogeneous pressure, swirling blood flow patterns and very high wall shear stress 3D maps highlight similarities and differences among patients. Steady and pulsatile tuning results are similar. This work demonstrates (1) how to use routine clinical data to define boundary conditions for patient-specific 3D models in pre-stage 2 single ventricle circulations and (2) how simulations can help to check the coherence of clinical data, or provide insights to clinicians that are otherwise difficult to measure, such as in the presence of kinks. Finally, the choice of steady vs. pulsatile tuning, limitations and possible extensions of this work are discussed.

[1]  N. Stergiopulos,et al.  Three- and four-element Windkessel models: Assessment of their fitting performance in a large cohort of healthy middle-aged individuals , 2008, Proceedings of the Institution of Mechanical Engineers. Part H, Journal of engineering in medicine.

[2]  A. Marsden,et al.  A Primer on Computational Simulation in Congenital Heart Disease for the Clinician , 2010, 1101.3726.

[3]  C A Dawson,et al.  Anatomic distribution of pulmonary vascular compliance. , 1998, Journal of applied physiology.

[4]  F. Liang,et al.  Assessment of cardiovascular function by combining clinical data with a computational model of the cardiovascular system. , 2013, The Journal of thoracic and cardiovascular surgery.

[5]  A B DuBois,et al.  Longitudinal distribution of vascular resistance in the pulmonary arteries, capillaries, and veins. , 1968, The Journal of clinical investigation.

[6]  J-F Gerbeau,et al.  A methodological paradigm for patient‐specific multi‐scale CFD simulations: from clinical measurements to parameter estimates for individual analysis , 2014, International journal for numerical methods in biomedical engineering.

[7]  Charles A. Taylor,et al.  Outflow boundary conditions for three-dimensional finite element modeling of blood flow and pressure in arteries , 2006 .

[8]  N. Stergiopulos,et al.  Total arterial inertance as the fourth element of the windkessel model. , 1999, American journal of physiology. Heart and circulatory physiology.

[9]  C. T. Kelley,et al.  Estimation and identification of parameters in a lumped cerebrovascular model. , 2008, Mathematical biosciences and engineering : MBE.

[10]  A. Kassab,et al.  Computational analysis of hybrid Norwood circulation with distal aortic arch obstruction and reverse Blalock-Taussig shunt. , 2012, The Annals of thoracic surgery.

[11]  R. Pettigrew,et al.  Turbulent Fluctuation Velocity: The Most Significant Determinant of Signal Loss in Stenotic Vessels , 1995, Magnetic resonance in medicine.

[12]  I. Vignon-Clementel,et al.  Three-dimensional simulations in Glenn patients: clinically based boundary conditions, hemodynamic results and sensitivity to input data. , 2011, Journal of biomechanical engineering.

[13]  Mahmoud Ismail,et al.  Adjoint-based inverse analysis of windkessel parameters for patient-specific vascular models , 2013, J. Comput. Phys..

[14]  A. Yoganathan,et al.  Imaging and patient-specific simulations for the Fontan surgery: current methodologies and clinical applications. , 2010, Progress in pediatric cardiology.

[15]  A. Marsden,et al.  Predictive modeling of the virtual Hemi-Fontan operation for second stage single ventricle palliation: two patient-specific cases. , 2013, Journal of biomechanics.

[16]  M. Umezu,et al.  Computational Hemodynamic Analysis in Congenital Heart Disease: Simulation of the Norwood Procedure , 2010, Annals of Biomedical Engineering.

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

[18]  Shu Takagi,et al.  Patient‐specific assessment of cardiovascular function by combination of clinical data and computational model with applications to patients undergoing Fontan operation , 2014, International journal for numerical methods in biomedical engineering.

[19]  C A Taylor,et al.  Outflow boundary conditions for 3D simulations of non-periodic blood flow and pressure fields in deformable arteries , 2010, Computer methods in biomechanics and biomedical engineering.

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

[21]  Giancarlo Pennati,et al.  Boundary conditions of patient-specific fluid dynamics modelling of cavopulmonary connections: possible adaptation of pulmonary resistances results in a critical issue for a virtual surgical planning , 2011, Interface Focus.

[22]  Christopher E. Hann,et al.  Unique parameter identification for cardiac diagnosis in critical care using minimal data sets , 2010, Comput. Methods Programs Biomed..

[23]  Sanjay Pant,et al.  A Multiscale Filtering-Based Parameter Estimation Method for Patient-Specific Coarctation Simulations in Rest and Exercise , 2013, STACOM.

[24]  F. Migliavacca,et al.  Computational fluid dynamics in the evaluation of hemodynamic performance of cavopulmonary connections after the Norwood procedure for hypoplastic left heart syndrome. , 2003, The Journal of thoracic and cardiovascular surgery.

[25]  A. Marsden,et al.  An integrated approach to patient-specific predictive modeling for single ventricle heart palliation , 2014, Computer methods in biomechanics and biomedical engineering.

[26]  Charles A. Taylor,et al.  Efficient anisotropic adaptive discretization of the cardiovascular system , 2006 .

[27]  Xiangrong Li,et al.  Anisotropic adaptive finite element method for modelling blood flow , 2005, Computer methods in biomechanics and biomedical engineering.

[28]  F. Migliavacca,et al.  Use of rapid prototyping models in the planning of percutaneous pulmonary valved stent implantation , 2007, Proceedings of the Institution of Mechanical Engineers. Part H, Journal of engineering in medicine.

[29]  Charles A. Taylor,et al.  Tuning Multidomain Hemodynamic Simulations to Match Physiological Measurements , 2010, Annals of Biomedical Engineering.

[30]  Charles A. Taylor,et al.  Morphometry-Based Impedance Boundary Conditions for Patient-Specific Modeling of Blood Flow in Pulmonary Arteries , 2007, Annals of Biomedical Engineering.

[31]  P. Moireau,et al.  Sequential parameter estimation for fluid–structure problems: Application to hemodynamics , 2012, International journal for numerical methods in biomedical engineering.

[32]  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.

[33]  Jean-Frédéric Gerbeau,et al.  Modified Navier-Stokes equations for the outflow boundary conditions in hemodynamics , 2016 .

[34]  C. O'leary,et al.  Perioperative Distribution of Pulmonary Vascular Resistance in Patients Undergoing Coronary Artery Surgery , 1996, Anesthesia and analgesia.

[35]  A. Kassab,et al.  Computational fluid dynamics in congenital heart disease* , 2012, Cardiology in the Young.

[36]  F. Migliavacca,et al.  Percutaneous pulmonary valve implantation based on rapid prototyping of right ventricular outflow tract and pulmonary trunk from MR data. , 2007, Radiology.