PIV-measured versus CFD-predicted flow dynamics in anatomically realistic cerebral aneurysm models

Design, manufacturing, and multi-modal imaging of stereolithography 3D printed flexible intracranial aneurysm phantoms

BACKGROUND

Physical vascular phantoms are instrumental in studying intracranial aneurysms and testing relevant imaging tools and training systems to provide improved clinical care. Current vascular phantom production methods have major limitations in capturing the biophysical and morphological characteristics of intracranial aneurysms with good fidelity and multi-modal imaging capacity. With stereolithography (SLA) 3D printing technology becoming more accessible, newer flexible and transparent printing materials with higher precision controls open the door for improving the efficiency and quality of producing anthropomorphic vascular phantoms but have rarely been explored for the application.

PURPOSE

This technical note intends to report the feasibility of using SLA 3D printing technology to manufacture flexible intracranial aneurysm phantoms with similar scales to the real anatomy, as well as their capacity for multi-modal flow imaging and analysis, including ultrasound flow imaging, high-speed filming, and particle image velocimetry analysis.

METHODS

We designed and 3D-printed two intracranial aneurysm phantoms with an SLA 3D printer using Formlabs Elastic 50A resin. By using a micropump to introduce cyclical flows in the phantoms, we first employed conventional Doppler and vector flow ultrasonography to observe and measure different fluidic properties. Then, a high-speed camera was used to record particles flowing within the phantom, and we further conducted a particle image velocimetry analysis, including the distribution of mean 2D velocity vectors, average velocity magnitudes, and the mean vorticity fields in the phantom for the high-speed imaging data.

RESULTS

We successfully 3D-printed flexible intracranial aneurysm phantoms with similar dimensions to the real anatomy. Additionally, we validated the phantoms' ability to allow visualization, measurement, and analysis of flow dynamics based on both real-time ultrasound and optical imaging.

CONCLUSIONS

Our proof-of-concept study illustrates that SLA 3D printing using commercial elastic resins can significantly contribute towards facilitating the fabrication of flexible intracranial aneurysms phantoms for training, research, and preoperative planning.

On flow fluctuations in ruptured and unruptured intracranial aneurysms: resolved numerical study

Flow fluctuations have emerged as a promising hemodynamic metric for understanding of hemodynamics in intracranial aneurysms. Several investigations have reported flow instabilities using numerical tools. In this study, the occurrence of flow fluctuations is investigated using either Newtonian or non-Newtonian fluid models in five patient-specific intracranial aneurysms using high-resolution lattice Boltzmann simulation methods. Flow instabilities are quantified by computing power spectral density, proper orthogonal decomposition, and fluctuating kinetic energy of velocity fluctuations. Our simulations reveal substantial flow instabilities in two of the ruptured aneurysms, where the pulsatile inflow through the neck leads to hydrodynamic instability, particularly around the rupture position, throughout the entire cardiac cycle. In other monitoring points, the flow instability is primarily observed during the deceleration phase; typically, the fluctuations begin just after peak systole, gradually decay, and the flow returns to its original, laminar pulsatile state during diastole. Additionally, we assess the rheological impact on flow dynamics. The disparity between Newtonian and non-Newtonian outcomes remains minimal in unruptured aneurysms, with less than a 5% difference in key metrics. However, in ruptured cases, adopting a non-Newtonian model yields a substantial increase in the fluctuations within the aneurysm sac, with up to a 30% higher fluctuating kinetic energy compared to the Newtonian model. The study highlights the importance of using appropriate high-resolution simulations and non-Newtonian models to capture flow fluctuation characteristics that may be critical for assessing aneurysm rupture risk.

Literature Survey for In-Vivo Reynolds and Womersley Numbers of Various Arteries and Implications for Compliant In-Vitro Modelling

PURPOSE

In-vitro modelling can be used to investigate haemodynamics of arterial geometry and stent implants. However, in-vitro model fidelity relies on precise matching of in-vivo conditions. In pulsatile flow, velocity distribution and wall shear stress depend on compliance, and the Reynolds and Womersley numbers. However, matching such values may lead to unachievable tolerances in phantom fabrication.

METHODS

Published Reynolds and Womersley numbers for 14 major arteries in the human body were determined via a literature search. Preference was given to in-vivo publications but in-vitro and in-silico values were presented when in-vivo values were not found. Subsequently ascending aorta and carotid artery case studies were presented to highlight the limitations dynamic matching would apply to phantom fabrication.

RESULTS

Seven studies reported the in-vivo Reynolds and Womersley numbers for the aorta and two for the carotid artery. However, only one study each reported in-vivo numbers for the remaining ten arteries. No in-vivo data could be found for the femoral, superior mesenteric and renal arteries. Thus, information derived in-vitro and in-silico were provided instead. The ascending aorta and carotid artery models required scaling to 1.5× and 3× life-scale, respectively, to achieve dimensional tolerance restrictions. Modelling the ascending aorta with the comparatively high viscosity water/glycerine solution will lead to high pump power demands. However, all the working fluids considered could be dynamically matched with low pump demand for the carotid model.

CONCLUSION

This paper compiles available human haemodynamic information, and highlights the paucity of information for some arteries. It also provides a method for optimal in-vitro experimental configuration.

Turbulent blood flow in a cerebral artery with an aneurysm

Unruptured intracranial aneurysms are common in the general population, and many uncertainties remain when predicting rupture risks and treatment outcomes. One of the cutting-edge tools used to investigate this condition is computational fluid dynamics (CFD). However, CFD is not yet mature enough to guide the clinical management of this disease. In addition, recent studies have reported significant flow instabilities when refined numerical methods are used. Questions remain as to how to properly simulate and evaluate this flow, and whether these instabilities are really turbulence. The purpose of the present study is to evaluate the impact of the simulation setup on the results and investigate the occurrence of turbulence in a cerebral artery with an aneurysm. For this purpose, direct numerical simulations were performed with up to 200 cardiac cycles and with data sampling rates of up to 100,000 times per cardiac cycle. Through phase-averaging or triple decomposition, the contributions of turbulence and of laminar pulsatile waves to the velocity, pressure and wall shear stress fluctuations were distinguished. For example, the commonly used oscillatory shear index was found to be closely related to the laminar waves introduced at the inlet, rather than turbulence. The turbulence energy cascade was evaluated through energy spectrum estimates, revealing that, despite the low flow rates and Reynolds number, the flow is turbulent near the aneurysm. Phase-averaging was shown to be an approach that can help researchers better understand this flow, although the results are highly dependent on simulation setup and post-processing choices.

Computational Fluid Dynamic Simulations of Cerebral Aneurysms

Computational fluid dynamics (CFD) simulations have been introduced to enable individualized risk prognosis for patients with unruptured cerebral aneurysms. The present contribution provides an overview of the biomechanical and physiological principles of aneurysm formation and rupture. It describes the computational steps of the CFD and the evaluated parameters. The clinical value of CFD is then discussed based on a recent literature review. Finally, we discuss current methodological limitations and possible future developments to overcome the actual drawbacks of CFD.

Optimization of 4D Flow MRI Spatial and Temporal Resolution for Examining Complex Hemodynamics in the Carotid Artery Bifurcation

BACKGROUND

