An Accelerated Thrombosis Model for Computational Fluid Dynamics Simulations in Rotary Blood Pumps

Numerical Study on the Impact of Central Venous Catheter Placement on Blood Flow in the Cavo-Atrial Junction

An in silico study is performed to investigate fluid dynamic effects of central venous catheter (CVC) placement within patient-specific cavo-atrial junctions. Prior studies show the CVC infusing a liquid, but this study focuses on the placement without any liquid emerging from the CVC. A 7 or 15-French double-lumen CVC is placed virtually in two patient-specific models; the CVC tip location is altered to understand its effect on the venous flow field. Results show that the CVC impact is trivial on flow in the superior vena cava when the catheter-to-vein ratio ranges from 0.15 to 0.33. Results further demonstrate that when the CVC tip is directly in the right atrium, flow vortices in the right atrium result in elevated wall shear stress near the tip hole. A recirculation region characterizes a spatially variable flow field inside the CVC side hole. Furthermore, flow stagnation is present near the internal side hole corners but an elevated wall shear stress near the curvature of the side hole's exit. These results suggest that optimal CVC tip location is within the superior vena cava, so as to lower the potential for platelet activation due to elevated shear stresses and that CVC geometry and location depth in the central vein significantly influences the local CVC fluid dynamics. A thrombosis model also shows thrombus formation at the side hole and tip hole. After modifying the catheter design, the hemodynamics change, which alter thrombus formation. Future studies are warranted to study CVC design and placement location in an effort to minimize CVC-induced thrombosis incidence.

[Progress in the analysis of hemolysis and coagulation models for interventional micro-axial flow blood pumps]

Interventional micro-axial flow blood pump is widely used as an effective treatment for patients with cardiogenic shock. Hemolysis and coagulation are vital concerns in the clinical application of interventional micro-axial flow pumps. This paper reviewed hemolysis and coagulation models for micro-axial flow blood pumps. Firstly, the structural characteristics of commercial interventional micro-axial flow blood pumps and issues related to clinical applications were introduced. Then the basic mechanisms of hemolysis and coagulation were used to study the factors affecting erythrocyte damage and platelet activation in interventional micro-axial flow blood pumps, focusing on the current models of hemolysis and coagulation on different scales (macroscopic, mesoscopic, and microscopic). Since models at different scales have different perspectives on the study of hemolysis and coagulation, a comprehensive analysis combined with multi-scale models is required to fully consider the influence of complex factors of interventional pumps on hemolysis and coagulation.

A multi-constituent model for assessing the effect of impeller shroud on the thrombosis potential of a centrifugal blood pump

Centrifugal blood pumps can be used for treating heart failure patients. However, pump thrombosis has remained one of the complications that trouble clinical treatment. This study analyzed the effect of impeller shroud on the thrombosis risk of the blood pump, and predicted areas prone to thrombosis. Multi-constituent transport equations were presented, considering mechanical activation and biochemical activation. It was found that activated platelets concentration can increase with shear stress and adenosine diphosphate(ADP) concentration increasing, and the highest risk of thrombosis inside the blood pump was under extracorporeal membrane oxygenation (ECMO) mode. Under the same condition, ADP concentration and thrombosis index of semi-shroud impeller can increase by 7.3% and 7.2% compared to the closed-shroud impeller. The main reason for the increase in thrombosis risk was owing to elevated scalar shear stress and more coagulation promoting factor-ADP released. The regions with higher thrombosis potential were in the center hole, top and bottom clearance. As a novelty, the findings revealed that impeller shroud can influence mechanical and biochemical activation factors. It is useful for identifying potential risk regions of thrombus formation based on relative comparisons.

Study on the optimal elastic modulus of flexible blades for right heart assist device supporting patients with single-ventricle physiologies

Background

Patients with single-ventricle physiologies continue to experience insufficient circulatory power after undergoing palliative surgeries. This paper proposed a right heart assist device equipped with flexible blades to provide circulatory assistance for these patients. The optimal elastic modulus of the flexible blades was investigated through numerical simulation.

Methods

