Vitaliy L. Rayz

Aneurysm Growth Occurs at Region of Low Wall Shear Stress Patient-Specific Correlation of Hemodynamics and Growth in a Longitudinal Study

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 (WSSTA) 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 WSSTA values obtained for the areas with a displacement smaller and greater than 0.3 mm were 2.55±3.6 and 0.76±1.5 Pa, respectively (P<0.001). A linear correlation analysis demonstrated a significant relationship between WSSTA and surface displacement (P<0.001). Conclusions— These results indicate that aneurysm growth is likely to occur in regions where the endothelial layer lining the vessel wall is exposed to abnormally low wall shear stress.

Phase-contrast magnetic resonance imaging measurements in intracranial aneurysms in vivo of flow patterns, velocity fields, and wall shear stress: comparison with computational fluid dynamics.

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 (MRI) using a time‐dependent 3D phase‐contrast sequence with encoding of each of the three components of the velocity vectors (7D‐MRV). To study the accuracy of 7D‐MRV in estimating these parameters in vivo, in comparison with CFD, 7D‐MRV and patient‐specific CFD modeling was performed for 3 patients who had intracranial aneurysms. Visual and quantitative analyses of the flow pattern and distribution of velocities, MSS, and WSS were performed using the two techniques. Spearmans coefficients of correlation between the two techniques were 0.56 for the velocity field, 0.48 for MSS, and 0.59 for WSS. Visual analysis and Bland–Altman plots showed good agreement for flow pattern and velocities but large discrepancies for MSS and WSS. These results indicate that 7D‐MRV can be used in vivo to measure velocity flow fields and for estimating MSS and WSS. Currently, however, this method cannot accurately quantify the latter two parameters. Magn Reson Med 61:409–417, 2009.

Giant Intracranial Aneurysms: Evolution of Management in a Contemporary Surgical Series

BACKGROUND Many significant microsurgical series of patients with giant aneurysms predate changes in practice during the endovascular era. OBJECTIVE A contemporary surgical experience is presented to examine changes in management relative to earlier reports, to establish the role of open microsurgery in the management strategy, and to quantify results for comparison with evolving endovascular therapies. METHODS During a 13-year period, 140 patients with 141 giant aneurysms were treated surgically. One hundred aneurysms (71%) were located in the anterior circulation, and 41 aneurysms were located in the posterior circulation. RESULTS One hundred eight aneurysms (77%) were completely occluded, 14 aneurysms (10%) had minimal residual aneurysm, and 16 aneurysms (11%) were incompletely occluded with reversed or diminished flow. Three patients with calcified aneurysms were coiled after unsuccessful clipping attempts. Eighteen patients died in the perioperative period (surgical mortality, 13%). Bypass-related complications resulted from bypass occlusion (7 patients), aneurysm hemorrhage due to incomplete aneurysm occlusion (4 patients), or aneurysm thrombosis with perforator or branch artery occlusion (4 patients). Thirteen patients were worse at late follow-up (permanent neurological morbidity, 9%; mean length of follow-up, 23 ± 1.9 months). Overall, good outcomes (Glasgow Outcome Score 5 or 4) were observed in 114 patients (81%), and 109 patients (78%) were improved or unchanged after therapy. CONCLUSION A heavy reliance on bypass techniques plus indirect giant aneurysm occlusion distinguishes this contemporary surgical experience from earlier ones, and obviates the need for hypothermic circulatory arrest. Experienced neurosurgeons can achieve excellent results with surgery as the “first-line” management approach and endovascular techniques as adjuncts to surgery.

Flow Residence Time and Regions of Intraluminal Thrombus Deposition in Intracranial Aneurysms

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. Non-Newtonian 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.

Numerical Simulations of Flow in Cerebral Aneurysms: Comparison of CFD Results and In Vivo MRI Measurements

