MRI and CFD studies of pulsatile flow in healthy and stenosed carotid bifurcation models

Pulsatile flow was studied in physiologically realistic models of a normal and a moderately stenosed (30% diameter reduction) human carotid bifurcation. Time-resolved velocity measurements were made using magnetic resonance imaging, from which wall shear stress (WSS) vectors were calculated. Velocity measurements in the inflow and outflow regions were also used as boundary conditions for a computational fluid dynamics (CFD) model. Experimental flow patterns and derived WSS vectors were compared qualitatively with the corresponding CFD predictions. In the stenosed phantom, flow in the bulb region of the "internal carotid artery" was concentrated along the outer wall, with a region of low and recirculating flow near the inner wall. In the normal phantom, the converse was found, with a low flow region near the outer wall of the bulb. Time-averaged WSS and oscillatory shear index were also markedly different for the two phantoms.     

Effect of common carotid artery inlet length on normal carotid bifurcation hemodynamics

Controversy exists regarding the suitability of fully developed versus measured inlet velocity profiles for image-based computational fluid dynamics (CFD) studies of carotid bifurcation hemodynamics. Here, we attempt to resolve this by investigating the impact of the reconstructed common carotid artery (CCA) inlet length on computed metrics of "disturbed" flow. Twelve normal carotid bifurcation geometries were reconstructed from contrast-enhanced angiograms acquired as part of the Vascular Aging--The Link That Bridges Age to Atherosclerosis study (VALIDATE). The right carotid artery lumen geometry was reconstructed from its brachiocephalic origin to well above the bifurcation, and the CCA was truncated objectively at locations one, three, five, and seven diameters proximal to where it flares into the bifurcation. Relative to the simulations carried out using the full CCA, models truncated at one CCA diameter strongly overestimated the amount of disturbed flow. Substantial improvement was offered by using three CCA diameters, with only minor further improvement using five CCA diameters. With seven CCA diameters, the amounts of disturbed flow agreed unambiguously with those predicted by the corresponding full-length models. Based on these findings, we recommend that image-based CFD models of the carotid bifurcation should incorporate at least three diameters of CCA length if fully developed velocity profiles are to be imposed at the inlet. The need for imposing measured inlet velocity profiles would seem to be relevant only for those cases where the CCA is severely truncated.     

Comparison of 4D Phase-Contrast MRI Flow Measurements to Computational Fluid Dynamics Simulations of Cerebrospinal Fluid Motion in the Cervical Spine

Cerebrospinal fluid (CSF) dynamics in the cervical spinal subarachnoid space (SSS) have been thought to be important to help diagnose and assess craniospinal disorders such as Chiari I malformation (CM). In this study we obtained time-resolved three directional velocity encoded phase-contrast MRI (4D PC MRI) in three healthy volunteers and four CM patients and compared the 4D PC MRI measurements to subject-specific 3D computational fluid dynamics (CFD) simulations. The CFD simulations considered the geometry to be rigid-walled and did not include small anatomical structures such as nerve roots, denticulate ligaments and arachnoid trabeculae. Results were compared at nine axial planes along the cervical SSS in terms of peak CSF velocities in both the cranial and caudal direction and visual interpretation of thru-plane velocity profiles. 4D PC MRI peak CSF velocities were consistently greater than the CFD peak velocities and these differences were more pronounced in CM patients than in healthy subjects. In the upper cervical SSS of CM patients the 4D PC MRI quantified stronger fluid jets than the CFD. Visual interpretation of the 4D PC MRI thru-plane velocity profiles showed greater pulsatile movement of CSF in the anterior SSS in comparison to the posterior and reduction in local CSF velocities near nerve roots. CFD velocity profiles were relatively uniform around the spinal cord for all subjects. This study represents the first comparison of 4D PC MRI measurements to CFD of CSF flow in the cervical SSS. The results highlight the utility of 4D PC MRI for evaluation of complex CSF dynamics and the need for improvement of CFD methodology. Future studies are needed to investigate whether integration of fine anatomical structures and gross motion of the brain and/or spinal cord into the computational model will lead to a better agreement between the two techniques.     

