Source author record

Detlef Lohse

Detlef Lohse appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.flu-dyn cond-mat.soft nlin.CD physics.app-ph physics.chem-ph physics.comp-ph chao-dyn cond-mat.mtrl-sci physics.ao-ph physics.geo-ph cond-mat.stat-mech astro-ph.GA cond-mat.dis-nn physics.hist-ph physics.ins-det physics.med-ph physics.optics physics.soc-ph

Catalog footprint

What is connected

210works

18topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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 "porous" 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.

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.

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.

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.

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.

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.

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.

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 "co-" and "counter-rotating" due to their location or as "broad" and "narrow" 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.

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.

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.

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.

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

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.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

preprint2020arXiv

News & Views: The distance rule in times of Corona

The paper discussed the recent findings of Lydia Bourouiba and her group (MIT) on disease spreading, in the context of the Corona pandemic.

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.

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.

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

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.

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.

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.'s work.

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.

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

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.

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

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.

preprint2020arXiv

Towards DNS of the Ultimate Regime of Rayleigh--Bénard Convection

In this contribution we have briefly introduced the problem of turbulent thermal convection with a particular look at its transition to the ultimate regime and the resolution requirements needed for the direct numerical simulation of this flow.

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.

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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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 `$+$' 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.

preprint2016arXiv

Deformation and orientation statistics of neutrally buoyant sub-Kolmogorov ellipsoidal droplets in turbulent Taylor-Couette flow

The influence of the underlying flow topology on the shape and size of sub-Kolmogorov droplets dispersed in a turbulent flow is of considerable interest in many industrial and scientific applications. In this work we study the deformation and orientation statistics of sub-Kolmogorov droplets dispersed into a turbulent Taylor-Couette flow. Along with Direct Numerical Simulations (DNS) of the carrier phase and Lagrangian tracking of the dispersed droplets, we solve a phenomenological equation proposed by Maffettone and Minale (\emph{J. Fluid Mech.} 78, 227-241 (1998)) to track the shape evolution and orientation of approximately $10^5$ ellipsoidal droplets. By varying the capillary number $Ca$ and viscosity ratio $\hat μ$ of the droplets we find that the droplets deform more with increasing capillary number $Ca$ and this effect is more pronounced in the boundary layer regions. This indicates that along with a capillary number effect there is also a strong correlation between spatial position and degree of deformation of the droplet. Regardless of the capillary number $Ca$, the major-axis of the ellipsoids tends to align with the stream-wise direction and the extensional strain rate eigen direction in the boundary layer region while the distribution is highly isotropic in the bulk. When the viscosity ratio between the droplet and the carrier fluid is increased we find that there is no preferential stretched axis which is due to the increased influence of rotation over stretching and relaxation. Droplets in high viscosity ratio systems are thus less deformed and oblate (disk-like) as compared to highly deformed prolate (cigar-like) droplets in low viscosity ratio systems.

preprint2016arXiv

Direct numerical simulation of Taylor-Couette flow with grooved walls: torque scaling and flow structure

We present direct numerical simulations of Taylor-Couette flow with grooved walls at a fixed radius ratio $η=r_i/r_o=0.714$ with inner cylinder Reynolds number up to $Re_i=3.76\times10^4$, corresponding to Taylor number up to $Ta=2.15\times10^9$. The grooves are axisymmetric V-shaped obstacles attached to the wall with a tip angle of $90^\circ$. Results are compared to the smooth wall case in order to investigate the effects of grooves on Taylor-Couette flow. We focus on the effective scaling laws for the torque, flow structures, and boundary layers. It is found that, when the groove height is smaller than the boundary layer thickness, the torque is the same as that of the smooth wall cases. With increasing $Ta$, the boundary layer thickness becomes smaller than the groove height. Plumes are ejected from the tips of the grooves and secondary circulations between the latter are formed. This is associated to a sharp increase of the torque and thus the effective scaling law for the torque vs. $Ta$ becomes much steeper. Further increasing $Ta$ does not result in an additional slope increase. Instead, the effective scaling law saturates to the "ultimate" regime effective exponents seen for smooth walls. It is found that even though after saturation the slope is the same as for the smooth wall case, the absolute value of torque is increased, and the more the larger size of the grooves.

preprint2016arXiv

Disentangling the origins of torque enhancement through wall roughness in Taylor-Couette turbulence

Direct numerical simulations (DNSs) are performed to analyze the global transport properties of turbulent Taylor-Couette flow with inner rough wall up to Taylor number $Ta=10^{10}$. The dimensionless torque $Nu_ω$ shows an effective scaling of $Nu_ω\propto Ta^{0.42\pm0.01}$, which is steeper than the ultimate regime effective scaling $Nu_ω\propto Ta^{0.38}$ seen for smooth inner and outer walls. It is found that at the inner rough wall, the dominant contribution to the torque comes from the pressure forces on the radial faces of the rough elements; while viscous shear stresses on the rough surfaces contribute little to $Nu_ω$. Thus, the log layer close to the rough wall depends on the roughness length scale, rather than on the viscous length scale. We then separate the torque contributed from the smooth inner wall and the rough outer wall. It is found that the smooth wall torque scaling follows $Nu_s \propto Ta_s^{0.38\pm0.01}$, in excellent agreement with the case where both walls are smooth. In contrast, the rough wall torque scaling follows $Nu_r \propto Ta_r^{0.47\pm0.03}$, very close to the pure ultimate regime scaling $Nu_ω\propto Ta^{1/2}$. The energy dissipation rate at the wall of inner rough cylinder decreases significantly as a consequence of the wall shear stress reduction caused by the flow separation at the rough elements. On the other hand, the latter shed vortices in the bulk that are transported towards the outer cylinder and dissipated. Compared to the purely smooth case, the inner wall roughness renders the system more bulk dominated and thus increases the effective scaling exponent.

preprint2016arXiv

Dynamic Leidenfrost effect: relevant time- and length-scales

When a liquid droplet impacts a hot solid surface, enough vapor may be generated under it as to prevent its contact with the solid. The minimum solid temperature for this so-called Leidenfrost effect to occur is termed the Leidenfrost temperature, or the dynamic Leidenfrost temperature when the droplet velocity is non-negligible. We observe the wetting/drying and the levitation dynamics of the droplet impacting on an (isothermal) smooth sapphire surface using high speed total internal reflection imaging, which enables us to observe the droplet base up to about 100 nm above the substrate surface. By this method we are able to reveal the processes responsible for the transitional regime between the fully wetting and the fully levitated droplet as the solid temperature increases, thus shedding light on the characteristic time- and length-scales setting the dynamic Leidenfrost temperature for droplet impact on an isothermal substrate.

preprint2016arXiv

Energy spectra in turbulent bubbly flows

We conduct experiments in a turbulent bubbly flow to study the nature of the transition between the classical $-$5/3 energy spectrum scaling for a single-phase turbulent flow and the $-$3 scaling for a swarm of bubbles rising in a quiescent liquid and of bubble-dominated turbulence. The bubblance parameter, which measures the ratio of the bubble-induced kinetic energy to the kinetic energy induced by the turbulent liquid fluctuations before bubble injection, is often used to characterise the bubbly flow. We vary the bubblance parameter from $b = \infty$ (pseudo-turbulence) to $b = 0$ (single-phase flow) over 2-3 orders of magnitude ($0.01 - 5$) to study its effect on the turbulent energy spectrum and liquid velocity fluctuations. The probability density functions (PDFs) of the liquid velocity fluctuations show deviations from the Gaussian profile for $b > 0$, i.e. when bubbles are present in the system. The PDFs are asymmetric with higher probability in the positive tails. The energy spectra are found to follow the $-$3 scaling at length scales smaller than the size of the bubbles for bubbly flows. This $-$3 spectrum scaling holds not only in the well-established case of pseudo-turbulence, but surprisingly in all cases where bubbles are present in the system ($b > 0$). Therefore, it is a generic feature of turbulent bubbly flows, and the bubblance parameter is probably not a suitable parameter to characterise the energy spectrum in bubbly turbulent flows. The physical reason is that the energy input by the bubbles passes over only to higher wave numbers, and the energy production due to the bubbles can be directly balanced by the viscous dissipation in the bubble wakes as suggested by Lance $\&$ Bataille (1991). In addition, we provide an alternative explanation by balancing the energy production of the bubbles with viscous dissipation in the Fourier space.

preprint2016arXiv

Exploring the phase space of multiple states in highly turbulent Taylor-Couette flow

We investigate the existence of multiple turbulent states in highly turbulent Taylor-Couette flow in the range of $\mathrm{Ta}=10^{11}$ to $9\cdot10^{12}$, by measuring the global torques and the local velocities while probing the phase space spanned by the rotation rates of the inner and outer cylinder. The multiple states are found to be very robust and are expected to persist beyond $\mathrm{Ta}=10^{13}$. The rotation ratio is the parameter that most strongly controls the transitions between the flow states; the transitional values only weakly depend on the Taylor number. However, complex paths in the phase space are necessary to unlock the full region of multiple states. Lastly, by mapping the flow structures for various rotation ratios in a Taylor-Couette setup with an equal radius ratio but a larger aspect ratio than before, multiple states were again observed. Here, they are characterized by even richer roll structure phenomena, including, for the first time observed in highly turbulent TC flow, an antisymmetrical roll state.

preprint2016arXiv

Exposing nanobubble-like objects to a degassed environment

The primary attribute of interest of surface nanobubbles is their unusual stability and a number of theories trying to explain this have been put forward. Interestingly, the dissolution of nanobubbles is a topic that did not receive a lot of attention yet. In this work we applied two different experimental procedures which should cause gaseous nanobubbles to completely dissolve. In our experiments we nucleated nanobubble-like objects by putting a drop of water on HOPG using a plastic syringe and disposable needle. In method A, the nanobubble-like objects were exposed to a flow of degassed water (1.17 mg/l) for 96 hours. In method B, the ambient pressure was lowered in order to degas the liquid and the nanobubble-like objects. Interestingly, the nanobubble-like objects remained stable after exposure to both methods. After thorough investigation of the procedures and materials used during our experiments, we found that the nanobubble-like object were induced by the use of disposable needles in which PDMS contaminated the water. It is very important for the nanobubble community to be aware of the fact that, although features look and behave like nanobubbles, in some cases they might in fact be or induced by contamination. The presence of contamination could also resolve some inconsistencies found in the nanobubble literature.

preprint2016arXiv

Laser impact on a drop

The impact of a laser pulse on a highly absorbing liquid drop can lead to a violent response: the drop is accelerated, strongly deforms, and eventually fragments. Shock waves, the ejection of matter, and even plasma formation can accompany this process ...

preprint2016arXiv

Latent heat induced rotation limited aggregation in 2D ice nanocrystals

The basic science responsible for the fascinating shapes of ice crystals and snowflakes is still not understood. Insufficient knowledge of the interaction potentials and the lack of relevant experimental access to the growth process are to blame for this failure. Here, we study the growth of fractal nanostructures in a two-dimensional (2D) system, intercalated between mica and graphene. Based on our Scanning Tunneling Spectroscopy (STS) data we provide compelling evidence that these fractals are 2D ice. They grow while they are in material contact with the atmosphere at 20 $^{\circ}$C and without significant thermal contact to the ambient. The growth is studied in-situ, in real time and space at the nanoscale. We find that the growing 2D ice nanocrystals assume a fractal shape, which is conventionally attributed to Diffusion Limited Aggregation (DLA). However, DLA requires a low mass density mother phase, in contrast to the actual currently present high mass density mother phase. Latent heat effects and consequent transport of heat and molecules are found to be key ingredients for understanding the evolution of the snow (ice) flakes. We conclude that not the local availability of water molecules (DLA), but rather them having the locally required orientation is the key factor for incorporation into the 2D ice nanocrystal. In combination with the transport of latent heat we attribute the evolution of fractal 2D ice nanocrystals to local temperature dependent Rotation Limited Aggregation (RLA). The ice growth occurs under extreme supersaturation, i.e. the conditions closely resemble the natural ones for the growth of complex 2D snow (ice) flakes and we consider our findings crucial for solving the 'perennial' snow (ice) flake enigma.

preprint2016arXiv

Microbubbles and microparticles are not faithful tracers of turbulent acceleration

We report on the Lagrangian statistics of acceleration of small (sub-Kolmogorov) bubbles and tracer particles with Stokes number St << 1 in turbulent flow. At decreasing Reynolds number, the bubble accelerations show deviations from that of tracer particles, i.e. they deviate from the Heisenberg-Yaglom prediction and show a quicker decorrelation despite their small size and minute St. Using direct numerical simulations, we show that these effects arise due the drift of these particles through the turbulent flow. We theoretically predict this gravity-driven effect for developed isotropic turbulence, with the ratio of Stokes to Froude number or equivalently the particle drift-velocity governing the enhancement of acceleration variance and the reductions in correlation time and intermittency. Our predictions are in good agreement with experimental and numerical results. The present findings are relevant to a range of scenarios encompassing tiny bubbles and droplets that drift through the turbulent oceans and the atmosphere. They also question the common usage of micro-bubbles and micro-droplets as tracers in turbulence research.

preprint2016arXiv

Multiple states in highly turbulent Taylor-Couette flow

The ubiquity of turbulent flows in nature and technology makes it of utmost importance to fundamentally understand turbulence. Kolmogorov's 1941 paradigm suggests that for strongly turbulent flows with many degrees of freedom and its large fluctuations, there would only be \emph{one} turbulent state as the large fluctuations would explore the entire higher-dimensional phase space. Here we report the first conclusive evidence of multiple turbulent states for large Reynolds number $\text{Re}=\mathcal{O}(10^6)$ (Taylor number $\text{Ta}=\mathcal{O}(10^{12})$) Taylor-Couette flow in the regime of ultimate turbulence, by probing the phase space spanned by the rotation rates of the inner and outer cylinder. The manifestation of multiple turbulent states is exemplified by providing combined global torque and local velocity measurements. This result verifies the notion that bifurcations can occur in high-dimensional flows i.e. very large $\text{Re}$) and questions Kolmogorov's paradigm.

preprint2016arXiv

On the spreading of impacting drops

The energy budget and dissipation mechanisms during droplet impact on solid surfaces are studied numerically and theoretically. We find that for high impact velocities and negligible surface friction at the solid surface (i.e. free-slip), about one half of the initial kinetic energy is transformed into surface energy, independent of the impact parameters and the detailed energy loss mechanism(s). We argue that this seemingly universal rule is related to the deformation mode of the droplet and is reminiscent of pipe flow undergoing a sudden expansion, for which the head loss can be calculated by multiplying the kinetic energy of the incoming flow by a geometrical factor. For impacts on a no-slip surface also dissipation in the shear boundary layer at the solid surface is important. In this case the geometric head loss acts as a lower bound on the total dissipation (i.e. the spreading on a no-slip surface approaches that on a free-slip surface when the droplet viscosity is send to zero). This new view on the impact problem allows for simple analytical estimates of the maximum spreading diameter of impacting drops as a function of the impact parameters and the properties of the solid surface. It bridges the gap between previous momentum balance approaches and energy balance approaches, which hitherto did not give consistent predictions in the low viscosity limit. Good agreement is found between our models and experiments, both for impacts on "slippery" or lubricated surfaces (e.g. Leidenfrost droplet impacts and head-on droplet-droplet collisions) and for impacts on no-slip surfaces.

preprint2016arXiv

Particle tracking around surface nanobubbles

The exceptionally long lifetime of surface nanobubbles remains one of the biggest questions in the field. One of the proposed mechanisms for the stability is the \emph{dynamic equilibrium} model, which describes a constant flux of gas in and out of the bubble. Here, we describe results from particle tracking experiments to measure this flow. The results are analysed by measuring the Voronoï cell size distribution, the diffusion, and speed of the tracer particles. We show that there is no detectable difference in the movement of particles above nanobubble-laden surfaces, as compared to nanobubble-free surfaces.

preprint2016arXiv

Scaling laws and flow structures of double diffusive convection in the finger regime

Direct numerical simulations are conducted for double diffusive convection (DDC) bounded by two parallel plates, with fluid properties similar to the values of seawater. The DDC flow is driven by an unstable salinity difference and stabilized at the same time by a temperature difference. For these conditions the flow can be in the finger regime. We develop scaling laws for three key response parameters of the system: The non-dimensional salinity flux $Nu_S$ mainly depends on the salinity Rayleigh number $Ra_S$, which measures the strength of the salinity difference, and exhibits a very weak dependence on the density ratio $Λ$, which is the ratio of the buoyancy forces induced by two scalar differences. The non-dimensional flow velocity $Re$ and the non-dimensional heat flux $Nu_T$ are dependent on both $Ra_S$ and $Λ$. However, the rescaled Reynolds number $ReΛ^{α^{\rm eff}_u}$ and the rescaled convective heat flux $(Nu_T-1)Λ^{α^{\rm eff}_T}$ depend only on $Ra_S$. The two exponents are dependent on the fluid properties and are determined from the numerical results. Moreover, the behaviors of $Nu_S$ and $ReΛ^{α^{\rm eff}_u}$ agree with the predictions of the Grossmann-Lohse theory which was originally developed for the Rayleigh-Bénard flow. The non-dimensional salt-finger width and the thickness of the velocity boundary layers, after being rescaled by $Λ^{α^{\rm eff}_u/2}$, collapse and obey a similar power-law scaling relation with $Ra_S$. When $Ra_S$ is large enough, salt fingers do not extend from one plate to the other and horizontal zonal flows emerge in the bulk region. We then show that the current scaling strategy can be successfully applied to the experimental results of a heat-copper-ion system~(Hage and Tilgner, Phys. Fluids, 22, 076603, 2010).

preprint2016arXiv

Stick-jump mode in surface droplet dissolution

The analogy between evaporating surface droplets in air to dissolving long-chain alcohol droplets in water is worked out. We show that next to the three known modi for surface droplet evaporation or dissolution (constant contact angle mode, constant contact radius mode, and stick-slide mode), a fourth mode exists for small droplets on supposedly smooth substrates, namely the stick-jump mode: intermittent contact line pinning causes the droplet to switch between sticking and jumping during the dissolution. We present experimental data and compare them to theory to predict the dissolution time in this stick-jump mode. We also explain why these jumps were easily observed for microscale droplets but not for larger droplets.

preprint2016arXiv

The boiling Twente Taylor-Couette (BTTC) facility: Temperature controlled turbulent flow between independently rotating, coaxial cylinders

A new Taylor-Couette system has been designed and constructed with precise temperature control. Two concentric independently rotating cylinders are able to rotate at maximum rates of $f_i = \pm20$ Hz for the inner cylinder and $f_o = \pm10$ Hz for the outer cylinder. The inner cylinder has an outside radius of $r_i = 75$ mm, and the outer cylinder has an inside radius of $r_o = 105$ mm, resulting in a gap of $d=30$ mm. The height of the gap $L =549$ mm, giving a volume of $V=9.3$l. The geometric parameters are $η= r_i/r_o = 0.714$ and $Γ= L/d = 18.3$. With water as working fluid at room temperature the Reynolds numbers that can be achieved are $\text{Re}_i = ω_i r_i (r_o-r_i)/ν= 2.8 \times 10^5$ and $\text{Re}_o = ω_o r_o (r_o-r_i)/ν= 2 \times 10^5$, or a combined Reynolds number of up to $\text{Re} = (ω_i r_i -ω_o r_o)(r_o-r_i)/ν= 4.8 \times 10^5$. If the working fluid is changed to the fluorinated liquid FC-3284 with kinematic viscosity 0.42 cSt, the combined Reynolds number can reach $\text{Re} = 1.1 \times 10^6$. The apparatus features precise temperature control of the outer and inner cylinder separately, and is fully optically accessible from the side and top. The new facility offers the possibility to accurately study the process of boiling inside a turbulent flow and its effect on the flow.

preprint2016arXiv

Tribonucleation of bubbles

We report on the nucleation of bubbles on solids that are gently rubbed against each other in a liquid. The phenomenon is found to depend strongly on the material and roughness of the solid surfaces. For a given surface, temperature, and gas content, a trail of growing bubbles is observed if the rubbing force and velocity exceed a certain threshold. Direct observation through a transparent solid shows that each bubble in the trail results from the early coalescence of several microscopic bubbles, themselves detaching from microscopic gas pockets forming between the solids. From a detailed study of the wear tracks, with atomic force and scanning electron microscopy imaging, we conclude that these microscopic gas pockets originate from a local fracturing of the surface asperities, possibly enhanced by chemical reactions at the freshly created surfaces. Our findings will be useful either for preventing undesired bubble formation or, on the contrary, for "writing with bubbles," i.e., creating controlled patterns of microscopic bubbles.

preprint2016arXiv

Vertically bounded double diffusive convection in the fingering regime: comparing no-slip vs free-slip boundary conditions

Vertically bounded fingering double diffusive convection (DDC) is numerically investigated, focusing on the influences of different velocity boundary conditions, i.e. the no-slip condition which is inevitable in the lab-scale experimental research, and the free-slip condition which is an approximation for the interfaces in many natural environments, such as the oceans. For both boundary conditions the flow is dominated by fingers and the global responses follow the same scaling laws, with enhanced prefactors for the free-slip cases. Therefore, the laboratory experiments with the no-slip boundaries serve as a good model for the finger layers in the ocean. Moreover, in the free-slip case although the tangential shear stress is eliminated at the boundaries, the local dissipation rate in the near-wall region may exceed the value found in the no-slip cases, which is caused by the stronger vertical motions of fingers and sheet structures near the free-slip boundaries. This counter intuitive result might be relevant for properly estimating and modelling the mixing and entrainment phenomena at free-surfaces and interfaces.

preprint2016arXiv

Volume entrained in the wake of a disc intruding into an oil-water interface

An object moving through a plane interface into a fluid deforms the interface in such a way that fluid from one side of the interface is entrained into the other side, a phenomenon known as Darwin's drift. We investigate this phenomenon experimentally using a disc which is started exactly at the interface of two immiscible fluids, namely oil and water. First, we observe that due to the density difference between the two fluids the deformation of the interface is influenced by gravity, and show that there exits a time window of universal behavior. Secondly, we show by comparing with boundary integral simulations that, even though the deformation is universal, our results cannot be fully explained by potential flow solutions. We attribute this difference to the starting vortex, which is created in the wake of the disc. Universal behavior is preserved, however, because the size and strength of the vortex shows the same universality as the potential flow solution.

preprint2015arXiv

A multiple--resolution strategy for Direct Numerical Simulation of scalar turbulence

In this paper a numerical procedure to simulate low diffusivity scalar turbulence is presented. The method consists of using a grid for the advected scalar with a higher spatial resolutions than that of the momentum. The latter usually requires a less refined mesh and integrating both fields on a single grid tailored to the most demanding variable, produces an unnecessary computational overhead. A multiple resolution approach is used also in the time integration in order to maintain the stability of the scalars on the finer grid. The method is the more advantageous the less diffusive the scalar is with respect to momentum, therefore it is particularly well suited for large Prandtl or Schmidt number flows. However, even in the case of equal diffusivities the present procedure gives CPU time and memory occupation savings. The reason is that the absence of the pressure term in the scalar equation leads to much steeper gradients in the scalar field as compared to the velocity field.

preprint2015arXiv

Azimuthal velocity profiles in Rayleigh-stable Taylor-Couette flow and implied axial angular momentum transport

We present azimuthal velocity profiles measured in a Taylor-Couette apparatus, which has been used as a model of stellar and planetary accretion disks. The apparatus has a cylinder radius ratio of $η= 0.716$, an aspect-ratio of $Γ= 11.74$, and the plates closing the cylinders in the axial direction are attached to the outer cylinder. We investigate angular momentum transport and Ekman pumping in the Rayleigh-stable regime. The regime is linearly stable and is characterized by radially increasing specific angular momentum. We present several Rayleigh-stable profiles for shear Reynolds numbers $Re_S \sim O(10^5) \,$, both for $Ω_i > Ω_o > 0$ (quasi-Keplerian regime) and $Ω_o > Ω_i > 0$ (sub-rotating regime) where $Ω_{i,o}$ is the inner/outer cylinder rotation rate. None of the velocity profiles matches the non-vortical laminar Taylor-Couette profile. The deviation from that profile increased as solid-body rotation is approached at fixed $Re_S$. Flow super-rotation, an angular velocity greater than that of both cylinders, is observed in the sub-rotating regime. The velocity profiles give lower bounds for the torques required to rotate the inner cylinder that were larger than the torques for the case of laminar Taylor-Couette flow. The quasi-Keplerian profiles are composed of a well mixed inner region, having approximately constant angular momentum, connected to an outer region in solid-body rotation with the outer cylinder and attached axial boundaries. These regions suggest that the angular momentum is transported axially to the axial boundaries. Therefore, Taylor-Couette flow with closing plates attached to the outer cylinder is an imperfect model for accretion disk flows, especially with regard to their stability.

preprint2015arXiv

Comparison between two and three dimensional Rayleigh-Bénard convection

Two dimensional (2D) and three dimensional (3D) Rayleigh-Bénard convection is compared using results from direct numerical simulations and prior experiments. The explored phase diagrams for both cases are reviewed. The differences and similarities between 2D and 3D are studied using Nu(Ra) for Pr = 4.38 and Pr = 0.7 and Nu(Pr) for Ra up to $10^8$. In the Nu(Ra) scaling at higher Pr, 2D and 3D are very similar; differing only by a constant factor up to $\text{Ra} = 10^{10}$. In contrast, the difference is large at lower Pr, due to the strong roll state dependence of Nu in 2D. The behaviour of Nu(Pr) is similar in 2D and 3D at large Pr. However, it differs significantly around Pr = 1. The Reynolds number values are consistently higher in 2D and additionally converge at large Pr. Finally, the thermal boundary layer profiles are compared in 2D and 3D.

preprint2015arXiv

Drop deformation by laser-pulse impact

A free-falling absorbing liquid drop hit by a nanosecond laser-pulse experiences a strong recoil-pressure kick. As a consequence, the drop propels forward and deforms into a thin sheet which eventually fragments. We study how the drop deformation depends on the pulse shape and drop properties. We first derive the velocity field inside the drop on the timescale of the pressure pulse, when the drop is still spherical. This yields the kinetic-energy partition inside the drop, which precisely measures the deformation rate with respect to the propulsion rate, before surface tension comes into play. On the timescale where surface tension is important the drop has evolved into a thin sheet. Its expansion dynamics is described with a slender-slope model, which uses the impulsive energy-partition as an initial condition. Completed with boundary integral simulations, this two-stage model explains the entire drop dynamics and its dependance on the pulse shape: for a given propulsion, a tightly focused pulse results in a thin curved sheet which maximizes the lateral expansion, while a uniform illumination yields a smaller expansion but a flat symmetric sheet, in good agreement with experimental observations.

preprint2015arXiv

Drop Shaping by Laser-Pulse Impact

We show how the deposition of laser energy induces propulsion and strong deformation of an absorbing liquid body. Combining high speed with stroboscopic imaging, we observe that a millimeter-sized dyed water drop hit by a millijoule nanosecond laser pulse propels forward at several meters per second and deforms until it eventually fragments. The drop motion results from the recoil momentum imparted at the drop surface by water vaporization. We measure the propulsion speed and the time-deformation law of the drop, complemented by boundary-integral simulations. The drop propulsion and shaping are explained in terms of the laser-pulse energy, the drop size, and the liquid properties. These findings are, for instance, crucial for the generation of extreme ultraviolet light in nanolithography machines.

preprint2015arXiv

From convection rolls to finger convection in double-diffusive turbulence

Double diffusive convection (DDC), which is the buoyancy driven flow with fluid density depending on two scalar components, is ubiquitous in many natural and engineering enviroments. Of great interests are scalers transfer rate and flow structures. Here we systematically investigate DDC flow between two horizontal plates, driven by an unstable salinity gradient and stabilized by a temperature gradient. Counterintuitively, when increasing the stabilizing temperature gradient, the salinity flux first increases, even though the velocity monotonically decreases, before it finally breaks down to the purely diffusive value. The enhanced salinity transport is traced back to a transition in the overall flow pattern, namely from large scale convection rolls to well-organised vertically-oriented salt fingers. We also show and explain that the unifying theory of thermal convection originally developed by Grossmann and Lohse for Rayleigh-Bénard convection can be directly applied to DDC flow for a wide range of control parameters (Lewis number and density ratio), including those which cover the common values relevant for ocean flows.

preprint2015arXiv

Logarithmic mean temperature profiles and their connection to plume emissions in turbulent Rayleigh-Bénard convection

Two-dimensional simulations of Rayleigh-Bénard convection at $Ra = 5\times10^{10}$ show that vertical logarithmic mean temperature profiles can be observed in regions of the boundary layer where thermal plumes are emitted. The profile is logarithmic only in these regions and not in the rest of the boundary layer where it is sheared by the large scale wind and impacted by plumes. In addition, the logarithmic behavior is not visible in the horizontal average. The findings reveal that the temperature profiles are strongly connected to thermal plume emission and support a perception that parts of the boundary layer can be turbulent, while others are not. The transition to the ultimate regime, in which the boundary layers are considered to be fully turbulent, can therefore be understood as a gradual increases in fraction of the plume-emitting ('turbulent') regions of the boundary layer.

preprint2015arXiv

Salinity transfer in bounded double diffusive convection

The double diffusive convection between two parallel plates is numerically studied for a series of parameters. The flow is driven by the salinity difference and stabilized by the thermal field. Our simulations are directly compared to experiments by Hage and Tilgner (\emph{Phys. Fluids} 22, 076603 (2010)) for several sets of parameters and reasonable agreement is found. This in particular holds for the salinity flux and its dependence on the salinity Rayleigh number. Salt fingers are present in all simulations and extend through the entire height. The thermal Rayleigh number seems to have minor influence on salinity flux but affects the Reynolds number and the morphology of the flow. Next to the numerical calculation, we apply the Grossmann-Lohse theory for Rayleigh-Bénard flow to the current problem without introducing any new coefficients. The theory successfully predicts the salinity flux both with respect to the scaling and even with respect to the absolute value for the numerical and experimental results.

preprint2015arXiv

Taylor-Couette turbulence at radius ratio $η=0.5$: scaling, flow structures and plumes

Using high-resolution particle image velocimetry we measure velocity profiles, the wind Reynolds number and characteristics of turbulent plumes in Taylor-Couette flow for a radius ratio of 0.5 and Taylor number of up to $6.2\cdot10^9$. The extracted angular velocity profiles follow a log-law more closely than the azimuthal velocity profiles due to the strong curvature of this $η=0.5$ setup. The scaling of the wind Reynolds number with the Taylor number agrees with the theoretically predicted 3/7-scaling for the classical turbulent regime, which is much more pronounced than for the well-explored $η=0.71$ case, for which the ultimate regime sets in at much lower Ta. By measuring at varying axial positions, roll structures are found for counter-rotation while no clear coherent structures are seen for pure inner cylinder rotation. In addition, turbulent plumes coming from the inner and outer cylinder are investigated. For pure inner cylinder rotation, the plumes in the radial velocity move away from the inner cylinder, while the plumes in the azimuthal velocity mainly move away from the outer cylinder. For counter-rotation, the mean radial flow in the roll structures strongly affects the direction and intensity of the turbulent plumes. Furthermore, it is experimentally confirmed that in regions where plumes are emitted, boundary layer profiles with a logarithmic signature are created.

preprint2015arXiv

The effect of velocity boundary conditions on the heat transfer and flow topology in two-dimensional Rayleigh-Bénard convection

The effect of various velocity boundary condition is studied in two-dimensional Rayleigh-Bénard convection. Combinations of no-slip, stress-free and periodic boundary conditions are used on both the sidewalls and the horizontal plates. For the studied Rayleigh numbers Ra between $10^8$ and $10^{11}$ the heat transport is lower for $Γ= 0.33$ than for $Γ= 1$ in case of no-slip sidewalls. This is surprisingly opposite for stress-free sidewalls, where the heat transport increases for lower aspect-ratio. In wider cells the aspect-ratio dependence is observed to disappear for $\text{Ra} \ge 10^{10}$. Two distinct flow types with very different dynamics can be seen, mostly dependent on the plate velocity boundary condition, namely roll-like flow and horizontal zonal flow, which have a substantial effect on the dynamics and heat transport in the system. The predominantly horizontal zonal flow suppresses heat flux and is observed for stress-free and asymmetric plates. Low aspect-ratio periodic sidewall simulations with a no-slip boundary condition on the plates also exhibit zonal flow. In all the other cases, the flow is roll-like. In two-dimensional Rayleigh-Bénard convection, the velocity boundary conditions thus have large implications on both roll-like and zonal flow that have to be taken into consideration before the boundary conditions are imposed.

preprint2015arXiv

Turbulent Rayleigh-Bénard convection described by projected dynamics in phase space

Rayleigh-Bénard convection, i.e. the flow of a fluid between two parallel plates that is driven by a temperature gradient, is an idealised setup to study thermal convection. Of special interest are the statistics of the turbulent temperature field, which we are investigating and comparing for three different geometries, namely convection with periodic horizontal boundary conditions in three and two dimensions as well as convection in a cylindrical vessel, in order to work out similarities and differences. To this end, we derive an exact evolution equation for the temperature probability density function (PDF). Unclosed terms are expressed as conditional averages of velocities and heat diffusion, which are estimated from direct numerical simulations. This framework lets us identify the average behaviour of a fluid particle by revealing the mean evolution of fluid of different temperatures in different parts of the convection cell. We connect the statistics to the dynamics of Rayleigh-Bénard convection, giving deeper insights into the temperature statistics and transport mechanisms. We find that the average behaviour is described by closed cycles in phase space that reconstruct the typical Rayleigh-Bénard cycle of fluid heating up at the bottom, rising up to the top plate, cooling down and falling down again. The detailed behaviour shows subtle differences between the three cases.

preprint2015arXiv

Universal mechanism for air entrainment during liquid impact

When a mm-sized liquid drop approaches a deep liquid pool, both the interface of the drop and the pool deform before the drop touches the pool. The build up of air pressure prior to coalescence is responsible for this deformation. Due to this deformation, air can be entrained at the bottom of the drop during the impact. We quantify the amount of entrained air numerically, using the Boundary Integral Method (BIM) for potential flow for the drop and the pool, coupled to viscous lubrication theory for the air film that has to be squeezed out during impact. We compare our results to various experimental data and find excellent agreement for the amount of air that is entrapped during impact onto a pool. Next, the impact of a rigid sphere onto a pool is numerically investigated and the air that is entrapped in this case also matches with available experimental data. In both cases of drop and sphere impact onto a pool the numerical air bubble volume V_b is found to be in agreement with the theoretical scaling V_b/V_{drop/sphere} ~ St^{-4/3}, where St is the Stokes number. This is the same scaling that has been found for drop impact onto a solid surface in previous research. This implies a universal mechanism for air entrainment for these different impact scenarios, which has been suggested in recent experimental work, but is now further elucidated with numerical results.

preprint2015arXiv

Wake-driven dynamics of finite-sized buoyant spheres in turbulence

Particles suspended in turbulent flows are affected by the turbulence and at the same time act back on the flow. The resulting coupling can give rise to rich variability in their dynamics. Here we report experimental results from an investigation of finite-sized buoyant spheres in turbulence. We find that even a marginal reduction in the particle's density from that of the fluid can result in strong modification of its dynamics. In contrast to classical spatial filtering arguments and predictions of particle models, we find that the particle acceleration variance increases with size. We trace this reversed trend back to the growing contribution from wake-induced forces, unaccounted for in current particle models in turbulence. Our findings highlight the need for improved multi-physics based models that account for particle wake effects for a faithful representation of buoyant-sphere dynamics in turbulence.

preprint2014arXiv

From antinode clusters to node clusters: The concentration dependent transition of floaters on a standing Faraday wave

A hydrophilic floating sphere that is denser than water drifts to an amplitude maximum (antinode) of a surface standing wave. A few identical floaters therefore organize into antinode clusters. However, beyond a transitional value of the floater concentration $ϕ$, we observe that the same spheres spontaneously accumulate at the nodal lines, completely inverting the self-organized particle pattern on the wave. From a potential energy estimate we show (i) that at low $ϕ$ antinode clusters are energetically favorable over nodal ones and (ii) how this situation reverses at high $ϕ$, in agreement with the experiment.

preprint2014arXiv

Heat transport in boiling turbulent Rayleigh-Bénard convection

Boiling is an extremely effective way to promote heat transfer from a hot surface to a liquid due to several mechanisms many of which are not understood in quantitative detail. An important component of the overall process is that the buoyancy of the bubbles compounds with that of the liquid to give rise to a much enhanced natural convection. In this paper we focus specifically on this enhancement and present a numerical study of the resulting two-phase Rayleigh-Bénard convection process. We make no attempt to model other aspects of the boiling process such as bubble nucleation and detachment. We consider a cylindrical cell with a diameter equal to its height. The cell base and top are held at temperatures above and below the boiling point of the liquid, respectively. By keeping the temperature difference constant and changing the liquid pressure we study the effect of the liquid superheat in a Rayleigh number range that, in the absence of boiling, would be between $2\times10^6$ and $5\times10^9$. We find a considerable enhancement of the heat transfer and study its dependence on the bubble number, the degree of superheat of the hot cell bottom and the Rayleigh number. The increased buoyancy provided by the bubbles leads to more energetic hot plumes detaching from the hot cell bottom and, as a consequence, the strength of the circulation in the cell is significantly increased. Our results are in general agreement with recent experimental results of Zhong et al., Phys. Rev. Lett. {\bf 102}, 124501 (2009) for boiling Rayleigh-Bénard convection.

preprint2014arXiv

Quantfying rich patterns in agglomeration of floating beads

Macroscopic spherical particles spontaneously form rich patterns on a standing Faraday wave. These patterns are found to follow a very systematic trend depending on the floater concentration $ϕ$: The same floaters that accumulate at amplitude maxima (antinodes) of the wave at low $ϕ$, surprisingly move towards the nodal lines when $ϕ$ is beyond a certain value. In more detail, circular irregularly packed antinode clusters at low $ϕ$ give way to loosely packed filamentary structures at intermediate $ϕ$, and are then followed by densely packed grid-shaped node clusters at high $ϕ$. Here, we successfully characterize the morphology of these rich patterns using a metric analysis, i.e., the Minkowski functionals. We modify the Minkowski functionals such that we are able to measure the physical quantities of the clusters such as area, perimeter, and aspect ratio.

preprint2014arXiv

Sidewall effects in Rayleigh-Bénard convection

We investigate the influence of the temperature boundary conditions at the sidewall on the heat transport in Rayleigh-Bénard (RB) convection using direct numerical simulations. For relatively low Rayleigh numbers Ra the heat transport is higher when the sidewall is isothermal, kept at a temperature $T_c+Δ/2$ (where $Δ$ is the temperature difference between the horizontal plates and $T_c$ the temperature of the cold plate), than when the sidewall is adiabatic. The reason is that in the former case part of the heat current avoids the thermal resistance of the fluid layer by escaping through the sidewall that acts as a short-circuit. For higher Ra the bulk becomes more isothermal and this reduces the heat current through the sidewall. Therefore the heat flux in a cell with an isothermal sidewall converges to the value obtained with an adiabatic sidewall for high enough Ra ($\simeq 10^{10}$). However, when the sidewall temperature deviates from $T_c+Δ/2$ the heat transport at the bottom and top plates is different from the value obtained using an adiabatic sidewall. In this case the difference does not decrease with increasing Ra thus indicating that the ambient temperature of the experimental apparatus can influence the heat transfer. A similar behavior is observed when only a very small sidewall region close to the horizontal plates is kept isothermal, while the rest of the sidewall is adiabatic. The reason is that in the region closest to the horizontal plates the temperature difference between the fluid and the sidewall is highest. This suggests that one should be careful with the placement of thermal shields outside the fluid sample to minimize spurious heat currents.

preprint2013arXiv

Air flow in a collapsing cavity

We experimentally study the airflow in a collapsing cavity created by the impact of a circular disk on a water surface. We measure the air velocity in the collapsing neck in two ways: Directly, by means of employing particle image velocimetry of smoke injected into the cavity and indirectly, by determining the time rate of change of the volume of the cavity at pinch-off and deducing the air flow in the neck under the assumption that the air is incompressible. We compare our experiments to boundary integral simulations and show that close to the moment of pinch-off, compressibility of the air starts to play a crucial role in the behavior of the cavity. Finally, we measure how the air flow rate at pinch-off depends on the Froude number and explain the observed dependence using a theoretical model of the cavity collapse.

preprint2013arXiv

Drop impact of shear thickening liquids

The impact of drops of concentrated non-Brownian suspensions (cornstarch and polystyrene spheres) onto a solid surface is investigated experimentally. The spreading dynamics and maxi- mal deformation of the droplet of such shear thickening liquids are found to be markedly different from the impact of Newtonian drops. A particularly striking observation is that the maximal de- formation is independent of the drop velocity and that the deformation suddenly stops during the impact phase. Both observations are due to the shear-thickening rheology of the suspensions, as is theoretically explained from a balance between the kinetic energy and the viscously-dissipated en- ergy, from which we establish a scaling relation between drop maximal deformation and rheological parameters of concentrated suspensions.

preprint2013arXiv

Fast & Scalded: Capillary Leidenfrost Droplets in micro-Ratches

In this Fluid Dynamics Videos submitted to the 31st Gallery of Fluid Motion, (find a different version here http://youtu.be/CS0c05WQ_js) we illustrate the special dynamics of capillary self-propelled Leidenfrost droplets in micrometric Ratchets. In order to be able to propel water droplets of sizes of the order of 1 mm, micro-ratchets were produced by direct material removal using a picosecond pulsed laser source. Surface micro-patterning with picosecond laser pulses allows creating a well controlled topography on a variety of substrates, with a resolution typically in the micron range. More information can be found in references.

preprint2013arXiv

Growing bubbles in a slightly supersaturated liquid solution

We have designed and constructed an experimental system to study gas bubble growth in slightly supersatu- rated liquids. This is achieved by working with carbon dioxide dissolved in water, pressurized at a maximum of 1 MPa and applying a small pressure drop from saturation conditions. Bubbles grow from hydrophobic cavities etched on silicon wafers, which allows us to control their number and position. Hence, the experiment can be used to investigate the interaction among bubbles growing in close proximity when the main mass transfer mechanism is diffusion and there is a limited availability of the dissolved species.

preprint2013arXiv

Heat transport and flow structure in rotating Rayleigh-Bénard convection

Here we summarize the results from our direct numerical simulations (DNS) and experimental measurements on rotating Rayleigh-Bénard (RB) convection. Our experiments and simulations are performed in cylindrical samples with an aspect ratio Γvarying from 1/2 to 2. Here Γ=D/L, where D and L are the diameter and height of the sample, respectively. When the rotation rate is increased, while a fixed temperature difference between the hot bottom and cold top plate is maintained, a sharp increase in the heat transfer is observed before the heat transfer drops drastically at stronger rotation rates. Here we focus on the question of how the heat transfer enhancement with respect to the non-rotating case depends on the Rayleigh number Ra, the Prandtl number Pr, and the rotation rate, indicated by the Rossby number Ro. Special attention will be given to the influence of the aspect ratio on the rotation rate that is required to get heat transport enhancement. In addition, we will discuss the relation between the heat transfer and the large scale flow structures that are formed in the different regimes of rotating RB convection and how the different regimes can be identified in experiments and simulations.

preprint2013arXiv

Interacting bubble clouds and their sonochemical production

Acoustically driven air pockets trapped in artificial crevices on a sur- face can emit bubbles which organize in (interacting) bubble clusters. With increasing driving power Fernandez Rivas et al. [Angew. Chem. Int. Ed., 2010] observed three different behaviors: clusters close to the very pits out of which they had been created, clusters pointing toward each other, and merging clusters. The latter behavior is highly undesired for technological purposes as it is associated with a reduction of the radical production and an enhancement of the erosion of the reactor walls. The dependence on the control parameters such as the distance of the pits and the conditions for cluster-merging are examined. The underlying mechanism, governed by the secondary Bjerknes forces, turns out to be strongly influenced by the nonlinearity of the bubble oscillations and not directly by the number of nucleated bubbles. The Bjerknes forces are found to dampen the bubble oscillations, thus reducing the radical production. Therefore, the increased number of bubbles at high power could be the key to understand the experimental observation that, after a certain power threshold, any further increase of the driving does not improve the sonochemical efficiency.

preprint2013arXiv

Leidenfrost temperature increase for impacting droplets on carbon-nanofiber surfaces

Droplets impacting on a superheated surface can either exhibit a contact boiling regime, in which they make direct contact with the surface and boil violently, or a film boiling regime, in which they remain separated from the surface by their own vapor. The transition from the contact to the film boiling regime depends not only on the temperature of the surface and kinetic energy of the droplet, but also on the size of the structures fabricated on the surface. Here we experimentally show that surfaces covered with carbon-nanofibers delay the transition to film boiling to much higher temperature compared to smooth surfaces. We present physical arguments showing that, because of the small scale of the carbon fibers, they are cooled by the vapor flow just before the liquid impact, thus permitting contact boiling up to much higher temperatures than on smooth surfaces. We also show that, as long as the impact is in the film boiling regime, the spreading factor of impacting droplets follows the same $\We^{3/10}$ scaling (with $\We$ the Weber number) found for smooth surfaces, which is caused by the vapor flow underneath the droplet.

preprint2013arXiv

Logarithmic boundary layers in highly turbulent Taylor-Couette flow

We provide direct measurements of the boundary layer properties in highly turbulent Taylor-Couette flow up to $\text{Ta}=6.2 \times 10^{12}$ using high-resolution particle image velocimetry (PIV). We find that the mean azimuthal velocity profile at the inner and outer cylinder can be fitted by the von Kármán log law $u^+ = \frac 1κ\ln y^+ +B$. The von Kármán constant $κ$ is found to depend on the driving strength $\text{Ta}$ and for large $\text{Ta}$ asymptotically approaches $κ\approx 0.40$. The variance profiles of the local azimuthal velocity have a universal peak around $y^+ \approx 12$ and collapse when rescaled with the driving velocity (and not with the friction velocity), displaying a log-dependence of $y^+$ as also found for channel and pipe flows [1,2].

preprint2013arXiv

Stability of viscous long liquid filaments

We study the collapse of an axisymmetric liquid filament both analytically and by means of a numerical model. The liquid filament, also known as ligament, may either collapse stably into a single droplet or break up into multiple droplets. The dynamics of the filament are governed by the viscosity and the aspect ratio, and the initial perturbations of its surface. We find that the instability of long viscous filaments can be completely explained by the Rayleigh-Plateau instability, whereas a low viscous filament can also break up due to end pinching. We analytically derive the transition between stable collapse and breakup in the Ohnesorge number versus aspect ratio phase space. Our result is confirmed by numerical simulations based on the slender jet approximation and explains recent experimental findings by Castrejon-Pita et al., PRL 108, 074506 (2012).

preprint2013arXiv

Statistics of turbulent fluctuations in counter-rotating Taylor-Couette flows

The statistics of velocity fluctuations of turbulent Taylor-Couette flow are examined. The rotation rate of the inner and outer cylinder are varied while keeping the Taylor number fixed to $1.49 \times 10^{12}$ ($\mathcal{O}(\text{Re})=10^6$). The azimuthal velocity component of the flow is measured using laser Doppler anemometry (LDA). For each experiment $5\times10^6$ datapoints are acquired and carefully analysed. Using extended self-similarity (ESS) \cite{ben93b} the longitudinal structure function exponents are extracted, and are found to weakly depend on the ratio of the rotation rates. For the case where only the inner cylinder rotates the results are in good agreement with results measured by Lewis and Swinney \cite{lew99} using hot-film anemometry. The power spectra shows clear -5/3 scaling for the intermediate angular velocity ratios $-ω_o/ω_i \in \{0.6, 0.8, 1.0\}$, roughly -5/3 scaling for $-ω_o/ω_i \in \{0.2, 0.3, 0.4, 2.0\}$, and no clear scaling law can be found for $-ω_0/ω_i = 0$ (inner cylinder rotation only); the local scaling exponent of the spectra has a strong frequency dependence. We relate these observations to the shape of the probability density function of the azimuthal velocity and the presence of a neutral line.

preprint2013arXiv

The importance of bubble deformability for strong drag reduction in bubbly turbulent Taylor-Couette flow

Bubbly turbulent Taylor-Couette (TC) flow is globally and locally studied at Reynolds numbers of Re = 5 x 10^5 to 2 x 10^6 with a stationary outer cylinder and a mean bubble diameter around 1 mm. We measure the drag reduction (DR) based on the global dimensional torque as a function of the global gas volume fraction a_global over the range 0% to 4%. We observe a moderate DR of up to 7% for Re = 5.1 x 10^5. Significantly stronger DR is achieved for Re = 1.0 x 10^6 and 2.0 x 10^6 with, remarkably, more than 40% of DR at Re = 2.0 x 10^6 and a_global = 4%. To shed light on the two apparently different regimes of moderate DR and strong DR, we investigate the local liquid flow velocity and the local bubble statistics, in particular the radial gas concentration profiles and the bubble size distribution, for the two different cases; Re = 5.1 x 10^5 in the moderate DR regime and Re = 1.0 x 10^6 in the strong DR regime, both at a_global = 3 +/- 0.5%. By defining and measuring a local bubble Weber number (We) in the TC gap close to the IC wall, we observe that the crossover from the moderate to the strong DR regime occurs roughly at the crossover of We ~ 1. In the strong DR regime at Re = 1.0 x 10^6 we find We > 1, reaching a value of 9 (+7, -2) when approaching the inner wall, indicating that the bubbles increasingly deform as they draw near the inner wall. In the moderate DR regime at Re = 5.1 x 10^5 we find We ~ 1, indicating more rigid bubbles, even though the mean bubble diameter is larger, namely 1.2 (+0.7, -0.1) mm, as compared to the Re = 1.0 x 10^6 case, where it is 0.9 (+0.6, -0.1) mm. We conclude that bubble deformability is a relevant mechanism behind the observed strong DR. These local results match and extend the conclusions from the global flow experiments as found by van den Berg et al. (2005) and from the numerical simulations by Lu, Fernandez & Tryggvason (2005).

preprint2013arXiv

The microfluidic Kelvin water dropper

The so-called "Kelvin water dropper" is a simple experiment demonstrating the spontaneous appearance of induced free charge in droplets emitted through a tube. As Lord Kelvin explained, water droplets spontaneously acquire a net charge during detachment from a faucet due to the presence of electrical fields in their surrounding created by any metallic object. In his experiment, two streams of droplets are allowed to drip from separated nozzles into separated buckets, which are at the same time interconnected through the dripping needles. In this paper we build a microfluidic water dropper and demonstrate that the droplets get charged and break-up due to electrohydrodynamic instabilities. A comparison with recent simulations shows the dependence of the acquired charge in the droplets on different parameters of the system. The phenomenon opens a door to cheap and accessible transformation of pneumatic pressure into electrical energy and to an enhanced control in microfluidic and biophysical manipulation of capsules, cells and droplets via self-induced charging of the elements.

preprint2013arXiv

The unifying theory of scaling in thermal convection: The updated prefactors

The unifying theory of scaling in thermal convection (Grossmann & Lohse (2000)) (henceforth the GL theory) suggests that there are no pure power laws for the Nusselt and Reynolds numbers as function of the Rayleigh and Prandtl numbers in the experimentally accessible parameter regime. In Grossmann & Lohse (2001) the dimensionless parameters of the theory were fitted to 155 experimental data points by Ahlers & Xu (2001) in the regime $3\times 10^7 \le Ra \le 3 \times 10^{9}$ and $4\le Pr \le 34$ and Grossmann & Lohse (2002) used the experimental data point from Qiu & Tong (2001) and the fact that Nu(Ra,Pr) is independent of the parameter a, which relates the dimensionless kinetic boundary thickness with the square root of the wind Reynolds number, to fix the Reynolds number dependence. Meanwhile the theory is on one hand well confirmed through various new experiments and numerical simulations. On the other hand these new data points provide the basis for an updated fit in a much larger parameter space. Here we pick four well established (and sufficiently distant) Nu(Ra,Pr) data points and show that the resulting Nu(Ra,Pr) function is in agreement with almost all established experimental and numerical data up to the ultimate regime of thermal convection, whose onset also follows from the theory. One extra Re(Ra,Pr) data point is used to fix Re(Ra,Pr). As Re can depend on the definition and the aspect ratio the transformation properties of the GL equations are discussed in order to show how the GL coefficients can easily be adapted to new Reynolds number data while keeping Nu(Ra,Pr) unchanged.

preprint2013arXiv

Velocity profiles in strongly turbulent Taylor-Couette flow

We derive the velocity profiles in strongly turbulent Taylor-Couette flow for the general case of independently rotating cylinders. The theory is based on the Navier-Stokes equations in the appropriate (cylinder) geometry. In particular, we derive the axial and the angular velocity profiles as functions of distance from the cylinder walls and find that both follow a logarithmic profile, with downwards-bending curvature corrections, which are more pronounced for the angular velocity profile as compared to the axial velocity profile, and which strongly increase with decreasing ratio $η$ between inner and outer cylinder radius. In contrast, the azimuthal velocity does not follow a log-law. We then compare the angular and azimuthal velocity profiles with the recently measured profiles in the ultimate state of (very) large Taylor numbers. Though the {\em qualitative} trends are the same -- down-bending for large wall distances and (properly shifted and non-dimensionalized) angular velocity profile $ω^+(r)$ being closer to a log-law than (properly shifted and non-dimensionalized) azimuthal velocity profile $u^+_φ(r)$ -- {\em quantitative} deviations are found for large wall distances. We attribute these differences to the Taylor rolls and the height dependence of the profiles, neither of which are considered in the theoretical approach.

preprint2012arXiv

Breakdown of the large-scale wind in Γ=1/2 rotating Rayleigh-Bénard flow

Experiments and simulations of rotating Rayleigh-Bénard convection in cylindrical samples have revealed an increase in heat transport with increasing rotation rate. This heat transport enhancement is intimately related to a transition in the turbulent flow structure from a regime dominated by a large-scale circulation (LSC), consisting of a single convection roll, at no or weak rotation to a regime dominated by vertically-aligned vortices at strong rotation. For a sample with an aspect ratio Γ= D/L = 1 (D is the sample diameter and L its height) the transition between the two regimes is indicated by a strong decrease in the LSC strength. In contrast, for Γ= 1/2 Weiss and Ahlers [J. Fluid Mech. {\bf{688}}, 461 (2011)] revealed the presence of a LSC-like sidewall temperature signature beyond the critical rotation rate. They suggested that this might be due to the formation of a two-vortex state, in which one vortex extends vertically from the bottom into the sample interior and brings up warm fluid, while another vortex brings down cold fluid from the top; this flow field would yield a sidewall temperature signature similar to that of the LSC. Here we show by direct numerical simulations for Γ= 1/2 and parameters that allow direct comparison with experiment that the spatial organization of the vertically-aligned vortical structures in the convection cell do indeed yield (for the time average) a sinusoidal variation of the temperature near the sidewall, as found in the experiment. This is also the essential and non-trivial difference with the Γ=1 sample, where the vertically-aligned vortices are distributed randomly.

preprint2012arXiv

Building micro-soccer-balls with evaporating colloidal fakir drops

Evaporation-driven particle self-assembly can be used to generate three-dimensional microstructures. We present a new method to create these colloidal microstructures, in which we can control the amount of particles and their packing fraction. To this end, we evaporate colloidal dispersion droplets on a special type of superhydrophobic micro-structured surface, on which the droplet re- mains in Cassie-Baxter state during the entire evaporative process. The remainders of the droplet consist of a massive spherical cluster of the microspheres, with diameters ranging from a few tens up to several hundreds of microns. We present scaling arguments to show how the final particle packing fraction of these balls depends on the dynamics of the droplet evaporation.

preprint2012arXiv

Capillary droplets on Leidenfrost micro-ratchets

Leidenfrost ratchets are structures with the ability of transporting liquid droplets when heated over the critical Leidenfrost temperature. Once this temperature is reached, the droplet levitates over the surface and moves in the direction marked by the slope of the ratchet at terminal velocities around 10 cm/s. Here we provide new experiments with micron-sized ratchets, which have been produced with picosecond pulse laser ablation. In the following work, we use a simple method to measure the thrust driving droplets of capillary size over the micro-ratchets. The mechanism responsible for the force acting on the drop on superheated ratchets has been recently under debate. We extend the recently proposed 'viscous mechanism' proposed by Dupeaux et al. [Europhys. Lett., 96, 58001 (2011)] to capillary droplets and find good agreement with our measurements.

preprint2012arXiv

Collapse of a non-axisymmetric, impact-created air cavity in water

The axisymmetric collapse of a cylindrical air cavity in water follows a universal power law with logarithmic corrections. Nonetheless, it has been suggested that the introduction of a small azimuthal disturbance induces a long term memory effect, reflecting in oscillations which are no longer universal but remember the initial condition. In this work, we create non-axisymmetric air cavities by driving a metal disc through an initially-quiescent water surface and observe their subsequent gravity-induced collapse. The cavities are characterized by azimuthal harmonic disturbances with a single mode number $m$ and amplitude $a_m$. For small initial distortion amplitude (1 or 2% of the mean disc radius), the cavity walls oscillate linearly during collapse, with nearly constant amplitude and increasing frequency. As the amplitude is increased, higher harmonics are triggered in the oscillations and we observe more complex pinch-off modes. For small amplitude disturbances we compare our experimental results with the model for the amplitude of the oscillations by Schmidt et al. (2009) and the model for the collapse of an axisymmetric impact-created cavity previously proposed by Bergmann et al. (2009b). By combining these two models we can reconstruct the three-dimensional shape of the cavity at any time before pinch-off.

preprint2012arXiv

Efficient Sonochemistry through Microbubbles Generated with Micromachined Surfaces

Sonochemical reactors are used in water treatment, the synthesis of fine chemicals, pharmaceutics and others. The low efficiency of sonoreactors have prevented its massive usage at industrial scales. Controlling the appearance of bubbles in place and time is the most limiting factor. A novel type of sonochemical reactor was designed making use of micro-fabrication techniques to control the nucleation sites of micro-bubbles. The efficiency was increased first by locating the nucleation sites in the most active region of a micro-chamber; additionally the desired chemical effect was significantly higher at the same powers than when not controlled. Silicon substrates were micromachined with "artificial nucleation sites" or pits, and placed at the bottom of the micro-chamber. The pits entrap gas which, upon ultrasonic excitation, sheds off a stream of microbubbles. The gas content of the pits is not depleted but is replenished by diffusion and the emission of microbubbles can continue for hours.

preprint2012arXiv

Flow states in two-dimensional Rayleigh-Bénard convection as a function of aspect-ratio and Rayleigh number

In this numerical study on two-dimensional Rayleigh-Bénard convection we consider $10^7 \leq Ra \leq 10^{12}$ in aspect ratio $0.23 \leq Γ\leq 13$ samples. We focus on several cases. First we consider small aspect ratio cells, where at high Ra number we find a sharp transition from a low Ra number branch towards a high Ra number branch, due to changes in the flow structure. Subsequently, we show that the influence of the aspect ratio on the heat transport decreases with increasing aspect ratio, although even at very large aspect ratio of $Γ\approx10$ variations up to 2.5% in the heat transport as a function of Γare observed. Finally, we observe long-lived transients up to at least $Ra=10^9$, as in certain aspect ratio cells we observe different flow states that are stable for thousands of turnover times.

preprint2012arXiv

Fluctuation theorems for an asymmetric rotor in a granular gas

We investigate the validity of fluctuation theorems for an asymmetric rotor experiment in a granular gas. A first state, with a Gaussian distribution of the angular velocity, is found to be well described by a first order Langevin equation. We show that fluctuation theorems are valid for the injected work and for the total entropy production. In a second state the angular velocity distribution is double-peaked due to a spontaneous symmetry breaking: A convection roll develops in the granular gas, which strongly couples to the rotor. Surprisingly, in this case similar symmetry relations hold, which lead to a good prediction for the height ratio of the two peaks.

preprint2012arXiv

Highly focused supersonic microjets

The paper describes the production of thin, focused microjets with velocities up to 850 m/s by the rapid vaporization of a small mass of liquid in an open liquid-filled capillary. The vaporization is caused by the absorption of a low-energy laser pulse. A likely explanation of the observed phenomenon is based on the impingement of the shock wave caused by the nearly-instantaneous vaporization on the free surface of the liquid. An experimental study of the dependence of the jet velocity on several parameters is conducted, and a semi-empirical relation for its prediction is developed. The coherence of the jets, their high velocity and good reproducibility and controllability are unique features of the system described. A possible application is to the development of needle-free drug injection systems which are of great importance for global health care.

preprint2012arXiv

Highly focused supersonic microjets: numerical simulations

By focusing a laser pulse inside a capillary partially filled with liquid, a vapour bubble is created which emits a pressure wave. This pressure wave travels through the liquid and creates a fast, focused axisymmetric microjet when it is reflected at the meniscus. We numerically investigate the formation of this microjet using axisymmetric boundary-integral simulations, where we model the pressure wave as a pressure pulse applied on the bubble. We find a good agreement between the simulations and experimental results in terms of the time evolution of the jet and on all parameters that can be compared directly. We present a simple analytical model that accurately predicts the velocity of the jet after the pressure pulse and its maximum velocity.

preprint2012arXiv

How gravity and size affect the acceleration statistics of bubbles in turbulence

We report results from the first systematic Lagrangian experimental investigation in the regime of very light (air bubbles in water) and large particles (D/η >> 1) in turbulence. Using a traversing camera setup and particle tracking, we study the Lagrangian acceleration statistics of ~3 mm diameter (D) bubbles in a water tunnel with nearly homogeneous and isotropic turbulence generated by an active-grid. The Reynolds number (Re_λ) is varied from 145 to 230, resulting in size ratios, D/η in the range of 7.3--12.5, where η is the Kolmogorov length scale. The experiments reveal that gravity increases the acceleration variance and reduces the intermittency of the PDF in the vertical direction. Once the gravity offset is subtracted, the variances of both the horizontal and vertical acceleration components are about 5\pm2 times larger than the one measured in the same flow for fluid tracers. Moreover, for these light particles, the experimental acceleration PDF shows a substantial reduction in intermittency at growing size ratios, in contrast to neutrally buoyant or heavy particles. All these results are closely matched by numerical simulations of finite-size bubbles with the Faxén corrections.

preprint2012arXiv

Lagrangian statistics of light particles in turbulence

We study the Lagrangian velocity and acceleration statistics of light particles (micro-bubbles in water) in homogeneous isotropic turbulence. Micro-bubbles with a diameter of 340 microns and Stokes number from 0.02 to 0.09 are dispersed in a turbulent water tunnel operated at Taylor-Reynolds numbers (Re) ranging from 160 to 265. We reconstruct the bubble trajectories by employing three-dimensional particle tracking velocimetry (PTV). It is found that the probability density functions (PDFs) of the micro-bubble acceleration show a highly non-Gaussian behavior with flatness values in the range 23-30. The acceleration flatness values show an increasing trend with Re, consistent with previous experiments (Voth et al., JFM, 2002) and numerics (Ishihara et al., JFM, 2007). These acceleration PDFs show a higher intermittency compared to tracers (Ayyalasomayajula et al., Phys. Fluids, 2008) and heavy particles (Ayyalasomayajula et al., Phys. Rev. Lett., 2006) in wind tunnel experiments. In addition, the micro-bubble acceleration autocorrelation function decorrelates slower with increasing Re. We also compare our results with experiments in von Karman flows and point-particle direct numerical simulations with periodic boundary conditions.

preprint2012arXiv

Logarithmic temperature profiles in the ultimate regime of thermal convection

We report on the theory of logarithmic temperature profiles in very strongly developed thermal convection in the geometry of a Rayleigh-Benard cell with aspect ratio one and discuss the degree of agreement with the recently measured profiles in the ultimate state of very large Rayleigh number flow. The parameters of the log-profile are calculated and compared with the measure ones. Their physical interpretation as well as their dependence on the radial position are discussed.

preprint2012arXiv

Logarithmic temperature profiles in turbulent Rayleigh-Bénard convection

We report results for the temperature profiles of turbulent Rayleigh-Bénard convection (RBC) in the interior of a cylindrical sample of aspect ratio $Γ\equiv D/L = 0.50$ ($D$ and $L$ are the diameter and height respectively). Results from experiment over the Rayleigh number range $4\times 10^{12} \alt Ra \alt 10^{15}$ for a Prandtl number $\Pra \simeq 0.8$ and from direct numerical simulation (DNS) at $Ra = 2 \times 10^{12}$ for $\Pra = 0.7$ are presented. We find that the temperature varies as $A*ln(z/L) + B$ where $z$ is the distance from the bottom or top plate. This is the case in the classical as well as in the ultimate state of RBC. From DNS we find that $A$ in the classical state decreases in the radial direction as the distance from the side wall increases and becomes small near the sample center.

preprint2012arXiv

Maximal air bubble entrainment at liquid drop impact

At impact of a liquid drop on a solid surface an air bubble can be entrapped. Here we show that two competing effects minimize the (relative) size of this entrained air bubble: For large drop impact velocity and large droplets the inertia of the liquid flattens the entrained bubble, whereas for small impact velocity and small droplets capillary forces minimize the entrained bubble. However, we demonstrate experimentally, theoretically, and numerically that in between there is an optimum, leading to maximal air bubble entrapment. Our results have a strong bearing on various applications in printing technology, microelectronics, immersion lithography, diagnostics, or agriculture.

preprint2012arXiv

Microdroplet impact at very high velocity

Water microdroplet impact at velocities up to 100 m/s for droplet diameters from 12 to 100 um is studied. This parameter range covers the transition from capillary-limited to viscosity-limited spreading of the impacting droplet. Splashing is absent for all measurements; the droplets always gently spread over the surface. The maximum spreading radius is compared to several existing models. The model by Pasandideh-Fard et al. agrees well with the measured data, indicating the importance of a thin boundary layer just above the surface, in which most of the viscous dissipation in the spreading droplet takes place. As explained by the initial air layer under the impacting droplet, a contact angle of 180 degrees is used as model input.

preprint2012arXiv

Needle-free injection into skin and soft matter with highly focused microjets

The development of needle-free drug injection systems is of great importance to global healthcare. However, in spite of its great potential and research history over many decades, these systems are not commonly used. One of the main problems is that existing methods use diffusive jets, which result in scattered penetration and severe deceleration of the jets, causing frequent pain and insufficient penetration. Another longstanding challenge is the development of accurate small volume injections. In this paper we employ a novel method of needle-free drug injection, using highly-focused high speed microjets, which aims to solve these challenges. We experimentally demonstrate that these unique jets are able to penetrate human skin: the focused nature of these microjets creates an injection spot smaller than a mosquito's proboscis and guarantees a high percentage of the liquid being injected. The liquid substances can be delivered to a much larger depth than conventional methods, and create a well-controlled dispersion pattern. Thanks to the excellent controllability of the microjet, small volume injections become feasible. Furthermore, the penetration dynamics is studied through experiments performed on gelatin mixtures (human soft tissue equivalent) and human skin, agreeing well with a viscous stress model which we develop. This model predicts the depth of the penetration into both human skin and soft tissue. The results presented here take needle-free injections a step closer to widespread use.

preprint2012arXiv

Optimal Taylor-Couette turbulence

Strongly turbulent Taylor-Couette flow with independently rotating inner and outer cylinders with a radius ratio of η= 0.716 is experimentally studied. From global torque measurements, we analyse the dimensionless angular velocity flux Nu_ω(Ta, a) as a function of the Taylor number Ta and the angular velocity ratio a = -ω_o/ω_i in the large-Taylor-number regime 10^{11} \lesssim Ta \lesssim 10^{13}. We analyse the data with the common power-law ansatz for the dimensionless angular velocity transport flux Nu_ω(Ta, a) = f(a)Ta^γ, with an amplitude f(a) and an exponent γ. The data are consistent with one effective exponent γ= 0.39\pm0.03 for all a. The amplitude of the angular velocity flux f(a) = Nu_ω(Ta, a)/Ta^0.39 is measured to be maximal at slight counter-rotation, namely at an angular velocity ratio of a_opt = 0.33\pm0.04. This value is theoretically interpreted as the result of a competition between the destabilizing inner cylinder rotation and the stabilizing but shear-enhancing outer cylinder counter-rotation. With the help of laser Doppler anemometry, we provide angular velocity profiles and identify the radial position r_n of the neutral line. While for moderate counter-rotation -0.40 ω_i \lesssim ω_o < 0, the neutral line still remains close to the outer cylinder and the probability distribution function (p.d.f.) of the bulk angular velocity is observed to be monomodal. For stronger counter-rotation the neutral line is pushed inwards towards the inner cylinder; in this regime the p.d.f. of the bulk angular velocity becomes bimodal, reflecting intermittent bursts of turbulent structures beyond the neutral line into the outer flow domain, which otherwise is stabilized by the counter-rotating outer cylinder. Finally, a hypothesis is offered allowing a unifying view for all these various results.

preprint2012arXiv

Sonoluminescence and sonochemiluminescence from a microreactor

Micromachined pits on a substrate can be used to nucleate and stabilize microbubbles in a liquid exposed to an ultrasonic field. Under suitable conditions, the collapse of these bubbles can result in light emission (sonoluminescence, SL). Hydroxyl radicals (OH*) generated during bubble collapse can react with luminol to produce light (sonochemiluminescence, SCL). SL and SCL intensities were recorded for several regimes related to the pressure amplitude (low and high acoustic power levels) at a given ultrasonic frequency (200 kHz) for pure water, and aqueous luminol and propanol solutions. Various arrangements of pits were studied, with the number of pits ranging from no pits (comparable to a classic ultrasound reactor), to three-pits. Where there was more than one pit present, in the high pressure regime the ejected microbubbles combined into linear (two-pits) or triangular (three-pits) bubble clouds (streamers). In all situations where a pit was present on the substrate, the SL was intensified and increased with the number of pits at both low and high power levels. For imaging SL emitting regions, Argon (Ar) saturated water was used under similar conditions. SL emission from aqueous propanol solution (50 mM) provided evidence of transient bubble cavitation. Solutions containing 0.1 mM luminol were also used to demonstrate the radical production by attaining the SCL emission regions.

preprint2012arXiv

Spatial distribution of heat flux and fluctuations in turbulent Rayleigh-Bénard convection

We numerically investigate the radial dependence of the velocity and temperature fluctuations and of the time-averaged heat flux $\bar{j}(r)$ in a cylindrical Rayleigh-Bénard cell with aspect ratio Γ=1 for Rayleigh numbers Ra between $2 \times 10^6$ and $2\times 10^{9}$ at a fixed Prandtl number Pr = 5.2. The numerical results reveal that the heat flux close to the sidewall is larger than in the center and that, just as the global heat transport, it has an effective power law dependence on the Rayleigh number, $\bar{j}(r)\propto Ra^{γ_j(r)}$. The scaling exponent $γ_j(r)$ decreases monotonically from 0.43 near the axis ($r \approx 0$) to 0.29 close to the side walls ($r \approx D/2$). The effective exponents near the axis and the side wall agree well with the measurements of Shang et al. (Phys.\ Rev.\ Lett.\ \textbf{100}, 244503, 2008) and the predictions of Grossmann and Lohse (Phys.\ Fluids \textbf{16}, 1070, 2004). Extrapolating our results to large Rayleigh number would imply a crossover at $Ra\approx 10^{15}$, where the heat flux near the axis would begin to dominate. In addition, we find that the local heat flux is more than twice as high at the location where warm or cold plumes go up or down, than in the plume depleted regions.

preprint2012arXiv

Surface nanobubbles: Seeing is believing

The existence of surface nanobubbles has been previously suggested using various experimental techniques, including attenuated total reflection spectroscopy, quartz crystal microbalance, neutron reflectometry, and x-ray reflectivity, but all of these techniques provide a sole number to quantify the existence of gas over (usually) hundreds of square microns. Thus `nanobubbles' are indistinguishable from a `uniform gassy layer' between surface and liquid. Atomic force microscopy, on the other hand, does show the existence of surface nanobubbles, but the highly intrusive nature of the technique means that a uniform gassy layer could break down into nanobubbles \textit{due to} the motion of the microscope's probe. Here we demonstrate \textit{optical} visualisation of surface nanobubbles, thus validating their individual existence non-intrusively.

