No Arabic abstract
Building on the recent theoretical work of Wray, Duffy and Wilson [J. Fluid Mech. 884, A45 (2020)] concerning the competitive diffusion-limited evaporation of multiple thin sessile droplets in proximity to each other, we obtain theoretical predictions for the spatially non-uniform densities of the contact-line deposits (often referred to as coffee stains or ring stains) left on the substrate after such droplets containing suspended solid particles have completely evaporated. Neighbouring droplets interact via their vapour fields, which results in a spatially non-uniform shielding effect. We give predictions for the deposits from a pair of identical droplets, which show that the deposit is reduced the most where the droplets are closest together, and demonstrate excellent quantitative agreement with experimental results of Pradhan and Panigrahi [Coll. Surf. A 482, 562-567 (2015)]. We also give corresponding predictions for a triplet of identical droplets arranged in an equilateral triangle, which show that the effect of shielding on the deposit is more subtle in this case.
We numerically investigate both single and multiple droplet dissolution with droplets consisting of lighter liquid dissolving in a denser host liquid. 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 $Shsim 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 regime, exhibit different flow morphologies: when Ra>=10, a buoyant plume is clearly visible which contrasts sharply to 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 Matter, vol. 12, 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, that leads to non-trivial phenomena, contrasting the shielding effect.
One-dimensional numerical simulations based on hybrid Eulerian-Lagrangian method are performed to study the interactions between propagating shocks and dispersed evaporating water droplets. Two-way coupling for exchanges of mass, momentum, energy and vapour species is adopted for the dilute two-phase gas-droplet flows. Interphase interactions and droplet breakup dynamics are investigated with initial droplet diameters of 30, 50, 70 and 90 {mu}m under an incident shock wave Mach number of 1.3. Novel two-phase flow phenomena are observed when droplet breakup occurs. First, droplets near the two-phase contact surface show obvious dispersed distribution because of the reflected pressure wave that propagates in the reverse direction of the leading shock. The reflected pressure wave grows stronger for larger droplets. Second, spatial oscillations of the gas phase pressure, droplet quantities (e.g., diameter and net force) and two-phase interactions (e.g., mass, momentum, and energy exchange), are observed in the post-shock region when droplet breakup occurs, which are caused by shock / droplet interactions. Third, the spatial distribution of droplets (i.e., number density, volume fraction) also shows strong oscillation in the post-shock region when droplet breakup occurs, which is caused by the oscillating force exerted on the droplets.
One-dimensional numerical simulations based on hybrid Eulerian-Lagrangian approach are performed to investigate the interactions between propagating shock waves and dispersed evaporating water droplets in two-phase gas-droplet flows. Two-way coupling for interphase exchanges of mass, momentum and energy is adopted. Parametric study on shock attenuation, droplet evaporation, motion and heating is conducted, through considering various initial droplet diameters (5-20 {mu}m), number densities (2.5 x 1011 - 2 x 1012 1/m3) and incident shock Mach numbers (1.17-1.9). It is found that the leading shock may be attenuated to sonic wave and even subsonic wave when droplet volume fraction is large and/or incident shock Mach number is low. Attenuation in both strength and propagation speed of the leading shock is mainly caused by momentum transfer to the droplets that interact at the shock front. Total pressure recovery is observed in the evaporation region, whereas pressure loss results from shock compression, droplet drag and pressure gradient force behind the shock front. Recompression of the region between the leading shock and two-phase contact surface is observed when the following compression wave is supersonic. After a critical point, this region gets stable in width and interphase exchanges in mass, momentum, and energy. However, the recompression phenomenon is sensitive to droplet volume fraction and may vanish with high droplet loading. For an incident shock Mach number of 1.6, recompression only occurs when the initial droplet volume fraction is below 3.28 x 10-5.
We consider self-propelled droplets which are driven by internal flow. Tracer particles, which are advected by the flow, in general follow chaotic trajectories, even though the motion of the autonomous swimmer is completely regular. The flow is mixing, and for P{e}clet and Batchelor numbers, which are realized e.g. in eucaryotic cells, advective mixing can substantially accelerate and even dominate transport by diffusion.
Active droplets swim as a result of the nonlinear advective coupling of the distribution of chemical species they consume or release with the Marangoni flows created by their non-uniform surface distribution. Most existing models focus on the self-propulsion of a single droplet in an unbounded fluid, which arises when diffusion is slow enough (i.e. beyond a critical Peclet number, $mbox{Pe}_c$). Despite its experimental relevance, the coupled dynamics of multiple droplets and/or collision with a wall remains mostly unexplored. Using a novel approach based on a moving fitted bispherical grid, the fully-coupled nonlinear dynamics of the chemical solute and flow fields are solved here to characterise in detail the axisymmetric collision of an active droplet with a rigid wall (or with a second droplet). The dynamics is strikingly different depending on the convective-to-diffusive transport ratio, $mbox{Pe}$: near the self-propulsion threshold (moderate $mbox{Pe}$), the rebound dynamics are set by chemical interactions and are well captured by asymptotic analysis; in contrast, for larger $mbox{Pe}$, a complex and nonlinear combination of hydrodynamic and chemical effects set the detailed dynamics, including a closer approach to the wall and a velocity plateau shortly after the rebound of the droplet. The rebound characteristics, i.e. minimum distance and duration, are finally fully characterised in terms of $mbox{Pe}$.