Research Papers

A Zero-Dimensional Model and Protocol for Simulating Patient-Specific Pulmonary Hemodynamics From Limited Clinical Data

[+] Author and Article Information
Vitaly O. Kheyfets

University of Colorado Anschutz Medical Campus,
Children's Hospital Colorado,
Aurora, CO 80045
e-mail: vitaly.kheyfets@ucdenver.edu

Jamie Dunning

University of Colorado Anschutz Medical Campus,
Children's Hospital Colorado,
Aurora, CO 80045
e-mail: jamie.dunning@ucdenver.edu

Uyen Truong

University of Colorado Anschutz Medical Campus,
Children's Hospital Colorado,
Aurora, CO 80045
e-mail: Uyen.Truong@childrenscolorado.org

Dunbar Ivy

University of Colorado Anschutz Medical Campus,
Children's Hospital Colorado,
Aurora, CO 80045
e-mail: Dunbar.Ivy@childrenscolorado.org

Kendall Hunter

University of Colorado Anschutz Medical Campus,
Children's Hospital Colorado,
Aurora, CO 80045
e-mail: Kendall.hunter@ucdenver.edu

Robin Shandas

University of Colorado Anschutz Medical Campus,
Children's Hospital Colorado,
Aurora, CO 80045
e-mail: robin.shandas@ucdenver.edu

Manuscript received December 5, 2014; final manuscript received September 16, 2016; published online November 3, 2016. Assoc. Editor: Naomi Chesler.

J Biomech Eng 138(12), 121001 (Nov 03, 2016) (8 pages) Paper No: BIO-14-1610; doi: 10.1115/1.4034830 History: Received December 05, 2014; Revised September 16, 2016

In pulmonary hypertension (PH) diagnosis and management, many useful functional markers have been proposed that are unfeasible for clinical implementation. For example, assessing right ventricular (RV) contractile response to a gradual increase in pulmonary arterial (PA) impedance requires simultaneously recording RV pressure and volume, and under different afterload/preload conditions. In addition to clinical applications, many research projects are hampered by limited retrospective clinical data and could greatly benefit from simulations that extrapolate unavailable hemodynamics. The objective of this study was to develop and validate a 0D computational model, along with a numerical implementation protocol, of the RV–PA axis. Model results are qualitatively compared with published clinical data and quantitatively validated against right heart catheterization (RHC) for 115 pediatric PH patients. The RV–PA circuit is represented using a general elastance function for the RV and a three-element Windkessel initial value problem for the PA. The circuit mathematically sits between two reservoirs of constant pressure, which represent the right and left atriums. We compared Pmax, Pmin, mPAP, cardiac output (CO), and stroke volume (SV) between the model and RHC. The model predicted between 96% and 98% of the variability in pressure and 98–99% in volumetric characteristics (CO and SV). However, Bland Altman plots showed the model to have a consistent bias for most pressure and volumetric parameters, and differences between model and RHC to have considerable error. Future studies will address this issue and compare specific waveforms, but these initial results are extremely promising as preliminary proof of concept of the modeling approach.

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


