BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Detlef Lohse

Detlef Lohse contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dyncond-mat.softphysics.app-phcond-mat.mtrl-sciphysics.chem-phphysics.comp-phnlin.CDphysics.ao-phphysics.geo-phphysics.ins-detphysics.med-phphysics.opticsphysics.soc-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
79works
0followers
13topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

79 published item(s)

preprint2026arXiv

Settling of chiral particles in a turbulent flow

Chiral particles are experimentally investigated while settling inwater with various turbulence intensity levels. The locations and orientations of the particles are tracked over time, allowing the close investigation of the particles' settling dynamics. The generated turbulent flow is measured using laser Doppler anemometry (LDA), and the turbulence strength varies between experiments in the range $0 \leq Re_λ\leq 250$. Starting with quiescent particle settling, the chiral particle's orientation dynamics are studied, revealing a preferred alignment and a strong translation-rotation coupling. The particle chirality determines the preferred rotation direction, though the alignment and translation-rotation coupling gradually vanish with increasing turbulence. We identify multiple settling modes for the chiral particles, which are characterised by the evolution of the rotation angles. Finally, a theoretical model assuming a simplified chiral particle in Stokes flow clarifies the emergence of each settling mode.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Transition from classical to ultimate melting

Melting is omnipresent in nature and technology, with applications ranging from metallurgy, biology, food science, and latent thermal energy storage to oceanography, geophysics, and climate science, and occurring on all scales from sub-millimeter to global scales. The key objective is to understand the rate at which an object melts as a function of its size and of the ambient conditions. To achieve this it is important to be able to extrapolate from small scale experiments and observations to large or even global scales. This is done by scaling laws. However, these are only meaningful if there is no transition from one scaling relation to another one. Here we show, however, that for both fixed and freely-advected melting objects immersed in a turbulent flow a melting transition does exist, namely from slow melting at the small scales to fast melting at the large scales. We do so by controlled melting experiments and corresponding direct numerical simulations, covering four orders of magnitude in scale. The transition corresponds to the transition from a laminar-type boundary layer around the melting object to a turbulent-type boundary layer, i.e., from so-called classical turbulence to ultimate turbulence, with its enhanced transport properties. Our results thus provide a quantitative understanding of the flow physics of the melting process and thereby enable a better extrapolation and prediction of melt rates on large scales such as relevant in geophysics, oceanography, and climate science.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Collective effects of neighbouring melting ice objects

We present a study on the melting dynamics of neighbouring ice bodies by means of idealised simulations, focusing on collective effects, with the goal of obtaining fundamental insight into how collective interactions influence the melting of ice. Two neighbouring (vertically or horizontally aligned), square-shaped, and equally sized ice objects (size on the order of centimetres) are immersed in quiescent fresh water at a temperature of 20°C. By performing two-dimensional direct numerical simulations, and using the phase-field method to model the phase change, the collective melting of these objects is studied. When the objects are horizontally aligned, no significant influence of the neighbouring object on the melting time is observed. On the other hand, when vertically aligned, though the melting of the upper object is mostly unaffected, the melting time and the morphology of the lower ice body strongly depends on the initial inter-object distance. We report that the melting of the bottom object can be enhanced by more than 10%, or delayed more than 20%, displaying a non-monotonic dependence on the initial object size. We show that this behaviour results from a non-trivial competition between layering of cold fluid, which lowers the heat transfer, and convective flows, which favour mixing and heat transfer. For this melting in mixed convection, we were able to collapse our data onto a single curve.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Lifetimes of metastable windy states in two-dimensional Rayleigh-Bénard convection with stress-free boundaries

Two-dimensional horizontally periodic Rayleigh-Bénard convection between stress-free boundaries displays two distinct types of states, depending on the initial conditions. Roll states are composed of pairs of counter-rotating convection rolls. Windy states are dominated by strong horizontal wind (also called zonal flow) that is vertically sheared, precludes convection rolls, and suppresses heat transport. Windy states occur only when the Rayleigh number $Ra$ is sufficiently above the onset of convection. At intermediate $Ra$ values, windy states can be induced by suitable initial conditions, but they undergo a transition to roll states after finite lifetimes. At larger $Ra$ values, where windy states have been observed for the full duration of simulations, it is unknown whether they represent chaotic attractors or only metastable states that would eventually undergo a transition to roll states. We study this question using direct numerical simulations of a fluid with a Prandtl number of 10 in a layer whose horizontal period is 8 times its height. At each of seven $Ra$ values between $9\times10^6$ and $2.25\times10^7$ we have carried out 200 or more simulations, all from initial conditions leading to windy convection with finite lifetimes. The lifetime statistics at each $Ra$ indicate a memoryless process with survival probability decreasing exponentially in time. The mean lifetimes grow with $Ra$ approximately as $Ra^4$. This analysis provides no $Ra$ value at which windy convection becomes stable; it might remain metastable at larger $Ra$ with extremely long lifetimes.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Rising and settling 2D cylinders with centre-of-mass offset

Rotational effects are commonly neglected when considering the dynamics of freely rising or settling isotropic particles. Here, we demonstrate that particle rotations play an important role for rising as well as for settling cylinders in situations when mass eccentricity, and thereby a new pendulum timescale, is introduced to the system. We employ two-dimensional simulations to study the motion of a single cylinder in a quiescent unbounded incompressible Newtonian fluid. This allows us to vary the Galileo number, density ratio, relative moment of inertia, and Centre-Of-Mass offset (COM) systematically and beyond what is feasible experimentally. For certain buoyant density ratios, the particle dynamics exhibit a resonance mode, during which the coupling via the Magnus lift force causes a positive feedback between translational and rotational motions. This mode results in vastly different trajectories with significantly larger rotational and translational amplitudes and an increase of the drag coefficient easily exceeding a factor two. We propose a simple model that captures how the occurrence of the COM offset induced resonance regime varies, depending on the other input parameters, specifically the density ratio, the Galileo number, and the relative moment of inertia. Remarkably, depending on the input parameters, resonance can be observed for centre-of-mass offsets as small as a few percent of the particle diameter, showing that the particle dynamics can be highly sensitive to this parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Double-diffusive transport in multicomponent vertical convection

Motivated by the ablation of vertical ice faces in salt water, we use three-dimensional direct numerical simulations to investigate the heat and salt fluxes in two-scalar vertical convection. For parameters relevant to ice-ocean interfaces in the convection-dominated regime, we observe that the salinity field drives the convection and that heat is essentially transported as a passive scalar. By varying the diffusivity ratio of heat and salt (i.e., the Lewis number $Le$), we identify how the different molecular diffusivities affect the scalar fluxes through the system. Away from the walls, we find that the heat transport is determined by a turbulent Prandtl number of $Pr_t\approx 1$ and that double-diffusive effects are practically negligible. However, the difference in molecular diffusivities plays an important role close to the boundaries. In the (unrealistic) case where salt diffused faster than heat, the ratio of salt-to-heat fluxes would scale as $Le^{1/3}$, consistent with classical nested scalar boundary layers. However, in the realistic case of faster heat diffusion (relative to salt), we observe a transition towards a $Le^{1/2}$ scaling of the ratio of the fluxes. This coincides with the thermal boundary layer width growing beyond the thickness of the viscous boundary layer. We find that this transition is not determined by a critical Lewis number, but rather by a critical Prandtl number $Pr\approx 10$, slightly below that for cold seawater where $Pr=14$. We compare our results to similar studies of sheared and double-diffusive flow under ice shelves, and discuss the implications for fluxes in large-scale ice-ocean models. By coupling our results to ice-ocean interface thermodynamics, we describe how the flux ratio impacts the interfacial salinity, and hence the strength of solutal convection and the ablation rate.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Elastocapillary Worthington jets

The retraction of an impacting droplet on a non-wetting substrate is often associated with the formation of a Worthington jet, which is fed by the retracting liquid. A non-Newtonian rheology of the liquid is known to affect the retraction of the impacting droplet. Here we present a novel phenomenon related to the impact of viscoelastic droplets on non-wettable substrates. We reveal that the viscoelasticity of the liquid results in an \emph{elastocapillary} regime in the stretching Worthington jet, distinguished by a pinned contact line and a slender jet that does not detach from the droplet. We identify the impact conditions, in the Weber number -- Deborah number phase space, for observing these \emph{elastocapillary} Worthington jets. Such jets exhibit an effectively nearly linear (in time) variation of the strain rate. Upon further extension, the jet exhibits beads-on-a-string structures, characteristic of the \emph{elastocapillary} thinning of slender viscoelastic liquid filaments. The \emph{elastocapillary} Worthington jet is not only relevant for a droplet impact on a solid substrate scenario, but can also be expected in other configurations where a Worthington jet is observed for viscoelastic liquids, such as drop impact on a liquid pool and bubble bursting at an interface.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Data-driven identification of the spatio-temporal structure of turbulent flows by streaming Dynamic Mode Decomposition

Streaming Dynamic Mode Decomposition (sDMD) (Hemati et al., Phys. Fluids 26(2014)) is a low-storage version of Dynamic Mode Decomposition (DMD) (Schmid, J. Fluid Mech. 656 (2010)), a data-driven method to extract spatio-temporal flow patterns. Streaming DMD avoids storing the entire data sequence in memory by approximating the dynamic modes through incremental updates with new available data. In this paper, we use sDMD to identify and extract dominant spatio-temporal structures of different turbulent flows, requiring the analysis of large datasets. First, the efficiency and accuracy of sDMD are compared to the classical DMD, using a publicly available test dataset that consists of velocity field snapshots obtained by direct numerical simulation of a wake flow behind a cylinder. Streaming DMD not only reliably reproduces the most important dynamical features of the flow; our calculations also highlight its advantage in terms of the required computational resources. We subsequently use sDMD to analyse three different turbulent flows that all show some degree of large-scale coherence: rapidly rotating Rayleigh--Bénard convection, horizontal convection and the asymptotic suction boundary layer. Structures of different frequencies and spatial extent can be clearly separated, and the prominent features of the dynamics are captured with just a few dynamic modes. In summary, we demonstrate that sDMD is a powerful tool for the identification of spatio-temporal structures in a wide range of turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Droplet dissolution driven by emerging thermal gradients and Marangoni flow

The lifetime $τ$ of an isothermal and purely diffusively dissolving droplet in a host liquid scales as $τ\sim R_0^2$ with its initial radius $R_0$ [Langmuir, Phys. Rev. 12, 368 (1919)]. For a droplet dissolving due to natural convection driven by density differences, its lifetime scales as $τ\sim R_0^{5/4}$ [Dietrich et al., J. Fluid Mech. 794, 45 (2016)]. In this paper we experimentally find and theoretically derive yet another droplet dissolution behavior, resulting in $τ\sim R_0^4$. It occurs when the dissolution dynamics is controlled by local heating of the liquid, leading to a modified solubility and a thermal Marangoni flow around the droplet. The thermal gradient is achieved by plasmonic heating of a gold nanoparticle decorated sample surface, on which a sessile water droplet immersed in water-saturated 1-butanol solution is sitting. The resulting off-wall thermal Marangoni flow and the temperature dependence of the solubility determine the droplet dissolution rate, resulting in a shrinkage $R(t) \sim (τ-t )^{1/4}$ of the droplet radius and thus in $τ\sim R_0^{4}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Enhancing Heat Transport in Multiphase Rayleigh-Bénard Turbulence by Changing the Plate-Liquid Contact Angles

