0
Research Papers: Fundamental Issues and Canonical Flows

# Numerical Study of Turbulent Helical Pipe Flow With Comparison to the Experimental Results

[+] Author and Article Information
Anup Kumer Datta

and Technology,
Okayama University,
3-1-1 Tsushima-naka, Kita-Ku,
Okayama 700-8530, Japan
e-mail: akd_math_02@yahoo.com

Yasutaka Hayamizu

National Institute of Technology,
Yonago College,
4448 Hikona-cho,
Yonago-shi 683-8502, Tottori, Japan
e-mail: hayamizu@yonago-k.ac.jp

Toshinori Kouchi

and Technology,
Okayama University,
3-1-1 Tsushima-naka, Kita-Ku,
Okayama 700-8530, Japan
e-mail: kouchi@mech.okayama-u.ac.jp

Yasunori Nagata

and Technology,
Okayama University,
3-1-1 Tsushima-naka, Kita-Ku,
Okayama 700-8530, Japan
e-mail: ynagata@okayama-u.ac.jp

Kyoji Yamamoto

and Technology,
Okayama University,
3-1-1 Tsushima-naka, Kita-Ku,
Okayama 700-8530, Japan
e-mail: tetsukyo8801@earth.ocn.ne.jp

Shinichiro Yanase

and Technology,
Okayama University,
3-1-1 Tsushima-naka, Kita-Ku,
Okayama 700-8530, Japan
e-mail: yanase@mech.okayama-u.ac.jp

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received August 7, 2016; final manuscript received April 5, 2017; published online June 20, 2017. Assoc. Editor: Satoshi Watanabe.

J. Fluids Eng 139(9), 091204 (Jun 20, 2017) (13 pages) Paper No: FE-16-1503; doi: 10.1115/1.4036477 History: Received August 07, 2016; Revised April 05, 2017

## Abstract

Turbulent flow through helical pipes with circular cross section is numerically investigated comparing with the experimental results obtained by our team. Numerical calculations are carried out for two helical circular pipes having different pitches and the same nondimensional curvature δ (=0.1) over a wide range of the Reynolds number from 3000 to 21,000 for torsion parameter $β$ (=torsion $/2δ$  = 0.02 and 0.45). We numerically obtained the secondary flow, the axial flow and the intensity of the turbulent kinetic energy by use of three turbulence models incorporated in OpenFOAM. We found that the change to fully developed turbulence is identified by comparing experimental data with the results of numerical simulations using turbulence models. We also found that renormalization group (RNG) $k−ε$ turbulence model can predict excellently the fully developed turbulent flow with comparison to the experimental data. It is found that the momentum transfer due to turbulence dominates the secondary flow pattern of the turbulent helical pipe flow. It is interesting that torsion effect is more remarkable for turbulent flows than laminar flows.

<>

## References

