Background and Purpose-Evolution of intracranial aneurysmal disease is known to be related to hemodynamic forces acting on the vessel wall. Low wall shear stress (WSS) has been reported to have a negative effect on endothelial cells normal physiology and may be an important contributor to local remodeling of the arterial wall and to aneurysm growth and rupture. Methods-Seven patient-specific models of intracranial aneurysms were constructed using MR angiography data acquired at two different time points (mean 16.4Ϯ7.4 months between the two time points). Numeric simulations of the flow in the baseline geometries were performed to compute WSS distributions. The lumenal geometries constructed from the two time points were manually coregistered, and the radial displacement of the wall was calculated on a pixel-by-pixel basis. This displacement, corresponding to the local growth of the aneurysm, was compared to the time-averaged wall shear stress (WSS TA ) through the cardiac cycle at that location. For statistical analysis, radial displacement was considered to be significant if it was larger than half of the MR pixel resolution (0.3 mm). Results-Mean

Evolution of intracranial aneurysms is known to be related to hemodynamic forces such as wall shear stress (WSS) and maximum shear stress (MSS). Estimation of these parameters can be performed using numerical simulations with computational fluid dynamics (CFD), but can also be directly measured with magnetic resonance imaging (

Abstract-Thrombus formation in intracranial aneurysms, while sometimes stabilizing lesion growth, can present additional risk of thrombo-embolism. The role of hemodynamics in the progression of aneurysmal disease can be elucidated by patient-specific computational modeling. In our previous work, patient-specific computational fluid dynamics (CFD) models were constructed from MRI data for three patients who had fusiform basilar aneurysms that were thrombus-free and then proceeded to develop intraluminal thrombus. In this study, we investigated the effect of increased flow residence time (RT) by modeling passive scalar advection in the same aneurysmal geometries. NonNewtonian pulsatile flow simulations were carried out in base-line geometries and a new postprocessing technique, referred to as ''virtual ink'' and based on the passive scalar distribution maps, was used to visualize the flow and estimate the flow RT. The virtual ink technique clearly depicted regions of flow separation. The flow RT at different locations adjacent to aneurysmal walls was calculated as the time the virtual ink scalar remained above a threshold value. The RT values obtained in different areas were then correlated with the location of intra-aneurysmal thrombus observed at a follow-up MR study. For each patient, the wall shear stress (WSS) distribution was also obtained from CFD simulations and correlated with thrombus location. The correlation analysis determined a significant relationship between regions where CFD predicted either an increased RT or low WSS and the regions where thrombus deposition was observed to occur in vivo. A model including both low WSS and increased RT predicted thrombus-prone regions significantly better than the models with RT or WSS alone.

The deposition of intralumenal thrombus in intracranial aneurysms adds a risk of thrombo-embolism over and above that posed by mass-effect and rupture. In addition to biochemical factors, hemodynamic factors that are governed by lumenal geometry and blood flow rates likely play an important role in the thrombus formation and deposition process. In this study, patient-specific computational fluid dynamics (CFD) models of blood flow were constructed from MRA data for three patients who had fusiform basilar aneurysms that were thrombus-free and then proceeded to develop intra-lumenal thrombus. In order to determine whether features of the flow fields could suggest which regions had an elevated potential for thrombus deposition, the flow was modeled in the baseline, thrombus-free geometries. Pulsatile flow simulations were carried out using patientspecific inlet flow conditions measured with MR velocimetry. Newtonian and non-Newtonian blood behavior was considered. A strong similarity was found between the intra-aneurysmal regions with CFD-predicted slow, recirculating flows and the regions of thrombus deposition observed in vivo in the follow-up MR studies. In two cases with larger aneurysms, the agreement between the low velocity zones and clotted off regions improved when non-Newtonian blood behavior was taken into account. A similarity was also found between the calculated low shear stress regions and the regions that were later observed to clot.