Newtonian and non-Newtonian blood flow in coiled cerebral aneurysms

Roadmap for Clinical Translation of Mobile Microrobotics

Medical microrobotics is an emerging field to revolutionize clinical applications in diagnostics and therapeutics of various diseases. On the other hand, the mobile microrobotics field has important obstacles to pass before clinical translation. This article focuses on these challenges and provides a roadmap of medical microrobots to enable their clinical use. From the concept of a "magic bullet" to the physicochemical interactions of microrobots in complex biological environments in medical applications, there are several translational steps to consider. Clinical translation of mobile microrobots is only possible with a close collaboration between clinical experts and microrobotics researchers to address the technical challenges in microfabrication, safety, and imaging. The clinical application potential can be materialized by designing microrobots that can solve the current main challenges, such as actuation limitations, material stability, and imaging constraints. The strengths and weaknesses of the current progress in the microrobotics field are discussed and a roadmap for their clinical applications in the near future is outlined.

Endovascular Treatment of Intracranial Aneurysm: The Importance of the Rheological Model in Blood Flow Simulations

Hemodynamics in intracranial aneurysm strongly depends on the non-Newtonian blood behavior due to the large number of suspended cells and the ability of red blood cells to deform and aggregate. However, most numerical investigations on intracranial hemodynamics adopt the Newtonian hypothesis to model blood flow and predict aneurysm occlusion. The aim of this study was to analyze the effect of the blood rheological model on the hemodynamics of intracranial aneurysms in the presence or absence of endovascular treatment. A numerical investigation was performed under pulsatile flow conditions in a patient-specific aneurysm with and without the insertion of an appropriately reconstructed flow diverter stent (FDS). The numerical simulations were performed using Newtonian and non-Newtonian assumptions for blood rheology. In all cases, FDS placement reduced the intra-aneurysmal velocity and increased the relative residence time (RRT) on the aneurysmal wall, indicating progressive thrombus formation and aneurysm occlusion. However, the Newtonian model largely overestimated RRT values and consequent aneurysm healing with respect to the non-Newtonian models. Due to the non-Newtonian blood properties and the large discrepancy between Newtonian and non-Newtonian simulations, the Newtonian hypothesis should not be used in the study of the hemodynamics of intracranial aneurysm, especially in the presence of endovascular treatment.

Standardized viscosity as a source of error in computational fluid dynamic simulations of cerebral aneurysms

BACKGROUND

Computational fluid dynamics (CFD) simulations are a powerful tool for studying cerebral aneurysms, capable of evaluating hemodynamics in a way that is infeasible with imaging alone. However, the difficulty of incorporating patient-specific information and inherent obstacles of in vivo validation have limited the clinical usefulness of CFD of cerebral aneurysms. In this work we investigate the effect of using standardized blood viscosity values in CFD simulations of cerebral aneurysms when compared to simulations of the same aneurysms using patient-specific viscosity values derived from hematocrit measurements.

PURPOSE

The objective of this work is to determine the level of error, on average, that is, caused by using standardized values of viscosity in CFD simulations of cerebral aneurysms. By quantifying this error, we demonstrate the need for incorporating patient-specific viscosity in future CFD investigations of cerebral aneurysms.

METHODS

CFD simulations of forty-one cerebral aneurysms were conducted using patient-specific boundary conditions. For each aneurysm two simulations were conducted, one utilizing patient-specific blood viscosity derived from hematocrit measurements and another using a standardized value for blood viscosity. Hemodynamic parameters such as wall shear stress (WSS), wall shear stress gradient (WSSG), and the oscillatory shear index (OSI) were calculated for each of the simulations for each aneurysm. Paired t-tests for differences in the time-averaged maps of these hemodynamic parameters between standardized and patient-specific viscosity simulations were conducted for each aneurysm. Bland-Altman analysis was used to examine the cohort-wide changes in the hemodynamic parameters. Subjects were broken into two groups, those with higher than standard viscosity and those with lower than standard viscosity. An unpaired t-test was used to compare the percent change in WSS, WSSG, and OSI between patient-specific and standardized viscosity simulations for the two cohorts. The percent changes in hemodynamic parameters were correlated against the direction and magnitude of percent change in viscosity, aneurysm size, and aneurysm location. For all t-tests, a Bonferroni-corrected significance level of 0.0167 was used.

RESULTS

63.2%, 41.5%, and 48.7% of aneurysms showed statistically significant differences between patient-specific and standardized viscosity simulations for WSS, WSSG, and OSI respectively. No statistically significant difference was found in the percent changes in WSS, WSSG, and OSI between the group with higher than standard viscosity and those with lower than standard viscosity, indicating an increase in viscosity can cause either an increase or decrease in each of the hemodynamic parameters. On a study-wide level no significant bias was found in either direction for WSS, WSSG, or OSI between the simulation groups due to the bidirectional effect of changing viscosity. No correlation was found between percent change of viscosity and percent change of WSS, WSSG, or OSI, meaning an after-the-fact correction for patient-specific viscosity is not feasible.

CONCLUSION

Standardizing viscosity values in CFD of cerebral aneurysms has a large and unpredictable impact on the calculated WSS, WSSG, and OSI when compared to CFD simulations of the same aneurysms using a patient-specific viscosity. We recommend implementing hematocrit-based patient-specific blood viscosity values for all CFD simulations of cerebral aneurysms.

Study of the mechanism of action of sand therapy on atherosclerosis based on the two-phase flow-Casson model

BACKGROUND

Sand therapy is a non-pharmacological physiotherapy method that uses the natural environment and resources of Xinjiang to treat through the heat transfer and magnetic effects of sand.

OBJECTIVE

Employing the two-phase flow-Casson blood flow model, we investigate the mechanism of atherosclerosis prevention via sand therapy, offering a biomechanical theoretical rationale for the prevention of atherosclerosis through sand therapy via the prism of computational fluid dynamics (CFD).

METHODS

Sand therapy experiments were conducted to obtain popliteal artery blood flow velocity, and blood was considered as a two-phase flow composed of plasma and red blood cells, and CFD method was applied to analyze the hemodynamic effects of Casson's blood viscosity model before and after sand therapy.

RESULTS

(1) The blood flow velocity increased by 0.24 m/s and 0.04 m/s at peak systolic and diastolic phases, respectively, after sand therapy; the axial velocity of blood vessels increased by 28.56% after sand therapy. (2) The average red blood cell viscosity decreased by 0.00014 Pa ⋅ s after sand therapy. (3) The low wall shear stress increased by 1.09 Pa and the high wall shear stress reached 41.47 Pa after sand therapy. (4) The time-averaged wall shear stress, shear oscillation index and relative retention time were reduced after sand therapy.

CONCLUSION

The increase of blood flow velocity after sand therapy can reduce the excessive deposition of cholesterol and other substances, the decrease of erythrocyte viscosity is beneficial to the migration of erythrocytes to the vascular center, the increase of low wall shear stress has a positive effect on the prevention of atherosclerosis, and the decrease of time-averaged wall shear stress, shear oscillation index and relative retention time can reduce the occurrence of thrombosis.

Computational Assessment of Unsteady Flow Effects on Magnetic Nanoparticle Targeting Efficiency in a Magnetic Stented Carotid Bifurcation Artery

PURPOSE

Worldwide, cardiovascular disease is the leading cause of hospitalization and death. Recently, the use of magnetizable nanoparticles for medical drug delivery has received much attention for potential treatment of both cancer and cardiovascular disease. However, proper understanding of the interacting magnetic field forces and the hydrodynamics of blood flow is needed for effective implementation. This paper presents the computational results of simulated implant assisted medical drug targeting (IA-MDT) via induced magnetism intended for administering patient specific doses of therapeutic agents to specific sites in the cardiovascular system. The drug delivery scheme presented in this paper functions via placement of a faintly magnetizable stent at a diseased location in the carotid artery, followed by delivery of magnetically susceptible drug carriers guided by the local magnetic field. Using this method, the magnetic stent can apply high localized magnetic field gradients within the diseased artery, while only exposing the neighboring tissues, arteries, and organs to a modest magnetic field. The localized field gradients also produce the forces needed to attract and hold drug-containing magnetic nanoparticles at the implant site for delivering therapeutic agents to treat in-stent restenosis.

METHODS

The multi-physics computational model used in this work is from our previous work and has been slightly modified for the case scenario presented in this paper. The computational model is used to analyze pulsatile blood flow, particle motion, and particle capture efficiency in a magnetic stented region using the magnetic properties of magnetite (Fe3O4) and equations describing the magnetic forces acting on particles produced by an external cylindrical electromagnetic coil. The electromagnetic coil produces a uniform magnetic field in the computational arterial flow model domain, while both the particles and the implanted stent are paramagnetic. A Eulerian-Lagrangian technique is adopted to resolve the hemodynamic flow and the motion of particles under the influence of a range of magnetic field strengths (Br = 2T, 4T, 6T, and 8T). Particle diameter sizes of 10 nm-4 µm in diameter were evaluated. Two dimensionless numbers were evaluated in this work to characterize relative effects of Brownian motion (BM), magnetic force induced particle motion, and convective blood flow on particle motion.

