Wall shear stress at the initiation site of cerebral aneurysms

Hemodynamics are believed to play an important role in the initiation of cerebral aneurysms. In particular, studies have focused on wall shear stress (WSS), which is a key regulator of vascular biology and pathology. In line with the observation that aneurysms predominantly occur at regions of high WSS, such as bifurcation apices or outer walls of vascular bends, correlations have been found between the aneurysm initiation site and high WSS. The aim of our study was to analyze the WSS field at an aneurysm initiation site that was neither a bifurcation apex nor the outer wall of a vascular bend. Ten cases with aneurysms on the A1 segment of the anterior cerebral artery were analyzed and compared with ten controls. Aneurysms were virtually removed from the vascular models of the cases to mimic the pre-aneurysm geometry. Computational fluid dynamics (CFD) simulations were created to assess the magnitude, gradient, multidirectionality, and pulsatility of the WSS. To aid the inter-subject comparison of hemodynamic variables, we mapped the branch surfaces onto a two-dimensional parametric space. This approach made it possible to view the whole branch at once for qualitative evaluation. It also allowed us to empirically define a patch for quantitative analysis, which was consistent among subjects and encapsulated the aneurysm initiation sites in our dataset. To test the sensitivity of our results, CFD simulations were repeated with a second independent observer virtually removing the aneurysms and with a 20 % higher flow rate at the inlet. We found that branches harboring aneurysms were characterized by high WSS and high WSS gradients. Among all assessed variables, the aneurysm initiation site most consistently coincided with peaks of temporal variation in the WSS magnitude.     

High-resolution computational fluid dynamics detects flow instabilities in the carotid siphon: Implications for aneurysm initiation and rupture?

The carotid siphon is by nature a tortuous vessel segment with sharp bends and large area variations, and of relevance to the study of intracranial aneurysm initiation and rupture. The aim of this paper was to determine whether the siphon might harbor flow instabilities, if care is taken to resolve them. This study focused on five consecutive internal carotid artery (ICA) aneurysm cases from the open-source Aneurisk dataset. The aneurysm, always downstream of the siphon, was digitally removed using previously developed and verified tools. Computational fluid dynamic (CFD) models included long cervical segments upstream, and middle and anterior cerebral arteries downstream. High-resolution pulsatile simulations were performed using the equivalent of ~24$~24$ million linear tetrahedra on average (range 16-32 M) and 30,000 time-steps/cycle. Two of the five cases were laminar with mild flow instabilities right after peak systole. One of the cases experienced strong periodic vortex shedding at a frequency of around 100 Hz. The remaining two cases harbored higher frequency flow instabilities and complex 3D vortical structures, extending to the cerebral arteries downstream. Our findings suggest that the carotid siphon, a conduit to the majority of anterior intracranial aneurysms, may experience flow instabilities, consistent with in vitro reports, but seemingly at odds with the majority of CFD studies, which have been done at lower resolutions. This has obvious implications for elucidating the forces involved in aneurysm initiation; and propagation of flow instabilities into ICA or downstream aneurysms could also impact understanding of the forces involved in aneurysm rupture.     

Flow patterns and wall shear stress distribution in human internal carotid arteries: the geometric effect on the risk for stenoses

It has been widely observed that atherosclerotic stenosis occurs at sites with complex hemodynamics, such as arteries with high curvature or bifurcations. These regions usually have very low or highly oscillatory wall shear stress (WSS). In the present study, 3D sinusoidally pulsatile blood flow through the models of internal carotid artery (ICA) with different geometries was investigated with computational simulation. Three preferred sites of stenoses were found along the carotid siphon with low and highly oscillatory WSS. The risk for stenoses at these sites was scaled with the values of time-averaged WSS and oscillating shear index (OSI). The local risk for stenoses at every preferred site of stenoses was found different between 3 types of ICA, indicating that the geometry of the blood vessel plays significant roles in the atherogenesis. Specifically, the large curvature and planarity of the vessel were found to increase the risk for stenoses, because they tend to lower WSS and elevate OSI. Therefore, the geometric study makes it possible to estimate the stenosis location in the ICA siphon as long as the shape of ICA was measured.     

