Department of Electronics and Telecommunications, Politecnico di Torino, Turin, Italy.
Department of Electrical & Computer Engineering, Institute of Biomedical Engineering, University of Toronto, Toronto, Canada.
Cardiovasc Eng Technol. 2023 Aug;14(4):505-525. doi: 10.1007/s13239-023-00669-z. Epub 2023 Jun 12.
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.
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.
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.
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.
选择适当的边界条件是血流模拟心血管模型开发的关键步骤。三元件 Windkessel 模型通常作为集中边界条件使用,提供外周循环的降阶表示。然而,Windkessel 参数的系统估计仍然是一个悬而未决的问题。此外,Windkessel 模型并不总是足以模拟血流动力学,后者通常需要更精细的边界条件。在这项研究中,我们提出了一种从截断点处的压力和流量波形估计高阶边界条件(包括 Windkessel 模型)参数的方法。此外,我们研究了采用更高阶边界条件对模型准确性的影响,这些边界条件对应于具有一个以上存储元件的等效电路。
所提出的技术基于时域矢量拟合,这是一种建模算法,给定系统输入和输出的样本,例如压力和流量波形,可以推导出一个近似它们关系的微分方程。
该方法的能力在由 55 个人体最大系统动脉组成的 1D 循环模型上进行了测试,以证明其准确性和在估计高于传统 Windkessel 模型的阶数的边界条件方面的有用性。将所提出的方法与其他常见的估计技术进行了比较,并在存在噪声数据和心理应激引起的主动脉流量生理变化的情况下验证了其在参数估计中的稳健性。
结果表明,所提出的方法能够准确估计任意阶数的边界条件。高阶边界条件可以提高心血管模拟的准确性,并且时域矢量拟合可以自动估计它们。