An official website of the United States government
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
Stochastic modelling of the instantaneous velocity profile in rough-wall turbulent boundary layers
The statistical properties of uniform momentum zones (UMZs) are extracted from laboratory and field measurements in rough wall turbulent boundary layers to formulate a set of stochastic models for the simulation of instantaneous velocity profiles. A spatiotemporally resolved velocity dataset, covering a field of view of$$8 \times 9\,{\rm m}^2$$, was obtained in the atmospheric surface layer using super-large-scale particle image velocimetry (SLPIV), as part of the Grand-scale Atmospheric Imaging Apparatus (GAIA). Wind tunnel data from a previous study are included for comparison (Heiselet al.,J. Fluid Mech., vol. 887, 2020, R1). The probability density function of UMZ attributes such as their thickness, modal velocity and averaged vertical velocity are built at varying elevations and modelled using log-normal and Gaussian distributions. Inverse transform sampling of the distributions is used to generate synthetic step-like velocity profiles that are spatially and temporally uncorrelated. Results show that in the wide range of wall-normal distances and$$Re_\tau$$up to$$\sim O(10^6)$$investigated here, shear velocity scaling is manifested in the velocity jump across shear interfaces between adjacent UMZs, and attached eddy behaviour is observed in the linear proportionality between UMZ thickness and their wall normal location. These very same characteristics are recovered in the generated instantaneous profiles, using both fully stochastic and data-driven hybrid stochastic (DHS) models, which address, in different ways, the coupling between modal velocities and UMZ thickness. Our method provides a stochastic approach for generating an ensemble of instantaneous velocity profiles, consistent with the structural organisation of UMZs, where the ensemble reproduces the logarithmic mean velocity profile and recovers significant portions of the Reynolds stresses and, thus, of the streamwise and vertical velocity variability.
more »
« less
- PAR ID:
- 10531953
- Publisher / Repository:
- Cambridge University Press
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 979
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We explored the settling dynamics of vertically aligned particles in a quiescent, stratified two-layer fluid using particle tracking velocimetry. Glass spheres of$$d=4\,{\rm mm}$$diameter were released at frequencies of 4, 6 and 8 Hz near the free surface, traversing through an upper ethanol layer ($$H_1$$), whereHis height or layer thickess, varying from$$10d$$to$$40d$$and a lower oil layer. Results reveal pronounced lateral particle motion in the ethanol layer, attributed to a higher Galileo number ($$Ga = 976$$, ratio of buoyancy–gravity to viscous effects), compared with the less active behaviour in the oil layer ($$Ga = 16$$). The ensemble vertical velocity of particles exhibited a minimum just past the density interface, becoming more pronounced with increasing$$H_1$$, and suggesting that enhanced entrainment from ethanol to oil resulted in an additional buoyancy force. This produced distinct patterns of particle acceleration near the density interface, which were marked by significant deceleration, indicating substantial resistance to particle motion. An increased drag coefficient occurred for$$H_1/d = 40$$compared with a single particle settling in oil; drag reduced as the particle-release frequency ($$\,f_p$$) increased, likely due to enhanced particle interactions at closer proximity. Particle pair dispersions, lateral ($$R^2_L$$) and vertical ($$R^2_z$$), were modulated by$$H_1$$, initial separation$$r_0$$and$$f_p$$. The$$R^2_L$$dispersion displayed ballistic scaling initially, Taylor scaling for$$r_0 < H_1$$and Richardson scaling for$$r_0 > H_1$$. In contrast,$$R^2_z$$followed a$$R^2_z \sim t^{5.5}$$scaling under$$r_0 < H_1$$. Both$$R^2_L$$and$$R^2_z$$plateaued at a distance from the interface, depending on$$H_1$$and$$f_p$$.more » « less
-
The attached-eddy model (AEM) predicts that the mean streamwise velocity and streamwise velocity variance profiles follow a logarithmic shape, while the vertical velocity variance remains invariant with height in the overlap region of high Reynolds number wall-bounded turbulent flows. Moreover, the AEM coefficients are presumed to attain asymptotically constant values at very high Reynolds numbers. Here, the AEM predictions are examined using sonic anemometer measurements in the near-neutral atmospheric surface layer, with a focus on the logarithmic behaviour of the streamwise velocity variance. Utilizing an extensive 210-day dataset collected from a 62 m meteorological tower located in the Eastern Snake River Plain, Idaho, USA, the inertial sublayer is first identified by analysing the measured momentum flux and mean velocity profiles. The logarithmic behaviour of the streamwise velocity variance and the associated ‘$$-1$$’ scaling of the streamwise velocity energy spectra are then investigated. The findings indicate that the Townsend–Perry coefficient ($$A_1$$) is influenced by mild non-stationarity that manifests itself as a Reynolds number dependence. After excluding non-stationary runs, and requiring the bulk Reynolds number defined using the atmospheric boundary layer height to be larger than$$4 \times 10^{7}$$, the inferred$$A_1$$converges to values ranging between 1 and 1.25, consistent with laboratory experiments. Furthermore, nine benchmark cases selected through a restrictive quality control reveal a close relation between the ‘$$-1$$’ scaling in the streamwise velocity energy spectrum and the logarithmic behaviour of streamwise velocity variance. However, additional data are required to determine whether the plateau value of the pre-multiplied streamwise velocity energy spectrum is identical to$$A_1$$.more » « less
-
Not all the information in a turbulent field is relevant for understanding particular regions or variables in the flow. Here, we present a method for decomposing a source field into its informative$$\boldsymbol {\varPhi }_{I}(\boldsymbol {x},t)$$and residual$$\boldsymbol {\varPhi }_{R}(\boldsymbol {x},t)$$components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All the necessary information for physical understanding, reduced-order modelling and control of the target variable is contained in$$\boldsymbol {\varPhi }_{I}(\boldsymbol {x},t)$$, whereas$$\boldsymbol {\varPhi }_{R}(\boldsymbol {x},t)$$offers no substantial utility in these contexts. The decomposition is formulated as an optimisation problem that seeks to maximise the time-lagged mutual information of the informative component with the target variable while minimising the mutual information with the residual component. The method is applied to extract the informative and residual components of the velocity field in a turbulent channel flow, using the wall shear stress as the target variable. We demonstrate the utility of IND in three scenarios: (i) physical insight into the effect of the velocity fluctuations on the wall shear stress; (ii) prediction of the wall shear stress using velocities far from the wall; and (iii) development of control strategies for drag reduction in a turbulent channel flow using opposition control. In case (i), IND reveals that the informative velocity related to wall shear stress consists of wall-attached high- and low-velocity streaks, collocated with regions of vertical motions and weak spanwise velocity. This informative structure is embedded within a larger-scale streak–roll structure of residual velocity, which bears no information about the wall shear stress. In case (ii), the best-performing model for predicting wall shear stress is a convolutional neural network that uses the informative component of the velocity as input, while the residual velocity component provides no predictive capabilities. Finally, in case (iii), we demonstrate that the informative component of the wall-normal velocity is closely linked to the observability of the target variable and holds the essential information needed to develop successful control strategies.more » « less
-
Abstract For any subset$$Z \subseteq {\mathbb {Q}}$$, consider the set$$S_Z$$of subfields$$L\subseteq {\overline {\mathbb {Q}}}$$which contain a co-infinite subset$$C \subseteq L$$that is universally definable inLsuch that$$C \cap {\mathbb {Q}}=Z$$. Placing a natural topology on the set$${\operatorname {Sub}({\overline {\mathbb {Q}}})}$$of subfields of$${\overline {\mathbb {Q}}}$$, we show that ifZis not thin in$${\mathbb {Q}}$$, then$$S_Z$$is meager in$${\operatorname {Sub}({\overline {\mathbb {Q}}})}$$. Here,thinandmeagerboth mean “small”, in terms of arithmetic geometry and topology, respectively. For example, this implies that only a meager set of fieldsLhave the property that the ring of algebraic integers$$\mathcal {O}_L$$is universally definable inL. The main tools are Hilbert’s Irreducibility Theorem and a new normal form theorem for existential definitions. The normal form theorem, which may be of independent interest, says roughly that every$$\exists $$-definable subset of an algebraic extension of$${\mathbb Q}$$is a finite union of single points and projections of hypersurfaces defined by absolutely irreducible polynomials.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2023.999
