Comprehensive validation of computational fluid dynamics simulations of in-vivo blood flow in patient-specific cerebral aneurysms

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.

Identification of intra-individual variation in intracranial arterial flow by MRI and the effect on computed hemodynamic descriptors

OBJECTIVES

To determine the intra-individual flow variation in serially acquired studies, and the influence of this variation on subsequent hemodynamic simulations using the inlet flow as a boundary condition. Author: Kindly check and confirm whether the corresponding authors are correctly identified.Confirmed.

MATERIALS AND METHODS

This prospective study included 51 patients (37 females and 14 males) with unruptured intracranial aneurysms who have received more than three times follow-up of 2D phase-contrast MR. The flow and velocity parameters were extracted to calculate the reproducibility and variation. Patient-specific computational fluid dynamics simulations were performed using the measured flows.

RESULTS

Intraclass correlation coefficients for mean and maximum velocity and flow parameters ranged from 0.77 to 0.90. A 10% CV of mean flow was identified. Variations of 10% in inlet flow resulted in hemodynamic changes including 41.41% of peak systolic wall shear stress; 39.13% of end-diastolic wall shear stress; 2.79% of low shear area at peak systole; 2.12% of low shear area at end diastole: 47.57% of time-averaged wall shear stress; and 0.17% of oscillatory shear index.

CONCLUSION

This study identified 10% of intra-individual mean flow variation on phase-contrast MR. Intra-individual flow variation resulted in a non-negligible variation in wall shear stress, but relatively small variation in low shear area in hemodynamic calculations.

Computational feasibility of calculating the steady-state blood flow rate through the vasculature of the entire human body

The human body contains approximately 20 billion individual blood vessels that deliver nutrients and oxygen to tissues. While blood flow is a well-developed field of research, no previous studies have calculated the blood flow rates through more than 5 million connected vessels. The goal of this study was to test if it is computationally feasible to calculate the blood flow rates through a vasculature equal in size to that of the human body. We designed and implemented a two-step algorithm to calculate the blood flow rates using principles of steady-state fluid dynamics. Steady-state fluid dynamics is an accurate approximation for the microvascular and venous structures in the human body. To determine the computational feasibility, we measured and evaluated the execution time, scalability, and memory usage to quantify the computational requirements. We demonstrated that it is computationally feasible to calculate the blood flow rate through 17 billion vessels in 6.5 hours using 256 compute nodes. The computational modeling of blood flow rate in entire organisms may find application in research on drug delivery, treatment of cancer metastases, and modulation of physiological performance.

Role of Morphological and Hemodynamic Factors in Predicting Intracranial Aneurysm Rupture: A Review

Intracranial aneurysms (IAs) carry the risk of rupture, which will lead to subarachnoid hemorrhage, which has a high mortality and morbidity risk. However, the treatment of IA's carries mortality and morbidity risks too. There are well-known risk factors for the rupture of IAs like age, size, and site. However, choosing patients with unruptured IAs for treatment is still a big challenge. This review article aimed to find out the relationship between morphological and hemodynamic characters of IAs with their rupture and incorporate these factors with well-known factors to yield an accurate module for predicting the rupture of IAs and decision-making in the treatment of unruptured IAs.

We searched in PubMed and Medline databases by using the following keywords: IAs, subarachnoid hemorrhage, and risk of rupture, morphology, and hemodynamic “mesh.” A total of 19 studies with 7269 patients and 9167 IAs, of which 1701 had ruptured, were reviewed thoroughly. Some modules like population, hypertension, age, size, earlier subarachnoid hemorrhage, and site (PHASES) score that involve well-known risk factors can be used to assess the risk of rupture of IAs. However, decision making for treating unruptured IA needs more detailed and more accurate modules. Studying morphological and hemodynamic factors and incorporation of them with well-known risk factors to yield a more comprehensive module will be very helpful in treating unruptured IA. Among morphological factors, aspect ratio (AR), size ratio (SR), aneurysm height, and bottle-neck factor showed significant effects on the growth and rupture of IA. Besides, wall shear stress (WSS), oscillatory shear index (OSI), and low wall shear stress area (LSA) as hemodynamic factors could have a substantial impact on the formation, shape, growth, and rupture of unruptured IA.