preprint2012arXiv

Taming Acoustic Cavitation

In this fluid dynamics video we show acoustic cavitation occurring from pits etched on a silicon surface. By immersing the surface in a liquid, gas pockets are entrapped in the pits which upon ultrasonic insonation, are observed to shed cavitation bubbles. Modulating the driving pressure it is possible to induce different behaviours based on the force balance that determines the interaction among bubbles and the silicon surface. This system can be used for several applications like sonochemical water treatment, cleaning of surfaces with deposited materials such as biofilms.

preprint2012arXiv

The clustering morphology of freely rising deformable bubbles

We investigate the clustering morphology of a swarm of freely rising deformable bubbles. A three-dimensional Voronoï analysis enables us to quantitatively distinguish between two typical clustering configurations: preferential clustering and a grid-like structure. The bubble data is obtained from direct numerical simulations (DNS) using the front-tracking method. It is found that the bubble deformation, represented by the aspect ratio χ, plays a significant role in determining which type of clustering is realized: Nearly spherical bubbles with χ<~ 1.015 form a grid-like structure, while more deformed bubbles show preferential clustering. Remarkably, this criteria for the clustering morphology holds for different diameters of the bubbles, surface tension, and viscosity of the liquid in the studied parameter regime. The mechanism of this clustering behavior is connected to the amount of vorticity generated at the bubble surfaces.

