Past and future of blood damage modelling in a view of translational research

Empirical and Computational Evaluation of Hemolysis in a Microfluidic Extracorporeal Membrane Oxygenator Prototype

Microfluidic devices promise to overcome the limitations of conventional hemodialysis and oxygenation technologies by incorporating novel membranes with ultra-high permeability into portable devices with low blood volume. However, the characteristically small dimensions of these devices contribute to both non-physiologic shear that could damage blood components and laminar flow that inhibits transport. While many studies have been performed to empirically and computationally study hemolysis in medical devices, such as valves and blood pumps, little is known about blood damage in microfluidic devices. In this study, four variants of a representative microfluidic membrane-based oxygenator and two controls (positive and negative) are introduced, and computational models are used to predict hemolysis. The simulations were performed in ANSYS Fluent for nine shear stress-based parameter sets for the power law hemolysis model. We found that three of the nine tested parameters overpredict (5 to 10×) hemolysis compared to empirical experiments. However, three parameter sets demonstrated higher predictive accuracy for hemolysis values in devices characterized by low shear conditions, while another three parameter sets exhibited better performance for devices operating under higher shear conditions. Empirical testing of the devices in a recirculating loop revealed levels of hemolysis significantly lower (<2 ppm) than the hemolysis ranges observed in conventional oxygenators (>10 ppm). Evaluating the model's ability to predict hemolysis across diverse shearing conditions, both through empirical experiments and computational validation, will provide valuable insights for future micro ECMO device development by directly relating geometric and shear stress with hemolysis levels. We propose that, with an informed selection of hemolysis parameters based on the shear ranges of the test device, computational modeling can complement empirical testing in the development of novel high-flow blood-contacting microfluidic devices, allowing for a more efficient iterative design process. Furthermore, the low device-induced hemolysis measured in our study at physiologically relevant flow rates is promising for the future development of microfluidic oxygenators and dialyzers.

A Fibrin-Thrombin Based In Vitro Perfusion System to Study Flow-Related Prosthetic Heart Valves Thrombosis

Prosthetic heart valve (PHV) replacement has increased the survival rate and quality of life for heart valve-diseased patients. However, PHV thrombosis remains a critical problem associated with these procedures. To better understand the PHV flow-related thrombosis problem, appropriate experimental models need to be developed. In this study, we present an in vitro fibrin clot model that mimics clot accumulation in PHVs under relevant hydrodynamic conditions while allowing real-time imaging. We created 3D-printed mechanical aortic valve models that were inserted into a transparent glass aorta model and connected to a system that simulates human aortic flow pulse and pressures. Thrombin was gradually injected into a circulating fibrinogen solution to induce fibrin clot formation, and clot accumulation was quantified via image analysis. The results of valves positioned in a normal versus a tilted configuration showed that clot accumulation correlated with the local flow features and was mainly present in areas of low shear and high residence time, where recirculating flows are dominant, as supported by computational fluid dynamic simulations. Overall, our work suggests that the developed method may provide data on flow-related clot accumulation in PHVs and may contribute to exploring new approaches and valve designs to reduce valve thrombosis.

Turbulence in surgical suction heads as detected by MRI

BACKGROUND

Blood loss is common during surgical procedures, especially in open cardiac surgery. Allogenic blood transfusion is associated with increased morbidity and mortality. Blood conservation programs in cardiac surgery recommend re-transfusion of shed blood directly or after processing, as this decreases transfusion rates of allogenic blood. But aspiration of blood from the wound area is often associated with increased hemolysis, due to flow induced forces, mainly through development of turbulence.

METHODS

We evaluated magnetic resonance imaging (MRI) as a qualitative tool for detection of turbulence. MRI is sensitive to flow; this study uses velocity-compensated T1-weighted 3D MRI for turbulence detection in four geometrically different cardiotomy suction heads under comparable flow conditions (0-1250 mL/min).

RESULTS

Our standard control suction head Model A showed pronounced signs of turbulence at all flow rates measured, while turbulence was only detectable in our modified Models 1-3 at higher flow rates (Models 1 and 3) or not at all (Model 2).

CONCLUSIONS

The comparison of flow performance of surgical suction heads with different geometries via acceleration-sensitized 3D MRI revealed significant differences in turbulence development between our standard control Model A and the modified alternatives (Models 1-3). As flow conditions during measurement have been comparable, the specific geometry of the respective suction heads must have been the main factor responsible. The underlying mechanisms and causative factors can only be speculated about, but as other investigations have shown, hemolytic activity is positively associated with degree of turbulence. The turbulence data measured in this study correlate with data from other investigations about hemolysis induced by surgical suction heads. The experimental MRI technique used showed added value for further elucidating the underlying physical phenomena causing blood damage due to non-physiological flow.

Flow simulation-based particle swarm optimization for developing improved hemolysis models

The improvement and development of blood-contacting devices, such as mechanical circulatory support systems, is a life saving endeavor. These devices must be designed in such a way that they ensure the highest hemocompatibility. Therefore, in-silico trials (flow simulations) offer a quick and cost-effective way to analyze and optimize the hemocompatibility and performance of medical devices. In that regard, the prediction of blood trauma, such as hemolysis, is the key element to ensure the hemocompatibility of a device. But, despite decades of research related to numerical hemolysis models, their accuracy and reliability leaves much to be desired. This study proposes a novel optimization path, which is capable of improving existing models and aid in the development of future hemolysis models. First, flow simulations of three, turbulent blood flow test cases (capillary tube, FDA nozzle, FDA pump) were performed and hemolysis was numerically predicted by the widely-applied stress-based hemolysis models. Afterward, a multiple-objective particles swarm optimization (MOPSO) was performed to tie the physiological stresses of the simulated flow field to the measured hemolysis using an equivalent of over one million numerically determined hemolysis predictions. The results show that our optimization is capable of improving upon existing hemolysis models. However, it also unveils some deficiencies and limits of hemolysis prediction with stress-based models, which will need to be addressed in order to improve its reliability.

