Particle swarm optimizer for arterial blood flow models

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.

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.

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.

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.