Kass, D. A. , and Kelly, R. P. , 1992, “ Ventriculo-Arterial Coupling: Concepts, Assumptions, and Applications,” Ann. Biomed. Eng., 20(1), pp. 41–62. [CrossRef] [PubMed]
Wang, Z. , and Chesler, N. C. , 2011, “ Pulmonary Vascular Wall Stiffness: An Important Contributor to the Increased Right Ventricular Afterload With Pulmonary Hypertension,” Pulm. Circ., 1(2), pp. 212–223. [CrossRef] [PubMed]
Fourie, P. R. , Coetzee, A. R. , and Bolliger, C. T. , 1992, “ Pulmonary Artery Compliance: Its Role in Right Ventricular-Arterial Coupling,” Cardiovasc. Res., 26(9), pp. 839–844. [CrossRef] [PubMed]
Wesseling, K. H. , Jansen, J. R. , Settels, J. J. , and Schreuder, J. J. , 1993, “ Computation of Aortic Flow From Pressure in Humans Using a Nonlinear, Three-Element Model,” J. Appl. Physiol., 74(5), pp. 2566–2573. [PubMed]
Lankhaar, J. W. , Rovekamp, F. A. , Steendijk, P. , Faes, T. J. , Westerhof, B. E. , Kind, T. , Vonk-Noordegraaf, A. , and Westerhof, N. , 2009, “ Modeling the Instantaneous Pressure-Volume Relation of the Left Ventricle: A Comparison of Six Models,” Ann. Biomed. Eng., 37(9), pp. 1710–1726. [CrossRef] [PubMed]
Kind, T. , Faes, T. J. , Vonk-Noordegraaf, A. , and Westerhof, N. , 2011, “ Proportional Relations Between Systolic, Diastolic and Mean Pulmonary Artery Pressure are Explained by Vascular Properties,” Cardiovasc. Eng. Technol., 2(1), pp. 15–23. [CrossRef] [PubMed]
Lanzarone, E. , and Ruggeri, F. , 2013, “ Inertance Estimation in a Lumped-Parameter Hydraulic Simulator of Human Circulation,” ASME J. Biomech. Eng., 135(6), pp. 61012–61017. [CrossRef]
Lankhaar, J. W. , Westerhof, N. , Faes, T. J. , Marques, K. M. , Marcus, J. T. , Postmus, P. E. , and Vonk-Noordegraaf, A. , 2006, “ Quantification of Right Ventricular Afterload in Patients With and Without Pulmonary Hypertension,” Am. J. Physiol.: Heart Circ. Physiol., 291(4), pp. H1731–1737. [CrossRef] [PubMed]
Burkhoff, D. , Alexander, J., Jr. , and Schipke, J. , 1988, “ Assessment of Windkessel as a Model of Aortic Input Impedance,” Am. J. Physiol., 255(4 Pt 2), pp. H742–753. [PubMed]
Nichols, W. W. , Nichols, W. W. , and McDonald, D. A. , 2011, McDonald's Blood Flow in Arteries: Theoretic, Experimental, and Clinical Principles, Hodder Arnold, London.
Zamir, M. , 2005, The Physics of Coronary Blood Flow, Springer, New York.
Tian, L. , Hunter, K. S. , Kirby, K. S. , Ivy, D. D. , and Shandas, R. , 2010, “ Measurement Uncertainty in Pulmonary Vascular Input Impedance and Characteristic Impedance Estimated From Pulsed-Wave Doppler Ultrasound and Pressure: Clinical Studies on 57 Pediatric Patients,” Physiol. Meas., 31(6), pp. 729–748. [CrossRef] [PubMed]
Reiter, G. , Reiter, U. , Kovacs, G. , Kainz, B. , Schmidt, K. , Maier, R. , Olschewski, H. , and Rienmueller, R. , 2008, “ Magnetic Resonance-Derived 3-Dimensional Blood Flow Patterns in the Main Pulmonary Artery as a Marker of Pulmonary Hypertension and a Measure of Elevated Mean Pulmonary Arterial Pressure,” Circ. Cardiovasc. Imaging, 1(1), pp. 23–30. [CrossRef] [PubMed]
Bachler, P. , Pinochet, N. , Sotelo, J. , Crelier, G. , Irarrazaval, P. , Tejos, C. , and Uribe, S. , 2013, “ Assessment of Normal Flow Patterns in the Pulmonary Circulation by Using 4D Magnetic Resonance Velocity Mapping,” Magn. Reson. Imaging, 31(2), pp. 178–188. [CrossRef] [PubMed]
Parlikar, T. A. , 2007, “ Modeling and Monitoring of Cardiovascular Dynamics for Patients in Critical Care,” Ph.D. thesis, Massachusetts Institute of Technology, Boston, MA.
Kung, E. , and Taylor, C. , 2011, “ Development of a Physical Windkessel Module to Re-Create in vivo Vascular Flow Impedance for in vitro Experiments,” Cardiovasc. Eng. Technol., 2(1), pp. 2–14. [CrossRef] [PubMed]
Stevenson, D. , Revie, J. , Chase, J. G. , Hann, C. E. , Shaw, G. M. , Lambermont, B. , Ghuysen, A. , Kolh, P. , and Desaive, T. , 2012, “ Beat-to-Beat Estimation of the Continuous Left and Right Cardiac Elastance From Metrics Commonly Available in Clinical Settings,” Biomed. Eng. Online, 11(1), p. 73. [CrossRef] [PubMed]
Haddad, F. , Hunt, S. A. , Rosenthal, D. N. , and Murphy, D. J. , 2008, “ Right Ventricular Function in Cardiovascular Disease, Part I: Anatomy, Physiology, Aging, and Functional Assessment of the Right Ventricle,” Circulation, 117(11), pp. 1436–1448. [CrossRef] [PubMed]
Bellofiore, A. , and Chesler, N. C. , 2013, “ Methods for Measuring Right Ventricular Function and Hemodynamic Coupling With the Pulmonary Vasculature,” Ann. Biomed. Eng., 41(7), pp. 1384–1398. [CrossRef] [PubMed]
Mukkamala, R. , 2000, “ A Forward Model-Based Analysis of Cardiovascular System Identification Methods,” Ph.D. thesis, Department of Electrical Engineering and Computer Science, MIT, Cambridge, MA.
Marak, K. P. , 2013, “ The Time Varying Elastance Model Used as a Boundary Condition in Arterial Network Simulations,” Master's thesis, Norwegian University of Science and Technology, Trondheim, Norway.
Chuang, P. P. , Wilson, R. F. , Homans, D. C. , Stone, K. , Bergman, T. , Bennett, T. D. , and Kubo, S. H. , 1996, “ Measurement of Pulmonary Artery Diastolic Pressure From a Right Ventricular Pressure Transducer in Patients With Heart Failure,” J. Card. Failure, 2(1), pp. 41–46. [CrossRef]
Graham, T. P., Jr. , Jarmakani, J. M. , Atwood, G. F. , and Canent, R. V., Jr. , 1973, “ Right Ventricular Volume Determinations in Children. Normal Values and Observations With Volume or Pressure Overload,” Circulation, 47(1), pp. 144–153. [CrossRef] [PubMed]
Horsfield, K. , and Woldenberg, M. J. , 1989, “ Diameters and Cross-Sectional Areas of Branches in the Human Pulmonary Arterial Tree,” Anat. Rec., 223(3), pp. 245–251. [CrossRef] [PubMed]
Akay, H. O. , Ozmen, C. A. , Bayrak, A. H. , Senturk, S. , Katar, S. , Nazaroglu, H. , and Taskesen, M. , 2009, “ Diameters of Normal Thoracic Vascular Structures in Pediatric Patients,” Surg. Radiol. Anat.: SRA, 31(10), pp. 801–807. [CrossRef]
Westerhof, N. , Stergiopulos, N. , and Noble, M. I. , 2010, Snapshots of Hemodynamics—An Aid for Clinical Research and Graduate Education, Springer, New York.
Muthurangu, V. , Atkinson, D. , Sermesant, M. , Miquel, M. E. , Hegde, S. , Johnson, R. , Andriantsimiavona, R. , Taylor, A. M. , Baker, E. , Tulloh, R. , Hill, D. , and Razavi, R. S. , 2005, “ Measurement of Total Pulmonary Arterial Compliance Using Invasive Pressure Monitoring and MR Flow Quantification During MR-Guided Cardiac Catheterization,” Am. J. Physiol.: Heart Circ. Physiol., 289(3), pp. H1301–1306. [CrossRef] [PubMed]
Taylor, C. A. , and Figueroa, C. A. , 2009, “ Patient-Specific Modeling of Cardiovascular Mechanics,” Annu. Rev. Biomed. Eng., 11(1), pp. 109–134. [CrossRef] [PubMed]
Hunter, K. S. , Lee, P. F. , Lanning, C. J. , Ivy, D. D. , Kirby, K. S. , Claussen, L. R. , Chan, K. C. , and Shandas, R. , 2008, “ Pulmonary Vascular Input Impedance is a Combined Measure of Pulmonary Vascular Resistance and Stiffness and Predicts Clinical Outcomes Better Than Pulmonary Vascular Resistance Alone in Pediatric Patients With Pulmonary Hypertension,” Am. Heart J., 155(1), pp. 166–174. [CrossRef] [PubMed]
Spilker, R. L. , Feinstein, J. A. , Parker, D. W. , Reddy, V . M. , and Taylor, C. A. , 2007, “ Morphometry-Based Impedance Boundary Conditions for Patient-Specific Modeling of Blood Flow in Pulmonary Arteries,” Ann. Biomed. Eng., 35(4), pp. 546–559. [CrossRef] [PubMed]
Brimioulle, S. , Wauthy, P. , Ewalenko, P. , Rondelet, B. , Vermeulen, F. , Kerbaul, F. , and Naeije, R. , 2003, “ Single-Beat Estimation of Right Ventricular End-Systolic Pressure-Volume Relationship,” Am. J. Physiol.: Heart Circ. Physiol., 284(5), pp. H1625–1630. [CrossRef] [PubMed]
Aguero, J. , Ishikawa, K. , Hadri, L. , Santos-Gallego, C. , Fish, K. , Hammoudi, N. , Chaanine, A. , Torquato, S. , Naim, C. , Ibanez, B. , Pereda, D. , Garcia-Alvarez, A. , Fuster, V. , Sengupta, P. P. , Leopold, J. A. , and Hajjar, R. J. , 2014, “ Characterization of Right Ventricular Remodeling and Failure in a Chronic Pulmonary Hypertension Model,” Am. J. Physiol.: Heart Circ. Physiol., 307(8), pp. H1204–1215. [CrossRef] [PubMed]
Champion, H. C. , Michelakis, E. D. , and Hassoun, P. M. , 2009, “ Comprehensive Invasive and Noninvasive Approach to the Right Ventricle-Pulmonary Circulation Unit: State of the Art and Clinical and Research Implications,” Circulation, 120(11), pp. 992–1007. [CrossRef] [PubMed]
Gehalot, P. , Zhang, R. , Mathew, A. , and Behbehani, K. , 2006, “ Efficacy of Using Mean Arterial Blood Pressure Sequence for Three-Element Windkessel Model Estimation,” 28th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS '06), New York, Aug. 30-Sept. 3, pp. 1379–1382.
Olufsen, M. S. , Peskin, C. S. , Kim, W. Y. , Pedersen, E. M. , Nadim, A. , and Larsen, J. , 2000, “ Numerical Simulation and Experimental Validation of Blood Flow in Arteries With Structured-Tree Outflow Conditions,” Ann. Biomed. Eng., 28(11), pp. 1281–1299. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 1