preprint2012arXiv

Thermal boundary layer profiles in turbulent Rayleigh-Bénard convection in a cylindrical sample

We numerically investigate the structures of the near-plate temperature profiles close to the bottom and top plates of turbulent Rayleigh-Bénard flow in a cylindrical sample at Rayleigh numbers Ra=10^8 to Ra=2\times10^{12} and Prandtl numbers Pr=6.4 and Pr=0.7 with the dynamical frame method [Q. Zhou and K.-Q. Xia, Phys. Rev. Lett. 104, 104301 (2010)] thus extending previous results for quasi-2-dimensional systems to 3D systems for the first time. The dynamical frame method shows that the measured temperature profiles in the spatially and temporally local frame are much closer to the temperature profile of a laminar, zero-pressure gradient boundary layer according to Pohlhausen than in the fixed reference frame. The deviation between the measured profiles in the dynamical reference frame and the laminar profiles increases with decreasing Pr, where the thermal BL is more exposed to the bulk fluctuations due to the thinner kinetic BL, and increasing Ra, where more plumes are passing the measurement location.

preprint2012arXiv

Three-dimensional Lagrangian Voronoi analysis for clustering of particles and bubbles in turbulence

Three-dimensional Voronoi analysis is used to quantify the clustering of inertial particles in homogeneous isotropic turbulence using data from numerics and experiments. We study the clustering behavior at different density ratios and particle response times (i.e. Stokes numbers St). The Probability Density Functions (PDFs) of the Voronoi cell volumes of light and heavy particles show a different behavior from that of randomly distributed particles -i.e. fluid tracers-implying that clustering is present. The standard deviation of the PDF normalized by that of randomly distributed particles is used to quantify the clustering. Light particles show maximum clustering for St around 1-2. The results are consistent with previous investigations employing other approaches to quantify the clustering. We also present the joint PDFs of enstrophy and Voronoi volumes and their Lagrangian autocorrelations. The small Voronoi volumes of light particles correspond to regions of higher enstrophy than those of heavy particles, indicating that light particles cluster in higher vorticity regions. The Lagrangian temporal autocorrelation function of Voronoi volumes shows that the clustering of light particles lasts much longer than that of heavy or neutrally buoyant particles. Due to inertial effects, the Lagrangian autocorrelation time-scale of clustered light particles is even longer than that of the enstrophy of the flow itself.

