Source author record

Chao Sun

Chao Sun 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 astro-ph.IM astro-ph.GA eess.SP hep-ph nlin.CD nlin.PS Operating Systems Performance physics.chem-ph physics.optics

Catalog footprint

What is connected

78works

12topics

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.

preprint2025arXiv

Freezing and ice aging dynamics in saline water under natural convection

Understanding the coupled dynamics of liquid-solid phase change and fluid flows is crucial in a wide range of geophysical and industrial applications. When freezing occurs in saline water, the newly formed ice is mushy, with a porous structure that traps the brine within the ice. In this work, which combines experiments and theoretical analyses, we investigate the long-term evolution of saline ice, comprehensively accounting for the coupled dynamics of multiscale fluid flow, heat and mass transfer, and phase change. We show that in a closed convective system the rapid formation of a mushy ice layer is followed by desalination (i.e, the expulsion of salt from the ice) processes that might lead to a slow asymptotic decrease of the ice thickness. Desalination of mushy ice reduces its porosity, which alters the dynamic thermal equilibrium and ice thickness by weakening buoyancy-driven convection within the mushy layer. In turn, changes in brine convection and ice thickness affect the further desalination of the ice. The long-term dynamics of the system can be accurately predicted by a one-dimensional model based on appropriate parameterizations of global heat and mass transfer properties. Furthermore, within the same theoretical model we explore the ice dynamics across a broader parameter space. Our findings advance the understanding of the coupled phase-change physics of saline solutions in the presence of convective fluid flows and provide a basis for explaining and predicting real-world phenomena such as the aging of sea ice.

preprint2025arXiv

Gravity drives the flow within the Stewartson layer in centrifugal convection

We conduct three-dimensional numerical simulations on centrifugal convection (CC) in a closed annular container, incorporating gravity and no-slip top and bottom boundaries, to systematically investigate rotation-induced secondary flow. The Stewartson layer, identified by an elongated circulation in mean vertical velocity plots, emerges near the inner and outer cylinders only beyond a critical gravitational forcing. Quantitative analyses confirm that the layer thickness scales as $δ_{st}\sim Ek^{1/3}$ due to rotational effects, consistent with results from rotating Rayleigh-Bénard convection, where $Ek$ represents the Ekman number. The internal circulation strength, however, is determined by both gravitational buoyancy and rotational effects. We propose that gravitational buoyancy drives the internal flow, which balances against viscous forces to establish a terminal velocity. Through theoretical analysis, the vertical velocity amplitude follows $W_{st}\sim Ek^{5/3}Ro^{-1}Ra_gPr^{-1}$, showing good agreement with simulation results across a wide parameter range. Here, $Ro^{-1}$ represents the inverse Rossby number, $Ra_g$ the gravitational Rayleigh number, and $Pr$ the Prandtl number. The theoretical predictions match simulations well, demonstrating that the Stewartson layer is gravity-induced and rotationally constrained through geostrophic balance in the CC system. These findings yield fundamental insights into turbulent flow structures and heat transfer mechanisms in the CC system, offering both theoretical advances and practical engineering applications.

preprint2025arXiv

Marangoni-driven freezing dynamics of supercooled binary droplets

Solidification of droplets is of great importance to various technological applications, drawing considerable attention from scientists aiming to unravel the fundamental physical mechanisms. In the case of multicomponent droplets undergoing solidification, the emergence of concentration gradients may trigger significant interfacial flows that dominate the freezing dynamics. Here, we experimentally investigate the fascinating interfacial freezing dynamics of supercooled ethanol-water droplets, accompanied with the migration and growth of massive ice particles. We reveal that these unique freezing dynamics are driven by solidification-induced solutal Marangoni flow within the droplets. Our model, which incorporates the temperature- and concentration-dependent properties of the ethanol-water mixture, quantitatively predicts both the migration velocity and the growth rate of the ice particles. The former is determined by the solutal Marangoni flow velocity, while the latter is governed by a balance between the latent heat release and the enhanced thermal dissipation by the Marangoni flow. Moreover, we show that the final wrapping state of droplets can be modulated by the concentration of ethanol. Our findings may pave the way for novel insights into the physicochemical hydrodynamics of multicomponent liquids undergoing phase transitions.

preprint2024arXiv

A Closed-loop Brain-Machine Interface SoC Featuring a 0.2$μ$J/class Multiplexer Based Neural Network