This numerical study presents a simple but extremely effective way to considerably enhance heat transport in turbulent multiphase flows, namely by using oleophilic walls. As a model system, we pick the Rayleigh-Bénard setup, filled with an oil-water mixture. For oleophilic walls, e.g. 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 then only about $20\%$ as compared to pure water. The physical explanation of the highly-efficient heat transport for oleophilic walls is that thermal plumes detach from the oil-rich boundary layer and are transported together with the oil phase. In the bulk, the oil-water interface prevents the plumes to mix with the turbulent water bulk. To confirm this physical picture, we show that the minimum amount of oil 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 thermal turbulence. Moreover, the physical insight of how coherent structures are coupled with one phase of a two-phase system has very general applicability for controlling transport properties in other turbulent multiphase flows.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Entrapment of Interfacial Nanobubbles on Nano Structured Surfaces

Spherical-cap-shaped interfacial nanobubbles (NBs) forming on hydrophobic surfaces in aqueous solutions have extensively been studied both from a fundamental point of view and due to their relevance for various practical applications. In this study, the nucleation mechanism of spontaneously generated NBs at solid-liquid interfaces of immersed nanostructured hydrophobic surfaces is studied. Depending on the size and density of the surface nanostructures, NBs with different size and density were reproducibly and deterministically obtained. A two-step process can explain the NB nucleation, based on the crevice model, i.e., entrapped air pockets in surface cavities which grow by diffusion. The results show direct evidence for the spontaneous formation of NBs on a surface at its immersion. Next, the influence of size and shape of the nanostructures on the nucleated NBs are revealed. In particular, on non-circular nanopits we obtain NBs with a non-circular footprint, demonstrating the strong pinning forces at the three-phase contact line.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Impact forces of water drops falling on superhydrophobic surfaces

A falling liquid drop, after impact on a rigid substrate, deforms and spreads, owing to the normal reaction force. Subsequently, if the substrate is non-wetting, the drop retracts and then jumps off. As we show here, not only is the impact itself associated with a distinct peak in the temporal evolution of the normal force, but also the jump-off, which was hitherto unknown. We characterize both peaks and elucidate how they relate to the different stages of the drop impact process. The time at which the second peak appears coincides with the formation of a Worthington jet, emerging through flow-focusing, and it is independent of the impact velocity. However, the magnitude of this peak is dictated by the drop's inertia and surface tension. We show that even low-velocity impacts can lead to a surprisingly high peak in the normal force, namely when a more pronounced singular Worthington jet occurs due to the collapse of an air cavity in the drop.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Micro-droplet nucleation through solvent exchange in a turbulent buoyant jet

Solvent exchange is a process involving mixing between a good solvent with dissolved solute and a poor solvent. The process creates local oversaturation which causes the nucleation of minute solute droplets. Such ternary systems on a macro-scale have remained unexplored in the turbulent regime. We experimentally study the solvent exchange process by injecting mixtures of ethanol and trans-anethole into water, forming a turbulent buoyant jet in upward direction. Locally, turbulent mixing causes oversaturation of the trans-anethole following turbulent entrainment. We optically measure the concentration of the nucleated droplets using a light attenuation technique and find that the radial concentration profile is sub-Gaussian. In contrast to the entrainment-based models, the spatial evolution of the oversaturation reveals continuous droplet nucleation downstream and radially across the jet, which we attribute to the limited mixing capacity of the jet. Though we are far from a full quantitative understanding, this work extends the knowledge on solvent exchange into the turbulent regime, and brings in a novel type of flow, broadening the scope of multicomponent, multiphase turbulent jets with phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Oscillatory droplet dissolution from competing Marangoni and gravitational flows

The dissolution or growth of a droplet in a host liquid is an important part for processes like chemical extraction, chromatography or emulsification. In this work we look at the dissolution of a pair of vertically aligned droplets immersed in water, both experimentally and with numerical simulations. The liquids used for the droplets are long chain alcohols with a low but finite solubility in water and a significantly lower density than that of the host liquid. Therefore, a solutal plume is formed above of the bottom droplet and natural convection dominates the dissolution process. We monitor the volume of the droplets and the velocity field around them over time. When the liquids of the two droplets are the same, our previously found scaling laws for the Sherwood and Reynolds numbers as functions of the Rayleigh number (Dietrich et al., 2016, J. Fluid Mech.) can be applied to the lower droplet. However, remarkably, when the liquid of the top droplet is different than that of the bottom droplet the volume as function of time becomes non-monotonic, and an oscillatory Marangoni flow at the top droplet is observed. We identify the competition between solutal Marangoni flow and density driven convection as the origin of the oscillation, and numerically model the process.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Physiochemical hydrodynamics of the phase segregation in an evaporating binary microdroplet

Phase segregation triggered by selective evaporation can emerge in multicomponent systems, leading to complex physiochemical hydrodynamics. Recently, Li et al. (Phys. Rev. Lett., vol. 120, 2018, 224501) and Kim & Stone (J. Fluid Mech., vol. 850, 2018, pp. 769-783) reported a segregative behavior (i.e., demixing) in an evaporating binary droplet. In this work, by means of experiments and theoretical analysis, we investigate the flow dynamics after the occurrence of the phase segregation. As example, we take the 1,2-hexanediol-water binary droplet system. First, we experimentally reveal the overall physiochemical hydrodynamics of the evaporation process, including the segregative behavior and the resulting flow structure close to the substrate. By quantifying the evolution of the radial flow, we identify three successive life stages of the evaporation process. At Stage I, a radially outward flow is observed. It is driven by the Marangoni effect. At the transition to Stage II, the radial flow partially reverses, starting from the contact line. This flow breaks the axial symmetry and remarkably is driven by the segregation itself. Finally at Stage III, the flow decays as the evaporation gradually ceases. At this stage the segregation has grown to the entire droplet, and the flow is again controlled by the Marangoni effect. The resulting Marangoni flow homogenizes the distribution of the entrapped volatile water over the whole droplet.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Surface Properties of Colloidal Particles Affect Colloidal Self-Assembly in Evaporating Self-Lubricating Ternary Droplets

In this work, we unravel the role of surface properties of colloidal particles on the formation of supraparticles (clusters of colloidal particles) in a colloidal Ouzo droplet. Self-lubricating colloidal Ouzo droplets are an efficient and simple approach to form supraparticles, overcoming the challenge of the coffee stain effect in situ. Supraparticles are an efficient route to high-performance materials in various fields, from catalysis to carriers for therapeutics. Yet, the role of the surface of colloidal particles in the formation of supraparticles using Ouzo droplets remains unknown. Therefore, we used silica particles as a model system and compared sterically stabilized versus electrostatically stabilized silica particles - positively and negatively charged. Additionally, we studied the effect of hydration. Hydrated negatively charged silica particles and sterically stabilized silica particles form supraparticles. Conversely, dehydrated negatively charged silica particles and positively charged amine-coated particles form flat film-like deposits. Notably, the assembly process is different for all the four types of particles. The surface modifications alter (a) the contact line motion of the Ouzo droplet and (b) the particle-oil and particle-substrate interactions. These alterations modify the particle accumulation at the various interfaces, which ultimately determines the shape of the final deposit. Thus, by modulating the surface properties of the colloidal particles, we can tune the shape of the final deposit, from a spheroidal supraparticle to a flat deposit. In the future, this approach can be used to tailor the supraparticles for applications such as optics and catalysis, where the shape affects the functionality.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Taylor-Culick retractions and the influence of the surroundings

When a freely suspended liquid film ruptures, it retracts spontaneously under the action of surface tension. If the film is surrounded by air, the retraction velocity is known to approach the constant Taylor-Culick velocity. However, when surrounded by an external viscous medium, the dissipation within that medium dictates the magnitude of the retraction velocity. In the present work, we study the retraction of a liquid (water) film in a viscous oil ambient (\emph{two-phase} Taylor-Culick retractions), and that sandwiched between air and a viscous oil (\emph{three-phase} Taylor-Culick retractions). In the latter case, the experimentally-measured retraction velocity is observed to have a weaker dependence on the viscosity of the oil phase as compared to the configuration where the water film is surrounded completely by oil. Numerical simulations indicate that this weaker dependence arises from the localization of viscous dissipation near the three-phase contact line. The speed of retraction only depends on the viscosity of the surrounding medium and not on that of the film. From the experiments and the numerical simulations, we reveal unprecedented regimes for the scaling of the Weber number $We_f$ of the film (based on its retraction velocity) or the capillary number $Ca_s$ of the surroundings vs. the Ohnesorge number $Oh_s$ of the surroundings in the regime of large viscosity of the surroundings ($Oh_s \gg 1$), namely $We_f \sim Oh_s^{-2}$ and $Ca_s \sim Oh_s^{0}$ for the two-phase Taylor-Culick configuration, and $We_f \sim Oh_s^{-1}$ and $Ca_s \sim Oh_s^{1/2}$ for the three-phase Taylor-Culick configuration.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The asymmetric coalescence of two droplets with different surface tensions is caused by capillary waves

When two droplets with different surface tensions collide, the shape evolution of the merging droplets is asymmetric. Using experimental and numerical techniques, we reveal that this asymmetry is caused by asymmetric capillary waves, which are the result of the different surface tensions of the droplets. We show that the asymmetry is enhanced by increasing the surface tension difference, and suppressed by increasing the inertia of the colliding droplets. Furthermore, we study capillary waves in the limit of no inertia. We reveal that the asymmetry is not directly caused by Marangoni forces. In fact, somehow counterintuitive, asymmetry is strongly reduced by the Marangoni effect. Rather, the different intrinsic capillary wave amplitudes and velocities associated with the different surface tensions of the droplets lie at the origin of the asymmetry during droplet coalescence.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Turbulent Rayleigh-Bénard convection with bubbles attached to the plate

We numerically investigate turbulent Rayleigh-Bénard convection with gas bubbles attached to the hot plate, mimicking a core feature in electrolysis, catalysis, or boiling. The existence of bubbles on the plate reduces the global heat transfer due to the much lower thermal conductivity of gases as compared to liquids and changes the structure of the boundary layers. The numerical simulations are performed in 3D at Prandtl number Pr=4.38 (water) and Rayleigh number $10^7\le Ra \le10^8$. For simplicity, we assume the bubbles to be equally-sized and having pinned contact lines. We vary the total gas-covered area fraction $0.18 \le S_0 \le 0.62$, the relative bubble height $0.02\le h/H \le0.05$ (where $H$ is the height of the Rayleigh-Bénard cell), the bubble number $40 \le n \le 144$, and their spatial distribution. In all cases, asymmetric temperature profiles are observed, which we quantitatively explain based on the heat flux conservation at each horizontal section. We further propose the idea of using an equivalent single-phase setup to mimic the system with attached bubbles. Based on this equivalence, we can calculate the heat transfer. Without introducing any free parameter, the predictions for the Nusselt number, the upper and lower thermal boundary layer thicknesses, and the mean centre temperature well agree with the numerical results. Finally, our predictions also work for the cases with much larger Pr (e.g. $400$), which indicates that our results can also be applied to predict the mass transfer in water electrolysis with bubbles attached to the electrode surface or in catalysis.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Vorticity-induced flow-focusing leads to bubble entrainment in an inkjet printhead: synchrotron X-ray and volume-of-fluid visualizations

