Nonasymptotic behavior of developing turbulent pipe flow

1976 ◽  
Vol 54 (3) ◽  
pp. 268-278 ◽  
Author(s):  
J. K. Reichert ◽  
R. S. Azad

Detailed measurements of mean velocity U profiles, in the inlet 70 diameters of a pipe, show that the development of turbulent pipe flow is nonasymptotic. Experiments were done at seven Reynolds numbers in the range 56 000–15 3000. Contours of U and V fields are presented for two representative Reynolds numbers. A U component peak exceeding the fully developed values has been found to occur along the pipe centerline. The Reynolds number behavior of the peak position has been determined. Hot film measurements of the mean wall shear stresses in the inlet region also show a nonasymptotic development consistent with the mean velocity results.

Author(s):  
Sowjanya Vijiapurapu ◽  
Jie Cui

Fully developed turbulent pipe flow is investigated by large eddy simulations (LES). The three-dimensional, unsteady, incompressible, filtered continuity and Navier-Stokes equations in cylindrical coordinates are discretized by a finite difference method. The spatial derivatives are approximated by second order conservative schemes. This scheme eliminates the numerical generation or dissipation of energy. The pressure Poisson equation is solved by FFT method and time is advanced through a third order Runge-Kutta method. The commonly used subgrid scale (SGS) models — the Smagorinsky model and the dynamic model are implemented and simulations are performed for fully developed turbulent pipe flow at two different Reynolds numbers. The flow features in terms of mean velocity as well as higher order turbulence intensities and correlations are presented and compared to experimental and DNS data available in literature. Extensive comparisons are made for cases using different grid resolution, different streamwise domain dimension, different sub-grid scale model, and, at two different Reynolds number. For two Reynolds numbers (5,000 and 30,000) tested in this study, the fine mesh (64 × 96 × 64, circumferential × radial × longitudinal) produces better results than the coarse mesh (32 × 48 × 32), indicating the significance of the grid resolution, especially near the pipe surface. On the fine mesh for the two Reynolds numbers, the results exhibit a slight Reynolds number effect, indicating the mesh needs to be further refined at higher Reynolds number. Simulations were performed for two domain sizes, namely 6D and 12D, where D is the pipe diameter. When the streamwise grid resolution remains unchanged, the two simulations show negligible difference. This ensures that a 6D domain is adequate to include the largest eddies in a fully developed turbulent pipe flow at the current Reynolds number. When the fine mesh is used, the subgrid scale models (Smagorinsky and Dynamic) provide limited contribution to the total turbulent kinetic energy. Although the current results agree quite well with other published LES simulations, when compared with the Law of the wall, benchmark experiments and DNS results, the simulated mean velocity in the log region is higher than the experimental and DNS data. Overall, it was observed that the numerical methods work satisfactorily well for turbulent pipe flows at low and high Reynolds numbers, and, the method has capability to be used in the simulation of flows with practical interest.


2015 ◽  
Vol 81 (826) ◽  
pp. 15-00091-15-00091 ◽  
Author(s):  
Yuki WADA ◽  
Noriyuki FURUICHII ◽  
Yoshiya TERAO ◽  
Yoshiyuki TSUJI

1973 ◽  
Vol 60 (4) ◽  
pp. 665-687 ◽  
Author(s):  
K. S. Yajnik ◽  
M. V. Subbaiah

The effects of swirl on internal turbulent flows are studied by conducting experiments on turbulent pipe flow with variable initial swirl. This first part of the study is primarily concerned with similarity laws. The mean velocity profiles, both away from and close to the wall, are found to admit similarity representations at sufficiently large Reynolds numbers, provided that flow reversal does not take place near the entrance. While the wall law is not sensibly dependent on swirl, the velocity defect law in its extended form is sensitive to swirl. Further, a logarithmic skin-friction law is obtained in which only the additive coefficient depends on swirl. This coefficient is found to vary linearly with the swirl angle in the range of the present experiments.


2008 ◽  
Vol 608 ◽  
pp. 81-112 ◽  
Author(s):  
XIAOHUA WU ◽  
PARVIZ MOIN

