Evaluation of the influence of inlet boundary conditions on computational fluid dynamics for intracranial aneurysms: A virtual experiment

Modelling lower-limb peripheral arterial disease using clinically available datasets: impact of inflow boundary conditions on hemodynamic indices for restenosis prediction

BACKGROUND AND OBJECTIVES

The integration of hemodynamic markers as risk factors in restenosis prediction models for lower-limb peripheral arteries is hindered by fragmented clinical datasets. Computed tomography (CT) scans enable vessel geometry reconstruction and can be obtained at different times than the Doppler ultrasound (DUS) images, which provide information on blood flow velocity. Computational fluid dynamics (CFD) simulations allow the computation of near-wall hemodynamic indices, whose accuracy depends on the prescribed inlet boundary condition (BC), derived from the DUS images. This study aims to: (i) investigate the impact of different DUS-derived velocity waveforms on CFD results; (ii) test whether the same vessel areas, subjected to altered hemodynamics, can be detected independently of the applied inlet BC; (iii) suggest suitable DUS images to obtain reliable CFD results.

METHODS

CFD simulations were conducted on three patients treated with bypass surgery, using patient-specific DUS-derived inlet BCs recorded at either the same or different time points than the CT scan. The impact of the chosen inflow condition on bypass hemodynamics was assessed in terms of wall shear stress (WSS)-derived quantities. Patient-specific critical thresholds for the hemodynamic indices were applied to identify critical luminal areas and compare the results with a reference obtained with a DUS image acquired in close temporal proximity to the CT scan.

RESULTS

The main findings indicate that: (i) DUS-derived inlet velocity waveforms acquired at different time points than the CT scan led to statistically significantly different CFD results (p<0.001); (ii) the same luminal surface areas, exposed to low time-averaged WSS, could be identified independently of the applied inlet BCs; (iii) similar outcomes were observed for the other hemodynamic indices if the prescribed inlet velocity waveform had the same shape and comparable systolic acceleration time to the one recorded in close temporal proximity to the CT scan.

CONCLUSIONS

Despite a lack of standardised data collection for diseased lower-limb peripheral arteries, an accurate estimation of luminal areas subjected to altered near-wall hemodynamics is possible independently of the applied inlet BC. This holds if the applied inlet waveform shares some characteristics - derivable from the DUS report - as one matching the acquisition time of the CT scan.

Deep-learning-based image segmentation for image-based computational hemodynamic analysis of abdominal aortic aneurysms: a comparison study

Computational hemodynamics is increasingly being used to quantify hemodynamic characteristics in and around abdominal aortic aneurysms (AAA) in a patient-specific fashion. However, the time-consuming manual annotation hinders the clinical translation of computational hemodynamic analysis. Thus, we investigate the feasibility of using deep-learning-based image segmentation methods to reduce the time required for manual segmentation. Two of the latest deep-learning-based image segmentation methods, ARU-Net and CACU-Net, were used to test the feasibility of automated computer model creation for computational hemodynamic analysis. Morphological features and hemodynamic metrics of 30 computed tomography angiography (CTA) scans were compared between pre-dictions and manual models. The DICE score for both networks was 0.916, and the correlation value was above 0.95, indicating their ability to generate models comparable to human segmentation. The Bland-Altman analysis shows a good agreement between deep learning and manual segmentation results. Compared with manual (computational hemodynamics) model recreation, the time for automated computer model generation was significantly reduced (from ∼2 h to ∼10 min). Automated image segmentation can significantly reduce time expenses on the recreation of patient-specific AAA models. Moreover, our study showed that both CACU-Net and ARU-Net could accomplish AAA segmentation, and CACU-Net outperformed ARU-Net in terms of accuracy and time-saving.

How patient specific are patient-specific computational models of cerebral aneurysms? An overview of sources of error and variability

Computational modeling of cerebral aneurysms, derived from clinical 3D angiography, has become widespread over the past 15 years. While such "image-based" or "patient-specific" models have shown promise for the assessment of rupture risk, much debate remains about their reliability in light of necessary modeling assumptions and incomplete or uncertain model input parameters derived from the clinic. The aims of this review were to walk through the various steps of this so-called patient-specific modeling pipeline and to highlight evidence supporting those steps that we can or cannot rely on. The relative importance of the different sources of error and variability on hemodynamic predictions is summarized, with recommendations to standardize for those that can be avoided and to pay closer attention those to that cannot.