The oscillatory flows present in an inkjet printhead can lead to strong deformations of the air-liquid interface at the nozzle exit. Such deformations may lead to an inward directed air jet with bubble pinch-off and the subsequent entrainment of an air bubble, which is highly detrimental to the stability of inkjet printing. Understanding the mechanisms of bubble entrainment is therefore crucial in improving print stability. In the present work, we use ultrafast X-ray phase-contrast imaging and direct numerical simulations based on the Volume-of-Fluid method to study the mechanisms underlying the bubble entrainment in a piezo-acoustic printhead. We first demonstrate good agreement between experiments and numerics. We then show the different classes of bubble pinch-off obtained in experiments, and that those were also captured numerically. The numerical results are then used to show that the baroclinic torque, which is generated at the gas-liquid interface due to the misalignment of density and pressure gradients, results in a flow-focusing effect that drives the formation of the air jet from which a bubble can pinch-off.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An efficient phase-field method for turbulent multiphase flows

With the aim of efficiently simulating three-dimensional multiphase turbulent flows with a phase-field method, we propose a new discretization scheme for the biharmonic term (the 4th-order derivative term) of the Cahn-Hilliard equation. This novel scheme can significantly reduce the computational cost while retaining the same accuracy as the original procedure. Our phase-field method is built on top of a direct numerical simulation solver, named AFiD (www.afid.eu) and open-sourced by our research group. It relies on a pencil distributed parallel strategy and a FFT-based Poisson solver. To deal with large density ratios between the two phases, a pressure split method [1] has been applied to the Poisson solver. To further reduce computational costs, we implement a multiple-resolution algorithm which decouples the discretizations for the Navier-Stokes equations and the scalar equation: while a stretched wall-resolving grid is used for the Navier-Stokes equations, for the Cahn-Hilliard equation we use a fine uniform mesh. The present method shows excellent computational performance for large-scale computation: on meshes up to 8 billion nodes and 3072 CPU cores, a multiphase flow needs only slightly less than 1.5 times the CPU time of the single-phase flow solver on the same grid. The present method is validated by comparing the results to previous studies for the cases of drop deformation in shear flow, including the convergence test with mesh refinement, and breakup of a rising buoyant bubble with density ratio up to 1000. Finally, we simulate the breakup of a big drop and the coalescence of O(10^3) drops in turbulent Rayleigh-Bénard convection at a Rayleigh number of $10^8$, observing good agreement with theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Emergence of bimodal motility in active droplets

To explore and react to their environment, living micro-swimmers have developed sophisticated strategies for locomotion - in particular, motility with multiple gaits. To understand the physical principles associated with such a behavioural variability,synthetic model systems capable of mimicking it are needed. Here, we demonstrate bimodal gait switching in autophoretic droplet swimmers. This minimal experimental system is isotropic at rest, a symmetry that can be spontaneously broken due to the nonlinear coupling between hydrodynamic and chemical fields, inducing a variety of flow patterns that lead to different propulsive modes. We report a dynamical transition from quasi-ballistic to bimodal chaotic motion, controlled by the viscosity of the swimming medium. By simultaneous visualisation of the chemical and hydrodynamic fields, supported quantitatively by an advection-diffusion model, we show that higher hydrodynamic modes become excitable with increasing viscosity, while the recurrent mode-switching is driven by the droplet's interaction with self-generated chemical gradients. We further demonstrate that this gradient interaction results in anomalous diffusive swimming akin to self-avoiding spatial exploration strategies observed in nature.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Extended lifetime of respiratory droplets in a turbulent vapour puff and its implications on airborne disease transmission

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 10um) 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 William F. 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 two meters advection range of the respiratory droplets within one second. 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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fabrication of free-standing Pt nanowires for use as thermal anemometry probes in turbulence measurements

We report a robust fabrication method for patterning free-standing Pt nanowires for the use as thermal anemometry probes for small-scale turbulence measurements. Using e-beam lithography, high aspect ratio Pt nanowires (~300 nm width, ~70 $μ$m length, ~100 nm thickness) were patterned on the surface of oxidized silicon (Si) wafers. Combining precise wet etching processes with dry etching processes, these Pt nanowires have been successfully released free-standing between two silicon dioxide (SiO2) beams supported on Si cantilevers. Moreover, the unique design of the bridge holding the device allowed to release the device gently without damaging the Pt nanowires. The total fabrication time was minimized by restricting the use of e-beam lithography to the patterning of the Pt nanowires while standard photolithography was employed for other parts of the devices. We demonstrate that the fabricated sensors are suitable for turbulence measurements when operated in a constant-current mode. A robust calibration between output voltage and fluid velocity was established over the velocity range from 0.5 m s-1 to 5 m s-1 in an SF6 atmosphere at a pressure of 2 bar and a temperature of 21°C. The sensing signal from the nanowires showed negligible drift over a period of several hours. Moreover, we confirmed that the nanowires are able to withstand high dynamic pressures by testing them in air at room temperature velocities up to 55 m/s.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Heat transfer in turbulent Rayleigh-Bénard convection within two immiscible fluid layers

We numerically investigate turbulent Rayleigh-Bénard convection within two immiscible fluid layers, aiming to understand how the layer thickness and fluid properties affect the heat transfer (characterized by the Nusselt number $Nu$) in two-layer systems. Both two- and three-dimensional simulations are performed at fixed global Rayleigh number $Ra=10^8$, Prandtl number $Pr=4.38$, and Weber number $We=5$. We vary the relative thickness of the upper layer between $0.01 \le α\le 0.99$ and the thermal conductivity coefficient ratio of the two liquids between $0.1 \le λ_k \le 10$. Two flow regimes are observed: In the first regime at $0.04\leα\le0.96$, convective flows appear in both layers and $Nu$ is not sensitive to $α$. In the second regime at $α\le0.02$ or $α\ge0.98$, convective flow only exists in the thicker layer, while the thinner one is dominated by pure conduction. In this regime, $Nu$ is sensitive to $α$. To predict $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 $Nu$ and for the temperature at the interface well agree with our numerical results and previous experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Layering and vertical transport in sheared double diffusive convection in the diffusive regime

A sequence of two and three-dimensional simulations is conducted for the double diffusive convection (DDC) flows in the diffusive regime subjected to an imposed shear. The flow is confined between two horizontal plates which are maintained at different constant temperature, salinity, and different velocity, thus setting up a shear across the flow. The lower plate is fixed at higher temperature and salinity, while the overall (unperturbed) density gradient is statically stable. For a wide range of control parameters, and for sufficiently strong perturbation of the conductive initial state, we find that staircase-like structures spontaneously develop, with relatively well-mixed layers separated by sharp interfaces of enhanced scalar gradient. Such staircases appear to be robust even in the presence of strong shear over very long times, although we typically observe early time coarsening of the number of observed layers. For the same set of control parameters, different asymptotic layered states, with markedly different vertical scalar fluxes, can arise for different initial perturbation structures. The imposed shear does significantly spatio-temporally modify the vertical transport of the various scalars. The flux ratio (i.e., the ratio between the density fluxes due to the total (convective and diffusive) salt flux and the total heat flux) is found, at steady state, to be essentially equal to the square root of the ratio of the salt diffusivity to the thermal diffusivity, consistently with the physical model originally proposed by Linden and Shirtcliffe (1978) and the variational arguments presented by Stern (1982) for unsheared double diffusive convection.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Lifting a sessile oil drop from a superamphiphobic surface with an impacting one

Colliding drops are widely encountered in everyday technologies and natural processes, from combustion engines and commodity sprays to raindrops and cloud formation. The outcome of a collision depends on many factors, including the impact velocity and the degree of head-on alignment, in addition to intrinsic properties like surface tension. Yet little is known on the binary impact dynamics of low surface tension oil drops on a low-wetting surface. We experimentally and numerically investigate the dynamics of an oil drop impacting an identical sessile drop sitting on a superamphiphobic surface. We observe five rebound scenarios, four of which do not involve coalescence. We describe two previously unexplored cases for sessile oil drop lift-off, resulting from a drop-on-drop impact event. The simulations quantitatively reproduce all rebound scenarios and enable quantification of the velocity profiles, the energy transfer, and the viscous dissipation. Our results illustrate how varying the relative offset and the impact velocity results in controllable rebound dynamics for low surface tension drop collisions on superamphiphobic surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Off-center gravity induces large-scale flow patterns in spherical Rayleigh-Bénard

Inspired by the hemispherical asymmetry observed in the Earth's inner core, we perform direct numerical simulations to study the effect of the gravity center offset in spherical Rayleigh-Bénard convection. We find that even a minimal shift of the gravity center has a pronounced influence on the flow structures. When the gravity center is shifted towards the South, the co-latitudinal buoyancy component creates an energetic jet on the Northern side of the inner sphere that is directed towards the outer sphere. As a result, a large-scale meridional circulation is formed. However, surprisingly, the global heat flux is not affected by the changes in the large-scale flow organization induced by the gravity center offset. Our results suggest that the hemispherical core asymmetry is key to model the flow phenomena in the Earth's outer core and mantle.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On explosive boiling of a multicomponent Leidenfrost drop

The gasification of multicomponent fuel drops is relevant in various energy-related technologies. An interesting phenomenon associated with this process is the self-induced explosion of the drop, producing a multitude of smaller secondary droplets, which promotes overall fuel atomization and, consequently, improves the combustion efficiency and reduces emissions of liquid-fueled engines. Here, we study a unique explosive gasification process of a tricomponent droplet consisting of water, ethanol, and oil ("ouzo"), by high-speed monitoring of the entire gasification event taking place in the well-controlled, levitated Leidenfrost state over a superheated plate. It is observed that the preferential evaporation of the most volatile component, ethanol, triggers nucleation of the oil microdroplets/nanodroplets in the remaining drop, which, consequently, becomes an opaque oil-in-water microemulsion. The tiny oil droplets subsequently coalesce into a large one, which, in turn, wraps around the remnant water. Because of the encapsulating oil layer, the droplet can no longer produce enough vapor for its levitation, and, thus, falls and contacts the superheated surface. The direct thermal contact leads to vapor bubble formation inside the drop and consequently drop explosion in the final stage.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Plasmonic Microbubble Dynamics in Binary Liquids

