Robustness of common hemodynamic indicators with respect to numerical resolution in 38 middle cerebral artery aneurysms

Background Using computational fluid dynamics (CFD) to compute the hemodynamics in cerebral aneurysms has received much attention in the last decade. The usability of these methods depends on the quality of the computations, highlighted in recent discussions. The purpose of this study is to investigate the convergence of common hemodynamic indicators with respect to numerical resolution. Methods 38 middle cerebral artery bifurcation aneurysms were studied at two different resolutions (one comparable to most studies, and one finer). Relevant hemodynamic indicators were collected from two of the most cited studies, and were compared at the two refinements. In addition, correlation to rupture was investigated. Results Most of the hemodynamic indicators were very well resolved at the coarser resolutions, correlating with the finest resolution with a correlation coefficient >0.95. The oscillatory shear index (OSI) had the lowest correlation coefficient of 0.83. A logarithmic Bland-Altman plot revealed noticeable variations in the proportion of the aneurysm under low shear, as well as in spatial and temporal gradients not captured by the correlation alone. Conclusion Statistically, hemodynamic indicators agree well across the different resolutions studied here. However, there are clear outliers visible in several of the hemodynamic indicators, which suggests that special care should be taken when considering individual assessment.

[1]  A. Algra,et al.  Prevalence and risk of rupture of intracranial aneurysms: a systematic review. , 1998, Stroke.

[2]  Jianping Xiang,et al.  Rupture Resemblance Score (RRS): toward risk stratification of unruptured intracranial aneurysms using hemodynamic–morphological discriminants , 2014, Journal of NeuroInterventional Surgery.

[3]  J. Womersley Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known , 1955, The Journal of physiology.

[4]  Makoto Yamamoto,et al.  Hemodynamic Differences Between Unruptured and Ruptured Intracranial Aneurysms During Observation , 2012, Stroke.

[5]  Anders Logg,et al.  Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book , 2012 .

[6]  H. Marquering,et al.  Rupture-Associated Changes of Cerebral Aneurysm Geometry: High-Resolution 3D Imaging before and after Rupture , 2014, American Journal of Neuroradiology.

[7]  C. Putman,et al.  Quantitative Characterization of the Hemodynamic Environment in Ruptured and Unruptured Brain Aneurysms , 2010, American Journal of Neuroradiology.

[8]  David A Steinman,et al.  Characterization of volumetric flow rate waveforms at the carotid bifurcations of older adults , 2010, Physiological measurement.

[9]  M. Markl,et al.  Dampening of Blood-Flow Pulsatility along the Carotid Siphon: Does Form Follow Function? , 2011, American Journal of Neuroradiology.

[10]  J. C. Simo,et al.  Unconditional stability and long-term behavior of transient algorithms for the incompressible Navier-Stokes and Euler equations , 1994 .

[11]  Alejandro F Frangi,et al.  @neurIST complex information processing toolchain for the integrated management of cerebral aneurysms , 2011, Interface Focus.

[12]  M Umezu,et al.  Risk Analysis of Unruptured Aneurysms Using Computational Fluid Dynamics Technology: Preliminary Results , 2011, American Journal of Neuroradiology.

[13]  D F Kallmes,et al.  Point: CFD—Computational Fluid Dynamics or Confounding Factor Dissemination , 2012, American Journal of Neuroradiology.

[14]  Kendall D Dennis,et al.  Grid convergence errors in hemodynamic solution of patient-specific cerebral aneurysms. , 2012, Journal of biomechanics.

[15]  Alejandro F Frangi,et al.  Morphological descriptors as rupture indicators in middle cerebral artery aneurysms , 2010, 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology.

[16]  T F Sherman,et al.  On connecting large vessels to small. The meaning of Murray's law , 1981, The Journal of general physiology.

[17]  F. Mut,et al.  Association of Hemodynamic Characteristics and Cerebral Aneurysm Rupture , 2011, American Journal of Neuroradiology.

[18]  D. Steinman,et al.  Estimation of Inlet Flow Rates for Image-Based Aneurysm CFD Models: Where and How to Begin? , 2015, Annals of Biomedical Engineering.