Computational fluid dynamics (CFD) methods can be used to compute the velocity field in patient-specific vascular geometries for pulsatile physiological flow. Those simulations require geometric and hemodynamic boundary values. The purpose of this study is to demonstrate that CFD models constructed from patient-specific magnetic resonance (MR) angiography and velocimetry data predict flow fields that are in good agreement with in vivo measurements and therefore can provide valuable information for clinicians. The effect of the inlet flow rate conditions on calculated velocity fields was investigated. We assessed the internal consistency of our approach by comparing CFD predictions of the in-plane velocity field to the corresponding in vivo MR velocimetry measurements. Patient-specific surface models of four basilar artery aneurysms were constructed from contrast-enhanced MR angiography data. CFD simulations were carried out in those models using patient-specific flow conditions extracted from MR velocity measurements of flow in the inlet vessels. The simulation results computed for slices through the vasculature of interest were compared with in-plane velocity measurements acquired with phase-contrast MR imaging in vivo. The sensitivity of the flow fields to inlet flow ratio variations was assessed by simulating five different inlet flow scenarios for each of the basilar aneurysm models. In the majority of cases, altering the inlet flow ratio caused major changes in the flow fields predicted in the aneurysm. A good agreement was found between the flow fields measured in vivo using the in-plane MR velocimetry technique and those predicted with CFD simulations. The study serves to demonstrate the consistency and reliability of both MR imaging and numerical modeling methods. The results demonstrate the clinical relevance of computational models and suggest that realistic patient-specific flow conditions are required for numerical simulations of the flow in aneurysmal blood vessels.

Carotid Atheroma Rupture Observed In Vivo and FSI-Predicted Stress Distribution Based on Pre-rupture Imaging

Atherosclerosis at the carotid bifurcation is a major risk factor for stroke. As mechanical forces may impact lesion stability, finite element studies have been conducted on models of diseased vessels to elucidate the effects of lesion characteristics on the stresses within plaque materials. It is hoped that patient-specific biomechanical analyses may serve clinically to assess the rupture potential for any particular lesion, allowing better stratification of patients into the most appropriate treatments. Due to a sparsity of in vivo plaque rupture data, the relationship between various mechanical descriptors such as stresses or strains and rupture vulnerability is incompletely known, and the patient-specific utility of biomechanical analyses is unclear. In this article, we present a comparison between carotid atheroma rupture observed in vivo and the plaque stress distribution from fluid–structure interaction analysis based on pre-rupture medical imaging. The effects of image resolution are explored and the calculated stress fields are shown to vary by as much as 50% with sub-pixel geometric uncertainty. Within these bounds, we find a region of pronounced elevation in stress within the fibrous plaque layer of the lesion with a location and extent corresponding to that of the observed site of plaque rupture.

Bypass Surgery for the Treatment of Dolichoectatic Basilar Trunk Aneurysms: A Work in Progress.

BACKGROUND The treatment of dolichoectatic basilar trunk aneurysms has been ineffectual or morbid due to nonsaccular morphology, deep location, and involvement of brainstem perforators. Treatment with bypass surgery has been advocated to eliminate malignant hemodynamics and to stabilize aneurysm growth. OBJECTIVE To validate that flow alteration with bypass and parent artery occlusion favorably impacts aneurysm progression. METHODS Surgical management evolved in 3 phases, each with different hemodynamic alterations. RESULTS During a 17-year period, 37 patients with dolichoectatic basilar trunk aneurysms were retrospectively identified, of whom 21 patients were observed, 12 treated immediately, and 4 selected for treatment after clinical progression. In phase 1, flow reversal was overly thrombogenic, despite heparin (N = 5, final mortality, 100%). In phase 2, flow reduction with intracranial-to-intracranial bypass was safer than flow reversal, but did not prevent progressive aneurysm enlargement (N = 3, final mortality 67%). In phase 3, distal clip occlusion of the basilar trunk aneurysm preserved anterograde flow in the aneurysm without rupture, but reduced flow threatened perforator patency, despite treatment with clopidogrel (N = 8, final mortality 62%). CONCLUSION Shifting treatment strategy for dolichoectatic basilar trunk aneurysms improved surgical (80% to 50%) and final mortalities (100% to 62%), with stabilization of aneurysms in the phase 3 survivors. Good outcomes are determined by perforator preservation and mitigating aneurysm thrombosis. Occlusion techniques with increased distal run-off seem to benefit perforators. The treatment of dolichoectatic basilar trunk aneurysms can advance through concentrated management in dedicated centers, concerted efforts to study morphology and hemodynamics with computational methods, and widespread collection of registry data. ABBREVIATIONS 4D PC-MRI, time-resolved phase-contrast MRIAICA, anterior inferior cerebellar arteryCE-MRA, high-resolution contrast-enhanced MR angiographyEC-IC, extracranial-to-intracranial bypassMCA, middle cerebral arteryMR, magnetic resonancemRS, modified Rankin ScalePCA, posterior cerebral arteryPICA, posterior inferior cerebellar arterySCA, superior cerebellar arterySTA, superficial temporal arteryVA, vertebral artery.