This work presents the first fabricated electrophysiology-optogenetic closed-loop bidirectional brain-machine interface (CL-BBMI) system-on-chip (SoC) with electrical neural signal recording, on-chip sleep staging and optogenetic stimulation. The first multiplexer with static assignment based table lookup solution (MUXnet) for multiplier-free NN processor was proposed. A state-of-the-art average accuracy of 82.4% was achieved with an energy consumption of only 0.2$μ$J/class in sleep staging task.

preprint2023arXiv

HLC2: a highly efficient cross-matching framework for large astronomical catalogues on heterogeneous computing environments

Cross-matching operation, which is to find corresponding data for the same celestial object or region from multiple catalogues,is indispensable to astronomical data analysis and research. Due to the large amount of astronomical catalogues generated by the ongoing and next-generation large-scale sky surveys, the time complexity of the cross-matching is increasing dramatically. Heterogeneous computing environments provide a theoretical possibility to accelerate the cross-matching, but the performance advantages of heterogeneous computing resources have not been fully utilized. To meet the challenge of cross-matching for substantial increasing amount of astronomical observation data, this paper proposes Heterogeneous-computing-enabled Large Catalogue Cross-matcher (HLC2), a high-performance cross-matching framework based on spherical position deviation on CPU-GPU heterogeneous computing platforms. It supports scalable and flexible cross-matching and can be directly applied to the fusion of large astronomical cataloguesfrom survey missions and astronomical data centres. A performance estimation model is proposed to locate the performance bottlenecks and guide the optimizations. A two-level partitioning strategy is designed to generate an optimized data placement according to the positions of celestial objects to increase throughput. To make HLC2 a more adaptive solution, the architecture-aware task splitting, thread parallelization, and concurrent scheduling strategies are designed and integrated. Moreover, a novel quad-direction strategy is proposed for the boundary problem to effectively balance performance and completeness. We have experimentally evaluated HLC2 using public released catalogue data. Experiments demonstrate that HLC2 scales well on different sizes of catalogues and the cross-matching speed is significantly improved compared to the state-of-the-art cross-matchers.

preprint2022arXiv

Accumulation and alignment of elongated gyrotactic swimmers in turbulence