Pressure and Flow Rate Changes During Contrast Injections in Cerebral Angiography: Correlation to Reflux Length

Cerebral angiography involves the antegrade injection of contrast media through a catheter into the vasculature to visualize the region of interest under X-ray imaging. Depending on the injection and blood flow parameters, the bolus of contrast can propagate in the upstream direction and proximal to the catheter tip, at which point contrast is said to have refluxed. In this in vitro study, we investigate the relationship of fundamental hemodynamic variables to this phenomenon. Contrast injections were carried out under steady and pulsatile flow using various vessel diameters, catheter sizes, working fluid flow rates, and injection rates. The distance from the catheter tip to the proximal edge of the contrast bolus, called reflux length, was measured on the angiograms; the relation of this reflux length to different hemodynamic parameters was evaluated. Results show that contrast reflux occurs when the pressure distal to the catheter tip increases to be greater than the pressure proximal to the catheter tip. The ratio of this pressure difference to the baseline flow rate, called reflux resistance here, was linearly correlated to the normalized reflux length (reflux length/vessel diameter). Further, the ratio of blood flow to contrast fluid momentums, called the Craya-Curtet number, was correlated to the normalized reflux length via a sigmoid function. A sigmoid function was also found to be representative of the relationship between the ratio of the Reynolds numbers of blood flow to contrast and the normalized reflux length. As described by previous reports, catheter based contrast injections cause substantial increases in local flow and pressure. Contrast reflux should generally be avoided during standard antegrade angiography. Our study shows two specific correlations between contrast reflux length and baseline and intra-injection parameters that have not been published previously. Further studies need to be conducted to fully characterize the phenomena and to extract reliable indicators of clinical utility. Parameters relevant to cerebral angiography are studied here, but the essential principles are applicable to all angiographic procedures involving antegrade catheter injections.

Comparison of phase-contrast MR and flow simulations for the study of CSF dynamics in the cervical spine

Background Investigators use phase-contrast magnetic resonance (PC-MR) and computational fluid dynamics (CFD) to assess cerebrospinal fluid dynamics. We compared qualitative and quantitative results from the two methods. Methods Four volunteers were imaged with a heavily T2-weighted volume gradient echo scan of the brain and cervical spine at 3T and with PC-MR. Velocities were calculated from PC-MR for each phase in the cardiac cycle. Mean pressure gradients in the PC-MR acquisition through the cardiac cycle were calculated with the Navier-Stokes equations. Volumetric MR images of the brain and upper spine were segmented and converted to meshes. Models of the subarachnoid space were created from volume images with the Vascular Modeling Toolkit. CFD simulations were performed with a previously verified flow solver. The flow patterns, velocities and pressures were compared in PC-MR and CFD flow images. Results PC-MR images consistently revealed more inhomogeneous flow patterns than CFD, especially in the anterolateral subarachnoid space where spinal nerve roots are located. On average, peak systolic and diastolic velocities in PC-MR exceeded those in CFD by 31% and 41%, respectively. On average, systolic and diastolic pressure gradients calculated from PC-MR exceeded those of CFD by 11% and 39%, respectively. Conclusions PC-MR shows local flow disturbances that are not evident in typical CFD. The velocities and pressure gradients calculated from PC-MR are systematically larger than those calculated from CFD.

Computational Fluid Dynamics modeling of contrast transport in basilar aneurysms following flow-altering surgeries