[19]  D. Altman,et al.  Measuring agreement in method comparison studies , 1999, Statistical methods in medical research.

[20]  David A Steinman,et al.  High-resolution CFD detects high-frequency velocity fluctuations in bifurcation, but not sidewall, aneurysms. , 2013, Journal of biomechanics.

[21]  Eric L. Miller,et al.  3D Shape Analysis of Intracranial Aneurysms Using the Writhe Number as a Discriminant for Rupture , 2011, Annals of Biomedical Engineering.

[22]  Alejandro F. Frangi,et al.  Automated segmentation of cerebral vasculature with aneurysms in 3DRA and TOF-MRA using geodesic active regions: an evaluation study. , 2010, Medical physics.

[23]  J. Mocco,et al.  Hemodynamic–Morphologic Discriminants for Intracranial Aneurysm Rupture , 2011, Stroke.

[24]  Kristian Valen-Sendstad,et al.  A study of wall shear stress in 12 aneurysms with respect to different viscosity models and flow conditions. , 2013, Journal of biomechanics.

[25]  L. Jou,et al.  Wall Shear Stress on Ruptured and Unruptured Intracranial Aneurysms at the Internal Carotid Artery , 2008, American Journal of Neuroradiology.

[26]  J. Mocco,et al.  MORPHOLOGY PARAMETERS FOR INTRACRANIAL ANEURYSM RUPTURE RISK ASSESSMENT , 2008, Neurosurgery.

[27]  Fujimaro Ishida,et al.  Low Wall Shear Stress Is Independently Associated With the Rupture Status of Middle Cerebral Artery Aneurysms , 2013, Stroke.

[28]  Sabine Roller,et al.  Transitional flow in intracranial aneurysms - a space and time refinement study below the Kolmogorov scales using Lattice Boltzmann Method , 2015, ArXiv.

[29]  D. Steinman,et al.  Narrowing the Expertise Gap for Predicting Intracranial Aneurysm Hemodynamics: Impact of Solver Numerics versus Mesh and Time-Step Resolution , 2015, American Journal of Neuroradiology.

[30]  Mikael Mortensen,et al.  Oasis: A high-level/high-performance open source Navier-Stokes solver , 2015, Comput. Phys. Commun..

[31]  Alejandro F. Frangi,et al.  Morphological Characterization of Intracranial Aneurysms Using 3-D Moment Invariants , 2007, IEEE Transactions on Medical Imaging.

[32]  H. Langtangen,et al.  Direct numerical simulation of transitional flow in a patient-specific intracranial aneurysm. , 2011, Journal of biomechanics.

[33]  Rohini Retarekar Hemodynamics and natural history outcome in unruptured intracranial aneurysms , 2012 .

[34]  Prahlad G. Menon,et al.  Variability of computational fluid dynamics solutions for pressure and flow in a giant aneurysm: the ASME 2012 Summer Bioengineering Conference CFD Challenge. , 2013, Journal of biomechanical engineering.

[35]  David A Steinman,et al.  High-resolution computational fluid dynamics detects flow instabilities in the carotid siphon: implications for aneurysm initiation and rupture? , 2014, Journal of biomechanics.

[36]  D. Steinman,et al.  Mind the Gap: Impact of Computational Fluid Dynamics Solution Strategy on Prediction of Intracranial Aneurysm Hemodynamics and Rupture Status Indicators , 2014, American Journal of Neuroradiology.

[37]  M. L. Raghavan,et al.  Quantified aneurysm shape and rupture risk. , 2005, Journal of neurosurgery.

[38]  J Max Findlay,et al.  The aspect ratio (dome/neck) of ruptured and unruptured aneurysms. , 2003, Journal of neurosurgery.

[39]  A. Algra,et al.  Case-fatality rates and functional outcome after subarachnoid hemorrhage: a systematic review. , 1997, Stroke.

[40]  F. Kajiya,et al.  Effects of size and shape (aspect ratio) on the hemodynamics of saccular aneurysms: a possible index for surgical treatment of intracranial aneurysms. , 1999, Neurosurgery.

[41]  Y Ventikos Resolving the Issue of Resolution , 2014, American Journal of Neuroradiology.