preprint2012arXiv

Why surface nanobubbles live for hours

We present a theoretical model for the experimentally found but counter-intuitive exceptionally long lifetime of surface nanobubbles. We can explain why, under normal experimental conditions, surface nanobubbles are stable for many hours or even up to days rather than the expected microseconds. The limited gas diffusion through the water in the far field, the cooperative effect of nanobubble clusters, and the pinned contact line of the nanobubbles lead to the slow dissolution rate.

preprint2011arXiv

Angular momentum transport and turbulence in laboratory models of Keplerian flows

We present angular momentum transport (torque) measurements in two recent experimental studies of the turbulent flow between independently rotating cylinders. In addition to these studies, we reanalyze prior torque measurements to expand the range of control parameters for the experimental Taylor-Couette flows. We find that the torque may be described as a product of functions that depend only on the Reynolds number, which describes the turbulent driving intensity, and the rotation number, which characterizes the effects of global rotation. For a given Reynolds number, the global angular momentum transport for Keplerian-like flow profiles is approximately 14% of the maximum achievable transport rate. We estimate that this level of transport would produce an accretion rate of $\dot{M}/\dot{M_0} \sim 10^{-3}$ in astrophysical disks. We argue that this level of transport from hydrodynamics alone could be significant.

preprint2011arXiv