In vivo measurement of blood velocity fields and flow descriptors remains challenging due to image artifacts and limited resolution of current imaging methods; however, in vivo imaging data can be used to inform and validate patient-specific computational fluid dynamics (CFD) models. Image-based CFD can be particularly useful for planning surgical interventions in complicated cases such as fusiform aneurysms of the basilar artery, where it is crucial to alter pathological hemodynamics while preserving flow to the distal vasculature. In this study, patient-specific CFD modeling was conducted for two basilar aneurysm patients considered for surgical treatment. In addition to velocity fields, transport of contrast agent was simulated for the preoperative and postoperative conditions using two approaches. The transport of a virtual contrast passively following the flow streamlines was simulated to predict post-surgical flow regions prone to thrombus deposition. In addition, the transport of a mixture of blood with an iodine-based contrast agent was modeled to compare and verify the CFD results with X-ray angiograms. The CFD-predicted patterns of contrast flow were qualitatively compared to in vivo X-ray angiograms acquired before and after the intervention. The results suggest that the mixture modeling approach, accounting for the flow rates and properties of the contrast injection, is in better agreement with the X-ray angiography data. The virtual contrast modeling assessed the residence time based on flow patterns unaffected by the injection procedure, which makes the virtual contrast modeling approach better suited for prediction of thrombus deposition, which is not limited to the peri-procedural state.

Noninvasive hemodynamic assessment, treatment outcome prediction and follow-up of aortic coarctation from MR imaging

PURPOSE

Coarctation of the aorta (CoA) is a congenital heart disease characterized by an abnormal narrowing of the proximal descending aorta. Severity of this pathology is quantified by the blood pressure drop (△P) across the stenotic coarctation lesion. In order to evaluate the physiological significance of the preoperative coarctation and to assess the postoperative results, the hemodynamic analysis is routinely performed by measuring the △P across the coarctation site via invasive cardiac catheterization. The focus of this work is to present an alternative, noninvasive measurement of blood pressure drop △P through the introduction of a fast, image-based workflow for personalized computational modeling of the CoA hemodynamics.

METHODS

The authors propose an end-to-end system comprised of shape and computational models, their personalization setup using MR imaging, and a fast, noninvasive method based on computational fluid dynamics (CFD) to estimate the pre- and postoperative hemodynamics for coarctation patients. A virtual treatment method is investigated to assess the predictive power of our approach.

RESULTS

Automatic thoracic aorta segmentation was applied on a population of 212 3D MR volumes, with mean symmetric point-to-mesh error of 3.00 ± 1.58 mm and average computation time of 8 s. Through quantitative evaluation of 6 CoA patients, good agreement between computed blood pressure drop and catheter measurements is shown: average differences are 2.38 ± 0.82 mm Hg (pre-), 1.10 ± 0.63 mm Hg (postoperative), and 4.99 ± 3.00 mm Hg (virtual stenting), respectively.

CONCLUSIONS

The complete workflow is realized in a fast, mostly-automated system that is integrable in the clinical setting. To the best of our knowledge, this is the first time that three different settings (preoperative--severity assessment, poststenting--follow-up, and virtual stenting--treatment outcome prediction) of CoA are investigated on multiple subjects. We believe that in future-given wider clinical validation-our noninvasive in-silico method could replace invasive pressure catheterization for CoA.

In vitro validation of endovascular Doppler-derived flow rates in models of the cerebral circulation

This study presents validation of endovascular Doppler velocimetry-based volumetric flow rate measurements conducted in a pulsatile flow loop simulating conditions in both the internal carotid and basilar artery. In vitro models of cerebral vessels, each containing an aneurysm, were fabricated from patient anatomies extracted from 3D rotational angiography. Flow velocity measurements were collected with three different experimental techniques: an endovascular Doppler wire, Particle Image Velocimetry, and a time-resolved ultrasonic flow meter. Womersley's theory of pulsatile flow in a cylindrical vessel was used to compute time-resolved volumetric flow rates from the endovascular Doppler velocity. The volumetric flow rates computed from the Doppler measurements were compared to those from the Particle Image Velocimetry profile measurements, and the direct measurements from the ultrasonic flow meter. The study establishes confidence intervals for any systematic or random errors associated with the wire-derived flow rates as benchmarked to the other two modalities. There is an approximately 10% random error in the Doppler-derived peak and time-averaged flow rates. There is a measurable uniform bias, about 15% too low, in the time-averaged Doppler-derived flow rates. There is also a small proportional bias in the peak systolic Doppler-derived flow rates. Potential sources of error are also discussed.