How patient-specific do internal carotid artery inflow rates need to be for computational fluid dynamics of cerebral aneurysms?

BACKGROUND

Computational fluid dynamics (CFD) has become a popular tool for studying 'patient-specific' blood flow dynamics in cerebral aneurysms; however, rarely are the inflow boundary conditions patient-specific. We aimed to test the impact of widespread reliance on generalized inflow rates.

METHODS

Internal carotid artery (ICA) flow rates were measured via 2D cine phase-contrast MRI for 24 patients scheduled for endovascular therapy of an ICA aneurysm. CFD models were constructed from 3D rotational angiography, and pulsatile inflow rates imposed as measured by MRI or estimated using an average older-adult ICA flow waveform shape scaled by a cycle-average flow rate (Q_avg) derived from the patient's ICA cross-sectional area via an assumed inlet velocity.

RESULTS

There was good overall qualitative agreement in the magnitudes and spatial distributions of time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), and spectral power index (SPI) using generalized versus patient-specific inflows. Sac-averaged quantities showed moderate to good correlations: R²=0.54 (TAWSS), 0.80 (OSI), and 0.68 (SPI). Using patient-specific Q_avg to scale the generalized waveform shape resulted in near-perfect agreement for TAWSS, and reduced bias, but not scatter, for SPI. Patient-specific waveform had an impact only on OSI correlations, which improved to R²=0.93.

CONCLUSIONS

Aneurysm CFD demonstrates the ability to stratify cases by nominal hemodynamic 'risk' factors when employing an age- and vascular-territory-specific recipe for generalized inflow rates. Q_avg has a greater influence than waveform shape, suggesting some improvement could be achieved by including measurement of patient-specific Q_avg into aneurysm imaging protocols.

Discrimination of intracranial aneurysm rupture status: patient-specific inflow boundary may not be a must-have condition in hemodynamic simulations

BACKGROUND

Computational fluid dynamics (CFD) are important in evaluating the hemodynamics of intracranial aneurysm rupture, and the setting of inflow boundary conditions is critical. We evaluated intracranial aneurysm hemodynamics based on generalized versus patient-specific inflow boundary conditions to examine the effect of different hemodynamic results on the discrimination of intracranial aneurysm rupture status.

METHODS

We enrolled 148 patients with 156 intracranial aneurysms. For each included aneurysm, we performed CFD simulation once based on patient-specific and once based on generalized inflow boundary conditions. First, we compared the hemodynamics of intracranial aneurysms based on different inflow boundary conditions. Then, we divided the included aneurysms into a ruptured and unruptured group and compared the hemodynamics between the two groups under patient-specific and generalized inflow boundary conditions.

RESULTS

For the hemodynamic parameters using specific inflow boundary conditions, more complex flow (p = 0.002), larger minimum WSS (p = 0.024), lower maximum low WSS area (LSA) (p = 0.038), and oscillatory shear index (p = 0.002) were found. Furthermore, we compared the hemodynamics between ruptured and unruptured groups based on different inflow boundary conditions. We found that the significant hemodynamic parameters associated with rupture status were the same, including the proportion of aneurysms with flow complex and unstable flow and the minimum and maximum of LSA (p = 0.011, p = 0.003, p = 0.001 and p = 0.004, respectively).

CONCLUSION

Patient-specific and generalized inflow boundary conditions of aneurysmal hemodynamics resulted in significant differences. However, the significant parameters associated with rupture status were the same in both conditions, indicating that patient-specific inflow boundary conditions may not be necessary for predicting rupture risk.

What does computational fluid dynamics tell us about intracranial aneurysms? A meta-analysis and critical review