Three-dimensional, ECG-gated, time-resolved, three-directional, velocity-encoded phase-contrast MRI (4D flow MRI) has been applied extensively to measure blood velocity in great vessels but has been much less used in diseased carotid arteries. Carotid artery webs (CaW) are non-inflammatory intraluminal shelf-like projections into the internal carotid artery (ICA) bulb that are associated with complex flow and cryptogenic stroke.

PURPOSE

Optimize 4D flow MRI for measuring the velocity field of complex flow in the carotid artery bifurcation model that contains a CaW.

METHODS

A 3D printed phantom model created from computed tomography angiography (CTA) of a subject with CaW was placed in a pulsatile flow loop within the MRI scanner. 4D Flow MRI images of the phantom were acquired with five different spatial resolutions (0.50-2.00  mm3) and four different temporal resolutions (23-96 ms) and compared to a computational fluid dynamics (CFD) solution of the flow field as a reference. We examined four planes perpendicular to the vessel centerline, one in the common carotid artery (CCA) and three in the internal carotid artery (ICA) where complex flow was expected. At these four planes pixel-by-pixel velocity values, flow, and time average wall shear stress (TAWSS) were compared between 4D flow MRI and CFD.

HYPOTHESIS

An optimized 4D flow MRI protocol will provide a good correlation with CFD velocity and TAWSS values in areas of complex flow within a clinically feasible scan time (~ 10 min).

RESULTS

Spatial resolution affected the velocity values, time average flow, and TAWSS measurements. Qualitatively, a spatial resolution of 0.50  mm3 resulted in higher noise, while a lower spatial resolution of 1.50-2.00  mm3 did not adequately resolve the velocity profile. Isotropic spatial resolutions of 0.50-1.00  mm3 showed no significant difference in total flow compared to CFD. Pixel-by-pixel velocity correlation coefficients between 4D flow MRI and CFD were > 0.75 for 0.50-1.00  mm3 but were < 0.5 for 1.50 and 2.00  mm3. Regional TAWSS values determined from 4D flow MRI were generally lower than CFD and decreased at lower spatial resolutions (larger pixel sizes). TAWSS differences between 4D flow and CFD were not statistically significant at spatial resolutions of 0.50-1.00  mm3 but were different at 1.50 and 2.00 mm3. Differences in temporal resolution only affected the flow values when temporal resolution was > 48.4 ms; temporal resolution did not affect TAWSS values.

CONCLUSION

A spatial resolution of 0.74-1.00  mm3 and a temporal resolution of 23-48 ms (1-2 k-space segments) provides a 4D flow MRI protocol capable of imaging velocity and TAWSS in regions of complex flow within the carotid bifurcation at a clinically acceptable scan time.

Hemodynamic and fluid flow analysis of a cerebral aneurysm: a CFD simulation

In this study, we investigate the hemodynamics parameters and their impact on the aneurysm rupture. The simulations are performed on an ideal (benchmark) and realistic model for the intracranial aneurysm that appears at the anterior communicating artery. The realistic geometry was reconstructed from patient-specific cerebral arteries. The computational fluid dynamics simulations are utilized to investigate the hemodynamic parameters such as flow recirculation, wall shear stress, and wall pressure. The boundary conditions are measured from the patient using ultrasonography. The solution of the governing equations is obtained by using the ANSYS-FLUENT 19.2 package. The CFD results indicate that the flow recirculation appears in the aneurysms zone. The effect of the flow recirculation on the bulge hemodynamics wall parameters is discussed to identify the rupture zone.

Method for Fabricating Transparent Patient-Specific Vocal Tract Replicas

INTRODUCTION

Transparent, patient-specific vocal tract replicas are helpful in research and educational endeavors but challenging to procure. An accessible method for fabricating these models, improving on previously suggested processes, would make them more widely available.

METHOD

Detailed instructions for fabricating a transparent, patient-specific vocal tract model were addressed. The broad steps were (1) digitally reconstructing (patient-specific) vocal tract geometry, (2) producing a vocal tract mold (using methods such as three-dimensional [3D] printing), and (3) casting transparent material (such as silicone) around the vocal tract mold and removing the mold. The cavities remaining within the cast represented the exact geometry of the vocal tract.

DISCUSSION

A combination of 3D printing and silicone casting can produce useful vocal tract replicas. Several simple changes to previous methods can improve consistency and reduce the labor and cost of production. Limitations and potential modifications to expand the applications of this method are discussed.

Fabrication of Compliant and Transparent Hollow Cerebral Vascular Phantoms for In Vitro Studies Using 3D Printing and Spin-Dip Coating

Endovascular surgery through flow diverters and coils is increasingly used for the minimally invasive treatment of intracranial aneurysms. To study the effectiveness of these devices, in vitro tests are performed in which synthetic vascular phantoms are typically used to reproduce in vivo conditions. In this paper, we propose a manufacturing process to obtain compliant and transparent hollow vessel replicas to assess the mechanical behaviour of endovascular devices and perform flow measurements. The vessel models were obtained in three main steps. First, a mould was 3D-printed in a water-soluble material; two techniques, fusion deposition modelling and stereolithography, were compared for this purpose. Then, the mould was covered with a thin layer of silicone through spin-dip coating, and finally, when the silicone layer solidified, it was dissolved in a hot water bath. The final models were tested in terms of the quality of the final results, the mechanical properties of the silicone, thickness uniformity, and transparency properties. The proposed approach makes it possible to produce models of different sizes and complexity whose transparency and mechanical properties are suitable for in vitro experiments. Its applicability is demonstrated through idealised and patient-specific cases.

Developing an in vitro validated 3D in silico internal carotid artery sidewall aneurysm model

Introduction: Direct quantification of hemodynamic factors applied to a cerebral aneurysm (CA) remains inaccessible due to the lack of technologies to measure the flow field within an aneurysm precisely. This study aimed to develop an in vitro validated 3D in silico patient-specific internal carotid artery sidewall aneurysm (ICASA) model which can be used to investigate hemodynamic factors on the CA pathophysiology. Methods: The validated ICASA model was developed by quantifying and comparing the flow field using particle image velocimetry (PIV) measurements and computational fluid dynamics (CFD) simulations. Specifically, the flow field characteristics, i.e., blood flowrates, normalized velocity profiles, flow streamlines, and vortex locations, have been compared at representative time instants in a cardiac pulsatile period in two designated regions of the ICASA model, respectively. One region is in the internal carotid artery (ICA) inlet close to the aneurysm sac, the other is across the middle of the aneurysmal sac. Results and Discussion: The results indicated that the developed computational fluid dynamics model presents good agreements with the results from the parallel particle image velocimetry and flowrate measurements, with relative differences smaller than 0.33% in volumetric flow rate in the ICA and relative errors smaller than 9.52% in averaged velocities in the complex aneurysmal sac. However, small differences between CFD and PIV in the near wall regions were observed due to the factors of slight differences in the 3D printed model, light reflection and refraction near arterial walls, and flow waveform uncertainties. The validated model not only can be further employed to investigate hemodynamic factors on the cerebral aneurysm pathophysiology statistically, but also provides a typical model and guidance for other professionals to evaluate the hemodynamic effects on cerebral aneurysms.