Errors in power-law estimations of inflow rates for intracranial aneurysm CFD

Patient-specific inflow rates are rarely available for computational fluid dynamics (CFD) studies of intracranial aneurysms. Instead, inflow rates are often estimated from parent artery diameters via power laws, i.e. Q ∝ Dⁿ, reflecting adaptation of conduit arteries to demanded flow. The present study aimed to validate the accuracy of these power laws. Internal carotid artery (ICA) flow rates were measured from 25 ICA aneurysm patients via 2D phase contrast MRI. ICA diameters, derived from 3D segmentation of rotational angiograms, were used to estimate inflow rates via power laws from the aneurysm CFD literature assuming the same inlet wall shear stress (WSS) (n = 3), velocity (n = 2) or flow rate (n = 0) for all cases. To illustrate the potential impact of errors in flow rate estimates, pulsatile CFD was carried out for four cases having large errors for at least one power law. Flow rates estimated by n = 3 and n = 0 power laws had significant (p < 0.01) mean biases of −22% to +32%, respectively, but with individual errors ranging from −78% to +120%. The n = 2 power law had no significant bias, but had non-negligible individual errors of −58% to +71%. CFD showed similarly large errors for time-averaged sac WSS; however, these were reduced after normalizing by parent artery WSS. High frequency WSS fluctuations, evident in 2/4 aneurysms, were also sensitive to inflow rate errors. Care should therefore be exercised in the interpretation of aneurysm CFD studies that rely on power law estimates of inflow rates, especially if absolute (vs. normalized) WSS, or WSS instabilities, are of interest.     

Can temporal fluctuation in spatial wall shear stress gradient initiate a cerebral aneurysm? A proposed novel hemodynamic index,the gradient oscillatory number (GON)

We propose a new hemodynamic index for the initiation of a cerebral aneurysm, defined by the temporal fluctuations of tension/compression forces acting on endothelial cells. We employed a patient-specific geometry of a human internal carotid artery (ICA) with an aneurysm, and reconstructed the geometry of the ICA before aneurysm formation by artificially removing the aneurysm. We calculated the proposed hemodynamic index and five other hemodynamic indices (wall shear stress (WSS) at peak systole, time-averaged WSS, time-averaged spatial WSS gradient, oscillatory shear index (OSI), and potential aneurysm formation indicator (AFI)) for the geometry before aneurysm formation using a computational fluid dynamics technique. By comparing the distribution of each index at the location of aneurysm formation, we discussed the validity of each. The results showed that only the proposed hemodynamic index had a significant correlation with the location of aneurysm formation. Our findings suggest that the proposed index may be useful as a hemodynamic index for the initiation of cerebral aneurysms.     

Automated and objective removal of bifurcation aneurysms: Incremental improvements,and validation against healthy controls

Abnormal hemodynamic stresses are thought to correlate with aneurysm initiation, growth, and rupture. We have previously investigated the role of wall shear stress (WSS) and WSS gradients (WSSG) in search for a mechanistic link to formation of sidewall aneurysms using an automated and objective tool for aneurysm removal and arterial reconstruction in combination with computational fluid dynamics (CFD). However, we warned against the use of the tool for bifurcation type aneurysms because of a potential unrealistic reconstruction of the apex. We hypothesized that inclusion of additional morphological features from the surrounding vasculature could overcome these constraints. We extended the previously published method for removal and reconstruction of the bifurcation vasculature based on diverging and converging points of the parent and daughter artery centerlines, to also include two new centerlines between the daughter vessels, one of them passed through the bifurcation center. Validation was performed by comparing the efficacy of the two algorithms, using ten healthy models of the internal carotid artery terminus as ground truth. Qualitative results showed that the bifurcation apexes became smoother relative to the original algorithm; more consistent with the reference models. This was reflected quantitatively by a reduced maximum distance between the reference and reconstructed surfaces, although not statistically significant. Furthermore, the modified algorithm also quantitatively improved CFD derived WSS and WSSG, especially the latter. In conclusion, the modified algorithm does not perfectly reconstruct the bifurcation apex, but provides an incremental improvement, especially important for the derived hemodynamic metrics of interest in vascular pathobiology.     

