A higher-order numerical procedure is applied to simulate typical transient phenomena in natural gas transportation. Reliable modeling and prediction of transients features in transmission pipelines are desirable for optimal control of gas deliverability, design and implementation of active controls, and modeling of operational behavior of network peripheral equipment (e.g., chokes, valves, compressors, etc.). As an alternative to the method of characteristics (MOC) that had been widely used for modeling these systems, higher-order total variation diminishing (TVD) methods are used to model some transient problems. This class of methods has the capability of capturing fine-scale phenomena, and they do provide a better resolution of frontal discontinuities. In this study, the TVD schemes are utilized in conjunction with upwind methods. Also, in order to ensure a stable time-stepping scheme over a wide range of Courant-Friedrich-Lewy (CFL) number, a special Runge-Kutta method is employed as the base solution algorithm to integrate the highly nonlinear, hyperbolic equations which govern the transportation of natural gas in pipelines. The overall procedure is stable, robust, and accurate when applied to solve practical problems with dynamic pressure waves.

1.
Chin, W. C., 1994, Wave Propagation in Petroleum Engineering, Gulf Pub. Co., Houston, TX.
2.
Fox, S. A., 1989, Transient Flow in Pipes, Open Channels and Sewers, Elis Horwood, Ltd., England.
3.
Harten
A.
,
1983
, “
High Resolution Schemes for Hyperbolic Conservation Laws
,”
Journal of Computational Physics
, Vol.
49
, p.
357
357
.
4.
Hirose, A., and Lonngren, K. E., 1985, Introduction to Wave Phenomena, John Wiley and Sons, New York, NY.
5.
Hirsch, C., 1990, Numerical Computation of Internal and External Flows, Vol. 2, John Wiley and Sons, New York, NY.
6.
Ibraheem
S. O.
, and
Demuren
A. O.
,
1996
, “
On Bi-grid Local Mode Analysis of Solution Techniques for 3-D Euler and Navier-Stokes Equations
,”
Journal of Computational Physics
, Vol.
125
, pp.
354
377
; also, NASA-TA 106749, 1995.
7.
Jameson, A., Schmidt, W., and Turkel, E., 1981, “Numerical Solutions of the Euler Equations by Finite Volume Methods Using Runge-Kutta Time-Stepping Schemes,” AIAA Paper No. 81-1259.
8.
Johnston
R. L.
, and
Pal
S. K.
,
1972
, “
The Numerical Solution of Hyperbolic Systems Using Bicharacteristics
,”
Mathematics of Computation
, Vol.
26
, p.
377
377
.
9.
Rao
C. V. K.
, and
Eswaran
S. K.
,
1993
, “
On the Analysis of Pressure Transients in Pipelines Carrying Compressible Fluids
,”
International Journal of Pressure Vessels & Piping
, Vol.
56
, pp.
107
129
.
10.
Roe
P. L.
,
1981
, “
Approximate Riemann Solvers, Parameter Vectors and Difference Schemes
,”
Journal of Computational Physics
, Vol.
43
, pp.
357
372
.
11.
Sod
G. A.
,
1978
, “
A Survey of Several Finite Difference Methods for Systems of Nonlinear Hyperbolic Conservation Laws
,”
Journal of Computational Physics
, Vol.
27
, p.
1
1
.
12.
Steger
J. L.
, and
Warming
R. F.
,
1980
, “
Flux Vector Splitting of the Inviscid Gasdynamic Equations with Application to Finite-Difference Methods
,”
Journal of Computational Physics
, Vol.
40
, pp.
263
293
.
13.
Tan
J. K.
,
Ng
K. C.
, and
Nathan
G. K.
,
1987
, “
Application of the Center Implicit Method for Investigation of Pressure Transients in Pipelines
,”
International Journal for Numerical Methods in Fluids
, Vol.
7
, pp.
395
406
.
14.
Taylor, T. D., Wood, N. E., and Powers, J. E., 1962, “A Computer Simulation of Gas Flow in Long Pipelines,” SPE Journal, Dec. pp. 297–302.
15.
Van Leer, B., 1982, “Flux-Vector Splitting for the Euler Equations,” ICASE Report 82–30.
16.
Wylie, E. B., Streeter, V. L., and Stoner, M. A., 1974, “Unsteady-State Natural-Gas Calculation in Complex Pipe Systems,” SPE Journal, Feb., pp. 35–43.
17.
Zhou, J., and Adewumi, M. A., 1995, “Simulation of Transient Flow in Natural Gas Pipelines,” Proceedings, PSIG 27th Annual Meeting, Albuquerque, NM
This content is only available via PDF.
You do not currently have access to this content.