We study the dynamics of gyrotactic swimmers in turbulence, whose orientation is governed by gravitational torque and local fluid velocity gradient. The gyrotaxis strength is measured by the ratio of the Kolmogorov time scale to the reorientation time scale due to gravity, and a large value of this ratio means the gyrotaxis is strong. By means of direct numerical simulations, we investigate the effects of swimming velocity and gyrotactic stability on spatial accumulation and alignment. Three-dimensional Vorono{\"ı} analysis is used to study the spatial distribution and time evolution of the particle concentration. We study spatial distribution by examing the overall preferential sampling and where clusters and voids (subsets of particles that have small and large Vorono{\"ı} volumes respectively) form. Compared with the ensemble particles, the preferential sampling of clusters and voids is found to be more pronounced. The clustering of fast swimmers lasts much longer than slower swimmers when the gyrotaxis is strong and intermediate, but an opposite trend emerges when the gyrotaxis is weak. In addition, we study the preferential alignment with the Lagrangian stretching direction, with which passive slender rods have been known to align. We show that the Lagrangian alignment is reduced by the swimming velocity when the gyrotaxis is weak, while the Lagrangian alignment is enhanced for the regime in which gyrotaxis is strong.

preprint2022arXiv

Characterization of Coupled Turbulent Wind-wave Flows Based on Large Eddy Simulation

Wind-wave interaction involves wind forcing on wave surface and wave effects on the turbulent wind structures, which essentially influences the wind and wave loading on structures. Existing research on wind-wave interaction modeling ignores the inherent strong turbulences of wind. The present study aims to characterize the turbulent airflow over wave surfaces and wave dynamics under wind driving force. A high-fidelity two-phase model is developed to simulate highly turbulent wind-wave fields. Instead of using uniform wind, inherent wind turbulences are prescribed at the inlet boundary using the turbulent spot method. The developed model is validated by comparing the simulated wind-wave flow characteristics with experimental data. With the validated model, a numerical case study is conducted on a 10^2 m scale under extreme wind and wave conditions. The result shows that when inherent wind turbulences are considered, the resultant turbulence is strengthened and is the summation of the inherent turbulence and the wave-induced turbulence. In addition, the wave coherent velocities and shelter effect are enhanced because of the presence of wind inherent turbulence. The regions of intense turbulence depend on the relative speed between wind velocity and wave phase speed. Higher wind velocities induce greater turbulence intensities, which can be increased by up to 17%. The different relative speed between wind and wave can induce opposite positive-negative patterns of wave coherent velocities. The wave-coherent velocity is approximately proportional to the wind velocity, and the influenced region mainly depends on the wave heights.

preprint2022arXiv

Decrypting material performance by wide-field femtosecond interferometric imaging of energy carrier evolution

Energy carrier evolution is crucial for material performance. Ultrafast microscopy has been widely applied to visualize the spatiotemporal evolution of energy carriers. However, direct imaging of small amounts of energy carriers on nanoscale remains difficult due to extremely weak transient signals. Here we present a method for ultrasensitive and high-throughput imaging of energy carrier evolution in space and time. This method combines femtosecond pump-probe techniques with interferometric scattering microscopy (iSCAT), named Femto-iSCAT. The interferometric principle and unique spatially-modulated contrast enhancement increase the transient image contrast by >2 orders of magnitude and enable the exploration of new science. We address three important and challenging problems: transport of different energy carriers at various interfaces, heterogeneous hot electron distribution and relaxation in single plasmonic resonators, and distinct structure-dependent edge state dynamics of carriers and excitons in optoelectronic semiconductors. Femto-iSCAT holds great potential as a universal tool for ultrasensitive imaging of energy carrier evolution in space and time.

preprint2022arXiv

Dynamics of finite-size spheroids in turbulent flow: the roles of flow structures and particle boundary layers

We study the translational and rotational dynamics of neutrally-buoyant finite-size spheroids in hydrodynamic turbulence by means of fully resolved numerical simulations. We examine axisymmetric shapes, from oblate to prolate, and the particle volume dependences. We show that the accelerations and rotations experienced by non-spherical inertial-scale particles result from volume filtered fluid forces and torques, similar to spherical particles. However, the particle orientations carry signatures of preferential alignments with the surrounding flow structures, which is reflected in distinct axial and lateral fluctuations for accelerations and rotation rates. The randomization of orientations does not occur even for particles with volume equivalent diameter size in the inertial range, here up to 60 $η$ at $Re_λ=120$. Additionally, we demonstrate that the role of fluid boundary layers around the particles cannot be neglected to reach a quantitative understanding of particle statistical dynamics, as they affect the intensities of angular velocities, and the relative importance of tumbling with respect to spinning rotations. This study brings to the fore the importance of inertial-scale flow structures in homogeneous and isotropic turbulence and their impacts on the transport of neutrally-buoyant bodies with size in the inertial range.

preprint2022arXiv

LeaFTL: A Learning-Based Flash Translation Layer for Solid-State Drives

In modern solid-state drives (SSDs), the indexing of flash pages is a critical component in their storage controllers. It not only affects the data access performance, but also determines the efficiency of the precious in-device DRAM resource. A variety of address mapping schemes and optimization techniques have been proposed. However, most of them were developed with human-driven heuristics. They cannot automatically capture diverse data access patterns at runtime in SSD controllers, which leaves a large room for improvement. In this paper, we present a learning-based flash translation layer (FTL), named LeaFTL, which learns the address mapping to tolerate dynamic data access patterns via linear regression at runtime. By grouping a large set of mapping entries into a learned segment, it significantly reduces the memory footprint of the address mapping table, which further benefits the data caching in SSD controllers. LeaFTL also employs various optimization techniques, including out-of-band metadata verification to tolerate mispredictions, optimized flash allocation, and dynamic compaction of learned index segments. We implement LeaFTL with an SSD simulator and evaluate it with various storage workloads. LeaFTL saves the memory consumption of the mapping table by 2.9x on average and improves the storage performance by 1.4x on average, in comparison with state-of-the-art FTL schemes.

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

Recent developments of turbulent emulsions in Taylor-Couette flow

Emulsions are common in many natural and industrial settings. Recently, much attention has been put on understanding the dynamics of turbulent emulsions. This paper reviews some recent studies of emulsions in turbulent Taylor-Couette flow, mainly focusing on the statistics of the dispersed phase and the global momentum transport of the system. We first study the size distribution and the breakup mechanism of the dispersed droplets for turbulent emulsions with a low volume-fraction (dilute) of the dispersed phase. For systems with a high volume-fraction of the dispersed phase (dense), we address the detailed response of the global transport (effective viscosity) of the turbulent emulsion and its connection to the droplet statistics. Finally, we will discuss catastrophic phase inversions, which can happen when the volume fraction of the dispersed phase exceeds a critical value during dynamic emulsification. We end the manuscript with a summary and an outlook including some open questions for future research. This article is part of the theme issue `Taylor-Couette and Related Flows on the Centennial of Taylor's Seminal Philosophical Transactions Paper'.

preprint2022arXiv

Supergravitational turbulent thermal convection

High-Rayleigh number convective turbulence is ubiquitous in many natural phenomena and in industries, such as atmospheric circulations, oceanic flows, flows in the fluid core of planets, and energy generations. In this work, we present a novel approach to boost the Rayleigh number in thermal convection by exploiting centrifugal acceleration and rapidly rotating a cylindrical annulus to reach an effective gravity of 60 times Earth's gravity. We show that in the regime where the Coriolis effect is strong, the scaling exponent of Nusselt number versus Rayleigh number exceeds one-third once the Rayleigh number is large enough. The convective rolls revolve in prograde direction, signifying the emergence of zonal flow. The present findings open a new avenue on the exploration of high-Rayleigh number turbulent thermal convection and will improve the understanding of the flow dynamics and heat transfer processes in geophysical and astrophysical flows and other strongly rotating systems.

preprint2022arXiv

The Decay constants of $B_c(nS)$ and $B^*_c(nS)$

The decay constants of the low lying S-wave $B_c$ mesons, i.e. $B_c(nS)$ and $B^*_c(nS)$ with $n\leq 3$, are calculated in the nonrelativistic quark model. The running coupling of the strong interaction is taken into account, and the uncertainties due to varying parameters and losing Lorentz covariance are considered carefully. As a byproduct, the decay constants of the low lying S-wave charmonium and bottomium states are given in the appendixes.

preprint2021arXiv

How do the finite-size particles modify the drag in Taylor-Couette turbulent flow

We experimentally investigate the drag modification by neutrally buoyant finite-size particles with various aspect ratios in a Taylor-Couette (TC) turbulent flow. The current Reynolds number, $Re$, ranges from $6.5\times10^3$ to $2.6\times10^4$, and the particle volume fraction, $Φ$, is up to $10\%$. Particles with three kinds of aspect ratio, $λ$, are used: $λ=1/3$ (oblate), $λ=1$ (spherical) and $λ=3$ (prolate). Unlike the case of bubbly TC flow, we find that the suspended finite-size particles increase the drag of the TC system regardless of their aspect ratios. In addition, the normalized friction coefficient, $c_{f,Φ}/c_{f,Φ=0}$, decreases with increasing $Re$, the reason could be that in the current low volume fractions the turbulent stress becomes dominant at higher $Re$. As $Re$ increases, the particles distribute more evenly in the entire system, which results from both the greater turbulence intensity and the more pronounced finite-size effects of the particles at higher $Re$. Moreover, it is found that the variation of the particle aspect ratios leads to different particle collective effects. The suspended spherical particles, which tend to cluster near the walls and form a particle layer, significantly affect the boundary layer and result in maximum drag modification. The minimal drag modification is found in the oblate case, where the particles preferentially cluster in the bulk region, and thus the particle layer is absent. Based on the optical measurement results, it can be concluded that, in the low volume fraction ranges ($Φ=0.5\%$ and $Φ= 2\%$ here), the larger drag modification is connected to the near-wall particle clustering. The present findings suggest that the particle shape plays a significant role in drag modification.

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.

preprint2020arXiv

A Redistribution Tool for Long-Term Archive of Astronomical Observation Data

Astronomical observation data require long-term preservation, and the rapid accumulation of observation data makes it necessary to consider the cost of long-term archive storage. In addition to low-speed disk-based online storage, optical disk or tape-based offline storage can be used to save costs. However, for astronomical research that requires historical data (particularly time-domain astronomy), the performance and energy consumption of data-accessing techniques cause problems because the requested data (which are organized according to observation time) may be located across multiple storage devices. In this study, we design and develop a tool referred to as AstroLayout to redistribute the observation data using spatial aggregation. The core algorithm uses graph partitioning to generate an optimized data placement according to the original observation data statistics and the target storage system. For the given observation data, AstroLayout can copy the long-term archive in the target storage system in accordance with this placement. An efficiency evaluation shows that AstroLayout can reduce the number of devices activated when responding to data-access requests in time-domain astronomy research. In addition to improving the performance of data-accessing techniques, AstroLayout can also reduce the storage systems power consumption. For enhanced adaptability, it supports storage systems of any media, including optical disks, tapes, and hard disks.

preprint2020arXiv

Anisotropic particles in two-dimensional convective turbulence

The orientational dynamics of inertialess anisotropic particles transported by two-dimensional convective turbulent flows display a coexistence of regular and chaotic features. We numerically demonstrate that very elongated particles (rods) align preferentially with the direction of the fluid flow, i.e., horizontally close to the isothermal walls and dominantly vertically in the bulk. This behaviour is due to the the presence of a persistent large scale circulation flow structure, which induces strong shear at wall boundaries and in up/down-welling regions. The near-wall horizontal alignment of rods persists at increasing the Rayleigh number, while the vertical orientation in the bulk is progressively weakened by the corresponding increase of turbulence intensity. Furthermore, we show that very elongated particles are nearly orthogonal to the orientation of the temperature gradient, an alignment independent of the system dimensionality and which becomes exact only in the limit of infinite Prandtl number. Tumbling rates are extremely vigorous adjacent to the walls, where particles roughly perform Jeffery orbits. This implies that the root-mean-square near-wall tumbling rates for spheres are much stronger than for rods, up to $\mathcal{O}(10)$ times at $Ra\simeq 10^9$. In the turbulent bulk the situation reverses and rods tumble slightly faster than isotropic particles, in agreement with earlier observations in two-dimensional turbulence.

preprint2020arXiv

AstroCatR: a Mechanism and Tool for Efficient Time Series Reconstruction of Large-Scale Astronomical Catalogues

Time series data of celestial objects are commonly used to study valuable and unexpected objects such as extrasolar planets and supernova in time domain astronomy. Due to the rapid growth of data volume, traditional manual methods are becoming extremely hard and infeasible for continuously analyzing accumulated observation data. To meet such demands, we designed and implemented a special tool named AstroCatR that can efficiently and flexibly reconstruct time series data from large-scale astronomical catalogues. AstroCatR can load original catalogue data from Flexible Image Transport System (FITS) files or databases, match each item to determine which object it belongs to, and finally produce time series datasets. To support the high-performance parallel processing of large-scale datasets, AstroCatR uses the extract-transform-load (ETL) preprocessing module to create sky zone files and balance the workload. The matching module uses the overlapped indexing method and an in-memory reference table to improve accuracy and performance. The output of AstroCatR can be stored in CSV files or be transformed other into formats as needed. Simultaneously, the module-based software architecture ensures the flexibility and scalability of AstroCatR. We evaluated AstroCatR with actual observation data from The three Antarctic Survey Telescopes (AST3). The experiments demonstrate that AstroCatR can efficiently and flexibly reconstruct all time series data by setting relevant parameters and configuration files. Furthermore, the tool is approximately 3X faster than methods using relational database management systems at matching massive catalogues.

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

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

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

Global and local statistics in turbulent emulsions

Turbulent emulsions are complex physical systems characterized by a strong and dynamical coupling between small-scale droplets and large-scale rheology. By using a specifically designed Taylor-Couette (TC) shear flow system, we are able to characterize the statistical properties of a turbulent emulsion made of oil droplets dispersed in an ethanol-water continuous solution, at the oil volume fraction up to 40%. We find that the dependence of the droplet size on the Reynolds number of the flow at the volume fraction of 1% can be well described by Hinze's criterion. The distribution of droplet sizes is found to follow a log-normal distribution, hinting at a fragmentation process as the possible mechanism dominating droplet formation. Additionally, the effective viscosity of the turbulent emulsion increases with the volume fraction of the dispersed oil phase, and decreases when the shear strength is increased. We find that the dependence of the effective viscosity on the shear rate can be described by the Herschel-Bulkley model, with a flow index monotonically decreasing with increasing the oil volume fraction. This finding indicates that the degree of shear thinning systematically increases with the volume fraction of the dispersed phase. The current findings have important implications for bridging the knowledge on turbulence and low-Reynolds-number emulsion flows to turbulent emulsion flows.

preprint2020arXiv

Leidenfrost drop impact on inclined superheated substrates

In real applications, drops always impact on solid walls with various inclinations. For the oblique impact of a Leidenfrost drop, which has a vapor layer under its bottom surface to prevent its direct contact with the superheated substrate, the drop can nearly frictionlessly slide along the substrate accompanied by the spreading and the retracting. To individually study these processes, we experimentally observe ethanol drops impact on superheated inclined substrates using high-speed imaging from two different views synchronously. We first study the dynamic Leidenfrost temperature, which mainly depends on the normal Weber number ${We}_\perp $. Then the substrate temperature is set to be high enough to study the Leidenfrost drop behaviors. During the spreading process, drops always keep uniform. And the maximum spreading factor $D_m/D_0$ follows a power-law dependence on the large normal Weber number ${We}_\perp $ as $D_m/D_0 = \sqrt{We_\perp /12+2}$ for $We_\perp \geq 30$. During the retracting process, drops with low impact velocities become non-uniform due to the gravity effect. For the sliding process, the residence time of all studied drops is nearly a constant, which is not affected by the inclination and $We$ number. The frictionless vapor layer results in the dimensionless sliding distance $L/D_0$ follows a power-law dependence on the parallel Weber number $We_\parallel$ as $L/D_0 \propto We_\parallel^{1/2}$. Without direct contact with the substrate, the behaviors of drops can be separately determined by ${We}_\perp $ and $We_\parallel$. When the impact velocity is too high, the drop fragments into many tiny droplets, which is called the splashing phenomenon. The critical splashing criterion is found to be $We_\perp ^*\simeq$ 120 or $K_\perp = We_\perp Re_\perp^{1/2} \simeq$ 5300 in the current parameter regime.

preprint2020arXiv

Rotation of anisotropic particles in Rayleigh-Bénard turbulence

Inertialess anisotropic particles in a Rayleigh-Bénard turbulent flow show maximal tumbling rates for weakly oblate shapes, in contrast with the universal behaviour observed in developed turbulence where the mean tumbling rate monotonically decreases with the particle aspect ratio. This is due to the concurrent effect of turbulent fluctuations and of a mean shear flow whose intensity, we show, is determined by the kinetic boundary layers. In Rayleigh-Bénard turbulence prolate particles align preferentially with the fluid velocity, while oblate ones orient with the temperature gradient. This analysis elucidates the link between particle angular dynamics and small-scale properties of convective turbulence and has implications for the wider class of sheared turbulent flows.

preprint2020arXiv

Rotational dynamics of bottom-heavy rods in turbulence from experiments and numerical simulations

We successfully perform the three-dimensional tracking in a turbulent fluid flow of small asymmetrical particles that are neutrally-buoyant and bottom-heavy, i.e., they have a non-homogeneous mass distribution along their symmetry axis. We experimentally show how a tiny mass inhomogeneity can affect the particle orientation along the preferred vertical direction and modify its tumbling rate. The experiment is complemented by a series of simulations based on realistic Navier-Stokes turbulence and on a point-like particle model that is capable to explore the full range of parameter space characterized by the gravitational torque stability number and by the particle aspect ratio. We propose a theoretical perturbative prediction valid in the high bottom-heaviness regime that agrees well with the observed preferential orientation and tumbling rate of the particles. We also show that the heavy-tail shape of the probability distribution function of the tumbling rate is weakly affected by the bottom-heaviness of the particles.

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.

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

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.

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

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

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

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

Translational and rotational dynamics of a large buoyant sphere in turbulence

We report experimental measurements of the translational and rotational dynamics of a large buoyant sphere in isotropic turbulence. We introduce an efficient method to simultaneously determine the position and (absolute) orientation of a spherical body from visual observation. The method employs a minimization algorithm to obtain the orientation from the 2D projection of a specific pattern drawn onto the surface of the sphere. This has the advantages that it does not require a database of reference images, is easily scalable using parallel processing, and enables accurate absolute orientation reference. Analysis of the sphere's translational dynamics reveals clear differences between the streamwise and transverse directions. The translational auto-correlations and PDFs provide evidence for periodicity in the particle's dynamics even under turbulent conditions. The angular autocorrelations show weak periodicity. The angular accelerations exhibit wide tails, however without a directional dependence.

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.

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

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

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

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.

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

Lagrangian single particle turbulent statistics through the Hilbert-Huang Transform

The Hilbert-Huang transform is applied to analyze single particle Lagrangian velocity data from numerical simulations of hydrodynamic turbulence. The velocity trajectory is described in terms of a set of intrinsic mode functions, C_{i}(t), and of their instantaneous frequency, ω_{i}(t). On the basis of this decomposition we define the ω-conditioned statistical moments of the C_{i} modes, named q-order Hilbert Spectra (HS). We show that such new quantities have enhanced scaling properties as compared to traditional Fourier transform- or correlation-based (Structure Functions) statistical indicators, thus providing better insights into the turbulent energy transfer process. We present a clear empirical evidence that the energy-like quantity, i.e. the second-order HS, displays a linear scaling in time in the inertial range, as expected from dimensional analysis and never observed before. We also measure high order moment scaling exponents in a direct way, without resorting the Extended Self Similarity (ESS) procedure. This leads to a new estimate of the Lagrangian structure functions exponents which are consistent with the multifractal prediction in the Lagrangian frame as proposed in [Biferale et al., Phys. Rev. Lett. vol. 93, 064502 (2004)].

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

Levitation of a drop over a moving surface

We investigate the levitation of a drop gently deposited onto the inner wall of a rotating hollow cylinder. For a sufficient velocity of the wall, the drop steadily levitates over a thin air film and reaches a stable angular position in the cylinder, where the drag and lift balance the weight of the drop. Interferometric measurement yields the three-dimensional (3D) air film thickness under the drop and reveals the asymmetry of the profile along the direction of the wall motion. A two-dimensional (2D) model is presented which explains the levitation mechanism, captures the main characteristics of the air film shape and predicts two asymptotic regimes for the film thickness $h_0$: For large drops $h_0\sim\ca^{2/3}κ_b^{-1}$, as in the Bretherton problem, where $\ca$ is the capillary number based on the air viscosity and $κ_b$ is the curvature at the bottom of the drop. For small drops $h_0\sim\ca^{4/5}(aκ_b)^{4/5}κ_b^{-1}$, where $a$ is the capillary length.

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

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

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

Applying Laser Doppler Anemometry inside a Taylor-Couette geometry - Using a ray-tracer to correct for curvature effects

In the present work it will be shown how the curvature of the outer cylinder affects Laser Doppler anemometry measurements inside a Taylor-Couette apparatus. The measurement position and the measured velocity are altered by curved surfaces. Conventional methods for curvature correction are not applicable to our setup, and it will be shown how a ray-tracer can be used to solve this complication. By using a ray-tracer the focal position can be calculated, and the velocity can be corrected. The results of the ray-tracer are verified by measuring an a priori known velocity field, and after applying refractive corrections good agreement with theoretical predictions are found. The methods described in this paper are applied to measure the azimuthal velocity profiles in high Reynolds number Taylor-Couette flow for the case of outer cylinder rotation.

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

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

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

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

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

Wall forces on a sphere in a rotating liquid-filled cylinder

We experimentally study the behavior of a particle slightly denser than the surrounding liquid in solid body rotating flow. Earlier work revealed that a heavy particle has an unstable equilibrium point in unbounded rotation flows. In the confinement of the rotational flow by a cylindrical wall a heavy sphere with density 1.05 g/cm$^3$ describes an orbital motion in our experiments. This is due to the effect of the wall near the sphere, i.e. a repulsive force ($F_w$). We model $F_w$ on the sphere as a function of the distance from the wall ($L$): $F_W \propto L^{-4}$ as proposed by Takemura and Magnaudet (2003). Remarkably, the path from the model including $F_w$ reproduce the experimentally measured trajectory. In addition during an orbital motion the particle does not spin around its axis, and we provide a possible explanation for this phenomenon.

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

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

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

The role of Stewartson and Ekman layers in turbulent rotating Rayleigh-Bénard convection

When the classical Rayleigh-Bénard (RB) system is rotated about its vertical axis roughly three regimes can be identified. In regime I (weak rotation) the large scale circulation (LSC) is the dominant feature of the flow. In regime II (moderate rotation) the LSC is replaced by vertically aligned vortices. Regime III (strong rotation) is characterized by suppression of the vertical velocity fluctuations. Using results from experiments and direct numerical simulations of RB convection for a cell with a diameter-to-height aspect ratio equal to one at $Ra \sim 10^8-10^9$ ($Pr=4-6$) and $0 \lesssim 1/Ro \lesssim 25$ we identified the characteristics of the azimuthal temperature profiles at the sidewall in the different regimes. In regime I the azimuthal wall temperature profile shows a cosine shape and a vertical temperature gradient due to plumes that travel with the LSC close to the sidewall. In regime II and III this cosine profile disappears, but the vertical wall temperature gradient is still observed. It turns out that the vertical wall temperature gradient in regimes II and III has a different origin than that observed in regime I. It is caused by boundary layer dynamics characteristic for rotating flows, which drives a secondary flow that transports hot fluid up the sidewall in the lower part of the container and cold fluid downwards along the sidewall in the top part.

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

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.

preprint2008arXiv

Oscillations of the large-scale circulation in turbulent Rayleigh-Bénard convection: the off-center mode and its relationship with the torsional mode

We report an experimental study of the large-scale circulation (LSC) in a turbulent Rayleigh-Bénard convection cell with aspect ratio unity. The temperature-extremum-extraction (TEE) method for obtaining the dynamic information of the LSC is presented. With this method, the azimuthal angular positions of the hot ascending and cold descending flows along the sidewall are identified from the measured instantaneous azimuthal temperature profile. The motion of the LSC is then decomposed into two different modes: the azimuthal mode and the translational or off-center mode. Comparing to the previous sinusoidal-fitting (SF) method, it is found that both methods give the same information about the azimuthal motion of the LSC, but the TEE method in addition can provide information about the off-center motion of the LSC, which is found to oscillate time-periodically around the cell's central vertical axis with an amplitude being nearly independent of the turbulent intensity. It is further found that the azimuthal angular positions of the hot ascending flow near the bottom plate and the cold descending flow near the top plate oscillate periodically out of phase by $π$, leading to the torsional mode of the LSC. These oscillations are then propagated vertically along the sidewall by the hot ascending and cold descending fluids. When they reach the mid-height plane, the azimuthal positions of the hottest and coldest fluids again oscillate out of phase by $π$. It is this out-of-phase horizontal positional oscillation of the hottest and coldest fluids at the same horizontal plane that produces the off-center oscillation of the LSC. A direct velocity measurement further confirms the existence of the bulk off-center mode of the flow field near cell center.

Chao Sun

What is connected

Connect this record

See the researcher in context

Building this map preview

78 published item(s)

Freezing and ice aging dynamics in saline water under natural convection

Gravity drives the flow within the Stewartson layer in centrifugal convection

Marangoni-driven freezing dynamics of supercooled binary droplets

A Closed-loop Brain-Machine Interface SoC Featuring a 0.2$μ$J/class Multiplexer Based Neural Network

HLC2: a highly efficient cross-matching framework for large astronomical catalogues on heterogeneous computing environments

Accumulation and alignment of elongated gyrotactic swimmers in turbulence

Characterization of Coupled Turbulent Wind-wave Flows Based on Large Eddy Simulation

Decrypting material performance by wide-field femtosecond interferometric imaging of energy carrier evolution

Dynamics of finite-size spheroids in turbulent flow: the roles of flow structures and particle boundary layers

LeaFTL: A Learning-Based Flash Translation Layer for Solid-State Drives

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

Recent developments of turbulent emulsions in Taylor-Couette flow

Supergravitational turbulent thermal convection

The Decay constants of $B_c(nS)$ and $B^*_c(nS)$

How do the finite-size particles modify the drag in Taylor-Couette turbulent flow

On explosive boiling of a multicomponent Leidenfrost drop

A Redistribution Tool for Long-Term Archive of Astronomical Observation Data

Anisotropic particles in two-dimensional convective turbulence

AstroCatR: a Mechanism and Tool for Efficient Time Series Reconstruction of Large-Scale Astronomical Catalogues

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

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

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

Global and local statistics in turbulent emulsions

Leidenfrost drop impact on inclined superheated substrates

Rotation of anisotropic particles in Rayleigh-Bénard turbulence

Rotational dynamics of bottom-heavy rods in turbulence from experiments and numerical simulations

Water entry of spheres into a rotating liquid

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

Drag reduction in boiling Taylor-Couette turbulence

Dynamic Leidenfrost effect: relevant time- and length-scales

Energy spectra in turbulent bubbly flows

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

Laser impact on a drop

Microbubbles and microparticles are not faithful tracers of turbulent acceleration

Multiple states in highly turbulent Taylor-Couette flow

On the spreading of impacting drops

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

Translational and rotational dynamics of a large buoyant sphere in turbulence

Tribonucleation of bubbles

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

Drop Shaping by Laser-Pulse Impact

Salinity transfer in bounded double diffusive convection

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

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

Growing bubbles in a slightly supersaturated liquid solution

Lagrangian single particle turbulent statistics through the Hilbert-Huang Transform

Leidenfrost temperature increase for impacting droplets on carbon-nanofiber surfaces

Levitation of a drop over a moving surface

Logarithmic boundary layers in highly turbulent Taylor-Couette flow

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

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

Velocity profiles in strongly turbulent Taylor-Couette flow

Applying Laser Doppler Anemometry inside a Taylor-Couette geometry - Using a ray-tracer to correct for curvature effects

Highly focused supersonic microjets

Highly focused supersonic microjets: numerical simulations

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

Lagrangian statistics of light particles in turbulence

Maximal air bubble entrainment at liquid drop impact

Microdroplet impact at very high velocity

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

Optimal Taylor-Couette turbulence

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

The clustering morphology of freely rising deformable bubbles

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

Wall forces on a sphere in a rotating liquid-filled cylinder

Angular momentum transport and turbulence in laboratory models of Keplerian flows

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

Drop impact on superheated surfaces

Flow reversals in thermally driven turbulence

The role of Stewartson and Ekman layers in turbulent rotating Rayleigh-Bénard convection

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