Search

We investigate drop break-up morphology, occurrence, time and size distribution, through large ensembles of high-fidelity direct-numerical simulations of drops in homogeneous isotropic turbulence, spanning a wide range of parameters in terms of the Weber number $We$ , viscosity ratio between the drop and the carrier flow $\mu _r=\mu _d/\mu _l$ , where d is the drop diameter, and Reynolds ($Re$) number. For $\mu _r \leq 20$ , we find a nearly constant critical $We$ , while it increases with $\mu _r$ (and $Re$) when $\mu _r > 20$ , and the transition can be described in terms of a drop Reynolds number. The break-up time is delayed when $\mu _r$ increases and is a function of distance to criticality. The first break-up child-size distributions for $\mu _r \leq 20$ transition from M to U shape when the distance to criticality is increased. At high $\mu _r$ , the shape of the distribution is modified. The first break-up child-size distribution gives only limited information on the fragmentation dynamics, as the subsequent break-up sequence is controlled by the drop geometry and viscosity. At high $We$ , a $d^{-3/2}$ size distribution is observed for $\mu _r \leq 20$ , which can be explained by capillary-driven processes, while for $\mu _r > 20$ , almost all drops formed by the fragmentation process are at the smallest scale, controlled by the diameter of the very extended filament, which exhibits a snake-like shape prior to break-up. [ABSTRACT FROM AUTHOR]

We investigate the statistics of turbulence in emulsions of two immiscible fluids of the same density. We compute velocity increments between points conditioned to be located in the same phase or in different phases, and examine their probability density functions (PDFs) and the associated structure functions (SFs). This enables us to demonstrate that the presence of the interface reduces the skewness of the PDF at small scales and therefore the magnitude of the energy flux towards the dissipative scales, which is quantified by the third-order SF. The analysis of the higher-order SFs shows that multiphase turbulence is more intermittent than single-phase turbulence. In particular, the local scaling exponents of the SFs display a saturation below the Kolmogorov–Hinze scale, which indicates the presence of large velocity gradients across the interface. Interestingly, the statistics of the velocity differences in the carrier phase recovers that of single-phase turbulence when the viscosity of the dispersed phase is high. [ABSTRACT FROM AUTHOR]

The dynamics of droplet fragmentation in turbulence is described by the Kolmogorov-Hinze framework. Yet, a quantitative theory is lacking at higher concentrations when strong interactions between the phases and coalescence become relevant, which is common in most flows. Here, we address this issue through a fully-coupled numerical study of the droplet dynamics in a turbulent flow at Rλ ≈ 140, the highest attained up to now. By means of time-space spectral statistics, not currently accessible to experiments, we demonstrate that the characteristic scale of the process, the Hinze scale, can be precisely identified as the scale at which the net energy exchange due to capillarity is zero. Droplets larger than this scale preferentially break up absorbing energy from the flow; smaller droplets, instead, undergo rapid oscillations and tend to coalesce releasing energy to the flow. Further, we link the droplet-size distribution with the probability distribution of the turbulent dissipation. This shows that key in the fragmentation process is the local flux of energy which dominates the process at large scales, vindicating its locality. Dynamics of droplet fragmentation in turbulence is described by the Kolmogorov-Hinze theory, but at higher concentrations common in most flows a quantitative theory is required. The authors use direct numerical simulations of turbulent multiphase flows finding that larger droplets break up absorbing energy from the flow, while smaller droplets undergo rapid oscillations and tend to coalesce releasing energy to the flow. [ABSTRACT FROM AUTHOR]

We present high-resolution three-dimensional (3-D) direct numerical simulations of breaking waves solving for the two-phase Navier-Stokes equations. We investigate the role of the Reynolds number (Re, wave inertia relative to viscous effects) and Bond number (Bo, wave scale over the capillary length) on the energy, bubble and droplet statistics of strong plunging breakers. We explore the asymptotic regimes at high Re and Bo, and compare with laboratory breaking waves. Energetically, the breaking wave transitions from laminar to 3-D turbulent flow on a time scale that depends on the turbulent Re up to a limiting value Reλ ∼ 100, consistent with the mixing transition in other canonical turbulent flows. We characterize the role of capillary effects on the impacting jet and ingested main cavity shape and subsequent fragmentation process, and extend the buoyant-energetic scaling from Deike et al. (J. Fluid Mech., vol. 801, 2016, pp. 91-129) to account for the cavity shape and its scale separation from the Hinze scale, rH. We confirm two regimes in the bubble size distribution, N(r/rH) ∝ (r/rH)-10/3 forr > rH, and ∝ r-2d. The evolution of the droplet statistics appears controlled by the details (r/rH)-3/2 for r < rH. Bubbles are resolved up to one order of magnitude below rH, and we observe a good collapse of the numerical data compared to laboratory breaking waves (Deane & Stokes, Nature, vol. 418 (6900), 2002, pp. 839-844). We resolve droplet statistics at high Bo in good agreement with recent experiments (Erinin et al., Geophys. Res. Lett., vol. 46 (14), 2019, pp. 8244-8251), with a distribution shape close to Nd(rd) ∝ r-2d. The evolution of the droplet statistics appears controlled by the details of the impact process and subsequent splash-up. We discuss velocity distributions for the droplets, finding ejection velocities up to four times the phase speed of the wave, which are produced during the most intense splashing events of the breaking process. [ABSTRACT FROM AUTHOR]

