Hemodynamics in a three-dimensional printed aortic model: a comparison of four-dimensional phase-contrast magnetic resonance and image-based computational fluid dynamics.

OBJECTIVE: This study aims to compare an electrocardiogram (ECG)-gated four-dimensional (4D) phase-contrast (PC) magnetic resonance imaging (MRI) technique and computational fluid dynamics (CFD) using variables controlled in a laboratory environment to minimize bias factors. MATERIALS AND METHODS: Data from 4D PC-MRI were compared with computational fluid dynamics using steady and pulsatile flows at various inlet velocities. Anatomically realistic models for a normal aorta, a penetrating atherosclerotic ulcer, and an abdominal aortic aneurysm were constructed using a three-dimensional printer. RESULTS: For the normal aorta model, the errors in the peak and the average velocities were within 5%. The peak velocities of the penetrating atherosclerotic ulcer and the abdominal aortic aneurysm models displayed a more extensive range of differences because of the high-speed and vortical fluid flows generated by the shape of the blood vessel. However, the average velocities revealed only relatively minor differences. CONCLUSIONS: This study compared the characteristics of PC-MRI and CFD through a phantom study that only included controllable experimental parameters. Based on these results, 4D PC-MRI and CFD are powerful tools for analyzing blood flow patterns in vivo. However, there is room for future developments to improve velocity measurement accuracy.

Analysis of the time-velocity curve in phase-contrast magnetic resonance imaging: a phantom study.

The aim of this study was to analyze the characteristics of time-velocity curve acquired by phase-contrast magnetic resonance imaging (PC-MRI) using an in-vitro flow model as a reference for hemodynamic studies. The time- velocity curves of the PC-MRI were compared with Doppler ultrasonography (US) and also compared with those obtained in the electromagnetic flowmeter. The correlation between techniques was analyzed using an electromagnetic flowmeter as a reference standard; the maximum, minimum, and average velocities, full-width at half-maximum (FWHM), and ascending gradient (AG) were measured from time-velocity curves. The correlations between an electromagnetic flowmeter and the respective measurement technique for the PC-MRI and Doppler US were found to be high (mean R2 > 0.9, p < 0.05). These results indicate that these measurement techniques are useful for measuring blood flow information and reflect actual flow. The PC-MRI was the best fit for the minimum velocity and FWHM, and the maximum velocity and AG were the best fit for Doppler US. The PC-MRI showed lower maximum velocity value and higher minimum velocity value than Doppler US. Therefore, PC-MRI demonstrates more obtuse time-velocity curve than Doppler US. In addition, the time- velocity curve of PC-MRI could be calibrated by introducing formulae that can convert each measurement value to a reference standard value within a 10% error. The PC-MRI can be used to estimate the Doppler US using this formula.

Hemodynamic Parameters and Early Intimal Thickening in Branching Blood Vessels.

Intimal thickening due to atherosclerotic lesions or intimal hyperplasia in medium to large blood vessels is a major contributor to heart disease, the leading cause of death in the Western World. Balloon angioplasty with stenting, bypass surgery, and endarterectomy (with or without patch reconstruction) are some of the techniques currently applied to occluded blood vessels. On the basis of the preponderance of clinical evidence that disturbed flow patterns play a key role in the onset and progression of atherosclerosis and intimal hyperplasia, it is of interest to analyze suitable hemodynamic wall parameters that indicate susceptible sites of intimal thickening and/or favorable conditions for thrombi formation. These parameters, based on the wall shear stress, wall pressure, or particle deposition, are applied to interpret experimental/clinical observations of intimal thickening. Utilizing the parameters as "indicator" functions, internal branching blood vessel geometries are analyzed and possibly altered for different purposes: early detection of possibly highly stenosed vessel segments, prediction of future disease progression, and vessel redesign to potentially improve long-term patency rates. At the present time, the focus is on the identification of susceptible sites in branching blood vessels and their subsequent redesign, employing hemodynamic wall parameters. Specifically, the time-averaged wall shear stress (WSS), its spatial gradient (WSSG), the oscillatory shear index (OSI), and the wall shear stress angle gradient (WSSAG) are compared with experimental data for an aortoceliac junction. Then, the OSI, wall particle density (WPD), and WSSAG are segmentally averaged for different carotid artery bifurcations and compared with clinical data of intimal thickening. The third branching blood vessel under consideration is the graft-to-vein anastomosis of a vascular access graft Suggested redesigns reduce several hemodynamic parameters (i.e., the WSSG, WSSAG, and normal pressure gradient [NPG]), thereby reducing the likelihood of restenosis, especially near the critical toe region.

Computational modeling of HHH therapy and impact of blood pressure and hematocrit.

