In this work, we compute numerically breakage/coalescence rates and size distribution of surfactant-laden droplets in turbulent flow. We use direct numerical simulation of turbulence coupled with a two-order-parameter phase-field method to describe droplets and surfactant dynamics. We consider two different values of the surface tension (i.e. two values for the Weber number, $We$, the ratio between inertial and surface tension forces) and four types of surfactant (i.e. four values of the elasticity number, $\unicode[STIX]{x1D6FD}_{s}$, which defines the strength of the surfactant). Stretching, breakage and merging of droplet interfaces are controlled by the complex interplay among shear stresses, surface tension and surfactant distribution, which are deeply intertwined. Shear stresses deform the interface, changing the local curvature and thus surface tension forces, but also advect surfactant over the interface. In turn, local increases of surfactant concentration reduce surface tension, changing the interface deformability and producing tangential (Marangoni) stresses. Finally, the interface feeds back to the local shear stresses via the capillary stresses, and changes the local surfactant distribution as it deforms, breaks and merges. We find that Marangoni stresses have a major role in restoring a uniform surfactant distribution over the interface, contrasting, in particular, the action of shear stresses: this restoring effect is proportional to the elasticity number and is stronger for smaller droplets. We also find that lower surface tension (higher $We$ or higher $\unicode[STIX]{x1D6FD}_{s}$) increases the number of breakage events, as expected, but also the number of coalescence events, more unexpected. The increase of the number of coalescence events can be traced back to two main factors: the higher probability of inter-droplet collisions, favoured by the larger number of available droplets, and the decreased deformability of smaller droplets. Finally, we show that, for all investigated cases, the steady-state droplet size distribution is in good agreement with the $-10/3$ power-law scaling (Garrett et al., J. Phys. Oceanogr., vol. 30 (9), 2000, pp. 2163–2171), conforming to previous experimental observations (Deane & Stokes, Nature, vol. 418 (6900), 2002, p. 839) and numerical simulations (Skartlien et al., J. Chem. Phys., vol. 139 (17), 2013).

In this work, we propose and test the validity of a modified Phase Field Method (PFM), which was specifically developed for large scale simulations of turbulent flows with large and deformable surfactant-laden droplets. The time evolution of the phase field, φ, and of the surfactant concentration field, ψ, are obtained from two Cahn-Hilliard-like equations together with a two-order-parameter Time-Dependent Ginzburg-Landau (TDGL) free energy functional. The modifications introduced circumvent existing limitations of current approaches based on PFM and improve the well-posedness of the model. The effect of surfactant on surface tension is modeled via an Equation Of State (EOS), further improving the flexibility of the approach. This method can efficiently handle topological changes, i.e. breakup and coalescence, and describe adsorption/desorption of surfactant. The capabilities of the proposed approach are tested in this paper against previous experimental results on the effects of surfactant on the deformation of a single droplet and on the interactions between two droplets. Finally, to appreciate the performances of the model on a large scale complex simulation, a qualitative analysis of the behavior of surfactant-laden droplets in a turbulent channel flow is presented and discussed.

After the Spanish flu pandemic, it was apparent that airborne transmission was crucial to spreading virus contagion, and research responded by producing several fundamental works like the experiments of Duguid [J. P. Duguid, J. Hyg. 44, 6 (1946)] and the model of Wells [W. F. Wells, Am. J. Hyg. 20, 611–618 (1934)]. These seminal works have been pillars of past and current guidelines published by health organizations. However, in about one century, understanding of turbulent aerosol transport by jets and plumes has enormously progressed, and it is now time to use this body of developed knowledge. In this work, we use detailed experiments and accurate computationally intensive numerical simulations of droplet-laden turbulent puffs emitted during sneezes in a wide range of environmental conditions. We consider the same emission—number of drops, drop size distribution, and initial velocity—and we change environmental parameters such as temperature and humidity, and we observe strong variation in droplets’ evaporation or condensation in accordance with their local temperature and humidity microenvironment. We assume that 3% of the initial droplet volume is made of nonvolatile matter. Our systematic analysis confirms that droplets’ lifetime is always about one order of magnitude larger compared to previous predictions, in some cases up to 200 times. Finally, we have been able to produce original virus exposure maps, which can be a useful instrument for health scientists and practitioners to calibrate new guidelines to prevent short-range airborne disease transmission.