LAUSR

Over the past decade, much of the development of computational models of device-related thrombosis has focused on platelet activity. While those models have been successful in predicting thrombus formation in medical devices operating at high shear rates (> 5000 s−1), they cannot be directly applied to low-shear devices, such as blood oxygenators and catheters, where emerging information suggest that fibrin formation is the predominant mechanism of clotting and platelet activity plays a secondary role. In the current work, we augment an existing platelet-based model of thrombosis with a partial model of the coagulation cascade that includes contact activation of factor XII and fibrin production. To calibrate the model, we simulate a backward-facing-step flow channel that has been extensively characterized in-vitro. Next, we perform blood perfusion experiments through a microfluidic chamber mimicking a hollow fiber membrane oxygenator and validate the model against these observations. The simulation results closely match the time evolution of the thrombus height and length in the backward-facing-step experiment. Application of the model to the microfluidic hollow fiber bundle chamber capture both gross features such as the increasing clotting trend towards the outlet of the chamber, as well as finer local features such as the structure of fibrin around individual hollow fibers. Our results are in line with recent findings that suggest fibrin production, through contact activation of factor XII, drives the thrombus formation in medical devices operating at low shear rates with large surface area to volume ratios. Author summary Patients treated with blood-contacting medical devices suffer from clotting complications. Over the past decades, a great effort has been made to develop computational tools to predict and prevent clot formation in these devices. However, most models have focused on platelet activity and neglected other important parts of the problem such as the coagulation cascade reactions that lead to fibrin formation. In the current work, we incorporate this missing element into a well-established and validated model for platelet activity. We then use this novel approach to predict thrombus formation in two experimental configurations. Our results confirm that to accurately predict the clotting process in devices where surface area to volume ratios are large and flow velocity and shear stresses remain low, coagulation reactions and subsequent fibrin formation must be considered. This new model could have great implications for the design and optimization of medical devices such as blood oxygenators. In the long term, the model could evolve into a functional tool to inform anticoagulation therapies for these patients.