Emulsions are paramount in various interdisciplinary topical areas, yet a satisfactory understanding of their behavior in buoyancy-driven thermal flows has not been established. In the present work, we unravel the dynamical regimes of thermal convection in emulsions by leveraging a large set of mesoscale numerical simulations. Emulsions are prepared with a given volume fraction of the initially dispersed phase, φ, ranging from dilute (low values of φ) to jammed emulsions (high values of φ), resulting in different rheological responses of the emulsion, i.e., from Newtonian to non-Newtonian yield-stress behaviors, respectively. We then characterize the dynamics of the emulsions in the paradigmatic setup of the Rayleigh-Bénard convection, i.e., when confined between two parallel walls at different temperatures under the effect of buoyancy forces, the latter encoded in the dimensionless Rayleigh number Ra. We thoroughly investigated the dynamics of the emulsion in the changing of φ and Ra. For a given φ, at increasing Ra, we observe that the emulsion exhibits convection states, where structural changes may appear (i.e., droplet breakup, coalescence, or phase inversion), which inevitably impact the emulsion rheology. For sufficiently high values of Ra, two states of convection are observed: for low/moderate values of φ (Newtonian emulsions), we observe breakup-dominated dynamics, whereas for high values of φ (non-Newtonian emulsions), we observe phase-inverted states. For both scenarios, the droplet size distribution depends on Ra, and scaling laws for the average droplet size are analyzed and quantified. Our results offer insights into the rich dynamics of emulsions under thermal convection, offering a detailed characterization of the various dynamical regimes to be expected and their relation with structural changes occurring in such complex fluids.

Based on mesoscale lattice Boltzmann numerical simulations, we characterize the Rayleigh-Bénard (RB) convective dynamics of dispersions of liquid droplets in another liquid phase. Our numerical methodology allows us to modify the droplets’ interfacial properties to mimic the presence of an emulsifier (e.g., a surfactant), resulting in a positive disjoining pressure which stabilizes the droplets against coalescence. To appreciate the effects of this interfacial stabilization on the RB convective dynamics, we carry out a comparative study between a proper emulsion, i.e., a system where the stabilization mech- anism is present (stabilized liquid-liquid dispersion), and a system where the stabilization mechanism is absent (nonstabilized liquid-liquid dispersion). The study is conducted by systematically changing both the volume fraction φ and the Rayleigh number Ra. We find that the morphology of the two systems is dramatically different due to the different inter- facial properties. However, the two systems exhibit similar global heat transfer properties, expressed via the Nusselt number Nu. Significant differences in heat transfer emerge at smaller scales, which we analyze via the Nusselt number defined at mesoscales Numes. In particular, stabilized systems exhibit more intense mesoscale heat flux fluctuations due to the persistence of fluid velocity fluctuations down to small scales, which are instead dissipated in the interfacial dynamics of nonstabilized dispersions. For fixed Ra, the difference in mesoscale heat-flux fluctuations depends nontrivially on φ, featuring a maximum in the range 0.1 < φ < 0.2. Taken all together, our results highlight the role of interfacial physics in mesoscale convective heat transfer of complex fluids.

The dynamics of stabilised concentrated emulsions presents a rich phenomenology including chaotic emulsification, non-Newtonian rheology and ageing dynamics at rest. Macroscopic rheology results from the complex droplet microdynamics and, in turn, droplet dynamics is influenced by macroscopic flows via the competing action of hydrodynamic and interfacial stresses, giving rise to a complex tangle of elastoplastic effects, diffusion, breakups and coalescence events. This tight multiscale coupling, together with the daunting challenge of experimentally investigating droplets under flow, has hindered the understanding of concentrated emulsions dynamics. We present results from three-dimensional numerical simulations of emulsions that resolve the shape and dynamics of individual droplets, along with the macroscopic flows. We investigate droplet dispersion statistics, measuring probability density functions (p.d.f.s) of droplet displacements and velocities, changing the concentration, in the stirred and ageing regimes. We provide the first measurements, in concentrated emulsions, of the relative droplet–droplet separations p.d.f. and of the droplet acceleration p.d.f., which becomes strongly non-Gaussian as the volume fraction is increased above the jamming point. Cooperative effects, arising when droplets are in contact, are argued to be responsible of the anomalous superdiffusive behaviour of the mean square displacement and of the pair separation at long times, in both the stirred and in the ageing regimes. This superdiffusive behaviour is reflected in a non-Gaussian pair separation p.d.f., whose analytical form is investigated, in the ageing regime, by means of theoretical arguments. This work paves the way to developing a connection between Lagrangian dynamics and rheology in concentrated emulsions.