Smartphone-based particle image velocimetry for cardiovascular flows applications: A focus on coronary arteries

An experimental set-up is presented for the in vitro characterization of the fluid dynamics in personalized phantoms of healthy and stenosed coronary arteries. The proposed set-up was fine-tuned with the aim of obtaining a compact, flexible, low-cost test-bench for biomedical applications. Technically, velocity vector fields were measured adopting a so-called smart-PIV approach, consisting of a smartphone camera and a low-power continuous laser (30 mW). Experiments were conducted in realistic healthy and stenosed 3D-printed phantoms of left anterior descending coronary artery reconstructed from angiographic images. Time resolved image acquisition was made possible by the combination of the image acquisition frame rate of last generation commercial smartphones and the flow regimes characterizing coronary hemodynamics (velocities in the order of 10 cm/s). Different flow regimes (Reynolds numbers ranging from 20 to 200) were analyzed. The smart-PIV approach was able to provide both qualitative flow visualizations and quantitative results. A comparison between smart-PIV and conventional PIV (i.e., the gold-standard experimental technique for bioflows characterization) measurements showed a good agreement in the measured velocity vector fields for both the healthy and the stenosed coronary phantoms. Displacement errors and uncertainties, estimated by applying the particle disparity method, confirmed the soundness of the proposed smart-PIV approach, as their values fell within the same range for both smart and conventional PIV measured data (≈5% for the normalized estimated displacement error and below 1.2 pixels for displacement uncertainty). In conclusion, smart-PIV represents an easy-to-implement, low-cost methodology for obtaining an adequately robust experimental characterization of cardiovascular flows. The proposed approach, to be intended as a proof of concept, candidates to become an easy-to-handle test bench suitable for use also outside of research labs, e.g., for educational or industrial purposes, or as first-line investigation to direct and guide subsequent conventional PIV measurements.

Hemodynamic characteristics in a cerebral aneurysm model using non-Newtonian blood analogues

This study aims to develop an experimentally validated computational fluid dynamics (CFD) model to estimate hemodynamic characteristics in cerebral aneurysms (CAs) using non-Newtonian blood analogues. Blood viscosities varying with shear rates were measured under four temperatures first, which serves as the reference for the generation of blood analogues. Using the blood analogue, particle image velocimetry (PIV) measurements were conducted to quantify flow characteristics in a CA model. Then, using the identical blood properties in the experiment, CFD simulations were executed to quantify the flow patterns, which were used to compare with the PIV counterpart. Additionally, hemodynamic characteristics in the simplified Newtonian and non-Newtonian models were quantified and compared using the experimentally validated CFD model. Results showed the proposed non-Newtonian viscosity model can predict blood shear-thinning properties accurately under varying temperatures and shear rates. Another developed viscosity model based on the blood analogue can well represent blood rheological properties. The comparisons in flow characteristics show good agreements between PIV and CFD, demonstrating the developed CFD model is qualified to investigate hemodynamic factors within CAs. Furthermore, results show the differences of absolute values were insignificant between Newtonian and non-Newtonian fluids in the distributions of wall shear stress (WSS) and oscillatory shear index (OSI) on arterial walls. However, not only does the simplified Newtonian model underestimate WSS and OSI in most regions of the aneurysmal sac, but it also makes mistakes in identifying the high OSI regions on the sac surface, which may mislead the hemodynamic assessment on the pathophysiology of CAs.

Effects of Pulsatile Flow Rate and Shunt Ratio in Bifurcated Distal Arteries on Hemodynamic Characteristics Involved in Two Patient-Specific Internal Carotid Artery Sidewall Aneurysms: A Numerical Study

The pulsatile flow rate (PFR) in the cerebral artery system and shunt ratios in bifurcated arteries are two patient-specific parameters that may affect the hemodynamic characteristics in the pathobiology of cerebral aneurysms, which needs to be identified comprehensively. Accordingly, a systematic study was employed to study the effects of pulsatile flow rate (i.e., PFR−I, PFR−II, and PFR−III) and shunt ratio (i.e., 75:25 and 64:36) in bifurcated distal arteries, and transient cardiac pulsatile waveform on hemodynamic patterns in two internal carotid artery sidewall aneurysm models using computational fluid dynamics (CFD) modeling. Numerical results indicate that larger PFRs can cause higher wall shear stress (WSS) in some local regions of the aneurysmal dome that may increase the probability of small/secondary aneurysm generation than under smaller PFRs. The low WSS and relatively high oscillatory shear index (OSI) could appear under a smaller PFR, increasing the potential risk of aneurysmal sac growth and rupture. However, the variances in PFRs and bifurcated shunt ratios have rare impacts on the time-average pressure (TAP) distributions on the aneurysmal sac, although a higher PFR can contribute more to the pressure increase in the ICASA−1 dome due to the relatively stronger impingement by the redirected bloodstream than in ICASA−2. CFD simulations also show that the variances of shunt ratios in bifurcated distal arteries have rare impacts on the hemodynamic characteristics in the sacs, mainly because the bifurcated location is not close enough to the sac in present models. Furthermore, it has been found that the vortex location plays a major role in the temporal and spatial distribution of the WSS on the luminal wall, varying significantly with the cardiac period.

Review of the Development of Hemodynamic Modeling Techniques to Capture Flow Behavior in Arteries Affected by Aneurysm, Atherosclerosis, and Stenting

Cardiovascular diseases (CVDs) are the leading cause of death in the developed world. CVD can include atherosclerosis, aneurysm, dissection, or occlusion of the main arteries. Many CVDs are caused by unhealthy hemodynamics. Some CVDs can be treated with the implantation of stents and stent grafts. Investigations have been carried out to understand the effects of stents and stent grafts have on arteries and the hemodynamic changes post-treatment. Numerous studies on stent hemodynamics have been carried out using computational fluid dynamics (CFD) which has yielded significant insight into the effect of stent mesh design on near-wall blood flow and improving hemodynamics. Particle image velocimetry (PIV) has also been used to capture behavior of fluids that mimic physiological hemodynamics. However, PIV studies have largely been restricted to unstented models or intra-aneurysmal flow rather than peri or distal stent flow behaviors. PIV has been used both as a standalone measurement method and as a comparison to validate the CFD studies. This article reviews the successes and limitations of CFD and PIV-based modeling methods used to investigate the hemodynamic effects of stents. The review includes an overview of physiology and relevant mechanics of arteries as well as consideration of boundary conditions and the working fluids used to simulate blood for each modeling method along with the benefits and limitations introduced.

A multi-modality approach for enhancing 4D flow magnetic resonance imaging via sparse representation

This work evaluates and applies a multi-modality approach to enhance blood flow measurements and haemodynamic analysis with phase-contrast magnetic resonance imaging (4D flow MRI) in cerebral aneurysms (CAs). Using a library of high-resolution velocity fields from patient-specific computational fluid dynamic simulations and in vitro particle tracking velocimetry measurements, the flow field of 4D flow MRI data is reconstructed as the sparse representation of the library. The method was evaluated with synthetic 4D flow MRI data in two CAs. The reconstruction enhanced the spatial resolution and velocity accuracy of the synthetic MRI data, leading to reliable pressure and wall shear stress (WSS) evaluation. The method was applied on in vivo 4D flow MRI data acquired in the same CAs. The reconstruction increased the velocity and WSS by 6-13% and 39-61%, respectively, suggesting that the accuracy of these quantities was improved since the raw MRI data underestimated the velocity and WSS by 10-20% and 40-50%, respectively. The computed pressure fields from the reconstructed data were consistent with the observed flow structures. The results suggest that using the sparse representation flow reconstruction with in vivo 4D flow MRI enhances blood flow measurement and haemodynamic analysis.

