Abstract
Abstract
Purpose
The choice of appropriate boundary conditions is a crucial step in the development of cardiovascular models for blood flow simulations. The three-element Windkessel model is usually employed as a lumped boundary condition, providing a reduced order representation of the peripheral circulation. However, the systematic estimation of the Windkessel parameters remains an open problem. Moreover, the Windkessel model is not always adequate to model blood flow dynamics, which often require more elaborate boundary conditions. In this study, we propose a method for the estimation of the parameters of high order boundary conditions, including the Windkessel model, from pressure and flow rate waveforms at the truncation point. Moreover, we investigate the effect of adopting higher order boundary conditions, corresponding to equivalent circuits with more than one storage element, on the accuracy of the model.
Method
The proposed technique is based on Time-Domain Vector Fitting, a modeling algorithm that, given samples of the input and output of a system, such as pressure and flow waveforms, can derive a differential equation approximating their relation.
Results
The capabilities of the proposed method are tested on a 1D circulation model consisting of the 55 largest human systemic arteries, to demonstrate its accuracy and its usefulness to estimate boundary conditions with order higher than the traditional Windkessel models. The proposed method is compared to other common estimation techniques, and its robustness in parameter estimation is verified in presence of noisy data and of physiological changes of aortic flow rate induced by mental stress.
Conclusion
Results suggest that the proposed method is able to accurately estimate boundary conditions of arbitrary order. Higher order boundary conditions can improve the accuracy of cardiovascular simulations, and Time-Domain Vector Fitting can automatically estimate them.
Publisher
Springer Science and Business Media LLC
Subject
Cardiology and Cardiovascular Medicine,Biomedical Engineering
Reference55 articles.
1. Formaggia, L., A. Quarteroni, and A. Veneziani. Cardiovascular Mathematics: Modeling and Simulation of the Circulatory System, Vol. 1, Milano: Springer, 2010.
2. Pirola, S., Z. Cheng, O. Jarral, D. O’Regan, J. Pepper, T. Athanasiou, and X. Xu. On the choice of outlet boundary conditions for patient-specific analysis of aortic flow using computational fluid dynamics. J. Biomech. 60:15–21, 2017.
3. Morbiducci, U., R. Ponzini, D. Gallo, C. Bignardi, and G. Rizzo. Inflow boundary conditions for image-based computational hemodynamics: impact of idealized versus measured velocity profiles in the human aorta. J. Biomech. 46(1):102–109, 2013.
4. Shi, Y., P. Lawford, and R. Hose. Review of zero-d and 1-d models of blood flow in the cardiovascular system. Biomed. Eng. Online. 10(1):1–38, 2011.
5. Parker, K.H., and C. Jones. Forward and backward running waves in the arteries: analysis using the method of characteristics. 1990