RESULTS

The computational simulations demonstrate that the greatest particle capture efficiency results for particle diameters within the micron range of 0.7-4 µm, specifically in regions where flow separation and vortices are at a minimum. Similar to our previous work (which did not involve the use of a magnetic stent), it was also observed that the capture efficiency of particles decreases substantially with particle diameter, especially in the superparamagnetic regime. Contrary to our previous work, using a magnetic stent tripled the capture efficiency of superparamagnetic particles. The highest capture efficiency observed for superparamagnetic particles was 78% with an 8 T magnetic field strength and 65% with a 2 T magnetic field strength when analyzing 100 nm particles. For 10 nm particles and an 8 T magnetic field strength, the particle capture efficiency was 55% and for a 2 T magnetic field strength the particle capture efficiency was observed to be 43%. Furthermore, it was found that larger magnetic field strengths, large particle diameter sizes (1 µm and above), and slower blood flow velocity improves the particle capture efficiency. The distribution of captured particles on the vessel wall along the axial and azimuthal directions is also discussed. Results for captured particles on the vessel wall along the axial flow direction showed that the particle density decreased along the axial direction, especially after the stented region. For the entrance section of the stented region, the captured particle density distribution along the axial direction is large, corresponding to the center-symmetrical distribution of the magnetic force in that section.

CONCLUSION

The simulation results presented in this work have shown to yield favorable capture efficiencies for micron range particles and superparamagnetic particles using magnetized implants such as the stent discussed in this work. The results presented in this work justify further investigation of MDT as a treatment technique for cardiovascular disease.

Is Accurate Lumen Segmentation More Important than Outlet Boundary Condition in Image-Based Blood Flow Simulations for Intracranial Aneurysms?

PURPOSE

Image-based blood flow simulations are increasingly used to investigate the hemodynamics in intracranial aneurysms (IAs). However, a strong variability in segmentation approaches as well as the absence of individualized boundary conditions (BCs) influence the quality of these simulation results leading to imprecision and decreased reliability. This study aims to analyze these influences on relevant hemodynamic parameters within IAs.

METHODS

As a follow-up study of an international multiple aneurysms challenge, the segmentation results of five IAs differing in size and location were investigated. Specifically, five possible outlet BCs were considered in each of the IAs. These are comprised of the zero-pressure condition (BC1), a flow distribution based on Murray's law with the exponents n = 2 (BC2) and n = 3 (BC3) as well as two advanced flow-splitting models considering the real vessels by including circular cross sections (BC4) or anatomical cross sections (BC5), respectively. In total, 120 time-dependent blood flow simulations were analyzed qualitatively and quantitatively, focusing on five representative intra-aneurysmal flow and five shear parameters such as vorticity and wall shear stress.

RESULTS

The outlet BC variation revealed substantial differences. Higher shear stresses (up to Δ9.69 Pa), intrasaccular velocities (up to Δ0.15 m/s) and vorticities (up to Δ629.22 1/s) were detected when advanced flow-splitting was applied compared to the widely used zero-pressure BC. The tendency of outlets BCs to over- or underestimate hemodynamic parameters is consistent across different segmentations of a single aneurysm model. Segmentation-induced variability reaches Δ19.58 Pa, Δ0.42 m/s and Δ957.27 1/s, respectively. Excluding low fidelity segmentations, however, (a) reduces the deviation drastically (>43%) and (b) leads to a lower impact of the outlet BC on hemodynamic predictions.

CONCLUSION

With a more realistic lumen segmentation, the influence of the BC on the resulting hemodynamics is decreased. A realistic lumen segmentation can be ensured, e.g., by using high-resolved 2D images. Furthermore, the selection of an advanced outflow-splitting model is advised and the use of a zero-pressure BC and BC based on Murray's law with exponent n = 3 should be avoided.

Cerebral multi-autoregulation model based enhanced external counterpulsation treatment planning for cerebral ischemic stroke

Enhanced external counterpulsation (EECP) treatment for cerebral ischemic stroke patients with differing severity of stenosis, is subject to uncertainties due to the varying effects of the cerebral autoregulation mechanism on haemodynamics. The current study reports the development of a cerebral multi-autoregulation (MR) mathematical model, based on cerebral arteriole regulation of neurogenic, vascular smooth muscle reflex and shear stress mechanisms which takes into account the severity of stenosis. The model was evaluated by comparison to authentic clinical measurements of cerebral autoregulatory efficiency. Then it was applied to a 0D/3D geometric multi-scale haemodynamic model of a cerebral artery. Haemodynamic indicators were calculated under different pressurization durations of EECP to evaluate the efficacy for different stenosis lesions. Moderate stenosis of 50% to 60% produced excessive time-averaged wall shear stress in the distal area of the stenosis (>7 Pa) during prolonged pressurization and may result in damage to vascular endothelial cells. However, prolonged pressurization did not result in haemodynamic risk for severe stenosis of 70% to 80%, indicating that the duration of pressurization may be extended with increasing severity of stenosis. The current MR model accurately simulated cerebral blood flow and has relevance to the simulation of cerebral haemodynamics in a clinical setting.

Do the stent blind insertion into the main portal vein (MPV) and stent diameter influence the surgical outcome of the transjugular intrahepatic portosystemic shunt (TIPS)?

Transjugular intrahepatic portosystemic shunt (TIPS) surgery is a clinical intervention to treat portal hypertension (PH) by deploying a covered stent to establish a shunt path for the portal vein (PV) system, and proper surgical strategy is of great importance to balance the shunt effect and the risk of complications. To understand the clinical strategies of the stent blind insertion and stent selection in clinic, this study investigated the effects of varying stent insertion positions and diameters on the PV hemodynamics and the shunt effect by computational fluid dynamics (CFD) analysis of five post-TIPS subjects. The results showed that the successful TIPS surgeries of the five PH subjects were confirmed by quantifying their pressure drops. The stent insertion positions at the main portal vein (MPV) slightly affected the clinically concerned hemodynamic indexes (i.e., MPV pressure, stent-outlet velocity) and the shunt index (SI). This indicated that the position of the stent going into the MPV may not need to be deliberately selected. Moreover, the covered stents with 6 mm and 8 mm diameters slightly influenced the hemodynamics as well, but the large-diameter stent better improved the shunt effect compared to the small-diameter one. Despite this, the 6 mm stent was suggested thanks to the higher risk of the hepatic encephalopathy (HE) observed in clinic, which indicated the excessive shunt of the 8 mm stent. The current work revealed the effects of different TIPS strategies on the surgical outcome, and could be useful for potential clinical practices.

Numerical Study of a Thrombus Migration Risk in Aneurysm After Coil Embolization in Patient Cases: FSI Modelling

PURPOSE

There are still many challenges for modelling a thrombus migration process in aneurysms. The main novelty of the present research lies in the modelling of aneurysm clot migration process in a realistic cerebral aneurysm, and the analysis of forces suffered by clots inside an aneurysm, through transient FSI simulations.

METHODS

The blood flow has been modelled using a Womersley velocity profile, and following the Carreau viscosity model. Hyperelastic Ogden model has been used for clot and isotropic linear elastic model for the artery walls. The FSI coupled model was implemented in ANSYS® software. The hemodynamic forces suffered by the clot have been quantified using eight different clot sizes and positions inside a real aneurysm.

RESULTS

The obtained results have shown that it is almost impossible for clots adjacent to aneurysm walls, to leave the aneurysm. Nevertheless, in clots positioned in the centre of the aneurysm, there is a real risk of clot migration. The risk of migration of a typical post-coiling intervention clot in an aneurysm, in contact with the wall and occupying a significant percentage of its volume is very low in the case studied, even in the presence of abnormally intense events, associated with sneezes or impacts.

CONCLUSIONS

The proposed methodology allows evaluating the clot migration risk, vital for evaluating the progress after endovascular interventions, it is a step forward in the personalized medicine, patient follow-up, and helping the medical team deciding the optimal treatment.

Assessment of the flow-diverter efficacy for intracranial aneurysm treatment considering pre- and post-interventional hemodynamics