Effect of inlet velocity profiles on patient-specific computational fluid dynamics simulations of the carotid bifurcation

Patient-specific computational fluid dynamics (CFD) is a powerful tool for researching the role of blood flow in disease processes. Modern clinical imaging technology such as MRI and CT can provide high resolution information about vessel geometry, but in many situations, patient-specific inlet velocity information is not available. In these situations, a simplified velocity profile must be selected. We studied how idealized inlet velocity profiles (blunt, parabolic, and Womersley flow) affect patient-specific CFD results when compared to simulations employing a "reference standard" of the patient's own measured velocity profile in the carotid bifurcation. To place the magnitude of these effects in context, we also investigated the effect of geometry and the use of subject-specific flow waveform on the CFD results. We quantified these differences by examining the pointwise percent error of the mean wall shear stress (WSS) and the oscillatory shear index (OSI) and by computing the intra-class correlation coefficient (ICC) between axial profiles of the mean WSS and OSI in the internal carotid artery bulb. The parabolic inlet velocity profile produced the most similar mean WSS and OSI to simulations employing the real patient-specific inlet velocity profile. However, anatomic variation in vessel geometry and the use of a nonpatient-specific flow waveform both affected the WSS and OSI results more than did the choice of inlet velocity profile. Although careful selection of boundary conditions is essential for all CFD analysis, accurate patient-specific geometry reconstruction and measurement of vessel flow rate waveform are more important than the choice of velocity profile. A parabolic velocity profile provided results most similar to the patient-specific velocity profile.     

Visualization and quantification of the human blood flow by magnetic resonance imaging.

Magnetic Resonance Imaging (MRI) offers new possibilities for the visualization and the noninvasive quantification of the blood flow in human vessels. By the application of conventional gradient echo sequences with electrocardiographic gating on a 1.5 Tesla whole body MRI system the flow induced phase shifts in the ascending and the abdominal aorta are analyzed. The instantaneous two-dimensional velocity profiles and the instantaneous flow rates are determined in a series of subsequent images with high temporal resolution throughout the cardiac cycle. For the flow analysis in further vessels and for the analysis of more complex flow patterns, as they occur in bifurcations or stenoses, a new MR flow imaging technique called FAcE with extremely short echo times is introduced and the first results of flow examinations in a bifurcation phantom and in the carotid artery are presented.     

Accuracy and reproducibility of CFD predicted wall shear stress using 3D ultrasound images

Computational fluid dynamics (CFD) flow simulation techniques have the potential to enhance our understanding of how haemodynamic factors are involved in atherosclerosis. Recently, 3D ultrasound has emerged as an alternative to other 3D imaging techniques, such as magnetic resonance angiography (MRA). The method can be used to generate realistic vascular geometry suitable for CFD simulations. In order to assess accuracy and reproducibility of the procedure from image acquisition to reconstruction to CFD simulation, a human carotid artery bifurcation phantom was scanned three times using 3D ultrasound. The geometry was reconstructed and flow simulations were carried out on the three sets as well as on a model generated using computer aided design (CAD) from the geometric information given by the manufacturer. It was found that the three reconstructed sets showed good reproducibility as well as satisfactory quantitative agreement with the CAD model. Analyzing two selected locations probably representing the 'worst cases,' accuracy comparing ultrasound and CAD reconstructed models was estimated to be between 7.2% and 7.7% of the maximum instantaneous WSS and reproducibility comparing the three scans to be between 8.2% and 10.7% of their average maximum.     

Validation of CFD predictions of flow in a 3D alveolated bend with experimental data