BACKGROUND: After an aneurysmal subarachnoid hemorrhage, cerebral microcirculatory changes occur as a result cerebral vasospasm. The objective of this study is to investigate, with a computational model, how various degrees of vasospasm are influenced by increasing the mean blood pressure and decreasing the blood viscosity. METHODS: Using ANSYS CFX software, a computational model was constructed to simulate steady-state fully developed laminar blood flow through a rigid wall system consisting of the internal carotid artery (ICA), anterior cerebral artery, posterior cerebral artery, and middle cerebral artery (MCA). The MCA was selected for the site of a single acute vasospasm. Five severities of vasospasm were studied: 3 mm (normal), 2.5, 2, 1.5, and 1 mm. The ICA was assumed to have a constant inlet flow rate of 315 mL/min. The anterior cerebral artery and posterior cerebral artery were assumed to have constant outlet flow rates of 105 mL/min and 30 mL/min, respectively. The MCA was assumed to have a constant outlet pressure of 92 mL/min. Two different hematocrits, 45% and 32%, were simulated using the models. RESULTS: For a hematocrit of 45, the mean ICA inlet pressure required to pump blood through the system was 104 mm Hg for the 3-mm diameter MCA and 105, 108, 116, and 158 mm Hg for vasospasm diameters of 2.5, 2, 1.5, and 1 mm, respectively. For a hematocrit of 32, the mean ICA inlet pressure required was 102, 103, 105, 113, and 152 mm Hg, respectively. CONCLUSIONS: The MCA required a large increase in mean ICA inlet pressure for vasospasm diameters less than 1.5 mm, which suggests that for vasospasms more than 50% diameter reduction, the blood pressure must be increased dramatically. Decreasing the hematocrit had minimal impact on blood flow in a constricted vessel.

Computational analysis of effects of external carotid artery flow and occlusion on adverse carotid bifurcation hemodynamics.

OBJECTIVE: This is a computational analysis of the effects of external carotid artery (ECA) flow, waveform, and occlusion geometry on two hemodynamic wall parameters associated with intimal hyperplasia and atherosclerosis. Study design Transient three-dimensional fluid mechanics analysis was applied to a standard carotid artery bifurcation. Mean internal carotid artery (ICA) flow was maintained at 236 mL/min with a normal waveform. ECA flow was increased from zero to 151 mL/min (64% of ICA flow) with both a normal biphasic waveform and a damped waveform. Geometry of five ECA occlusions was studied: distal, proximal stump, smooth, smooth without carotid sinus, and optimal reconstruction.Primary outcome measures Two time-averaged and area-averaged hemodynamic wall parameters were computed from the velocity and wall shear stress (WSS) solutions, ie, wall shear stress angle gradient (WSSAG) and oscillatory shear index (OSI). Both local and area-averaged hemodynamic wall parameters were computed for the distal common carotid artery (CCA) and the proximal ICA. RESULTS: When ECA flow with a normal waveform is increased from zero to 151 mL/min, area-averaged WSS values increase in the CCA, from 3.0 to 4.4 dynes/cm(2) (46%), and in the ICA, from 16.5 to 17.1 dynes/cm(2) (4%); minimum local WSS values in the carotid sinus remain less than 1 dyne/cm(2); maximum local values of WSSAG and OSI are observed in the carotid sinus and increase from 3.5 to 9.1 radian/cm (160%) and 0.23 to 0.46 (100%), respectively; CCA plus ICA area-averaged WSSAG increases by 52%, and OSI increases by 144%; and damping of the ECA waveform has little effect on local or area-averaged WSSAG but reduces OSI to 68%. When the ECA is occluded, the minimum local WSS in the carotid sinus is less than 1 dyne/cm(2). However, if the carotid sinus is removed or the CCA-ICA geometry hemodynamically optimized, the minimum WSS is approximately 4 dynes/cm(2). Similarly, eliminating the carotid sinus markedly reduces local maximum WSSAG, from 3.0-3.5 radian/cm to 0.3 radian/cm, and reduces local maximum OSI from 0.22-0.49 to 0.04. Area-averaged WSSAG and OSI over the CCA and ICA are reduced by approximately 50% with elimination of the carotid sinus. CONCLUSIONS: The degree of adverse carotid bifurcation hemodynamics as measured with WSSAG and OSI is directly proportional to ECA flow. The marked difference in normal ICA and ECA flow waveforms does not contribute to adverse wall hemodynamics. Location of an ECA occlusion (distal, proximal, stump, smooth) does not affect adverse carotid hemodynamics; however, marked improvement is obtained with elimination of the carotid sinus.

Acoustic streaming induced by ultrasonic flexural vibrations and associated enhancement of convective heat transfer.

Acoustic streaming induced by ultrasonic flexural vibrations and the associated convection enhancement are investigated. Acoustic streaming pattern, streaming velocity, and associated heat transfer characteristics are experimentally observed. Moreover, analytical analysis based on Nyborg's formulation is performed along with computational fluid dynamics (CFD) simulation using a numerical solver CFX 4.3. Two distinctive acoustic streaming patterns in half-wavelength of the flexural vibrations are observed, which agree well with the theory. However, acoustic streaming velocities obtained from CFD simulation, based on the incompressible flow assumption, exceed the theoretically estimated velocity by a factor ranging from 10 to 100, depending upon the location along the beam. Both CFD simulation and analytical analysis reveal that the acoustic streaming velocity is proportional to the square of the vibration amplitude and the wavelength of the vibrating beam that decreases with the excitation frequency. It is observed that the streaming velocity decreases with the excitation frequency. Also, with an open-ended channel, a substantial increase in streaming velocity is observed from CFD simulations. Using acoustic streaming, a temperature drop of 40 degrees C with a vibration amplitude of 25 microm at 28.4 kHz is experimentally achieved.