Endovascular treatment of intracranial aneurysms with flow diverters (FD) has become one of the most promising interventions. Due to its woven high-density structure they are particularly applicable for challenging lesions. Although several studies have already conducted realistic hemodynamic quantification of the FD efficacy, a comparison with morphologic post-interventional data is still missing. This study analyses the hemodynamics of ten intracranial aneurysm patients treated with a novel FD device. Based on pre- and post-interventional 3D digital subtraction angiography image data, patient-specific 3D models of both treatment states are generated applying open source threshold-based segmentation methods. Using a fast virtual stenting approach, the real stent positions available in the post-interventional data are virtually replicated and both treatment scenarios were characterized using image-based blood flow simulations. The results show FD-induced flow reductions at the ostium by a decrease in mean neck flow rate (51%), inflow concentration index (56%) and mean inflow velocity (53%). Intraluminal reductions in flow activity for time-averaged wall shear stress (47%) and kinetic energy (71%) are present as well. However, an intra-aneurysmal increase in flow pulsatility (16%) for the post-interventional cases can be observed. Patient-specific FD simulations demonstrate the desired flow redirection and activity reduction inside the aneurysm beneficial for thrombosis formation. Differences in the magnitude of hemodynamic reduction exist over the cardiac cycle which may be addressed in a clinical setting by anti-hypertensive treatment in selected cases.

Classification of Blood Rheological Models through an Idealized Symmetrical Bifurcation

The assumed rheological behavior of blood influences the hemodynamic characteristics of numerical blood flow simulations. Until now, alternative rheological specifications have been utilized, with uncertain implications for the results obtained. This work aims to group sixteen blood rheological models in homogeneous clusters, by exploiting data generated from numerical simulations on an idealized symmetrical arterial bifurcation. Blood flow is assumed to be pulsatile and is simulated using a commercial finite volume solver. An appropriate mesh convergence study is performed, and all results are collected at three different time instants throughout the cardiac cycle: at peak systole, early diastole, and late diastole. Six hemodynamic variables are computed: the time average wall shear stress, oscillatory shear index, relative residence time, global and local non-Newtonian importance factor, and non-Newtonian effect factor. The resulting data are analyzed using hierarchical agglomerative clustering algorithms, which constitute typical unsupervised classification methods. Interestingly, the rheological models can be partitioned into three homogeneous groups, whereas three specifications appear as outliers which do not belong in any partition. Our findings suggest that models which are defined in a similar manner from a mathematical perspective may behave substantially differently in terms of the data they produce. On the other hand, models characterized by different mathematical formulations may belong to the same statistical group (cluster) and can thus be considered interchangeably.

EFFECTS OF SAND THERAPY ON FLUID DYNAMICS OF BLOOD IN THE POPLITEAL ARTERY BASED ON THE EULER MULTIPHASE MODEL

The effect of sand therapy on the hemodynamics of the popliteal artery (POA) was investigated to elucidate its mechanism in atherosclerosis physiotherapy. The hemodynamic parameters of the POA before and after sand therapy were obtained from 58 subjects and statistically analyzed. A three-dimensional finite element model of POA bifurcation was established. Blood was regarded as a triple-phase flow composed of plasma, red blood cells (RBCs), and white blood cells (WBCs). Using computational fluid dynamics (CFD), blood flow velocity, wall shear stress (WSS), and blood cell volume distribution before and after sand therapy were calculated. Sand therapy could reduce the vortex and reflux in blood vessels. The maximum increase in blood axial velocity was 34.38%. The maximum WSS increased by 22.82 Pa and the blood cells gradually moved away from the vessel wall and toward the axis of the blood vessel. The reduction in vortex and reflux, and the increase in blood axial velocity after sand therapy are beneficial in reducing the accumulation of cholesterol and other substances, thereby alleviating the formation of atherosclerotic plaques and lipid streaks. The rise in WSS slows the tendency of arterial vascular thickening in regions of low WSS. The movement of blood cells to the axis of blood vessels can reduce the possibility of thrombosis caused by blood cells squeezing on the blood vessel walls. These findings suggest that sand therapy may slow atherosclerosis progression.

Cerebral Aneurysm Occlusion at 12-Month Follow-Up After Flow-Diverter Treatment: Statistical Modeling for V&V With Real-World Data

Background: Flow-Diverter (FD) porosity has been pointed as a critical factor in the occlusion of cerebral aneurysms after treatment.

Objective: Verification and Validation of computational models in terms of predictive capacity, relating FD porosity and occlusion after cerebral aneurysms treatment.

Methods: Sixty-four aneurysms, with pre-treatment and follow-up images, were considered. Patient demographics and aneurysm morphological information were collected. The computational simulation provided by ANKYRAS provided FD porosity, expansion, and mesh angle. FD occlusion was assessed and recorded from follow-up images. Multiple regression Logit and analysis of covariance (ANCOVA) models were used to model the data with both categorical and continuous models.

Results: Occlusion of the aneurysm after 12 months was affected by aneurysm morphology but not by FD mesh morphology. A Time-To-Occlusion (TTO) of 6.92 months on average was observed with an SE of 0.24 months in the aneurysm population surveyed. TTO was estimated with statistical significance from the resulting model for the data examined and was capable of explaining 92% of the data variation.

Conclusions: Porosity was found to have the most correction power when assessing TTO, proving its importance in the process of aneurysm occlusion. Still, further Verification and Validation (V&V) of treatment simulation in more extensive, multi-center, and randomized databases is required.

Role of patient-specific blood properties in computational fluid dynamics simulation of flow diverter deployed cerebral aneurysms

BACKGROUND:

Hemodynamics and their clinical outcome of cerebral aneurysms treated with flow diverter (FD) stents have thus far been investigated using computational fluid dynamics (CFD) simulations. Although human blood is characterized as a non-Newtonian patientspecific fluid, non-patient-specific blood properties (PSBP) were applied in most extant studies.

OBJECTIVE:

To investigate the hemodynamic effects caused by PSBPs in aneurysms treated with FD stents.

METHODS:

We measured blood properties (density and viscosity) for 12 patients who underwent FD stent deployment. We conducted CFD simulations with the measured PSBPs and non-PSBPs quoted from previous studies. The average blood flow velocity and wall shear stress within the aneurysms were calculated and two simulation patterns were compared.

RESULTS:

The velocity and wall shear stress changed by 2.93% and 3.16% on average, respectively, without an FD stent deployed. Conversely, with the FD stents deployed, the change rates increased to 11.1% and 9.06% on average, respectively.

CONCLUSIONS:

The change in hemodynamic parameters if PSBPs are considered, may not be negligible when conducting CFD simulations of FD stent deployed aneurysms To obtain an adequate hemodynamic environment for cerebral aneurysms with FD stents deployed, it is recommended to use PSBPs for CFD simulations.

A nonlinear multi-scale model for blood circulation in a realistic vascular system

In the last decade, numerical models have become an increasingly important tool in biological and medical science. Numerical simulations contribute to a deeper understanding of physiology and are a powerful tool for better diagnostics and treatment. In this paper, a nonlinear multi-scale model framework is developed for blood flow distribution in the full vascular system of an organ. We couple a quasi one-dimensional vascular graph model to represent blood flow in larger vessels and a porous media model to describe flow in smaller vessels and capillary bed. The vascular model is based on Poiseuille's Law, with pressure correction by elasticity and pressure drop estimation at vessels' junctions. The porous capillary bed is modelled as a two-compartment domain (artery and venous) using Darcy's Law. The fluid exchange between the artery and venous capillary bed compartments is defined as blood perfusion. The numerical experiments show that the proposed model for blood circulation: (i) is closely dependent on the structure and parameters of both the larger vessels and of the capillary bed, and (ii) provides a realistic blood circulation in the organ. The advantage of the proposed model is that it is complex enough to reliably capture the main underlying physiological function, yet highly flexible as it offers the possibility of incorporating various local effects. Furthermore, the numerical implementation of the model is straightforward and allows for simulations on a regular desktop computer.

Cerebral aneurysm flow diverter modeled as a thin inhomogeneous porous medium in hemodynamic simulations

Rapid and accurate simulation of cerebral aneurysm flow modifications by flow diverters (FDs) can help improving patient-specific intervention and predicting treatment outcome. However, when FD devices are explicitly represented in computational fluid dynamics (CFD) simulations, flow around the stent wires must be resolved, leading to high computational cost. Classic porous medium (PM) methods can reduce computational expense but cannot capture the inhomogeneous FD wire distribution once implanted on a cerebral artery and thus cannot accurately model the post-stenting aneurysmal flow. We report a novel approach that models the FD flow modification as a thin inhomogeneous porous medium (iPM). It improves over the classic PM approaches in two ways. First, the FD is more appropriately treated as a thin screen, which is more accurate than the classic 3D-PM-based Darcy-Forchheimer relation. Second, pressure drop is calculated cell-by-cell using the local FD geometric parameters across an inhomogeneous PM. We applied the iPM technique to simulating the post-stenting hemodynamics of three patient-specific aneurysms. To test its accuracy and speed, we compared the results from the iPM technique against CFD simulations with explicit FD devices. The iPM CFD ran 500% faster than the explicit CFD while achieving 94%-99% accuracy; thus, iPM is a promising clinical bedside modeling tool to assist endovascular interventions with FD and stents.

