Newtonian and non-Newtonian blood flow in coiled cerebral aneurysms

A Patient-Specific Three-Dimensional Hemodynamic Model of the Circle of Willis

Circle of Willis (CoW) is one of the most important cerebral arteries in the human body and various attempts have been made to study the hemodynamic of blood flow in this vital part of the brain. In the present study, blood flow in a patient specific CoW is numerically modeled to predict disease-prone regions of the CoW. Medical images and computer aided design software are used to construct a realistic three-dimensional model of the CoW for this particular case. The arteries are considered as elastic conduits and the interactions between arterial walls and the blood flow are taken into account. Mooney-Rivlin hyperelastic model is used to describe the behavior of arterial walls and blood is considered as a non-Newtonian fluid obeying the Carreau model. An available experimental-based pulsatile velocity profile is used at the entrance of the CoW. The finite element-based commercial software, ADINA, is used to solve the governing equations. Blood pressure and velocity and arterial wall shear stress are calculated in different regions of the CoW. A simplified form of the model is also compared with the available published data. Results affirmed that the proposed computational model has the potential to capture the hemodynamic characteristics of the CoW. The computational results can be used to determine disease-prone locations for a given CoW.

Virtual endovascular treatment of intracranial aneurysms: models and uncertainty

Virtual endovascular treatment models (VETMs) have been developed with the view to aid interventional neuroradiologists and neurosurgeons to pre-operatively analyze the comparative efficacy and safety of endovascular treatments for intracranial aneurysms. Based on the current state of VETMs in aneurysm rupture risk stratification and in patient-specific prediction of treatment outcomes, we argue there is a need to go beyond personalized biomechanical flow modeling assuming deterministic parameters and error-free measurements. The mechanobiological effects associated with blood clot formation are important factors in therapeutic decision making and models of post-treatment intra-aneurysmal biology and biochemistry should be linked to the purely hemodynamic models to improve the predictive power of current VETMs. The influence of model and parameter uncertainties associated to each component of a VETM is, where feasible, quantified via a random-effects meta-analysis of the literature. This allows estimating the pooled effect size of these uncertainties on aneurysmal wall shear stress. From such meta-analyses, two main sources of uncertainty emerge where research efforts have so far been limited: (1) vascular wall distensibility, and (2) intra/intersubject systemic flow variations. In the future, we suggest that current deterministic computational simulations need to be extended with strategies for uncertainty mitigation, uncertainty exploration, and sensitivity reduction techniques. WIREs Syst Biol Med 2017, 9:e1385. doi: 10.1002/wsbm.1385 For further resources related to this article, please visit the WIREs website.

Does the DSA reconstruction kernel affect hemodynamic predictions in intracranial aneurysms? An analysis of geometry and blood flow variations

BACKGROUND

Computational fluid dynamics (CFD) blood flow predictions in intracranial aneurysms promise great potential to reveal patient-specific flow structures. Since the workflow from image acquisition to the final result includes various processing steps, quantifications of the individual introduced potential error sources are required.

METHODS

Three-dimensional (3D) reconstruction of the acquired imaging data as input to 3D model generation was evaluated. Six different reconstruction modes for 3D digital subtraction angiography (DSA) acquisitions were applied to eight patient-specific aneurysms. Segmentations were extracted to compare the 3D luminal surfaces. Time-dependent CFD simulations were carried out in all 48 configurations to assess the velocity and wall shear stress (WSS) variability due to the choice of reconstruction kernel.

RESULTS

All kernels yielded good segmentation agreement in the parent artery; deviations of the luminal surface were present at the aneurysm neck (up to 34.18%) and in distal or perforating arteries. Observations included pseudostenoses as well as noisy surfaces, depending on the selected reconstruction kernel. Consequently, the hemodynamic predictions show a mean SD of 11.09% for the aneurysm neck inflow rate, 5.07% for the centerline-based velocity magnitude, and 17.83%/9.53% for the mean/max aneurysmal WSS, respectively. In particular, vessel sections distal to the aneurysms yielded stronger variations of the CFD values.

CONCLUSIONS

The choice of reconstruction kernel for DSA data influences the segmentation result, especially for small arteries. Therefore, if precise morphology measurements or blood flow descriptions are desired, a specific reconstruction setting is required. Furthermore, research groups should be encouraged to denominate the kernel types used in future hemodynamic studies.

