Concentrated emulsions flowing through channels of varying widths are omnipresent in daily life, from dispensing mayonnaise in our kitchens to large-scale industrial processing of food, pharmaceuticals, etc. Local changes in channel geometry affect the stability of emulsions over length scales far beyond the droplet magnitude, for example through propagation of coalescence events called a coalescence avalanche. The underlying mechanisms are not well understood. In this work, we investigated the stability of concentrated emulsions flowing through microchannels featuring a constriction. We found that in this model geometry, the acceleration of the droplets induced near the entrance of the constriction triggers a coalescence event between the leading and the trailing droplet, but only above a critical droplet velocity. This separation-induced coalescence event, in turn, was found to trigger a coalescence avalanche in the upstream direction. Analysis of the flow behavior through particle image velocimetry and particle tracking velocimetry revealed that the propagation also follows a separation-induced coalescence mechanism, due to the retraction of the interface of the trailing droplet upon coalescence and the corresponding acceleration of the liquid inside the coalesced fluid thread. The constriction ratio was found to enhance the coalescence occurrence but did not affect the speed of coalescence propagation.

This study reports the distribution of gas and liquid velocity fluctuations and flow structures in N₂-water and CO₂-water bubbly flows within a pseudo-2D bubble column reactor (BCR). PIV and BIV techniques, operating at a resolution of 750 Hz, were synchronized, and ensemble averaging of the recorded velocity fields was performed to mitigate distortions caused by light interference, reflections, and shadows. Probability density functions (PDF), power spectral density (PSD), and cross-correlation values of velocity fluctuations were calculated to analyze gas bubble interaction and dissolution in dilute to dense bubbly flow regimes. The singular value decomposition (SVD) technique was applied to correlate the flow structure distribution and energy content of velocity fluctuation. The results indicate that CO₂-water bubbly flow exhibits a bimodal distribution of velocity fluctuation at low gas superficial velocity (U_g = 5.25 ± 1e-4 mm/s) due to bubble dissolution, whereas N₂-water bubbly flow at U_g = 8.67 ± 0.035 mm/s shows unimodal distribution. At high gas superficial velocities (U_g> 20.91 ± 0.075 mm/s for N₂ and U_g> 15.75 ± 0.05 mm/s for CO₂), bubbly flows exhibit unimodal distribution of velocity fluctuations. The bimodal distribution is attributed to bubble coalescence and breakup. Experimental findings also suggest that recirculation zones near the BCR walls primarily contain small flow structures with the highest energy density and rapid structural decay. In contrast, larger flow structures are found in the center of the BCR, exhibiting slower decay. Detailed insights into CO₂ bubble dissolution could enhance BCR design, thereby improving the efficiency of relevant industrial applications.

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.