Fully developed incompressible turbulent pipe flow at bulk-velocity- and pipe-diameter-based Reynolds number ReD=44000 was simulated with second-order finite-difference methods on 630 million grid points. The corresponding Kármán number R+, based on pipe radius R, is 1142, and the computational domain length is 15R. The computed mean flow statistics agree well with Princeton Superpipe data at ReD=41727 and at ReD=74000. Second-order turbulence statistics show good agreement with experimental data at ReD=38000. Near the wall the gradient of $\mbox{ln}\overline{u}_{z}^{+}$ with respect to ln(1−r)+ varies with radius except for a narrow region, 70 < (1−r)+ < 120, within which the gradient is approximately 0.149. The gradient of $\overline{u}_{z}^{+}$ with respect to ln{(1−r)++a+} at the present relatively low Reynolds number of ReD=44000 is not consistent with the proposition that the mean axial velocity $\overline{u}_{z}^{+}$ is logarithmic with respect to the sum of the wall distance (1−r)+ and an additive constant a+ within a mesolayer below 300 wall units. For the standard case of a+=0 within the narrow region from (1−r)+=50 to 90, the gradient of $\overline{u}_{z}^{+}$ with respect to ln{(1−r)++a+} is approximately 2.35. Computational results at the lower Reynolds number ReD=5300 also agree well with existing data. The gradient of $\overline{u}_{z}$ with respect to 1−r at ReD=44000 is approximately equal to that at ReD=5300 for the region of 1−r > 0.4. For 5300 < ReD < 44000, bulk-velocity-normalized mean velocity defect profiles from the present DNS and from previous experiments collapse within the same radial range of 1−r > 0.4. A rationale based on the curvature of mean velocity gradient profile is proposed to understand the perplexing existence of logarithmic mean velocity profile in very-low-Reynolds-number pipe flows. Beyond ReD=44000, axial turbulence intensity varies linearly with radius within the range of 0.15 < 1−r < 0.7. Flow visualizations and two-point correlations reveal large-scale structures with comparable near-wall azimuthal dimensions at ReD=44000 and 5300 when measured in wall units. When normalized in outer units, streamwise coherence and azimuthal dimension of the large-scale structures in the pipe core away from the wall are also comparable at these two Reynolds numbers.


1979 ◽  
Vol 44 (3) ◽  
pp. 700-710 ◽  
Author(s):  
Ivan Fořt ◽  
Hans-Otto Möckel ◽  
Jan Drbohlav ◽  
Miroslav Hrach

Profiles of the mean velocity have been analyzed in the stream streaking from the region of rotating standard six-blade disc turbine impeller. The profiles were obtained experimentally using a hot film thermoanemometer probe. The results of the analysis is the determination of the effect of relative size of the impeller and vessel and the kinematic viscosity of the charge on three parameters of the axial profile of the mean velocity in the examined stream. No significant change of the parameter of width of the examined stream and the momentum flux in the stream has been found in the range of parameters d/D ##m <0.25; 0.50> and the Reynolds number for mixing ReM ##m <2.90 . 101; 1 . 105>. However, a significant influence has been found of ReM (at negligible effect of d/D) on the size of the hypothetical source of motion - the radius of the tangential cylindrical jet - a. The proposed phenomenological model of the turbulent stream in region of turbine impeller has been found adequate for values of ReM exceeding 1.0 . 103.


Author(s):  
B.J McKeon ◽  
J.F Morrison

The streamwise velocity component in turbulent pipe flow is assessed to determine whether it exhibits asymptotic behaviour that is indicative of high Reynolds numbers. The asymptotic behaviour of both the mean velocity (in the form of the log law) and that of the second moment of the streamwise component of velocity in the outer and overlap regions is consistent with the development of spectral regions which indicate inertial scaling. It is shown that an ‘inertial sublayer’ in physical space may be considered as a spatial analogue of the inertial subrange in the velocity spectrum and such behaviour only appears for Reynolds numbers R + >5×10 3 , approximately, much higher than was generally thought.


2014 ◽  
Vol 749 ◽  
pp. 79-98 ◽  
Author(s):  
S. C. C. Bailey ◽  
M. Vallikivi ◽  
M. Hultmark ◽  
A. J. Smits

AbstractFive separate data sets on the mean velocity distributions in the Princeton University/ONR Superpipe are used to establish the best estimate for the value of von Kármán’s constant for the flow in a fully developed, hydraulically smooth pipe. The profiles were taken using Pitot tubes, conventional hot wires and nanoscale thermal anemometry probes. The value of the constant was found to vary significantly due to measurement uncertainties in the mean velocity, friction velocity and the wall distance, and the number of data points included in the analysis. The best estimate for the von Kármán constant in turbulent pipe flow is found to be $0.40 \pm 0.02$. A more precise estimate will require improved instrumentation.


2013 ◽  
Vol 734 ◽  
pp. 275-316 ◽  
Author(s):  
Rashad Moarref ◽  
Ati S. Sharma ◽  
Joel A. Tropp ◽  
Beverley J. McKeon