Despite the plethora of published studies on intracranial aneurysms (IAs) hemodynamic using computational fluid dynamics (CFD), limited progress has been made towards understanding the complex physics and biology underlying IA pathophysiology. Guided by 1733 published papers, we review and discuss the contemporary IA hemodynamics paradigm established through two decades of IA CFD simulations. We have traced the historical origins of simplified CFD models which impede the progress of comprehending IA pathology. We also delve into the debate concerning the Newtonian fluid assumption used to represent blood flow computationally. We evidently demonstrate that the Newtonian assumption, used in almost 90% of studies, might be insufficient to describe IA hemodynamics. In addition, some fundamental properties of the Navier-Stokes equation are revisited in supplementary material to highlight some widely spread misconceptions regarding wall shear stress (WSS) and its derivatives. Conclusively, our study draws a roadmap for next-generation IA CFD models to help researchers investigate the pathophysiology of IAs.

Hemodynamics of Cerebral Aneurysms: Connecting Medical Imaging and Biomechanical Analysis

In the last two decades, numerous studies have conducted patient-specific computations of blood flow dynamics in cerebral aneurysms and reported correlations between various hemodynamic metrics and aneurysmal disease progression or treatment outcomes. Nevertheless, intra-aneurysmal flow analysis has not been adopted in current clinical practice, and hemodynamic factors usually are not considered in clinical decision making. This review presents the state of the art in cerebral aneurysm imaging and image-based modeling, discussing the advantages and limitations of each approach and focusing on the translational value of hemodynamic analysis. Combining imaging and modeling data obtained from different flow modalities can improve the accuracy and fidelity of resulting velocity fields and flow-derived factors that are thought to affect aneurysmal disease progression. It is expected that predictive models utilizing hemodynamic factors in combination with patient medical history and morphological data will outperform current risk scores and treatment guidelines. Possible future directions include novel approaches enabling data assimilation and multimodality analysis of cerebral aneurysm hemodynamics.

COMPUTATIONAL STUDY ON THE RUPTURE RISK IN REAL CEREBRAL ANEURYSMS WITH GEOMETRICAL AND FLUID-MECHANICAL PARAMETERS USING FSI SIMULATIONS AND MACHINE LEARNING ALGORITHMS

Fluid-mechanical and morphological parameters are recognized as major factors in the rupture risk of human aneurysms. On the other hand, it is well known that a lot of machine learning tools are available to study a variety of problems in many fields. In this work, fluid–structure interaction (FSI) simulations were carried out to examine a database of 60 real saccular cerebral aneurysms (30 ruptured and 30 unruptured) using reconstructions by angiography images. With the results of the simulations and geometric analyses, we studied the analysis of variance (ANOVA) statistic test in many variables and we obtained that aspect ratio (AR), bottleneck factor (BNF), maximum height of the aneurysms (MH), relative residence time (RRT), Womersley number (WN) and Von-Mises strain (VMS) are statically significant and good predictors for the models. In consequence, these ones were used in five machine learning algorithms to determine the rupture risk predictions of the aneurysms, where the adaptative boosting (AdaBoost) was calculated with the highest area under the curve (AUC) in the receiver operating characteristic (ROC) curve (AUC 0.944).

The effect of endothelial nitric oxide synthase on the hemodynamics and wall mechanics in murine arteriovenous fistulas

Creation of a hemodialysis arteriovenous fistula (AVF) causes aberrant vascular mechanics at and near the AVF anastomosis. When inadequately regulated, these aberrant mechanical factors may impede AVF lumen expansion to cause AVF maturation failure, a significant clinical problem with no effective treatments. The endothelial nitric oxide synthase (NOS3) system is crucial for vascular health and function, but its effect on AVF maturation has not been fully characterized. We hypothesize that NOS3 promotes AVF maturation by regulating local vascular mechanics following AVF creation. Here we report the first MRI-based fluid-structure interaction (FSI) study in a murine AVF model using three mouse strains: NOS3 overexpression (NOS3 OE) and knockout (NOS3−/−) on C57BL/6 background, with C57BL/6 as the wild-type control (NOS3+/+). When compared to NOS3+/+ and NOS3−/−, AVFs in the OE mice had larger lumen area. AVFs in the OE mice also had smoother blood flow streamlines, as well as lower blood shear stress at the wall, blood vorticity, inner wall circumferential stretch, and radial wall thinning at the anastomosis. Our results demonstrate that overexpression of NOS3 resulted in distinct hemodynamic and wall mechanical profiles associated with favorable AVF remodeling. Enhancing NOS3 expression may be a potential therapeutic approach for promoting AVF maturation.