The growth of surface plasmonic microbubbles in binary water/ethanol solutions is experimentally studied. The microbubbles are generated by illuminating a gold nanoparticle array with a continuous wave laser. Plasmonic bubbles exhibit ethanol concentration-dependent behaviors. For low ethanol concentrations (f_e) of < 67.5%, bubbles do not exist at the solid-liquid interface. For high f_e values of >80%, the bubbles behave as in pure ethanol. Only in an intermediate window of 67.5% < f_e < 80% do we find sessile plasmonic bubbles with a highly nontrivial temporal evolution, in which as a function of time three phases can be discerned. (1) In the first phase, the microbubbles grow, while wiggling. (2) As soon as the wiggling stops, the microbubbles enter the second phase in which they suddenly shrink, followed by (3) a steady reentrant growth phase. Our experiments reveal that the sudden shrinkage of the microbubbles in the second regime is caused by a depinning event of the three phase contact line. We systematically vary the ethanol concentration, laser power, and laser spot size to unravel water recondensation as the underlying mechanism of the sudden bubble shrinkage in phase 2.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Strong alignment of prolate ellipsoids in Taylor-Couette flow

We report on the mobility and orientation of finite-size, neutrally buoyant prolate ellipsoids (of aspect ratio $Λ=4$) in Taylor-Couette flow, using interface resolved numerical simulations. The setup consists of a particle-laden flow in between a rotating inner and a stationary outer cylinder. We simulate two particle sizes $\ell/d=0.1$ and $\ell/d=0.2$, $\ell$ denoting the particle major axis and $d$ the gap-width between the cylinders. The volume fractions are $0.01\%$ and $0.07\%$, respectively. The particles, which are initially randomly positioned, ultimately display characteristic spatial distributions which can be categorised into four modes. Modes $(i)$ to $(iii)$ are observed in the Taylor vortex flow regime, while mode ($iv$) encompasses both the wavy vortex, and turbulent Taylor vortex flow regimes. Mode $(i)$ corresponds to stable orbits away from the vortex cores. Remarkably, in a narrow $\textit{Ta}$ range, particles get trapped in the Taylor vortex cores (mode ($ii$)). Mode $(iii)$ is the transition when both modes $(i)$ and $(ii)$ are observed. For mode $(iv)$, particles distribute throughout the domain due to flow instabilities. All four modes show characteristic orientational statistics. We find the particle clustering for mode ($ii$) to be size-dependent, with two main observations. Firstly, particle agglomeration at the core is much higher for $\ell/d=0.2$ compared to $\ell/d=0.1$. Secondly, the $\textit{Ta}$ range for which clustering is observed depends on the particle size. For this mode $(ii)$ we observe particles to align strongly with the local cylinder tangent. The most pronounced particle alignment is observed for $\ell/d=0.2$ around $\textit{Ta}=4.2\times10^5$. This observation is found to closely correspond to a minimum of axial vorticity at the Taylor vortex core ($\textit{Ta}=6\times10^5$) and we explain why.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Universality in microdroplet nucleation during solvent exchange in Hele-Shaw like channels

Micro and nanodroplets have many important applications such as in drug delivery, liquid-liquid extraction, nanomaterial synthesis and cosmetics. A commonly used method to generate a large number of micro or nanodroplets in one simple step is solvent exchange (also called nanoprecipitation), in which a good solvent of the droplet phase is displaced by a poor one, generating an oversaturation pulse that leads to droplet nucleation. Despite its crucial importance, the droplet growth resulting from the oversaturation pulse in this ternary system is still poorly understood. We experimentally and theoretically study this growth in Hele-Shaw like channels by measuring the total volume of the oil droplets that nucleates out of it. In order to prevent the oversaturated oil from exiting the channel, we decorated some of the channels with a porous region in the middle. Solvent exchange is performed with various solution compositions, flow rates and channel geometries, and the measured droplets volume is found to increase with the Péclet number $Pe$ with an approximate effective power law $V\propto Pe^{0.50}$. A theoretical model is developed to account for this finding. With this model we can indeed explain the $V\propto Pe^{1/2}$ scaling, including the prefactor, which can collapse all data of the &#34;porous&#34; channels onto one universal curve, irrespective of channel geometry and composition of the mixtures. Our work provides a macroscopic approach to this bottom-up method of droplet generation and may guide further studies on oversaturation and nucleation in ternary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bubble puzzles: From fundamentals to applications

For centuries, bubbles have fascinated artists, engineers, and scientists alike. In spite of century-long research on them, new and often surprising bubble phenomena, features, and applications keep popping up. In this paper I sketch my personal scientific bubble journey, starting with single bubble sonoluminescence, continuing with sound emission and scattering of bubbles, cavitation, snapping shrimp, impact events, air entrainment, surface micro- and nanobubbles, and finally coming to effective force models for bubbles and dispersed bubbly two-phase flow. In particular, I also cover various applications of bubbles, namely in ultrasound diagnostics, drug and gene delivery, piezo-acoustic inkjet printing, immersion lithography, sonochemistry, electrolysis, catalysis, acoustic marine geophysical survey, and bubble drag reduction for naval vessels, and show how these applications crossed my way. I also try to show that good and interesting fundamental science and relevant applications are not a contradiction, but mutually stimulate each other in both directions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bubbly and Buoyant Particle-Laden Turbulent Flows

Fluid turbulence is commonly associated with stronger drag, greater heat transfer, and more efficient mixing than in laminar flows. In many natural and industrial settings, turbulent liquid flows contain suspensions of dispersed bubbles and light particles. Recently, much attention has been devoted to understanding the behavior and underlying physics of such flows by use of both experiments and high-resolution direct numerical simulations. This review summarizes our present understanding of various phenomenological aspects of bubbly and buoyant particle-laden turbulent flows. We begin by discussing different dynamical regimes, including those of crossing trajectories and wake-induced oscillations of rising particles, and regimes in which bubbles and particles preferentially accumulate near walls or within vortical structures. We then address how certain paradigmatic turbulent flows, such as homogeneous isotropic turbulence, channel flow, Taylor-Couette turbulence, and thermally driven turbulence, are modified by the presence of these dispersed bubbles and buoyant particles. We end with a list of summary points and future research questions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bubbly drag reduction using a hydrophobic inner cylinder in Taylor-Couette turbulence

In this study we experimentally investigate bubbly drag reduction in a highly turbulent flow of water with dispersed air at $5.0 \times 10^{5} \leq \text{Re} \leq 1.7 \times 10^{6}$ over a non-wetting surface containing micro-scale roughness. To do so, the Taylor-Couette geometry is used, allowing for both accurate global drag and local flow measurements. The inner cylinder - coated with a rough, hydrophobic material - is rotating, whereas the smooth outer cylinder is kept stationary. The crucial control parameter is the air volume fraction $α$ present in the working fluid. For small volume fractions ($α< {4}\,\%$), we observe that the surface roughness from the coating increases the drag. For large volume fractions of air ($α\geq 4\,\%$), the drag decreases compared to the case with both the inner and outer cylinders uncoated, i.e. smooth and hydrophilic, using the same volume fraction of air. This suggests that two competing mechanisms are at place: on the one hand the roughness invokes an extension of the log-layer - resulting in an increase in drag - and on the other hand there is a drag-reducing mechanism of the hydrophobic surface interacting with the bubbly liquid. The balance between these two effects determines whether there is overall drag reduction or drag enhancement. For further increased bubble concentration $α= {6}\,\%$ we find a saturation of the drag reduction effect. Our study gives guidelines for industrial applications of bubbly drag reduction in hydrophobic wall-bounded turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Catastrophic phase inversion in high-Reynolds number turbulent Taylor--Couette flow

Emulsions are omnipresent in the food industry, health care, and chemical synthesis. In this Letter the dynamics of meta-stable oil-water emulsions in highly turbulent ($10^{11}\leq\text{Ta}\leq 3\times 10^{13}$) Taylor--Couette flow, far from equilibrium, is investigated. By varying the oil-in-water void fraction, catastrophic phase inversion between oil-in-water and water-in-oil emulsions can be triggered, changing the morphology, including droplet sizes, and rheological properties of the mixture, dramatically. The manifestation of these different states is exemplified by combining global torque measurements and local in-situ laser induced fluorescence (LIF) microscopy imaging. Despite the turbulent state of the flow and the dynamic equilibrium of the oil-water mixture, the global torque response of the system is found to be as if the fluid were Newtonian, and the effective viscosity of the mixture was found to be several times bigger or smaller than either of its constituents.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Competing Marangoni and Rayleigh convection in evaporating binary droplets

For a small sessile or pendant droplet it is generally assumed that gravity does not play any role once the Bond number is small. This is even assumed for evaporating binary sessile or pendant droplets, in which convective flows can be driven due to selective evaporation of one component and the resulting concentration and thus surface tension differences at the air-liquid interface. However, recent studies have shown that in such droplets gravity indeed can play a role and that natural convection can be the dominant driving mechanism for the flow inside evaporating binary droplets (Edwards et al., Phys. Rev. Lett. 121, 184501 (2018); Li et al., Phys. Rev. Lett. 122, 114501 (2019)). In this study, we derive and validate a quasi-stationary model for the flow inside evaporating binary sessile and pendant droplets, which successfully allows to predict the prevalence and the intriguing interaction of Rayleigh and/or Marangoni convection on the basis of a phase diagram for the flow field expressed in terms of the Rayleigh and Marangoni numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Diffusive interaction of multiple surface nanobubbles and nanodroplets: shrinkage, growth, and coarsening

Surface nanobubbles are nanoscopic spherical-cap shaped gaseous domains on immersed substrates which are stable, even for days. After the stability of a {\it single} surface nanobubble has been theoretically explained, i.e. contact line pinning and gas oversaturation are required to stabilize them against diffusive dissolution [Lohse and Zhang, Phys.\ Rev.\ E 91, 031003 (R) (2015)], here we focus on the {\it collective} diffusive interaction of {\it multiple} nanobubbles. For that purpose we develop a finite difference scheme for the diffusion equation with the appropriate boundary conditions and with the immersed boundary method used to represent the growing or shrinking bubbles. After validation of the scheme against the exact results of Epstein and Plesset for a bulk bubble [J. Chem. Phys. 18, 1505 (1950)] and of Lohse and Zhang for a surface bubble, the framework of these simulations is used to describe the coarsening process of competitively growing nanobubbles. The coarsening process for such diffusively interacting nanobubbles slows down with advancing time and thus increasing bubble distance. The present results for surface nanobubbles are also applicable for immersed surface nanodroplets, for which better controlled experimental results of the coarsening process exist.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Double maxima of angular momentum transport in $η=0.91$ TC turbulence