We study the process of thermal convection in jammed emulsions with a yield-stress rheology. We find that heat transfer occurs via an intermittent mechanism, whereby intense short-lived convective “heat bursts” are spaced out by long-lasting conductive periods. This behavior is the result of a sequence of fluidization-rigidity transitions, rooted in a nontrivial interplay between emulsion yield-stress rheology and plastic activity, which we characterize via a statistical analysis of the dynamics at the droplet scale. We also show that droplets’ coalescence induced during heat bursts leads to a spatially heterogeneous phase inversion of the emulsion which eventually supports a sustained convective state.

We present a mathematical model describing the evolution of sea ice and meltwater during summer. The system is described by two coupled partial differential equations for the ice thickness h(x,t) and pond depth w(x,t) fields. The model is similar, in principle, to the one put forward by Luthije et al. (2006), but it features i) a modified melting term, ii) a non-uniform seepage rate of meltwater through the porous ice medium and a minimal coupling with the atmosphere via a surface wind shear term, ?s (Scagliarini et al. 2020). We test, in particular, the sensitivity of the model to variations of parameters controlling fluid- dynamic processes at the pond level, namely the variation of turbulent heat flux with pond depth and the lateral melting of ice enclosing a pond. We observe that different heat flux scalings determine different rates of total surface ablations, while the system is relatively robust in terms of probability distributions of pond surface areas. Finally, we study pond morphology in terms of fractal dimensions, showing that the role of lateral melting is minor, whereas there is evidence of an impact from the initial sea ice topography.

Switchable and adaptive substrates emerged as valuable tools for controlling wetting and actuation of droplet motion. Here, we report a computational study of the dynamics of an unstable thin liquid film deposited on a switchable substrate, modeled with a space- and time-varying contact angle. For a sufficiently large rate of wettability variation, a topological transition appears. Instead of breaking up into droplets, as expected for a substrate with multiple wetting minima, a metastable rivulet state emerges. A criterion discriminating whether or not rivulets occur is identified in terms of a single dimensionless parameter. Finally, we show and derive theoretically how the film rupture times, droplet shape, and rivulet lifetime depend on the pattern wavelength and speed.

