Standardized viscosity as a source of error in computational fluid dynamic simulations of cerebral aneurysms

Abstract

BackgroundComputational fluid dynamics (CFD) simulations are a powerful tool for studying cerebral aneurysms, capable of evaluating hemodynamics in a way that is infeasible with imaging alone. However, the difficulty of incorporating patient‐specific information and inherent obstacles of in vivo validation have limited the clinical usefulness of CFD of cerebral aneurysms. In this work we investigate the effect of using standardized blood viscosity values in CFD simulations of cerebral aneurysms when compared to simulations of the same aneurysms using patient‐specific viscosity values derived from hematocrit measurements.PurposeThe objective of this work is to determine the level of error, on average, that is, caused by using standardized values of viscosity in CFD simulations of cerebral aneurysms. By quantifying this error, we demonstrate the need for incorporating patient‐specific viscosity in future CFD investigations of cerebral aneurysms.MethodsCFD simulations of forty‐one cerebral aneurysms were conducted using patient‐specific boundary conditions. For each aneurysm two simulations were conducted, one utilizing patient‐specific blood viscosity derived from hematocrit measurements and another using a standardized value for blood viscosity. Hemodynamic parameters such as wall shear stress (WSS), wall shear stress gradient (WSSG), and the oscillatory shear index (OSI) were calculated for each of the simulations for each aneurysm. Paired t‐tests for differences in the time‐averaged maps of these hemodynamic parameters between standardized and patient‐specific viscosity simulations were conducted for each aneurysm. Bland–Altman analysis was used to examine the cohort‐wide changes in the hemodynamic parameters. Subjects were broken into two groups, those with higher than standard viscosity and those with lower than standard viscosity. An unpaired t‐test was used to compare the percent change in WSS, WSSG, and OSI between patient‐specific and standardized viscosity simulations for the two cohorts. The percent changes in hemodynamic parameters were correlated against the direction and magnitude of percent change in viscosity, aneurysm size, and aneurysm location. For all t‐tests, a Bonferroni‐corrected significance level of 0.0167 was used.Results63.2%, 41.5%, and 48.7% of aneurysms showed statistically significant differences between patient‐specific and standardized viscosity simulations for WSS, WSSG, and OSI respectively. No statistically significant difference was found in the percent changes in WSS, WSSG, and OSI between the group with higher than standard viscosity and those with lower than standard viscosity, indicating an increase in viscosity can cause either an increase or decrease in each of the hemodynamic parameters. On a study‐wide level no significant bias was found in either direction for WSS, WSSG, or OSI between the simulation groups due to the bidirectional effect of changing viscosity. No correlation was found between percent change of viscosity and percent change of WSS, WSSG, or OSI, meaning an after‐the‐fact correction for patient‐specific viscosity is not feasible.ConclusionStandardizing viscosity values in CFD of cerebral aneurysms has a large and unpredictable impact on the calculated WSS, WSSG, and OSI when compared to CFD simulations of the same aneurysms using a patient‐specific viscosity. We recommend implementing hematocrit‐based patient‐specific blood viscosity values for all CFD simulations of cerebral aneurysms.

A practical strategy for data assimilation of cerebral intra-aneurysmal flows using a variational method with boundary control of velocity

“…More efficient gradient-based optimization techniques, such as the Broyden–Fletcher– Goldfarb–Shanno method [11], will mitigate this issue. Fifth, standardized blood viscosity has a large and unpredictable impact on hemodynamic factors in CFD analyses [29]. Uncertainty analyses or extension of the optimization problem to include blood viscosity will be needed to correctly quantify patient-specific hemodynamic features in DA analyses.…”