This work focuses on the advanced modeling of the thixotropic nature of blood, coupled with an elasto-visco-plastic formulation by invoking a consistent and validated model for TEVP materials. The proposed model has been verified for the adequate description of the rheological behavior of suspensions, introducing a scalar variable that describes dynamically the level of internal microstructure of rouleaux at any instance, capturing accurately the aggregation and disaggregation mechanisms of the RBCs. Also, a non-linear fitting is adopted for the definition of the model’s parameters on limited available experimental data of steady and transient rheometric flows of blood samples. We present the predictability of the new model in various steady and transient rheometric flows, including startup shear, rectangular shear steps, shear cessation, triangular shear steps and LAOS tests. Our model provides predictions for the elasto-thixotropic mechanism in startup shear flows, demonstrating a non-monotonic relationship of the thixotropic index on the shear-rate. The intermittent shear step test reveals the dynamics of the structural reconstruction, which in turn is associated with the aggregation process. Moreover, our model offers robust predictions for less examined tests such as uniaxial elongation, in which normal stress was found to have considerable contribution. Apart from the integrated modeling of blood rheological complexity, our implementation is adequate for multi-dimensional simulations due to its tensorial formalism accomplished with a single time scale for the thixotropic effects, resulting in a low computational cost compared to other TEVP models.

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 Fåhraeus and the Fåhraeus-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.