Better Than Nothing: A Rational Approach for Minimizing the Impact of Outflow Strategy on Cerebrovascular Simulations

BACKGROUND AND PURPOSE: Computational fluid dynamics simulations of neurovascular diseases are impacted by various modeling assumptions and uncertainties, including outlet boundary conditions. Many studies of intracranial aneurysms, for example, assume zero pressure at all outlets, often the default (“do-nothing”) strategy, with no physiological basis. Others divide outflow according to the outlet diameters cubed, nominally based on the more physiological Murray's law but still susceptible to subjective choices about the segmented model extent. Here we demonstrate the limitations and impact of these outflow strategies, against a novel “splitting” method introduced here. MATERIALS AND METHODS: With our method, the segmented lumen is split into its constituent bifurcations, where flow divisions are estimated locally using a power law. Together these provide the global outflow rate boundary conditions. The impact of outflow strategy on flow rates was tested for 70 cases of MCA aneurysm with 0D simulations. The impact on hemodynamic indices used for rupture status assessment was tested for 10 cases with 3D simulations. RESULTS: Differences in flow rates among the various strategies were up to 70%, with a non-negligible impact on average and oscillatory wall shear stresses in some cases. Murray-law and splitting methods gave flow rates closest to physiological values reported in the literature; however, only the splitting method was insensitive to arbitrary truncation of the model extent. CONCLUSIONS: Cerebrovascular simulations can depend strongly on the outflow strategy. The default zero-pressure method should be avoided in favor of Murray-law or splitting methods, the latter being released as an open-source tool to encourage the standardization of outflow strategies.

[1]  Yiqian,et al.  Hemodynamic Differences Between Unruptured and Ruptured Intracranial Aneurysms During Observation , 2012 .

[2]  D F Kallmes,et al.  Artery Length Sensitivity in Patient-Specific Cerebral Aneurysm Simulations , 2015, American Journal of Neuroradiology.

[3]  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.

[4]  Ivo G. H. Jansen,et al.  Generalized versus Patient-Specific Inflow Boundary Conditions in Computational Fluid Dynamics Simulations of Cerebral Aneurysmal Hemodynamics , 2014, American Journal of Neuroradiology.

[5]  T. Wonnacott,et al.  Relation between diameter and flow in major branches of the arch of the aorta. , 1992, Journal of biomechanics.

[6]  G. Duckwiler,et al.  Sensitivity of patient-specific numerical simulation of cerebal aneurysm hemodynamics to inflow boundary conditions. , 2007, Journal of neurosurgery.

[7]  Bongjae Chung,et al.  CFD for Evaluation and Treatment Planning of Aneurysms: Review of Proposed Clinical Uses and Their Challenges , 2014, Annals of Biomedical Engineering.

[8]  L. Antiga,et al.  Intracranial Aneurysm Neck Size Overestimation with 3D Rotational Angiography: The Impact on Intra-Aneurysmal Hemodynamics Simulated with Computational Fluid Dynamics , 2013, American Journal of Neuroradiology.

[9]  D A Steinman,et al.  Non‐Newtonian versus numerical rheology: Practical impact of shear‐thinning on the prediction of stable and unstable flows in intracranial aneurysms , 2017, International journal for numerical methods in biomedical engineering.

[10]  G. Karniadakis,et al.  Outflow Boundary Conditions for Arterial Networks with Multiple Outlets , 2008, Annals of Biomedical Engineering.

[11]  David A Steinman,et al.  A dual-pressure boundary condition for use in simulations of bifurcating conduits. , 2002, Journal of biomechanical engineering.

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

[13]  C Chnafa,et al.  Improved reduced-order modelling of cerebrovascular flow distribution by accounting for arterial bifurcation pressure drops. , 2017, Journal of biomechanics.

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

[15]  Simone Manini,et al.  pyNS: An Open-Source Framework for 0D Haemodynamic Modelling , 2014, Annals of Biomedical Engineering.

[16]  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.

[17]  Alejandro F. Frangi,et al.  Computational Hemodynamics in Cerebral Aneurysms: The Effects of Modeled Versus Measured Boundary Conditions , 2010, Annals of Biomedical Engineering.

[18]  A. Moura,et al.  Sensitivity to outflow boundary conditions and level of geometry description for a cerebral aneurysm , 2012, International journal for numerical methods in biomedical engineering.

[19]  Bo Carlberg,et al.  Aging alters the dampening of pulsatile blood flow in cerebral arteries , 2016, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[20]  三浦 洋一,et al.  Low wall shear stress is independently associated with the rupture status of middle cerebral artery aneurysms , 2013 .

[21]  C Chnafa,et al.  Vessel calibre and flow splitting relationships at the internal carotid artery terminal bifurcation , 2017, Physiological measurement.

[22]  Alberto Avolio,et al.  Image segmentation methods for intracranial aneurysm haemodynamic research. , 2014, Journal of biomechanics.

[23]  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.

[24]  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.

[25]  D A Steinman,et al.  The Computational Fluid Dynamics Rupture Challenge 2013—Phase I: Prediction of Rupture Status in Intracranial Aneurysms , 2015, American Journal of Neuroradiology.

[26]  C D Murray,et al.  The Physiological Principle of Minimum Work: I. The Vascular System and the Cost of Blood Volume. , 1926, Proceedings of the National Academy of Sciences of the United States of America.

[27]  David A. Steinman,et al.  Automatic Neck Plane Detection and 3D Geometric Characterization of Aneurysmal Sacs , 2012, Annals of Biomedical Engineering.

[28]  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.

[29]  H. Schirmer,et al.  Bifurcation geometry and the presence of cerebral artery aneurysms. , 2004, Journal of neurosurgery.

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

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