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.

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 $R{e}_{\tau }\approx 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.

Using direct numerical simulation (DNS) in combination with the volume of fluid method (VoF), we investigate the influence of the density ratio between the carrier and dispersed phase on emulsions, where the baseline simulation approximately corresponds to the ratio of water-in-gasoline emulsions. For this purpose, homogeneous isotropic turbulence (HIT) is generated using a linear forcing method, enhanced by a proportional–integral–derivative (PID) controller, ensuring a constant turbulent kinetic energy (TKE) for two-phase flows, where the TKE balance equation contains an additional term due to surface tension. Then, the forcing is stopped, and gravitational acceleration is activated. The proposed computational setup represents a unique and well-controlled configuration to study emulsification and segregation. We consider four different density ratios, which are applied in industrial processes, to investigate the influence of the density ratio on the statistically steady state of the emulsions, and their segregation under decaying turbulence and constant gravitational acceleration. At the statistically steady state, we hold the turbulence constant and study the effects of the density ratio ${\rho }_d/{\rho }_c$, on the interface area, the Sauter mean diameter (SMD), and the statistical droplet size distribution. We find that all are affected by the density ratio, and we observe a relation between the SMD and ${\rho }_d/{\rho }_c$. Furthermore, we assume a dependence of the critical Weber number on the density ratio. In the second part of our work, we study the segregation process. To this end, we consider the change in the center of mass of the disperse phase and the energy release, to analyze the dependence of segregation on the density difference $\Delta \rho /{\rho }_d$. We show that segregation scales with the density difference and the droplet size, and a segregation time scale has been suggested that collapses the height of the center of mass for different density ratios.

This study is concerned with Large Eddy Simulation of liquid jet atomization using a two-way coupled Eulerian-Lagrangian multiscale approach. The proposed framework combines Volume-of-Fluid interface capturing with Lagrangian Particle Tracking. The former is used to compute the core jet and large liquid elements in the near-nozzle region, whereas the latter is used to track the large number of small droplets in the dilute downstream region of the spray. The convective and surface tension sub-grid scale terms arising in the context of two-phase flow LES are closed using suitable models, and secondary atomization is considered by employing a modified version of the Taylor Analogy Breakup model. The introduced framework is used to simulate an oil-in-air atomization as well as the Diesel-like Spray A test case of the Engine Combustion Network. Compared to previous studies based on Eulerian-Lagrangian methods, the present work stands out for the high-fidelity numerical approach, the complex test cases and the detailed comparison of the results to experimental data, which indicates a promising performance.