We study turbulent emulsions and the emulsification process in homogeneous isotropic turbulence (HIT) using direct numerical simulations (DNS) in combination with the volume of fluid method (VOF). For generating a turbulent flow field, we employ a linear forcing approach augmented by a proportional‐integral‐derivative (PID) controller, which ensures a constant turbulent kinetic energy for two phase flow scenarios and accelerates the emulsification process. For the simulations, the density ratio of dispersed and carrier phases is chosen to be similar to that of oil and water (0.9), representing a typical application. We vary the turbulence intensity and the surface tension coefficient. Thus, we modulate those parameters that directly affect the Hinze scale, which is expected to be the most stable maximum droplet diameter in emulsions in HIT. The considered configurations can be characterized with Taylor Reynolds numbers in the range of 100–140 and Weber numbers, evaluated with the velocity fluctuations and the integral length scale, of 4–70. Using the 3‐D simulation results, we study the emulsification process as well as the emulsions at a statistically stationary state. For the latter, droplet size distributions are evaluated and compared. We observe a Hinze scale similarity of the size distributions considering a fixed integral length scale, that is, similar Hinze scales obtained at different turbulence intensities or for different fluid properties result in similar distributions.

We study the segregation of emulsions in decaying turbulence using direct numerical simulations in combination with the volume of fluid method. To this end, we generate emulsions in forced homogeneous isotropic turbulence and then turn the forcing off and activate the gravitational acceleration. This allows us to study the segregation process in decaying turbulence and under gravity. We consider non-iso-density emulsions, where the dispersed phase is the lighter one. The segregation process is driven by both the minimization of the potential energy achieved by the sinking of the heavier phase as well as the minimization of the surface energy achieved by coalescence. To study these two processes and their impacts on the segregation progress in detail, we consider different buoyancy forces and surface tension coefficients in our investigation, resulting in five different configurations. The surface tension coefficient also alters the droplet size distribution of the emulsion. Using the three-dimensional simulation results and the monitored data, we analyze the driving mechanisms and their impact on the segregation progress in detail. We propose a dimensionless number that reflects the energy release dominating the segregation. Moreover, we evaluate the time required for the rise of the lighter phase and study correlations with the varied parameters: gravitational acceleration and surface tension coefficient.

The influence of non-Newtonian fluid behavior on the flow statistics of turbulent bubble-laden downflow in a vertical channel is investigated. A Direct Numerical Simulation (DNS) study is conducted for power-law fluids with power-law indexes of 0.7 (shear-thinning), 1 (Newtonian) and 1.3 (shear-thickening) in the liquid phase at a gas volume fraction of 6%. The flow is driven downward by a constant volumetric flow rate corresponding to a friction Reynolds number of Reτ≈127.3. The Eötvös number is varied between Eo=0.3125 and Eo=3.75 in order to investigate the influence of quasi-spherical as well as wobbling bubbles and thus the interplay of the bubble deformability with the power-law behavior of the liquid bulk. The resulting first- and second-order fluid statistics, i.e., the gas fraction, mean velocity and velocity fluctuation profiles across the channel, show clear trends in reply to varying power-law indexes. In addition, it was observed that the bubble oscillations increase with decreasing power-law index. In the channel core, the bubbles significantly increase the dissipation rate, which, in contrast to its behavior at the wall, shows similar orders of magnitude for all power-law indexes.

The influence of non-Newtonian fluid behavior and the Eötvös number on conditional and unconditional second-order structure functions of bubbly channel flows has been investigated by conducting a series of direct numerical simulations at a friction Reynolds number of 127.3. Two Eötvös numbers have been considered (Eo = 0.3125 and Eo = 3.75) together with three different power-law indexes representing shear-thinning (n = 0.7), Newtonian (n = 1.0), and shear-thickening (n = 1.3) fluid behavior. The scaling of the second-order structure functions (SFs) can be translated into an inertial range scaling of the turbulent kinetic energy spectrum. However, because of the discontinuous character of the fluid properties in bubbly flows, SFs are more easily accessible than turbulence spectra, which are based on Fourier transform. It has been found that the different parameters (i.e., Eo, n) have an influence on the energy content as well as the peak location of the compensated second-order SFs (i.e., the dimensions of the large scales). However, after appropriate scaling, the curves nearly collapse. To confirm and further explain the above findings, directional length scales have been evaluated and discussed in detail. Finally, the anisotropy of the Reynolds stress tensor and dissipation tensor has been analyzed in terms of the Lumley triangle, showing that bubbly channel flows are less isotropic than their single-phase counterpart, although they are more homogeneous in the channel center. While the dissipation tensor is slightly more isotropic than the Reynolds stress tensor in the bulk region of the channel flow, overall, a very similar behavior is observed.

No abstract