Computational fluid dynamics study of common stent models inside idealised curved coronary arteries

The haemodynamic behaviour of blood inside a coronary artery after stenting is greatly affected by individual stent features as well as complex geometrical properties of the artery including tortuosity and curvature. Regions at higher risk of restenosis, as measured by low wall shear stress (WSS < 0.5 Pa), have not yet been studied in detail in curved stented arteries. In this study, three-dimensional computational modelling and computational fluid dynamics methodologies were used to analyse the haemodynamic characteristics in curved stented arteries using several common stent models. Results in this study showed that stent strut thickness was one major factor influencing the distribution of WSS in curved arteries. Regions of low WSS were found behind struts, particularly those oriented at a large angle relative to the streamwise flow direction. These findings were similar to those obtained in studies of straight arteries. An uneven distribution of WSS at the inner and outer bends of curved arteries was observed where the WSS was lower at the inner bend. In this study, it was also shown that stents with a helical configuration generated an extra swirling component of the flow based on the helical direction; however, this extra swirl in the flow field did not cause significant changes on the distribution of WSS under the current setup.     

Numerical investigations of the haemodynamic changes associated with stent malapposition in an idealised coronary artery

The deployment of a coronary stent near complex lesions can sometimes lead to incomplete stent apposition (ISA), an undesirable side effect of coronary stent implantation. Three-dimensional computational fluid dynamics (CFD) calculations are performed on simplified stent models (with either square or circular cross-section struts) inside an idealised coronary artery to analyse the effect of different levels of ISA to the change in haemodynamics inside the artery. The clinical significance of ISA is reported using haemodynamic metrics like wall shear stress (WSS) and wall shear stress gradient (WSSG). A coronary stent with square cross-sectional strut shows different levels of reverse flow for malapposition distance (MD) between 0 mm and 0.12 mm. Chaotic blood flow is usually observed at late diastole and early systole for MD=0 mm and 0.12 mm but are suppressed for MD=0.06 mm. The struts with circular cross section delay the flow chaotic process as compared to square cross-sectional struts at the same MD and also reduce the level of fluctuations found in the flow field. However, further increase in MD can lead to chaotic flow not only at late diastole and early systole, but it also leads to chaotic flow at the end of systole. In all cases, WSS increases above the threshold value (0.5 Pa) as MD increases due to the diminishing reverse flow near the artery wall. Increasing MD also results in an elevated WSSG as flow becomes more chaotic, except for square struts at MD=0.06 mm.     

Approximating hemodynamics of cerebral aneurysms with steady flow simulations

Computational fluid dynamics (CFD) simulations can be employed to gain a better understanding of hemodynamics in cerebral aneurysms and improve diagnosis and treatment. However, introduction of CFD techniques into clinical practice would require faster simulation times. The aim of this study was to evaluate the use of computationally inexpensive steady flow simulations to approximate the aneurysm's wall shear stress (WSS) field. Two experiments were conducted. Experiment 1 compared for two cases the time-averaged (TA), peak systole (PS) and end diastole (ED) WSS field between steady and pulsatile flow simulations. The flow rate waveform imposed at the inlet was varied to account for variations in heart rate, pulsatility index, and TA flow rate. Consistently across all flow rate waveforms, steady flow simulations accurately approximated the TA, but not the PS and ED, WSS field. Following up on experiment 1, experiment 2 tested the result for the TA WSS field in a larger population of 20 cases covering a wide range of aneurysm volumes and shapes. Steady flow simulations approximated the space-averaged WSS with a mean error of 4.3%. WSS fields were locally compared by calculating the absolute error per node of the surface mesh. The coefficient of variation of the root-mean-square error over these nodes was on average 7.1%. In conclusion, steady flow simulations can accurately approximate the TA WSS field of an aneurysm. The fast computation time of 6 min per simulation (on 64 processors) could help facilitate the introduction of CFD into clinical practice.     