The Effect of Blood Viscosity on Shear-Induced Hemolysis using a Magnetically Levitated Shearing Device

INTRODUCTION

Blood contacting medical devices, including rotary blood pumps, can cause shear-induced blood damage that may lead to adverse effects in patients. Due in part to an inadequate understanding of how cell-scale fluid mechanics impact red blood cell membrane deformation and damage, there is currently not a uniformly accepted engineering model for predicting blood damage caused by complex flow fields within ventricular assist devices (VADs).

METHODS

We empirically investigated hemolysis in a magnetically levitated axial Couette flow device typical of a rotary VAD. The device is able to accurately control the shear rate and exposure time experienced by blood and to minimize the effects of other uncharacterized stresses. Using this device, we explored the effects of both hematocrit and plasma viscosity on shear-induced hemolysis to characterize blood damage based on the viscosity-independent shear rate, rather than on shear stress.

RESULTS

Over a shear rate range of 20,000-80,000 1/s, the Index of Hemolysis (IH) was found to be dependent upon and well-predicted by shear rate alone. IH was independent of hematocrit, bulk viscosity, or the suspension media viscosity, and less correlated to shear stress (MSE=0.46-0.75) than to shear rate (MSE=0.06-0.09).

CONCLUSION

This study recommends that future investigations of shear-induced blood damage report findings with respect to the viscosity-neutral term of shear rate, in addition to the bulk whole blood viscosity measured at an appropriate shear rate relevant to the flow conditions of the device.

Equivalent Scalar Stress Formulation Taking into Account Non-Resolved Turbulent Scales

PURPOSE

Cardiovascular engineering includes flows with fluid-dynamical stresses as a parameter of interest. Mechanical stresses are high-risk factors for blood damage and can be assessed by computational fluid dynamics. By now, it is not described how to calculate an adequate scalar stress out of turbulent flow regimes when the whole share of turbulence is not resolved by the simulation method and how this impacts the stress calculation.

METHODS

We conducted direct numerical simulations (DNS) of test cases (a turbulent channel flow and the FDA nozzle) in order to access all scales of flow movement. After validation of both DNS with literature und experimental data using magnetic resonance imaging, the mechanical stress is calculated as a baseline. Afterwards, same flows are calculated using state-of-the-art turbulence models. The stresses are computed for every result using our definition of an equivalent scalar stress, which includes the influence from respective turbulence model, by using the parameter dissipation. Afterwards, the results are compared with the baseline data.

RESULTS

The results show a good agreement regarding the computed stress. Even when no turbulence is resolved by the simulation method, the results agree well with DNS data. When the influence of non-resolved motion is neglected in the stress calculation, it is underpredicted in all cases.

CONCLUSION

With the used scalar stress formulation, it is possible to include information about the turbulence of the flow into the mechanical stress calculation even when the used simulation method does not resolve any turbulence.

Turbulence and turbulent flow structures in a ventricular assist device-A numerical study using the large-eddy simulation

Numerical flow simulations that analyze the turbulent flow characteristics within a turbopump are important for optimizing the efficiency of such machines. In the case of ventricular assist devices (VADs), turbulent flow characteristics must be also examined in order to improve hemocompatibility. Turbulence increases the shear stresses in the VAD flow, which can lead to an increased damage to the transported blood components. Therefore, an understanding of the turbulent flow patterns and their significance for the numerical blood damage prediction is particularly important for flow optimizations in VADs in order to identify and thus minimize flow regions where blood could be damaged due to high turbulent stresses. Nevertheless, the turbulence occurring in VADs and the local turbulent structures that lead to increased turbulent stresses have not yet been analyzed in detail in these machines. Therefore, this study aims to investigate the turbulence in an axial VAD in a comprehensive and double tracked way. First, the flow in an axial VAD was computed using the large-eddy simulation method, and it was verified that the majority of the turbulence was directly resolved by the simulation. Then, the turbulent flow state of the VAD was quantified globally. For this purpose, a self-designed evaluation method, the power loss analysis, was used. Subsequently, local flow regions and flow structures were identified where significant turbulent stresses prevail. It will be shown that the identified regions are universal and will also occur in other axial blood pumps as well, for example, in the HeartMate II.

Influence of turbulent shear stresses on the numerical blood damage prediction in a ventricular assist device

The blood damage prediction in rotary blood pumps is an important procedure to evaluate the hemocompatibility of such systems. Blood damage is caused by shear stresses to the blood cells and their exposure times. The total impact of an equivalent shear stress can only be taken into account when turbulent stresses are included in the blood damage prediction. The aim of this article was to analyze the influence of the turbulent stresses on the damage prediction in a rotary blood pump's flow. Therefore, the flow in a research blood pump was computed using large eddy simulations. A highly turbulence-resolving setup was used in order to directly resolve most of the computed stresses. The simulations were performed at the design point and an operation point with lower flow rate. Blood damage was predicted using three damage models (volumetric analysis of exceeded stress thresholds, hemolysis transport equation, and hemolysis approximation via volume integral) and two shear stress definitions (with and without turbulent stresses). For both simulations, turbulent stresses are the dominant stresses away from the walls. Here, they act in a range between 9 and 50 Pa. Nonetheless, the mean stresses in the proximity of the walls reach levels, which are one order of magnitude higher. Due to this, the turbulent stresses have a small impact on the results of the hemolysis prediction. Yet, turbulent stresses should be included in the damage prediction, since they belong to the total equivalent stress definition and could impact the damage on proteins or platelets.