Can sites prone to flow induced vascular complications in a-v fistulas be assessed using computational fluid dynamics?

Arterio-venous fistulas (shunts between arteries and veins) are the preferred vascular access for hemodialysis. Despite their superior patency, compared with synthetic tubes and grafts, functional problems and inadequate flow rates are the common complications. Local flow conditions, in particular low and oscillating wall shear stresses (WSS), are central to vascular problems and a robust framework for analyzing flow conditions in vascular structures could provide an understanding of the mechanisms leading to vascular complications, such as stenoses, aneurisms, and thromboses. We hypothesize that a validated computational fluid dynamics (CFD) framework can be used to identify critical fistula configurations with elevated risk of complications. Therefore, the aim of the present study was to develop a CFD framework for analyzing fluid flow in complex vascular structures, such as arterio-venous fistulas validated by comparisons of in vitro volume flows with CFD results and flow fields from ultrasound scans with CFD simulations. Volume flows measured in vitro and CFD data differed quantitatively. However, good relative correlations exist between the data using logarithmic scales. Qualitatively, visual comparisons between ultrasound and CFD images showed good agreement between the two methods. In addition, WSS levels and the oscillatory shear index (OSI) were calculated and visualized on the model surface. The method was successfully validated and the method is deemed suitable for more thorough investigations into the field of vascular complications in a-v fistulas.

In-vivo examples of flow patterns with the fast vector velocity ultrasound method.

PURPOSE: Conventional ultrasound methods for acquiring color flow images of the blood motion are limited by a relatively low frame rate and are restricted to only giving velocity estimates along the ultrasound beam direction. To circumvent these limitations, the Plane Wave Excitation (PWE) method has been proposed. MATERIAL AND METHODS: The PWE method can estimate the 2D vector velocity of the blood with a high frame rate. Vector velocity estimates are acquired by using the following approach: The ultrasound is not focused during the ultrasound transmission, and a full speckle image of the blood can be acquired for each pulse emission. The pulse is a 13 bit Barker code transmitted simultaneously from each transducer element. The 2D vector velocity of the blood is found using 2D speckle tracking between segments in consecutive speckle images. Implemented on the experimental scanner RASMUS and using a 100 CPU linux cluster for post processing, PWE can achieve a frame of 100 Hz where one vector velocity sequence of approximately 3 sec, takes 10 h to store and 48 h to process. In this paper a case study is presented of in-vivo vector velocity estimates in different complex vessel geometries. RESULTS: The flow patterns of six bifurcations and two veins were investigated. It was shown: 1. that a stable vortex in the carotid bulb was present opposed to other examined bifurcations, 2. that retrograde flow was present in the superficial branch of the femoral artery during diastole, 3. that retrograde flow was present in the subclavian artery and antegrade in the common carotid artery during diastole, 4. that vortices were formed in the sinus pockets behind the venous valves in both antegrade and retrograde flow, and 5. that secondary flow was present in various vessels. CONCLUSION: Using a fast vector velocity ultrasound method, in-vivo scans have been recorded where complex flow patterns were visualized in greater detail than previously visualized by conventional color flow imaging techniques.

In vivo comparison of three ultrasound vector velocity techniques to MR phase contrast angiography.

The objective of this paper is to validate angle independent vector velocity methods for blood velocity estimation. Conventional Doppler ultrasound (US) only estimates the blood velocity along the US beam direction where the estimate is angle corrected assuming laminar flow parallel to vessel boundaries. This results in incorrect blood velocity estimates, when angle of insonation approaches 90 degrees or when blood flow is non-laminar. Three angle independent vector velocity methods are evaluated in this paper: directional beamforming (DB), synthetic aperture flow imaging (STA) and transverse oscillation (TO). The performances of the three methods were investigated by measuring the stroke volume in the right common carotid artery of 11 healthy volunteers with magnetic resonance phase contrast angiography (MRA) as reference. The correlation with confidence intervals (CI) between the three vector velocity methods and MRA were: DB vs. MRA: R=0.84 (p<0.01, 95% CI: 0.49-0.96); STA vs. MRA: R=0.71 (p<0.05, 95% CI: 0.19-0.92) and TO vs. MRA: R=0.91 (p<0.01, 95% CI: 0.69-0.98). No significant differences were observed for any of the three comparisons (DB vs. MRA: p=0.65; STA vs. MRA: p=0.24; TO vs. MRA: p=0.36). Bland-Altman plots were additionally constructed, and mean differences with limits of agreements (LoA) for the three comparisons were: DB vs. MRA=0.17 ml (95% CI: -0.61-0.95) with LoA=-2.11-2.44 ml; STA vs. MRA=-0.55 ml (95% CI: -1.54-0.43) with LoA=-3.42-2.32 ml; TO vs. MRA=0.24 ml (95% CI: -0.32-0.81) with LoA=-1.41-1.90 ml. According to the results, reliable volume flow estimates can be obtained with all three methods. The three US vector velocity techniques can yield quantitative insight into flow dynamics and visualize complex flow patterns, which potentially can give the clinician a novel tool for cardiovascular disease assessment.

Pattern dynamics of vortex ripples in sand: nonlinear modeling and experimental validation.

Vortex ripples in sand are studied experimentally in a one-dimensional setup with periodic boundary conditions. The nonlinear evolution, far from the onset of instability, is analyzed in the framework of a simple model developed for homogeneous patterns. The interaction function describing the mass transport between neighboring ripples is extracted from experimental runs using a recently proposed method for data analysis, and the predictions of the model are compared to the experiment. An analytic explanation of the wavelength selection mechanism in the model is provided, and the width of the stable band of ripples is measured.