Axially-homogeneous Rayleigh-Benard convection in a cylindrical cell

Previous numerical studies have shown that the "ultimate regime of thermal convection" can be attained in a Rayleigh-Benard cell when the kinetic and thermal boundary layers are eliminated by replacing the walls with periodic boundary conditions (homogeneous Rayleigh-Benard convection). Then, the heat transfer scales like Nu ~ Ra^{1/2} and turbulence intensity as Re ~ Ra^{1/2}, where the Rayleigh number Ra indicates the strength of the driving force. However, experiments never operate in unbounded domains and it is important to understand how confinement might alter the approach to this ultimate regime. Here we consider homogeneous Rayleigh-Benard convection in a laterally confined geometry - a small aspect-ratio vertical cylindrical cell - and show evidence of the ultimate regime as Ra is increased: In spite of the confinement and the resulting kinetic boundary layers, we still find Nu ~ Re ~ Ra^{1/2}. The system supports exact solutions composed of modes of exponentially growing vertical velocity and temperature fields, with Ra as the critical parameter determining the properties of these modes. Counterintuitively, in the low Ra regime, or for very narrow cylinders, the numerical simulations are susceptible to these solutions which can dominate the dynamics and lead to very high and unsteady heat transfer. As Ra is increased, interaction between modes stabilizes the system, evidenced by the increasing homogeneity and reduced fluctuations in the r.m.s. velocity and temperature fields. We also test that physical results become independent of the periodicity length of the cylinder, a purely numerical parameter, as the aspect ratio is increased.

preprint2011arXiv

Boundary layer structure in turbulent thermal convection and its consequences for the required numerical resolution

Results on the Prandtl-Blasius type kinetic and thermal boundary layer thicknesses in turbulent Rayleigh-Bénard convection in a broad range of Prandtl numbers are presented. By solving the laminar Prandtl-Blasius boundary layer equations, we calculate the ratio of the thermal and kinetic boundary layer thicknesses, which depends on the Prandtl number Pr only. It is approximated as $0.588Pr^{-1/2}$ for $Pr\ll Pr^*$ and as $0.982 Pr^{-1/3}$ for $Pr^*\ll\Pr$, with $Pr^*= 0.046$. Comparison of the Prandtl--Blasius velocity boundary layer thickness with that evaluated in the direct numerical simulations by Stevens, Verzicco, and Lohse (J. Fluid Mech. 643, 495 (2010)) gives very good agreement. Based on the Prandtl--Blasius type considerations, we derive a lower-bound estimate for the minimum number of the computational mesh nodes, required to conduct accurate numerical simulations of moderately high (boundary layer dominated) turbulent Rayleigh-Bénard convection, in the thermal and kinetic boundary layers close to bottom and top plates. It is shown that the number of required nodes within each boundary layer depends on Nu and Pr and grows with the Rayleigh number Ra not slower than $\sim\Ra^{0.15}$. This estimate agrees excellently with empirical results, which were based on the convergence of the Nusselt number in numerical simulations.

preprint2011arXiv

Bubble size prediction in co-flowing streams

