Thrombus Clotting Simulation Method Based on the Gillespie Method

Clotting of thrombus is a reflection of blood physiological properties and it plays an important role in virtual surgery and physiological simulation of a human body.In this paper,we propose a thrombus simulation method based on the Gillespie method.Firstly,the blood is modeled in view of thrombogenesis and physiology reaction rate are set according to reaction components.Then we determine the reactions and relevant time durations among particles based on the Gillespie method’s Monte Carlo steps within the kernel boundary before fluid particles are transformed to thrombus particles.At last,we fix the position of fibrous particles generated from the previous step and thrombus is clotted.In order to show the result in real time,CUDA parallel computing toolkit was utilized to speed up the computation.We analyzed thrombus clotting,clotting speed,and runtime efficiency through experiments.The results show that our method can simulate the clotting of thrombus efficiently,which avoids the high computational price of the Euler method.Meanwhile,the simulated thrombus would restrict the blood flow path,which shows us the coupling effect between blood and thrombus.We also rendered some realistic thrombus clotting pictures.