Hemodynamic Investigation of the Effectiveness of a Two Overlapping Flow Diverter Configuration for Cerebral Aneurysm Treatment

Flow diverters (FDs) are widely employed as endovascular treatment devices for large or wide-neck cerebral aneurysms. Occasionally, overlapped FDs are deployed to enhance the flow diversion effect. In this study, we investigated the hemodynamics of overlapping FDs via computational fluid dynamics (CFD) simulations. We reproduced the arterial geometry of a patient who had experienced the deployment of two overlapping FDs. We utilized two stent patterns, namely the patterns for one FD and two overlapping FDs. We calculated the velocity, mass flow rate, wall shear stress, and pressure loss coefficient as well as their change rates for each pattern relative to the no-FD pattern results. The CFD simulation results indicated that the characteristics of the blood flow inside the aneurysm were minimally affected by the deployment of a single FD; in contrast, the overlapping FD pattern results revealed significant changes in the flow. Further, the velocity at an inspection plane within the aneurysm sac decreased by up to 92.2% and 31.0% in the cases of the overlapping and single FD patterns, respectively, relative to the no-FD pattern. The simulations successfully reproduced the hemodynamics, and the qualitative and quantitative investigations are meaningful with regard to the clinical outcomes of overlapped FD deployment.

Comparison of Newtonian and Non-newtonian Fluid Models in Blood Flow Simulation in Patients With Intracranial Arterial Stenosis

Background

Newtonian fluid model has been commonly applied in simulating cerebral blood flow in intracranial atherosclerotic stenosis (ICAS) cases using computational fluid dynamics (CFD) modeling, while blood is a shear-thinning non-Newtonian fluid. We aimed to investigate the differences of cerebral hemodynamic metrics quantified in CFD models built with Newtonian and non-Newtonian fluid assumptions, in patients with ICAS.

Methods

We built a virtual artery model with an eccentric 75% stenosis and performed static CFD simulation. We also constructed CFD models in three patients with ICAS of different severities in the luminal stenosis. We performed static simulations on these models with Newtonian and two non-Newtonian (Casson and Carreau-Yasuda) fluid models. We also performed transient simulations on another patient-specific model. We measured translesional pressure ratio (PR) and wall shear stress (WSS) values in all CFD models, to reflect the changes in pressure and WSS across a stenotic lesion. In all the simulations, we compared the PR and WSS values in CFD models derived with Newtonian, Casson, and Carreau-Yasuda fluid assumptions.

Results

In all the static and transient simulations, the Newtonian/non-Newtonian difference on PR value was negligible. As to WSS, in static models (virtual and patient-specific), the rheological difference was not obvious in areas with high WSS, but observable in low WSS areas. In the transient model, the rheological difference of WSS areas with low WSS was enhanced, especially during diastolic period.

Conclusion

Newtonian fluid model could be applicable for PR calculation, but caution needs to be taken when using the Newtonian assumption in simulating WSS especially in severe ICAS cases.

Hemodynamic and Morphological Factors Related to Coil Compaction in Basilar Artery Tip Aneurysms

PURPOSE

Coil compaction is directly related to the degree of cerebral aneurysmal recanalization. The Degree of Recanalization (DoR) was quantified by measuring the volume vacated by coil deformation. The purpose of this study was to clarify the hemodynamic and morphological factors associated with coil compaction.

METHODS

Computational fluid dynamics (CFD) simulations were performed on 28 middle size (5-10 mm), unruptured basilar artery tip aneurysms. The DoR was measured by comparing the coil mass shape obtained from three-dimensional digital subtraction angiography data immediately after coil embolization and again within 1 to 2 years of follow-up. Deployed coils were modeled using a virtual coiling technique for CFD simulations. Hemodynamic and morphological factors to predict the DoR were derived using multiple linear regression.

RESULTS

Aneurysmal neck area, the maximum pressure generated on the neck surface after coil embolization, and the high-pressure position on the neck surface predicted DoR with statistical significance (p<0.001, p<0.001, p=0.004, respectively). DoR tended to increase when the neck area was large, the pressure generated on the coils was high, and the high-pressure position was close to the center of the neck surface. The volume embolization ratio was not statistically relevant for the DoR in the cases of this study.

CONCLUSIONS

Coil compaction occurs in cerebral aneurysms with a wide neck, high pressure generated on the coils, and high pressure in the center of the neck surface. Establishing the DoR can contribute to the prediction of recanalization after coil embolization.

Analysis of Morphological-Hemodynamic Risk Factors for Aneurysm Rupture Including a Newly Introduced Total Volume Ratio

The purpose of this study was to evaluate morphological and hemodynamic factors, including the newly developed total volume ratio (TVR), in evaluating rupture risk of cerebral aneurysms using ≥7 mm sized aneurysms. Twenty-three aneurysms (11 unruptured and 12 ruptured) ≥ 7 mm were analyzed from 3-dimensional rotational cerebral angiography and computational fluid dynamics (CFD). Ten morphological and eleven hemodynamic factors of the aneurysms were qualitatively and quantitatively compared. Correlation analysis between morphological and hemodynamic factors was performed, and the relationship among the hemodynamic factors was analyzed. Morphological factors (ostium diameter, ostium area, aspect ratio, and bottleneck ratio) and hemodynamic factors (TVR, minimal wall shear stress of aneurysms, time-averaged wall shear stress of aneurysms, oscillatory shear index, relative residence time, low wall shear stress area, and ratio of low wall stress area) were statistically different between ruptured and unruptured aneurysms (p < 0.05). By simple regression analysis, the morphological factor aspect ratio and the hemodynamic factor TVR were significantly correlated (r² = 0.602, p = 0.001). Ruptured aneurysms had complex and unstable flow. In ≥7 mm ruptured aneurysms, high aspect ratio, bottleneck ratio, complex flow, unstable flow, low TVR, wall shear stress at aneurysm, high oscillatory shear index, relative resistance time, low wall shear stress area, and ratio of low wall stress area were significant in determining the risk of aneurysm rupture.

Intra-saccular device modeling for treatment planning of intracranial aneurysms: from morphology to hemodynamics

MOTIVATION

Intra-saccular devices (ID), developed for the treatment of bifurcation aneurysms, offer new alternatives for treating complex terminal and bifurcation aneurysms. In this work, a complete workflow going from medical images to post-treatment CFD analysis is described and used in the assessment of a concrete clinical problem.

MATERIALS AND METHODS

Two different intra-saccular device sizes were virtually implanted in 3D models of the patient vasculature using the ID-Fit method. After deployment, the local porosity at the closed end of the device in contact with the blood flow was computed. This porosity was then used to produce a CFD porous medium model of the device. Velocities and wall shear stress were assessed for each model.

RESULTS

Six patients treated with intra-saccular devices were included in this work. For each case, 2 different device sizes were virtually implanted and 3 CFD simulations were performed: after deployment simulation with each size and before deployment simulation (untreated). A visible reduction in velocities was observed after device implantation. Velocity and WSS reduction was statistically significant (K-S statistics, [Formula: see text]).

CONCLUSIONS

Placement of different device size can lead to a partial filling of the aneurysm, either at the dome or at the neck, depending on the particular positioning by the interventionist. The methodology used in this work can have a strong clinical impact, since it provides additional information in the process of device selection using preoperative data.

Fluid-structure interaction modeling of lactating breast: Newtonian vs. non-Newtonian milk

Breastfeeding is a highly dynamic and complex mechanism. The suckling process by the infant involves compression and intra-oral vacuum pressure, leading to milk expression from breast. The accumulated milk from the nipple varies depending on the milk properties and transient flow rate during the suckling cycle. Rheological studies on raw human milk indicate that milk has a non-Newtonian shear-thinning flow behavior. This study aims to investigate the effect of non-Newtonian milk on flow behavior through the breast ductal system using fluid-structure interaction (FSI) simulation. The results of the non-Newtonian effects on flow velocity and the volumetric flow rate of expressed milk are presented. The results show that non-Newtonian Carreau model is promising for the simulation of human milk flow through the breast ductal systems. Also, the results show that the non-Newtonian effects on the milk flow behavior appear for 30-35% of the suckling cycle. Therefore, the Newtonian model is acceptable for the purpose of numerical simulation.

Assessment and visualization of hemodynamic loading in aneurysm sac and neck: Effect of foam insertion

