Purpose: Pediatric and adult patients with sickle cell anemia (SCA) are at increased risk of stroke and cerebrovascular accident. In the general adult population, there is a relationship between arterial hemodynamics and pathology; however, this relationship in SCA patients remains to be elucidated. The aim of this work was to characterize circle of Willis hemodynamics in patients with SCA and quantify the impact of viscosity choice on pathophysiologically-relevant hemodynamics measures. Methods: Based on measured vascular geometries, time-varying flow rates, and blood parameters, detailed patient-specific simulations of the circle of Willis were conducted for SCA patients (n = 6). Simulations quantified the impact of patient-specific and standard blood viscosities on wall shear stress (WSS). Results: These results demonstrated that use of a standard blood viscosity introduces large errors into the estimation of pathophysiologically-relevant hemodynamic parameters. Standard viscosity models overpredicted peak WSS by 55% and 49% for steady and pulsatile flow, respectively. Moreover, these results demonstrated non-uniform, spatial patterns of positive and negative WSS errors related to viscosity, and standard viscosity simulations overpredicted the time-averaged WSS by 32% (standard deviation = 7.1%). Finally, differences in shear rate demonstrated that the viscosity choice alters the simulated near-wall flow field, impacting hemodynamics measures. Conclusions: This work presents simulations of circle of Willis arterial flow in SCA patients and demonstrates the importance and feasibility of using a patient-specific viscosity in these simulations. Accurately characterizing cerebrovascular hemodynamics in SCA populations has potential for elucidating the pathophysiology of large-vessel occlusion, aneurysms, and tissue damage in these patients.
Total cavopulmonary connection is the result of a series of palliative surgical repairs performed on patients with single ventricle heart defects. The resulting anatomy has complex and unsteady hemodynamics characterized by flow mixing and flow separation. Although varying degrees of flow pulsatility have been observed in vivo, non-pulsatile (time-averaged) boundary conditions have traditionally been assumed in hemodynamic modeling, and only recently have pulsatile conditions been incorporated without completely characterizing their effect or importance. In this study, 3D numerical simulations with both pulsatile and non-pulsatile boundary conditions were performed for 24 patients with different anatomies and flow boundary conditions from Georgia Tech database. Flow structures, energy dissipation rates and pressure drops were compared under rest and simulated exercise conditions. It was found that flow pulsatility is the primary factor in determining the appropriate choice of boundary conditions, whereas the anatomic configuration and cardiac output had secondary effects. Results show that the hemodynamics can be strongly influenced by the presence of pulsatile flow. However, there was a minimum pulsatility threshold, identified by defining a weighted pulsatility index (wPI), above which the influence was significant. It was shown that when wPI < 30%, the relative error in hemodynamic predictions using time-averaged boundary conditions was less than 10% compared to pulsatile simulations. In addition, when wPI < 50, the relative error was less than 20%. A correlation was introduced to relate wPI to the relative error in predicting the flow metrics with non-pulsatile flow conditions.