A parameter estimation framework for patient-specific hemodynamic computations

We propose a fully automated parameter estimation framework for performing patient-specific hemodynamic computations in arterial models. To determine the personalized values of the windkessel models, which are used as part of the geometrical multiscale circulation model, a parameter estimation problem is formulated. Clinical measurements of pressure and/or flow-rate are imposed as constraints to formulate a nonlinear system of equations, whose fixed point solution is sought. A key feature of the proposed method is a warm-start to the optimization procedure, with better initial solution for the nonlinear system of equations, to reduce the number of iterations needed for the calibration of the geometrical multiscale models. To achieve these goals, the initial solution, computed with a lumped parameter model, is adapted before solving the parameter estimation problem for the geometrical multiscale circulation model: the resistance and the compliance of the circulation model are estimated and compensated.The proposed framework is evaluated on a patient-specific aortic model, a full body arterial model, and multiple idealized anatomical models representing different arterial segments. For each case it leads to the best performance in terms of number of iterations required for the computational model to be in close agreement with the clinical measurements.

[1]  Thomas J. R. Hughes,et al.  In vivo validation of a one-dimensional finite-element method for predicting blood flow in cardiovascular bypass grafts , 2003, IEEE Transactions on Biomedical Engineering.

[2]  T. Wonnacott,et al.  Relation between diameter and flow in major branches of the arch of the aorta. , 1992, Journal of biomechanics.

[3]  L Kadem,et al.  Modeling the impact of concomitant aortic stenosis and coarctation of the aorta on left ventricular workload. , 2011, Journal of biomechanics.

[4]  Charles A. Taylor,et al.  Comparative study of viscoelastic arterial wall models in nonlinear one-dimensional finite element simulations of blood flow. , 2011, Journal of biomechanical engineering.

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

[6]  Jarek Rossignac,et al.  Simulating hemodynamics of the Fontan Y-graft based on patient-specific in vivo connections. , 2013, The Journal of thoracic and cardiovascular surgery.

[7]  N R Cholvin,et al.  Pressure Drop across Artificially Induced Stenoses in the Femoral Arteries of Dogs , 1975, Circulation research.

[8]  B.N. Steele Using one-dimensional finite element analysis to estimate differential pressure of renal artery stenoses , 2007, 2007 Computers in Cardiology.

[9]  M. Olufsen,et al.  Numerical Simulation and Experimental Validation of Blood Flow in Arteries with Structured-Tree Outflow Conditions , 2000, Annals of Biomedical Engineering.

[10]  M. Olufsen,et al.  Dynamics of cerebral blood flow regulation explained using a lumped parameter model. , 2002, American journal of physiology. Regulatory, integrative and comparative physiology.

[11]  F. Lazeyras,et al.  Validation of a patient-specific one-dimensional model of the systemic arterial tree. , 2011, American journal of physiology. Heart and circulatory physiology.

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

[13]  S. Sherwin,et al.  Pulse wave propagation in a model human arterial network: Assessment of 1-D visco-elastic simulations against in vitro measurements , 2011, Journal of biomechanics.

[14]  Puneet Sharma,et al.  A patient-specific reduced-order model for coronary circulation , 2012, 2012 9th IEEE International Symposium on Biomedical Imaging (ISBI).

[15]  Ghassan S. Kassab,et al.  A validated predictive model of coronary fractional flow reserve , 2012, Journal of The Royal Society Interface.

[16]  C. A. Figueroa,et al.  Sequential identification of boundary support parameters in a fluid-structure vascular model using patient image data , 2012, Biomechanics and Modeling in Mechanobiology.

[17]  D. Bessems,et al.  On the propagation of pressure and flow waves through the patient specific arterial system , 2003 .

[18]  P. Nithiarasu,et al.  A 1D arterial blood flow model incorporating ventricular pressure, aortic valve and regional coronary flow using the locally conservative Galerkin (LCG) method , 2008 .

[19]  Dorin Comaniciu,et al.  Non-Invasive Hemodynamic Assessment of Aortic Coarctation: Validation with In Vivo Measurements , 2013, Annals of Biomedical Engineering.

[20]  J. Womersley Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known , 1955, The Journal of physiology.

[21]  Sansuke M. Watanabe,et al.  Identification of vascular territory resistances in one-dimensional hemodynamics simulations. , 2012, Journal of biomechanics.

[22]  David A. Steinman,et al.  Image-Based Modeling of Blood Flow and Vessel Wall Dynamics: Applications, Methods and Future Directions , 2010, Annals of Biomedical Engineering.

[23]  G. Karniadakis,et al.  Outflow Boundary Conditions for Arterial Networks with Multiple Outlets , 2008, Annals of Biomedical Engineering.

[24]  K Low,et al.  An improved baseline model for a human arterial network to study the impact of aneurysms on pressure‐flow waveforms , 2012, International journal for numerical methods in biomedical engineering.

[25]  Pablo J. Blanco,et al.  A two-level time step technique for the partitioned solution of one-dimensional arterial networks , 2012 .

[26]  N. Stergiopulos,et al.  Simple and accurate way for estimating total and segmental arterial compliance: The pulse pressure method , 1994, Annals of Biomedical Engineering.

[27]  C. D. Murray THE PHYSIOLOGICAL PRINCIPLE OF MINIMUM WORK , 1931, The Journal of general physiology.

[28]  Mahmoud Ismail,et al.  CFD Challenge: Hemodynamic Simulation of a Patient-Specific Aortic Coarctation Model with Adjoint-Based Calibrated Windkessel Elements , 2012, STACOM.

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

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

[31]  F. Mut,et al.  Association of Hemodynamic Characteristics and Cerebral Aneurysm Rupture , 2011, American Journal of Neuroradiology.

[32]  J P Mynard,et al.  A simple, versatile valve model for use in lumped parameter and one‐dimensional cardiovascular models , 2012, International journal for numerical methods in biomedical engineering.

[33]  Berend E. Westerhof,et al.  The arterial Windkessel , 2009, Medical & Biological Engineering & Computing.

[34]  M. Chial,et al.  in simple , 2003 .

[35]  Charles A. Taylor,et al.  Computational simulations for aortic coarctation: representative results from a sampling of patients. , 2011, Journal of biomechanical engineering.

[36]  C D Murray,et al.  The Physiological Principle of Minimum Work: I. The Vascular System and the Cost of Blood Volume. , 1926, Proceedings of the National Academy of Sciences of the United States of America.

[37]  A. Quarteroni,et al.  One-dimensional models for blood flow in arteries , 2003 .

[38]  D. F. Young,et al.  Computer simulation of arterial flow with applications to arterial and aortic stenoses. , 1992, Journal of biomechanics.

[39]  Alfio Quarteroni,et al.  Computational vascular fluid dynamics: problems, models and methods , 2000 .