Incorporating Unresolved Stresses in Blood Damage Modeling: Energy Dissipation More Accurate Than Reynolds Stress Formulation

OBJECTIVE

Reynolds Averaged Navier Stokes (RANS) models are often used as the basis for modeling blood damage in turbulent flows. To predict blood damage by turbulence stresses that are not resolved in RANS, a stress formulation that represents the corresponding scales is required. Here, we compare two commonly employed stress formulations: a scalar stress representation that uses Reynolds stresses as a surrogate for unresolved fluid stresses, and an effective stress formulation based on energy dissipation.

METHODS

We conducted unsteady RANS simulations of the CentriMag blood pump with three different closure models and a Large Eddy Simulation (LES) for reference. We implemented both stress representations in all models and compared the resulting total stress distributions in Eulerian and Lagrangian frameworks.

RESULTS

The Reynolds-stress-based approach overestimated the contribution of unresolved stresses in RANS, with differences between closure models of up to several orders of magnitude. With the dissipation-based approach, the total stresses predicted with RANS deviated by about 50% from the LES reference, which was more accurate than only considering resolved stresses.

CONCLUSION

The Reynolds-stress-based formulation proved unreliable for estimating scalar stresses in our RANS simulations, while the dissipation-based approach provided an accuracy improvement over simply neglecting unresolved stresses.

SIGNIFICANCE

Our results suggest that dissipation-based inclusion of unresolved stresses should be the preferred choice for blood damage modeling in RANS.

Computational Fluid Dynamics Turbulence Model and Experimental Study for a Fontan Cavopulmonary Assist Device

Head-flow HQ curves for a Fontan cavopulmonary assist device (CPAD) were measured using a blood surrogate in a mock circulatory loop and simulated with various computational fluid dynamics (CFD) models. The tests benchmarked the CFD tools for further enhancement of the CPAD design. Recommended Reynolds-Averaged Navier-Stokes (RANS) CFD approaches for the development of conventional ventricular assist devices (VAD) were found to have shortcomings when applied to the Fontan CPAD, which is designed to neutralize off-condition obstruction risks that could contribute to a major adverse event. The no-obstruction condition is achieved with a von Karman pump, utilizing large clearances and small blade heights, which challenge conventional VAD RANS-based CFD hemodynamic simulations. High-fidelity large eddy simulation (LES) is always recommended; however, this may be cost-inhibitive for optimization studies in commercial settings, thus the reliance on RANS models. This study compares head and power predictions of various RANS turbulence models, employing experimental measurements and LES results as a basis for comparison. The models include standard k-ϵ, re-normalization group k-ϵ, realizable k-ϵ, shear stress transport (SST) k-ω, SST with transitional turbulence, and Generalized k-ω. For the pressure head predictions, it was observed that the standard k-ϵ model provided far better agreement with experiment. For the rotor torque, k-ϵ predictions were 30% lower than LES, while the SST and LES torque values were near identical. For the Fontan CPAD, the findings support using LES for the final design simulations, k-ϵ model for head and general flow simulation, and SST for power, shear stress, hemolysis, and thrombogenicity predictions.