Research Papers: Fundamental Issues and Canonical Flows

Energy Dissipation Effect in the One-Dimensional Limit of the Energy Equation in Turbulent Compressible Flow

[+] Author and Article Information
Tor Ytrehus

Department of Energy and Process Engineering,
Norwegian University of Science
and Technology,
7491 Trondheim, Norway e-mail: tor.ytrehus@ntnu.no

Jan Fredrik Helgaker

Polytec Research Institute,
5527 Haugesund, Norway;
Department of Energy and Process Engineering,
Norwegian University of Science
and Technology,
7491 Trondheim, Norway
e-mail: jan.fredrik.helgaker@polytec.no

1Corresponding author.

Manuscript received July 6, 2012; final manuscript received January 23, 2013; published online April 8, 2013. Assoc. Editor: Ye Zhou.

J. Fluids Eng 135(6), 061201 (Apr 08, 2013) (8 pages) Paper No: FE-12-1318; doi: 10.1115/1.4023656 History: Received July 06, 2012; Revised January 23, 2013

The transportation of natural gas through high pressure transmission pipelines has been modeled by numerically solving the conservation equations for mass, momentum, and energy for one-dimensional compressible viscous heat conducting flow. Since the one-dimensional version is a result of averages over the pipe cross-section and the flow is normally turbulent, the order of averaging in space and time is an issue; in particular, for the dissipation term. The Reynolds decomposition and time averaging should be performed first, followed by the contraction to the one-dimensional version by the cross-sectional averaging. The result is a correction factor, which is close to unity, on the usual expression of the dissipation term in the energy equation. This factor will, to some extent, affect the temperature distribution along the pipeline. For low Reynolds numbers (Re104) it reduces the dissipation by as much as 7%, irrespective of roughness. For high Reynolds numbers (Re ≥ 107) and roughness in the high range of the micron decade, the dissipation is increased by 10%. If the pipeline is also thermally isolated such that the flow can be considered adiabatic, the effect of turbulent dissipation gains further importance.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


Ramsen, J., Losnegaard, S., Langelandsvik, L., Simonsen, A., and Postvoll, W., 2009, “Important Aspects of Gas Temperature Modeling in Long Subsea Pipelines,” Proceedings of the 40th PSIG Annual Meeting 2009, Galveston, TX.
Thorley, A., and Tiley, C., 1987, “Unsteady and Transient Flow of Compressible Fluids in Pipelines—A Review of Theoretical and Some Experimental Studies,” Int. J. Heat Fluid Flow, 8(1), pp. 3–15. [CrossRef]
Issa, R. I., and Spalding, D. B., 1972, “Unsteady One-Dimensional Compressible Frictional Flow With Heat Transfer,” J. Mech. Eng. Sci., 14(6), pp. 365–369. [CrossRef]
Wylie, E. B., Streeter, V. L., and Stoner, M. A., 1974, “Unsteady-State Natural-Gas Calculations in Complex Pipe Systems,” SPE J., 14, pp. 35–43. [CrossRef]
Poloni, M., Winterbone, D. E., and Nichols, J. R., 1987, “Comparison of Unsteady Flow Calculations in a Pipe by the Method of Characteristics and the Two-Step Differential Lax-Wendroff Method,” Int. J. Mech. Sci., 29(5), pp. 367–378. [CrossRef]
Kiuchi, T., 1994, “An Implicit Method for Transient Gas Flows in Pipe Networks,” Int. J. Heat Fluid Flow, 15(5), pp. 378–383. [CrossRef]
Abbaspour, M., and Chapman, K. S., 2008, “Nonisothermal Transient Flow in Natural Gas Pipeline,” ASME J. Appl. Mech., 75(3), p. 031018. [CrossRef]
Luskin, M., 1979, “An Approximate Procedure for Nonsymmetric Nonlinear Hyperbolic Systems With Integral Boundary Conditions,” SIAM (Soc. Ind. Appl. Math.) J. Numer. Anal., 16(1), pp. 145–164. [CrossRef]
Chaczykowski, M., 2010, “Transient Flow in Natural Gas Pipeline—The Effect of Pipeline Thermal Model,” Appl. Math. Model., 34, pp. 1051–1067. [CrossRef]
Osiadacz, A. J., and Chaczykowski, M., 2001, “Comparison of Isothermal and Non-Isothermal Pipeline Gas Flow Models,” Chem. Eng. J., 81, pp. 41–51. [CrossRef]
Chaczykowski, M., 2009 “Sensitivity of Pipeline Gas Flow Model to the Selection of the Equation of State,” Chem. Eng. Res. Des., 87, pp. 1596–1603. [CrossRef]
Starling, K. E., 1973, Fluid Thermodynamic Properties for Light Petroleum Systems, Gulf Pub., Houston, TX.
Colebrook, C. F., 1939, “Turbulent Flow in Pipes, With Particular Reference to the Transition Region Between the Smooth and Rough Pipe Laws,” J. Inst. Civil Eng., 11, pp. 133–156.
White, F. M., 2006, Viscous Fluid Flow, 3rd ed., McGraw-Hill, New York.
Hinze, J. O., 1975, Turbulence, 2nd ed., McGraw-Hill, New York.
Zagarola, M. V., and Smits, A. J., 1998, “Mean-Flow Scaling of Turbulent Pipe Flow,” J. Fluid Mech., 373, pp. 33–79. [CrossRef]
Spalding, D. B., 1961, “A Single Formula for the ‘Law of the Wall’,” ASME J. Appl. Mech., 28(3), pp. 455–458. [CrossRef]


Grahic Jump Location
Fig. 1

Velocity profile of the inner wall in the inner variables u+ and y+. The figure taken from White [14] (p. 420).

Grahic Jump Location
Fig. 2

Coordinate system used in integrating Eq. (17) over the pipe cross-section

Grahic Jump Location
Fig. 3

Correction factor F as a function of the Reynolds number, surface roughness, and pipe diameter. Results are given for values of ε/D * 106 equal to 20, 10, 5, 3, and 1. Values for F have been computed from Eq. (27).

Grahic Jump Location
Fig. 4

Local errors for p, m·, and T as a function of the grid points N

Grahic Jump Location
Fig. 5

Outlet boundary condition for the mass flow for the smooth pipe

Grahic Jump Location
Fig. 6

Results for the smooth pipe: (top) outlet pressure, and (bottom) inlet mass flow

Grahic Jump Location
Fig. 7

Outlet temperature for a smooth pipe with and without the correction factor at Re = 3 × 107 and ε/D = 10−6. The difference in the outlet temperature is approximately 0.1 °C.

Grahic Jump Location
Fig. 8

Boundary conditions: (left) inlet mass flow (f1(t)), and (right) inlet temperature (f2(t))

Grahic Jump Location
Fig. 9

The Reynolds number and turbulent dissipation factor at the pipe outlet

Grahic Jump Location
Fig. 10

Outlet temperature as a function of time with and without correcting the dissipation term

Grahic Jump Location
Fig. 11

Temperature profile along the pipeline at t = 20 h with and without correcting the dissipation term

Grahic Jump Location
Fig. 12

The difference between temperatures for the corrected and noncorrected dissipation term for the heat conducting (U = 2.5) and adiabatic (U = 0) cases

Grahic Jump Location
Fig. 13

(left) Modeled inlet pressure, and (right) modeled outlet mass flow



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 eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In