Quantification of Discretization Error in Wall Shear Stress Calculations for Aneurismal Flows Using CFD

Many recent studies suggest that hemodynamic factors such as wall shear stress (WSS) and pressure contribute to the genesis and growth of intracranial aneurysms. Recently there have been a number of computational hemodynamics studies that calculate the values of wall shear stress in arterial and aneurismal flows. However there is a lack of comprehensive error analysis in many of the computational hemodynamics studies. This is perhaps the reason for speculative and ambiguous conclusions drawn by various studies as to the nature of wall shear stress responsible for aneurysm growth. In the current study, geometry involving an actual aneurysm is built from angiogram images. Another geometry consisting of the primary artery where the aneurysm formed is also built by removing the aneurysm volume. The two geometries are meshed using three different grid densities. Second order schemes are used to simulate the pulsatile hemodynamics through each of the geometries. Various representative planes along the geometries are considered and the major flow variables and WSS are plotted as a function of grid densities. The procedure for estimation of discretization error, suggested by ASME Journal of Fluids Engineering, is applied at various representative locations along the aneurysm and arterial geometry. The results suggest high dependence of calculated WSS on local grid density. The contours of WSS in the arterial geometry suggest that high WSS does not necessarily occur at the location where the aneurysm originated. Possible remedies are suggested so that this uncertainty could be eliminated from future studies.