Quantifying turbulent wall shear stress in a subject specific human aorta using large eddy simulation.

In this study, large-eddy simulation (LES) is employed to calculate the disturbed flow field and the wall shear stress (WSS) in a subject specific human aorta. Velocity and geometry measurements using magnetic resonance imaging (MRI) are taken as input to the model to provide accurate boundary conditions and to assure the physiological relevance. In total, 50 consecutive cardiac cycles were simulated from which a phase average was computed to get a statistically reliable result. A decomposition similar to Reynolds decomposition is introduced, where the WSS signal is divided into a pulsating part (due to the mass flow rate) and a fluctuating part (originating from the disturbed flow). Oscillatory shear index (OSI) is plotted against time-averaged WSS in a novel way, and locations on the aortic wall where elevated values existed could easily be found. In general, high and oscillating WSS values were found in the vicinity of the branches in the aortic arch, while low and oscillating WSS were present in the inner curvature of the descending aorta. The decomposition of WSS into a pulsating and a fluctuating part increases the understanding of how WSS affects the aortic wall, which enables both qualitative and quantitative comparisons.

Evaluation of aortic geometries created by magnetic resonance imaging data in healthy volunteers.

INTRODUCTION: The development of atherosclerotic plaques has been associated with the patterns of wall shear stress (WSS). However, much is still uncertain with the methods used to calculate WSS. Correct vessel geometries are mandatory to get reliable estimations, and the purpose of this study was to evaluate an in vivo method for creating aortic 3D geometry in human based on data from magnetic resonance imaging (MRI) with ultrasound as reference. METHODS: The aortas of ten healthy men, 23·4 ± 1·6 years of age, were examined with a 1·5 T MRI system using a 3D gadolinium-enhanced gradient-echo sequence. Three-dimensional geometries were created using manual segmentation of images. Lumen diameters (LD) were measured in the abdominal aorta (AA) and the thoracic aorta (TA) with non-invasive B-mode ultrasound as a reference. RESULTS: The anteroposterior diameter of the AA was 13·6 ± 1·1 mm for the MRI and 13·8 ± 1·3 mm for the ultrasound (NS). Intraobserver variability (CV) for MRI and ultrasound was <0·92% and <0·40%, respectively. Interobserver variability for MRI and ultrasound was 0·96% and 0·56%, respectively. The diameter of the TA was 19·2 ± 1·4 mm for the MRI, and the intraobserver variability (CV) was <0·78% and interobserver variability (CV) was 0·92%. CONCLUSION: Specific arterial geometries can be constructed with a high degree of accuracy using MRI. This indicates that the MRI geometries may be used to create realistic and correct geometries in the calculation of WSS in the aorta of human.

Quantifying turbulent wall shear stress in a stenosed pipe using large eddy simulation.

Large eddy simulation was applied for flow of Re=2000 in a stenosed pipe in order to undertake a thorough investigation of the wall shear stress (WSS) in turbulent flow. A decomposition of the WSS into time averaged and fluctuating components is proposed. It was concluded that a scale resolving technique is required to completely describe the WSS pattern in a subject specific vessel model, since the poststenotic region was dominated by large axial and circumferential fluctuations. Three poststenotic regions of different WSS characteristics were identified. The recirculation zone was subject to a time averaged WSS in the retrograde direction and large fluctuations. After reattachment there was an antegrade shear and smaller fluctuations than in the recirculation zone. At the reattachment the fluctuations were the largest, but no direction dominated over time. Due to symmetry the circumferential time average was always zero. Thus, in a blood vessel, the axial fluctuations would affect endothelial cells in a stretched state, whereas the circumferential fluctuations would act in a relaxed direction.

Simulation of phase contrast MRI of turbulent flow.

Phase contrast MRI is a powerful tool for the assessment of blood flow. However, especially in the highly complex and turbulent flow that accompanies many cardiovascular diseases, phase contrast MRI may suffer from artifacts. Simulation of phase contrast MRI of turbulent flow could increase our understanding of phase contrast MRI artifacts in turbulent flows and facilitate the development of phase contrast MRI methods for the assessment of turbulent blood flow. We present a method for the simulation of phase contrast MRI measurements of turbulent flow. The method uses an Eulerian-Lagrangian approach, in which spin particle trajectories are computed from time-resolved large eddy simulations. The Bloch equations are solved for each spin for a frame of reference moving along the spins trajectory. The method was validated by comparison with phase contrast MRI measurements of velocity and intravoxel velocity standard deviation (IVSD) on a flow phantom consisting of a straight rigid pipe with a stenosis. Turbulence related artifacts, such as signal drop and ghosting, could be recognized in the measurements as well as in the simulations. The velocity and the IVSD obtained from the magnitude of the phase contrast MRI simulations agreed well with the measurements.

On MRI turbulence quantification.

Turbulent flow, characterized by velocity fluctuations, accompanies many forms of cardiovascular disease and may contribute to their progression and hemodynamic consequences. Several studies have investigated the effects of turbulence on the magnetic resonance imaging (MRI) signal. Quantitative MRI turbulence measurements have recently been shown to have great potential for application both in human cardiovascular flow and in engineering flow. In this article, potential pitfalls and sources of error in MRI turbulence measurements are theoretically and numerically investigated. Data acquisition strategies suitable for turbulence quantification are outlined. The results show that the sensitivity of MRI turbulence measurements to intravoxel mean velocity variations is negligible, but that noise may degrade the estimates if the turbulence encoding parameter is set improperly. Different approaches for utilizing a given amount of scan time were shown to influence the dynamic range and the uncertainty in the turbulence estimates due to noise. The findings reported in this work may be valuable for both in vitro and in vivo studies employing MRI methods for turbulence quantification.

Feasibility of patient specific aortic blood flow CFD simulation.

Patient specific modelling of the blood flow through the human aorta is performed using computational fluid dynamics (CFD) and magnetic resonance imaging (MRI). Velocity patterns are compared between computer simulations and measurements. The workflow includes several steps: MRI measurement to obtain both geometry and velocity, an automatic levelset segmentation followed by meshing of the geometrical model and CFD setup to perform the simulations follwed by the actual simulations. The computational results agree well with the measured data.