Biomechanical behaviour of cerebral aneurysm and its relation with the formation of intraluminal thrombus: a patient-specific modelling study

Influence of vortical structures on fibrin clot formation in cerebral aneurysms: A two-dimensional computational study

Thrombosis is an important contributor to cerebral aneurysm growth and progression. A number of sophisticated multiscale and multiphase in silico models have been developed with a view towards interventional planning. Many of these models are able to account for clotting outcomes, but do not provide detailed insight into the role of flow during clot development. In this study, we present idealised, two-dimensional in silico cerebral fibrin clot model based on computational fluid dynamics (CFD), biochemical modelling and variable porosity, permeability, and diffusivity. The model captures fibrin clot growth in cerebral aneurysms over a period at least 1000 s in five different geometries. The fibrin clot growth results were compared to an experiment presented in literature. The biochemistry was found to be more sensitive to mesh size compared to the haemodynamics, while larger timesteps overpredicted clot size in pulsatile flow. When variable diffusivity was used, the predicted clot size was 25.4% lesser than that with constant diffusivity. The predicted clot size in pulsatile flow was 14.6% greater than in plug flow. Different vortex modes were observed in plug and pulsatile flow; the latter presented smaller intermediate modes where the main vortex was smaller and less likely to disrupt the growing fibrin clot. Furthermore, smaller vortex modes were seen to support fibrin clot propagation across geometries. The model clearly demonstrates how the growing fibrin clot alters vortical structures within the aneurysm sac and how this changing flow, in turn, shapes the growing fibrin clot.

Numerical investigation of unruptured middle cerebral artery bifurcation aneurysms: influence of aspect ratio

An aneurysm is a disease condition, which is due to the pathological weakening of an arterial wall. These aneurysms are often found in various branch points and bifurcations of an artery in the cerebral circulation. Most aneurysms come to medical attention, either due to brain hemorrhages caused by rupture or found unruptured. To consider surgically invasive treatment modalities, clinicians need scientific methods such as, hemodynamic analysis to assess rupture risk. The arterial wall loses its structural integrity when wall shear stress (WSS) and other hemodynamic parameters exceed a certain threshold. In the present study, numerical simulations are carried out for unruptured middle cerebral artery (MCA) aneurysms. Three distinct representative sizes are chosen from a larger patient pool of 26 MCA aneurysms. Logically, these aneurysms represent three growth stages of any patient with similar anatomical structure. Simulations are performed to compare the three growth phases (with different aspect ratios) of an aneurysm and correlate their hemodynamic parameters. Simulations with patient specific boundary conditions reveal that, aneurysms with a higher aspect ratio (AR) correspond to an attendant decrease in both time-averaged wall shear stress (TAWSS) and spatial wall shear stress gradients (WSSG). Smaller MCAs were observed to have higher positive wall shear stress divergence (WSSD), exemplifying the tensile nature of arterial wall stretching. Present study identifies positive wall shear stress divergence (PWSSD) to be a potential biomarker for evaluating the growth of an aneurysm.

Structural modelling of the cardiovascular system

Computational modelling of the cardiovascular system offers much promise, but represents a truly interdisciplinary challenge, requiring knowledge of physiology, mechanics of materials, fluid dynamics and biochemistry. This paper aims to provide a summary of the recent advances in cardiovascular structural modelling, including the numerical methods, main constitutive models and modelling procedures developed to represent cardiovascular structures and pathologies across a broad range of length and timescales; serving as an accessible point of reference to newcomers to the field. The class of so-called hyperelastic materials provides the theoretical foundation for the modelling of how these materials deform under load, and so an overview of these models is provided; comparing classical to application-specific phenomenological models. The physiology is split into components and pathologies of the cardiovascular system and linked back to constitutive modelling developments, identifying current state of the art in modelling procedures from both clinical and engineering sources. Models which have originally been derived for one application and scale are shown to be used for an increasing range and for similar applications. The trend for such approaches is discussed in the context of increasing availability of high performance computing resources, where in some cases computer hardware can impact the choice of modelling approach used.

In vitro measurement of platelet adhesion to intact endothelial cells under low shear conditions

BACKGROUND

Prediction of thrombus formation at intact arterial walls under low shear flow conditions is clinically important particularly for better prognoses of embolisation in cerebral aneurysms. Although a new mathematical model for this purpose is necessary, little quantitative information has been known about platelet adhesion to intact endothelial cells.

OBJECTIVE

The objective of this study is to measure the number of platelets adhering to intact endothelial cells with a focus upon the influence of the shear rate.

METHODS

Endothelial cells disseminated in μ-slides were exposed to swine whole blood at different shear rates. Adenosine diphosphate (ADP) was used as an agonist. Adherent platelets were counted by means of scanning electron microscopy.

RESULTS

At an ADP concentration of 1 µM, 20.8 ± 3.1 platelets per 900 µm2 were observed after 30-minute perfusion at a shear rate of 0.8 s-1 whereas only 3.0 ± 1.4 per 900 µm2 at 16.8 s-1.

