Parameters calculation of turbulent fluid flow in a pipe of a circular cross section

International Journal of Innovation Engineering and Science Research
www.ijiesr.com
Volume 2 Issue 6 November-December 2018 38|P a g e
ABSTRACT
Parameters calculation of turbulent fluid flow in a
pipe of a circular cross section
Denis Chemezov, Svetlana Tyurina, Irina Medvedeva, Lyudmila Smirnova, Elena Bogomolova,
Margarita Bakhmeteva, Alexandra Strunina, Nina Melenteva
Vladimir Industrial College
Vladimir, Russian Federation
The dependencies of total pressure, velocity, vorticity, turbulent length, turbulent dissipation, turbulent viscosity,
turbulent energy and turbulent time of moving fluid from a straight pipe length of a circular cross section are
presented in graphical and mathematical forms. Changing analysis of considered parameters was performed at
mass flow rates of 0.45, 1.0 and 1.5 kg/s. A transition boundary of laminar flow of fluid to turbulent flow is at the
distance of 2
/5 of length from the inlet of the pipe (at accepted total length of the pipe of 1000 mm).
Keywords – fluid, turbulent flow, a pipe, a model, a section.
I. INTRODUCTION
At present, a question of sustainability of laminar/transient flow regimes of fluid in pipelines has not
been fully studied. According to numerous experimental data, it is determined that even in small
straight sections of the pipeline, fluid flow changes from laminar to turbulent [1 – 7].
Flow pattern of fluid is determined by viscosity, flow velocity, a cross-sectional area of the pipe,
composition of fluid mixture and other parameters. Also, flow pattern of fluid depends on the ratio of
accelerating force and viscous friction force.
Conducting of experiments on a research of fluid flow regimes in the production, and even in the
laboratory conditions, is difficult. A computer simulation of hydrodynamic processes of fluid flow in the
pipelines allows not only to present a visual model of fluid flow, but also to obtain accurate or
approximate mathematical equations describing pattern of each flow regime.
II. MATERIAL AND METHOD
Intensity changing of the parameters of turbulent fluid flow in the pipe of the circular cross section was
determined by calculations in the Flow Simulation special computer program. Three solid-state pipe
models with the same overall dimensions (internal diameter is 30 mm and total length is 1000 mm)
were built for the computer simulation.
Different mass of the incompressible fluid model (water) at temperature of 293.2 K was supplied per
unit of time (second) in the inlet of the pipe models. Specified mass flow rates for three pipe models
are presented in the table I.
TABLE I. MASS FLOW RATES OF FLUID.
The pipe model The first The second The third
Mass flow rate of fluid QM, kg/s 0.45 1.0 1.5
Fluid flow direction was taken normal to a face of the inlet of the pipe model. Turbulence intensity It
(2%) and turbulence length Lt (0.0008 m) were taken as the basic parameters of developed turbulent
fluid motion in the pipe model. Turbulence intensity of fluid flows depends on characteristics of the
pipe and the Reynolds number Re. Fluid pressure at the outlet of the pipe model was taken of 101325
Pa. Laminar and turbulent fluid flows were taken into account in the calculation. The cavitation
process in the calculation was not accepted.
An inner wall of the pipe model had the properties of the adiabatic wall with the perfectly smooth
surface (the surface roughness is 0 µm). Initial temperature of the inner surface of the pipe model was
taken of 293.2 K.