We present a numerical study of emulsions in homogeneous and isotropic turbulence (HIT) at Reλ = 137. The problem is addressed via direct numerical simulations, where the volume of fluid is used to represent the complex features of the liquid–liquid interface. We consider a mixture of two iso-density fluids, where fluid properties are varied with the goal of understanding their role in turbulence modulation, in particular the volume fraction (0.03 < α < 0.5), viscosity ratio (0.01 < μd/μc < 100) and large-scale Weber number (10.6 < WeL < 106.5). The analysis, performed by studying integral quantities and spectral scale-by-scale analysis, reveals that energy is transported consistently from large to small scales by the interface, and no inverse cascade is observed. Furthermore, the total surface is found to be directly proportional to the amount of energy transported, while viscosity and surface tension alter the dynamic that regulates energy transport. We also observe the −10/3 and −3/2 scaling on droplet size distributions, suggesting that the dimensional arguments that led to their derivation are verified in HIT conditions. [ABSTRACT FROM AUTHOR]

We investigate the modes of deformation of an initially spherical bubble immersed in a homogeneous and isotropic turbulent background flow. We perform direct numerical simulations of the two-phase incompressible Navier–Stokes equations, considering a low-density bubble in the high-density turbulent flow at various Weber numbers (the ratio of turbulent and surface tension forces) using the air–water density ratio. We discuss a theoretical framework for the bubble deformation in a turbulent flow using a spherical harmonic decomposition. We propose, for each mode of bubble deformation, a forcing term given by the statistics of velocity and pressure fluctuations, evaluated on a sphere of the same radius. This approach formally relates the bubble deformation and the background turbulent velocity fluctuations, in the limit of small deformations. The growth of the total surface deformation and of each individual mode is computed from the direct numerical simulations using an appropriate Voronoi decomposition of the bubble surface. We show that two successive temporal regimes occur: the first regime corresponds to deformations driven only by inertial forces, with the interface deformation growing linearly in time, in agreement with the model predictions, whereas the second regime results from a balance between inertial forces and surface tension. The transition time between the two regimes is given by the period of the first Rayleigh mode of bubble oscillation. We discuss how our approach can be used to relate the bubble lifetime to the turbulence statistics and eventually show that at high Weber numbers, bubble lifetime can be deduced from the statistics of turbulent fluctuations at the bubble scale. [ABSTRACT FROM AUTHOR]

We study bubble break-up in homogeneous and isotropic turbulence by direct numerical simulations of the two-phase incompressible Navier–Stokes equations. We create the turbulence by forcing in physical space and introduce the bubble once a statistically stationary state is reached. We perform a large ensemble of simulations to investigate the effect of the Weber number (the ratio of turbulent and surface tension forces) on bubble break-up dynamics and statistics, including the child bubble size distribution, and discuss the numerical requirements to obtain results independent of grid size. We characterize the critical Weber number below which no break-up occurs and the associated Hinze scale dh. At Weber number close to stable conditions (initial bubble sizes d0≈dh), we observe binary and tertiary break-ups, leading to bubbles mostly between 0.5dh and dh, a signature of a production process local in scale. For large Weber numbers (d0>3dh), we observe the creation of a wide range of bubble radii, with numerous child bubbles between 0.1dh and 0.3dh, an order of magnitude smaller than the parent bubble. The separation of scales between the parent and child bubble is a signature of a production process non-local in scale. The formation mechanism of these sub-Hinze scale bubbles relates to rapid large deformation and successive break-ups: the first break-up in a sequence leaves highly deformed bubbles which will break again, without recovering a spherical shape and creating an array of much smaller bubbles. We discuss the application of this scenario to the production of sub-Hinze bubbles under breaking waves. [ABSTRACT FROM AUTHOR]