We use experiments and direct numerical simulations to probe the phase-space of low-curvature Taylor--Couette (TC) flow in the vicinity of the ultimate regime. The cylinder radius ratio is fixed at $η=r_i/r_o=0.91$. Non-dimensional shear drivings (Taylor numbers $\text{Ta}$) in the range $10^7\leq\text{Ta}\leq10^{11}$ are explored for both co- and counter-rotating configurations. In the $\text{Ta}$ range $10^8\leq\text{Ta}\leq10^{10}$, we observe two local maxima of the angular momentum transport as a function of the cylinder rotation ratio, which can be described as either as &#34;co-&#34; and &#34;counter-rotating&#34; due to their location or as &#34;broad&#34; and &#34;narrow&#34; due to their shape. We confirm that the broad peak is accompanied by the strengthening of the large-scale structures, and that the narrow peak appears once the driving (Ta) is strong enough. As first evidenced in numerical simulations by Brauckmann \emph{et al.}~(2016), the broad peak is produced by centrifugal instabilities and that the narrow peak is a consequence of shear instabilities. We describe how the peaks change with $\text{Ta}$ as the flow becomes more turbulent. Close to the transition to the ultimate regime when the boundary layers (BLs) become turbulent, the usual structure of counter-rotating Taylor vortex pairs breaks down and stable unpaired rolls appear locally. We attribute this state to changes in the underlying roll characteristics during the transition to the ultimate regime. Further changes in the flow structure around $\text{Ta}\approx10^{10}$ cause the broad peak to disappear completely and the narrow peak to move. This second transition is caused when the regions inside the BLs which are locally smooth regions disappear and the whole boundary layer becomes active.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effect of axially varying sandpaper roughness on bubbly drag reduction in Taylor-Couette turbulence

We experimentally investigate the influence of alternating rough and smooth walls on bubbly drag reduction (DR). We apply rough sandpaper bands of width $s$ between $48.4\,mm$ and $148.5\,mm$, and roughness height $k = 695\,μm$, around the smooth inner cylinder (IC) of the Twente Turbulent Taylor-Couette facility. Between sandpaper bands, the IC is left uncovered over similar width $s$, resulting in alternating rough and smooth bands, a constant pattern in axial direction. We measure the DR in water that originates from introducing air bubbles to the fluid at (shear) Reynolds numbers $\textit{Re}_s$ ranging from $0.5 \times 10^6$ to $1.8 \times 10^6$. Results are compared to bubbly DR measurements with a completely smooth IC and an IC that is completely covered with sandpaper of the same roughness $k$. The outer cylinder is left smooth for all variations. Results are also compared to bubbly DR measurements where a smooth outer cylinder is rotating in opposite direction to the smooth IC. This counter rotation induces secondary flow structures that are very similar to those observed when the IC is composed of alternating rough and smooth bands. For the measurements with roughness, the bubbly DR is found to initially increase more strongly with $\textit{Re}_s$, before levelling off to reach a value that no longer depends on $\textit{Re}_s$. This is attributed to a more even axial distribution of the air bubbles, resulting from the increased turbulence intensity of the flow compared to flow over a completely smooth wall at the same $\textit{Re}_s$. The air bubbles are seen to accumulate at the rough wall sections in the flow. Here, locally, the drag is largest and so the drag reducing effect of the bubbles is felt strongest. Therefore, a larger maximum value of bubbly DR is found for the alternating rough and smooth walls compared to the completely rough wall.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electrochemically Induced pH Change: Time-Resolved Confocal Fluorescence Microscopy Measurements and Comparison with Numerical Model

Confocal fluorescence microscopy is a proven technique, which can image near-electrode pH changes. For a complete understanding of electrode processes, time-resolved measurements are required, which have not yet been provided. Here we present the first measurements of time-resolved pH profiles with confocal fluorescence microscopy. The experimental results compare favorably with a one-dimensional reaction-diffusion model; this holds up to the point where the measurements reveal three-dimensionality in the pH distribution. Specific factors affecting the pH measurement such as attenuation of light and the role of dye migration are also discussed in detail. The method is further applied to reveal the buffer effects observed in sulfate-containing electrolytes. The work presented here is paving the way toward the use of confocal fluorescence microscopy in the measurement of 3D time-resolved pH changes in numerous electrochemical settings, for example in the vicinity of bubbles.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Entrapment and Dissolution of Microbubbles Inside Microwells

The formation and evolution of immersed surface micro- and nanobubbles are essential in various practical applications, such as the usage of superhydrophobic rematerials, drug delivery, and mineral flotation. In this work, we investigate the entrapment of microbubbles on a hydrophobic surface, structured with microwells, when water flow passes along, and the subsequent microbubble dissolution. At entrapment, the microbubble is initially pinned at the edge of the microwell. At some point, the three-phase contact line detaches from one side of the edge and separates from the wall, after which it further recedes. We systematically investigate the evolution of the footprint diameter and the contact angle of the entrapped microbubbles, which reveals that the dissolution process is in the constant contact angle mode. By varying the gas undersaturation level, we quantify how a high gas undersaturation enhances the dissolution process, and compare with simplified theoretical predictions for dissolving bubbles on a plane surface. We find that geometric partial blockage effects of the diffusive flux out of the microbubble trapped in the microwell lead to reduced dissolution rates.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Evaporation-Induced Crystallization of Surfactants in Sessile Multicomponent Droplets

Surfactants have been widely studied and used in controlling droplet evaporation. In this work, we observe and study the crystallization of sodium dodecyl sulfate (SDS) within an evaporating glycerol-water mixture droplet. The crystallization is induced by the preferential evaporation of water, which decreases the solubility of SDS in the mixture. As a consequence, the crystals shield the droplet surface and cease the evaporation. The universality of the evaporation characteristics for a range of droplet sizes is revealed by applying a diffusion model, extended by Raoult&#39;s law. To describe the nucleation and growth of the crystals, we employ the 2-dimensional crystallization model of Weinberg [J. Non-Cryst. Solids 1991, 134, 116]. The results of this model compare favorably to our experimental results. Our findings may inspire the community to reconsider the role of high concentration of surfactants in a multicomponent evaporation system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fast-freezing kinetics inside a droplet impacting on a cold surface

Freezing or solidification of impacting droplets is omnipresent in nature and technology, be it a rain droplet falling on a supercooled surface, be it in inkjet printing where often molten wax is used, be it in added manufacturing or in metal production processes or in extreme ultraviolet lithography (EUV) for the chip production where molten tin is used to generate the EUV radiation. For many of these industrial applications, a detailed understanding of the solidification process is essential. Here, by adopting a totally new optical technique in the context of freezing, namely TIR (Total-Internal-Reflection), we elucidate the freezing kinetics during the solidification of a droplet while it impacts on an undercooled surface. We show for the first time that at sufficiently high undercooling a peculiar freezing morphology exists that involves sequential advection of frozen fronts from the centre of the droplet to its boundaries. This phenomenon is examined by combining elements of classical nucleation theory to the large scale hydrodynamics on the droplet scale, bringing together two subfields which traditionally have been quite separated. Furthermore, we report a peculiar self-peeling phenomenon of a frozen splat that is driven by the existence of a transient crystalline state during solidification.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Flow organization in laterally unconfined Rayleigh-Bénard turbulence

We investigate the large-scale circulation (LSC) of turbulent Rayleigh-Bénard convection in a large box of aspect ratio $Γ=32$ for Rayleigh numbers up to $Ra=10^9$ and at a fixed Prandtl number $Pr=1$. A conditional averaging technique allows us to extract statistics of the LSC even though the number and the orientation of the structures vary throughout the domain. We find that various properties of the LSC obtained here, such as the wall-shear stress distribution, the boundary layer thicknesses and the wind Reynolds number, do not differ significantly from results in confined domains ($Γ\approx 1$). This is remarkable given that the size of the structures (as measured by the width of a single convection roll) more than doubles at the highest $Ra$ as the confinement is removed. An extrapolation towards the critical shear Reynolds number of $Re_s^{\textrm{crit}} \approx 420$, at which the boundary layer (BL) typically becomes turbulent, predicts that the transition to the ultimate regime is expected at $Ra_{\textrm{crit}} \approx \mathcal{O}(10^{15})$ in unconfined geometries. This result is in line with the Göttingen experimental observations. Furthermore, we confirm that the local heat transport close to the wall is highest in the plume impacting region, where the thermal BL is thinnest, and lowest in the plume emitting region, where the thermal BL is thickest. This trend, however, weakens with increasing $Ra$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

From Rayleigh-Bénard convection to porous-media convection: how porosity affects heat transfer and flow structure

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 $ϕ$ in the flow properties during the transition process from the traditional RB convection with $ϕ=1$ (so no obstacles included) to Darcy-type porous-media convection with $ϕ$ approaching 0. Simulations are carried out in a cell with unity aspect ratio, for the Rayleigh number $Ra$ from $10^5$ to $10^{10}$ and varying porosities $ϕ$, 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-monotonously as a function of $ϕ$; namely, with decreasing $ϕ$, it first increases, before it decreases for $ϕ$ approaching 0. The non-monotonous behaviour of $Nu(ϕ)$ 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 $ϕ$, depending on $Ra$, two different heat transfer regimes are identified, with different effective power-law behaviours of $Nu$ vs $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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

From zonal flow to convection rolls in Rayleigh-Bénard convection with free-slip plates

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-dimension (2D) RB convection and the other one three-dimension (3D) RB convection with a rotating axis parallel to the plate. We explore the parameter range of Rayleigh numbers Ra from $10^7 to $10^9$ and Prandtl numbers $Pr$ from $1$ to $100$. We show that zonal flow, which was observed, for example, by Goluskin \emph{et al}. \emph{J. Fluid. Mech.} 759, 360-385 (2014) for $Γ=2$, is only stable when $Γ$ is smaller than a critical value, which depends on $Ra$ and $Pr$. With increasing $Γ$, we find a second regime in which both zonal flow and different convection roll states can be statistically stable. For even larger $Γ$, in a third regime, only convection roll states are statistically stable and zonal flow is not sustained. For the 3D simulations, we fix $Ra=10^7$ and $Pr=0.71$, and compare the flow for $Γ=8$ and $Γ= 16$. We demonstrate that with increasing aspect ratio $Γ$, zonal flow, which was observed for small $Γ=2π$ by von Hardenberg \emph{et al}. \emph{Phys. Rev. Lett.} 15, 134501 (2015), completely disappears for $Γ=16$. For such large $Γ$ only convection roll states are statistically stable. In between, here for medium aspect ratio $Γ= 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 2D case.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gas-Vapor Interplay in Plasmonic Bubble Shrinkage

The understanding of the shrinkage dynamics of plasmonic bubbles formed around metallic nanoparticles immersed in liquid and irradiated by a resonant light source is crucial for the usage of these bubbles in numerous applications. In this paper we experimentally show and theoretically explain that a plasmonic bubble during its shrinkage undergoes two different phases: first, a rapid partial bubble shrinkage governed by vapor condensation and, second, a slow diffusion-controlled bubble dissolution. The history of the bubble formation plays an important role in the shrinkage dynamics during the first phase, as it determines the gas-vapor ratio in the bubble composition. Higher laser powers lead to more vaporous bubbles, while longer pulses and higher dissolved air concentrations lead to more gaseous bubbles. The dynamics of the second phase barely depends on the history of bubble formation, i.e. laser power and pulse duration, but strongly on the dissolved air concentration, which defines the concentration gradient at the bubble interface. Finally, for the bubble dissolution in the second phase, with decreasing dissolved air concentration, we observe a gradual transition from a $R(t) \propto (t_0 - t) ^{1/3}$ scaling law to a $R(t) \propto (t_0 - t) ^{1/2}$ scaling law, where $t_0$ is the lifetime of the bubble and theoretically explain this transition.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Giant plasmonic bubbles nucleation under different ambient pressures