Chemezov, D. et al. “International Journal of Innovation Engineering and Science Research”
Accuracy of the calculation results was average when a number of finite elements of the fluid model
equal to 3264 and the pipe model equal to 528.
III. RESULT AND DISCUSSION
Flow pattern of fluid in the pipe model of the circular cross section is presented in the form of vectors
(the Fig. 1).
Fig. 1. Water flow at the straight section of the three-dimensional pipe model of the circular cross section.
Vortex formation of fluid flows is observed in a center part of a calculated field of the pipe model,
laminar regime prevails at the inlet and the outlet.
The dependencies of the changing values of the turbulent fluid flow parameters from length of the
model pipe of the circular cross section are presented in the Fig. 2 – 5. The data were obtained from
the axial line of fluid motion in the pipe model.
a) b) c)
Fig. 2. The dependencies of flow velocity and total pressure of fluid from the pipe length: a) at QM = 0.45 kg/s; b) at QM = 1.0 kg/s;
c) at QM = 1.5 kg/s.
a) b) c)
Fig. 3. The dependencies of turbulent length and vorticity of fluid from the pipe length: a) at QM = 0.45 kg/s; b) at QM = 1.0 kg/s; c)
at QM = 1.5 kg/s.

a) b) c)
Fig. 4. The dependencies of turbulent viscosity and turbulent dissipation of fluid from the pipe length: a) at QM = 0.45 kg/s; b) at
QM = 1.0 kg/s; c) at QM = 1.5 kg/s.
a) b) c)
Fig. 5. The dependencies of turbulent energy and turbulent time of fluid from the pipe length: a) at QM = 0.45 kg/s; b) at QM = 1.0
kg/s; c) at QM = 1.5 kg/s.
Based on the analysis of the dependencies graphs of flow velocity and total fluid pressure, turbulent
length and vorticity of fluid, turbulent viscosity and turbulent dissipation of fluid, turbulent energy and
turbulent time of fluid flow from length of the pipe model, it could be argued that:
1. The values of the considered parameters of turbulent fluid flow increase with increasing of mass flow
rate at the distance from the inlet to the outlet of the pipe model. Herewith, turbulent time of fluid in the
each section of the pipe model decreases.
2. Pressure, flow velocity, turbulent dissipation and turbulent viscosity of moving fluid have almost the
same pattern of changing over entire length of the pipe model. The other parameters of turbulent fluid
flow in the pipe model have different regularities.
Let us consider flow pattern of fluid in the second pipe model. Total length of the pipe model was
divided into 10 equal sections. Length of the each section was 100 mm (0.1 m). Herewith, 0 was taken
by the inlet of the pipe model, and 1 was taken by the outlet. Let us write down the parameters
changing of turbulent fluid flow in the each section of the pipe model in the functions form (the tables II
– IV): Ptot(Lp), εt(Lp), Lt(Lp), u(Lp), tt(Lp), νt(Lp), kt(Lp) and ω(Lp).
TABLE II. THE PARAMETERS CHANGING OF TURBULENT FLUID FLOW IN THE FIRST TO THE FOURTH SECTIONS OF THE PIPE MODEL.
1
101095579
2029501
Lp.
etotP

 2
91001511
2029502
Lp.
etotP

 3
91003821
2029503
Lp.
etotP

 4
91013461
2029504
Lp.
etotP


1
00020
0.2983=1
Lp.
et 2
00010
29840=2
Lp.
e.t 3
00010
298903
Lp.
e.=tε 4
51016371
300804
Lp.
e.=tε


1
61037629
001701
Lp.
e.tL

 2
61097637
001702
Lp.
e.tL

 3
61023876
001703
Lp.
e.tL

 4
61091883
001704
Lp.
e.tL


1
00290
003601
Lp.
e.u  2
00280
003702
Lp.
e.u  3
00270
003703
Lp.
e.u  4
00240
003804
Lp.
e.u 
1
51060336
071101
Lp.
e.tt

 2
51039325
071102
Lp.
e.tt

 3
51090383
071103
Lp.
e.tt

 4
61063392
070904
Lp.
e.tt


1
51062083
147501
Lp.
e.tv

 2
51098792
147502
Lp.
e.tv

 3
51021732
147503
Lp.
e.tv

 4
71047528
147804
Lp.
e.tv


1
00010
02201
Lp.
e.tk  2
51039048
02202
Lp.
e.tk

 3
51023246
02203
Lp.
e.tk

 4
61066436
022104
Lp.
e.tk