Mis-sizing of stent promotes intimal hyperplasia: impact of endothelial shear and intramural stress

Stent can cause flow disturbances on the endothelium and compliance mismatch and increased stress on the vessel wall. These effects can cause low wall shear stress (WSS), high wall shear stress gradient (WSSG), oscillatory shear index (OSI), and circumferential wall stress (CWS), which may promote neointimal hyperplasia (IH). The hypothesis is that stent-induced abnormal fluid and solid mechanics contribute to IH. To vary the range of WSS, WSSG, OSI, and CWS, we intentionally mismatched the size of stents to that of the vessel lumen. Stents were implanted in coronary arteries of 10 swine. Intravascular ultrasound (IVUS) was used to size the coronary arteries and stents. After 4 wk of stent implantation, IVUS was performed again to determine the extent of IH. In conjunction, computational models of actual stents, the artery, and non-Newtonian blood were created in a computer simulation to yield the distribution of WSS, WSSG, OSI, and CWS in the stented vessel wall. An inverse relation (R(2) = 0.59, P < 0.005) between WSS and IH was found based on a linear regression analysis. Linear relations between WSSG, OSI, and IH were observed (R(2) = 0.48 and 0.50, respectively, P < 0.005). A linear relation (R(2) = 0.58, P < 0.005) between CWS and IH was also found. More statistically significant linear relations between the ratio of CWS to WSS (CWS/WSS), the products CWS × WSSG and CWS × OSI, and IH were observed (R(2) = 0.67, 0.54, and 0.56, respectively, P < 0.005), suggesting that both fluid and solid mechanics influence the extent of IH. Stents create endothelial flow disturbances and intramural wall stress concentrations, which correlate with the extent of IH formation, and these effects were exaggerated with mismatch of stent/vessel size. These findings reveal the importance of reliable vessel and stent sizing to improve the mechanics on the vessel wall and minimize IH.     

MRI-based quantification of outflow boundary conditions for computational fluid dynamics of stenosed human carotid arteries

Accurate assessment of wall shear stress (WSS) is vital for studies on the pathogenesis of atherosclerosis. WSS distributions can be obtained by computational fluid dynamics (CFD) using patient-specific geometries and flow measurements. If patient-specific flow measurements are unavailable, in- and outflow have to be estimated, for instance by using Murray’s Law. It is currently unknown to what extent this law holds for carotid bifurcations, especially in cases where stenoses are involved. We performed flow measurements in the carotid bifurcation using phase-contrast MRI in patients with varying degrees of stenosis. An empirical relation between outflow and degree of area stenosis was determined and the outflow measurements were compared to estimations based on Murray’s Law. Furthermore, the influence of outflow conditions on the WSS distribution was studied.For bifurcations with an area stenosis smaller than 65%, the outflow ratio of the internal carotid artery (ICA) to the common carotid artery (CCA) was 0.62±0.12 while the outflow ratio of the external carotid artery (ECA) was 0.35±0.13. If the area stenosis was larger than 65%, the flow to the ICA decreased linearly to zero at 100% area stenosis. The empirical relation fitted the flow data well (R²=0.69), whereas Murray’s Law overestimated the flow to the ICA substantially for larger stenosis, resulting in an overestimation of the WSS. If patient-specific flow measurements of the carotid bifurcation are unavailable, estimation of the outflow ratio by the presented empirical relation will result in a good approximation of calculated WSS using CFD.     

Deviation from optimal vascular caliber control at middle cerebral artery bifurcations harboring aneurysms

