Computer simulation of vertebral artery occlusion in endovascular procedures

The aim of this work is to establish a computational pipeline for the simulation of blood flow in vasculatures and apply this pipeline to endovascular interventional scenarios, e.g. angioplasty in vertebral arteries.

A patient-specific supra-aortal vasculature is digitized from a 3D CT angiography image. By coupling a reduced formulation of the governing Navier-Stokes equations with a wall constitutive equation we are able to solve the transient flow in elastic vessels. By further incorporating a bifurcation model the blood flow across vascular branches can be evaluated, thus flow in a large vasculature can be modeled. Vascular diseases are simulated by modifying the arterial tree configurations, e.g. the effective diameters, schematic connectivity, etc. Occlusion in an artery is simulated by removing that artery from the arterial tree.

It takes about 2 min per cardiac cycle to compute blood flow in an arterial tree consisting of 38 vessels and 18 bifurcations on a laptop PC. The simulation results show that blood supply in the posterior region is compensated from the contralateral vertebral artery and the anterior cerebral arteries if one of the vertebral arteries is occluded.

The computational pipeline is computationally efficient and can capture main flow patterns at any point in the arterial tree. With further improvement it can serve as a powerful tool for the haemodynamic analysis in patient-specific vascular structures.