Validation of computational fluid dynamics methods with anatomically exact, 3D printed MRI phantoms and 4D pcMRI

A new concept of rapid 3D prototyping was implemented using cost-effective 3D printing for creating anatomically correct replica of cerebral aneurysms. With a dedicated flow loop set-up in a full body human MRI scanner, flow measurements were performed using 4D phase contrast magnetic resonance imaging to visualize and quantify intra-aneurysmal flow patterns. Ultrashort TE sequences were employed to obtain high-resolution 3D image data to visualize the lumen inside the plastic replica. In-vitro results were compared with retrospectively obtained in-vivo data and results from computational fluid dynamics simulations (CFD). Rapid prototyping of anatomically realistic 3D models may have future impact in treatment planning, design of image acquisition methods for MRI and angiographic systems and for the design and testing of advanced image post-processing technologies.

Accuracy of computational cerebral aneurysm hemodynamics using patient-specific endovascular measurements

Computational hemodynamic simulations of cerebral aneurysms have traditionally relied on stereotypical boundary conditions (such as blood flow velocity and blood pressure) derived from published values as patient-specific measurements are unavailable or difficult to collect. However, controversy persists over the necessity of incorporating such patient-specific conditions into computational analyses. We perform simulations using both endovascularly-derived patient-specific and typical literature-derived inflow and outflow boundary conditions. Detailed three-dimensional anatomical models of the cerebral vasculature are developed from rotational angiography data, and blood flow velocity and pressure are measured in situ by a dual-sensor pressure and velocity endovascular guidewire at multiple peri-aneurysmal locations in 10 unruptured cerebral aneurysms. These measurements are used to define inflow and outflow boundary conditions for computational hemodynamic models of the aneurysms. The additional in situ measurements which are not prescribed in the simulation are then used to assess the accuracy of the simulated flow velocity and pressure drop. Simulated velocities using patient-specific boundary conditions show good agreement with the guidewire measurements at measurement locations inside the domain, with no bias in the agreement and a random scatter of ≈25%. Simulated velocities using the simplified, literature-derived values show a systematic bias and over-predicted velocity by ≈30% with a random scatter of ≈40%. Computational hemodynamics using endovascularly measured patient-specific boundary conditions have the potential to improve treatment predictions as they provide more accurate and precise results of the aneurysmal hemodynamics than those based on commonly accepted reference values for boundary conditions.

Variability of computational fluid dynamics solutions for pressure and flow in a giant aneurysm: the ASME 2012 Summer Bioengineering Conference CFD Challenge

Stimulated by a recent controversy regarding pressure drops predicted in a giant aneurysm with a proximal stenosis, the present study sought to assess variability in the prediction of pressures and flow by a wide variety of research groups. In phase I, lumen geometry, flow rates, and fluid properties were specified, leaving each research group to choose their solver, discretization, and solution strategies. Variability was assessed by having each group interpolate their results onto a standardized mesh and centerline. For phase II, a physical model of the geometry was constructed, from which pressure and flow rates were measured. Groups repeated their simulations using a geometry reconstructed from a micro-computed tomography (CT) scan of the physical model with the measured flow rates and fluid properties. Phase I results from 25 groups demonstrated remarkable consistency in the pressure patterns, with the majority predicting peak systolic pressure drops within 8% of each other. Aneurysm sac flow patterns were more variable with only a few groups reporting peak systolic flow instabilities owing to their use of high temporal resolutions. Variability for phase II was comparable, and the median predicted pressure drops were within a few millimeters of mercury of the measured values but only after accounting for submillimeter errors in the reconstruction of the life-sized flow model from micro-CT. In summary, pressure can be predicted with consistency by CFD across a wide range of solvers and solution strategies, but this may not hold true for specific flow patterns or derived quantities. Future challenges are needed and should focus on hemodynamic quantities thought to be of clinical interest.