Shape memory polymer (SMP) foam is often proposed as the future alternative of coils in aneurysm treatment devices. Present work numerically investigates the unsteady, three-dimensional simulation of blood flow in a cerebral aneurysm filled with SMP foam. Simulations are conducted on patient-specific geometries with realistic blood velocity waveform imposed at the inlet while SMP foam is treated as a porous medium. The present study introduces a "loading risk map" that helps to visualize the hemodynamic effect of foam insertion on the aneurysm sac and neck. The loading risk maps suggest that while the SMP foam subdues the flow and wall shear pulsations in the aneurysm sac, the pressure distribution is minimally affected. The maps suggest that while the downstream lip is the most risk-prone site for both geometries, downstream vascular anatomy significantly influences foam efficiency in reducing pressure and wall shear stress loading.

Numerical study on the energy cascade of pulsatile Newtonian and power-law flow models in an ICA bifurcation

The complex physics and biology underlying intracranial hemodynamics are yet to be fully revealed. A fully resolved direct numerical simulation (DNS) study has been performed to identify the intrinsic flow dynamics in an idealized carotid bifurcation model. To shed the light on the significance of considering blood shear-thinning properties, the power-law model is compared to the commonly used Newtonian viscosity hypothesis. We scrutinize the kinetic energy cascade (KEC) rates in the Fourier domain and the vortex structure of both fluid models and examine the impact of the power-law viscosity model. The flow intrinsically contains coherent structures which has frequencies corresponding to the boundary frequency, which could be associated with the regulation of endothelial cells. From the proposed comparative study, it is found that KEC rates and the vortex-identification are significantly influenced by the shear-thinning blood properties. Conclusively, from the obtained results, it is found that neglecting the non-Newtonian behavior could lead to underestimation of the hemodynamic parameters at low Reynolds number and overestimation of the hemodynamic parameters by increasing the Reynolds number. In addition, we provide physical insight and discussion onto the hemodynamics associated with endothelial dysfunction which plays significant role in the pathogenesis of intracranial aneurysms.

Stenting as porous media in anatomically accurate geometries. A comparison of models and spatial heterogeneity

Modelling intracranial aneurysm blood flow after flow diverter treatment has proven to be of great scientific and clinical interest. One of the reasons for not having CFD as an everyday clinical tool yet is the time required to set-up such simulations plus the required computational time. The speed-up of these simulations can have a considerable impact during treatment planning and device selection. Modelling flow diverters as a porous medium (PM) can considerably improve the computational time. Many models have been presented in literature, but quantitative comparisons between models are scarce. In this study, the untreated case, the explicit definition of the flow diverter wires as no-slip boundary condition and five different porous medium models were chosen for comparison, and evaluated on intracranial aneurysm of 14 patients with different shapes, sizes, and locations. CFD simulations were made using finite volume method on steady flow conditions. Velocities, kinetic energy, wall shear stress, and computational time were assessed for each model. Then, all models are compared against the no-slip boundary condition using non parametric Kolmogorov-Smirnov test. The model with least performance showed a mean K-S statistic of 0.31 and deviance of 0.2, while the model with best values always gave K-S statistics below 0.2. Kinetic energy between PM models varied between an over estimation of 218.3% and an under estimation of 73.06%. Also, speedups were between 4.75x and 5.3x (stdev: 0.38x and 0.15x) when using PM models. Flow diverters can be simulated with PM with a good agreement to standard CFD simulations were FD wires are represented with no-slip boundary condition in less than a quarter of the time. Best results were obtained on PM models based on geometrical properties, in particular, when using a heterogeneous medium based on equations for flat rhomboidal wire frames.

Hemodynamic and morphological characteristics of a growing cerebral aneurysm

The growth of cerebral aneurysms is linked to local hemodynamic conditions, but the driving mechanisms of the growth are poorly understood. The goal of this study was to examine the association between intraaneurysmal hemodynamic features and areas of aneurysm growth, to present the key hemodynamic parameters essential for an accurate prediction of the growth, and to gain a deeper understanding of the underlying mechanisms. Patient-specific images of a growing cerebral aneurysm in 3 different growth stages acquired over a period of 40 months were segmented and reconstructed. A unique aspect of this patient-specific case study was that while one side of the aneurysm stayed stable, the other side continued to grow. This unique case enabled the authors to examine their aims in the same patient with parent and daughter arteries under the same inlet flow conditions. Pulsatile flow in the aneurysm models was simulated using computational fluid dynamics and was validated with in vitro experiments using particle image velocimetry measurements. The authors' detailed analysis of intrasaccular hemodynamics linked the growing regions of aneurysms to flow instabilities and complex vortex structures. Extremely low velocities were observed at or around the center of the unstable vortex structure, which matched well with the growing regions of the studied cerebral aneurysm. Furthermore, the authors observed that the aneurysm wall regions with a growth greater than 0.5 mm coincided with wall regions of lower (< 0.5 Pa) time-averaged wall shear stress (TAWSS), lower instantaneous (< 0.5 Pa) wall shear stress (WSS), and high (> 0.1) oscillatory shear index (OSI). To determine which set of parameters can best identify growing and nongrowing aneurysms, the authors performed statistical analysis for consecutive stages of the growing CA. The results demonstrated that the combination of TAWSS and the distance from the center of the vortical structure has the highest sensitivity and positive predictive value, and relatively high specificity and negative predictive value. These findings suggest that an unstable, recirculating flow structure within the aneurysm sac created in the region adjacent to the aneurysm wall with low TAWSS may be introduced as an accurate criterion to explain the hemodynamic conditions predisposing the aneurysm to growth. The authors' findings are based on one patient's data set, but the study lays out the justification for future large-scale verification. The authors' findings can assist clinicians in differentiating stable and growing aneurysms during preinterventional planning.

A review on the reliability of hemodynamic modeling in intracranial aneurysms: why computational fluid dynamics alone cannot solve the equation

Computational blood flow modeling in intracranial aneurysms (IAs) has enormous potential for the assessment of highly resolved hemodynamics and derived wall stresses. This results in an improved knowledge in important research fields, such as rupture risk assessment and treatment optimization. However, due to the requirement of assumptions and simplifications, its applicability in a clinical context remains limited.This review article focuses on the main aspects along the interdisciplinary modeling chain and highlights the circumstance that computational fluid dynamics (CFD) simulations are embedded in a multiprocess workflow. These aspects include imaging-related steps, the setup of realistic hemodynamic simulations, and the analysis of multidimensional computational results. To condense the broad knowledge, specific recommendations are provided at the end of each subsection.Overall, various individual substudies exist in the literature that have evaluated relevant technical aspects. In this regard, the importance of precise vessel segmentations for the simulation outcome is emphasized. Furthermore, the accuracy of the computational model strongly depends on the specific research question. Additionally, standardization in the context of flow analysis is required to enable an objective comparison of research findings and to avoid confusion within the medical community. Finally, uncertainty quantification and validation studies should always accompany numerical investigations.In conclusion, this review aims for an improved awareness among physicians regarding potential sources of error in hemodynamic modeling for IAs. Although CFD is a powerful methodology, it cannot provide reliable information, if pre- and postsimulation steps are inaccurately carried out. From this, future studies can be critically evaluated and real benefits can be differentiated from results that have been acquired based on technically inaccurate procedures.

Computational and experimental investigation of particulate matter deposition in cerebral side aneurysms

Intracranial aneurysms frequently develop blood clots, plaque and inflammations, which are linked to enhanced particulate mass deposition. In this work, we propose a computational model for particulate deposition, that accounts for the influence of field forces, such as gravity and electrostatics, which produce an additional flux of particles perpendicular to the fluid motion and towards the wall. This field-mediated flux can significantly enhance particle deposition in low-shear environments, such as in aneurysm cavities. Experimental investigation of particle deposition patterns in in vitro models of side aneurysms, demonstrated the ability of the model to predict enhanced particle adhesion at these sites. Our results showed a significant influence of gravity and electrostatic forces (greater than 10%), indicating that the additional terms presented in our models may be necessary for modelling a wide range of physiological flow conditions and not only for ultra-low shear regions. Spatial differences between the computational model and the experimental results suggested that additional transport and fluidic mechanisms affect the deposition pattern within aneurysms. Taken together, the presented findings may enhance our understanding of pathological deposition processes at cardiovascular disease sites, and facilitate rational design and optimization of cardiovascular particulate drug carriers.

Long-Term Venovenous Connection for Extracorporeal Carbon Dioxide Removal (ECCO2R)-Numerical Investigation of the Connection to the Common Iliac Veins

Purpose

Currently used cannulae for extracorporeal carbon dioxide removal (ECCO₂R) are associated with complications such as thrombosis and distal limb ischemia, especially for long-term use. We hypothesize that the risk of these complications is reducible by attaching hemodynamically optimized grafts to the patient’s vessels. In this study, as a first step towards a long-term stable ECCO₂R connection, we investigated the feasibility of a venovenous connection to the common iliac veins. To ensure its applicability, the drainage of reinfused blood (recirculation) and high wall shear stress (WSS) must be avoided.

Methods

