Computational Simulation of Aneurysms Using Smoothed Particle Hydrodynamics

Modeling and simulation of aneurysm formation, growth, and rupture plays an essential role in a wide spectrum of application scenarios, ranging from risk stratification to stability prediction, and from clinical decision-making to treatment innovation. Unfortunately, it remains a non-trivial task due to the difficulties imposed by the complex and under-researched pathophysiological mechanisms behind the different development stages of various aneurysms. In this paper, we present a novel computational method for aneurysm simulation using smoothed particle hydrodynamics (SPH). Firstly, we consider blood in a vessel as a kind of incompressible fluid and model its flow dynamics using the SPH method; and then, to simulate aneurysm growth and rupture, the relationship between the aneurysm development and the properties of fluid particles is established by solving the motion control equation. In view of the prevalence of aneurysms in bifurcation vessels, we further enhance the capability of the model by introducing a solution for bifurcation aneurysms simulation according to Murray’s law. In addition, a CUDA parallel computing scheme is also designed to speed up the simulation process. To evaluate the performance of the proposed method, we conduct extensive experiments with different physical parameters associated with morphological characteristics of an aneurysm. The experimental results demonstrate the effectiveness and efficiency of proposed method in modeling and simulating aneurysm formation, growth, and rupture.