Water-immersed gold nanoparticles irradiated by a laser can trigger the nucleation of plasmonic bubbles after a delay time of a few microseconds [Wang et al., Proc. Natl. Acad. Sci. USA 122, 9253,(2018)]. Here we systematically investigated the light-vapor conversion efficiency, η, of these plasmonic bubbles as a function of the ambient pressure. The efficiency of the formation of these initial-phase and mainly water-vapor containing bubbles, which is defined as the ratio of the energy that is required to form the vapor bubbles and the total energy dumped in the gold nanoparticles before nucleation of the bubble by the laser, can be as high as 25%. The amount of vaporized water first scales linearly with the total laser energy dumped in the gold nanoparticles before nucleation, but for larger energies the amount of vaporized water levels off. The efficiency η decreases with increasing ambient pressure. The experimental observations can be quantitatively understood within a theoretical framework based on the thermal diffusion equation and the thermal dynamics of the phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

In-air microfluidics enables rapid fabrication of emulsions, suspensions, and 3D modular (bio)materials

Microfluidic chips provide unparalleled control over droplets and jets, which have advanced all natural sciences. However, microfluidic applications could be vastly expanded by increasing the per-channel throughput and directly exploiting the output of chips for rapid additive manufacturing. We unlock these features with in-air microfluidics, a new chip-free platform to manipulate microscale liquid streams in the air. By controlling the composition and in-air impact of liquid microjets by surface tension-driven encapsulation, we fabricate monodisperse emulsions, particles, and fibers with diameters of 20 to 300um at rates that are 10 to 100 times higher than chip-based droplet microfluidics. Furthermore, in-air microfluidics uniquely enables module-based production of three-dimensional (3D) multiscale (bio)materials in one step because droplets are partially solidified in-flight and can immediately be printed onto a substrate. In-air microfluidics is cytocompatible, as demonstrated by additive manufacturing of 3D modular constructs with tailored micro-environments for multiple cell types. Its in-line control, high throughput and resolution, and cytocompatibility make in-air microfluidics a versatile platform technology for science, industry, and health care.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Initial solidification dynamics of spreading droplets

When a droplet is brought in contact with an undercooled surface, it wets the substrate and solidifies at the same time. The interplay between the phase transition effects and the contact-line motion, leading to its arrest, remains poorly understood. Here we reveal the early solidification patterns and dynamics of spreading hexadecane droplets. Total internal reflection (TIR) imaging is employed to temporally and spatially resolve the early solidification behaviour. With this, we determine the conditions leading to the contact-line arrest. We quantify the overall nucleation behaviour, \textit{i.e.} the nucleation rate and the crystal growth speed, and show its sensitivity to the applied undercooling of the substrate. By combining the Johnson-Mehl-Avrami-Kolmogorov nucleation theory and scaling relations for the spreading, we can calculate the temporal evolution of the solid area fraction, which is in good agreement with our observations. We also show that for strong enough undercooling it is the rapid growth of the crystals which determines the eventual arrest of the spreading contact line.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modelling large scale airgun-bubble dynamics with highly non-spherical features

A thorough understanding of the dynamics of meter-sized airgun-bubbles is very crucial to seabed geophysical exploration. In this study, we use the boundary integral method to investigate the highly non-spherical airgun-bubble dynamics and its corresponding pressure wave emission. Moreover, a model is proposed to also consider the process of air release from the airgun port, which is found to be the most crucial factor to estimate the initial peak of the pressure wave. The numerical simulations show good agreement with experiments, in terms of non-spherical bubble shapes and pressure waves. Thereafter, the effects of the port opening time $T\rm_{open}$, airgun firing depth, heat transfer, and gravity are numerically investigated. We find that a smaller $T\rm_{open}$ leads to a more violent air release that consequently causes stronger high-frequency pressure wave emissions; however, the low-frequency pressure waves are little affected. Additionally, the non-spherical bubble dynamics is highly dependent on the Froude number $Fr$. Starting from $Fr=2$, as $Fr$ increases, the jet contains lower kinetic energy, resulting in a stronger energy focusing of the bubble collapse itself and thus a larger pressure peak during the bubble collapse phase. For $Fr \ge 7$, the spherical bubble theory becomes an appropriate description of the airgun-bubble. The new findings of this study may provide a reference for practical operations and designing environmentally friendly airguns in the near future.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Morphology of Evaporating Sessile Microdroplets on Lyophilic Elliptical Patches

The evaporation of droplets occurs in a large variety of natural and technological processes such as medical diagnostics, agriculture, food industry, printing, and catalytic reactions. We study the different droplet morphologies adopted by an evaporating droplet on a surface with an elliptical patch with a different contact angle. We perform experiments to observe these morphologies and use numerical calculations to predict the effects of the patched surfaces. We observe that tuning the geometry of the patches offers control over the shape of the droplet. In the experiments, the drops of various volumes are placed on elliptical chemical patches of different aspect ratios and imaged in 3D using laser scanning confocal microscopy, extracting the droplets shape. In the corresponding numerical simulations, we minimize the interfacial free energy of the droplet, by employing Surface Evolver. The numerical results are in good qualitative agreement with our experimental data and can be used for the design of micropatterned structures, potentially suggesting or excluding certain morphologies for particular applications. However, the experimental results show the effects of pinning and contact angle hysteresis, which are obviously absent in the numerical energy minimization. The work culminates with a morphology diagram in the aspect ratio vs relative volume parameter space, comparing the predictions with the measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multiple states in turbulent large-aspect ratio thermal convection: What determines the number of convection rolls?

Recent findings suggest that wall-bounded turbulent flow can take different statistically stationary turbulent states, with different transport properties, even for the very same values of the control parameters. What state the system takes depends on the initial conditions. Here we analyze the multiple states in large-aspect ratio ($Γ$) two-dimensional turbulent Rayleigh--Bénard flow with no-slip plates and horizontally periodic boundary conditions as model system. We determine the number $n$ of convection rolls, their mean aspect ratios $Γ_r = Γ/n$, and the corresponding transport properties of the flow (i.e., the Nusselt number $Nu$), as function of the control parameters Rayleigh ($Ra$) and Prandtl number. The effective scaling exponent $β$ in $Nu \sim Ra^β$ is found to depend on the realized state and thus $Γ_r$, with a larger value for the smaller $Γ_r$. By making use of a generalized Friedrichs inequality, we show that the elliptical instability and viscous damping determine the $Γ_r$-window for the realizable turbulent states. The theoretical results are in excellent agreement with our numerical finding $2/3 \le Γ_r \le 4/3$, where the lower threshold is approached for the larger $Ra$. Finally, we show that the theoretical approach to frame $Γ_r$ also works for free-slip boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Physicochemical Hydrodynamics of Droplets out of Equilibrium: A Perspective Review

Droplets abound in nature and technology. In general, they are multicomponent, and, when out of equilibrium, with gradients in concentration, implying flow and mass transport. Moreover, phase transitions can occur, with either evaporation, solidification, dissolution, or nucleation of a new phase. The droplets and their surrounding liquid can be binary, ternary, or contain even more components, and with several even in different phases. In the last two decades the rapid advances in experimental and numerical fluid dynamical techniques has enabled major progress in our understanding of the physicochemical hydrodynamics of such droplets, further narrowing the gap from fluid dynamics to chemical engineering and colloid & interfacial science and arriving at a quantitative understanding of multicomponent and multiphase droplet systems far from equilibrium, and aiming towards a one-to-one comparison between experiments and theory or numerics. This review will discuss various examples of the physicochemical hydrodynamics of droplet systems far from equilibrium and our present understanding of them. These include immiscible droplets in a concentration gradient, coalescence of droplets of different liquids, droplets in concentration gradients emerging from chemical reactions (including droplets as microswimmers) and phase transitions such as evaporation, solidification, dissolution, or nucleation, and droplets in ternary liquids, including solvent exchange, nano-precipitation, and the so-called ouzo effect. We will also discuss the relevance of the physicochemical hydrodynamics of such droplet systems for many important applications, including in chemical analysis and diagnostics, microanalysis, pharmaceutics, synthetic chemistry and biology, chemical and environmental engineering, the oil and remediation industries, inkjet-printing, for micro- and nano-materials, and in nanotechnolgoy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Plasmonic Bubble Nucleation and Growth in Water: Effect of Dissolved Air

Under continuous laser irradiation, noble metal nanoparticles immersed in water can quickly heat up, leading to the nucleation of so-called plasmonic bubbles. In this work, we want to further understand the bubble nucleation and growth mechanism. In particular, we quantitatively study the effect of the amount of dissolved air on the bubble nucleation and growth dynamics, both for the initial giant bubble, which forms shortly after switching on the laser and is mainly composed of vapor, and for the final life phase of the bubble, during which it mainly contains air expelled from water. We found that the bubble nucleation temperature depends on the gas concentration: the higher the gas concentration, the lower the bubble nucleation temperature. Also, the long-term diffusiondominated bubble growth is governed by the gas concentration. The radius of the bubbles grows as R(t)~t^1/3 for airequilibrated and air-oversaturated water. In contrast, in partially degassed water, the growth is much slower since, even for the highest temperature we achieve, the water remains undersaturated.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rayleigh-Taylor instability by segregation in an evaporating multi-component microdroplet

The evaporation of multi-component droplets is relevant to various applications but challenging to study due to the complex physicochemical dynamics. Recently, Li (2018) reported evaporation-triggered segregation in 1,2-hexanediol-water binary droplets. In this present work, we added 0.5 wt% silicone oil into the 1,2-hexanediol-water binary solution. This minute silicone oil concentration dramatically modifies the evaporation process as it triggers an early extraction of the 1,2-hexanediol from the mixture. Surprisingly, we observe that the segregation of 1,2-hexanediol forms plumes, rising up from the rim of the sessile droplet towards the apex during the droplet evaporation. By orientating the droplet upside down, i.e., by studying a pendant droplet, the absence of the plumes indicates that the flow structure is induced by buoyancy, which drives a Rayleigh-Taylor instability (i.e., driven by density differences & gravitational acceleration). From micro-PIV measurement, we further prove that the segregation of the non-volatile component (1,2-hexanediol) hinders the evaporation near the contact line, which leads to a suppression of the Marangoni flow in this region. Hence, on long time scales, gravitational effects play the dominant role in the flow structure, rather than Marangoni flows. We compare the measurement of the evaporation rate with the diffusion model of Popov (2005), coupled with Raoult&#39;s law and the activity coefficient. This comparison indeed confirms that the silicone-oil-triggered segregation of the non-volatile 1,2-hexanediol significantly delays the evaporation. With an extended diffusion model, in which the influence of the segregation has been implemented, the evaporation can be well described.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Scaling in internally heated convection: a unifying theory