Properties and Applications of PDMS for Biomedical Engineering: A Review

Polydimethylsiloxane (PDMS) is an elastomer with excellent optical, electrical and mechanical properties, which makes it well-suited for several engineering applications. Due to its biocompatibility, PDMS is widely used for biomedical purposes. This widespread use has also led to the massification of the soft-lithography technique, introduced for facilitating the rapid prototyping of micro and nanostructures using elastomeric materials, most notably PDMS. This technique has allowed advances in microfluidic, electronic and biomedical fields. In this review, an overview of the properties of PDMS and some of its commonly used treatments, aiming at the suitability to those fields' needs, are presented. Applications such as microchips in the biomedical field, replication of cardiovascular flow and medical implants are also reviewed.

Hemodynamic Comparison of Stent-Grafts for the Treatment of Aortoiliac Occlusive Disease

Purpose:

To compare the flow patterns and hemodynamics of the AFX stent-graft and the covered endovascular reconstruction of aortic bifurcation (CERAB) configuration using laser particle image velocimetry (PIV) experiments.

Materials and Methods:

Two anatomically realistic aortoiliac phantoms were constructed using polydimethylsiloxane polymer. An AFX stent-graft with a transparent cover made with a new method was inserted into one phantom. A CERAB configuration using Atrium’s Avanta V12 with transparent covers made with a previously established method was inserted into the other phantom, both modified stent-grafts were suitable for laser PIV, enabling visualization of the flow fields and quantification of time average wall shear stress (TAWSS), oscillatory shear index (OSI), and relative residence time (RRT).

Results:

Disturbed flow was observed at the bifurcation region of the AFX, especially at the end systolic velocity (ESV) time-point where recirculation was noticeable due to vortical flow. In contrast, predominantly unidirectional flow was observed at the CERAB bifurcation. These observations were confirmed by the quantified hemodynamic results from PIV analysis where mean TAWSS of 0.078 Pa (range: 0.009–0.242 Pa) was significantly lower in AFX as compared with 0.229 Pa (range: 0.013–0.906 Pa) for CERAB (p<0.001). Mean OSI of 0.318 (range: 0.123–0.496) in AFX was significantly higher than 0.252 (range: 0.055–0.472) in CERAB (p<0.001). Likewise, mean RRT of 180 Pa−1 (range: 9–3603 Pa−1) in AFX was also significantly higher than 88 Pa−1 (range: 2–840 Pa−1) in CERAB (p=0.0086).

Conclusion:

In this in vitro study, the flow pattern of a modified AFX stent-graft was found to be more disturbed especially at the end systolic phase, its hemodynamic outcomes less desirable than CERAB configuration.

Clinical Relevance:

While the AFX stent-graft has an advantage over the CERAB configuration in eliminating radial mismatch, and maintaining the anatomical bifurcation for future endovascular intervention, this in vitro study revealed that the associated lower TAWSS, higher OSI and RRT may predispose to thrombosis and are, thus, less desirable as compared to a CERAB configuration. Further investigation is warranted to confirm whether these findings translate into the clinical setting.

Higher oscillatory shear index is related to aneurysm recanalization after coil embolization in posterior communicating artery aneurysms

BACKGROUND

The recurrence rate of posterior communicating artery (Pcom) aneurysms after endovascular treatment (EVT) is higher than that for aneurysms located in other sites. However, it is still unclear what mechanisms are responsible for the recanalization of cerebral aneurysms. In this investigation, we compared hemodynamic factors related with recanalization of Pcom aneurysms treated by endoluminal coiling using computational fluid dynamics (CFD) with high-resolution three-dimensional digital subtraction angiography images.

METHODS

Twenty patients were enrolled. A double-sinogram acquisition was performed with and without contrast injection after coil embolization to get true blood vessel lumen by relatively complementing the first sinogram with the second. Adaptive Cartesian meshing was performed to produce vascular wall objects for CFD simulation. The boundary condition for inlet (ICA) was set for dynamic velocity according to the cardiac cycle (0.8 s). Hemodynamic parameters were recorded at two specific points (branching point of Pcom and residual sac). The peak pressure, peak WSS, and oscillatory shear index (OSI) were recorded and analyzed.

RESULTS

The median age was 61.0 years, and 18 patients (90%) were female. During a median follow-up of 12 months, seven (35%) treated aneurysms showed recanalization. The median aneurysm volume was significantly higher, and aneurysm height and neck sizes were significantly longer in the recanalization group than those in the stable group. At the branching point of the Pcom, the peak pressure, peak WSS, or OSI did not significantly differ between the two groups. The only statistically significant hemodynamic parameter related with recanalization was the OSI at the aneurysm point. Multivariate logistic regression showed that with an increase of 0.01 OSI at the aneurysm point, the odds ratio for the aneurysm recanalization was 1.19.

CONCLUSIONS

A higher OSI is related with recanalization after coil embolization for a Pcom aneurysm.

FSI modeling on the effect of artery-aneurysm thickness and coil embolization in patient cases

BACKGROUND AND OBJECTIVE

The attainment of a methodology to simulate the hemodynamic in patient-specific cerebral vessels with aneurysms is still a challenge. The novelty of this work is focused on the effect of coil embolization in a realistic cerebral aneurysm, according to the vessel wall thickness and aneurysm thickness, through transient FSI simulations.

METHODS

The quality of the mesh for simulations was checked with a specific mesh convergence study; and the numerical methodology was validated using numerical research data of the literature. The model was implemented in ANSYS® software. The total deformation and equivalent stress evolution in the studied cases, before and after coil embolization, were compared. More than 20 different models were employed due to different arterial wall thickness and aneurysm wall thickness combinations.

RESULTS

The obtained results have showed that deformation and stress values are highly influenced with the sac thickness. The thinner sac aneurysm thickness is, the greater deformation and stress are. The results after coil embolization process have highlighted that considering typical values of arterial wall thickness and aneurysm thickness 0.3 mm and 0.15 mm respectively, a deformation reduction around 50% and a stress reduction around 70% can be achieved.

CONCLUSIONS

The proposed methodology is a step forward in the personalized medicine, quantifying the aneurysm rupture risk reduction, and helping the medical team in the preoperative planning, or to deciding the optimal treatment.

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.

Hemodynamic Comparison of Stent-Grafts for the Treatment of Aortoiliac Occlusive Disease

PURPOSE

To compare the flow patterns and hemodynamics of the AFX stent-graft and the covered endovascular reconstruction of aortic bifurcation (CERAB) configuration using laser particle image velocimetry (PIV) experiments.

MATERIALS AND METHODS

