Advanced Constitutive Modeling of the Thixotropic Elasto-Visco-Plastic Behavior of Blood: Steady-State Blood Flow in Microtubes

The present work focuses on the in-silico investigation of the steady-state blood flow in straight microtubes, incorporating advanced constitutive modeling for human blood and blood plasma. The blood constitutive model accounts for the interplay between thixotropy and elasto-visco-plasticity via a scalar variable that describes the level of the local blood structure at any instance. The constitutive model is enhanced by the non-Newtonian modeling of the plasma phase, which features bulk viscoelasticity. Incorporating microcirculation phenomena such as the cell-free layer (CFL) formation or the Fahraeus and the Fahraeus-Lindqvist effects is an indispensable part of the blood flow investigation. The coupling between them and the momentum balance is achieved through correlations based on experimental observations. Notably, we propose a new simplified form for the dependence of the apparent viscosity on the hematocrit that predicts the CFL thickness correctly. Our investigation focuses on the impact of the microtube diameter and the pressure-gradient on velocity profiles, normal and shear viscoelastic stresses, and thixotropic properties. We demonstrate the microstructural configuration of blood in steady-state conditions, revealing that blood is highly aggregated in narrow tubes, promoting a flat velocity profile. Additionally, the proper accounting of the CFL thickness shows that for narrow microtubes, the reduction of discharged hematocrit is significant, which in some cases is up to 70%. At high pressure-gradients, the plasmatic proteins in both regions are extended in the flow direction, developing large axial normal stresses, which are more significant in the core region. We also provide normal stress predictions at both the blood/plasma interface (INS) and the tube wall (WNS), which are difficult to measure experimentally. Both decrease with the tube radius; however, they exhibit significant differences in magnitude and type of variation. INS varies linearly from 4.5 to 2 Pa, while WNS exhibits an exponential decrease taking values from 50 mPa to zero.

Automated summary

This study focuses on the in-silico investigation of steady-state blood flow in straight microtubes, incorporating advanced constitutive modeling for human blood and blood plasma. The research demonstrates the microstructural configuration of blood in steady-state conditions, revealing highly aggregated blood in narrow tubes, and the extension of plasmatic proteins in the flow direction at high pressure-gradients generating large axial normal stresses. Predictions of normal stress at both the blood/plasma interface and the tube wall show significant differences in magnitude and type of variation, with INS linearly decreasing from 4.5 to 2 Pa and WNS exponentially decreasing from 50 mPa to zero with decreasing tube radius.