We offer a unifying theory for turbulent purely internally heated convection, generalizing the unifying theories of Grossmann and Lohse (2000, 2001) for Rayleigh--Bénard turbulence and of Shishkina, Grossmann and Lohse (2016) for turbulent horizontal convection, which are both based on the splitting of the kinetic and thermal dissipation rates in respective boundary and bulk contributions. We obtain the mean temperature of the system and the Reynolds number (which are the response parameters) as function of the control parameters, namely the internal thermal driving strength (called, when nondimensionalized, the Rayleigh--Roberts number) and the Prandtl number. The results of the theory are consistent with our direct numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Self-propelling Microdroplets Generated and Sustained by Liquid-liquid Phase Separation in Confined Spaces

Flow transport in confined spaces is ubiquitous in technological processes, ranging from separation and purification of pharmaceutical ingredients by microporous membranes and drug delivery in biomedical treatment to chemical and biomass conversion in catalyst-packed reactors and carbon dioxide sequestration. In this work, we suggest a distinct pathway for enhanced liquid transport in a confined space via self-propelling microdroplets. These microdroplets can form spontaneously from localized liquid-liquid phase separation as a ternary mixture is diluted by a diffusing poor solvent. High speed images reveal how the microdroplets grow, break up and propel rapidly along the solid surface, with a maximal velocity up to ~160 um/s, in response to a sharp concentration gradient resulting from phase separation. The microdroplet self-propulsion induces a replenishing flow between the walls of the confined space towards the location of phase separation, which in turn drives the mixture out of equilibrium and leads to a repeating cascade of events. Our findings on the complex and rich phenomena of self-propelling droplets suggest an effective approach to enhanced flow motion of multicomponent liquid mixtures within confined spaces for time effective separation and smart transport processes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Solvent Exchange in a Hele-Shaw Cell Universality of Surface Nanodroplet Nucleation

Solvent exchange (also called solvent shifting or Ouzo effect) is a generally used bottom-up process to mass-produce nanoscale droplets. In this process, a good solvent for some oil is displaced by a poor one, leading to oil nanodroplet nucleation and subsequent growth. Here we perform this process on a hydrophobic substrate so that sessile droplets so-called surface nanodroplets-develop, following the work of Zhang et al. [Zhang, X.; Lu, Z.; Tan, H.; Bao, L.; He, Y.; Sun, C.; Lohse, D. Proc. Natl. Acad. Sci. U.S.A. 2015, 122, 9253-9257]. In contrast to what was done in that paper, we chose a very well-controlled Hele-Shaw geometry with negligible gravitational effects, injecting the poor solvent in the center of the Hele-Shaw cell, and characterize the emerging nanodroplets as a function of radial distance and flow rates. We find that the mean droplet volume per area <Vol>_area strongly depends on the local Peclet number Pe and follows a universal scaling law <Vol>_area~Pe^(3/4). Moreover, the probability distribution function of the droplet volume strongly depends on the local Pe as well, regardless of the flow rates and radial distance, giving strong support to the theoretical model of the solvent exchange process developed in Zhang et al.&#39;s work.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Speeding up biphasic reactions with surface nanodroplets

Biphasic chemical reactions compartmentalized in small droplets offer advantages, such as streamlined procedures for chemical analysis, enhanced chemical reaction efficiency and high specificity of conversion. In this work, we experimentally and theoretically investigate the rate for biphasic chemical reactions between acidic nanodroplets on a substrate surface and basic reactants in a surrounding bulk flow. The reaction rate is measured by droplet shrinkage as the product is removed from the droplets by the flow. In our experiments, we determine the dependence of the reaction rate on the flow rate and the solution concentration. The theoretical analysis predicts that the life time $τ$ of the droplets scales with Peclet number $Pe$ and the reactant concentration in the bulk flow $c_{re,bulk}$ as $τ\propto Pe^{-3/2}c_{re,bulk}^{-1}$, in good agreement with our experimental results. Furthermore, we found that the product from the reaction on an upstream surface can postpone the droplet reaction on a downstream surface, possibly due to the adsorption of interface-active products on the droplets in the downstream. The time of the delay decreases with increasing $Pe$ of the flow and also with increasing reactant concentration in the flow, following the scaling same as that of the reaction rate with these two parameters. Our findings provide insight for the ultimate aim to enhance droplet reactions under flow conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Splitting droplet through coalescence of two different three-phase contact lines

Moving contact lines of more than two phases dictate a large number of interfacial phenomena. Despite its significance to fundamental and applied processes, the contact lines at a junction of four-phases (two immiscible liquids, solid and gas) have been addressed only in a few investigations. Here, we report an intriguing phenomenon that follows after the four phases of oil, water, solid and gas make contact through the coalescence of two different three-phase contact lines. We combine experimental study and theoretical analysis to reveal and rationalize the dynamics exhibited upon the coalescence between the contact line of a micron-sized oil droplet and the receding contact line of a millimetre-sized water drop that covers the oil droplet on the substrate. We find that after the coalescence a four-phase contact line is formed for a brief period. However, this quadruple contact line is not stable, leading to a `droplet splitting&#39; effect and eventual expulsion of the oil droplet from the water drop. We then show that the interfacial tension between the different phases and the viscosity of oil droplet dictate the splitting dynamics. More viscous oils display higher resistance to the extreme deformations of the droplet induced by the instability of the quadruple contact line and no droplet expulsion is observed for such cases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Subcritical behaviour in double diffusive convection within the diffusive regime

We conduct two- and three-dimensional simulations for double diffusive convection in the diffusive regime, where the fluid flow is driven by a destabilizing temperature gradient and stabilized by a stably stratified salinity gradient. We study how the heat flux, Reynolds number, and flow structures change with the density ratio $Λ$, which is the ratio of the buoyancy force induced by the salinity gradient to that by the temperature gradient. When $Λ$ increases from zero, the flow first behaves similarly as in pure Rayleigh-Bénard (RB) convection, both with respect to flow structure and to heat transport. The linear stability analysis of Baines & Gill (J. Fluid Mech., vol. 37, 1969, pp. 289-306) had estimated the critical density ratio $Λ_c$, above which the flow becomes stable. However, here we show that by using a large-scale circulation as initial condition (rather than the linear profiles assumed in the linear stability analysis), DDC in the diffusive regime can exhibit subcritical behaviour when $Λ> Λ_c$, i.e., coexistence of states at the same control parameters. Even though the density ratio becomes thousands times that of the critical value $Λ_c$, there is still convection with strongly enhanced heat transfer properties compared to the pure conduction case. We reveal the corresponding flow structures and find an unstably-stratified region sandwiched between two stably-stratified layers. Our results demonstrate the importance of the initial condition for DDC in the diffusive regime, especially in the situation of a large density ratio, which occurs in high-latitude ocean regions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The effect of Prandtl number on turbulent sheared thermal convection

In turbulent wall sheared thermal convection, there are three different flow regimes, depending on the relative relevance of thermal forcing and wall shear. In this paper we report the results of direct numerical simulations of such sheared Rayleigh-Bénard convection, at fixed Rayleigh number $Ra=10^6$, varying the wall Reynolds number in the range $0 \leq Re_w \leq 4000$ and Prandtl number $0.22 \leq Pr \leq 4.6$, extending our prior work by Blass et al. (2020), where $Pr$ was kept constant at unity and the thermal forcing ($Ra$) varied. We cover a wide span of bulk Richardson numbers $0.014 \leq Ri \leq 100$ and show that the Prandtl number strongly influences the morphology and dynamics of the flow structures. In particular, at fixed $Ra$ and $Re_w$, a high Prandtl number causes stronger momentum transport from the walls and therefore yields a greater impact of the wall shear on the flow structures, resulting in an increased effect of $Re_w$ on the Nusselt number. Furthermore, we analyse the thermal and kinetic boundary layer thicknesses and relate their behaviour to the resulting flow regimes. For the largest shear rates and $Pr$ numbers, we observe the emergence of a Prandtl-von Karman log-layer, signalling the onset of turbulent dynamics in the boundary layer. Finally, our results allow to extend the Grossmann-Lohse theory for heat transport in Rayleigh-Bénard convection to the sheared case, universally describing $Nu(Ra,Pr,Re_w)$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Time-Resolved In Situ Liquid-Phase Atomic Force Microscopy and Infrared Nanospectroscopy during the Formation of Metal-Organic Framework Thin Films

Metal-organic framework (MOF) thin-films show unmatched promise as smart membranes and photo-catalytic coatings. However, their nucleation and growth resulting from intricate molecular assembly processes are not well understood, yet crucial to control the thin-film properties. Here, we directly observe the nucleation and growth behaviour of HKUST-1 thin-films by real-time in-situ AFM at different temperatures in the Cu-BTC solution. In combination with ex-situ infrared (nano-)spectroscopy, synthesis at 25°C reveals initial nucleation of rapidly growing HKUST-1 islands surrounded by a continuously nucleating, but slowly growing HKUST-1 carpet. Monitoring at 13°C and 50°C shows the strong impact of temperature on thin-film formation, resulting in (partial) nucleation and growth inhibition. The nucleation and growth mechanisms as well as their kinetics provide insights to aid in future rational design of MOF thin-films.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Two-layer Thermally Driven Turbulence: Mechanisms for Interface Breakup

It is commonly accepted that the breakup criteria of drops or bubbles in turbulence is governed by surface tension and inertia. However, also {\it{buoyancy}} can play an important role at breakup. In order to better understand this role, here we numerically study Rayleigh-Bénard convection for two immiscible fluid layers, in order to identify the effects of buoyancy on interface breakup. We explore the parameter space spanned by the Weber number $5\leq We \leq 5000$ (the ratio of inertia to surface tension) and the density ratio between the two fluids $0.001 \leq Λ\leq 1$, at fixed Rayleigh number $Ra=10^8$ and Prandtl number $Pr=1$. At low $We$, the interface undulates due to plumes. When $We$ is larger than a critical value, the interface eventually breaks up. Depending on $Λ$, two breakup types are observed: The first type occurs at small $Λ\ll 1$ (e.g. air-water systems) when local filament thicknesses exceed the Hinze length scale. The second, strikingly different, type occurs at large $Λ$ with roughly $0.5 < Λ\le 1$ (e.g. oil-water systems): The layers undergo a periodic overturning caused by buoyancy overwhelming surface tension. For both types the breakup criteria can be derived from force balance arguments and show good agreement with the numerical results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ultrasound-enhanced mass transfer during single-bubble diffusive growth

Ultrasound is known to enhance surface bubble growth and removal in catalytic and microfluidic applications, yet the contributions of rectified diffusion and microstreaming phenomena towards mass transfer remain unclear. We quantify the effect of ultrasound on the diffusive growth of a single spherical CO$_2$ bubble growing on a substrate in supersaturated water. The time dependent bubble size, shape, oscillation amplitude and microstreaming flow field are resolved. We show and explain how ultrasound can enhance the diffusive growth of surface bubbles by up to two orders of magnitude during volumetric resonance. The proximity of the wall forces the bubble to oscillate non-spherically, thereby generating vigorous streaming during resonance that results in convection-dominated growth.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Water entry of spheres into a rotating liquid

The transient cavity dynamics during water entry of a heavy, non-rotating sphere impacting a rotating pool of liquid is studied experimentally, numerically, and theoretically. We show that the pool rotation advances the transition of the cavity type - from deep seal to surface seal - marked by a reduction in the transitional Froude number. The role of the dimensionless rotational number $\mathcal{S} \equiv ωR_0/U_0$ on the transient cavity dynamics is unveiled, where $R_0$ is the sphere radius, $ω$ the angular speed of the liquid, and $U_0$ the impact velocity. The rotating background liquid has two discernible effects on the cavity evolution. Firstly, an increase in the underwater pressure field due to centripetal effects, and secondly a reduction in the pressure of airflow in the cavity neck near the water surface. The non-dimensional pinch-off time of the deep seal shows a robust 1/2 power-law dependence on the Froude number, but with a reducing prefactor for increasing $ω$. Our findings reveal that the effects of a rotating background liquid on the water entry can be traced back to the subtle differences in the initial stage splash and the near-surface cavity dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

What rotation rate maximizes heat transport in rotating Rayleigh-Bénard convection with Prandtl number larger than one?

The heat transfer and flow structure in rotating Rayleigh-Bénard convection are strongly influenced by the Rayleigh ($Ra$), Prandtl ($Pr$), and Rossby ($Ro$) number. For $Pr\gtrsim 1$ and intermediate rotation rates, the heat transfer is increased compared to the non-rotating case. We find that the regime of increased heat transfer is subdivided into a low and a high $Ra$ number regime. For $Ra\lesssim 5\times10^8$ the heat transfer at a given $Ra$ and $Pr$ is highest at an optimal rotation rate, at which the thickness of the viscous and thermal boundary layer is about equal. From the scaling relations of the thermal and viscous boundary layer thicknesses, we derive that the optimal rotation rate scales as $1/Ro_\mathrm{opt} \approx 0.12 Pr^{1/2}Ra^{1/6}$. In the low $Ra$ regime the heat transfer is similar in a periodic domain and cylindrical cells with different aspect ratios, i.e.\ the ratio of diameter to height. This is consistent with the view that the vertically aligned vortices are the dominant flow structure. For $Ra\gtrsim 5\times10^8$ the above scaling for the optimal rotation rate does not hold anymore. It turns out that in the high $Ra$ regime, the flow structures at the optimal rotation rate are very different than for lower $Ra$. Surprisingly, the heat transfer in the high $Ra$ regime differs significantly for a periodic domain and cylindrical cells with different aspect ratios, which originates from the sidewall boundary layer dynamics and the corresponding secondary circulation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Coherence of temperature and velocity superstructures in turbulent Rayleigh-Bénard flow

We investigate the interplay between large-scale patterns, so-called superstructures, in the fluctuation fields of temperature $θ$ and vertical velocity $w$ in turbulent Rayleigh-Bénard convection at large aspect ratios. Earlier studies suggested that velocity superstructures were smaller than their thermal counterparts in the center of the domain. However, a scale-by-scale analysis of the correlation between the two fields employing the linear coherence spectrum reveals that superstructures of the same size exist in both fields, which are almost perfectly correlated. The issue is further clarified by the observation that in contrast to the temperature, and unlike assumed previously, superstructures in the vertical velocity field do not result in a peak in the power spectrum of $w$. The origin of this difference is traced back to the production terms of the $θ$- and $w$-variance. These results are confirmed for a range of Rayleigh numbers $Ra = 10^5$--$10^9$, the superstructure size is seen to increase monotonically with $Ra$. Furthermore, the scale distribution of particularly the temperature fluctuations is pronouncedly bimodal. In addition to the large-scale peak caused by the superstructures, there exists a strong small-scale peak. This `inner peak&#39; is most intense at a distance of $δ_θ$ from the wall and associated with structures of size $\approx 10 δ_θ$, where $δ_θ$ is the thermal boundary layer thickness. Finally, based on the vertical coherence relative to a reference height of $δ_θ$, a self-similar structure is identified in the velocity field (vertical and horizontal components) but not in the temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Controlling secondary flow in Taylor-Couette turbulence through spanwise-varying roughness

