Abstract
Despite recent advances in the development of computational methods of modeling thrombosis, relatively little effort has been made in developing methods of modeling blood clot embolization. Such a model would provide substantially greater understanding of the mechanics of embolization, as in-vitro and in-vivo characterization of embolization is difficult. Here, a method of computationally simulating embolization is developed. Experiments are performed of blood clots formed in a polycarbonate tube, where phosphate-buffered saline is run through the tube at increasing flow rates until the clot embolizes. The experiments revealed embolization can be initiated by leading edge and trailing edge detachment or by non-uniform detachment. Stress-relaxation experiments are also performed to establish values of constitutive parameters for subsequent simulations. The embolization in the tube is reproduced in silico using a multiphase volume-of-fluid approach, where the clot is modeled as viscoelastic. By varying the constitutive parameters at the wall, embolization can be reproduced in-silico at varying flow rates, and a range of constitutive parameters fitting the experiments is reported. Here, the leading edge embolization is simulated at flow rates consistent with the experiments demonstrating excellent agreement in this specific behavior.