Researcher profile

Sai Peng

Sai Peng contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
1topics
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

5 published item(s)

preprint2022arXiv

Comparison of three numerical stabilization techniques of viscoelastic flows: vortex shedding behind a confined cylinder

In this study, the OpenFOAM platform, based on the finite volume method, is applied to investigate the two-dimensional viscoelastic flow past a circular cylinder. The FENE-P model, which considers the bounded elongation of polymer molecules, is chosen to describe the elastic constitutive relationship of the polymer solution. The maximum molecular chain lengths of L = 10, 50, 100, and 200 are considered, which describe the molecular conformation characteristics of the polymer solution. To improve the numerical instability of the viscoelastic flow simulation, three different methods, i.e., the traditional method (Td) with the addition of artificial viscosity, the logarithmic reconstruction method (Log), and the square root tensor method (Sqrt), are evaluated. The results show that the artificial viscosity has a little effect on the accuracy for the simulation with a small molecular chain length (L = 10). However, for long molecular chain lengths such as L = 100 and L = 200, the addition of artificial dissipation tends to overestimate the drag, which indicates that special caution is needed to incorporate the artificial dissipation in the simulation. Moreover, the logarithmic reconstruction method shows a strong grid-dependent characteristics, which may produce unphysical results.

preprint2022arXiv

Effect of wall slip on laminar flow past a circular cylinder

A numerical study of two-dimensional flow past a confined circular cylinder with slip wall is performed. A dimensionless number, Knudsen number ($Kn$) is used to describe the slip length of cylinder wall. The Reynolds number ($Re$) and Knudsen number ($Kn$) ranges considered are $Re = [1, 180]$ and $Kn = [0, \infty)$, respectively. Time-averaged flow separation angle ($\bar{θ_s}$), dimensionless recirculation length ($\bar{L_s}$) and the tangential velocity ($\bar{u_τ}$) distributed on the cylinder's wall, drag coefficient ($\bar{C_d}$) and drag reduction ($DR$) are investigated. The time-averaged tangential velocity distribution on the cylinder's wall fit well with the formula $\bar{u_τ} = [\fracα{1+βe^{-γ(π-θ)}}+δ]sin(θ) $, where the coefficients ($α$, $β$, $γ$, $δ$) are related with $Re$ and $Kn$. Several scaling-laws are found, $log(\bar{u_{τmax}})\sim{log(Re)}$ and $\bar{u_{τmax}}\sim{Kn}$ for low $Kn$, ($\bar{u_{τmax}}$ is the maximum tangential velocity on the cylinder's wall), $log(DR)\sim{log(Re)}$ ($Re\leq45$ and $Kn\leq0.1$), $log(DR)\sim{log(Kn)}$ ($Kn\leq0.05$). At low $Re$, $DR_v$ (the friction drag reduction) is the main source of $DR$. However, $DR_p$ (the differential pressure drag reduction) contributes the most to $DR$ at high $Re$ ($Re>\sim60$) and $Kn$ over a critical number. $DR_v$ is found almost independent to $Re$.

preprint2022arXiv

Subcritical insability of viscoelastic flow over a circular cylinder: A numerical study

In this paper, we discuss whether the instability of viscoelastic flow around a circular cylinder is subcritical or supercritical by numerical simulation. The Oldroyd-B model is selected to describe the viscoelastic constitutive relationship. The Log-conformation reformulation is employed to stabilize numerical simulation. The parameter ranges investigated are the Reynolds numbers ($Re$) spanning from 5 to 100 and the Weissenberg number ($Wi$) spanning from 0 to 10, with a fixed viscosity ratio of $β= 0.9$. Simulations are performed under two paths, i.e., 1) increasing $Re$ (or $Wi$) and 2) decreasing $Re$ (or $Wi$) slowly and gradually from one state to the next. The results show that the statistical solutions such as the time-averaged velocity obtained along the two paths are not identical over certain parameter range, which is around the transition point from the steady to unsteady flow. This loading path dependence behaviour indicates that the flow instability is subcritical.

preprint2022arXiv

Wake asymmetry weakening in viscoelastic fluids: Numerical discovery and mechanism exploration

Viscoelasticity weakens the asymmetry of laminar shedding flow behind a blunt body in a free domain. In the present study, this finding is confirmed by four unsteady viscoelastic flows with asymmetric flow configuration, i.e., flow over an inclined flat plate with various angles of incidence, flow over a rotating circular cylinder, flow over a circular cylinder with asymmetric slip boundary distribution, and flow over an inclined row of eight equally closely spaced circular cylinders (which can be considered as a single large blunt body) through direct numerical simulation combined with the Peterlin approximation of the finitely extensible nonlinear elastic (FENE-P) model. At high Weissenberg number, an arc shape region with high elastic stress, which is similar to shock wave, forms in the frontal area of the blunt body. This region acts as a stationary shield to separate the flow into different regions. Thus, the free stream resembles to pass this shield instead of the original blunt body. As this shield has symmetric feature, the wake flow restores symmetry.

preprint2022arXiv

Wall-induced translation of a rotating particle in viscoelastic fluid

Shear-thinning and viscoelasticity are two non-Newtonian fluid properties widely existing in biological fluids. In this study, we found that the translation motion of a rotating particle near a wall speed up firstly, and then slows down with enhancement of fluid viscoelasticity, which is different from the behavior reported in shear thinning fluid (Chen et al. J. Fluid Mech. 2021, 927). Our research is carried out by numerical simulation of Navier-Stokes equations combined with Oldroyd-B constitutive model. This work is expected to be helpful to understand the movement of a rotating sphere near a wall in complex fluids comprehensively.