A reference model was selected for computational fluid dynamics, on the basis of the analysis of imaging data. Initially, a sensitivity analysis regarding recirculation was conducted using as variables: blood flow, the distance of drainage and return to the iliocaval junction, as well as the diameter and position of the grafts. Subsequently, the connection was optimized regarding recirculation and the WSS was evaluated. We validated the simulations in a silicone model traversed by dyed fluid.

Results

The simulations were in good agreement with the validation measurements (mean deviation 1.64%). The recirculation ranged from 32.1 to 0%. The maximum WSS did not exceed 5.57 Pa. The position and diameter of the return graft show the highest influence on recirculation. A correlation was ascertained between recirculation and WSS. Overall, an inflow jet directed at a vessel wall entails not only high WSS, but also a flow separation and thereby an increased recirculation. Therefore, return grafts aligned to the vena cava are crucial.

Conclusion

In conclusion, a connection without recirculation could be feasible and therefore provides a promising option for a long-term ECCO₂R connection.

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.

Recanalization of Embolized Endovascular Intracranial Aneurysms and Changes in the Blood Viscosity: A Pilot Study

Background

The purpose of our research was to evaluate the relationships between blood viscosity and recanalization of coiled intracranial aneurysms.

Material/Methods

The study included consecutives patients treated endovascularly by a team of experienced neurosurgeons and neuroradiologists due to brain aneurysm. A total of 50 patients (the average age was 57.48 years, SD=13.71) were assigned to 2 groups: group A with recanalization (4 male and 8 female patients) and group B without recanalization (10 male and 28 female patients) were examined. All patients underwent a 6-month follow-up of the whole-blood viscosity test with a Brookfield DV III+pro cone-plate viscometer using the Rheocalc program. Differences between groups were assessed using the Statistica 12 computer program (StatSoft Inc., Tulsa, OK, USA).

Results

Studies have shown no significant difference in the age range between group A and B (P=0.31). In group A, higher viscosity values were found for whole blood [median: 4.14 dyn×sec/cm² (mPa×sec) quartile range 0.42], compared to group B [median: 3.92 dyn×sec/cm² (mPa×sec); quartile range 0.40; (P=0.04)]. This difference was significant (P=0.04). Additionally, the level of hematocrit was positively related with recanalization, the higher the hematocrit, the more frequent recanalization. A very strong and statistically significant relationship occurred between the frequency of recanalization and smoking (P<0.001).

Conclusions

The occurrence of higher values of whole blood viscosity which increase turbulent flow through the vessels may be a risk for recanalization of the coiled intracranial aneurysm.

Size-Dependent Distribution of Patient-Specific Hemodynamic Factors in Unruptured Cerebral Aneurysms Using Computational Fluid Dynamics

PURPOSE

To analyze size-dependent hemodynamic factors [velocity, shear rate, blood viscosity, wall shear stress (WSS)] in unruptured cerebral aneurysms using computational fluid dynamics (CFD) based on the measured non-Newtonian model of viscosity.

METHODS

Twenty-one patients with unruptured aneurysms formed the study cohort. Patient-specific geometric models were reconstructed for CFD analyses. Aneurysms were divided into small and large groups based on a cutoff size of 5 mm. For comparison between small and large aneurysms, 5 morphologic variables were measured. Patient-specific non-Newtonian blood viscosity was applied for more detailed CFD simulation. Quantitative and qualitative analyses of velocity, shear rate, blood viscosity, and WSS were conducted to compare small and large aneurysms.

RESULTS

Complex flow patterns were found in large aneurysms. Large aneurysms had a significantly lower shear rate (235 ± 341 s^-1)) than small aneurysms (915 ± 432 s^-1) at peak-systole. Two times higher blood viscosity was observed in large aneurysms compared with small aneurysms. Lower WSS was found in large aneurysms (1.38 ± 1.36 Pa) than in small aneurysms (3.53 ± 1.22 Pa). All the differences in hemodynamic factors between small and large aneurysms were statistically significant.

CONCLUSIONS

Large aneurysms tended to have complex flow patterns, low shear rate, high blood viscosity, and low WSS. The hemodynamic factors that we analyzed might be useful for decision making before surgical treatment of aneurysms.

Role of distal cerebral vasculature in vessel constriction after aneurysm treatment with flow diverter stents

BACKGROUND

Treatment of intracranial aneurysms with flow diverter stent (FDS) procedures can lead to caliber changes of jailed vessels. The reason some branches remain unchanged and others are affected by narrowing remains unknown.

OBJECTIVE

To investigate the influence of resistance to flow from distal vasculature on stent-induced hemodynamic modifications affecting bifurcating vessels.

MATERIALS AND METHODS

Radiological images and demographic data were acquired for 142 aneurysms treated with a FDS. Vascular resistance was estimated from patient-specific anatomic data. Correlation analysis was used to identify correspondence between anatomic data and clinical outcome. Computational Fluid Dynamics was performed on a typical patient-specific model to evaluate the influence of FDS on flow. Relevant hemodynamic variables along the bifurcating vessels were quantitatively analyzed and validated with in vitro data obtained using power Doppler ultrasound.

RESULTS

Statistical analysis showed a correlation between clinical outcome and FDS resistance to flow considering overall jailed vessel vascular resistance (r=0.5, P<0.001). Computational predictions of blood flow showed that hemodynamics is minimally affected by FDS treatment in the ophthalmic artery.

CONCLUSIONS

Jailed vessels are affected by narrowing when resistance to flow from the FDS constitutes a larger proportion of the overall vessel resistance to flow. This knowledge may contribute to better understanding of intracranial hemodynamics after a FDS procedure and reinforce indications for flow diversion in the treatment of intracranial aneurysms.

Advances in Biomaterials and Technologies for Vascular Embolization

Minimally invasive transcatheter embolization is a common nonsurgical procedure in interventional radiology used for the deliberate occlusion of blood vessels for the treatment of diseased or injured vasculature. A wide variety of embolic agents including metallic coils, calibrated microspheres, and liquids are available for clinical practice. Additionally, advances in biomaterials, such as shape-memory foams, biodegradable polymers, and in situ gelling solutions have led to the development of novel preclinical embolic agents. The aim here is to provide a comprehensive overview of current and emerging technologies in endovascular embolization with respect to devices, materials, mechanisms, and design guidelines. Limitations and challenges in embolic materials are also discussed to promote advancement in the field.

Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH)-Phase Ib: Effect of morphology on hemodynamics

Background

Image-based blood flow simulations have been increasingly applied to investigate intracranial aneurysm (IA) hemodynamics. However, the acceptance among physicians remains limited due to the high variability in the underlying assumptions and quality of results.

Methods

To evaluate the vessel segmentation as one of the most important sources of error, the international Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH) was announced. 26 research groups from 13 different countries segmented three datasets, which contained five IAs in total. Based on these segmentations, 73 time-dependent blood flow simulations under consistent conditions were carried out. Afterwards, relevant flow and shear parameters (e.g., neck inflow rate, parent vessel flow rate, spatial mean velocity, and wall shear stress) were analyzed both qualitatively and quantitatively.

Results

Regarding the entire vasculature, the variability of the segmented vessel radius is 0.13 mm, consistent and independent of the local vessel radius. However, the centerline velocity shows increased variability in more distal vessels. Focusing on the aneurysms, clear differences in morphological and hemodynamic parameters were observed. The quantification of the segmentation-induced variability showed approximately a 14% difference among the groups for the parent vessel flow rate. Regarding the mean aneurysmal velocity and the neck inflow rate, a variation of 30% and 46% was observed, respectively. Finally, time-averaged wall shear stresses varied between 28% and 51%, depending on the aneurysm in question.

Conclusions

MATCH reveals the effect of state-of-the-art segmentation algorithms on subsequent hemodynamic simulations for IA research. The observed variations may lead to an inappropriate interpretation of the simulation results and thus, can lead to inappropriate conclusions by physicians. Therefore, accurate segmentation of the region of interest is necessary to obtain reliable and clinically helpful flow information.

Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH)-phase II: rupture risk assessment

PURPOSE

Assessing the rupture probability of intracranial aneurysms (IAs) remains challenging. Therefore, hemodynamic simulations are increasingly applied toward supporting physicians during treatment planning. However, due to several assumptions, the clinical acceptance of these methods remains limited.

METHODS

To provide an overview of state-of-the-art blood flow simulation capabilities, the Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH) was conducted. Seventeen research groups from all over the world performed segmentations and hemodynamic simulations to identify the ruptured aneurysm in a patient harboring five IAs. Although simulation setups revealed good similarity, clear differences exist with respect to the analysis of aneurysm shape and blood flow results. Most groups (12/71%) included morphological and hemodynamic parameters in their analysis, with aspect ratio and wall shear stress as the most popular candidates, respectively.

RESULTS