Cerebral aneurysms treated with flow-diverting stents: computational models with intravascular blood flow measurements

BACKGROUND AND PURPOSE

Computational fluid dynamics modeling is useful in the study of the hemodynamic environment of cerebral aneurysms, but patient-specific measurements of boundary conditions, such as blood flow velocity and pressure, have not been previously applied to the study of flow-diverting stents. We integrated patient-specific intravascular blood flow velocity and pressure measurements into computational models of aneurysms before and after treatment with flow-diverting stents to determine stent effects on aneurysm hemodynamics.

MATERIALS AND METHODS

Blood flow velocity and pressure were measured in peri-aneurysmal locations by use of an intravascular dual-sensor pressure and Doppler velocity guidewire before and after flow-diverting stent treatment of 4 unruptured cerebral aneurysms. These measurements defined inflow and outflow boundary conditions for computational models. Intra-aneurysmal flow rates, wall shear stress, and wall shear stress gradient were calculated.

RESULTS

Measurements of inflow velocity and outflow pressure were successful in all 4 patients. Computational models incorporating these measurements demonstrated significant reductions in intra-aneurysmal wall shear stress and wall shear stress gradient and a trend in reduced intra-aneurysmal blood flow.

CONCLUSIONS

Integration of intravascular dual-sensor guidewire measurements of blood flow velocity and blood pressure provided patient-specific computational models of cerebral aneurysms. Aneurysm treatment with flow-diverting stents reduces blood flow and hemodynamic shear stress in the aneurysm dome.

Hemodynamic characteristics at the rupture site of cerebral aneurysms: a case study

BACKGROUND AND IMPORTANCE

Hemodynamics play an important role in the mechanisms of aneurysm formation, growth, and rupture. However, little is known about the hemodynamics of rupture sites.

CLINICAL PRESENTATION

We incidentally acquired 3-dimensional images before and at the moment of rebleeding of a cerebral aneurysm in a patient. Comparison of these 2 images enabled precise identification of the rupture site. On the basis of computational fluid dynamics simulation, we propose that there are characteristic hemodynamic parameters of the rupture site in cerebral aneurysms. We evaluated flow velocity, wall shear stress (WSS), pressure, and the oscillatory shear index to determine characteristic parameters at the rupture site. Among the hemodynamic parameters in the cardiac cycle, the rupture site was most markedly distinguished by a combination of low WSS at end diastole and high pressure at peak systole. The flow patterns around the rupture site uniquely changed in the cardiac cycle. The rupture site was an impingement zone at peak systole. Flow separation at the rupture site was observed at end diastole.

CONCLUSION

In this case, a region with low WSS at end diastole and high pressure at peak systole was at the rupture site. A possible mechanism of rupture in this particular aneurysm is that low WSS at end diastole caused degeneration and thinning of the aneurysm wall and that high pressure at peak systole (impingement zone) resulted in rupture of the thinning wall.

A workflow for patient-individualized virtual angiogram generation based on CFD simulation

Increasing interest is drawn on hemodynamic parameters for classifying the risk of rupture as well as treatment planning of cerebral aneurysms. A proposed method to obtain quantities such as wall shear stress, pressure, and blood flow velocity is to numerically simulate the blood flow using computational fluid dynamics (CFD) methods. For the validation of those calculated quantities, virtually generated angiograms, based on the CFD results, are increasingly used for a subsequent comparison with real, acquired angiograms. For the generation of virtual angiograms, several patient-specific parameters have to be incorporated to obtain virtual angiograms which match the acquired angiograms as best as possible. For this purpose, a workflow is presented and demonstrated involving multiple phantom and patient cases.