Cerebral aneurysms form preferentially at arterial bifurcations. The vascular optimality principle (VOP) decrees that minimal energy loss across bifurcations requires optimal caliber control between radii of parent (r₀) and daughter branches (r₁ and r₂): r₀ⁿ=r₁ⁿ+r₂ⁿ, with n approximating three. VOP entails constant wall shear stress (WSS), an endothelial phenotype regulator. We sought to determine if caliber control is maintained in aneurysmal intracranial bifurcations. Three-dimensional rotational angiographic volumes of 159 middle cerebral artery (MCA) bifurcations (62 aneurysmal) were processed using 3D gradient edge-detection filtering, enabling threshold-insensitive radius measurement. Radius ratio (RR)=r₀³/(r₁³+r₂³) and estimated junction exponent (n) were compared between aneurysmal and non-aneurysmal bifurcations using Student t-test and Wilcoxon rank-sum analysis. The results show that non-aneurysmal bifurcations display optimal caliber control with mean RR of 1.05 and median n of 2.84. In contrast, aneurysmal bifurcations had significantly lower RR (0.76, p<.0001) and higher n (4.28, p<.0001). Unexpectedly, 37% of aneurysmal bifurcations revealed a daughter branch larger than its parent vessel, an absolute violation of optimality, not witnessed in non-aneurysmal bifurcations. The aneurysms originated more often off the smaller daughter (52%) vs. larger daughter branch (16%). Aneurysm size was not statistically correlated to RR or n. Aneurysmal males showed higher deviation from VOP. Non-aneurysmal MCA bifurcations contralateral to aneurysmal ones showed optimal caliber control. Aneurysmal bifurcations, in contrast to non-aneurysmal counterparts, disobey the VOP and may exhibit dysregulation in WSS-mediated caliber control. The mechanism of this focal divergence from optimality may underlie aneurysm pathogenesis and requires further study.     

Comparison of the hemodynamics in 6mm and 4-7 mm hemodialysis grafts by means of CFD

The aim of our study is to investigate with computational fluid dynamics (CFD) whether different arterial anastomotic geometries result in a different hemodynamics at the arterial (AA) and venous anastomosis (VA) of hemodialysis vascular access grafts. We have studied a 6mm graft (CD) and a 4-7 mm graft (TG). A validated three-dimensional CFD model is developed to simulate flow in the two graft types. Only the arterial anastomosis (AA) geometry differs. The boundary conditions applied are a periodic velocity signal at the arterial inlet and a periodic pressure wave at the venous outlet. Flow rate is set to 1,000 ml/min. The time dependent Navier-Stokes equations are solved. Wall shear stress (WSS), wall shear stress gradient (WSSG) and pressure gradient (PG) are calculated. Anastomotic flow is asymmetric although the anastomosis geometry is symmetric. The hemodynamic parameters, WSS, WSSG and PG, values at the suture line of the arterial anastomosis of the TG are at least twice as much as in the CD. Comparing the parameters at the two AA indicate that little flow rate increase introduces the risk of hemolysis in the TG whereas the CD is completely free of hemolysis. The hemodynamic parameter values at the venous anastomosis of the CD are 24 till 35% higher compared to the values of the TG. WSS values (> 3 Pa) in the VA are in the critical range for stenosis development in both graft geometries. The zones where the parameters reach extreme values correspond to the locations where intimal hyperplasia formation is reported in literature. In all anastomoses, the hemodynamic parameter levels are in the range where leucocytes and platelets get activated. Our simulations confirm clinical results where TG did not show a better outcome when compared to the CD.     

Influence of stent-induced vessel deformation on hemodynamic feature of bloodstream inside ICA aneurysms