A one-way fluid-structure interaction (FSI) simulation was employed to study the deformation of flexible blades during rotation and its impact on device performance. The process began with a computational fluid dynamics (CFD) simulation to calculate the blood pressure rise and the pressure on the blades' surface. Subsequently, these pressure data were exported for finite element analysis (FEA) to compute the deformation of the blades. The fluid domain was then recreated based on the deformed blades' shape. Iterative CFD and FEA simulations were performed until both the blood pressure rise and the blades' shape stabilized. The blood pressure rise, hemolysis risk, and thrombosis risk corresponding to blades with different elastic moduli were exhaustively evaluated to determine the optimal elastic modulus.

Results

Except for the case at 8,000 rpm with a blade elastic modulus of 40 MPa, the pressure rise associated with flexible blades within the studied range (rotational speeds of 4,000 rpm and 8,000 rpm, elastic modulus between 10 MPa and 200 MPa) was lower than that of rigid blades. It was observed that the pressure rise corresponding to flexible blades increased as the elastic modulus increased. Additionally, no significant difference was found in the hemolysis risk and thrombus risk between flexible blades of various elastic moduli and rigid blades.

Conclusion

Except for one specific case, deformation of the flexible blades within the studied range led to a decrease in the impeller's functionality. Notably, rotational speed had a more significant impact on hemolysis risk and thrombus risk compared to blade deformation. After a comprehensive analysis of blade compressibility, blood pressure rise, hemolysis risk, and thrombus risk, the optimal elastic modulus for the flexible blades was determined to be between 40 MPa and 50 MPa.

A fibrin enhanced thrombosis model for medical devices operating at low shear regimes or large surface areas

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⁻¹), 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.

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.

A New Mathematical Numerical Model to Evaluate the Risk of Thrombosis in Three Clinical Ventricular Assist Devices

(1) Background: Thrombosis is the main complication in patients supported with ventricular assist devices (VAD). Models that accurately predict the risk of thrombus formation in VADs are still lacking. When VADs are clinically assisted, their complex geometric configuration and high rotating speed inevitably generate complex flow fields and high shear stress. These non-physiological factors can damage blood cells and proteins, release coagulant factors and trigger thrombosis. In this study, a more accurate model for thrombus assessment was constructed by integrating parameters such as shear stress, residence time and coagulant factors, so as to accurately assess the probability of thrombosis in three clinical VADs. (2) Methods: A mathematical model was constructed to assess platelet activation and thrombosis within VADs. By solving the transport equation, the influence of various factors such as shear stress, residence time and coagulation factors on platelet activation was considered. The diffusion equation was applied to determine the role of activated platelets and substance deposition on thrombus formation. The momentum equation was introduced to describe the obstruction to blood flow when thrombus is formed, and finally a more comprehensive and accurate model for thrombus assessment in patients with VAD was obtained. Numerical simulations of three clinically VADs (CH-VAD, HVAD and HMII) were performed using this model. The simulation results were compared with experimental data on platelet activation caused by the three VADs. The simulated thrombogenic potential in different regions of MHII was compared with the frequency of thrombosis occurring in the regions in clinic. The regions of high thrombotic risk for HVAD and HMII observed in experiments were compared with the regions predicted by simulation. (3) Results: It was found that the percentage of activated platelets within the VAD obtained by solving the thrombosis model developed in this study was in high agreement with the experimental data (r² = 0.984), the likelihood of thrombosis in the regions of the simulation showed excellent correlation with the clinical statistics (r² = 0.994), and the regions of high thrombotic risk predicted by the simulation were consistent with the experimental results. Further study revealed that the three clinical VADs (CH-VAD, HVAD and HMII) were prone to thrombus formation in the inner side of the secondary flow passage, the clearance between cone and impeller, and the corner region of the inlet pipe, respectively. The risk of platelet activation and thrombus formation for the three VADs was low to high for CH-VAD, HVAD, and HM II, respectively. (4) Conclusions: In this study, a more comprehensive and accurate thrombosis model was constructed by combining parameters such as shear stress, residence time, and coagulation factors. Simulation results of thrombotic risk received with this model showed excellent correlation with experimental and clinical data. It is important for determining the degree of platelet activation in VAD and identifying regions prone to thrombus formation, as well as guiding the optimal design of VAD and clinical treatment.