Newtonian viscosity model could overestimate wall shear stress in intracranial aneurysm domes and underestimate rupture risk

Objective Computational fluid dynamics (CFD) simulations of intracranial aneurysm hemodynamics usually adopt the simplification of the Newtonian blood rheology model. A study was undertaken to examine whether such a model affects the predicted hemodynamics in realistic intracranial aneurysm geometries. Methods Pulsatile CFD simulations were carried out using the Newtonian viscosity model and two non-Newtonian models (Casson and Herschel-Bulkley) in three typical internal carotid artery saccular aneurysms (A, sidewall, oblong-shaped with a daughter sac; B, sidewall, quasi-spherical; C, near-spherical bifurcation). For each aneurysm model the surface distributions of shear rate, blood viscosity and wall shear stress (WSS) predicted by the three rheology models were compared. Results All three rheology models produced similar intra-aneurysmal flow patterns: aneurysm A had a slowly recirculating secondary vortex near the dome whereas aneurysms B and C contained only a large single vortex. All models predicted similar shear rate, blood viscosity and WSS in parent vessels of all aneurysms and in the sacs of B and C. However, large discrepancies in shear rate, viscosity and WSS among predictions by the various rheology models were found in the dome area of A where the flow was relatively stagnant. Here the Newtonian model predicted higher shear rate and WSS values and lower blood viscosity than the two non-Newtonian models. Conclusions The Newtonian fluid assumption can underestimate viscosity and overestimate shear rate and WSS in regions of stasis or slowly recirculating secondary vortices, typically found at the dome in elongated or complex-shaped saccular aneurysms as well as in aneurysms following endovascular treatment. Because low shear rates and low WSS in such flow conditions indicate a high propensity for thrombus formation and rupture, CFD based on the Newtonian assumption may underestimate the propensity of these events.

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

[2]  Aichi Chien,et al.  Patient-specific flow analysis of brain aneurysms at a single location: comparison of hemodynamic characteristics in small aneurysms , 2008, Medical & Biological Engineering & Computing.

[3]  C. Derdeyn,et al.  Mortality rates after subarachnoid hemorrhage: variations according to hospital case volume in 18 states. , 2003, Journal of neurosurgery.

[4]  R. Heros,et al.  Intracranial arterial aneurysm‐‐an update. , 1983, Stroke.

[5]  J. Rossmann,et al.  Effect of non-newtonian behavior on hemodynamics of cerebral aneurysms. , 2009, Journal of biomechanical engineering.

[6]  Adnan Siddiqui,et al.  Size Ratio Correlates With Intracranial Aneurysm Rupture Status: A Prospective Study , 2010, Stroke.

[7]  K Takakura,et al.  Clinical Analysis of Incidentally Discovered Unruptured Aneurysms , 1993, Stroke.

[8]  Zhijie Wang,et al.  Complex Hemodynamics at the Apex of an Arterial Bifurcation Induces Vascular Remodeling Resembling Cerebral Aneurysm Initiation , 2007, Stroke.

[9]  C. Putman,et al.  Characterization of cerebral aneurysms for assessing risk of rupture by using patient-specific computational hemodynamics models. , 2005, AJNR. American journal of neuroradiology.

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

[11]  P. Chaturani,et al.  A study of non-Newtonian aspects of blood flow through stenosed arteries and its applications in arterial diseases. , 1985, Biorheology.

[12]  Barbara M. Johnston,et al.  Non-Newtonian blood flow in human right coronary arteries: transient simulations. , 2006, Journal of biomechanics.

[13]  J. Mocco,et al.  Mapping vascular response to in vivo Hemodynamics: application to increased flow at the basilar terminus , 2010, Biomechanics and modeling in mechanobiology.

[14]  S. Renowden,et al.  Intra-Aneurysmal Thrombosis as a Possible Cause of Delayed Aneurysm Rupture after Flow-Diversion Treatment , 2010, American Journal of Neuroradiology.

[15]  J. Xiang,et al.  Alteration of intra-aneurysmal hemodynamics for flow diversion using enterprise and vision stents. , 2010, World neurosurgery.

[16]  T. Hori,et al.  Is the Aspect Ratio a Reliable Index for Predicting the Rupture of a Saccular Aneurysm? , 2001, Neurosurgery.

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

[18]  Ender A Finol,et al.  Blood flow dynamics in saccular aneurysm models of the basilar artery. , 2006, Journal of biomechanical engineering.

[19]  S Oka,et al.  Optimality principle in vascular bifurcation. , 1987, Biorheology.

[20]  J. Broderick,et al.  Initial and Recurrent Bleeding Are the Major Causes of Death Following Subarachnoid Hemorrhage , 1994, Stroke.

[21]  Sujan Dhar,et al.  INFLUENCE OF INTRACRANIAL ANEURYSM‐TO‐PARENT VESSEL SIZE RATIO ON HEMODYNAMICS AND IMPLICATION FOR RUPTURE: RESULTS FROM A VIRTUAL EXPERIMENTAL STUDY , 2009, Neurosurgery.

[22]  Elad I Levy,et al.  Size ratio for clinical assessment of intracranial aneurysm rupture risk , 2010, Neurological research.

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

[24]  L. Ge,et al.  Numerical Modeling of the Flow in Intracranial Aneurysms: Prediction of Regions Prone to Thrombus Formation , 2008, Annals of Biomedical Engineering.

[25]  K. Katada,et al.  Magnitude and Role of Wall Shear Stress on Cerebral Aneurysm: Computational Fluid Dynamic Study of 20 Middle Cerebral Artery Aneurysms , 2004, Stroke.

[26]  Elad I Levy,et al.  QUANTIFICATION OF HEMODYNAMIC CHANGES INDUCED BY VIRTUAL PLACEMENT OF MULTIPLE STENTS ACROSS A WIDE‐NECKED BASILAR TRUNK ANEURYSM , 2007, Neurosurgery.

[27]  Bu-Lang Gao,et al.  Aneurysm Inflow-Angle as a Discriminant for Rupture in Sidewall Cerebral Aneurysms: Morphometric and Computational Fluid Dynamic Analysis , 2010, Stroke.

[28]  H. Meng,et al.  Effects of arterial geometry on aneurysm growth: three-dimensional computational fluid dynamics study. , 2004, Journal of neurosurgery.

[29]  Alejandro F. Frangi,et al.  Efficient pipeline for image-based patient-specific analysis of cerebral aneurysm hemodynamics: technique and sensitivity , 2005, IEEE Transactions on Medical Imaging.

[30]  D. Holdsworth,et al.  A three-dimensional cerebrovascular flow phantom. , 1999, Medical physics.

[31]  F. N. van de Vosse,et al.  The influence of the non-Newtonian properties of blood on the flow in large arteries: steady flow in a carotid bifurcation model. , 1999, Journal of biomechanics.

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

[33]  D. Saloner,et al.  Numerical analysis of flow through a severely stenotic carotid artery bifurcation. , 2002, Journal of biomechanical engineering.

[34]  K. Takayama,et al.  A proposed parent vessel geometry-based categorization of saccular intracranial aneurysms: computational flow dynamics analysis of the risk factors for lesion rupture. , 2005, Journal of neurosurgery.