Abstract
Thrombosis and thromboembolism are deadly risk factors in blood-contacting biomedical devices, and in-silico models of thrombosis are attractive tools to understand the mechanics of these processes, though the simulation of thromboembolism remains underdeveloped. The purpose of this study is to modify an existing computational thrombosis model to allow for thromboembolism and to investigate the behavior of the modified model at a range of flow rates. The new and existing models are observed to lead to similar predictions of thrombosis in a canonical backward-facing step geometry across flow rates, and neither model predicts thrombosis in a turbulent flow. Simulations are performed by increasing flow rates in the case of a clot formed at lower flow to induce embolization. While embolization is observed, most of the clot breakdown is by shear rather than by breakup and subsequent transport of clotted material, and further work is required in the formulation and validation of embolization. This model provides a framework to further investigate thromboembolization.