Verifying numerical predictions with experimental data is an important aspect of any modeling studies. In the case of the lung, the absence of direct in vivo flow measurements makes such verification almost impossible. We performed computational fluid dynamics (CFD) simulations in a 3D scaled-up model of an alveolated bend with rigid walls that incorporated essential geometrical characteristics of human alveolar structures and compared numerical predictions with experimental flow measurements made in the same model by particle image velocimetry (PIV). Flow in both models was representative of acinar flow during normal breathing (0.82ml/s). The experimental model was built in silicone and silicone oil was used as the carrier fluid. Flow measurements were obtained by an ensemble averaging procedure. CFD simulation was performed with STAR-CCM+ (CD-Adapco) using a polyhedral unstructured mesh. Velocity profiles in the central duct were parabolic and no bulk convection existed between the central duct and the alveoli. Velocities inside the alveoli were approximately 2 orders of magnitude smaller than the mean velocity in the central duct. CFD data agreed well with those obtained by PIV. In the central duct, data agreed within 1%. The maximum simulated velocity along the centerline of the model was 0.5% larger than measured experimentally. In the alveolar cavities, data agreed within 15% on average. This suggests that CFD techniques can satisfactorily predict acinar-type flow. Such a validation ensure a great degree of confidence in the accuracy of predictions made in more complex models of the alveolar region of the lung using similar CFD techniques.     

Flow pattern analysis in a highly stenotic patient-specific carotid bifurcation model using a turbulence model

The aim of this study is to investigate the blood flow pattern in carotid bifurcation with a high degree of luminal stenosis, combining in vivo magnetic resonance imaging (MRI) and computational fluid dynamics (CFD). A newly developed two-equation transitional model was employed to evaluate wall shear stress (WSS) distribution and pressure drop across the stenosis, which are closely related to plaque vulnerability. A patient with an 80% left carotid stenosis was imaged using high resolution MRI, from which a patient-specific geometry was reconstructed and flow boundary conditions were acquired for CFD simulation. A transitional model was implemented to investigate the flow velocity and WSS distribution in the patient-specific model. The peak time-averaged WSS value of approximately 73 Pa was predicted by the transitional flow model, and the regions of high WSS occurred at the throat of the stenosis. High oscillatory shear index values up to 0.50 were present in a helical flow pattern from the outer wall of the internal carotid artery immediately after the throat. This study shows the potential suitability of a transitional turbulent flow model in capturing the flow phenomena in severely stenosed carotid arteries using patient-specific MRI data and provides the basis for further investigation of the links between haemodynamic variables and plaque vulnerability. It may be useful in the future for risk assessment of patients with carotid disease.     

MRI-based quantification of outflow boundary conditions for computational fluid dynamics of stenosed human carotid arteries

Accurate assessment of wall shear stress (WSS) is vital for studies on the pathogenesis of atherosclerosis. WSS distributions can be obtained by computational fluid dynamics (CFD) using patient-specific geometries and flow measurements. If patient-specific flow measurements are unavailable, in- and outflow have to be estimated, for instance by using Murray’s Law. It is currently unknown to what extent this law holds for carotid bifurcations, especially in cases where stenoses are involved. We performed flow measurements in the carotid bifurcation using phase-contrast MRI in patients with varying degrees of stenosis. An empirical relation between outflow and degree of area stenosis was determined and the outflow measurements were compared to estimations based on Murray’s Law. Furthermore, the influence of outflow conditions on the WSS distribution was studied.For bifurcations with an area stenosis smaller than 65%, the outflow ratio of the internal carotid artery (ICA) to the common carotid artery (CCA) was 0.62±0.12 while the outflow ratio of the external carotid artery (ECA) was 0.35±0.13. If the area stenosis was larger than 65%, the flow to the ICA decreased linearly to zero at 100% area stenosis. The empirical relation fitted the flow data well (R²=0.69), whereas Murray’s Law overestimated the flow to the ICA substantially for larger stenosis, resulting in an overestimation of the WSS. If patient-specific flow measurements of the carotid bifurcation are unavailable, estimation of the outflow ratio by the presented empirical relation will result in a good approximation of calculated WSS using CFD.     

Hemodynamics in the carotid artery bifurcation: a comparison between numerical simulations and in vitro MRI measurements