Highly turbulent Taylor-Couette flow with spanwise-varying roughness is investigated experimentally and numerically (direct numerical simulations (DNS) with an immersed boundary method (IBM)) to determine the effects of the spacing and axial width $s$ of the spanwise varying roughness on the total drag and {on} the flow structures. We apply sandgrain roughness, in the form of alternating {rough and smooth} bands to the inner cylinder. Numerically, the Taylor number is $\mathcal{O}(10^9)$ and the roughness width is varied between $0.47\leq \tilde{s}=s/d \leq 1.23$, where $d$ is the gap width. Experimentally, we explore $\text{Ta}=\mathcal{O}(10^{12})$ and $0.61\leq \tilde s \leq 3.74$. For both approaches the radius ratio is fixed at $η=r_i/r_o = 0.716$, with $r_i$ and $r_o$ the radius of the inner and outer cylinder respectively. We present how the global transport properties and the local flow structures depend on the boundary conditions set by the roughness spacing $\tilde{s}$. Both numerically and experimentally, we find a maximum in the angular momentum transport as function of $\tilde s$. This can be atributed to the re-arrangement of the large-scale structures triggered by the presence of the rough stripes, leading to correspondingly large-scale turbulent vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Convection-dominated dissolution for single and multiple immersed sessile droplets

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 $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 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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Drag reduction in boiling Taylor-Couette turbulence

We create a highly controlled lab environment-accessible to both global and local monitoring-to analyse turbulent boiling flows and in particular their shear stress in a statistically stationary state. Namely, by precisely monitoring the drag of strongly turbulent Taylor-Couette flow (the flow in between two co-axially rotating cylinders, Reynolds number $\textrm{Re}\approx 10^6$) during its transition from non-boiling to boiling, we show that the intuitive expectation, namely that a few volume percent of vapor bubbles would correspondingly change the global drag by a few percent, is wrong. Rather, we find that for these conditions a dramatic global drag reduction of up to 45% occurs. We connect this global result to our local observations, showing that for major drag reduction the vapor bubble deformability is crucial, corresponding to Weber numbers larger than one. We compare our findings with those for turbulent flows with gas bubbles, which obey very different physics than vapor bubbles. Nonetheless, we find remarkable similarities and explain these.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Drop fragmentation by laser-pulse impact

We study the fragmentation of a liquid drop that is hit by a laser pulse. The drop expands into a thin sheet that breaks by the radial expulsion of ligaments from its rim and the nucleation and growth of holes on the sheet. By combining experimental data from two liquid systems with vastly different time- and length scales we show how the early-time laser-matter interaction affects the late-time fragmentation. We identify two Rayleigh--Taylor instabilities of different origins as the prime cause of the fragmentation and derive scaling laws for the characteristic breakup time and wavenumber. The final web of ligaments results from a subtle interplay between these instabilities and deterministic modulations of the local sheet thickness, which originate from the drop deformation dynamics and spatial variations in the laser-beam profile.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Flow organization and heat transfer in turbulent wall sheared thermal convection

We perform direct numerical simulations of wall sheared Rayleigh-Bénard (RB) convection for Rayleigh numbers up to $Ra=10^8$, Prandtl number unity, and wall shear Reynolds numbers up to $Re_w=10000$. Using the Monin-Obukhov length $L_{MO}$ we identify three different flow states, a buoyancy dominated regime ($L_{MO} \lesssim λ_θ$; with $λ_θ$ the thermal boundary layer thickness), a transitional regime ($0.5H \gtrsim L_{MO} \gtrsim λ_θ$; with $H$ the height of the domain), and a shear dominated regime ($L_{MO} \gtrsim 0.5H$). In the buoyancy dominated regime the flow dynamics are similar to that of turbulent thermal convection. The transitional regime is characterized by rolls that are increasingly elongated with increasing shear. The flow in the shear dominated regime consists of very large-scale meandering rolls, similar to the ones found in conventional Couette flow. As a consequence of these different flow regimes, for fixed $Ra$ and with increasing shear, the heat transfer first decreases, due to the breakup of the thermal rolls, and then increases at the beginning of the shear dominated regime. For $L_{MO} \gtrsim 0.5H$ the Nusselt number $Nu$ effectively scales as $Nu \sim Ra^α$, with $α\ll 1/3$ while we find $α\simeq 0.31$ in the buoyancy dominated regime. In the transitional regime the effective scaling exponent is $α> 1/3$, but the temperature and velocity profiles in this regime are not logarithmic yet, thus indicating transient dynamics and not the ultimate regime of thermal convection.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Solidification of liquid metal drops during impact

Hot liquid metal drops impacting onto a cold substrate solidify during their subsequent spreading. Here we experimentally study the influence of solidification on the outcome of an impact event. Liquid tin drops are impacted onto sapphire substrates of varying temperature. The impact is visualised both from the side and from below, which provides a unique view on the solidification process. During spreading an intriguing pattern of radial ligaments rapidly solidifies from the centre of the drop. This pattern determines the late-time morphology of the splat. A quantitative analysis of the drop spreading and ligament formation is supported by scaling arguments. Finally, a phase diagram for drop bouncing, deposition and splashing as a function of substrate temperature and impact velocity is provided.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Direct numerical simulations of Taylor--Couette turbulence: the effect of sand grain roughness

Progress in roughness research, mapping any given roughness geometry to its fluid dynamic behaviour, has been hampered by the lack of accurate and direct measurements of skin-friction drag, especially in open systems. The Taylor--Couette (TC) system has the benefit of being a closed system, but its potential for characterizing irregular, realistic, 3-D roughness has not been previously considered in depth. Here, we present direct numerical simulations (DNSs) of TC turbulence with sand grain roughness mounted on the inner cylinder. The model proposed by Scotti (\textit{Phys. Fluids}, vol. 18, 031701, 2006) has been improved to simulate a random rough surface of monodisperse sand grains, which is characterized by the equivalent sand grain height $k_s$. Taylor numbers range from $Ta = 1.0\times 10^7$(corresponding to $Re_τ= 82$) to $Ta = 1.0\times 10^9$($Re_τ= 635$). We focus on the influence of the roughness height $k_s^+$ in the transitionally rough regime, through simulations of TC with rough surfaces, ranging from $k_s^+=5$ up to $k_s^+ = 92$, where the superscript `$+$&#39; indicates non-dimensionalization in viscous units. We find that the downwards shift of the logarithmic layer, due to transitionally rough sand grains exhibits remarkably similar behavior to that of the Nikuradse (\textit{VDI-Forschungsheft} 361, 1933) data of sand grain roughness in pipe flow, regardless of the Taylor number dependent constants of the logarithmic layer.

0 citations0 reviewsNo code linkedOpen access