Physical model with circuit diagram of the simulated RV–PA axis. In parts of this manuscript, R = PVR − Zc.

Grahic Jump Location
Fig. 2

Numerical protocol for simulating RV–PA hemodynamics. EDV = end diastolic volume; EDPVR = end diastolic pressure–volume relationship.

Grahic Jump Location
Fig. 3

(a) Typical RV and PA pressure waveforms computed using RV–PA axis model. (b) and (c) Ventricular volume and ventricular pressure–volume loop, respectively. (d) Pulmonary vascular impedance in the frequency domain, computed using simulated PA pressure and flow waveforms as outlined in Ref. [2]. Note: PPAand PRVare pulmonary and RV pressure, respectively.

Grahic Jump Location
Fig. 4

Max (left column), min (middle column), and mean (right column) PA pressure comparison between RV and PA axis model and measured RHC hemodynamics. In each column, the top row shows a correlation between measured and simulated values. The slope (m) and y-intercept (b) are coefficients for the fitted line: CVSIM=m⋅RHC+b. The bottom rows show Bland Altman plots, where the middle and outer lines represent the consistent bias and 1.96SD, respectively.

Grahic Jump Location
Fig. 5

(Left) cardiac output (CO) computed by integrating the flow waveform simulated with the RV–PA model, compared with CO measured using a thermodilution catheter. (Right) stroke volume (SV) computed according the difference between the maximum and minimum volume measured by implementing the elastance function, compared with SV = CO/HR. The slope (m) and y-intercept (b) are coefficients for the fitted line: CVSIM=m⋅RHC+b. The bottom rows show Bland Altman plots, where the middle and outer lines represent the consistent bias and 1.96SD, respectively.



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