This study reports the effect of N₂ and CO₂ bubbles on dilute to dense gas-liquid two-phase bubbly flow. A shadowgraph imaging technique captured bubble images at a high spatiotemporal resolution. The recordings of bubble images allow us to compute gas fraction distribution. It requires challenging segmentation and gas-liquid interface detection approaches in image processing. Hence a novel gas contour characterization technique has been introduced in this study that analyses light intensity per pixel for quantifying the effect of local gas volume fraction. The dominant gas structure and repetitive gas pattern have also been determined here using Fourier transform-based power spectral density and 2D cross-correlation functions, respectively. Gas-liquid flow regimes of dissolved CO₂ bubbles are found quite different than that of N₂ bubbles. The plausible reasons are that gas fraction distribution at the sparger region may inhibit bubble coalescence and the positive surface charge of CO₂ bubbles acts as a barrier to the interface deformation.

Introduction: Wall shear stress (WSS) is associated with the growth and rupture of an intracranial aneurysm. To reveal their underlying connections, many image-based computational fluid dynamics (CFD) studies have been conducted. However, the methodological validations using both in vivo medical imaging and in vitro optical flow measurements were rarely accompanied in such studies. Methods: In the present study, we performed a comparative assessment on the hemodynamics of a patient-specific intracranial saccular aneurysm using in vivo 4D Flow MRI, in silico CFD, in vitro stereoscopic and tomographic particle imaging velocimetry (Stereo-PIV and Tomo-PIV) techniques. PIV experiments and CFD were conducted under steady state corresponding to the peak systole of 4D Flow MRI. Results: The results showed that all modalities provided similar flow features and overall surface distribution of WSS. However, a large variation in the absolute WSS values was found. 4D Flow MRI estimated a 2- to 4-fold lower peak WSS (3.99 Pa) and a 1.6- to 2-fold lower mean WSS (0.94 Pa) than Tomo-PIV, Stereo-PIV, and CFD. Bland-Altman plots of WSS showed that the differences between PIV-/CFD-based WSS and 4D Flow MRI-based WSS increase with higher WSS magnitude. Such proportional trend was absent in the Bland-Altman comparison of velocity where the resolutions of PIV and CFD datasets were matched to 4D Flow MRI. We also found that because of superior resolution in the out-of-plane direction, WSS estimation by Tomo-PIV was higher than Stereo-PIV. Conclusions: Our results indicated that the differences in spatial resolution could be the main contributor to the discrepancies between each modality. The findings of this study suggest that with current techniques, care should be taken when using absolute WSS values to perform a quantitative risk analysis of aneurysm rupture.

Left ventricular (LV) blood flow is an inherently complex time-varying 3-D phenomenon, where 2-D quantification often ignores the effect of out-of-plane motion. In this study, we describe high frame rate 4-D echocardiographic particle image velocimetry (echo-PIV) using a prototype matrix transesophageal transducer and a dynamic LV phantom for testing the accuracy of echo-PIV in the presence of complex flow patterns. Optical time-resolved tomographic PIV (tomo-PIV) was used as a reference standard for comparison. Echo-PIV and tomo-PIV agreed on the general profile of the LV flow patterns, but echo-PIV smoothed out the smaller flow structures. Echo-PIV also underestimated the flow rates at greater imaging depths, where the PIV kernel size and transducer point spread function were large relative to the velocity gradients. We demonstrate that 4-D echo-PIV could be performed in just four heart cycles, which would require only a short breath-hold, providing promising results. However, methods for resolving high velocity gradients in regions of poor spatial resolution are required before clinical translation.

Left ventricular flow is intrinsically complex, three-dimensional and unsteady. Its features are susceptible to cardiovascular pathology and treatment, in particular to surgical interventions involving the valves (mitral valve replacement). To improve our understanding of intraventricular fluid mechanics and the impact of various types of prosthetic valves thereon, we have developed a custom-designed versatile left ventricular phantom with anatomically realistic moving left ventricular membrane. A biological, a tilting disc and a bileaflet valve (in two different orientations) were mounted in the mitral position and tested under the same settings. To investigate 3D flow within the phantom, a four-view tomographic particle image velocimetry setup has been implemented. The results compare side-by-side the evolution of the 3D flow topology, vortical structures and kinetic energy in the left ventricle domain during the cardiac cycle. Except for the tilting disc valve, all tested prosthetic valves induced a crossed flow path, where the outflow crosses the inflow path, passing under the mitral valve. The biological valve shows a strong jet with a peak velocity about twice as high compared to all mechanical heart valves, which makes it easier to penetrate deeply into the cavity. Accordingly, the peak kinetic energy in the left ventricle in case of the biological valve is about four times higher than the mechanical heart valves. We conclude that the tomographic particle imaging velocimetry setup provides a useful ground truth measurement of flow features and allows a comparison of the effects of different valve types on left ventricular flow patterns.