CONCLUSIONS

The number of adherent platelets is determined by a balance between the shear and the degree of stimulation by the agonist. At an ADP concentration of 1 µM, a limit to the shear rate at which platelets can adhere to intact endothelial cells is considered to be slightly higher than 16.8 s-1.

Phase-contrast MRI versus numerical simulation to quantify hemodynamical changes in cerebral aneurysms after flow diverter treatment

Cerebral aneurysms are a major risk factor for intracranial bleeding with devastating consequences for the patient. One recently established treatment is the implantation of flow-diverters (FD). Methods to predict their treatment success before or directly after implantation are not well investigated yet. The aim of this work was to quantitatively study hemodynamic parameters in patient-specific models of treated cerebral aneurysms and its correlation with the clinical outcome. Hemodynamics were evaluated using both computational fluid dynamics (CFD) and phase contrast (PC) MRI. CFD simulations and in vitro MRI measurements were done under similar flow conditions and results of both methods were comparatively analyzed. For preoperative and postoperative distribution of hemodynamic parameters, CFD simulations and PC-MRI velocity measurements showed similar results. In both cases where no occlusion of the aneurysm was observed after six months, a flow reduction of about 30-50% was found, while in the clinically successful case with complete occlusion of the aneurysm after 6 months, the flow reduction was about 80%. No vortex was observed in any of the three models after treatment. The results are in agreement with recent studies suggesting that CFD simulations can predict post-treatment aneurysm flow alteration already before implantation of a FD and PC-MRI could validate the predicted hemodynamic changes right after implantation of a FD.

Experimental and CFD flow studies in an intracranial aneurysm model with Newtonian and non-Newtonian fluids

BACKGROUND

According to the clinical data, flow conditions play a major role in the genesis of intracranial aneurysms. The disorder of the flow structure is the cause of damage of the inner layer of the vessel wall, which leads to the development of cerebral aneurysms. Knowledge of the alteration of the flow field in the aneurysm region is important for treatment.

OBJECTIVE

The aim is to study quantitatively the flow structure in an patient-specific aneurysm model of the internal carotid artery using both experimental and computational fluid dynamics (CFD) methods with Newtonian and non-Newtonian fluids.

METHODS

A patient-specific geometry of aneurysm of the internal carotid artery was used. Patient data was segmented and smoothed to obtain geometrical model. An elastic true-to-scale silicone model was created with stereolithography. For initial investigation of the blood flow, the flow was visualized by adding particles into the silicone model. The precise flow velocity measurements were done using 1D Laser Doppler Anemometer with a spatial resolution of 50 μ m and a temporal resolution of 1 ms. The local velocity measurements were done at a distance of 4 mm to each other. A fluid with non-Newtonian properties was used in the experiment. The CFD simulations for unsteady-state problem were done using constructed hexahedral mesh for Newtonian and non-Newtonian fluids.

RESULTS

Using 1D laser Doppler Anemometer the minimum velocity magnitude at the end of systole -0.01 m/s was obtained in the aneurysm dome while the maximum velocity 1 m/s was at the center of the outlet segment. On central cross section of the aneurysm the maximum velocity value is only 20% of the average inlet velocity. The average velocity on the cross-section is only 11% of the inlet axial velocity. Using the CFD simulation the wall shear stresses for Newtonian and non-Newtonian fluid at the end of systolic phase (t= 0.25 s) were computed. The wall shear stress varies from 3.52 mPa (minimum value) to 10.21 Pa (maximum value) for the Newtonian fluid. For the non-Newtonian fluid the wall shear stress minimum is 2.94 mPa; the maximum is 9.14 Pa. The lowest value of the wall shear stress for both fluids was obtained at the dome of the aneurysm while the highest wall shear stress was at the beginning of the outlet segment. The vortex in the aneurysm region is unstable during the cardiac cycle. The clockwise rotation of the streamlines at the inlet segment for Newtonian and non-Newtonian fluid is shown.

CONCLUSION

The results of the present study are in agreement with the hemodynamics theory of aneurysm genesis. Low value of wall shear stress is observed at the aneurysm dome which can cause a rupture of an aneurysm.

Understanding the role of hemodynamics in the initiation, progression, rupture, and treatment outcome of cerebral aneurysm from medical image-based computational studies

About a decade ago, the first image-based computational hemodynamic studies of cerebral aneurysms were presented. Their potential for clinical applications was the result of a right combination of medical image processing, vascular reconstruction, and grid generation techniques used to reconstruct personalized domains for computational fluid and solid dynamics solvers and data analysis and visualization techniques. A considerable number of studies have captivated the attention of clinicians, neurosurgeons, and neuroradiologists, who realized the ability of those tools to help in understanding the role played by hemodynamics in the natural history and management of intracranial aneurysms. This paper intends to summarize the most relevant results in the field reported during the last years.