One of the effective treatment options for intracranial aneurysms is stent-assisted coiling. Though, previous works have demonstrated that stent usage would result in the deformation of the local vasculature. The effect of simple stent on the blood hemodynamics is still uncertain. In this work, hemodynamic features of the blood stream on four different ICA aneurysm with/without interventional are investigated. To estimate the relative impacts of vessel deformation, four distinctive ICA aneurysm is simulated by the one-way FSI technique. Four hemodynamic factors of aneurysm blood velocity, wall pressure and WSS are compared in the peak systolic stage to disclose the impact of defamation by the stent in two conditions. The stent usage would decrease almost all of the mentioned parameters, except for OSI. Stenting reduces neck inflow rate, while the effect of interventional was not consistent among the aneurysms. The deformation of an aneurysm has a strong influence on the hemodynamics of an aneurysm. This outcome is ignored by most of the preceding investigations, which focused on the pre-interventional state for studying the relationship between hemodynamics and stents. Present results show that the application of stent without coiling would improve most hemodynamic factors, especially when the deformation of the aneurysm is high enough.

     

An inelastic multi-mechanism constitutive equation for cerebral arterial tissue

Intracranial aneurysms (ICA) are abnormal saccular dilations of cerebral arteries, commonly found at apices of arterial bifurcations and outer walls of curved arterial segments. Histological evidence suggests the stages in ICA development include the deformation of a segment of arterial wall into a “bleb” with no identifiable neck region followed by the development of an aneurysm with a clear neck. Afterwards, the aneurysm may undergo further enlargement, possibly with significant biological response including calcification and thrombosis. Past studies of the biomechanics of cerebral aneurysm tissue have been directed at modeling elastic deformations of pre-existing aneurysms. Taking this approach, the aneurysm wall is treated as a different entity than the arterial tissue from which it developed. In the current work, a nonlinear, inelastic, dual-mechanism constitutive equation for cerebral arterial tissue is developed. It is the first to model the recruitment of collagen fibers and degradation of the internal elastic lamina, two important characteristics of early stage aneurysm formation.     

A theoretical model for fibroblast-controlled growth of saccular cerebral aneurysms

A new theoretical model for the growth of saccular cerebral aneurysms is proposed by extending the recent constitutive framework of Kroon and Holzapfel [2007a. A model for saccular cerebral aneurysm growth by collagen fibre remodelling. J. Theor. Biol. 247, 775-787]. The continuous turnover of collagen is taken to be the driving mechanism in aneurysmal growth. The collagen production rate depends on the magnitude of the cyclic deformation of fibroblasts, caused by the pulsating blood pressure during the cardiac cycle. The volume density of fibroblasts in the aneurysmal tissue is taken to be constant throughout the growth process. The growth model is assessed by considering the inflation of an axisymmetric membranous piece of aneurysmal tissue, with material characteristics representative of a cerebral aneurysm. The diastolic and systolic states of the aneurysm are computed, together with its load-free state. It turns out that the value of collagen pre-stretch, that determines growth speed and stability of the aneurysm, is of pivotal importance. The model is able to predict aneurysms with typical berry-like shapes observed clinically, and the predicted wall stresses correlate well with the experimentally obtained ultimate stresses of this type of tissue. The model predicts that aneurysms should fail when reaching a size of about 1.2-3.6 mm, which is smaller than what has been clinically observed. With some refinements, the model may, however, be used to predict future growth of diagnosed aneurysms.     

Effects of vessel compliance on flow pattern in porcine epicardial right coronary arterial tree

The compliance of the vessel wall affects hemodynamic parameters which may alter the permeability of the vessel wall. Based on experimental measurements, the present study established a finite element (FE) model in the proximal elastic vessel segments of epicardial right coronary arterial (RCA) tree obtained from computed tomography. The motion of elastic vessel wall was measured by an impedance catheter and the inlet boundary condition was measured by an ultrasound flow probe. The Galerkin FE method was used to solve the Navier–Stokes and Continuity equations, where the convective term in the Navier–Stokes equation was changed in the arbitrary Lagrangian–Eulerian (ALE) framework to incorporate the motion due to vessel compliance. Various hemodynamic parameters (e.g., wall shear stress—WSS, WSS spatial gradient—WSSG, oscillatory shear index—OSI) were analyzed in the model. The motion due to vessel compliance affects the time-averaged WSSG more strongly than WSS at bifurcations. The decrease of WSSG at flow divider in elastic bifurcations, as compared to rigid bifurcations, implies that the vessel compliance decreases the permeability of vessel wall and may be atheroprotective. The model can be used to predict coronary flow pattern in subject-specific anatomy as determined by noninvasive imaging.     

