Source author record

P. K. Yeung

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

physics.flu-dyn physics.comp-ph Computation nlin.CD physics.data-an

Catalog footprint

Research overview

6works

5topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

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

preprint2016arXiv

Refined similarity hypothesis using 3D local averages

The refined similarity hypotheses of Kolmogorov, regarded as an important ingredient of intermittent turbulence, has been tested in the past using one-dimensional data and plausible surrogates of energy dissipation. We employ data from direct numerical simulations, at the microscale Reynolds number $R_λ\sim 650$ , on a periodic box of $4096^3$ grid points to test the hypotheses using 3D averages. In particular, we study the small-scale properties of the stochastic variable $V = Δu(r)/(r ε_r)^{1/3}$ , where $Δu(r)$ is the longitudinal velocity increment and $ε_r$ is the dissipation rate averaged over a three-dimensional volume of linear size $r$ . We show that $V$ is universal in the inertial subrange. In the dissipation range, the statistics of $V$ are shown to depend solely on a local Reynolds number.

preprint2008arXiv

Universal intermittent properties of particle trajectories in highly turbulent flows

We present a collection of eight data sets, from state-of-the-art experiments and numerical simulations on turbulent velocity statistics along particle trajectories obtained in different flows with Reynolds numbers in the range $R_λ\in [120:740]$ . Lagrangian structure functions from all data sets are found to collapse onto each other on a wide range of time lags, revealing a universal statistics, and calling for a unified theoretical description. Parisi-Frisch Multifractal theory, suitable extended to the dissipative scales and to the Lagrangian domain, is found to capture intermittency of velocity statistics over the whole three decades of temporal scales here investigated.

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.