Two anatomically realistic aortoiliac phantoms were constructed using polydimethylsiloxane polymer. An AFX stent-graft with a transparent cover made with a new method was inserted into one phantom. A CERAB configuration using Atrium's Avanta V12 with transparent covers made with a previously established method was inserted into the other phantom, both modified stent-grafts were suitable for laser PIV, enabling visualization of the flow fields and quantification of time average wall shear stress (TAWSS), oscillatory shear index (OSI), and relative residence time (RRT).

RESULTS

Disturbed flow was observed at the bifurcation region of the AFX, especially at the end systolic velocity (ESV) time-point where recirculation was noticeable due to vortical flow. In contrast, predominantly unidirectional flow was observed at the CERAB bifurcation. These observations were confirmed by the quantified hemodynamic results from PIV analysis where mean TAWSS of 0.078 Pa (range: 0.009-0.242 Pa) was significantly lower in AFX as compared with 0.229 Pa (range: 0.013-0.906 Pa) for CERAB (p<0.001). Mean OSI of 0.318 (range: 0.123-0.496) in AFX was significantly higher than 0.252 (range: 0.055-0.472) in CERAB (p<0.001). Likewise, mean RRT of 180 Pa^-1 (range: 9-3603 Pa^-1) in AFX was also significantly higher than 88 Pa^-1 (range: 2-840 Pa^-1) in CERAB (p=0.0086).

CONCLUSION

In this in vitro study, the flow pattern of a modified AFX stent-graft was found to be more disturbed especially at the end systolic phase, its hemodynamic outcomes less desirable than CERAB configuration.

CLINICAL RELEVANCE

While the AFX stent-graft has an advantage over the CERAB configuration in eliminating radial mismatch, and maintaining the anatomical bifurcation for future endovascular intervention, this in vitro study revealed that the associated lower TAWSS, higher OSI and RRT may predispose to thrombosis and are, thus, less desirable as compared to a CERAB configuration. Further investigation is warranted to confirm whether these findings translate into the clinical setting.

Comparison of Hemodynamic Visualization in Cerebral Arteries: Can Magnetic Resonance Imaging Replace Computational Fluid Dynamics?

A multimodality approach was applied using four-dimensional flow magnetic resonance imaging (4D flow MRI), time-of-flight magnetic resonance angiography (TOF-MRA) signal intensity gradient (SIG), and computational fluid dynamics (CFD) to investigate the 3D blood flow characteristics and wall shear stress (WSS) of the cerebral arteries. TOF-MRA and 4D flow MRI were performed on the major cerebral arteries in 16 healthy volunteers (mean age 34.7 ± 7.6 years). The flow rate measured with 4D flow MRI in the internal carotid artery, middle cerebral artery, and anterior cerebral artery were 3.8, 2.5, and 1.2 mL/s, respectively. The 3D blood flow pattern obtained through CFD and 4D flow MRI on the cerebral arteries showed reasonable consensus. CFD delivered much greater resolution than 4D flow MRI. TOF-MRA SIG and CFD WSS of the major cerebral arteries showed reasonable consensus with the locations where the WSS was relatively high. However, the visualizations were very different between TOF-MRA SIG and CFD WSS at the internal carotid artery bifurcations, the anterior cerebral arteries, and the anterior communicating arteries. 4D flow MRI, TOF-MRA SIG, and CFD are complementary methods that can provide additional insight into the hemodynamics of the human cerebral artery.

Central moments multiple relaxation time LBM for hemodynamic simulations in intracranial aneurysms: An in-vitro validation study using PIV and PC-MRI

The lattice Boltzmann method (LBM) has recently emerged as an efficient alternative to classical Navier-Stokes solvers. This is particularly true for hemodynamics in complex geometries. However, in its most basic formulation, i.e. with the so-called single relaxation time (SRT) collision operator, it has been observed to have a limited stability domain in the Courant/Fourier space, strongly constraining the minimum time-step and grid size. The development of improved collision models such as the multiple relaxation time (MRT) operator in central moments space has tremendously widened the stability domain, while allowing to overcome a number of other well-documented artifacts, therefore opening the door for simulations over a wider range of grid and time-step sizes. The present work focuses on implementing and validating a specific collision operator, the central Hermite moments multiple relaxation time model with the full expansion of the equilibrium distribution function, to simulate blood flows in intracranial aneurysms. The study further proceeds with a validation of the numerical model through different test-cases and against experimental measurements obtained via stereoscopic particle image velocimetry (PIV) and phase-contrast magnetic resonance imaging (PC-MRI). For a patient-specific aneurysm both PIV and PC-MRI agree fairly well with the simulation. Finally, low-resolution simulations were shown to be able to capture blood flow information with sufficient accuracy, as demonstrated through both qualitative and quantitative analysis of the flow field while leading to strongly reduced computation times. For instance in the case of the patient-specific configuration, increasing the grid-size by a factor of two led to a reduction of computation time by a factor of 14 with very good similarity indices still ranging from 0.83 to 0.88.

In Vitro Assessment of Flow Variability in an Intracranial Aneurysm Model Using 4D Flow MRI and Tomographic PIV

Aneurysm rupture has been suggested to be related to aneurysm geometry, morphology, and complex flow activity; therefore, understanding aneurysm-specific hemodynamics is crucial. 4D Flow MRI has been shown to be a feasible tool for assessing hemodynamics in intracranial aneurysms with high spatial resolution. However, it requires averaging over multiple heartbeats and cannot account for cycle-to-cycle hemodynamics variations. This study aimed to assess cycle-to-cycle flow dynamics variations in a patient-specific intracranial aneurysm model using tomographic particle image velocimetry (tomo-PIV) at a high image rate under pulsatile flow conditions. Time-resolved and time-averaged velocity flow fields within the aneurysm sac and estimations of wall shear stress (WSS) were compared with those from 4D Flow MRI. A one-way ANOVA showed a significant difference between cardiac cycles (p value < 0.0001); however, differences were not significant after PIV temporal and spatial resolution was matched to that of MRI (p value 0.9727). This comparison showed the spatial resolution to be the main contributor to assess cycle-to-cycle variability. Furthermore, the comparison with 4D Flow MRI between velocity components, streamlines, and estimated WSS showed good qualitative and quantitative agreement. This study showed the feasibility of patient-specific in-vitro experiments using tomo-PIV to assess 4D Flow MRI with high repeatability in the measurements.

A novel roller pump for physiological flow

Having physiological correct flow waveforms is a key feature for experimental studies of blood flow, especially in the process of developing and testing a new medical device such as stent, mechanical heart valve, or any implantable medical device that involves circulation of blood through the device. It is also a critical part of a perfusion system for cardiopulmonary bypass and extracorporeal membrane oxygenation procedures. This study investigated the feasibility of a novel roller pump for use in experimental flow phantoms. Flow rates of carotid flow profile measured directly with the ultrasonic flow meter matched well with the reference flow rates programmed into the machine with similarity index of 0.97 and measured versus programmed flow rates at specific time‐points of peak systolic velocity (PSV): 0.894 vs 0.880, end systolic velocity (ESV): 0.333 vs 0.319, and peak diastolic velocity (PDV): 0.514 vs 0.520 L/min. Flow rates derived from video analysis of the pump motion for carotid, suprarenal, and infrarenal flows also matched well with references with similarity indices of 0.99, 0.99, and 0.96, respectively. Measured flow rates (mean/standard deviation) at PSV, ESV, and PDV time‐points for carotid: 0.883/0.016 vs 0.880, 0.342/0.007 vs 0.319, and 0.485/0.009 vs 0.520; suprarenal: 3.497/0.014 vs 3.500, 0.004/0.003 vs 0, and 1.656/0.073 vs 1.453; infrarenal: 4.179/0.024 vs 4.250, −1.147/0.015 vs −1.213, and 0.339/0.017 vs 0.391 L/min, respectively. The novel roller pump is suitable for benchtop testing of physiological flow.