STUDY ON THE RELATIONSHIP BETWEEN WALL SHEAR STRESS AND ASPECT RATIO OF CEREBRAL ANEURYSMS WITH DIFFERENT PRESSURE DIFFERENCES USING CFD SIMULATIONS

CFD simulations were performed for 60 human cerebral aneurysms (30 previously ruptured and 30 previously unruptured) to study the behavior of the time-averaged wall shear stress (TAWSS) with respect to the aspect ratio (AR), implementing a set of low, normal, and high-pressure differences between the inlet and the outlets of each artery. It is well known that there exists a direct relationship between TAWSS and the rupture. In this investigation, we presented an important result because the condition of the pressure among the branches and the AR may be measured in any patient, then a slope may be associated, and finally a TAWSS may be estimated. We found that when the pressure difference increased, the absolute slopes between TAWSS and AR increased as well. Also, the magnitude of the slope in the previously unruptured aneurysms was 4.7 times the slope in the previously ruptured aneurysms. On the other hand, TAWSS was higher in the previously unruptured aneurysm than previously ruptured aneurysms due to the unruptured aneurysms that have a smaller surface area. Furthermore, we analyzed the relationship between TAWSS and other geometric parameters of the aneurysm, such as bottleneck and non-sphericity index; however, no correlation was found for either cases.

Assessment of boundary conditions for CFD simulation in human carotid artery

Computational fluid dynamics (CFD) is an increasingly used method for investigation of hemodynamic parameters and their alterations under pathological conditions, which are important indicators for diagnosis of cardiovascular disease. In hemodynamic simulation models, the employment of appropriate boundary conditions (BCs) determines the computational accuracy of the CFD simulation in comparison with pressure and velocity measurements. In this study, we have first assessed the influence of inlet boundary conditions on hemodynamic CFD simulations. We selected two typical patients suspected of carotid artery disease, with mild stenosis and severe stenosis. Both patients underwent digital subtraction angiography (DSA), magnetic resonance angiography, and the invasive pressure guide wire measured pressure profile. We have performed computational experiments to (1) study the hemodynamic simulation outcomes of distributions of wall shear stress, pressure, pressure gradient and (2) determine the differences in hemodynamic performances caused by inlet BCs derived from DSA and Womersley analytical solution. Our study has found that the difference is related to the severity of the stenosis; the greater the stenosis, the more the difference ensues. Further, in our study, the two typical subjects with invasively measured pressure profile and thirty subjects with ultrasound Doppler velocimeter (UDV) measurement served as the criteria to evaluate the hemodynamic outcomes of wall shear stress, pressure, pressure gradient and velocity due to different outlet BCs based on the Windkessel model, structured-tree model, and fully developed flow model. According to the pressure profiles, the fully developed model appeared to have more fluctuations compared with the other two models. The Windkessel model had more singularities before convergence. The three outlet BCs models also showed good correlation with the UDV measurement, while the Windkessel model appeared to be slightly better ([Formula: see text]). The structured-tree model was seen to have the best performance in terms of available computational cost and accuracy. The results of our numerical simulation and the good correlation with the computed pressure and velocity with their measurements have highlighted the effectiveness of CFD simulation in patient-specific human carotid artery with suspected stenosis.

The effect of inlet and outlet boundary conditions in image-based CFD modeling of aortic flow

Background

Computational modeling of cardiovascular flow is a growing and useful field, but such simulations usually require the researcher to guess the flow’s inlet and outlet conditions since they are difficult and expensive to measure. It is critical to determine the amount of uncertainty introduced by these assumptions in order to evaluate the degree to which cardiovascular flow simulations are accurate. Our work begins to address this question by examining the sensitivity of flow to several different assumed velocity inlet and outlet conditions in a patient-specific aorta model.

Methods