AbstractWe study the Reynolds-number scaling and the geometric self-similarity of a gain-based, low-rank approximation to turbulent channel flows, determined by the resolvent formulation of McKeon & Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382), in order to obtain a description of the streamwise turbulence intensity from direct consideration of the Navier–Stokes equations. Under this formulation, the velocity field is decomposed into propagating waves (with single streamwise and spanwise wavelengths and wave speed) whose wall-normal shapes are determined from the principal singular function of the corresponding resolvent operator. Using the accepted scalings of the mean velocity in wall-bounded turbulent flows, we establish that the resolvent operator admits three classes of wave parameters that induce universal behaviour with Reynolds number in the low-rank model, and which are consistent with scalings proposed throughout the wall turbulence literature. In addition, it is shown that a necessary condition for geometrically self-similar resolvent modes is the presence of a logarithmic turbulent mean velocity. Under the practical assumption that the mean velocity consists of a logarithmic region, we identify the scalings that constitute hierarchies of self-similar modes that are parameterized by the critical wall-normal location where the speed of the mode equals the local turbulent mean velocity. For the rank-1 model subject to broadband forcing, the integrated streamwise energy density takes a universal form which is consistent with the dominant near-wall turbulent motions. When the shape of the forcing is optimized to enforce matching with results from direct numerical simulations at low turbulent Reynolds numbers, further similarity appears. Representation of these weight functions using similarity laws enables prediction of the Reynolds number and wall-normal variations of the streamwise energy intensity at high Reynolds numbers (${Re}_{\tau } \approx 1{0}^{3} {\unicode{x2013}} 1{0}^{10} $). Results from this low-rank model of the Navier–Stokes equations compare favourably with experimental results in the literature.


1977 ◽  
Vol 99 (3) ◽  
pp. 556-560 ◽  
Author(s):  
E. M. Sparrow ◽  
C. E. Anderson

Consideration is given to the developing laminar flow in a parallel plate channel, with the fluid being drawn from a large upstream space. The flow fields upstream and downstream of the channel inlet were solved simultaneously. A finite-difference technique was employed which was facilitated by a coordinate transformation that telescoped the broadly extended flow domain into a more compact size. For the solutions, the Reynolds number was assigned values from 1 to 1000, covering the range from viscous-dominated flows to those where both viscous and inertia effects are relevant. Streamline maps indicate that whereas a low Reynolds number flow glides smoothly into the channel, a high Reynolds number flow has to turn sharply to enter the channel, with the result that the sharply turning fluid tends to overshoot at first and then readjust. A significant amount of upstream predevelopment occurs at low and intermediate Reynolds numbers. Thus, for example, at Re = 1 and 100, the center-line velocities at inlet are, respectively, 1.37 and 1.13 times the mean velocity (the fully developed center-line velocity is 1.5 times the mean). The upstream pressure drop, measured in terms of the velocity head, is substantially increased by viscous effects at low and intermediate Reynolds numbers.


2008 ◽  
Vol 614 ◽  
pp. 425-446 ◽  
Author(s):  
MINA NISHI ◽  
BÜLENT ÜNSAL ◽  
FRANZ DURST ◽  
GAUTAM BISWAS

Laminar-to-turbulent transition of pipe flows occurs, for sufficiently high Reynolds numbers, in the form of slugs. These are initiated by disturbances in the entrance region of a pipe flow, and grow in length in the axial direction as they move downstream. Sequences of slugs merge at some distance from the pipe inlet to finally form the state of fully developed turbulent pipe flow. This formation process is generally known, but the randomness in time of naturally occurring slug formation does not permit detailed study of slug flows. For this reason, a special test facility was developed and built for detailed investigation of deterministically generated slugs in pipe flows. It is also employed to generate the puff flows at lower Reynolds numbers. The results reveal a high degree of reproducibility with which the triggering device is able to produce puffs. With increasing Reynolds number, ‘puff splitting’ is observed and the split puffs develop into slugs. Thereafter, the laminar-to-turbulent transition occurs in the same way as found for slug flows. The ring-type obstacle height, h, required to trigger fully developed laminar flows to form first slugs or puffs is determined to show its dependence on the Reynolds number, Re = DU/ν (where D is the pipe diameter, U is the mean velocity in the axial direction and ν is the kinematic viscosity of the fluid). When correctly normalized, h+ turns out to be independent of Reτ (where h+ = hUτ/ν, Reτ = DUτ/ν and $U_{\tau}\,{=}\,\sqrt{\tau_{w}/ \rho}$; τw is the wall shear stress and ρ is the density of the fluid).


Sign in / Sign up

Export Citation Format

Share Document