A novel roller pump for physiological flow

Having physiological correct flow waveforms is a key feature for experimental studies of blood flow, especially in the process of developing and testing a new medical device such as stent, mechanical heart valve, or any implantable medical device that involves circulation of blood through the device. It is also a critical part of a perfusion system for cardiopulmonary bypass and extracorporeal membrane oxygenation procedures. This study investigated the feasibility of a novel roller pump for use in experimental flow phantoms. Flow rates of carotid flow profile measured directly with the ultrasonic flow meter matched well with the reference flow rates programmed into the machine with similarity index of 0.97 and measured versus programmed flow rates at specific time‐points of peak systolic velocity (PSV): 0.894 vs 0.880, end systolic velocity (ESV): 0.333 vs 0.319, and peak diastolic velocity (PDV): 0.514 vs 0.520 L/min. Flow rates derived from video analysis of the pump motion for carotid, suprarenal, and infrarenal flows also matched well with references with similarity indices of 0.99, 0.99, and 0.96, respectively. Measured flow rates (mean/standard deviation) at PSV, ESV, and PDV time‐points for carotid: 0.883/0.016 vs 0.880, 0.342/0.007 vs 0.319, and 0.485/0.009 vs 0.520; suprarenal: 3.497/0.014 vs 3.500, 0.004/0.003 vs 0, and 1.656/0.073 vs 1.453; infrarenal: 4.179/0.024 vs 4.250, −1.147/0.015 vs −1.213, and 0.339/0.017 vs 0.391 L/min, respectively. The novel roller pump is suitable for benchtop testing of physiological flow.

Hemodynamics of Cerebral Aneurysms: Connecting Medical Imaging and Biomechanical Analysis

In the last two decades, numerous studies have conducted patient-specific computations of blood flow dynamics in cerebral aneurysms and reported correlations between various hemodynamic metrics and aneurysmal disease progression or treatment outcomes. Nevertheless, intra-aneurysmal flow analysis has not been adopted in current clinical practice, and hemodynamic factors usually are not considered in clinical decision making. This review presents the state of the art in cerebral aneurysm imaging and image-based modeling, discussing the advantages and limitations of each approach and focusing on the translational value of hemodynamic analysis. Combining imaging and modeling data obtained from different flow modalities can improve the accuracy and fidelity of resulting velocity fields and flow-derived factors that are thought to affect aneurysmal disease progression. It is expected that predictive models utilizing hemodynamic factors in combination with patient medical history and morphological data will outperform current risk scores and treatment guidelines. Possible future directions include novel approaches enabling data assimilation and multimodality analysis of cerebral aneurysm hemodynamics.

Evaluation of a Desktop 3D Printed Rigid Refractive-Indexed-Matched Flow Phantom for PIV Measurements on Cerebral Aneurysms

PURPOSE

Fabrication of a suitable flow model or phantom is critical to the study of biomedical fluid dynamics using optical flow visualization and measurement methods. The main difficulties arise from the optical properties of the model material, accuracy of the geometry and ease of fabrication.

METHODS

Conventionally an investment casting method has been used, but recently advancements in additive manufacturing techniques such as 3D printing have allowed the flow model to be printed directly with minimal post-processing steps. This study presents results of an investigation into the feasibility of fabrication of such models suitable for particle image velocimetry (PIV) using a common 3D printing Stereolithography process and photopolymer resin.

RESULTS

An idealised geometry of a cerebral aneurysm was printed to demonstrate its applicability for PIV experimentation. The material was shown to have a refractive index of 1.51, which can be refractive matched with a mixture of de-ionised water with ammonium thiocyanate (NH₄SCN). The images were of a quality that after applying common PIV pre-processing techniques and a PIV cross-correlation algorithm, the results produced were consistent within the aneurysm when compared to previous studies.

CONCLUSIONS

This study presents an alternative low-cost option for 3D printing of a flow phantom suitable for flow visualization simulations. The use of 3D printed flow phantoms reduces the complexity, time and effort required compared to conventional investment casting methods by removing the necessity of a multi-part process required with investment casting techniques.

Development of Phantom Material that Resembles Compression Properties of Human Brain Tissue for Training Models

There is a need to quantify and reproduce the mechanical behavior of brain tissue for a variety of applications from designing proper training models for surgeons to enabling research on the effectiveness of personal protective gear, such as football helmets. The mechanical response of several candidate phantom materials, including hydrogels and emulsions, was characterized and compared to porcine brain tissue under similar strains and strain rates. Some candidate materials were selected since their compositions were similar to brain tissue, such as emulsions that mimic the high content of lipids. Others, like silicone, were included since these are currently used as phantom materials. The mechanical response of the emulsion was closer to that of the native porcine brain tissue than the other candidates. The emulsions, created by addition of oil to a hydrogel, were able to withstand compressive strain greater than 40%. The addition of lipids in the emulsions also prevented the syneresis typically seen with hydrogel materials. This allowed the emulsion material to undergo freeze-thaw cycles with no significant change in their mechanical properties.

Multi-modality cerebral aneurysm haemodynamic analysis: in vivo 4D flow MRI, in vitro volumetric particle velocimetry and in silico computational fluid dynamics

Typical approaches to patient-specific haemodynamic studies of cerebral aneurysms use image-based computational fluid dynamics (CFD) and seek to statistically correlate parameters such as wall shear stress (WSS) and oscillatory shear index (OSI) to risk of growth and rupture. However, such studies have reported contradictory results, emphasizing the need for in-depth multi-modality haemodynamic metric evaluation. In this work, we used in vivo 4D flow MRI data to inform in vitro particle velocimetry and CFD modalities in two patient-specific cerebral aneurysm models (basilar tip and internal carotid artery). Pulsatile volumetric particle velocimetry experiments were conducted, and the particle images were processed using Shake-the-Box, a particle tracking method. Distributions of normalized WSS and relative residence time were shown to be highly yet inconsistently affected by minor flow field and spatial resolution variations across modalities, and specific relationships among these should be explored in future work. Conversely, OSI, a non-dimensional parameter, was shown to be more robust to the varying assumptions, limitations and spatial resolutions of each subject and modality. These results suggest a need for further multi-modality analysis as well as development of non-dimensional haemodynamic parameters and correlation of such metrics to aneurysm risk of growth and rupture.

Fabrication of Low-Cost Patient-Specific Vascular Models for Particle Image Velocimetry

PURPOSE

Particle image velocimetry (PIV), an in vitro experimentation technique that optically measures velocity components to analyze fluid velocity fields, has become increasingly popular to study flow dynamics in various vascular territories. However, it can be difficult and expensive to create patient-specific clear models for PIV due to the importance of refractive index matching of the model and the fluid. We aim to implement and test the use of poly-vinyl alcohol (PVA) in a lost-core casting technique to create low-cost, patient-specific models for PIV.

