Effects of hematocrit and non-Newtonian blood rheology on pulsatile wall shear stress distributions in vascular anomalies: A multiple relaxation time lattice Boltzmann approach

In the present study, we investigate the flow dynamics of non-Newtonian blood, focusing on the distribution of wall shear stress (WSS) and hematocrit levels, which is the volume percentage of red blood cells in whole blood. We analyze these factors under pulsatile conditions, in vascular anomalies such as stent channels and intracranial aneurysms. To achieve this, a three-dimensional computational approach based on the lattice Boltzmann method (LBM) with a multiple relaxation time (MRT) collision operator is employed. To represent the blood's shear-thinning properties, we developed a constitutive model inspired by the Carreau–Yasuda model. This model considers the variability in blood viscosity with shear rate correlated with hematocrit levels based on experimental data documented in the literature. The accuracy of the employed MRT-LBM is demonstrated by the consistency of results with analytical solutions for steady state and experimental data for pulsatile WSS distributions in non-Newtonian and Newtonian fluids. Results indicate that, in areas narrowed by stenosis or expanded by aneurysms, hematocrit levels affect flow dynamics. Higher hematocrit levels intensify pulsatile flow through stenotic regions, increasing WSS cyclic variations. We derived a density distribution function to demonstrate how shear rates vary in vascular anomalies, revealing blood viscosity changes and non-Newtonian properties. These properties complicate flow patterns, resulting in non-linear WSS distributions, which are essential for understanding endothelial cell reactions and disease pathways. Pulsatile blood flow and altered rheological properties due to increased hematocrit affect saccular aneurysm fluid dynamics over time and space, causing vorticities to change shape, size, and intensity.