The presence of atherosclerotic plaques has been shown to be closely related to the vessel geometry. Studies on postmortem human arteries and on the experimental animal show positive correlation between the presence of plaque thickness and low shear stress, departure of unidirectional flow and regions of flow separation and recirculation. Numerical simulations of arterial blood flow and direct blood flow velocity measurements by magnetic resonance imaging (MRI) are two approaches for the assessment of arterial blood flow patterns. In order to verify that both approaches give equivalent results magnetic resonance velocity data measured in a compliant anatomical carotid bifurcation model were compared to the results of numerical simulations performed for a corresponding computational vessel model. Cross sectional axial velocity profiles were calculated and measured for the midsinus and endsinus internal carotid artery. At both locations a skewed velocity profile with slow velocities at the outer vessel wall, medium velocities at the side walls and high velocities at the flow divider (inner) wall were observed. Qualitative comparison of the axial velocity patterns revealed no significant differences between simulations and in vitro measurements. Even quantitative differences such as for axial peak flow velocities were less than 10%. Secondary flow patterns revealed some minor differences concerning the form of the vortices but maximum circumferential velocities were in the same range for both methods.     

Validation of computational fluid dynamics methodology used for human upper airway flow simulations

An anatomically accurate human upper airway model was constructed from multiple magnetic resonance imaging axial scans. This model was used to conduct detailed Computational Fluid Dynamics (CFD) simulations during expiration, to investigate the fluid flow in the airway regions where obstruction could occur. An identical physical model of the same airway was built using stereo lithography. Pressure and velocity measurements were conducted in the physical model. Both simulations and experiments were performed at a peak expiratory flow rate of 200 L/min. Several different numerical approaches within the FLUENT commercial software framework were used in the simulations; unsteady Large Eddy Simulation (LES), steady Reynolds-Averaged Navier-Stokes (RANS) with two-equation turbulence models (i.e. k?ε, standard k?ω, and k?ω Shear Stress Transport (SST)) and with one-equation Spalart–Allmaras model. The CFD predictions of the average wall static pressures at different locations along the airway wall were favorably compared with the experimental data. Among all the approaches, standard k?ω turbulence model resulted in the best agreement with the static pressure measurements, with an average error of ～20% over all ports. The highest positive pressures were observed in the retroglossal regions below the epiglottis, while the lowest negative pressures were recorded in the retropalatal region. The latter is a result of the airflow acceleration in the narrow retropalatal region. The largest pressure drop was observed at the tip of the soft palate. This location has the smallest cross section of the airway. The good agreement between the computations and the experimental results suggest that CFD simulations can be used to accurately compute aerodynamic flow characteristics of the upper airway.     

Model studies of nonsteady flow using magnetic resonance imaging

A bolus-tracking magnetic resonance imaging (MRI) method has been employed to measure velocity profiles for oscillatory flow with and without a steady flow component as well as pulsatile flow in an axisymmetric tube model. A range of flow conditions within normal physiological limits was tested. The imaged velocity profiles were observed to be generally in accord with theoretical predictions. Instantaneous flow rates calculated from the MR images agreed well with those assessed using an ultrasonic flowmeter. Because MRI is noninvasive and poses few risks to subjects, this technique is potentially useful for studying vascular hemodynamics in vivo.     

Proposition of an outflow boundary approach for carotid artery stenosis CFD simulation

The purpose of this study was to propose an innovative approach of setting outlet boundary conditions for the computational fluid dynamics (CFD) simulation of human common carotid arteries (CCAs) bifurcation based on the concept of energy loss minimisation at flow bifurcation. Comparisons between this new approach and previously reported boundary conditions were also made. The results showed that CFD simulation based on the proposed boundary conditions gave an accurate prediction of the critical stenosis ratio of carotid arteries (at around 65%). Other boundary conditions, such as the constant external pressure (P = 0) and constant outflow ratio, either overestimated or underestimated the critical stenosis ratio of carotid arteries. The patient-specific simulation results furthermore indicated that the calculated internal carotid artery flow ratio at CCA bifurcation (61%) coincided with the result obtained by clinical measurements through the use of Colour Doppler ultrasound.     

Novel method to correct displacement profile images and calculate flow vector fields using saturation-tagging MR imaging.