Active turbulence is a paradigmatic and fascinating example of self-organized motion at large scales occurring in active matter. We employ massive hydrodynamic simulations of suspensions of resolved model microswimmers to tackle the phenomenon in semi-diluted conditions at a mesoscopic level. We measure the kinetic energy spectrum and find that it decays as k-3 over a range of interme- diate wavenumbers. The velocity distributions are of L ́evy type, a distinct difference with inertial turbulence. Furthermore, we propose a reduced order dynamical deterministic model for active turbulence, inspired to shell models for classical turbulence, whose numerical and analytical study confirms the spectrum power-law observed in the simulations and reveals hints of a non-Gaussian, intermittent, physics of active turbulence. Direct numerical simulations and modelling also agree in pointing to a phenomenological picture whereby, in the absence of an energy cascade `a la Richardson forbidden by the low Reynolds number regime, it is the coupling between fluid velocity gradients and bacterial orientation that gives rise to a multiscale dynamics.

Stabilised dense emulsions display a rich phenomenology connecting microstructure and rheology. In this work, we study how an emulsion with a finite yield stress can be built via large-scale stirring. By gradually increasing the volume fraction of the dispersed minority phase, under the constant action of a stirring force, we are able to achieve a volume fraction close to 80%. Despite the fact that our system is highly concentrated and not yet turbulent we observe a droplet size distribution consistent with the -10/3 scaling, often associated with inertial range droplets breakup. We report that the polydispersity of droplet sizes correlates with the dynamics of the emulsion formation process. Additionally, we quantify the visco-elastic properties of the dense emulsion finally obtained and we demonstrate the presence of a finite yield stress. The approach reported can pave the way to a quantitative understanding of the complex interplay between the dynamics of mesoscale constituents and the large-scale flow properties of yield stress fluids.

We introduce a novel mesoscopic computational model based on a multiphase-multicomponent lattice Boltzmann method for the simulation of self-phoretic particles in the presence of liquid-liquid interfaces. Our model features fully resolved solvent hydrodynamics, and, thanks to its versatility, it can handle important aspects of the multiphysics of the problem, including particle wettability and differential solubility of the product in the two liquid phases. The method is extensively validated in simple numerical experiments, whose outcome is theoretically predictable, and then applied to the study of the behavior of active particles next to and trapped at interfaces. We show that their motion can be variously steered by tuning relevant control parameters, such as the phoretic mobilities, the contact angle, and the product solubility.

The sedimentation process in an active suspension is the result of the competition between gravity and the autonomous motion of particles. We carry out simulations of run-and-tumble squirmers that move in a fluid medium, focusing on the dependence of the non-equilibrium steady state on the swimming properties. We find that for large enough activity, the density profiles are no longer simple exponentials; we recover the numerical results through the introduction of a local effective temperature, suggesting that the breakdown of the Perrin-like exponential form is a collective effect due to fluid-mediated dynamic correlations among particles. We show that analogous concepts can also fit the case of active non-motile particles, for which we report the first study of this kind. Moreover, we provide evidence of scenarios where the solvent hydrodynamics induces non-local effects which require the full three-dimensional dynamics to be taken into account in order to understand sedimentation in active suspensions. Finally, analyzing the statistics of the orientations of microswimmers, the emergence of a height-dependent polar order in the system is discussed.

Switchable and adaptive substrates emerged as valuable tools for the control of wetting and actuation of droplet motion. Here we report a computational study of the dynamics of an unstable thin liquid film deposited on a switchable substrate, modelled with a space and time varying contact angle. With a static pattern, all the fluid is drained into droplets located around contact angle minima, whereas for a sufficiently large rate of wettability variation a state consisting of metastable rivulets is observed. A criterion discriminating whether rivulets can be observed or not is identified in terms of a single dimensionless parameter. Finally, we show and explain theoretically how the film rupture times, droplet shape and rivulet life time depend on the pattern wavelength and speed.

We study numerically the effect of thermal fluctuations and of variable fluid-substrate interactions on the spontaneous dewetting of thin liquid films. To this aim, we use a recently developed lattice Boltzmann method for thin liquid film flows, equipped with a properly devised stochastic term. While it is known that thermal fluctuations yield shorter rupture times, we show that this is a general feature of hydrophilic substrates, irrespective of the contact angle $\theta$. The ratio between deterministic and stochastic rupture times, though, decreases with $\theta$. Finally, we discuss the case of fluctuating thin film dewetting on chemically patterned substrates and its dependence on the form of the wettability gradients.

We present mesoscale numerical simulations of Rayleigh-Bénard (RB) convection in a two-dimensional model emulsion. The systems under study are constituted of finite-size droplets, whose concentration is systematically varied from small (Newtonian emulsions) to large values (non-Newtonian emulsions). We focus on the characterisation of the heat transfer properties close to the transition from conductive to convective states, where it is well known that a homogeneous Newtonian system exhibits a steady flow and a time-independent heat flux. In marked contrast, emulsions exhibit non-steady dynamics with fluctuations in the heat flux. In this paper, we aim at the characterisation of such non-steady dynamics via detailed studies on the time-averaged heat flux and its fluctuations. To quantitatively understand the time-averaged heat flux, we propose a side-by-side comparison between the emulsion system and a single-phase (SP) system, whose viscosity is suitably constructed from the shear rheology of the emulsion. We show that such local closure works well only when a suitable degree of coarse-graining (at the droplet scale) is introduced in the local viscosity. To delve deeper into the fluctuations in the heat flux, we furthermore propose a side-by-side comparison between a Newtonian emulsion (i.e., with a small droplet concentration) and a non-Newtonian emulsion (i.e., with a large droplet concentration), at fixed time-averaged heat flux. This comparison elucidates that finite-size droplets and the non-Newtonian rheology cooperate to trigger enhanced heat-flux fluctuations at the droplet scales. These enhanced fluctuations are rooted in the emergence of space correlations among distant droplets, which we highlight via direct measurements of the droplets displacement and the characterisation of the associated correlation function. The observed findings offer insights on heat transfer properties for confined systems possessing finite-size constituents.

We investigate the rheology of strain-hardening spherical capsules, from the dilute tothe concentrated regime under a confined shear flow using three-dimensional numericalsimulations. We consider the effect of capillary number, volume fraction and membrane inextensibility on the particle deformation and on the effective suspension viscosity andnormal stress differences of the suspension. The suspension displays a shear-thinningbehaviour that is a characteristic of soft particles such as emulsion droplets, vesicles,strain-softening capsules and red blood cells. We find that the membrane inextensibilityplays a significant role on the rheology and can almost suppress the shear-thinning. Forconcentrated suspensions a non-monotonic dependence of the normal stress differenceson the membrane inextensibility is observed, reflecting a similar behaviour in the particleshape. The effective suspension viscosity, instead, grows and eventually saturates, for verylarge inextensibilities, approaching the solid particle limit. In essence, our results revealthat strain-hardening capsules share rheological features with both soft and solid particlesdepending on the ratio of the area dilatation to shear elastic modulus. Furthermore, thesuspension viscosity exhibits a universal behaviour for the parameter space defined by the capillary number and the membrane inextensibility, when introducing the particlegeometrical changes at the steady state in the definition of the volume fraction.

The rheology of pressure-driven flows of two-dimensional dense monodisperse emulsions in neutral wetting microchannels is investigated by means of mesoscopic lattice Boltzmann simulations, capable of handling large collections of droplets, in the order of several hundreds. The simulations reveal that the fluidization of the emulsion proceeds through a sequence of discrete steps, characterized by yielding events whereby layers of droplets start rolling over each other, thus leading to sudden drops of the relative effective viscosity. It is shown that such discrete fluidization is robust against loss of confinement, namely it persists also in the regime of small ratios of the droplet diameter over the microchannel width. We also develop a simple phenomenological model which predicts a linear relation between the relative effective viscosity of the emulsion and the product of the confinement parameter (global size of the device over droplet radius) and the viscosity ratio between the disperse and continuous phases. The model shows excellent agreement with the numerical simulations. The present work offers new insights to enable the design of microfluidic scaffolds for tissue engineering applications and paves the way to detailed rheological studies of soft-glassy materials in complex geometries.

We present a mathematical model describing the evolution ofsea ice andmeltwater during summer. The system is described by two coupled partial differential equations for the ice thickness h and pond depth w fields. We test the sensitivity of the model to variations of parameters controlling fluid-dynamic processes at the pond level, namely the variation of turbulent heat flux with pond depth and the lateral melting of ice enclosing a pond. We observe that different heat flux scalings determine different rates of total surface ablations, while the system is relatively robust in terms of probability distributions of pond surface areas. Finally, we study pond morphology in terms of fractal dimensions, showing that the role of lateral melting is minor, whereas there is evidence of an impact from the initial sea ice topography.

Active fluids comprise a variety of systems composed of elements immersed in a fluid environment which can convert some form of energy into directed motion; as such they are intrinsically out-of-equilibrium in the absence of any external force. A fundamental problem in the physics of active matter concerns the understanding of how the characteristics of autonomous propulsion and agent-agent interactions determine the collective dynamics of the system. We study numerically the suspensions of self-propelled diffusiophoretic colloids, in (quasi)-2d configurations, accounting for both dynamically resolved solute-mediated phoretic interactions and solvent-mediated hydrodynamic interactions. Our results show that the system displays different scenarios at changing the colloid-solute affinity and it develops a cluster phase in the chemoattractive case. We study the statistics of cluster sizes and cluster morphologies for different magnitudes of colloidal activity. Finally, we provide evidences that hydrodynamics plays a relevant role in the aggregation kinetics and cluster morphology, significantly hindering cluster growth.

The dynamics of active colloids is very sensitive to the presence of boundaries and interfaces which therefore can be used to control their motion. Here we analyze the dynamics of active colloids adsorbed at a fluid-fluid interface. By using a mesoscopic numerical approach which relies on an approximated numerical solution of the Navier-Stokes equation, we show that when adsorbed at a fluid interface, an active colloid experiences a net torque even in the absence of a viscosity contrast between the two adjacent fluids. In particular, we study the dependence of this torque on the contact angle of the colloid with the fluid-fluid interface and on its surface properties. We rationalize our results via an approximate approach which accounts for the appearance of a local friction coefficient. By providing insight into the dynamics of active colloids adsorbed at fluid interfaces, our results are relevant for two-dimensional self assembly and emulsion stabilization by means of active colloids.

We explore the impact of geometrical corrugations on the near-wall flow properties of a soft material driven in a confined rough microchannel. By means of numerical simulations, we perform a quantitative analysis of the relation between the flow rate ? and the wall stress ?w for a number of setups, by changing both the roughness values as well as the roughness shape. Roughness suppresses the flow, with the existence of a characteristic value of ?w at which flow sets in. Just above the onset of flow, we quantitatively analyze the relation between ? and ?w. While for smooth walls a linear dependence is observed, steeper behaviours are found to set in by increasing wall roughness. The variation of the steepness, in turn, depends on the shape of the wall roughness, wherein gentle steepness changes are promoted by a variable space localization of the roughness.