Sreenivasan, K. R. , and Strykowski, P. J. , 1983, “ Stabilization Effects in Flow Through Helically Coiled Pipes,” Exp. Fluids, 1(1), pp. 31–36.
Webster, D. R. , and Humphrey, J. A. C. , 1993, “ Experimental Observations of Flow Instability in a Helical Coil,” ASME J. Fluids Eng., 115(3), pp. 436–443.
Webster, D. R. , and Humphrey, J. A. C. , 1997, “ Traveling Wave Instability in Helical Coil Flow,” Phys. Fluids, 9(2), pp. 407–418.
Hüttl, T. J. , and Friedrich, R. , 2001, “ Direct Numerical Simulation of Turbulent Flows in Curved and Helically Coiled Pipes,” Comput. Fluids, 30(5), pp. 591–605.
Yamamoto, K. , Hayamizu, Y. , Matsuoka, Y. , and Ohta, K. , 2005, “ Measurement of Flow Velocity and Turbulence in a Helically Curved Pipe,” Trans. Jpn. Soc. Mech. Eng., Ser. B, 71(705), pp. 1278–1285.
Cioncolini, A. , and Santini, L. , 2006, “ An Experimental Investigation Regarding the Laminar to Turbulent Flow Transition in Helically Coiled Pipes,” Exp. Therm. Fluid Sci., 30(4), pp. 367–380.
Ito, H. , 1959, “ Friction Factors for Turbulent Flow in Curved Pipes,” J. Basic Eng., Trans ASME D., 81, pp. 123–134.
Hayamizu, Y. , Yamamoto, K. , Yanase, S. , Hyakutake, T. , Shinohara, T. , and Morita, S. , 2008, “ Experimental Study of the Flow in Helical Circular Pipes: Torsion Effect on the Flow Velocity and Turbulence,” J. Therm. Sci., 17(3), pp. 193–198.
Datta, A. K. , Kouchi, T. , Hayamizu, Y. , Nagata, Y. , Yamamoto, K. , and Yanase, Y. , 2016, “ Numerical Study on the Instability of Helical Pipe Flow,” Fluid Dyn. Res. (submitted).
Wang, C. Y. , 1981, “ On the Low-Reynolds-Number Flow in a Helical Pipe,” J. Fluid Mech., 108, pp. 185–194.
Germano, M. , 1982, “ On the Effect of Torsion on a Helical Pipe Flow,” J. Fluid Mech., 125, pp. 1–8.
Kao, H. C. , 1987, “ Torsion Effect on Fully Developed Flow in a Helical Pipe,” J. Fluid Mech., 184, pp. 335–356.
Xie, D. G. , 1990, “ Torsion Effect on Secondary Flow in a Helical Pipe,” Int. J. Heat Fluid Flow, 11(2), pp. 114–119.
Tuttle, E. R. , 1990, “ Laminar Flow in Twisted Pipes,” J. Fluid Mech., 219, pp. 545–570.
Liu, S. , and Masliyah, J. H. , 1993, “ Axially Invariant Laminar Flow in Helical Pipes With a Finite Pitch,” J. Fluid Mech., 251, pp. 315–353.
Yamamoto, K. , Yanase, S. , and Yoshida, T. , 1994, “ Torsion Effect on the Flow in a Helical Pipe,” Fluid Dyn. Res., 14(5), pp. 259–273.
Yamamoto, K. , Akita, T. , Ikeuchi, H. , and Kita, Y. , 1995, “ Experiment Study of the Flow in a Helical Circular Tube,” Fluid Dyn. Res., 16(4), pp. 237–249.
Yamamoto, K. , Yanase, S. , and Jiang, R. , 1998, “ Stability of the Flow in a Helical Tube,” Fluid Dyn. Res., 22(3), pp. 153–170.
OpenFOAM, 2014, “ OpenFOAM 2.3.0 User Guide,” OpenCFD Ltd., Bracknell, UK,
Launder, B. E. , and Sharma, B. I. , 1974, “ Application of the Energy-Dissipation Model of Turbulence to the Calculation of Flow Near a Spinning Disc,” Lett. Heat Mass Transfer, 1(2), pp. 131–137.
Lien, F. S. , Chen, W. L. , and Leschziner, M. A. , 1996, “ Low Reynolds-Number Eddy-Viscosity Modelling Based on Non-Linear Stress-Strain/Vorticity Relations,” Engineering Turbulence Modelling and Experiments, Vol. 3, W. Rodi and G. Bergeles, eds., Elsevier Science B. V., Amsterdam, The Netherlands, pp. 91–100.
Etemad, S. , Sunden, B. , and Daunius, O. , 2006, “ Turbulent Flow and Heat Transfer in a Square-Sectioned U-Bend,” Prog. Comput. Fluid Dyn., 6(1–3), pp. 89–100.
Yakhot, V. , and Orszag, S. A. , 1986, “ Renormalization Group Analysis of Turbulence—I: Basic Theory,” J. Sci. Comput., 1(1), pp. 3–51.
Patankar, S. V. , and Spalding, D. B. , 1972, “ A Calculation Procedure for Heat, Mass and Momentum Transfer in Three-Dimensional Parabolic Flows,” Int. J. Heat Mass Transfer, 15(10), pp. 1787–1806.
Wada, K. , 2010, “ Turbulence Measurement of the Flow in a Helical Circular Pipes,” Master thesis, Okayama University, Okayama, Japan (in Japanese).
Lesieur, M. , Metais, O. , and Comte, P. , 2005, Large-Eddy Simulation of Turbulence, Cambridge University Press, Cambridge, UK.