Saturation-tagging Magnetic Resonance (MR) imaging provides a simple and robust means to directly visualize displacement profiles within fluid flow fields. Although useful for velocity quantitation as well as for qualitative depiction of flow patterns in certain well-defined flow fields, the technique is prone to distortions due to oblique flow (misregistration artifact) and ambiguity of fluid vector trajectories in complex flow situations. A novel method is proposed whereby two images are acquired, differing in the temporal position of the phase encoding gradient. Theoretical analysis shows that from the paired images, distortion of the two-dimensional displacement profile can be corrected and fluid velocity vectors extracted, even if the flow directions are unknown. In the simpler case of flow oblique to the gradient principal axes, but with a known trajectory, only one image is necessary to correct the displacement profile distortion and extract the velocity information. MRI experiments in a straight tube model have been carried out to evaluate the feasibility of this method. Good agreement is achieved between the results from MR imaging and those predicted via computer simulation.     

Measurement of hemodynamics in human carotid artery using ultrasound and computational fluid dynamics.

The objective of the study was to investigate the feasibility of using computational fluid dynamic modeling (CFD) with noninvasive ultrasound measurements to determine time-variant three-dimensional (3D) carotid arterial hemodynamics in humans in vivo. The effects of hyperoxia and hypoxic hypercapnia on carotid artery local hemodynamics were examined by use of this approach. Six normotensive volunteers followed a double-blind randomized crossover design. Blood pressure, heart rate, and carotid blood flow were measured while subjects breathed normal air, a mixture of 5% CO(2) and 15% O(2) (hypoxic hypercapnia), and 100% O(2) (hyperoxia). Carotid artery geometry was reconstructed on the basis of B-mode ultrasound images by using purpose-built image processing software. Time-variant 3D carotid hemodynamics were estimated by using finite volume-based CFD. Systemic blood pressure was not significantly affected by hyperoxia or hypoxic hypercapnia, but heart rate decreased significantly with hyperoxia. There was an increase in diastolic flow velocity in the external carotid artery after hypoxic hypercapnia, but otherwise carotid blood flow velocities did not change significantly. Compared with normal air, hyperoxic conditions were associated with a decrease in the width of the region of flow separation in the external carotid artery. During hyperoxia, there was also an increase in the minimum and a decrease in maximum shear stress in the bifurcation and hence a reduction in cyclic variation in shear stress. Hypoxic hypercapnia was associated with a reduced duration of flow separation in the external carotid artery and an increase in the minimum shear stress without affecting the cyclic variation in shear stress. This study demonstrates the feasibility of using noninvasive ultrasound techniques in conjunction with CFD to describe time-variant 3D hemodynamics in the human carotid arterial bifurcation in vivo.     

Validation of CFD simulations of cerebral aneurysms with implication of geometric variations

BACKGROUND: Computational fluid dynamics (CFD) simulations using medical-image-based anatomical vascular geometry are now gaining clinical relevance. This study aimed at validating the CFD methodology for studying cerebral aneurysms by using particle image velocimetry (PIV) measurements, with a focus on the effects of small geometric variations in aneurysm models on the flow dynamics obtained with CFD. METHOD OF APPROACH: An experimental phantom was fabricated out of silicone elastomer to best mimic a spherical aneurysm model. PIV measurements were obtained from the phantom and compared with the CFD results from an ideal spherical aneurysm model (S1). These measurements were also compared with CFD results, based on the geometry reconstructed from three-dimensional images of the experimental phantom. We further performed CFD analysis on two geometric variations, S2 and S3, of the phantom to investigate the effects of small geometric variations on the aneurysmal flow field. Results. We found poor agreement between the CFD results from the ideal spherical aneurysm model and the PIV measurements from the phantom, including inconsistent secondary flow patterns. The CFD results based on the actual phantom geometry, however, matched well with the PIV measurements. CFD of models S2 and S3 produced qualitatively similar flow fields to that of the phantom but quantitatively significant changes in key hemodynamic parameters such as vorticity, positive circulation, and wall shear stress. CONCLUSION: CFD simulation results can closely match experimental measurements as long as both are performed on the same model geometry. Small geometric variations on the aneurysm model can significantly alter the flow-field and key hemodynamic parameters. Since medical images are subjected to geometric uncertainties, image-based patient-specific CFD results must be carefully scrutinized before providing clinical feedback.     