The majority of groups (7/41%) selected the largest aneurysm as being the ruptured one. Four (24%) of the participating groups were able to correctly select the ruptured aneurysm, while three groups (18%) ranked the ruptured aneurysm as the second most probable. Successful selections were based on the integration of clinically relevant information such as the aneurysm site, as well as advanced rupture probability models considering multiple parameters. Additionally, flow characteristics such as the quantification of inflow jets and the identification of multiple vortices led to correct predictions.

CONCLUSIONS

MATCH compares state-of-the-art image-based blood flow simulation approaches to assess the rupture risk of IAs. Furthermore, this challenge highlights the importance of multivariate analyses by combining clinically relevant metadata with advanced morphological and hemodynamic quantification.

Patient-Specific Computational Fluid Dynamics in Ruptured Posterior Communicating Aneurysms Using Measured Non-Newtonian Viscosity : A Preliminary Study

Objective

The objective of this study was to analyze patient-specific blood flow in ruptured aneurysms using obtained non-Newtonian viscosity and to observe associated hemodynamic features and morphological effects.

Methods

Five patients with acute subarachnoid hemorrhage caused by ruptured posterior communicating artery aneurysms were included in the study. Patients’ blood samples were measured immediately after enrollment. Computational fluid dynamics (CFD) was conducted to evaluate viscosity distributions and wall shear stress (WSS) distributions using a patient-specific geometric model and shear-thinning viscosity properties.

Results

Substantial viscosity change was found at the dome of the aneurysms studied when applying non-Newtonian blood viscosity measured at peak-systole and end-diastole. The maximal WSS of the non-Newtonian model on an aneurysm at peaksystole was approximately 16% lower compared to Newtonian fluid, and most of the hemodynamic features of Newtonian flow at the aneurysms were higher, except for minimal WSS value. However, the differences between the Newtonian and non-Newtonian flow were not statistically significant. Rupture point of an aneurysm showed low WSS regardless of Newtonian or non-Newtonian CFD analyses.

Conclusion

By using measured non-Newtonian viscosity and geometry on patient-specific CFD analysis, morphologic differences in hemodynamic features, such as changes in whole blood viscosity and WSS, were observed. Therefore, measured non-Newtonian viscosity might be possibly useful to obtain patient-specific hemodynamic and morphologic result.

Evidence for non-Newtonian behavior of intracranial blood flow from Doppler ultrasonography measurements

Computational fluid dynamics (CFD) studies of intracranial hemodynamics often use Newtonian viscosity model to close the shear rate term in the Navier-Stokes equation. This is based on a commonly accepted hypothesis which state that non-Newtonian effects can be neglected in intracranial blood flow. This study aims to examine the validity of such hypothesis to guide future CFD studies of intracranial hemodynamics. Doppler ultrasonography (DUS) measurements of systolic and diastolic vessel diameter and blood velocity were conducted on 16 subjects (mean age 50.6). The measurements were conducted on the internal carotid (ICA), middle cerebral (MCA), and anterior communicating (AComA) arteries. Systolic and diastolic wall shear stress (WSS) values were calculated via the Hagen-Poiseuille exact solution using Newtonian and three different non-Newtonian models: namely Carreau, power-law and Herschel-Bulkley models. The Weissenberg-Rabinowitsch correction for blood shear-thinning viscosity was applied to the non-Newtonian models. The error percentage between the two sets of models was calculated and discussed. The Newtonian hypothesis was tested statistically and discussed using paired t tests. Significant differences (P < 0.0001) were found between the Newtonian and non-Newtonian WSS in ICA. In MCA and AComA, similar differences were found except in the systole and diastole for the Herschel-Bulkley and power-law models (P = 0.0669, P = 0.7298), respectively. The error between the Newtonian and non-Newtonian models ranged from - 27 to 30% (0.2 to 2.2 Pa). These values could affect the physical interpretation of IA CFD studies. Evidence suggests that the Newtonian assumption may be inappropriate to investigate intracranial hemodynamics. Graphical abstract The WSS estimation error resulting from using the Newtonian assumption compared to three non-Newtonian models for ICA, MCA, and AComA in systole and diastole conditions, based on TCCD measurements of 16 subjects. The error due to the Newtonian assumption ranged from 0.2 to 2.2 Pa (- 27 to 30%). These values could affect the physical interpretation of IA CFD studies.

A comparative study of porous medium CFD models for flow diverter stents: Advantages and shortcomings

In computational fluid dynamics, there is a high interest in modeling flow diverter stents as porous media due to its reduced computational loads. One of the main difficulties of such models is proper parameter setup. Most authors assume flow diverter's wire screen as an isotropic and homogeneous medium, while others proposes anisotropic configurations, yet very little is discussed about the effect of these assumptions on model's accuracy. In this paper, we compare the effect of different models on hemodynamics in relation to their parameters. The fidelity and efficiency of the different models to capture wire screen effect on fluid flow are quantitatively analyzed and compared.

3DRA Reconstruction of Intracranial Aneurysms - How does Voxel Size Influences Morphologic and Hemodynamic Parameters

Three-dimensional shape analysis and imagebased hemodynamic simulations are widely used to assess the individual rupture risk of intracranial aneurysms. However, the quality of those results highly depends on pre-simulative working steps including image reconstruction and segmentation. Within this study, three patient-specific aneurysms were reconstructed using three different voxel sizes (0.1 mm, 0.3 mm, 0.5 mm). Afterwards, 3D segmentations and time-dependent blood flow simulations were carried out to evaluate the impact of the reconstruction size. The results indicate that overall all voxel sizes lead to a qualitatively good agreement with respect to the aneurysm surfaces. However, deviations occur regarding the neck representation as well as the consideration of perforating arteries. Further, morphological differences lead to clear hemodynamic variations, especially for shear force predictions. The findings indicate that depending on the desired analysis, careful reconstruction parameter selection is required. Particularly, for quantitative morphology and blood flow studies, the early step of reconstruction can have a crucial effect on subsequent results.

Software-based simulation for preprocedural assessment of braided stent sizing: a validation study

OBJECTIVE

The authors sought to validate the use of a software-based simulation for preassessment of braided self-expanding stents in the treatment of wide-necked intracranial aneurysms.

METHODS

This was a retrospective, observational, single-center study of 13 unruptured and ruptured intracranial aneurysms treated with braided self-expanding stents. Pre- and postprocedural angiographic studies were analyzed. ANKYRAS software was used to compare the following 3 variables: the manufacturer-given nominal length (NL), software-calculated simulated length (SL), and the actual measured length (ML) of the stent. Appropriate statistical methods were used to draw correlations among the 3 lengths.

RESULTS

In this study, data obtained in 13 patients treated with braided self-expanding stents were analyzed. Data for the 3 lengths were collected for all patients. Error discrepancy was calculated by mean squared error (NL to ML -22.2; SL to ML -6.14, p < 0.05), mean absolute error (NL to ML 3.88; SL to ML -1.84, p < 0.05), and mean error (NL to ML -3.81; SL to ML -1.22, p < 0.05).

CONCLUSIONS

The ML was usually less than the NL given by the manufacturer, indicating significant change in length in most cases. Computational software-based simulation for preassessment of the braided self-expanding stents is a safe and effective way for accurately calculating the change in length to aid in choosing the right-sized stent for optimal placement in complex intracranial vasculature.

Wall Shear Stress and Flow Patterns in Unruptured and Ruptured Anterior Communicating Artery Aneurysms Using Computational Fluid Dynamics

Objective

The goal of this study was to compare several parameters, including wall shear stress (WSS) and flow pattern, between unruptured and ruptured anterior communicating artery (ACoA) aneurysms using patient-specific aneurysm geometry.

Methods

In total, 18 unruptured and 24 ruptured aneurysms were analyzed using computational fluid dynamics (CFD) models. Minimal, average, and maximal wall shear stress were calculated based on CFD simulations. Aneurysm height, ostium diameter, aspect ratio, and area of aneurysm were measured. Aneurysms were classified according to flow complexity (simple or complex) and inflow jet (concentrated or diffused). Statistical analyses were performed to ascertain differences between the aneurysm groups.

Results

Average wall shear stress of the ruptured group was greater than that of the unruptured group (9.42% for aneurysm and 10.38% for ostium). The average area of ruptured aneurysms was 31.22% larger than unruptured aneurysms. Simple flow was observed in 14 of 18 (78%) unruptured aneurysms, while all ruptured aneurysms had complex flow (p<0.001). Ruptured aneurysms were more likely to have a concentrated inflow jet (63%), while unruptured aneurysms predominantly had a diffused inflow jet (83%, p=0.004).

Conclusion

Ruptured aneurysms tended to have a larger geometric size and greater WSS compared to unruptured aneurysms, but the difference was not statistically significant. Flow complexity and inflow jet were significantly different between unruptured and ruptured ACoA aneurysms.