1
00170
194201
Lp.
e.ω  2
00170
194402
Lp.
e.ω  3
00160
195103
Lp.
e.ω  4
00170
19304
Lp.
e.ω 
TABLE III. THE PARAMETERS CHANGING OF TURBULENT FLUID FLOW IN THE FIFTH TO THE SEVENTH SECTIONS OF THE PIPE MODEL.
5
91034251
2029505
Lp.
etotP

 6
51059816
2066306
Lp.
etotP

 7
00260
7356407
Lp.
etotP


5
00020
306705
Lp.
e.=tε

6
01810
003406
Lp.
e.=tε 7
00880
502707
Lp.
e.=tε
5
71060222
001705
Lp.
e.tL

 6
01670510759126
Lp.
e.tL  7
00260
05907
Lp.
e.tL 
5
00120
004305
Lp.
e.u  6
01310
000206
Lp.
e.u  7
00890
003307
Lp.
e.u 
5
51067117
070405
Lp.
e.tt

 6
00510
020206
Lp.
e.tt  7
00120
635507
Lp.
e.tt


5
51043555
148505
Lp.
e.tv

 6
02790
000106
Lp.
e.tv  7
00630
214187
Lp.
e.tv 
5
00010
022405
Lp.
e.tk

 6
0230510385676
Lp.
e.tk  7
00760
319407
Lp.
e.tk 
5
00170
196105
Lp.
e. 6
00310
456906
Lp.
e.

 7
00240
087407
Lp.
e.
TABLE IV. THE PARAMETERS CHANGING OF TURBULENT FLUID FLOW IN THE EIGHTH TO THE TENTH SECTIONS OF THE PIPE MODEL.
8
00380
141950008
Lp.
etotP

 9
00460
21623009
Lp.
etotP

 10
00670
762440010
Lp.
etotP


8
0040
645868
Lp.
e.=tε 9
00190
802219
Lp.
e.=tε 10
00110
51483410
Lp.
e.=tε
8
00140
108208
Lp.
e.tL  9
00040
187909
Lp.
e.tL  10
00010
2273010
Lp.
e.tL 
8
010
1834938
Lp.
e.u

 9
00360
044209
Lp.
e.u  10
00380
0396010
Lp.
e.u 
8
00040
402808
Lp.
e.tt

 9
00030
391509
Lp.
e.tt

 10
00030
3813010
Lp.
e.tt


8
00320
7102968
Lp.
e.tv  9
00120
01353009
Lp.
e.tv  10
00050
54245010
Lp.
e.tv 
8
00360
675528
Lp.
e.tk  9
00160
535589
Lp.
e.tk  10
00080
16071310
Lp.
e.tk 
8
02970
32363008
Lp.
e

 9
00780
604579
Lp.
e.

 10
00920
0003010
Lp.
e.
In the presented functions are conditionally designated: Ptot1, Ptot2, Ptot3, …, Ptot10 – total pressure of
fluid in the first, the second, the third, …, the tenth sections of the pipe model, Pa; Lp1 = (0;0.1], Lp2 =
(0.1;0.2], Lp3 = (0.2;0.3], Lp4 = (0.3;0.4], Lp5 = (0.4;0.5], Lp6 = (0.5;0.6], Lp7 = (0.6;0.7], Lp8 = (0.7;0.8], Lp9
= (0.8;0.9], Lp10 = (0.9;1.0] – the length ranges of the each section of the pipe model, m (the values in
the range are location coordinates of the corresponding section, starting from the inlet of the pipe
model); εt1, εt2, εt3, …, εt10 – turbulent dissipation of fluid in the first, the second, the third, …, the tenth
sections of the pipe model, W/kg; Lt1, Lt2, Lt3, …, Lt10 – turbulent length of fluid flow in the first, the