Targeted particle tracking in computational models of human carotid bifurcations

A significant and largely unsolved problem of computational fluid dynamics (CFD) simulation of flow in anatomically relevant geometries is that very few calculated pathlines pass through regions of complex flow. This in turn limits the ability of CFD-based simulations of imaging techniques (such as MRI) to correctly predict in vivo performance. In this work, I present two methods designed to overcome this filling problem, firstly, by releasing additional particles from areas of the flow inlet that lead directly to the complex flow region ("preferential seeding") and, secondly, by tracking particles both "downstream" and "upstream" from seed points within the complex flow region itself. I use the human carotid bifurcation as an example of complex blood flow that is of great clinical interest. Both idealized and healthy volunteer geometries are investigated. With uniform seeding in the inlet plane (in the common carotid artery (CCA)) of an idealized bifurcation geometry, approximately half the particles passed through the internal carotid artery (ICA) and half through the external carotid artery. However, of those particles entering the ICA, only 16% passed directly through the carotid bulb region. Preferential seeding from selected regions of the CCA was able to increase this figure to 47%. In the second method, seeding of particles within the carotid bulb region itself led to a very high proportion (97%) of pathlines running from CCA to ICA. Seeding of particles in the bulb plane of three healthy volunteer carotid bifurcation geometries led to much better filling of the bulb regions than by particles seeded at the inlet alone. In all cases, visualization of the origin and behavior of recirculating particles led to useful insights into the complex flow patterns. Both seeding methods produced significant improvements in filling the carotid bulb region with particle tracks compared with uniform seeding at the inlet and led to an improved understanding of the complex flow patterns. The methods described may be combined and are generally applicable to CFD studies of fluid and gas flow and are, therefore, of relevance in hemodynamics, respiratory mechanics, and medical imaging science.     

MRI strain imaging of the carotid artery: Present limitations and future challenges

Rupture of atherosclerotic plaques in the carotid artery is a main cause of stroke. Current diagnostics are not sufficient to identify all rupture-prone plaques, and studies have shown that biomechanical factors improve current plaque risk assessment. Strain imaging may be a valuable contribution to this risk assessment. MRI is a versatile imaging technique that offers various methods that are capable of measuring tissue strain. In this review, MR imaging techniques with displacement (DENSE), velocity (PC MRI), or strain (SENC) encoding protocols are discussed, together with post-processing techniques based on time-resolved MRI data. Although several MRI techniques are being developed to improve time-resolved MR imaging, current technical limitations related to spatial and temporal resolutions render MRI strain imaging currently unfit for carotid plaque strain evaluation.     

Numerical simulation of blood flow in the left ventricle and aortic sinus using magnetic resonance imaging and computational fluid dynamics

Understanding cardiac blood flow patterns has many applications in analysing haemodynamics and for the clinical assessment of heart function. In this study, numerical simulations of blood flow in a patient-specific anatomical model of the left ventricle (LV) and the aortic sinus are presented. The realistic 3D geometry of both LV and aortic sinus is extracted from the processing of magnetic resonance imaging (MRI). Furthermore, motion of inner walls of LV and aortic sinus is obtained from cine-MR image analysis and is used as a constraint to a numerical computational fluid dynamics (CFD) model based on the moving boundary approach. Arbitrary Lagrangian–Eulerian finite element method formulation is used for the numerical solution of the transient dynamic equations of the fluid domain. Simulation results include detailed flow characteristics such as velocity, pressure and wall shear stress for the whole domain. The aortic outflow is compared with data obtained by phase-contrast MRI. Good agreement was found between simulation results and these measurements.     

In vivo velocity vector imaging and time-resolved strain rate measurements in the wall of blood vessels using MRI

In this paper, we present a new approach for velocity vector imaging and time-resolved measurements of strain rates in the wall of human arteries using MRI and we prove its feasibility on two examples: in vitro on a phantom and in vivo on the carotid artery of a human subject. Results point out the promising potential of this approach for investigating the mechanics of arterial tissues in vivo.