Non-Newtonian Patient-specific Simulations of Left Atrial Hemodynamics
Atrial fibrillation (AF) is the most common arrhythmia, affecting ~35M people worldwide. The irregular beating of the left atrial (LA) caused by AF impacts the LA hemodynamics increasing the risk of thrombosis and ischemic stroke. Most LA thrombi appear in its appendage (LAA), a narrow sac of varied morphology where blood is prone to stagnate. In the LAA, the combination of slow blood flow and low shear rates (<100 [1/s]) promotes the formation of red blood cell aggregations called rouleaux. Blood experiences a non-Newtonian behavior when rouleaux formed that has not been considered in previous CFD analysis of the LA. We model the anatomy and motion of the LA from 4D-CT images and solve the blood flow inside the LA geometry with our CFD in-house code, which models Non-Newtonian rheology with the shear-hematocrit-dependent Carreau-Yasuda equation. We cover a wide range of non-Newtonian effects considering a small and a large hematocrit, including an additional constitutive relation to account for themrouleaux formation time, and we compare our results with Newtonian simulations. Blood rheology influence in LAA hemostasis is studied in 6 patient-specific anatomies. Two subjects had an LAA thrombus (digitally removed before running the simulations), another had a history of TIAs, and the remaining three had normal atrial function. In our simulations, the shear rate remains below 50 [1/s] in the LAA for all non-Newtonian models considered. This triggers an increase of viscosity that alters the flow behavior in that site, which exhibits different flow patterns than Newtonian simulations. These hemodynamic changes translate into differences in the LAA hemostasis, calculated with the residence time.