Accounting for residence-time in blood rheology models: do we really need non-Newtonian blood flow modelling in large arteries?

Patient-specific computational fluid dynamics (CFD) is a promising tool that provides highly resolved haemodynamics information. The choice of blood rheology is an assumption in CFD models that has been subject to extensive debate. Blood is known to exhibit shear-thinning behaviour, and non-Newtonian modelling has been recommended for aneurysmal flows. Current non-Newtonian models ignore rouleaux formation, which is the key player in blood's shear-thinning behaviour. Experimental data suggest that red blood cell aggregation and rouleaux formation require notable red blood cell residence-time (RT) in a low shear rate regime. This study proposes a novel hybrid Newtonian and non-Newtonian rheology model where the shear-thinning behaviour is activated in high RT regions based on experimental data. Image-based abdominal aortic and cerebral aneurysm models are considered and highly resolved CFD simulations are performed using a minimally dissipative solver. Lagrangian particle tracking is used to define a backward particle RT measure and detect stagnant regions with increased rouleaux formation likelihood. Our novel RT-based non-Newtonian model shows a significant reduction in shear-thinning effects and provides haemodynamic results qualitatively identical and quantitatively close to the Newtonian model. Our results have important implications in patient-specific CFD modelling and suggest that non-Newtonian models should be revisited in large artery flows.

Effects on Aortoiliac Fluid Dynamics After Endovascular Sealing of Abdominal Aneurysms

OBJECTIVES:

To evaluate the effects on aortoiliac fluid dynamics after the implantation of an endograft based on endovascular aneurysm sealing (EVAS) versus endovascular aneurysm repair (EVAR) strategy.

METHODS:

An adaptive geometrical deformable model was used for aortic lumen segmentation in 8 patients before and after the surgery. Abdominal aneurysms were treated with an endograft based on the EVAS system (Nellix, n = 4) and with a device based on an anatomical fixation technology (n = 4). Pressure, blood velocity, and wall shear stress (WSS) were estimated at different aortic regions using computational fluid dynamics methods. Physiologic inlet/outlet flow values at the abdominal aorta, the celiac trunk, and the mesenteric and the renal arteries were set. Pressure references were set at iliac arteries outlet.

RESULTS:

Maximum aneurysm sizes were similar for both groups in the preoperative scans. The lumen area was lower after EVAR ( P < .05) and EVAS ( P < .01) compared to preoperative aortic lumen sizes. Pressure increase was higher in the proximal abdominal aorta after EVAS compared to EVAR (2.3 ± 0.3 mm Hg vs 0.9 ± 0.3 mm Hg, P < .001). Peak blood velocities inside the endografts were 3-fold higher for EVAS compared to EVAR (54 ± 5 cm/s vs 17 ± 4 cm/s, P < .01). Velocities at the iliac arteries also remained higher for EVAS (38 ± 4 cm/s vs 24 ± 4 cm/s, P < .05). Peak WSS at the iliac arteries remained higher for EVAS compared to EVAR group ( P < .05).

CONCLUSION:

The significant modification of the aortic bifurcation anatomy after EVAS alters aortoiliac fluid dynamics, showing a pressure impact at the renal arteries level and an acceleration of the blood velocity at the iliac region with a concomitant increase in peak WSS.

Virtual stenting of intracranial aneurysms: A pilot study for the prediction of treatment success based on hemodynamic simulations

Endovascular treatment of intracranial aneurysms using flow-diverting devices has revolutionized the treatment of large and complex lesions due to its minimally invasive nature and potential clinical outcomes. However, incomplete or delayed occlusion and persistent intracranial aneurysm growth are still an issue for up to one-third of the patients. We evaluated two patients with intracranial aneurysm located at the internal carotid artery who were treated with flow-diverting devices and had opposite outcomes. Both patients presented with similar aneurysms and were treated with the same device, but after a 1-year follow-up, one case presented with complete occlusion (Case 1) and the other required further treatment (Case 2). To reproduce the interventions, virtual stents were deployed and blood flow simulations were carried out using the respective patient-specific geometries. Afterward, hemodynamic metrics such as aneurysmal inflow reduction, wall shear stresses, oscillatory shear, and inflow concentration indices were quantified. The hemodynamic simulations reveal that for both cases, the neck inflow was clearly reduced due to the therapy (Case 1: 19%, Case 2: 35%). In addition, relevant hemodynamic parameters such as time-averaged wall shear stress (Case 1: 35.6%, Case 2: 57%) and oscillatory shear (Case 1: 33.1%, Case 2: 26.7%) were decreased considerably. However, although stronger relative reductions occurred in the unsuccessful case, the absolute flow values in the successful case were approximately halved. The findings demonstrate that a high relative effect of endovascular devices is not necessarily associated with the desired treatment outcome. Instead, it appears that a successful intracranial aneurysm therapy requires a certain patient-specific inflow threshold.

Effect of Local Coil Density on Blood Flow Stagnation in Densely Coiled Cerebral Aneurysms: A Computational Study Using a Cartesian Grid Method

Aneurysm recurrence is the most critical concern following coil embolization of a cerebral aneurysm. Adequate packing density (PD) and coil uniformity are believed necessary to achieve sufficient flow stagnation, which decreases the risk of aneurysm recurrence. The effect of coil distribution on the extent of flow stagnation, however, especially in cases of dense packing (high PD), has received less attention. Thus, the cause of aneurysm recurrence despite dense packing is still an open question. The primary aim of this study is to evaluate the effect of local coil density on the extent of blood flow stagnation in densely coiled aneurysms. For this purpose, we developed a robust computational framework to determine blood flow using a Cartesian grid method, by which the complex fluid pathways in coiled aneurysms could be flexibly treated using an implicit function. This tool allowed us to conduct blood flow analyses in two patient-specific geometries with 50 coil distribution patterns in each aneurysm at clinically adequate PD. The results demonstrated that dense packing in the aneurysm may not necessarily block completely the inflow into the aneurysm and local flow that formed in the neck region, whose strength was inversely related to this local PD. This finding suggests that local coil density in the neck region still plays an important role in disturbing the remaining local flow, which possibly prevents thrombus formation in a whole aneurysm sac, increasing the risk of aneurysm regrowth and subsequent recurrence.

Geometrically enhanced sensor surfaces for the selective capture of cell-like particles in a laminar flow field

Medical wires inserted into the blood stream of patients offer an attractive perspective to capture rare cells such as circulating tumor cells in vivo. A major challenge in such systems is to achieve an efficient interaction of the desired cells with the sensing surface and avoid those cells that simply flow by the wire without any contact while floating in a laminar flow field at some small distance to the sensor surface. We describe a new strategy to increase the interaction of cells or cell-like particles to such wire-shaped sensor surfaces both from an experimental and a theoretical point of view. For model experiments, we use cell-like particles that are flowing past the profile wire in a blood-like liquid stream. In the fluid dynamics simulations, this sensor is inserted into small capillaries. The influence of geometry and orientation of the wire with respect to the surrounding capillary onto the capture behavior is studied. Parameters, such as wire diameter, profile shape, wire torsion, and orientation of it with respect to the liquid stream, induce in some cases quite strong crossflows. These crossflows enhance the contact probability compared to a straight line wire of the same length by factors of up to about 80. A general model connecting the wire geometry with the crossflow intensity and the particle capture behavior is developed. Particle capture experiments demonstrate that the identified geometric factors can improve the capture of cell-like particles in laminar fluid flows and enhance the performance of such cell sensors.

Variability of hemodynamic parameters using the common viscosity assumption in a computational fluid dynamics analysis of intracranial aneurysms

BACKGROUND

In most simulations of intracranial aneurysm hemodynamics, blood is assumed to be a Newtonian fluid. However, it is a non-Newtonian fluid, and its viscosity profile differs among individuals. Therefore, the common viscosity assumption may not be valid for all patients.

OBJECTIVE

This study aims to test the suitability of the common viscosity assumption.

METHODS

Blood viscosity datasets were obtained from two healthy volunteers. Three simulations were performed for three different-sized aneurysms, two using measured value-based non-Newtonian models and one using a Newtonian model. The parameters proposed to predict an aneurysmal rupture obtained using the non-Newtonian models were compared with those obtained using the Newtonian model.

RESULTS

The largest difference (25%) in the normalized wall shear stress (NWSS) was observed in the smallest aneurysm. Comparing the difference ratio to the NWSS with the Newtonian model between the two Non-Newtonian models, the difference of the ratio was 17.3%.

CONCLUSIONS

Irrespective of the aneurysmal size, computational fluid dynamics simulations with either the common Newtonian or non-Newtonian viscosity assumption could lead to values different from those of the patient-specific viscosity model for hemodynamic parameters such as NWSS.