skip to main content

Title: Direct numerical simulations of turbulent pipe flow up to

Well-resolved direct numerical simulations (DNS) have been performed of the flow in a smooth circular pipe of radius$R$and axial length$10{\rm \pi} R$at friction Reynolds numbers up to$Re_\tau =5200$using the pseudo-spectral code OPENPIPEFLOW. Various turbulence statistics are documented and compared with other DNS and experimental data in pipes as well as channels. Small but distinct differences between various datasets are identified. The friction factor$\lambda$overshoots by$2\,\%$and undershoots by$0.6\,\%$the Prandtl friction law at low and high$Re$ranges, respectively. In addition,$\lambda$in our results is slightly higher than in Pirozzoliet al.(J. Fluid Mech., vol. 926, 2021, A28), but matches well the experiments in Furuichiet al.(Phys. Fluids, vol. 27, issue 9, 2015, 095108). The log-law indicator function, which is nearly indistinguishable between pipe and channel up to$y^+=250$, has not yet developed a plateau farther away from the wall in the pipes even for the$Re_\tau =5200$cases. The wall shear stress fluctuations and the inner peak of the axial turbulence intensity – which grow monotonically with$Re_\tau$– are lower in the pipe than in the channel, but the difference decreases with increasing$Re_\tau$. While the wall value is slightly lower in the channel than in the pipe at the same$Re_\tau$, the inner peak of the pressure fluctuation shows negligible differences between them. The Reynolds number scaling of all these quantities agrees with both the logarithmic and defect-power laws if the coefficients are properly chosen. The one-dimensional spectrum of the axial velocity fluctuation exhibits a$k^{-1}$dependence at an intermediate distance from the wall – also seen in the channel. In summary, these high-fidelity data enable us to provide better insights into the flow physics in the pipes as well as the similarity/difference among different types of wall turbulence.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Cambridge University Press
Date Published:
Journal Name:
Journal of Fluid Mechanics
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The equation for the fluid velocity gradient along a Lagrangian trajectory immediately follows from the Navier–Stokes equation. However, such an equation involves two terms that cannot be determined from the velocity gradient along the chosen Lagrangian path: the pressure Hessian and the viscous Laplacian. A recent model handles these unclosed terms using a multi-level version of the recent deformation of Gaussian fields (RDGF) closure (Johnson & Meneveau,Phys. Rev. Fluids, vol. 2 (7), 2017, 072601). This model is in remarkable agreement with direct numerical simulations (DNS) data and works for arbitrary Taylor Reynolds numbers$\textit {Re}_\lambda$. Inspired by this, we develop a Lagrangian model for passive scalar gradients in isotropic turbulence. The equation for passive scalar gradients also involves an unclosed term in the Lagrangian frame, namely the scalar gradient diffusion term, which we model using the RDGF approach. However, comparisons of the statistics obtained from this model with DNS data reveal substantial errors due to erroneously large fluctuations generated by the model. We address this defect by incorporating into the closure approximation information regarding the scalar gradient production along the local trajectory history of the particle. This modified model makes predictions for the scalar gradients, their production rates, and alignments with the strain-rate eigenvectors that are in very good agreement with DNS data. However, while the model yields valid predictions up to$\textit {Re}_\lambda \approx 500$, beyond this, the model breaks down.

    more » « less
  2. Drag for wall-bounded flows is directly related to the spatial flux of spanwise vorticity outward from the wall. In turbulent flows a key contribution to this wall-normal flux arises from nonlinear advection and stretching of vorticity, interpretable as a cascade. We study this process using numerical simulation data of turbulent channel flow at friction Reynolds number$Re_\tau =1000$. The net transfer from the wall of spanwise vorticity created by downstream pressure drop is due to two large opposing fluxes, one which is ‘down-gradient’ or outward from the wall, where most vorticity concentrates, and the other which is ‘up-gradient’ or toward the wall and acting against strong viscous diffusion in the near-wall region. We present evidence that the up-gradient/down-gradient transport occurs by a mechanism of correlated inflow/outflow and spanwise vortex stretching/contraction that was proposed by Lighthill. This mechanism is essentially Lagrangian, but we explicate its relation to the Eulerian anti-symmetric vorticity flux tensor. As evidence for the mechanism, we study (i) statistical correlations of the wall-normal velocity and of wall-normal flux of spanwise vorticity, (ii) vorticity flux cospectra identifying eddies involved in nonlinear vorticity transport in the two opposing directions and (iii) visualizations of coherent vortex structures which contribute to the transport. The ‘D-type’ vortices contributing to down-gradient transport in the log layer are found to be attached, hairpin-type vortices. However, the ‘U-type’ vortices contributing to up-gradient transport are detached, wall-parallel, pancake-shaped vortices with strong spanwise vorticity, as expected by Lighthill's mechanism. We discuss modifications to the attached eddy model and implications for turbulent drag reduction.

    more » « less
  3. Abstract

    We investigate a novel geometric Iwasawa theory for${\mathbf Z}_p$-extensions of function fields over a perfect fieldkof characteristic$p>0$by replacing the usual study ofp-torsion in class groups with the study ofp-torsion class groupschemes. That is, if$\cdots \to X_2 \to X_1 \to X_0$is the tower of curves overkassociated with a${\mathbf Z}_p$-extension of function fields totally ramified over a finite nonempty set of places, we investigate the growth of thep-torsion group scheme in the Jacobian of$X_n$as$n\rightarrow \infty $. By Dieudonné theory, this amounts to studying the first de Rham cohomology groups of$X_n$equipped with natural actions of Frobenius and of the Cartier operatorV. We formulate and test a number of conjectures which predict striking regularity in the$k[V]$-module structure of the space$M_n:=H^0(X_n, \Omega ^1_{X_n/k})$of global regular differential forms as$n\rightarrow \infty .$For example, for each tower in a basic class of${\mathbf Z}_p$-towers, we conjecture that the dimension of the kernel of$V^r$on$M_n$is given by$a_r p^{2n} + \lambda _r n + c_r(n)$for allnsufficiently large, where$a_r, \lambda _r$are rational constants and$c_r : {\mathbf Z}/m_r {\mathbf Z} \to {\mathbf Q}$is a periodic function, depending onrand the tower. To provide evidence for these conjectures, we collect extensive experimental data based on new and more efficient algorithms for working with differentials on${\mathbf Z}_p$-towers of curves, and we prove our conjectures in the case$p=2$and$r=1$.

    more » « less
  4. The effects of surfactants on a mechanically generated plunging breaker are studied experimentally in a laboratory wave tank. Waves are generated using a dispersively focused wave packet with a characteristic wavelength of$\lambda _0 = 1.18$m. Experiments are performed with two sets of surfactant solutions. In the first set, increasing amounts of the soluble surfactant Triton X-100 are mixed into the tank water, while in the second set filtered tap water is left undisturbed in the tank for wait times ranging from 15 min to 21 h. Increasing Triton X-100 concentrations and longer wait times lead to surfactant-induced changes in the dynamic properties of the free surface in the tank. It is found that low surface concentrations of surfactants can dramatically change the wave breaking process by changing the shape of the jet and breaking up the entrained air cavity at the time of jet impact. Direct numerical simulations (DNS) of plunging breakers with constant surface tension are used to show that there is significant compression of the free surface near the plunging jet tip and dilatation elsewhere. To explore the effect of this compression/dilatation, the surface tension isotherm is measured in all experimental cases. The effects of surfactants on the plunging jet are shown to be primarily controlled by the surface tension gradient ($\Delta \mathcal {E}$) while the ambient surface tension of the undisturbed wave tank ($\sigma _0$) plays a secondary role.

    more » « less
  5. Based on a generalized local Kolmogorov–Hill equation expressing the evolution of kinetic energy integrated over spheres of size$\ell$in the inertial range of fluid turbulence, we examine a possible definition of entropy and entropy generation for turbulence. Its measurement from direct numerical simulations in isotropic turbulence leads to confirmation of the validity of the fluctuation relation (FR) from non-equilibrium thermodynamics in the inertial range of turbulent flows. Specifically, the ratio of probability densities of forward and inverse cascade at scale$\ell$is shown to follow exponential behaviour with the entropy generation rate if the latter is defined by including an appropriately defined notion of ‘temperature of turbulence’ proportional to the kinetic energy at scale$\ell$.

    more » « less