In this paper, the size of bubbles formed through the breakup of a gaseous jet in a co-axial microfluidic device is derived. The gaseous jet surrounded by a co-flowing liquid stream breaks up into monodisperse microbubbles and the size of the bubbles is determined by the radius of the inner gas jet and the bubble formation frequency. We obtain the radius of the gas jet by solving the Navier-Stokes equations for low Reynolds number flows and by minimization of the dissipation energy. The prediction of the bubble size is based on the system's control parameters only, i.e. the inner gas flow rate $Q_i$, the outer liquid flow rate $Q_o$, and the tube radius $R$. For a very low gas-to-liquid flow rate ratio ($Q_i / Q_o \rightarrow 0$) the bubble radius scales as $r_b / R \propto \sqrt{Q_i / Q_o}$, independently of the inner to outer viscosity ratio $η_i/η_o$ and of the type of the velocity profile in the gas, which can be either flat or parabolic, depending on whether high-molecular-weight surfactants cover the gas-liquid interface or not. However, in the case in which the gas velocity profiles are parabolic and the viscosity ratio is sufficiently low, i.e. $η_i/η_o \ll 1$, the bubble diameter scales as $r_b \propto (Q_i/Q_o)^β$, with $β$ smaller than 1/2.

preprint2011arXiv

Connecting flow structures and heat flux in turbulent Rayleigh-Bénard convection

The aspect ratio (Γ) dependence of the heat transfer (Nusselt number Nu in dimensionless form) in turbulent (two-dimensional) Rayleigh-Bénard convection is numerically studied in the regime $0.4 \le Γ\le 1.25$ for Rayleigh numbers $10^7 \le Ra \le Ra^{9}$ and Prandtl numbers Pr =0.7 (gas) and 4.3 (water). Nu (Γ) shows a very rich structure with sudden jumps and sharp transitions. We connect these structures to the way the flow organizes itself in the sample and explain why the aspect ratio dependence of Nu is more pronounced for small Pr. Even for fixed Γdifferent turbulent states (with different resulting Nu) can exist, between which the flow can or cannot switch. In the latter case the heat transfer thus depends on the initial conditions.

preprint2011arXiv

Direct measurements of air layer profiles under impacting droplets using high-speed color interferometry

A drop impacting on a solid surface deforms before the liquid makes contact with the surface. We directly measure the time evolution of the air layer profile under the droplet using high-speed color interferometry, obtaining the air layer thickness before and during the wetting process. Based on the time evolution of the extracted profiles obtained at multiple times, we measure the velocity of air exiting from the gap between the liquid and the solid, and account for the wetting mechanism and bubble entrapment. The present work offers a tool to accurately measure the air layer profile and quantitatively study the impact dynamics at a short time scale before impact.

preprint2011arXiv

Drop impact on superheated surfaces

At impact of a liquid droplet on a smooth surface heated above the liquid's boiling point, the droplet either immediately boils when it contacts the surfaces (``contact boiling''), or without any surface contact forms a Leidenfrost vapor layer towards the hot surface and bounces back (``gentle film boiling''), or both forms the Leidenfrost layer and ejects tiny droplets upward (``spraying film boiling''). We experimentally determine conditions under which impact behaviors in each regime can be realized. We show that the dimensionless maximum spreading $γ$ of impacting droplets on the heated surfaces in both gentle and spraying film boiling regimes shows a universal scaling with the Weber number $\We$ ($γ\sim\We^{2/5}$) -- regardless of surface temperature and of liquid properties -- which is much steeper than for the impact on non-heated (hydrophilic or hydrophobic) surfaces ($γ\sim\We^{1/4}$). We also intereferometrically measure the vapor thickness under the droplet.

preprint2011arXiv

Effect of aspect-ratio on vortex distribution and heat transfer in rotating Rayleigh-Bénard convection

Numerical and experimental data for the heat transfer as function of the Rossby number Ro in turbulent rotating Rayleigh-Bénard convection are presented for Prandtl number Pr=4.38 and Rayleigh number $Ra=2.91\times10^8$ up to $Ra=4.52\times10^9$. The aspect ratio Γ= D/L, where L is the height and D the diameter of the cylindrical sample, is varied between Γ=0.5 and Γ=2.0. Without rotation, where the aspect ratio influences the global large scale circulation, we see a small aspect-ratio dependence in the Nusselt number for $Ra=2.91\times10^8$. However, for stronger rotation, i.e. $1/Ro \gg 1/Ro_c$, the heat transport becomes independent of the aspect-ratio. We interpret this finding as follows: In the rotating regime the heat is mainly transported by vertically-aligned vortices. Since the vertically-aligned vortices are local, the aspect ratio has a negligible effect on the heat transport in the rotating regime. Indeed, a detailed analysis of vortex statistics shows that the fraction of the horizontal area that is covered by vortices is independent of the aspect ratio when $1/Ro \gg 1/Ro_c$. In agreement with the results of Weiss et al. Phys. Rev. Lett., vol 105, 224501 (2010) we find a vortex-depleted area close to the sidewall. Here, we in addition show that there is also an area with enhanced vortex concentration next to the vortex-depleted edge region and that the absolute widths of both regions are independent of the aspect ratio.

preprint2011arXiv

Effect of Impurities in Description of Surface Nanobubbles

Surface nanobubbles emerging at solid-liquid interfaces of submerged hydrophobic surfaces show extreme stability and very small (gas-side) contact angles. In a recent study Ducker (W. A. Ducker, Langmuir 25, 8907 (2009).) conjectured that these effects may arise from the presence of impurities at the air-water interface of the nanobubbles. In this paper we present a quantitative analysis of this hypothesis by estimating the dependence of the contact angle and the Laplace pressure on the fraction of impurity coverage at the liquid-gas interface. We first develop a general analytical framework to estimate the effect of impurities (ionic or non-ionic) in lowering the surface tension of a given air-water interface. We then employ this model to show that the (gas-side) contact angle and the Laplace pressure across the nanobubbles indeed decrease considerably with an increase in the fractional coverage of the impurities, though still not sufficiently small to account for the observed surface nanobubble stability. The proposed model also suggests the dependencies of the Laplace pressure and the contact angle on the type of impurity.

preprint2011arXiv

Effect of Plumes on Measuring the Large Scale Circulation in Turbulent Rayleigh-Bénard Convection

We studied the properties of the large-scale circulation (LSC) in turbulent Rayleigh-Bénard (RB) convection by using results from direct numerical simulations in which we placed a large number of numerical probes close to the sidewall. The LSC orientation is determined by either a cosine or a polynomial fit to the azimuthal temperature or azimuthal vertical velocity profile measured with the probes. We study the LSC in Γ=D/L=1/2 and Γ=1 samples, where D is the diameter and L the height. For Pr=6.4 in an aspect ratio Γ=1 sample at $Ra=1\times10^8$ and $5\times10^8$ the obtained LSC orientation is the same, irrespective of whether the data of only 8 or all 64 probes per horizontal plane are considered. In a Γ=1/2 sample with $Pr=0.7$ at $Ra=1\times10^8$ the influence of plumes on the azimuthal temperature and azimuthal vertical velocity profiles is stronger. Due to passing plumes and/or the corner flow the apparent LSC orientation obtained using a cosine fit can result in a misinterpretation of the character of the large-scale flow. We introduce the relative LSC strength, which we define as the ratio between the energy in the first Fourier mode and the energy in all modes that can be determined from the azimuthal temperature and azimuthal vertical velocity profiles, to further quantify the large-scale flow. For $Ra=1\times10^8$ we find that this relative LSC strength is significantly lower in a Γ=1/2 sample than in a Γ=1 sample, reflecting that the LSC is much more pronounced in a Γ=1 sample than in a Γ=1/2 sample. The determination of the relative LSC strength can be applied directly to available experimental data to study high Rayleigh number thermal convection and rotating RB convection.

preprint2011arXiv

Finite-size effects lead to supercritical bifurcations in turbulent rotating Rayleigh-Bénard convection

In turbulent thermal convection in cylindrical samples of aspect ratio Γ= D/L (D is the diameter and L the height) the Nusselt number Nu is enhanced when the sample is rotated about its vertical axis, because of the formation of Ekman vortices that extract additional fluid out of thermal boundary layers at the top and bottom. We show from experiments and direct numerical simulations that the enhancement occurs only above a bifurcation point at a critical inverse Rossby number $1/\Ro_c$, with $1/\Ro_c \propto 1/Γ$. We present a Ginzburg-Landau like model that explains the existence of a bifurcation at finite $1/\Ro_c$ as a finite-size effect. The model yields the proportionality between $1/\Ro_c$ and $1/Γ$ and is consistent with several other measured or computed system properties.

preprint2011arXiv

Flow reversals in thermally driven turbulence

We analyze the reversals of the large scale flow in Rayleigh-Bénard convection both through particle image velocimetry flow visualization and direct numerical simulations (DNS) of the underlying Boussinesq equations in a (quasi) two-dimensional, rectangular geometry of aspect ratio 1. For medium Prandtl number there is a diagonal large scale convection roll and two smaller secondary rolls in the two remaining corners diagonally opposing each other. These corner flow rolls play a crucial role for the large scale wind reversal: They grow in kinetic energy and thus also in size thanks to plume detachments from the boundary layers up to the time that they take over the main, large scale diagonal flow, thus leading to reversal. Based on this mechanism we identify a typical time scale for the reversals. We map out the Rayleigh number vs Prandtl number phase space and find that the occurrence of reversals very sensitively depends on these parameters.

preprint2011arXiv

Horizontal Structures of Velocity and Temperature Boundary Layers in 2D Numerical Turbulent Rayleigh-Bénard Convection

