BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

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

Network

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

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

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

Researcher profile

Chao Sun

Chao Sun contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dyncond-mat.softastro-ph.IMeess.SPhep-phOperating SystemsPerformancephysics.chem-phphysics.optics

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
32works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

32 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On explosive boiling of a multicomponent Leidenfrost drop

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bubbly and Buoyant Particle-Laden Turbulent Flows

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

We perform a numerical study of the heat transfer and flow structure of Rayleigh-Bénard (RB) convection in (in most cases regular) porous media, which are comprised of circular, solid obstacles located on a square lattice. This study is focused on the role of porosity $ϕ$ in the flow properties during the transition process from the traditional RB convection with $ϕ=1$ (so no obstacles included) to Darcy-type porous-media convection with $ϕ$ approaching 0. Simulations are carried out in a cell with unity aspect ratio, for the Rayleigh number $Ra$ from $10^5$ to $10^{10}$ and varying porosities $ϕ$, at a fixed Prandtl number $Pr=4.3$, and we restrict ourselves to the two dimensional case. For fixed $Ra$, the Nusselt number $Nu$ is found to vary non-monotonously as a function of $ϕ$; namely, with decreasing $ϕ$, it first increases, before it decreases for $ϕ$ approaching 0. The non-monotonous behaviour of $Nu(ϕ)$ originates from two competing effects of the porous structure on the heat transfer. On the one hand, the flow coherence is enhanced in the porous media, which is beneficial for the heat transfer. On the other hand, the convection is slowed down by the enhanced resistance due to the porous structure, leading to heat transfer reduction. For fixed $ϕ$, depending on $Ra$, two different heat transfer regimes are identified, with different effective power-law behaviours of $Nu$ vs $Ra$, namely, a steep one for low $Ra$ when viscosity dominates, and the standard classical one for large $Ra$. The scaling crossover occurs when the thermal boundary layer thickness and the pore scale are comparable. The influences of the porous structure on the temperature and velocity fluctuations, convective heat flux, and energy dissipation rates are analysed, further demonstrating the competing effects of the porous structure to enhance or reduce the heat transfer.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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&#39;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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Water entry of spheres into a rotating liquid

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Drag reduction in boiling Taylor-Couette turbulence

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

0 citations0 reviewsNo code linkedOpen access