Direct numerical simulation of transitional flow in a patient-specific intracranial aneurysm

In experiments turbulence has previously been shown to occur in intracranial aneurysms. The effects of turbulence induced oscillatory wall stresses could be of great importance in understanding aneurysm rupture. To investigate the effects of turbulence on blood flow in an intracranial aneurysm, we performed a high resolution computational fluid dynamics (CFD) simulation in a patient specific middle cerebral artery (MCA) aneurysm using a realistic, pulsatile inflow velocity. The flow showed transition to turbulence just after peak systole, before relaminarization occurred during diastole. The turbulent structures greatly affected both the frequency of change of wall shear stress (WSS) direction and WSS magnitude, which reached a maximum value of 41.5Pa. The recorded frequencies were predominantly in the range of 1-500Hz. The current study confirms, through properly resolved CFD simulations that turbulence can occur in intracranial aneurysms.     

Flow instability detected in ruptured versus unruptured cerebral aneurysms at the internal carotid artery

Flow instability has emerged as a new hemodynamic metric hypothesized to have potential value in assessing the rupture risk of cerebral aneurysms. However, diverse findings have been reported in the literature. In the present study, high-resolution hemodynamic simulations were performed retrospectively on 35 aneurysms (10 ruptured & 25 unruptured) located at the internal carotid artery (ICA). Simulated hemodynamic parameters were statistically compared between the ruptured and unruptured aneurysms, with emphasis on examining the correlation of flow instability with the status of aneurysm rupture. Pronounced flow instability was detected in 20% (2 out of 10) of the ruptured aneurysms, whereas in 44% (11 out of 25) of the unruptured aneurysms. Statistically, the flow instability metric (quantified by the temporally and spatially averaged fluctuating kinetic energy over the aneurysm sac) did not differ significantly between the ruptured and unruptured aneurysms. In contrast, low wall shear stress area (LSA) and pressure loss coefficient (PLC) exhibited significant correlations with the status of aneurysm rupture. In conclusion, the present study suggests that the presence of flow instability may not correlate closely with the status of aneurysm rupture, at least for ICA aneurysms. On the other hand, the retrospective nature of the study and the small sample size may have to some extent compromised the reliability of the conclusion, and therefore large-scale prospective studies would be needed to further address the issue.     

Tetrahedral vs. polyhedral mesh size evaluation on flow velocity and wall shear stress for cerebral hemodynamic simulation

Haemodynamic factors, in particular wall shear stresses (WSSs) may have significant impact on growth and rupture of cerebral aneurysms. Without a means to measure WSS reliably in vivo, computational fluid dynamic (CFD) simulations are frequently employed to visualise and quantify blood flow from patient-specific computational models. With increasing interest in integrating these CFD simulations into pretreatment planning, a better understanding of the validity of the calculations in respect to computation parameters such as volume element type, mesh size and mesh composition is needed. In this study, CFD results for the two most common aneurysm types (saccular and terminal) are compared for polyhedral- vs. tetrahedral-based meshes and discussed regarding future clinical applications. For this purpose, a set of models were constructed for each aneurysm with spatially varying surface and volume mesh configurations (mesh size range: 5119-258, 481 volume elements). WSS distribution on the model wall and point-based velocity measurements were compared for each configuration model. Our results indicate a benefit of polyhedral meshes in respect to convergence speed and more homogeneous WSS patterns. Computational variations of WSS values and blood velocities are between 0.84 and 6.3% from the most simple mesh (tetrahedral elements only) and the most advanced mesh design investigated (polyhedral mesh with boundary layer).