We investigate the structures of the near-plate velocity and temperature profiles at different horizontal positions along the conducting bottom (and top) plate of a Rayleigh-Bénard convection cell, using two-dimensional (2D) numerical data obtained at the Rayleigh number Ra=10^8 and the Prandtl number Pr=4.4 of an Oberbeck-Boussinesq flow with constant material parameters. The results show that most of the time, and for both velocity and temperature, the instantaneous profiles scaled by the dynamical frame method [Q. Zhou and K.-Q. Xia, Phys. Rev. Lett. 104, 104301 (2010) agree well with the classical Prandtl-Blasius laminar boundary layer (BL) profiles. Therefore, when averaging in the dynamical reference frames, which fluctuate with the respective instantaneous kinematic and thermal BL thicknesses, the obtained mean velocity and temperature profiles are also of Prandtl-Blasius type for nearly all horizontal positions. We further show that in certain situations the traditional definitions based on the time-averaged profiles can lead to unphysical BL thicknesses, while the dynamical method also in such cases can provide a well-defined BL thickness for both the kinematic and the thermal BLs.

preprint2011arXiv

How water droplets evaporate on a superhydrophobic substrate

Evaporation of water droplets on a superhydrophobic substrate, on which the contact line is pinned, is investigated. While previous studies mainly focused on droplets with contact angles smaller than 90^\circ, here we analyze almost the full range of possible contact angles (10^\circ -150^\circ). The greater contact angles and pinned contact lines can be achieved by the use of superhydrophobic Carbon Nanofiber substrates. The time-evolutions of the contact angle and the droplet mass are examined. The experimental data is in good quantitative agreement with the model presented by Popov (Physical Review E 71, 2005), demonstrating that the evaporation process is quasi-static, diffusion-driven, and that thermal effects play no role. Furthermore, we show that the experimental data for the evolution of both the contact angle and the droplet mass can be collapsed onto one respective universal curve for all droplet sizes and initial contact angles.

preprint2011arXiv

Knudsen gas provides nanobubble stability

We provide a model for the remarkable stability of surface nanobubbles to bulk dissolution. The key to the solution is that the gas in a nanobubble is of Knudsen type. This leads to the generation of a bulk liquid flow which effectively forces the diffusive gas to remain local. Our model predicts the presence of a vertical water jet immediately above a nanobubble, with an estimated speed of $\sim3.3\,\mathrm{m/s}$, in good agreement with our experimental atomic force microscopy measurement of $\sim2.7\,\mathrm{m/s}$. In addition, our model also predicts an upper bound for the size of nanobubbles, which is consistent with the available experimental data.

preprint2011arXiv

Microbubble formation and pinch-off scaling exponent in flow-focusing devices

We investigate the gas jet breakup and the resulting microbubble formation in a microfluidic flow-focusing device using ultra high-speed imaging at 1 million frames/s. In recent experiments [Dollet et al., Phys. Rev. Lett. 100, 034504 (2008)] it was found that in the final stage of the collapse the radius of the neck scales with time with a 1/3 power-law exponent, which suggested that gas inertia and the Bernoulli suction effect become important. Here, ultra high-speed imaging was used to capture the complete bubble contour and quantify the gas flow through the neck. It revealed that the resulting decrease in pressure, due to Bernoulli suction, is too low to account for an accelerated pinch-off. The high temporal resolution images enable us to approach the final moment of pinch-off to within 1 μs. We observe that the final moment of bubble pinch-off is characterized by a scaling exponent of 0.41 +/- 0.01. This exponent is approximately 2/5, which can be derived, based on the observation that during the collapse the neck becomes less slender, due to the exclusive driving through liquid inertia.

preprint2011arXiv

Order-to-disorder transition in ring-shaped colloidal stains

A colloidal dispersion droplet evaporating from a surface, such as a drying coffee drop, leaves a distinct ring-shaped stain. Although this mechanism is frequently used for particle self-assembly, the conditions for crystallization have remained unclear. Our experiments with monodisperse colloidal particles reveal a structural transition in the stain, from ordered crystals to disordered packings. We show that this sharp transition originates from a temporal singularity of the flow velocity inside the evaporating droplet at the end of its life. When the deposition speed is low, particles have time to arrange by Brownian motion, while at the end, high-speed particles are jammed into a disordered phase.

preprint2011arXiv

Origin of line tension for a Lennard-Jones nanodroplet

The existence and origin of line tension has remained controversial in literature. To address this issue we compute the shape of Lennard-Jones nanodrops using molecular dynamics and compare them to density functional theory in the approximation of the sharp kink interface. We show that the deviation from Young's law is very small and would correspond to a typical line tension length scale (defined as line tension divided by surface tension) similar to the molecular size and decreasing with Young's angle. We propose an alternative interpretation based on the geometry of the interface at the molecular scale.

preprint2011arXiv

Prandtl and Rayleigh number dependence of heat transport in high Rayleigh number thermal convection

Results from direct numerical simulation for three-dimensional Rayleigh-Bénard convection in samples of aspect ratio $Γ=0.23$ and $Γ=0.5$ up to Rayleigh number $Ra=2\times10^{12}$ are presented. The broad range of Prandtl numbers $0.5<Pr<10$ is considered. In contrast to some experiments, we do not see any increase in $Nu/Ra^{1/3}$, neither due to $Pr$ number effects, nor due to a constant heat flux boundary condition at the bottom plate instead of constant temperature boundary conditions. Even at these very high $Ra$, both the thermal and kinetic boundary layer thicknesses obey Prandtl-Blasius scaling.

preprint2011arXiv

The Twente turbulent Taylor-Couette (T3C) facility: Strongly turbulent (multiphase) flow between two independently rotating cylinders

A new turbulent Taylor-Couette system consisting of two independently rotating cylinders has been constructed. The gap between the cylinders has a height of 0.927 m, an inner radius of 0.200 m, and a variable outer radius (from 0.279 to 0.220 m). The maximum angular rotation rates of the inner and outer cylinder are 20 and 10 Hz, respectively, resulting in Reynolds numbers up to 3.4 x 10^6 with water as working fluid. With this Taylor-Couette system, the parameter space (Re_i, Re_o, η) extends to (2.0 x 10^6, {\pm}1.4 x 10^6, 0.716-0.909). The system is equipped with bubble injectors, temperature control, skin-friction drag sensors, and several local sensors for studying turbulent single-phase and two-phase flows. Inner cylinder load cells detect skin-friction drag via torque measurements. The clear acrylic outer cylinder allows the dynamics of the liquid flow and the dispersed phase (bubbles, particles, fibers, etc.) inside the gap to be investigated with specialized local sensors and nonintrusive optical imaging techniques. The system allows study of both Taylor-Couette flow in a high-Reynolds-number regime, and the mechanisms behind skin-friction drag alterations due to bubble injection, polymer injection, and surface hydrophobicity and roughness.

preprint2011arXiv

Torque Scaling in Turbulent Taylor-Couette Flow with Co- and Counterrotating Cylinders

We analyze the global transport properties of turbulent Taylor-Couette flow in the strongly turbulent regime for independently rotating outer and inner cylinders, reaching Reynolds numbers of the inner and outer cylinders of Re_i = 2 x 10^6 and Re_o = 1.4 x 10^6, respectively. For all Re_i, Re_o, the dimensionless torque G scales as a function of the Taylor number Ta (which is proportional to the square of the difference between the angular velocities of the inner and outer cylinders) with a universal effective scaling law G \propto Ta^{0.88}, corresponding to Nu_omega \propto Ta^{0.38} for the Nusselt number characterizing the angular velocity transport between the inner and outer cylinders. The exponent 0.38 corresponds to the ultimate regime scaling for the analogous Rayleigh-Benard system. The transport is most efficient for the counterrotating case along the diagonal in phase space with omega_o \approx -0.4 omega_i.

preprint2011arXiv

Ultimate Turbulent Taylor-Couette Flow

The flow structure of strongly turbulent Taylor-Couette flow with Reynolds numbers up to Re_i = 2*10^6 of the inner cylinder is experimentally examined with high-speed particle image velocimetry (PIV). The wind Reynolds numbers Re_w of the turbulent Taylor-vortex flow is found to scale as Re_w ~ Ta^(1/2), exactly as predicted for the ultimate turbulence regime, in which the boundary layers are turbulent. The dimensionless angular velocity flux has an effective scaling of Nu_ω ~ Ta^0.38, also in correspondence with turbulence in the ultimate regime. The scaling of Nu_ω is confirmed by local angular velocity flux measurements extracted from high-speed PIV measurements: though the flux shows huge fluctuations, its spatial and temporal average nicely agrees with the result from the global torque measurements.

preprint2010arXiv

Avalanche of particles in evaporating coffee drops

The pioneering work of Deegan et al. [Nature 389, (1997)] showed how a drying sessile droplet suspension of particles presents a maximum evaporating flux at its contact line which drags liquid and particles creating the well known coffee stain ring. In this Fluid Dynamics Video, measurements using micro Particle Image Velocimetry and Particle Tracking clearly show an avalanche of particles being dragged in the last moments, for vanishing contact angles and droplet height. This explains the different characteristic packing of the particles in the layers of the ring: the outer one resembles a crystalline array, while the inner one looks more like a jammed granular fluid. Using the basic hydrodynamic model used by Deegan et al. [Phys. Rev. E 62, (2000)] it will be shown how the liquid radial velocity diverges as the droplet life comes to an end, yielding a good comparison with the experimental data.

preprint2010arXiv

Boundary layers in rotating weakly turbulent Rayleigh-Benard convection

The effect of rotation on the boundary layers (BLs) in a Rayleigh-Benard (RB) system at a relatively low Rayleigh number, i.e. $Ra = 4\times10^7$, is studied for different Pr by direct numerical simulations and the results are compared with laminar BL theory. In this regime we find a smooth onset of the heat transfer enhancement as function of increasing rotation rate. We study this regime in detail and introduce a model based on the Grossmann-Lohse theory to describe the heat transfer enhancement as function of the rotation rate for this relatively low Ra number regime and weak background rotation $Ro\gtrsim 1$. The smooth onset of heat transfer enhancement observed here is in contrast to the sharp onset observed at larger $Ra \gtrsim 10^8$ by Stevens {\it{et al.}} [Phys. Rev. Lett. {\bf{103}}, 024503, 2009], although only a small shift in the Ra-Ro-Pr phase space is involved.

preprint2010arXiv

Breakup of diminutive Rayleigh jets

Discharging a liquid from a nozzle at sufficient large velocity leads to a continuous jet that due to capillary forces breaks up into droplets. Here we investigate the formation of microdroplets from the breakup of micron-sized jets with ultra high-speed imaging. The diminutive size of the jet implies a fast breakup time scale $τ_\mathrm{c} = \sqrt{ρr^3 / γ}$ of the order of 100\,ns{}, and requires imaging at 14 million frames per second. We directly compare these experiments with a numerical lubrication approximation model that incorporates inertia, surface tension, and viscosity [Eggers and Dupont, J. Fluid Mech. 262, 205 (1994); Shi, Brenner, and Nagel, Science 265, 219 (1994)]. The lubrication model allows to efficiently explore the parameter space to investigate the effect of jet velocity and liquid viscosity on the formation of satellite droplets. In the phase diagram we identify regions where the formation of satellite droplets is suppressed. We compare the shape of the droplet at pinch-off between the lubrication approximation model and a boundary integral (BI) calculation, showing deviations at the final moment of the pinch-off. Inspite of this discrepancy, the results on pinch-off times and droplet and satellite droplet velocity obtained from the lubrication approximation agree with the high-speed imaging results.

preprint2010arXiv

Building water bridges in air: Electrohydrodynamics of the Floating Water Bridge

The interaction of electrical fields and liquids can lead to phenomena that defies intuition. Some famous examples can be found in Electrohydrodynamics as Taylor cones, whipping jets or non-coalescing drops. A less famous example is the Floating Water Bridge: a slender thread of water held between two glass beakers in which a high voltage difference is applied. Surprisingly, the water bridge defies gravity even when the beakers are separated at distances up to 2 cm. In the presentation, experimental measurements and simple models are proposed and discussed for the stability of the bridge and the source of the flow, revealing an important role of polarization forces on the stability of the water bridge. On the other hand, the observed flow can only be explained due to the non negligible free charge present in the surface. In this sense, the Floating Water Bridge can be considered as an extreme case of a leaky dielectric liquid (J. R. Melcher and G. I. Taylor, Annu. Rev. Fluid Mech., 1:111, 1969).

preprint2010arXiv

Drag and lift forces on a counter-rotating cylinder in rotating flow

Results are reported of an experimental investigation into the motion of a heavy cylinder free to move inside a water-filled drum rotating around a horizontal axis. The cylinder is observed to either co- or, counter intuitively, counter-rotate with respect to the rotating drum. The flow was measured with particle image velocimetry (PIV), and it was found that the inner cylinder significantly altered the bulk flow field from the solid-body rotation found for a fluid filled drum. In the counter-rotation case, the generated lift force allowed the cylinder to freely rotate without contact with the drum wall. Drag and lift coefficients of the freely counter-rotating cylinder were measured over a wide range of Reynolds numbers, 2,500 $<$ Re $<$ 25,000, dimensionless rotation rates, 0.0$ < α< $1.2, and gap to cylinder diameter ratios 0.003 $< G/2a <$ 0.5. Drag coefficients were consistent with previous measurements on a cylinder in a uniform flow. However, for the lift coefficient considerable larger values were observed in the present measurements. We found the enhancement of the lift force to be mainly caused by the vicinity of the wall.

preprint2010arXiv

How micropatterns and air pressure affect splashing on surfaces

We experimentally investigate the splashing mechanism of a millimeter-sized ethanol drop impinging on a structured solid surface, comprised of micro-pillars, through side-view and top-view high speed imaging. By increasing the impact velocity we can tune the impact outcome from a gentle deposition to a violent splash, at which tiny droplets are emitted as the liquid sheet spreads laterally. We measure the splashing threshold for different micropatterns and find that the arrangement of the pillars significantly affects the splashing outcome. In particular, directional splashing in direction in which air flow through pattern is possible. Our top-view observations of impact dynamics reveal that an trapped air is responsible for the splashing. Indeed by lowering the pressure of the surrounding air we show that we can suppress the splashing in the explored parameter regime.

preprint2010arXiv

Prandtl-Blasius temperature and velocity boundary layer profiles in turbulent Rayleigh-Bénard convection

The shape of velocity and temperature profiles near the horizontal conducting plates in turbulent Rayleigh-Bénard convection are studied numerically and experimentally over the Rayleigh number range $10^8\lesssim Ra\lesssim3\times10^{11}$ and the Prandtl number range $0.7\lesssim Pr\lesssim5.4$. The results show that both the temperature and velocity profiles well agree with the classical Prandtl-Blasius laminar boundary-layer profiles, if they are re-sampled in the respective dynamical reference frames that fluctuate with the instantaneous thermal and velocity boundary-layer thicknesses.

preprint2010arXiv

Supersonic Air Flow due to Solid-Liquid Impact

A solid object impacting on liquid creates a liquid jet due to the collapse of the impact cavity. Using visualization experiments with smoke particles and multiscale simulations we show that in addition a high-speed air-jet is pushed out of the cavity. Despite an impact velocity of only 1 m/s, this air-jet attains \emph{supersonic} speeds already when the cavity is slightly larger than 1 mm in diameter. The structure of the air flow resembles closely that of compressible flow through a nozzle -- with the key difference that here the "nozzle" is a \emph{liquid} cavity shrinking rapidly in time.

preprint2010arXiv

Surface bubble nucleation phase space

Recent research has revealed several different techniques for nanoscopic gas nucleation on submerged surfaces, with findings seemingly in contradiction with each other. In response to this, we have systematically investigated the occurrence of surface nanobubbles on a hydrophobised silicon substrate for various different liquid temperatures and gas concentrations, which we controlled independently. We found that nanobubbles occupy a distinct region of this phase space, occurring for gas concentrations of approximately 100-110%. Below the nanobubble phase we did not detect any gaseous formations on the substrate, whereas micropancakes (micron wide, nanometer high gaseous domains) were found at higher temperatures and gas concentrations. We moreover find that supersaturation of dissolved gases is not a requirement for nucleation of bubbles.

preprint2010arXiv

Wetting Splashing

We present fluid dynamics videos illustrating wetting splashing-produced by water drop impact onto hydrophobic microstructures at high impact velocity ($\sim 3$ ms$^{-1}$). The substrate consists of regular and transparent microtextures in square or hexagonal lattice, yielding a high contact angle of $\sim 150 ^\circ$. Our high speed top-or-bottom view recordings through the transparent surface shed light on the solid-liquid-air interfaces at impact. Despite the superhydrophobicity of the latticed micropillars (of a periodicity of $10 μm), water droplet wets a certain central area and moreover entraps an air bubble beneath the droplet. Besides the central wet area, lamella surf on air splashing outward at high impinging velocity. The effects of micropatterns and air pressure on the impact outcome are also examined. We show that microscopic boundary condition, imposed by the solid texture, profoundly influences the macroscopic flow dynamics upon superhydrophobic surfaces at high impinging velocity. In addition, the intervening air between the liquid and the solid plays a crucial role in directional splash, which can be eliminated by a reduced air pressure.

preprint2009arXiv

Collapse of Non-Axisymmetric Cavities

A round disk with a harmonic disturbance impacts on a water surface and creates a non-axisymmetric cavity which collapses under the influence of hydrostatic pressure. We use disks deformed with mode m=2 to m=6. For all mode numbers we find clear evidence for a phase inversion of the cavity wall during the collapse. We present a fluid dynamics video showing high speed imaging of different modes, pointing out the characteristic features during collapse.

preprint2009arXiv

Drop impact upon micro- and nanostructured superhydrophobic surfaces

We experimentally investigate drop impact dynamics onto different superhydrophobic surfaces, consisting of regular polymeric micropatterns and rough carbon nanofibers, with similar static contact angles. The main control parameters are the Weber number \We and the roughness of the surface. At small \We, i.e. small impact velocity, the impact evolutions are similar for both types of substrates, exhibiting Fakir state, complete bouncing, partial rebouncing, trapping of an air bubble, jetting, and sticky vibrating water balls. At large \We, splashing impacts emerge forming several satellite droplets, which are more pronounced for the multiscale rough carbon nanofiber jungles. The results imply that the multiscale surface roughness at nanoscale plays a minor role in the impact events for small \We $\apprle 120$ but an important one for large \We $\apprge 120$. Finally, we find the effect of ambient air pressure to be negligible in the explored parameter regime \We $\apprle 150$

preprint2009arXiv

Evaporation-triggered Wetting Transition for Water Droplets upon Hydrophobic Microstructures

When placed on rough hydrophobic surfaces, water droplets of diameter larger than a few millimeters can easily form pearls, as they are in the Cassie-Baxter state with air pockets trapped underneath the droplet. Intriguingly, a natural evaporating process can drive such a Fakir drop into a completely wetting (Wenzel) state. Our microscopic observations with simultaneous side and bottom views of evaporating droplets upon transparent hydrophobic microstructures elucidate the water-filling dynamics and the mechanism of this evaporation-triggered transition. For the present material the wetting transition occurs when the water droplet size decreases to a few hundreds of micrometers in radius. We present a general global energy argument which estimates the interfacial energies depending on the drop size and can account for the critical radius for the transition.

preprint2009arXiv

Nucleation threshold and deactivation mechanisms of nanoscopic cavitation nuclei

The acoustic nucleation threshold for bubbles trapped in cavities has theoretically been predicted within the crevice theory by Atchley & Prosperetti [J. Acoust. Soc. Am. 86, 1065-1084 (1989)]. Here, we determine this threshold experimentally, by applying a single pressure pulse to bubbles trapped in cylindrical nanoscopic pits ("artificial crevices") with radii down to 50 nm. By decreasing the minimum pressure stepwise, we observe the threshold for which the bubbles start to nucleate. The experimental results are quantitatively in excellent agreement with the theoretical predictions of Atchley & Prosperetti. In addition, we provide the mechanism which explains the deactivation of cavitation nuclei: gas diffusion together with an aspherical bubble collapse. Finally, we present superhydrophobic nuclei which cannot be deactivated, unless with a high-speed liquid jet directed into the pit.

preprint2009arXiv

Optimal Prandtl number for heat transfer in rotating Rayleigh-Benard convection

Numerical data for the heat transfer as a function of the Prandtl (Pr) and Rossby (Ro) numbers in turbulent rotating Rayleigh-Benard convection are presented for Rayleigh number Ra = 10^8. When Ro is fixed the heat transfer enhancement with respect to the non-rotating value shows a maximum as function of Pr. This maximum is due to the reduced efficiency of Ekman pumping when Pr becomes too small or too large. When Pr becomes small, i.e. for large thermal diffusivity, the heat that is carried by the vertical vortices spreads out in the middle of the cell, and Ekman pumping thus becomes less efficient. For higher Pr the thermal boundary layers (BLs) are thinner than the kinetic BLs and therefore the Ekman vortices do not reach the thermal BL. This means that the fluid that is sucked into the vertical vortices is colder than for lower Pr which limits the efficiency of the upwards heat transfer.

preprint2009arXiv

Preferred sizes and ordering in surface nanobubble populations

Two types of homogeneous surface nanobubble populations, created by different means, are analyzed statistically on both their sizes and spatial positions. In the first type (created by droplet-deposition, case A) the bubble size R is found to be distributed according to a generalized gamma law with a preferred radius R*=20 nm. The radial distribution function shows a preferred spacing at ~5.5 R*. These characteristics do not show up in comparable Monte-Carlo simulations of random packings of hard disks with the same size distribution and the same density, suggesting a structuring effect in the nanobubble formation process. The nanobubble size distribution of the second population type (created by ethanol-water exchange, case B) is a mixture of two clearly separated distributions, hence, with two preferred radii. The local ordering is less significant, due to the looser packing of the nanobubbles.

preprint2009arXiv

Quantifying effective slip length over micropatterned hydrophobic surfaces

We employ micro-particle image velocimetry ($μ$-PIV) to investigate laminar micro-flows in hydrophobic microstructured channels, in particular the slip length. These microchannels consist of longitudinal micro-grooves, which can trap air and prompt a shear-free boundary condition and thus slippage enhancement. Our measurements reveal an increase of the slip length when the width of the micro-grooves is enlarged. The result of the slip length is smaller than the analytical prediction by Philip et al. [1] for an infinitely large and textured channel comprised of alternating shear-free and no-slip boundary conditions. The smaller slip length (as compared to the prediction) can be attributed to the confinement of the microchannel and the bending of the meniscus (liquid-gas interface). Our experimental studies suggest that the curvature of the meniscus plays an important role in microflows over hydrophobic micro-ridges.

preprint2009arXiv

Radial boundary layer structure and Nusselt number in Rayleigh-Benard convection

Results from direct numerical simulations for three dimensional Rayleigh-Benard convection in a cylindrical cell of aspect ratio 1/2 and Pr=0.7 are presented. They span five decades of Ra from $2\times 10^6$ to $2 \times10^{11}$. Good numerical resolution with grid spacing $\sim$ Kolmogorov scale turns out to be crucial to accurately calculate the Nusselt number, which is in good agreement with the experimental data by Niemela et al., Nature, 404, 837 (2000). In underresolved simulations the hot (cold) plumes travel further from the bottom (top) plate than in the fully resolved case, because the thermal dissipation close to the sidewall (where the grid cells are largest) is insufficient. We compared the fully resolved thermal boundary layer profile with the Prandtl-Blasius profile. We find that the boundary layer profile is closer to the Prandtl Blasius profile at the cylinder axis than close to the sidewall, due to rising plumes in that region.

preprint2009arXiv

The effect of finite container size on granular jet formation

When an object is dropped into a bed of fine, loosely packed sand, a surprisingly energetic jet shoots out of the bed. In this work we study the effect that boundaries have on the granular jet formation. We did this by (i) decreasing the depth of the sand bed and (ii) reducing the container diameter to only a few ball diameters. These confinements change the behavior of the ball inside the bed, the void collapse, and the resulting jet height and shape. We map the parameter space of impact with Froude number, ambient pressure, and container dimensions as parameters. From these results we propose a new explanation for the thick-thin structure of the jet reported by several groups [J.R.Royer \textit{et al}., Nature Phys. \textbf{1}, 164 (2005)], [G.Caballero \textit{et al}., Phys. Rev. Lett. \textbf{99}, 018001 (2007)] and [J.O. Marston, \textit{et al}., Physics of Fluids \textbf{20}, 023301 (2008)].

preprint2009arXiv

Transitions between turbulent states in rotating Rayleigh-Benard convection

Weakly-rotating turbulent Rayleigh-Benard convection was studied experimentally and numerically. With increasing rotation and large enough Rayleigh number an abrupt transition from a turbulent state with nearly rotation-independent heat transport to another turbulent state with enhanced heat transfer is observed at a critical inverse Rossby number $1/Ro_c \simeq 0.4$. Whereas for $1/Ro < 1/Ro_c$ the strength of the large-scale convection-roll is either enhanced or essentially unmodified depending on parameters, its strength is increasingly diminished beyond $1/Ro_c$ where it competes with Ekman vortices that cause vertical fluid transport and thus heat-transfer enhancement.

preprint2008arXiv

Controlled impact of a disk on a water surface: Cavity dynamics

In this paper we study the transient surface cavity which is created by the controlled impact of a disk of radius h0 on a water surface at Froude numbers below 200. The dynamics of the transient free surface is recorded by high speed imaging and compared to boundary integral simulations. An excellent agreement is found between both. The flow surrounding the cavity is measured with high speed particle image velocimetry and is found to also agree perfectly with the flow field obtained from the simulations. We present a simple model for the radial dynamics of the cavity based on the collapse of an infinite cylinder. This model accounts for the observed asymmetry of the radial dynamics between the expansion and contraction phase of the cavity. It reproduces the scaling of the closure depth and total depth of the cavity which are both found to scale roughly proportional to Fr^{1/2} with a weakly Froude number dependent prefactor. In addition, the model accurately captures the dynamics of the minimal radius of the cavity, the scaling of the volume Vbubble of air entrained by the process, namely Vbubble/h0^3 proportional (1 + 0.26Fr^{1/2})Fr^{1/2}, and gives insight into the axial asymmetry of the pinch-off process.

preprint2008arXiv

Dimensionality and morphology of particle and bubble clusters in turbulent flow

We conduct numerical experiments to investigate the spatial clustering of particles and bubbles in simulations of homogeneous and isotropic turbulence. Varying the Stokes parameter and the densities, striking differences in the clustering of the particles can be observed. To quantify these visual findings we use the Kaplan--Yorke dimension. This local scaling analysis shows a dimension of approximately 1.4 for the light bubble distribution, whereas the distribution of very heavy particles shows a dimension of approximately 2.4. However, clearly separate parameter combinations yield the same dimensions. To overcome this degeneracy and to further develop the understanding of clustering, we perform a morphological (geometrical and topological) analysis of the particle distribution. For such an analysis, Minkowski functionals have been successfully employed in cosmology, in order to quantify the global geometry and topology of the large-scale distribution of galaxies. In the context of dispersed multiphase flow, these Minkowski functionals -- being morphological order parameters -- allow us to discern the filamentary structure of the light particle distribution from the wall-like distribution of heavy particles around empty interconnected tunnels.

preprint2008arXiv

Electrolytically Generated Nanobubbles on HOPG Surfaces

Electrolysis of water is employed to produce surface nanobubbles on highly orientated pyrolytic graphite (HOPG) surfaces. Hydrogen (oxygen) nanobubbles are formed when the HOPG surface acts as negative (positive) electrode. Coverage and volume of the nanobubbles enhance with increasing voltage. The yield of hydrogen nanobubbles is much larger than the yield of oxygen nanobubbles. The growth of the individual nanobubbles during the electrolysis process is recorded in time with the help of AFM measurements and correlated with the total current. Both the size of the individual nanobubbles and the total current saturate after typical 1 minute; then the nanobubbles are in a dynamic equilibrium, meaning that they do not further grow, in spite of ongoing gas production and nonzero current. The surface area of nanobubbles shows a good correlation with the nanobubble volume growth rate, suggesting that either the electrolytic gas emerges directly at the nanobubbles' surface, or it emerges at the electrode's surface and then diffuses through the nanobubbles' surface. Moreover, the experiments reveal that the time constants of the current and the aspect ratio of nanobubbles are the same under all conditions. Replacement of pure water by water containing a small amount of sodium chloride (0.01 M) allows for larger currents, but qualitatively gives the same results.

preprint2008arXiv

Flow organization in non-Oberbeck-Boussinesq Rayleigh-Benard convection in water

Non-Oberbeck-Boussinesq (NOB) effects on the flow organization in two-dimensional Rayleigh-Benard turbulence are numerically analyzed. The working fluid is water. We focus on the temperature profiles, the center temperature, the Nusselt number, and on the analysis of the velocity field. Several velocity amplitudes (or Reynolds numbers) and several kinetic profiles are introduced and studied; these together describe the various features of the rather complex flow organization. The results are presented both as functions of the Rayleigh number Ra (with Ra up to 10^8) for fixed temperature difference (Delta) between top and bottom plates and as functions of Delta ("non-Oberbeck-Boussinesqness") for fixed Ra with Delta up to 60 K. All results are consistent with the available experimental NOB data for the center temperature Tc and the Nusselt number ratio Nu_{NOB}/Nu_{OB} (the label OB meaning that the Oberbeck-Boussinesq conditions are valid). Beyond Ra ~ 10^6 the flow consists of a large diagonal center convection roll and two smaller rolls in the upper and lower corners. In the NOB case the center convection roll is still characterized by only one velocity scale.

preprint2008arXiv

Heat transfer mechanisms in bubbly Rayleigh-Benard convection

The heat transfer mechanism in Rayleigh-Benard convection in a liquid with a mean temperature close to its boiling point is studied through numerical simulations with point-like vapor bubbles, which are allowed to grow or shrink through evaporation and condensation and which act back on the flow both thermally and mechanically. It is shown that the effect of the bubbles is strongly dependent on the ratio of the sensible heat to the latent heat as embodied in the Jacob number Ja. For very small Ja the bubbles stabilize the flow by absorbing heat in the warmer regions and releasing it in the colder regions. With an increase in Ja, the added buoyancy due to the bubble growth destabilizes the flow with respect to single-phase convection and considerably increases the Nusselt number.

preprint2008arXiv

Large scale dynamics in turbulent Rayleigh-Benard convection

The progress in our understanding of several aspects of turbulent Rayleigh-Benard convection is reviewed. The focus is on the question of how the Nusselt number and the Reynolds number depend on the Rayleigh number Ra and the Prandtl number Pr, and on how the thicknesses of the thermal and the kinetic boundary layers scale with Ra and Pr. Non-Oberbeck-Boussinesq effects and the dynamics of the large-scale convection-roll are addressed as well. The review ends with a list of challenges for future research on the turbulent Rayleigh-Benard system.

preprint2008arXiv

Prandtl-, Rayleigh-, and Rossby-number dependence of heat transport in turbulent rotating Rayleigh-Benard convection

Experimental and numerical data for the heat transfer as a function of the Rayleigh-, Prandtl-, and Rossby numbers in turbulent rotating Rayleigh-Benard convection are presented. For relatively small $Ra ~ 10^8$ and large Pr modest rotation can enhance the heat transfer by up to 30%. At larger Ra there is less heat-transfer enhancement, and at small Pr = 0.7. there is no heat-transfer enhancement at all. We suggest that the small-Pr behavior is due to the breakdown of the heat-transfer-enhancing Ekman pumping because of larger thermal diffusion.

preprint2007arXiv

Microbubbly drag reduction in Taylor-Couette flow in the wavy vortex regime

We investigate the effect of microbubbles on Taylor-Couette flow by means of direct numerical simulations. We employ an Eulerian-Lagrangian approach with a gas-fluid coupling based on the point-force approximation. Added mass, drag, lift, and gravity are taken into account in the modeling of the motion of the individual bubble. We find that very dilute suspensions of small non-deformable bubbles (volume void fraction below 1%, zero Weber number and bubble Reynolds number <10) induce a robust statistically steady drag reduction (up to 20%) in the so called wavy vortex flow regime (Re = 600-2500). The Reynolds number dependence of the normalized torque (the so-called Torque Reduction Ratio (TRR) which corresponds to the drag reduction) is consistent with a recent series of experimental measurements performed by Murai et al. (J. Phys. 14, 143 (2005)). Our analysis suggests that the physical mechanism for the torque reduction in this regime is due to the local axial forcing, induced by rising bubbles, that is able to break the highly dissipative Taylor wavy vortices in the system. We finally show that the lift force acting on the bubble is crucial in this process. When neglecting it, the bubbles preferentially accumulate near the inner cylinder and the bulk flow is less efficiently modified.

preprint2006arXiv

Leaping shampoo and the stable Kaye effect

Shear-thinning fluids exhibit surprisingly rich behaviour. One example is the Kaye effect which occurs when a thin stream of a solution of polyisobutylene in Decalin is poured into a dish of the fluid. As pouring proceeds, a small stream of liquid occasionally leaps upward from the heap. This surprising effect, which lasts only a second or so, is named after its first observer A. Kaye, who could offer no explanation for this behaviour. Later, Collyer and Fischer suggested from 250 frames per second cine recordings that the fluid must be highly shear thinning as well as elastic and 'pituitous'. In addition, they concluded that a rigid surface is required to back the reflected liquid stream. While the words bouncing and reflection are associated with non-continuous and elastic effects, we will show here that the Kaye effect is in fact a continuous flow phenomenon. We show that the Kaye effect works for many common fluids, including shampoos and liquid soaps. We reveal its physical mechanism (formation, stability and disruption) through high-speed imaging. The measurements are interpreted with a simple theoretical model including only the shear thinning behaviour of the liquid; elastic properties of the liquid play no role. We show that the Kaye effect can be stable and that it can be directed. We even demonstrate a stable Kaye effect on a thin soap film excluding the necessity of a rigid backing surface.

preprint2004arXiv

Spontaneous Ratchet Effect in a Granular Gas

The spontaneous clustering of a vibrofluidized granular gas is employed to generate directed transport in two different compartmentalized systems: a "granular fountain" in which the transport takes the form of convection rolls, and a "granular ratchet" with a spontaneous particle current perpendicular to the direction of energy input. In both instances, transport is not due to any system-intrinsic anisotropy, but arises as a spontaneous collective symmetry breaking effect of many interacting granular particles. The experimental and numerical results are quantitatively accounted for within a flux model.

preprint1996arXiv

Acoustic Energy Storage in Single Bubble Sonoluminescence

Single bubble sonoluminescence is understood in terms of a shock focusing towards the bubble center. We present a mechanism for significantly enhancing the effect of shock focusing, arising from the storage of energy in the acoustic modes of the gas. The modes with strongest coupling are not spherically symmetric. The storage of acoustic energy gives a framework for understanding how light intensities depend so strongly on ambient gases and liquids and suggests that the light intensities of successive flashes are highly correlated.

preprint1996arXiv

Scaling and Linear Response in the GOY Turbulence model

The GOY model is a model for turbulence in which two conserved quantities cascade up and down a linear array of shells. When the viscosity parameter, $ν$, is small the model has a qualitative behavior which is similar to the Kolmogorov theories of turbulence. Here a static solution to the model is examined, and a linear stability analysis is performed to obtain response eigenvalues and eigenfunctions. Both the static behavior and the linear response show an inertial range with a relatively simple scaling structure. Our main results are: (i) The response frequencies cover a wide range of scales, with ratios which can be understood in terms of the frequency scaling properties of the model. (ii) Even small viscosities play a crucial role in determining the model's eigenvalue spectrum. (iii) As a parameter within the model is varied, it shows a ``phase transition'' in which there is an abrupt change in many eigenvalues from stable to unstable values. (iv) The abrupt change is determined by the model's conservation laws and symmetries. This work is thus intended to add to our knowledge of the linear response of a stiff dynamical systems and at the same time to help illuminate scaling within a class of turbulence models.

preprint1995arXiv

Bubble Shape Oscillations and the Onset of Sonoluminescence

An air bubble trapped in water by an oscillating acoustic field undergoes either radial or nonspherical pulsations depending on the strength of the forcing pressure. Two different instability mechanisms (the Rayleigh--Taylor instability and parametric instability) cause deviations from sphericity. Distinguishing these mechanisms allows explanation of many features of recent experiments on sonoluminescence, and suggests methods for finding sonoluminescence in different parameter regimes.

preprint1995arXiv

How the viscous subrange determines inertial range properties in turbulence shell models

We calculate static solutions of the 'GOY' shell model of turbulence and do a linear stability analysis. The asymptotic limit of large Reynolds numbers is analyzed. A phase diagram is presented which shows the range of stability of the static solution. We see an unexpected oscillatory dependence of the stability range upon $\lg ν$, where $ν$ is the viscosity. This effect depends upon the discrete structure of the shell model and goes to zero as the separation between the shells is brought to zero. These findings show how viscous effects play a role in determining inertial properties of shell models and give some hints for understanding the effects of viscous dissipation upon real turbulence.

Detlef Lohse

What is connected

Connect this record

See the researcher in context

Building this map preview

210 published item(s)

Settling of chiral particles in a turbulent flow

Transition from classical to ultimate melting

Collective effects of neighbouring melting ice objects

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

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

Double-diffusive transport in multicomponent vertical convection

Elastocapillary Worthington jets

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

Droplet dissolution driven by emerging thermal gradients and Marangoni flow

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

Entrapment of Interfacial Nanobubbles on Nano Structured Surfaces

Impact forces of water drops falling on superhydrophobic surfaces

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

Oscillatory droplet dissolution from competing Marangoni and gravitational flows

Physiochemical hydrodynamics of the phase segregation in an evaporating binary microdroplet

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

Taylor-Culick retractions and the influence of the surroundings

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

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

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

An efficient phase-field method for turbulent multiphase flows

Emergence of bimodal motility in active droplets

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

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

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

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

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

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

On explosive boiling of a multicomponent Leidenfrost drop

Plasmonic Microbubble Dynamics in Binary Liquids

Strong alignment of prolate ellipsoids in Taylor-Couette flow

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

Bubble puzzles: From fundamentals to applications

Bubbly and Buoyant Particle-Laden Turbulent Flows

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

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

Competing Marangoni and Rayleigh convection in evaporating binary droplets

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

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

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

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

Entrapment and Dissolution of Microbubbles Inside Microwells

Evaporation-Induced Crystallization of Surfactants in Sessile Multicomponent Droplets

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

Flow organization in laterally unconfined Rayleigh-Bénard turbulence

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

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

Gas-Vapor Interplay in Plasmonic Bubble Shrinkage

Giant plasmonic bubbles nucleation under different ambient pressures

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

Initial solidification dynamics of spreading droplets

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

Morphology of Evaporating Sessile Microdroplets on Lyophilic Elliptical Patches

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

News & Views: The distance rule in times of Corona

Physicochemical Hydrodynamics of Droplets out of Equilibrium: A Perspective Review

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

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

Scaling in internally heated convection: a unifying theory

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

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

Speeding up biphasic reactions with surface nanodroplets

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

Subcritical behaviour in double diffusive convection within the diffusive regime

The effect of Prandtl number on turbulent sheared thermal convection

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

Towards DNS of the Ultimate Regime of Rayleigh--Bénard Convection

Two-layer Thermally Driven Turbulence: Mechanisms for Interface Breakup

Ultrasound-enhanced mass transfer during single-bubble diffusive growth

Water entry of spheres into a rotating liquid

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

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

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

Convection-dominated dissolution for single and multiple immersed sessile droplets