An efficient two-stage approach for image-based FSI analysis of atherosclerotic arteries

Patient-specific biomechanical modeling of atherosclerotic arteries has the potential to aid clinicians in characterizing lesions and determining optimal treatment plans. To attain high levels of accuracy, recent models use medical imaging data to determine plaque component boundaries in three dimensions, and fluid–structure interaction is used to capture mechanical loading of the diseased vessel. As the plaque components and vessel wall are often highly complex in shape, constructing a suitable structured computational mesh is very challenging and can require a great deal of time. Models based on unstructured computational meshes require relatively less time to construct and are capable of accurately representing plaque components in three dimensions. These models unfortunately require additional computational resources and computing time for accurate and meaningful results. A two-stage modeling strategy based on unstructured computational meshes is proposed to achieve a reasonable balance between meshing difficulty and computational resource and time demand. In this method, a coarsegrained simulation of the full arterial domain is used to guide and constrain a fine-scale simulation of a smaller region of interest within the full domain. Results for a patient-specific carotid bifurcation model demonstrate that the two-stage approach can afford a large savings in both time for mesh generation and time and resources needed for computation. The effects of solid and fluid domain truncation were explored, and were shown to minimally affect accuracy of the stress fields predicted with the two-stage approach.

Numerical Simulation of Pre- and Postsurgical Flow in a Giant Basilar Aneurysm

Computational modeling of the flow in cerebral aneurysms is an evolving technique that may play an important role in surgical planning. In this study, we simulated the flow in a giant basilar aneurysm before and after surgical takedown of one vertebral artery. Patient-specific geometry and flowrates obtained from magnetic resonance (MR) angiography and velocimetry were used to simulate the flow prior to and after the surgery. Numerical solutions for steady and pulsatile flows were obtained. Highly three-dimensional flows, with strong secondary flows, were computed in the aneurysm in the presurgical and postsurgical conditions. The computational results predicted that occlusion of a vertebral artery would result in a significant increase of the slow flow region formed in the bulge of the aneurysm, where increased particle residence time and velocities lower than 2.5 cms were computed. The region of slow flow was found to have filled with thrombus following surgery. Predictions of numerical simulation methods are consistent with the observed outcome following surgical treatment of an aneurysm. The study demonstrates that computational models may provide hypotheses to test in future studies, and might offer guidance for the interventional treatment of cerebral aneurysms.

Computational Modeling of Vascular Hemodynamics

In addition to biochemical factors, hemodynamic factors that are governed by lumenal geometry and blood flow rates likely play an important role in the pathogenesis of cardiovascular disease. Numerous computational and experimental studies indicated correlation of certain hemodynamic parameters with initiation and progression of atherosclerotic plaques, while flow variables possibly affecting aneurysmal disease are still disputed. This chapter presents a review of publications on current state of the art for the Computational Fluid Dynamics (CFD) methods used in patient-specific flow modeling. Typical modeling assumptions and boundary conditions are described and discussed as well as some of the post-processing and visualization techniques. It is hoped that with the current advances in medical imaging and numerical methods, computational modeling will evolve into a clinical tool providing guidance for cardiovascular disease treatment on a patient-by-patient basis.