## Figures

Fig. 1

(a) Helical pipe with circular cross section and (b) coordinate system

Fig. 2

Mesh on the cross section of the helical pipe

Fig. 3

Axial flow distribution on the x-axis for β = 0.02 and δ = 0.1

Fig. 6

Axial velocity distribution on the y-axis for β = 0.02 at δ = 0.1. The upper wall is on the right-hand side: (a) Re = 5320, (b) Re = 6880, (c) Re = 12,140, and (d) Re = 20,850.

Fig. 4

A view of the 3.2 pitch helical pipe showing the three-pitch downstream cross section for β = 0.02 and δ = 0.1

Fig. 5

Axial velocity distribution on the x-axis for β = 0.02 at δ = 0.1. The outer wall is on the right-hand side: (a) Re = 5320, (b) Re = 6880, (c) Re = 12,140, and (d) Re = 20,850.

Fig. 7

Vector plots of the secondary flow for β = 0.02 and δ = 0.1. Two figures on the left-hand side are the results of numerical simulations and one on the right-hand side is the experimental result: (a) is for Re = 12,140 and (b) for Re = 20,850. The outer wall is on the right-hand side.

Fig. 8

Numerical result of vector plots of the secondary flow for β = 0.02 and δ = 0.1 at Re = 200

Fig. 9

Contours of the axial velocity for β = 0.02 and δ = 0.1. Two figures on the left-hand side are the results of numerical simulations and one on the right-hand side is the experimental result: (a) is for Re = 12,140 and (b) for Re = 20,850.

Fig. 14

Turbulence intensity, (2k)1/2/〈w¯〉, for β = 0.45 and δ = 0.1 at Re = 9340 and 19,910, where k is the intensity of turbulent kinetic energy and 〈w¯〉 the mean axial velocity. Two figures on the left-hand side are the results of numerical simulations and one on the right-hand side is the experimental result. (a) is for Re = 9340 and (b) for Re = 19,910.

Fig. 15

Friction factor of the helical pipe for β = 0.45 and δ = 0.1

Fig. 10

Turbulent intensity, (2k)1/2/〈w¯〉 for β = 0.02 and δ = 0.1 at Re = 12,140 and 20,850, where k is the intensity of turbulent kinetic energy and 〈w¯〉 the mean axial velocity. Two figures on the left-hand side are the results of numerical simulations and one on the right-hand side is the experimental result: (a) is for Re = 12,140 and (b) for Re = 20,850.

Fig. 17

Equi-Q-value surface for β = 0.45 and δ = 0.1: (a) shows the pipe wall, (b) the case for Re = 200, where Q = 0.5, and (c) the case for Re = 19,910, where Q = 70

Fig. 16

Equi-Q-value surface for β = 0.02 and δ = 0.1: (a) shows the pipe wall, (b) the case for Re = 200, where Q = 0.5, and (c) the case for Re = 20,850, where Q = 70

Fig. 11

Axial velocity distribution on the x-axis for β = 0.45 at δ = 0.1. The outer wall is on the right-hand side: (a) Re = 3060, (b) Re = 5550, (c) Re = 9340, and (d) Re = 19,910.

Fig. 12

Axial velocity distribution on the y-axis for β = 0.45 at δ = 0.1. The upper wall is on the right-hand side: (a) Re = 3060, (b) Re = 5550, (c) Re = 9340, and (d) Re = 19,910.

Fig. 13

Contours of the axial velocity for β = 0.45 and δ = 0.1. Two figures on the left-hand side are the results of numerical simulations and one on the right-hand side is the experimental result: (a) is for Re = 9340 and (b) for Re = 19,910.

## Errata

Some tools below are only available to our subscribers or users with an online account.

### Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related Proceedings Articles
Related eBook Content
Topic Collections