METHODS

Anonymized patient vascular anatomies were segmented and processed in Mimics/3Matic to create patient-specific cores from 3D digital subtraction angiographies. The cores were 3D-printed with PVA and post-processed with a 80:20 water:glue mixture to smooth the surface. Two silicones, Sylgard 184 and Solaris, were used to encapsulate the model and the PVA core was dissolved using warm water. Computed tomography scans were used to evaluate geometric accuracy using circumferences and surface differences in the model.

RESULTS

Mean geometric differences in circumference along the inlet centerline and the mean surface difference in the aneurysm between the final Silicone Model and the desired STL Print geometry were statistically insignificant (0.6 mm, 95% CI [- 1.4, 2.8] and 0.3 mm 95% CI [- 0.1, 0.7], respectively). Particle illumination within each model was successful. The cost of one 10 cm × 10 cm × 5 cm model was $69.

CONCLUSION

This technique was successful to implement and test the use of PVA in a lost-core casting technique to create low-cost, patient-specific in vitro models for PIV experimentation.

Comparison of intracranial aneurysm flow quantification techniques: standard PIV vs stereoscopic PIV vs tomographic PIV vs phase-contrast MRI vs CFD

Image-based hemodynamic simulations to assess the rupture risk or improve the treatment planning of intracranial aneurysms have become popular recently. However, due to strong modeling assumptions and limitations, the acceptance of numerical approaches remains limited. Therefore, validation using experimental methods is mandatory.

In this study, a unique compilation of four in-vitro flow measurements (three particle image velocimetry approaches using a standard (PIV), stereoscopic (sPIV), and tomographic (tPIV) setup, as well as a phase-contrast magnetic resonance imaging (PC-MRI) measurement) were compared with a computational fluid dynamics (CFD) simulation. This was carried out in a patient-specific silicone phantom model of an internal carotid artery aneurysm under steady flow conditions. To evaluate differences between each technique, a similarity index (SI) with respect to the velocity vectors and the average velocity magnitude differences among all involved modalities were computed.

The qualitative comparison reveals that all techniques are able to provide a reasonable description of the global flow structures. High quantitative agreement in terms of SI and velocity magnitude differences was found between all PIV methods and CFD. However, quantitative differences were observed between PC-MRI and the other techniques. Deeper analysis revealed that the limited resolution of the PC-MRI technique is a major contributor to the experienced differences and leads to a systematic underestimation of overall velocity magnitude levels inside the vessel. This confirms the necessity of using highly resolving flow measurement techniques, such as PIV, in an in-vitro environment to individually verify the validity of the numerically obtained hemodynamic results.

Computational Fluid Dynamics Modeling of the Human Pulmonary Arteries with Experimental Validation

Pulmonary hypertension (PH) is a chronic progressive disease characterized by elevated pulmonary arterial pressure, caused by an increase in pulmonary arterial impedance. Computational fluid dynamics (CFD) can be used to identify metrics representative of the stage of PH disease. However, experimental validation of CFD models is often not pursued due to the geometric complexity of the model or uncertainties in the reproduction of the required flow conditions. The goal of this work is to validate experimentally a CFD model of a pulmonary artery phantom using a particle image velocimetry (PIV) technique. Rapid prototyping was used for the construction of the patient-specific pulmonary geometry, derived from chest computed tomography angiography images. CFD simulations were performed with the pulmonary model with a Reynolds number matching those of the experiments. Flow rates, the velocity field, and shear stress distributions obtained with the CFD simulations were compared to their counterparts from the PIV flow visualization experiments. Computationally predicted flow rates were within 1% of the experimental measurements for three of the four branches of the CFD model. The mean velocities in four transversal planes of study were within 5.9 to 13.1% of the experimental mean velocities. Shear stresses were qualitatively similar between the two methods with some discrepancies in the regions of high velocity gradients. The fluid flow differences between the CFD model and the PIV phantom are attributed to experimental inaccuracies and the relative compliance of the phantom. This comparative analysis yielded valuable information on the accuracy of CFD predicted hemodynamics in pulmonary circulation models.

Identification of vortex structures in a cohort of 204 intracranial aneurysms

An intracranial aneurysm (IA) is a cerebrovascular pathology that can lead to death or disability if ruptured. Abnormal wall shear stress (WSS) has been associated with IA growth and rupture, but little is known about the underlying flow physics related to rupture-prone IAs. Previous studies, based on analysis of a few aneurysms or partial views of three-dimensional vortex structures, suggest that rupture is associated with complex vortical flow inside IAs. To further elucidate the relevance of vortical flow in aneurysm pathophysiology, we studied 204 patient IAs (56 ruptured and 148 unruptured). Using objective quantities to identify three-dimensional vortex structures, we investigated the characteristics associated with aneurysm rupture and if these features correlate with previously proposed WSS and morphological characteristics indicative of IA rupture. Based on the Q-criterion definition of a vortex, we quantified the degree of the aneurysmal region occupied by vortex structures using the volume vortex fraction (vVF) and the surface vortex fraction (sVF). Computational fluid dynamics simulations showed that the sVF, but not the vVF, discriminated ruptured from unruptured aneurysms. Furthermore, we found that the near-wall vortex structures co-localized with regions of inflow jet breakdown, and significantly correlated to previously proposed haemodynamic and morphologic characteristics of ruptured IAs.

Hemodynamics in a giant intracranial aneurysm characterized by in vitro 4D flow MRI

Experimental and computational data suggest that hemodynamics play a critical role in the development, growth, and rupture of cerebral aneurysms. The flow structure, especially in aneurysms with a large sac, is highly complex and three-dimensional. Therefore, volumetric and time-resolved measurements of the flow properties are crucial to fully characterize the hemodynamics. In this study, phase-contrast Magnetic Resonance Imaging is used to assess the fluid dynamics inside a 3D-printed replica of a giant intracranial aneurysm, whose hemodynamics was previously simulated by multiple research groups. The physiological inflow waveform is imposed in a flow circuit with realistic cardiovascular impedance. Measurements are acquired with sub-millimeter spatial resolution for 16 time steps over a cardiac cycle, allowing for the detailed reconstruction of the flow evolution. Moreover, the three-dimensional and time-resolved pressure distribution is calculated from the velocity field by integrating the fluid dynamics equations, and is validated against differential pressure measurements using precision transducers. The flow structure is characterized by vortical motions that persist within the aneurysm sac for most of the cardiac cycle. All the main flow statistics including velocity, vorticity, pressure, and wall shear stress suggest that the flow pattern is dictated by the aneurysm morphology and is largely independent of the pulsatility of the inflow, at least for the flow regimes investigated here. Comparisons are carried out with previous computational simulations that used the same geometry and inflow conditions, both in terms of cycle-averaged and systolic quantities.

Methodology for Computational Fluid Dynamic Validation for Medical Use: Application to Intracranial Aneurysm