We examined the differences between plug flow, parabolic flow, linear shear flows, skewed cubic flow profiles, and Womersley flow at the inlet. Only the shape of the inlet velocity profile was varied—all other parameters were identical among these simulations. Secondary flow in the form of a counter-rotating pair of vortices was also added to parabolic axial flow to study its effect on the solution. In addition, we examined the differences between two-element Windkessel, three element Windkessel and the outflow boundary conditions. In these simulations, only the outlet boundary condition was varied.

Results

The results show axial and in-plane velocities are considerably different close to the inlet for the cases with different inlet velocity profile shapes. However, the solutions are qualitatively similar beyond 1.75D, where D is the inlet diameter. This trend is also observed in other quantities such as pressure and wall shear stress. Normalized root-mean-square deviation, a measure of axial velocity magnitude differences between the different cases, generally decreases along the streamwise coordinate. The linear shear inlet velocity boundary condition and plug velocity boundary condition solution exhibit the highest time-averaged wall shear stress, approximately \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$8\%$$\end{document}8% higher than the parabolic inlet velocity boundary condition. Upstream of 1D from the inlet, adding secondary flow has a significant impact on temporal wall shear stress distributions. This is especially observable during diastole, when integrated wall shear stress magnitude varies about \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$26\%$$\end{document}26% between simulations with and without secondary flow. The results from the outlet boundary condition study show the Windkessel models differ from the outflow boundary condition by as much as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$18\%$$\end{document}18% in terms of time-averaged wall shear stress. Furthermore, normalized root-mean-square deviation of axial velocity magnitude, a measure of deviation between Windkessel and the outflow boundary condition, increases along the streamwise coordinate indicating larger variations near outlets.

Conclusion

It was found that the selection of inlet velocity conditions significantly affects only the flow region close to the inlet of the aorta. Beyond two diameters distal to the inlet, differences in flow solution are small. Although additional studies must be performed to verify this result, the data suggest that it is important to use patient-specific inlet conditions primarily if the researcher is concerned with the details of the flow very close to the inlet. Similarly, the selection of outlet conditions significantly affects the flow in the vicinity of the outlets. Upstream of five diameters proximal to the outlet, deviations between the outlet boundary conditions examined are insignificant. Although the inlet and outlet conditions only affect the flow significantly in their respective neighborhoods, our study indicates that outlet conditions influence a larger percentage of the solution domain.

Virtual endovascular treatment of intracranial aneurysms: models and uncertainty

Virtual endovascular treatment models (VETMs) have been developed with the view to aid interventional neuroradiologists and neurosurgeons to pre-operatively analyze the comparative efficacy and safety of endovascular treatments for intracranial aneurysms. Based on the current state of VETMs in aneurysm rupture risk stratification and in patient-specific prediction of treatment outcomes, we argue there is a need to go beyond personalized biomechanical flow modeling assuming deterministic parameters and error-free measurements. The mechanobiological effects associated with blood clot formation are important factors in therapeutic decision making and models of post-treatment intra-aneurysmal biology and biochemistry should be linked to the purely hemodynamic models to improve the predictive power of current VETMs. The influence of model and parameter uncertainties associated to each component of a VETM is, where feasible, quantified via a random-effects meta-analysis of the literature. This allows estimating the pooled effect size of these uncertainties on aneurysmal wall shear stress. From such meta-analyses, two main sources of uncertainty emerge where research efforts have so far been limited: (1) vascular wall distensibility, and (2) intra/intersubject systemic flow variations. In the future, we suggest that current deterministic computational simulations need to be extended with strategies for uncertainty mitigation, uncertainty exploration, and sensitivity reduction techniques. WIREs Syst Biol Med 2017, 9:e1385. doi: 10.1002/wsbm.1385 For further resources related to this article, please visit the WIREs website.

High resolution hemodynamic profiling of murine arteriovenous fistula using magnetic resonance imaging and computational fluid dynamics

Background

Arteriovenous fistula (AVF) maturation failure remains a major cause of morbidity and mortality in hemodialysis patients. The two major etiologies of AVF maturation failure are early neointimal hyperplasia development and persistent inadequate outward remodeling. Although hemodynamic changes following AVF creation may impact AVF remodeling and contribute to neointimal hyperplasia development and impaired outward remodeling, detailed AVF hemodynamics are not yet fully known. Since murine AVF models are valuable tools for investigating the pathophysiology of AVF maturation failure, there is a need for a new approach that allows the hemodynamic characterization of murine AVF at high resolutions.

