Abstract
The non-Newtonian properties of blood flow have been widely debated in hemodynamic research, particularly for congenital heart defects. Many studies comparing Newtonian and non-Newtonian models have overlooked dimensional group consistency, resulting in comparisons influenced by inconsistent Reynolds numbers rather than viscosity effects. In this study, we address this issue by applying a generalized Reynolds number formulation to ensure consistent dimensionless group comparisons. We compare flow structures and hemodynamic metrics in 20 pediatric Fontan circulations using the non-Newtonian Casson model against both conventional and generalized Reynolds number-corrected Newtonian models. Our results show that the conventional Newtonian model significantly overestimates flow rotation and underestimates stagnation regions, potentially misrepresenting thrombosis risk. The generalized Reynolds number method, however, predicts flow structures, wall shear stress, and energy-based metrics more in line with the non-Newtonian model. Percentage of power loss estimates from the generalized method (17.7 [10.1, 22.7]; p < 0.05 ) align more closely with the non-Newtonian model (12.9 [7.0, 17.1]) than with the conventional approach (8.5 [4.3, 10.2]; p < 0.001 ), offering a more clinically relevant prediction. Additionally, indexed viscous dissipation from the generalized method (2.14 [1.17, 3.69] n.d.) is statistically indistinguishable (p=0.97) from the non-Newtonian model (2.42 [1.07, 3.60] n.d.; p < 0.05 ). Our analysis highlights that while the generalized Reynolds number method cannot fully replicate local shear-thinning effects, it substantially improves upon the conventional Newtonian approach by correcting for viscosity mismatch. We emphasize the importance of dimensionless group consistency before drawing conclusions in hemodynamic studies and advocate for broader adoption of non-Newtonian models to obtain critical clinical insights.