Particle swarm optimizer for arterial blood flow models

BACKGROUND AND OBJECTIVE Mathematical modeling and computational simulations of arterial blood flow network can offer an insilico platform for both diagnostics and therapeutic phases of patients that suffer from cardiac diseases. These models are normally complex and involve many unknown parameters. For physiological relevance, these parameters should be optimized using in-vivo human/animal data sets. The main goal of this work is to develop an efficient, yet an accurate optimization algorithm to compute parameters in the arterial blood flow models. METHODS The particle swarm optimization (PSO) method is proposed herein for the first time, as an accurate algorithm that applies to computing parameters in the Windkessel type model of blood flow in the arterial system. We begin by defining a 6-element Windkessel (WK6) arterial flow model, which is then implemented and validated using multiple flow rate and aortic pressure measurements obtained from different subjects including dogs, pigs and humans. The parameters in the model are obtained using the PSO technique which minimizes the pressure root mean square (P-RMS) error between the computed and the measured aortic pressure waveform. RESULTS Model parameters obtained using the proposed PSO method were able to recover the pressure waveform in the aorta during the cardiac cycle for both healthy and diseased species (animals/humans). The PSO method provides an accurate approach to solve this challenging multi-dimensional parameter identification problem. The results obtained by PSO algorithm was compared with the classical gradient-based, namely the non-linear square fit (NLSF) algorithm. CONCLUSIONS The results indicate that the PSO method offers alternative and accurate method to find optimal physiological parameters involved in the Windkessel model for the study of arterial blood flow network. The PSO method has performed better than the NLSF approach as depicted from the P-RMS calculations. Finally, we believe that the PSO method offers a great potential and could be used for many other biomedicine optimization problems.

[1]  Y. Aboelkassem,et al.  A hybrid Windkessel-Womersley model for blood flow in arteries. , 2019, Journal of theoretical biology.

[2]  Xiangdong Wang,et al.  Parameters identification of chaotic systems via chaotic ant swarm , 2006 .

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

[4]  J B Charles,et al.  Mathematical modeling of human cardiovascular system for simulation of orthostatic response. , 1992, The American journal of physiology.

[5]  T. G. Coleman,et al.  Circulation: overall regulation. , 1972, Annual review of physiology.

[6]  Karim Azer,et al.  A One-dimensional Model of Blood Flow in Arteries with Friction and Convection Based on the Womersley Velocity Profile , 2007, Cardiovascular engineering.

[7]  Luca Ridolfi,et al.  Impact of atrial fibrillation on the cardiovascular system through a lumped-parameter approach , 2014, Medical & Biological Engineering & Computing.

[8]  N. Westerhof,et al.  Modeling pressure-flow relations in cardiac muscle in diastole and systole. , 1997, The American journal of physiology.

[9]  Anselmo Cardoso de Paiva,et al.  Convolutional neural network-based PSO for lung nodule false positive reduction on CT images , 2018, Comput. Methods Programs Biomed..

[10]  Roberto Burattini,et al.  Development of systemic arterial mechanical properties from infancy to adulthood interpreted by four-element windkessel models. , 2007, Journal of applied physiology.

[11]  Bram W. Smith,et al.  Velocity profile method for time varying resistance in minimal cardiovascular system models. , 2003, Physics in medicine and biology.

[12]  M. Ursino,et al.  Role of short-term cardiovascular regulation in heart period variability: a modeling study. , 2003, American journal of physiology. Heart and circulatory physiology.

[13]  Eun Bo Shim,et al.  Mathematical modeling of cardiovascular system dynamics using a lumped parameter method. , 2004, The Japanese journal of physiology.

[14]  John W. Clark,et al.  A Closed-Loop Model of the Canine Cardiovascular System That Includes Ventricular Interaction , 2000, Comput. Biomed. Res..

[15]  A. Quarteroni,et al.  On the coupling of 3D and 1D Navier-Stokes equations for flow problems in compliant vessels , 2001 .

[16]  Yubing Shi,et al.  Review of Zero-D and 1-D Models of Blood Flow in the Cardiovascular System , 2011, Biomedical engineering online.

[17]  S Scarsoglio,et al.  Impaired coronary blood flow at higher heart rates during atrial fibrillation: Investigation via multiscale modelling , 2019, Comput. Methods Programs Biomed..

[18]  R. Burattini,et al.  Complex and frequency-dependent compliance of viscoelastic windkessel resolves contradictions in elastic windkessels. , 1998, Medical engineering & physics.

[19]  Qing Lv,et al.  A novel heterogeneous feature ant colony optimization and its application on robot path planning , 2015, 2015 IEEE Congress on Evolutionary Computation (CEC).