Methods

This methods paper presents a magnetic resonance imaging (MRI)-based computational fluid dynamic (CFD) method that we developed to rigorously quantify the evolving hemodynamic environment in murine AVF. The lumen geometry of the entire murine AVF was reconstructed from high resolution, non-contrast 2D T2-weighted fast spin echo MRI sequence, and the flow rates of the AVF inflow and outflow were extracted from a gradient echo velocity mapping sequence. Using these MRI-obtained lumen geometry and inflow information, CFD modeling was performed and used to calculate blood flow velocity and hemodynamic factors at high resolutions (on the order of 0.5 μm spatially and 0.1 ms temporally) throughout the entire AVF lumen. We investigated both the wall properties (including wall shear stress (WSS), wall shear stress spatial gradient, and oscillatory shear index (OSI)) and the volumetric properties (including vorticity, helicity, and Q-criterion).

Results

Our results demonstrate increases in AVF flow velocity, WSS, spatial WSS gradient, and OSI within 3 weeks post-AVF creation when compared to pre-surgery. We also observed post-operative increases in flow disturbances and vortices, as indicated by increased vorticity, helicity, and Q-criterion.

Conclusions

This novel protocol will enable us to undertake future mechanistic studies to delineate the relationship between hemodynamics and AVF development and characterize biological mechanisms that regulate local hemodynamic factors in transgenic murine AVF models.

Electronic supplementary material

The online version of this article (doi:10.1186/s12976-017-0053-x) contains supplementary material, which is available to authorized users.

EFFECTS OF HYPERTENSION AND PRESSURE GRADIENT IN A HUMAN CEREBRAL ANEURYSM USING FLUID STRUCTURE INTERACTION SIMULATIONS

Fluid–structure interaction (FSI) simulations were carried out in a human cerebral aneurysm model with the objective of quantifying the effects of hypertension and pressure gradient on the behavior of fluid and solid mechanics. Six FSI simulations were conducted using a hyperelastic Mooney–Rivlin model. Important differences in wall shear stress (WSS), wall displacements, and effective von Mises stress are reported. The hypertension increases wall stress and displacements in the aneurysm region; however, the effects of hypertension on the hemodynamics in the aneurysm region were small. The pressure gradient affects the WSS in the aneurysm and also the displacement and wall stress on the aneurysm. Maximum wall stress with hypertension in the range of rupture strength was found.

4D Flow MRI in Neuroradiology: Techniques and Applications

Assessment of the intracranial flow is important for the understanding and management of cerebral vascular diseases. From brain aneurysms and arteriovenous malformations lesions to intracranial and cervical stenosis, the appraisal of the blood flow can be crucial and influence positively on patients' management. The determination of the intracranial hemodynamics and the collateral pattern seems to play to a major role in the management of these lesions. 4D flow magnetic resonance imaging is a noninvasive phase contrast derived method that has been developed and applied in neurovascular diseases. It has a great potential if followed by further technical improvements and comprehensive and systematic clinical studies.

Uncertainty quantification of wall shear stress in intracranial aneurysms using a data-driven statistical model of systemic blood flow variability

Adverse wall shear stress (WSS) patterns are known to play a key role in the localisation, formation, and progression of intracranial aneurysms (IAs). Complex region-specific and time-varying aneurysmal WSS patterns depend both on vascular morphology as well as on variable systemic flow conditions. Computational fluid dynamics (CFD) has been proposed for characterising WSS patterns in IAs; however, CFD simulations often rely on deterministic boundary conditions that are not representative of the actual variations in blood flow. We develop a data-driven statistical model of internal carotid artery (ICA) flow, which is used to generate a virtual population of waveforms used as inlet boundary conditions in CFD simulations. This allows the statistics of the resulting aneurysmal WSS distributions to be computed. It is observed that ICA waveform variations have limited influence on the time-averaged WSS (TAWSS) on the IA surface. In contrast, in regions where the flow is locally highly multidirectional, WSS directionality and harmonic content are strongly affected by the ICA flow waveform. As a consequence, we argue that the effect of blood flow variability should be explicitly considered in CFD-based IA rupture assessment to prevent confounding the conclusions.