Computational fluid dynamics (CFD) is a promising tool to aid in clinical diagnoses of cardiovascular diseases. However, it uses assumptions that simplify the complexities of the real cardiovascular flow. Due to high-stakes in the clinical setting, it is critical to calculate the effect of these assumptions in the CFD simulation results. However, existing CFD validation approaches do not quantify error in the simulation results due to the CFD solver's modeling assumptions. Instead, they directly compare CFD simulation results against validation data. Thus, to quantify the accuracy of a CFD solver, we developed a validation methodology that calculates the CFD model error (arising from modeling assumptions). Our methodology identifies independent error sources in CFD and validation experiments, and calculates the model error by parsing out other sources of error inherent in simulation and experiments. To demonstrate the method, we simulated the flow field of a patient-specific intracranial aneurysm (IA) in the commercial CFD software star-ccm+. Particle image velocimetry (PIV) provided validation datasets for the flow field on two orthogonal planes. The average model error in the star-ccm+ solver was 5.63 ± 5.49% along the intersecting validation line of the orthogonal planes. Furthermore, we demonstrated that our validation method is superior to existing validation approaches by applying three representative existing validation techniques to our CFD and experimental dataset, and comparing the validation results. Our validation methodology offers a streamlined workflow to extract the "true" accuracy of a CFD solver.

Computational Fluid Dynamics and Additive Manufacturing to Diagnose and Treat Cardiovascular Disease

Noninvasive engineering models are now being used for diagnosing and planning the treatment of cardiovascular disease. Techniques in computational modeling and additive manufacturing have matured concurrently, and results from simulations can inform and enable the design and optimization of therapeutic devices and treatment strategies. The emerging synergy between large-scale simulations and 3D printing is having a two-fold benefit: first, 3D printing can be used to validate the complex simulations, and second, the flow models can be used to improve treatment planning for cardiovascular disease. In this review, we summarize and discuss recent methods and findings for leveraging advances in both additive manufacturing and patient-specific computational modeling, with an emphasis on new directions in these fields and remaining open questions.

Computational and experimental assessment of influences of hemodynamic shear stress on carotid plaque

Background

Studies have identified hemodynamic shear stress as an important determinant of endothelial function and atherosclerosis. In this study, we assess the influences of hemodynamic shear stress on carotid plaques.

Methods

Carotid stenosis phantoms with three severity (30, 50, 70%) were made from 10% polyvinyl alcohol (PVA) cryogel. The phantoms were placed in a pulsatile flow loop with the same systolic/diastolic phase (35/65) and inlet flow rate (16 L/h). Ultrasonic particle imaging velocimetry (Echo PIV) and computational fluid dynamics (CFD) were used to calculate the velocity profile and shear stress distribution in the carotid stenosis phantoms. Inlet/outlet boundary conditions used in CFD were extracted from Echo PIV experiments to make sure that the results were comparable.

Results

Echo PIV and CFD results showed that velocity was largest in 70% than those in 30 and 50% at peak systole. Echo PIV results indicated that shear stress was larger in the upper wall and the surface of plaque than in the center of vessel. CFD results demonstrated that wall shear stress in the upstream was larger than in downstream of plaque. There was no significant difference in average velocity obtained by CFD and Echo PIV in 30% (p = 0.25). Velocities measured by CFD in 50% (93.01 cm/s) and in 70% (115.07 cm/s) were larger than those by Echo PIV in 50% (60.26 ± 5.36 cm/s) and in 70% (89.11 ± 7.21 cm/s).

Conclusions

The results suggested that Echo PIV and CFD could obtain hemodynamic shear stress on carotid plaques. Higher WSS occurred in narrower arteries, and the shoulder of plaque bore higher WSS than in bottom part.

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

BACKGROUND

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Transitional hemodynamics in intracranial aneurysms - Comparative velocity investigations with high resolution lattice Boltzmann simulations, normal resolution ANSYS simulations, and MR imaging

PURPOSE

Blood flow in intracranial aneurysms has, until recently, been considered to be disturbed but still laminar. Recent high resolution computational studies have demonstrated, in some situations, however, that the flow may exhibit high frequency fluctuations that resemble weakly turbulent or transitional flow. Due to numerous assumptions required for simplification in computational fluid dynamics (CFD) studies, the occurrence of these events, in vivo, remains unsettled. The detection of these fluctuations in aneurysmal blood flow, i.e., hemodynamics by CFD, poses additional challenges as such phenomena cannot be captured in clinical data acquisition with magnetic resonance (MR) due to inadequate temporal and spatial resolutions. The authors' purpose was to address this issue by comparing results from highly resolved simulations, conventional resolution laminar simulations, and MR measurements, identify the differences, and identify their causes.

METHODS

Two aneurysms in the basilar artery, one with disturbed yet laminar flow and the other with transitional flow, were chosen. One set of highly resolved direct numerical simulations using the lattice Boltzmann method (LBM) and another with adequate resolutions under laminar flow assumption were conducted using a commercially available ANSYS Fluent solver. The velocity fields obtained from simulation results were qualitatively and statistically compared against each other and with MR acquisition.

RESULTS

Results from LBM, ANSYS Fluent, and MR agree well qualitatively and quantitatively for one of the aneurysms with laminar flow in which fluctuations were <80 Hz. The comparisons for the second aneurysm with high fluctuations of > ∼ 600 Hz showed vivid differences between LBM, ANSYS Fluent, and magnetic resonance imaging. After ensemble averaging and down-sampling to coarser space and time scales, these differences became minimal.

CONCLUSIONS

A combination of MR derived data and CFD can be helpful in estimating the hemodynamic environment of intracranial aneurysms. Adequately resolved CFD would suffice gross assessment of hemodynamics, potentially in a clinical setting, and highly resolved CFD could be helpful in a detailed and retrospective understanding of the physiological mechanisms.

Flow Instability Detected by High-Resolution Computational Fluid Dynamics in Fifty-Six Middle Cerebral Artery Aneurysms

Recent high-resolution computational fluid dynamics (CFD) studies have detected persistent flow instability in intracranial aneurysms (IAs) that was not observed in previous in silico studies. These flow fluctuations have shown incidental association with rupture in a small aneurysm dataset. The aims of this study are to explore the capabilities and limitations of a commercial cfd solver in capturing such velocity fluctuations, whether fluctuation kinetic energy (fKE) as a marker to quantify such instability could be a potential parameter to predict aneurysm rupture, and what geometric parameters might be associated with such fluctuations. First, we confirmed that the second-order discretization schemes and high spatial and temporal resolutions are required to capture these aneurysmal flow fluctuations. Next, we analyzed 56 patient-specific middle cerebral artery (MCA) aneurysms (12 ruptured) by transient, high-resolution CFD simulations with a cycle-averaged, constant inflow boundary condition. Finally, to explore the mechanism by which such flow instabilities might arise, we investigated correlations between fKE and several aneurysm geometrical parameters. Our results show that flow instabilities were present in 8 of 56 MCA aneurysms, all of which were unruptured bifurcation aneurysms. Statistical analysis revealed that fKE could not differentiate ruptured from unruptured aneurysms. Thus, our study does not lend support to these flow instabilities (based on a cycle-averaged constant inflow as opposed to peak velocity) being a marker for rupture. We found a positive correlation between fKE and aneurysm size as well as size ratio. This suggests that the intrinsic flow instability may be associated with the breakdown of an inflow jet penetrating the aneurysm space.