To quantify the fate of respiratory droplets under different ambient relative humidities, direct numerical simulations of a typical respiratory event are performed. We found that, because small ...droplets (with initial diameter of 10 μm) are swept by turbulent eddies in the expelled humid puff, their lifetime gets extended by a factor of more than 30 times as compared to what is suggested by the classical picture by Wells, for 50% relative humidity. With increasing ambient relative humidity the extension of the lifetimes of the small droplets further increases and goes up to around 150 times for 90% relative humidity, implying more than 2 m advection range of the respiratory droplets within 1 sec. Employing Lagrangian statistics, we demonstrate that the turbulent humid respiratory puff engulfs the small droplets, leading to many orders of magnitude increase in their lifetimes, implying that they can be transported much further during the respiratory events than the large ones. Our findings provide the starting points for larger parameter studies and may be instructive for developing strategies on optimizing ventilation and indoor humidity control. Such strategies are key in mitigating the COVID-19 pandemic in the present autumn and upcoming winter.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UL, UM
We theoretically and numerically investigate the instabilities driven by diffusiophoretic flow, caused by a solutal concentration gradient along a reacting surface. The important control parameters ...are the Péclet number $Pe$, which quantifies the ratio of the solutal advection rate to the diffusion rate, and the Schmidt number $Sc$, which is the ratio of viscosity and diffusivity. First, we study the diffusiophoretic flow on a catalytic plane in two dimensions. From a linear stability analysis, we obtain that for $Pe$ larger than $8{\rm \pi}$ mass transport by convection overtakes that by diffusion, and a symmetry-breaking mode arises, which is consistent with numerical results. For even larger $Pe$, nonlinear terms become important. For $Pe > 16{\rm \pi}$, multiple concentration plumes are emitted from the catalytic plane, which eventually merge into a single larger one. When $Pe$ is even larger ($Pe \gtrsim 603$ for Schmidt number $Sc=1$), there are continuous emissions and merging events of the concentration plumes. The newly found flow states have different flow structures for different $Sc$: for $Sc\geqslant 1$, we observe the chaotic emission of plumes, but the fluctuations of concentration are only present in the region near the catalytic plane. In contrast, for $Sc<1$, chaotic flow motion occurs also in the bulk. In the second part of the paper, we conduct three-dimensional simulations for spherical catalytic particles, and beyond a critical Péclet number again find continuous plume emission and plume merging, now leading to a chaotic motion of the phoretic particle. Our results thus help us to understand the experimentally observed chaotic motion of catalytic particles in the high $Pe$ regime.
We numerically investigate both single and multiple droplet dissolution with droplets consisting of less dense liquid dissolving in a denser host liquid. In this situation, buoyancy can lead to ...convection and thus plays an important role in the dissolution process. The significance of buoyancy is quantified by the Rayleigh number
$Ra$
, which is the buoyancy force over the viscous damping force. In this study,
$Ra$
spans almost four decades from 0.1 to 400. We focus on how the mass flux, characterized by the Sherwood number
$Sh$
, and the flow morphologies depend on
$Ra$
. For single droplet dissolution, we first show the transition of the
$Sh(Ra)$
scaling from a constant value to
$Sh\sim Ra^{1/4}$
, which confirms the experimental results by Dietrich
et al.
(
J. Fluid Mech.
, vol. 794, 2016, pp. 45–67). The two distinct regimes, namely the diffusively and the convectively dominated regimes, exhibit different flow morphologies: when
$Ra\geqslant 10$
, a buoyant plume is clearly visible, which contrasts sharply with the pure diffusion case at low
$Ra$
. For multiple droplet dissolution, the well-known shielding effect comes into play at low
$Ra$
, so that the dissolution rate is slower as compared to the single droplet case. However, at high
$Ra$
, convection becomes more and more dominant so that a collective plume enhances the mass flux, and remarkably the multiple droplets dissolve faster than a single droplet. This has also been found in the experiments by Laghezza
et al.
(
Soft Matt.
, vol. 12 (26), 2016, pp. 5787–5796). We explain this enhancement by the formation of a single, larger plume rather than several individual plumes. Moreover, there is an optimal
$Ra$
at which the enhancement is maximized, because the single plume is narrower at larger
$Ra$
, which thus hinders the enhancement. Our findings demonstrate a new mechanism in collective droplet dissolution, which is the merging of the plumes, which leads to non-trivial phenomena, contrasting the shielding effect.
Indoor ventilation is essential for a healthy and comfortable living environment. A key issue is to discharge anthropogenic air contamination such as CO$_2$ gas or, of potentially more direct ...consequence, airborne respiratory droplets. Here, by employing direct numerical simulations, we study mechanical displacement ventilation with a wide range of ventilation rates $Q$ from 0.01 to 0.1 m$^3$ s$^{-1}$ person$^{-1}$. For this ventilation scheme, a cool lower zone is established beneath a warm upper zone with interface height $h$, which depends on $Q$. For weak ventilation, we find the scaling relation $h\sim Q^{3/5}$, as suggested by Hunt & Linden (Build. Environ., vol. 34, 1999, pp. 707–720). Also, the CO$_{2}$ concentration decreases with $Q$ within this regime. However, for too strong ventilation, the interface height $h$ becomes insensitive to $Q$, and the ambient averaged CO$_2$ concentration decreases towards the ambient value. At these values of $Q$, the concentrations of pollutants are very low and so further dilution has little effect. We suggest that such scenarios arise when the vertical kinetic energy associated with the ventilation flow is significant compared with the potential energy of the thermal stratification.
Rayleigh–Bénard (RB) convection with free-slip plates and horizontally periodic boundary conditions is investigated using direct numerical simulations. Two configurations are considered, one is ...two-dimensional (2-D) RB convection and the other one three-dimensional (3-D) RB convection with a rotating axis parallel to the plate, which for strong rotation mimics 2-D RB convection. For the 2-D simulations, we explore the parameter range of Rayleigh numbers $Ra$ from $10^{7}$ to $10^{9}$ and Prandtl numbers $Pr$ from $1$ to $100$. The effect of the width-to-height aspect ratio $\varGamma$ is investigated for $1\leqslant \varGamma \leqslant 128$. We show that zonal flow, which was observed, for example, by Goluskin et al. (J. Fluid. Mech., vol. 759, 2014, pp. 360–385) for $\varGamma =2$, is only stable when $\varGamma$ is smaller than a critical value, which depends on $Ra$ and $Pr$. The regime in which only zonal flow can exist is called the first regime in this study. With increasing $\varGamma$, we find a second regime in which both zonal flow and different convection roll states can be statistically stable. For even larger $\varGamma$, in a third regime, only convection roll states are statistically stable and zonal flow is not sustained. How many convection rolls form (or in other words, what the mean aspect ratio of an individual roll is), depends on the initial conditions and on $Ra$ and $Pr$. For instance, for $Ra=10^{8}$ and $Pr=10$, the aspect ratio $\varGamma _r$ of an individual, statistically stable convection roll can vary in a large range between $16/11$ and $64$. A convection roll with a large aspect ratio of $\varGamma _r = 64$, or more generally already with $\varGamma _r \gg 10$, can be seen as ‘localized’ zonal flow, and indeed carries over various properties of the global zonal flow. For the 3-D simulations, we fix $Ra=10^{7}$ and $Pr=0.71$, and compare the flow for $\varGamma =8$ and $\varGamma = 16$. We first show that with increasing rotation rate both the flow structures and global quantities like the Nusselt number $Nu$ and the Reynolds number $Re$ increasingly behave like in the 2-D case. We then demonstrate that with increasing aspect ratio $\varGamma$, zonal flow, which was observed for small $\varGamma =2{\rm \pi}$ by von Hardenberg et al. (Phys. Rev. Lett., vol. 15, 2015, 134501), completely disappears for $\varGamma =16$. For such large $\varGamma$, only convection roll states are statistically stable. In-between, here for medium aspect ratio $\varGamma = 8$, the convection roll state and the zonal flow state are both statistically stable. What state is taken depends on the initial conditions, similarly as we found for the 2-D case.
For dissolving active oil droplets in an ambient liquid, it is generally assumed that the Marangoni effect results in repulsive interactions, while the buoyancy effects caused by the density ...difference between the droplets, diffusing product and the ambient fluid are usually neglected. However, it has been observed in recent experiments that active droplets can form clusters due to buoyancy-driven convection (Krüger et al., Eur. Phys. J. E, vol. 39, 2016, pp. 1–9). In this study we numerically analyse the buoyancy effect, in addition to the propulsion caused by Marangoni flow (with its strength characterized by the Péclet number $Pe$). The buoyancy effects have their origin in (i) the density difference between the droplet and the ambient liquid, which is characterized by the Galileo number $Ga$; and (ii) the density difference between the diffusing product (i.e. filled micelles) and the ambient liquid, which can be quantified by a solutal Rayleigh number $Ra$. We analyse how the attracting and repulsing behaviour of neighbouring droplets depends on the control parameters $Pe$, $Ga$ and $Ra$. We find that while the Marangoni effect leads to the well-known repulsion between the interacting droplets, the buoyancy effect of the reaction product leads to buoyancy-driven attraction. At sufficiently large $Ra$, even collisions between the droplets can take place. Our study on the effect of $Ga$ further shows that with increasing $Ga$, the collision becomes delayed. Moreover, we derive that the attracting velocity of the droplets, which is characterized by a Reynolds number $Re_d$, is proportional to $Ra^{1/4}/( \ell /R)$, where $\ell /R$ is the distance between the neighbouring droplets normalized by the droplet radius. Finally, we numerically obtain the repulsive velocity of the droplets, characterized by a Reynolds number $Re_{rep}$, which is proportional to $PeRa^{-0.38}$. The balance of attractive and repulsive effect leads to $Pe\sim Ra^{0.63}$, which agrees well with the transition curve between the regimes with and without collision.
This numerical study presents a simple but extremely effective way to considerably enhance heat transport in turbulent wall-bounded multiphase flows, namely by using oleophilic walls. As a model ...system, we pick the Rayleigh–Bénard set-up, filled with an oil–water mixture. For oleophilic walls, using only $10\,\%$ volume fraction of oil in water, we observe a remarkable heat transport enhancement of more than $100\,\%$ as compared to the pure water case. In contrast, for oleophobic walls, the enhancement is only of about $20\,\%$ as compared to pure water. The physical explanation of the heat transport increment for oleophilic walls is that thermal plumes detach from the oil-rich boundary layer and carry the heat with them. In the bulk, the oil–water interface prevents the plumes from mixing with the turbulent water bulk and to diffuse their heat. To confirm this physical picture, we show that the minimum amount of oil necessary to achieve the maximum heat transport is set by the volume fraction of the thermal plumes. Our findings provide guidelines of how to optimize heat transport in wall-bounded thermal turbulence. Moreover, the physical insight of how coherent structures are coupled with one of the phases of a two-phase system has very general applicability for controlling transport properties in other turbulent wall-bounded multiphase flows.
In the preparation of café latte, spectacular layer formation can occur between the espresso shot in a glass of milk and the milk itself. Xue et al. (Nat. Commun., vol. 8, 2017, pp. 1–6) showed that ...the injection velocity of espresso determines the depth of coffee–milk mixture. After a while, when a stable stratification forms in the mixture, the layering process can be modelled as a double diffusive convection system with a stably stratified coffee–milk mixture cooled from the side. More specifically, we perform (two-dimensional) direct numerical simulations of laterally cooled double diffusive convection for a wide parameter range, where the convective flow is driven by a lateral temperature gradient while stabilized by a vertical concentration gradient. Depending on the strength of stabilization as compared to the thermal driving, the system exhibits different flow regimes. When the thermal driving force dominates over the stabilizing force, the flow behaves like vertical convection in which a large-scale circulation develops. However, with increasing strength of the stabilizing force, a meta-stable layered regime emerges. Initially, several vertically-stacked convection rolls develop, and these well-mixed layers are separated by sharp interfaces with large concentration gradients. The initial thickness of these emerging layers can be estimated by balancing the work exerted by thermal driving and the required potential energy to bring fluid out of its equilibrium position in the stably stratified fluid. In the layered regime, we further observe successive layer merging, and eventually only a single convection roll remains. We elucidate the following merging mechanism: as weakened circulation leads to accumulation of hot fluid adjacent to the hot sidewall, larger buoyancy forces associated with hotter fluid eventually break the layer interface. Then two layers merge into a larger layer, and circulation establishes again within the merged structure.
We numerically investigate turbulent Rayleigh–Bénard convection through two immiscible fluid layers, aiming to understand how the layer thickness and fluid properties affect the heat transfer ...(characterized by the Nusselt number $\mbox {Nu}$) in two-layer systems. Both two- and three-dimensional simulations are performed at fixed global Rayleigh number $\mbox {Ra}=10^8$, Prandtl number $\mbox {Pr}=4.38$ and Weber number $\mbox {We}=5$. We vary the relative thickness of the upper layer between $0.01 \le \alpha \le 0.99$ and the thermal conductivity coefficient ratio of the two liquids between $0.1 \le \lambda _k \le 10$. Two flow regimes are observed. In the first regime at $0.04\le \alpha \le 0.96$, convective flows appear in both layers and $\mbox {Nu}$ is not sensitive to $\alpha$. In the second regime at $\alpha \le 0.02$ or $\alpha \ge 0.98$, convective flow only exists in the thicker layer, while the thinner one is dominated by pure conduction. In this regime, $\mbox {Nu}$ is sensitive to $\alpha$. To predict $\mbox {Nu}$ in the system in which the two layers are separated by a unique interface, we apply the Grossmann–Lohse theory for both individual layers and impose heat flux conservation at the interface. Without introducing any free parameter, the predictions for $\mbox {Nu}$ and for the temperature at the interface agree well with our numerical results and previous experimental data.
We perform a numerical study of the heat transfer and flow structure of Rayleigh–Bénard (RB) convection in (in most cases regular) porous media, which are comprised of circular, solid obstacles ...located on a square lattice. This study is focused on the role of porosity
$\unicodeSTIX{x1D719}$
in the flow properties during the transition process from the traditional RB convection with
$\unicodeSTIX{x1D719}=1$
(so no obstacles included) to Darcy-type porous-media convection with
$\unicodeSTIX{x1D719}$
approaching 0. Simulations are carried out in a cell with unity aspect ratio, for Rayleigh number
$Ra$
from
$10^{5}$
to
$10^{10}$
and varying porosities
$\unicodeSTIX{x1D719}$
, at a fixed Prandtl number
$Pr=4.3$
, and we restrict ourselves to the two-dimensional case. For fixed
$Ra$
, the Nusselt number
$Nu$
is found to vary non-monotonically as a function of
$\unicodeSTIX{x1D719}$
; namely, with decreasing
$\unicodeSTIX{x1D719}$
, it first increases, before it decreases for
$\unicodeSTIX{x1D719}$
approaching 0. The non-monotonic behaviour of
$Nu(\unicodeSTIX{x1D719})$
originates from two competing effects of the porous structure on the heat transfer. On the one hand, the flow coherence is enhanced in the porous media, which is beneficial for the heat transfer. On the other hand, the convection is slowed down by the enhanced resistance due to the porous structure, leading to heat transfer reduction. For fixed
$\unicodeSTIX{x1D719}$
, depending on
$Ra$
, two different heat transfer regimes are identified, with different effective power-law behaviours of
$Nu$
versus
$Ra$
, namely a steep one for low
$Ra$
when viscosity dominates, and the standard classical one for large
$Ra$
. The scaling crossover occurs when the thermal boundary layer thickness and the pore scale are comparable. The influences of the porous structure on the temperature and velocity fluctuations, convective heat flux and energy dissipation rates are analysed, further demonstrating the competing effects of the porous structure to enhance or reduce the heat transfer.