Sensitivity of flow patterns in aneurysms on the anterior communicating artery to anatomic variations of the cerebral arterial network

Recent studies raised increasing concern about the reliability of computer models in reproducing in vivo hemodynamics in cerebral aneurysms. Boundary condition problem is among the most frequently addressed issues since three-dimensional (3-D) modeling is usually restricted to local arterial segments. The present study focused on aneurysms on the anterior communicating artery (ACoA) which represent a large subgroup of detected cerebral aneurysms and, in particular, have a relatively high risk of rupture compared to aneurysms located in other regions. The sensitivity of blood flows in three ACoA aneurysms to boundary conditions was investigated using 3-D hemodynamic models. The boundary conditions of the 3-D models were predicted by a one-dimensional (1-D) model of the cerebral arterial network. The parameters of the 1-D model were assigned based respectively on population-averaged data and patient-specific data derived from medical images, yielding a population-generic model and a patient-specific model. In addition, particle image velocimetry (PIV) experiments were performed to validate the code used to simulate intra-aneurysmal blood flows. Obtained results showed that switching the boundary conditions of the aneurysm models from population-generic ones to patient-specific ones led to pronounced changes in simulated intra-aneurysmal flow patterns in terms of vortex structure, impingement region and the magnitude and spatial distribution of wall shear stress and oscillatory shear index. In particular, the way and the degree in which hemodynamic quantities are influenced by boundary conditions exhibited pronounced inter-patient variability. In summary, our study underlines the importance of patient-specific treatment of boundary conditions in model studies focusing on ACoA aneurysms.

Computational fluid dynamics analysis of flow reduction induced by flow-diverting stents in intracranial aneurysms: a patient-unspecific hemodynamics change perspective

Background and purpose

Flow-diverter stents (FDSs) have been used effectively to treat large neck and complex saccular aneurysms on the anterior carotid circulation. Intra-aneurysmal flow reduction induces progressive aneurysm thrombosis in most patients. Understanding the degree of flow modification necessary to induce complete aneurysm occlusion among patients with considerable hemodynamics variability may be important for treatment planning.

Materials and methods

Patients with incidental intracranial saccular aneurysms who underwent FDS endovascular procedures were included and studied for a 12 months’ follow-up period. We used computational fluid dynamics on patient-specific geometries from 3D rotational angiography without and with virtual stent placement and thus compared intra-aneurysmal hemodynamic problems. Receiver operating characteristic analysis was used to estimate the stent:no-stent minimum hemodynamic ratio thresholds that significantly (p≤0.05) determined the condition necessary for long-term (12 months) aneurysm occlusion.

Results

We included 12 consecutive patients with sidewall aneurysms located in the internal carotid or vertebral artery. The measured porosity of the 12 deployed virtual FDSs was 83±3% (mean±SD). Nine aneurysms were occluded during the 12 months’ follow-up, whereas three were not. A significant (p=0.05) area under the curve (AUC) was found for spatiotemporal mean velocity reduction in the aneurysms: AUC=0.889±0.113 (mean±SD) corresponding to a minimum velocity reduction threshold of 0.353 for occlusion to occur. The 95% CI of the AUC was 0.66 to 1.00. The sensitivity and specificity of the method were ∼99% and ∼67%, respectively. For both wall shear stress and pressure reductions in aneurysms no thresholds could be determined: AUC=0.63±0.16 (p=0.518) and 0.67±0.165 (p=0.405), respectively.

Conclusions

For successful FDS treatment the post-stent average velocity in sidewall intracranial aneurysms must be reduced by at least one-third from the initial pre-stent conditions.

Numerical Study of Cerebroarterial Hemodynamic Changes Following Carotid Artery Operation: A Comparison Between Multiscale Modeling and Stand-Alone Three-Dimensional Modeling

