Researcher profile

P. K. Yeung

P. K. Yeung contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - UnverifiedVerification L1Unclaimed author
4works
0followers
4topics
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

4 published item(s)

preprint2020arXiv

A highly scalable particle tracking algorithm using partitioned global address space (PGAS) programming for extreme-scale turbulence simulations

A new parallel algorithm utilizing partitioned global address space (PGAS) programming model to achieve high scalability is reported for particle tracking in direct numerical simulations of turbulent flow. The work is motivated by the desire to obtain Lagrangian information necessary for the study of turbulent dispersion at the largest problem sizes feasible on current and next-generation multi-petaflop supercomputers. A large population of fluid particles is distributed among parallel processes dynamically, based on instantaneous particle positions such that all of the interpolation information needed for each particle is available either locally on its host process or neighboring processes holding adjacent sub-domains of the velocity field. With cubic splines as the preferred interpolation method, the new algorithm is designed to minimize the need for communication, by transferring between adjacent processes only those spline coefficients determined to be necessary for specific particles. This transfer is implemented very efficiently as a one-sided communication, using Co-Array Fortran (CAF) features which facilitate small data movements between different local partitions of a large global array. Detailed benchmarks are obtained on the Cray petascale supercomputer Blue Waters at the University of Illinois, Urbana-Champaign. For operations on the particles in a $8192^3$ simulation ($0.55$ trillion grid points) on $262,144$ Cray XE6 cores, the new algorithm is found to be orders of magnitude faster relative to a prior algorithm in which each particle is tracked by the same parallel process at all times. Improving support of PGAS models on major compilers suggests that this algorithm will be of wider applicability on most upcoming supercomputers.

preprint2020arXiv

Scaling exponents saturate in three-dimensional isotropic turbulence

From a database of direct numerical simulations of homogeneous and isotropic turbulence, generated in periodic boxes of various sizes, we extract the spherically symmetric part of moments of velocity increments and first verify the following (somewhat contested) results: the $4/5$-ths law holds in an intermediate range of scales and that the second order exponent over the same range of scales is {\it{anomalous}}, departing from the self-similar value of $2/3$ and approaching a constant of $0.72$ at high Reynolds numbers. We compare with some typical theories the dependence of longitudinal exponents as well as their derivatives with respect to the moment order $n$, and estimate the most probable value of the Hölder exponent. We demonstrate that the transverse scaling exponents saturate for large $n$, and trace this trend to the presence of large localized jumps in the signal. The saturation value of about $2$ at the highest Reynolds number suggests, when interpreted in the spirit of fractals, the presence of vortex sheets rather than more complex singularities. In general, the scaling concept in hydrodynamic turbulence appears to be more complex than even the multifractal description.

preprint2019arXiv

Extreme velocity gradients in turbulent flows

Fully turbulent flows are characterized by intermittent formation of very localized and intense velocity gradients. These gradients can be orders of magnitude larger than their typical value and lead to many unique properties of turbulence. Using direct numerical simulations of the Navier-Stokes equations with unprecedented small-scale resolution, we characterize such extreme events over a significant range of turbulence intensities, parameterized by the Taylor-scale Reynolds number ($R_λ$). Remarkably, we find the strongest velocity gradients to empirically scale as $τ_K^{-1} R_λ^β$, with $β\approx 0.775 \pm 0.025$, where $τ_K$ is the Kolmogorov time scale (with its inverse, $τ_K^{-1}$, being the {r.m.s.} of velocity gradient fluctuations). Additionally, we observe velocity increments across very small distances $r \le η$, where $η$ is the Kolmogorov length scale, to be as large as the {r.m.s.} of the velocity fluctuations. Both observations suggest that the smallest length scale in the flow behaves as $ηR_λ^{-α}$, with $α= β- \frac{1}{2}$, which is at odds with predictions from existing phenomenological theories. We find that extreme gradients are arranged in vortex tubes, such that strain conditioned on vorticity grows on average slower than vorticity, approximately as a power law with an exponent $γ< 1$, which weakly increases with $R_λ$. Using scaling arguments, we get $β=(2-γ)^{-1}$, which suggests that $β$ would also slowly increase with $R_λ$. We conjecture that approaching the limit of infinite $R_λ$, the flow is overall smooth, with intense velocity gradients over scale $ ηR_λ^{-1/2}$, corresponding to $β= 1$.

preprint2019arXiv

Fractal iso-level sets in high-Reynolds-number scalar turbulence

We study the fractal scaling of iso-levels sets of a passive scalar mixed by three-dimensional homogeneous and isotropic turbulence at high Reynolds numbers. The Schmidt number is unity. A fractal box-counting dimension $D_F$ can be obtained for iso-levels below about $3$ standard deviations of the scalar fluctuation on either side of its mean value. The dimension varies systematically with the iso-level, with a maximum of about $8/3$ for the iso-level at the mean; this maximum dimension also follows as an upper bound from the geometric measure theory. We interpret this result to mean that mixing in turbulence is always incomplete. A unique box-counting dimension for all iso-levels results when we consider the spatial support of the steep cliffs of the scalar conditioned on local strain; that unique dimension is about $4/3$.