Does Arterial Flow Rate Affect the Assessment of Flow-Diverter Stent Performance?

BACKGROUND AND PURPOSE

Our aim was to assess the performance of flow-diverter stents. The pre- and end-of-treatment angiographies are commonly compared. However, the arterial flow rate may change between acquisitions; therefore, a better understanding of its influence on the local intra-aneurysmal hemodynamics before and after flow-diverter stent use is required.

MATERIALS AND METHODS

Twenty-five image-based aneurysm models extracted from 3D rotational angiograms were conditioned for computational fluid dynamics simulations. Pulsatile simulations were performed at different arterial flow rates, covering a wide possible range of physiologic flows among 1-5 mL/s. The effect of flow-diverter stents on intra-aneurysmal hemodynamics was numerically simulated with a porous medium model. Spatiotemporal-averaged intra-aneurysmal flow velocity and flow rate were calculated for each case to quantify the hemodynamics after treatment. The short-term flow-diverter stent performance was characterized by the relative velocity reduction inside the aneurysm.

RESULTS

Spatiotemporal-averaged intra-aneurysmal flow velocity before and after flow-diverter stent use is linearly proportional to the mean arterial flow rate (minimum R² > 0.983 of the linear regression models for untreated and stented models). Relative velocity reduction asymptotically decreases with increasing mean arterial flow rate. When the most probable range of arterial flow rate was considered (3-5 mL/s), instead of the wide possible flow range, the mean SD of relative velocity reduction was reduced from 3.6% to 0.48%.

CONCLUSIONS

Both intra-aneurysmal aneurysm velocity and flow-diverter stent performance depend on the arterial flow rate. The performance could be considered independent of the arterial flow rates within the most probable range of physiologic flows.

Blood flow in intracranial aneurysms treated with Pipeline embolization devices: computational simulation and verification with Doppler ultrasonography on phantom models

Purpose:

The aim of this study was to validate a computational fluid dynamics (CFD) simulation of flow-diverter treatment through Doppler ultrasonography measurements in patient-specific models of intracranial bifurcation and side-wall aneurysms.

Methods:

Computational and physical models of patient-specific bifurcation and sidewall aneurysms were constructed from computed tomography angiography with use of stereolithography, a three-dimensional printing technology. Flow dynamics parameters before and after flow-diverter treatment were measured with pulse-wave and color Doppler ultrasonography, and then compared with CFD simulations.

Results:

CFD simulations showed drastic flow reduction after flow-diverter treatment in both aneurysms. The mean volume flow rate decreased by 90% and 85% for the bifurcation aneurysm and the side-wall aneurysm, respectively. Velocity contour plots from computer simulations before and after flow diversion closely resembled the patterns obtained by color Doppler ultrasonography.

Conclusion:

The CFD estimation of flow reduction in aneurysms treated with a flow-diverting stent was verified by Doppler ultrasonography in patient-specific phantom models of bifurcation and side-wall aneurysms. The combination of CFD and ultrasonography may constitute a feasible and reliable technique in studying the treatment of intracranial aneurysms with flow-diverting stents.

URINE FLOW DYNAMICS THROUGH THE URETHRA IN PATIENTS WITH BLADDER OUTLET OBSTRUCTION

Voiding dysfunction is a common disease among the elderly male population. However, few studies have elucidated the flow dynamics of the disease mechanism. We attempted to clarify the fluid dynamics of urine in the prostatic urethra (PU) in patients with bladder outlet obstruction (BOO), a common disease associated with voiding dysfunction. A model of normal PU and three phase models following the BOO progression were prepared using 3D CAD software tracing cystourethroscopic images. In fluid dynamics computation, the normal detrusor pressure was applied to the bladder side in each model. Results of fluid simulation were compared for flow trajectory line, fluid energy dynamics using Bernoulli's principle and pressure loss coefficient. In the BOO model, a large vortex was observed adjacent to the bladder outlet (BO) with increased hydraulic energy loss during the passage through the tract. In contrast, opening of the BO eliminated the vortex with reduction of hydraulic energy loss. These results corresponded with clinical data measured by catheterized pressure sensors inserted into the urethra, and those of urethral diameter measured by urethrography. Fluid dynamics simulation using CAD models can clarify the mechanism of voiding dysfunction in a less invasive and more acceptable way in patients with BOO.