Free outflow boundary conditions have been widely adopted in hemodynamic model studies, they, however, intrinsically lack the ability to account for the regulatory mechanisms of systemic hemodynamics and hence carry a risk of producing incorrect results when applied to vascular segments with multiple outlets. In the present study, we developed a multiscale model capable of incorporating global cardiovascular properties into the simulation of blood flows in local vascular segments. The multiscale model was constructed by coupling a three-dimensional (3D) model of local arterial segments with a zero-one-dimensional (0-1-D) model of the cardiovascular system. Numerical validation based on an idealized model demonstrated the ability of the multiscale model to preserve reasonable pressure/flow wave transmission among different models. The multiscale model was further calibrated with clinical data to simulate cerebroarterial hemodynamics in a patient undergoing carotid artery operation. The results showed pronounced hemodynamic changes in the cerebral circulation following the operation. Additional numerical experiments revealed that a stand-alone 3D model with free outflow conditions failed to reproduce the results obtained by the multiscale model. These results demonstrated the potential advantage of multiscale modeling over single-scale modeling in patient-specific hemodynamic studies. Due to the fact that the present study was limited to a single patient, studies on more patients would be required to further confirm the findings.

Influence of Parent Artery Segmentation and Boundary Conditions on Hemodynamic Characteristics of Intracranial Aneurysms

PURPOSE

The purpose of this study is to explore the influence of segmentation of the upstream and downstream parent artery and hemodynamic boundary conditions (BCs) on the evaluated hemodynamic factors for the computational fluid dynamics (CFD) analysis of intracranial aneurysms.

MATERIALS AND METHODS

Three dimensional patient-specific aneurysm models were analyzed by applying various combinations of inlet and outlet BCs. Hemodynamic factors such as velocity pattern, streamline, wall shear stress, and oscillatory shear index at the systolic time were visualized and compared among the different cases.

RESULTS

Hemodynamic factors were significantly affected by the inlet BCs while there was little influence of the outlet BCs. When the inlet length was relatively short, different inlet BCs showed different hemodynamic factors and the calculated hemodynamic factors were also dependent on the inlet length. However, when the inlet length (L) was long enough (L>20D, where D is the diameter of inlet section), the hemodynamic factors became similar regardless of the inlet BCs and lengths. The error due to different inlet BCs was negligible. The effect of the outlet length on the hemodynamic factors was similar to that of the inlet length.

CONCLUSION

Simulated hemodynamic factors are highly sensitive to inlet BCs and upstream parent artery segmentation. The results of this work can provide an insight into how to build models and to apply BCs for more accurate estimation of hemodynamic factors from CFD simulations of intracranial aneurysms.

Cerebral Blood Flow in a Healthy Circle of Willis and Two Intracranial Aneurysms: Computational Fluid Dynamics Versus Four-Dimensional Phase-Contrast Magnetic Resonance Imaging

Computational fluid dynamics (CFD) opens up multiple opportunities to investigate the hemodynamics of the human vascular system. However, due to numerous assumptions the acceptance of CFD among physicians is still limited in practice and validation through comparison is mandatory. Time-dependent quantitative phase-contrast magnetic resonance imaging PC-MRI measurements in a healthy volunteer and two intracranial aneurysms were carried out at 3 and 7 Tesla. Based on the acquired images, three-dimensional (3D) models of the aneurysms were reconstructed and used for the numerical simulations. Flow information from the MR measurements were applied as boundary conditions. The four-dimensional (4D) velocity fields obtained by CFD and MRI were qualitatively as well as quantitatively compared including cut planes and vector analyses. For all cases a high similarity of the velocity patterns was observed. Additionally, the quantitative analysis revealed a good agreement between CFD and MRI. Deviations were caused by minor differences between the reconstructed vessel models and the actual lumen. The comparisons between diastole and systole indicate that relative differences between MRI and CFD are intensified with increasing velocity. The findings of this study lead to the conclusion that CFD and MRI agree well in predicting intracranial velocities when realistic geometries and boundary conditions are provided. Due to the considerably higher temporal and spatial resolution of CFD compared to MRI, complex flow patterns can be further investigated in order to evaluate their role with respect to aneurysm formation or rupture. Nevertheless, special care is required regarding the vessel reconstruction since the geometry has a major impact on the subsequent numerical results.