second, the third, …, the tenth sections of the pipe model, m; u1, u2, u3, …, u10 – flow velocity of fluid in
the first, the second, the third, …, the tenth sections of the pipe model, m/s; tt1, tt2, tt3, …, tt10 – turbulent
time of fluid flow in the first, the second, the third, …, the tenth sections of the pipe model, s; νt1, νt2, νt3,
…, νt10 – turbulent viscosity of fluid in the first, the second, the third, …, the tenth sections of the pipe
model, Pa·s; kt1, kt2, kt3, …, kt10 – turbulent energy of fluid in the first, the second, the third, …, the tenth
sections of the pipe model, J/kg; ω1, ω2, ω3, …, ω10 – vorticity of fluid in the first, the second, the third,
…, the tenth sections of the pipe model, 1/s.
In the first four sections, non-gradient flow is observed, i.e. fluid in these sections does not completely
fill the cross section of the pipe model. Thickness of a boundary layer (the thin near-wall layer)
increases at non-gradient flow of fluid. In the remaining sections of the pipe model, an area with
positive pressure gradient was determined in which fluid flow was slowing down. Fluid flow slows in the
boundary layer and in the flow core (external flow). Fluid slowing is more quickly in the boundary layer
by increasing its thickness.
The general analytical formula was obtained after the analysis of the calculated mathematical
equations for determining of the parameters of turbulent fluid flow in the pipe model of the circular
cross section:
  )1( kxkx
e
k
A
dxAе ,
where А is numerical coefficient of the function; k is power coefficient of the function.
IV. CONCLUSION
1. Changing pattern of flow velocity and total pressure, and, consequently, turbulent dissipation and
viscosity of moving fluid does not depend on mass flow rate.
2. The more mass flow rate of fluid, the more strongly signs of turbulence are observed (the Fig. 3, c
and the Fig. 5, c).
3. At the distance of
2
/5 of length from the inlet of the pipe, fluid flow is closer to laminar. At the
distance of
3
/5 of length from the outlet of the pipe, interaction of direct and reverse fluid flows occurs,
which causes vortices, and, consequently, laminar regime turns into turbulent.
REFERENCES
[1] D. Chemezov, Calculation of pressure losses of liquid at a cylindrical straight pipeline section, ISJ Theoretical & Applied
Science, 12 (56), 2017, 19 – 22. DOI 10.15863/TAS.2017.12.56.5
[2] D.A. Chemezov, Hydrodynamic characteristics of water flow in straight and curved sections of a pipeline, Modern
materials, equipment and technology [Text]: Materials of the 4
th
International scientific and practical conference, South-
West State University, CJSC «University book», Kursk, 2015, 468 – 471. ISBN 978-5-9906195-4-8.
[3] D. Chemezov, N. Palev, Analytical models of the turbulent fluid flow in a circular pipe, ISJ Theoretical & Applied
Science, 09 (41), 2016, 77 – 84. DOI 10.15863/TAS.2016.09.41.12
[4] V. Oshovskiy, A. Dyubanov, Computer modelling of the hydrodynamic effects arising in the narrowing device, Scientific
Works of Donetsk National Technical University, 2(21), 2013, 168 – 169. ISSN 2074-6652.
[5] I. Tsukanov, V. Shapiro, S. Zhang, A Mesh free Method for Incompressible Fluid Dynamics Problems, Int. J. Numer.
Meth. Eng., 2003, Vol. 58.
[6] D.S. Henningson, P.J. Schmid, Stability and transition in shear flows, New-York, 2001.
[7] A.V. Proskurin, A.M. Sagalakov, Viscous Fluid Flow in a Coaxial Pipe, Izvestiya of Altai State University Journal, 2016,
62 – 68. DOI 10.14258/izvasu(2016)1-10

Parameters calculation of turbulent fluid flow in a pipe of a circular cross section

More Related Content

What's hot (18)

Similar to Parameters calculation of turbulent fluid flow in a pipe of a circular cross section (20)

More from International Journal of Innovation Engineering and Science Research (20)

Recently uploaded (20)

Parameters calculation of turbulent fluid flow in a pipe of a circular cross section