[20]  A P Avolio,et al.  Pulsatile Flow and Pressure in Human Systemic Arteries: Studies in Man and in a Multibranched Model of the Human Systemic Arterial Tree , 1980, Circulation research.

[21]  Modeling of aortic valve dynamics in a lumped parameter model of left ventricular-arterial coupling , 2008 .

[22]  Padmavathi Kora,et al.  ECG based Myocardial Infarction detection using Hybrid Firefly Algorithm , 2017, Comput. Methods Programs Biomed..

[23]  J R Jansen,et al.  Computation of aortic flow from pressure in humans using a nonlinear, three-element model. , 1993, Journal of applied physiology.

[24]  R Beyar,et al.  Interaction between cardiac chambers and thoracic pressure in intact circulation. , 1987, The American journal of physiology.

[25]  L. Paradowski,et al.  Characterization of pulmonary arterial input impedance with lumped parameter models. , 1987, The American journal of physiology.

[26]  Saïd Hanafi,et al.  A hybrid quantum particle swarm optimization for the Multidimensional Knapsack Problem , 2016, Eng. Appl. Artif. Intell..

[27]  J Geoffrey Chase,et al.  Minimal haemodynamic system model including ventricular interaction and valve dynamics. , 2004, Medical engineering & physics.

[28]  James Kennedy,et al.  Particle swarm optimization , 2002, Proceedings of ICNN'95 - International Conference on Neural Networks.

[29]  Nigel G Shrive,et al.  Time-domain representation of ventricular-arterial coupling as a windkessel and wave system. , 2003, American journal of physiology. Heart and circulatory physiology.

[30]  Carlos A. Coello Coello,et al.  Constraint-handling in nature-inspired numerical optimization: Past, present and future , 2011, Swarm Evol. Comput..

[31]  Tao Li,et al.  Particle swarm optimizer with crossover operation , 2018, Eng. Appl. Artif. Intell..

[32]  Ying Han,et al.  Clinical flexible needle puncture path planning based on particle swarm optimization , 2020, Comput. Methods Programs Biomed..

[33]  T G Coleman,et al.  HUMAN. A comprehensive physiological model. , 1983, The Physiologist.

[34]  S. Campbell,et al.  Mathematical modeling of aortic valve dynamics during systole. , 2015, Journal of theoretical biology.

[35]  H. H. Hardy,et al.  A digital computer model of the human circulatory system , 1982, Medical and Biological Engineering and Computing.

[36]  J. Clark,et al.  A dynamic model of ventricular interaction and pericardial influence. , 1997, The American journal of physiology.

[37]  N Westerhof,et al.  Total arterial inertance as the fourth element of the windkessel model. , 1999, The American journal of physiology.

[38]  D. Burkhoff,et al.  Assessment of Windkessel as a model of aortic input impedance. , 1988, The American journal of physiology.

[39]  A. Quarteroni,et al.  Analysis of lumped parameter models for blood flow simulations and their relation with 1D models , 2004 .

[40]  R. Mark,et al.  Computational modeling of cardiovascular response to orthostatic stress. , 2002, Journal of applied physiology.

[41]  M. Ursino,et al.  A mathematical model of the carotid baroregulation in pulsating conditions , 1999, IEEE Transactions on Biomedical Engineering.

[42]  P. Kolh,et al.  Systemic and pulmonary hemodynamics assessed with a lumped-parameter heart-arterial interaction model , 2003 .

[43]  L. Ridolfi,et al.  Fluid dynamics of heart valves during atrial fibrillation: a lumped parameter-based approach , 2015, Computer methods in biomechanics and biomedical engineering.

[44]  S. Chiaramida,et al.  A comprehensive model for right-left heart interaction under the influence of pericardium and baroreflex. , 1997, The American journal of physiology.

[45]  T. G. Coleman,et al.  Quantitative Circulatory Physiology: An integrative mathematical model of human physiology for medical education , 2007, Advances in physiology education.

[46]  J Beschi Raja,et al.  PSO-FCM based data mining model to predict diabetic disease , 2020, Comput. Methods Programs Biomed..

[47]  A. Capello,et al.  Tracking time-varying properties of the systemic vascular bed , 1989, IEEE Transactions on Biomedical Engineering.

[48]  R. Burattini,et al.  Computer identification of models for the arterial tree input impedance: Comparison between two new simple models and first experimental results , 1982, Medical and Biological Engineering and Computing.

[49]  Martin Hexamer,et al.  Simulation and prediction of cardiotherapeutical phenomena from a pulsatile model coupled to the Guyton circulatory model , 2002, IEEE Transactions on Biomedical Engineering.

[50]  Aleksander S. Popel,et al.  Modeling the circulation with three-terminal electrical networks containing special nonlinear capacitors , 2006, Annals of Biomedical Engineering.