Researcher profile

Chong Shen Ng

Chong Shen Ng contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
6topics
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

7 published item(s)

preprint2024arXiv

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

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

preprint2022arXiv

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

This numerical study presents a simple but extremely effective way to considerably enhance heat transport in turbulent multiphase flows, namely by using oleophilic walls. As a model system, we pick the Rayleigh-Bénard setup, filled with an oil-water mixture. For oleophilic walls, e.g. using only $10\%$ volume fraction of oil in water, we observe a remarkable heat transport enhancement of more than $100\%$ as compared to the pure water case. In contrast, for oleophobic walls, the enhancement is then only about $20\%$ as compared to pure water. The physical explanation of the highly-efficient heat transport for oleophilic walls is that thermal plumes detach from the oil-rich boundary layer and are transported together with the oil phase. In the bulk, the oil-water interface prevents the plumes to mix with the turbulent water bulk. To confirm this physical picture, we show that the minimum amount of oil to achieve the maximum heat transport is set by the volume fraction of the thermal plumes. Our findings provide guidelines of how to optimize heat transport in thermal turbulence. Moreover, the physical insight of how coherent structures are coupled with one phase of a two-phase system has very general applicability for controlling transport properties in other turbulent multiphase flows.

preprint2021arXiv

An efficient phase-field method for turbulent multiphase flows

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

preprint2021arXiv

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

To quantify the fate of respiratory droplets under different ambient relative humidities, direct numerical simulations of a typical respiratory event are performed. We found that, because small droplets (with initial diameter of 10um) are swept by turbulent eddies in the expelled humid puff, their lifetime gets extended by a factor of more than 30 times as compared to what is suggested by the classical picture by William F. Wells, for 50% relative humidity. With increasing ambient relative humidity the extension of the lifetimes of the small droplets further increases and goes up to around 150 times for 90% relative humidity, implying more than two meters advection range of the respiratory droplets within one second. Employing Lagrangian statistics, we demonstrate that the turbulent humid respiratory puff engulfs the small droplets, leading to many orders of magnitude increase in their lifetimes, implying that they can be transported much further during the respiratory events than the large ones. Our findings provide the starting points for larger parameter studies and may be instructive for developing strategies on optimizing ventilation and indoor humidity control. Such strategies are key in mitigating the COVID-19 pandemic in the present autumn and upcoming winter.

preprint2021arXiv

Strong alignment of prolate ellipsoids in Taylor-Couette flow

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

preprint2020arXiv

Two-layer Thermally Driven Turbulence: Mechanisms for Interface Breakup

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

preprint2019arXiv

Convection-dominated dissolution for single and multiple immersed sessile droplets

We numerically investigate both single and multiple droplet dissolution with droplets consisting of lighter liquid dissolving in a denser host liquid. The significance of buoyancy is quantified by the Rayleigh number Ra which is the buoyancy force over the viscous damping force. In this study, Ra spans almost four decades from 0.1 to 400. We focus on how the mass flux, characterized by the Sherwood number Sh, and the flow morphologies depend on Ra. For single droplet dissolution, we first show the transition of the Sh(Ra) scaling from a constant value to $Sh\sim Ra^{1/4}$, which confirms the experimental results by Dietrich et al. (J. Fluid Mech., vol. 794, 2016, pp. 45--67). The two distinct regimes, namely the diffusively- and the convectively-dominated regime, exhibit different flow morphologies: when Ra>=10, a buoyant plume is clearly visible which contrasts sharply to the pure diffusion case at low Ra. For multiple droplet dissolution, the well-known shielding effect comes into play at low Ra so that the dissolution rate is slower as compared to the single droplet case. However, at high Ra, convection becomes more and more dominant so that a collective plume enhances the mass flux, and remarkably the multiple droplets dissolve faster than a single droplet. This has also been found in the experiments by Laghezza et al. (Soft Matter, vol. 12, 2016, pp. 5787--5796). We explain this enhancement by the formation of a single, larger plume rather than several individual plumes. Moreover, there is an optimal Ra at which the enhancement is maximized, because the single plume is narrower at larger Ra, which thus hinders the enhancement. Our findings demonstrate a new mechanism in collective droplet dissolution, which is the merging of the plumes, that leads to non-trivial phenomena, contrasting the shielding effect.