Transmission of natural gas through high pressure pipelines has been modeled by numerically solving the governing equations for one-dimensional compressible flow using implicit finite difference methods. In the first case the backward Euler method is considered using both standard first-order upwind and second-order centered differences for the spatial derivatives. The first-order upwind approximation, which is a one-sided approximation, is found to be unstable for CFL numbers less than $1$, while the centered difference approximation is stable for any CFL number. In the second case a cell centered method is considered where flow values are calculated at the midpoint between grid points. This method is also stable for any CFL number. However, for a discontinuous change in inlet temperature, the method is observed to introduce unphysical oscillations in the temperature profile along the pipeline. A solution strategy where the hydraulic and thermal models are solved separately using different discretization techniques is suggested. Such a solution strategy does not introduce unphysical oscillations for discontinuous changes in inlet boundary conditions and is found to be stable for any CFL number. The one-dimensional flow model is validated using operational data from a high pressure natural gas pipeline.

## References

References
1.
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
.10.1016/0142-727X(87)90044-0
2.
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
.10.1243/JMES_JOUR_1972_014_045_02
3.
Wylie
,
E. B.
,
Streeter
,
V. L.
, and
Stoner
,
M. A.
,
1974
, “
Unsteady-State Natural-Gas Calculations in Complex Pipe Systems
,”
Soc. Pet. Eng.
,
14
(1), pp.
35
43
.10.2118/4004-PA
4.
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
.10.1016/0020-7403(87)90118-4
5.
Kiuchi
,
T.
,
1994
, “
An Implicit Method for Transient Gas Flows in Pipe Networks
,”
Int. J. Heat Fluid Flow
,
15
(
5
), pp.
378
383
.10.1016/0142-727X(94)90051-5
6.
Abbaspour
,
M.
, and
Chapman
,
K. S.
,
2008
, “
Nonisothermal Transient Flow in Natural Gas Pipeline
,”
ASME J. Appl. Mech.
,
75
(
3
), p.
031018
.10.1115/1.2840046
7.
Chaczykowski
,
M.
,
2010
, “
Transient Flow in Natural Gas Pipeline—The Effect of Pipeline Thermal Model
,”
Appl. Math. Model.
,
34
(4), pp.
1051
1067
.10.1016/j.apm.2009.07.017
8.
Ytrehus
,
T.
, and
Helgaker
,
J. F.
,
2013
, “
Energy Dissipation Effect in the One-Dimensional Limit of the Energy Equation in Turbulent Compressible Flow
,”
ASME J. Fluids Eng.
,
135
(6), p.
061201
.10.1115/1.4023656
9.
Luskin
,
M.
,
1979
, “
An Approximate Procedure for Nonsymmetric Nonlinear Hyperbolic Systems With Integral Boundary Conditions
,”
SIAM J. Numer. Anal.
,
16
(
1
), pp.
145
164
.10.1137/0716011
10.
Bisgaard
,
C.
,
Sørensen
,
H. H.
, and
Spangenberg
,
S.
,
1987
, “
A Finite Element Method for Transient Compressible Flow in Pipelines
,”
Int. J. Numer. Methods Fluids
,
7
(3), pp.
291
303
.10.1002/fld.1650070308
11.
Chaiko
,
M. A.
,
2006
, “
A Finite-Volume Approach for Simulation of Liquid-Column Separation in Pipelines
,”
ASME J. Fluids Eng.
,
128
(6), pp.
1324
1335
.10.1115/1.2353271
12.
Gato
,
L. M. C.
, and
Henriques
,
J. C. C.
,
2005
, “
Dynamic Behaviour of High-Pressure Natural-Gas Flow in Pipelines
,”
Int. J. Heat Fluid Flow
,
26
(5), pp.
817
825
.10.1016/j.ijheatfluidflow.2005.03.011
13.
,
A. J.
, and
Yedroudj
,
M.
,
1987
, “
A Comparison of a Finite Element Method and a Finite Difference Method for Transient Simulation of a Gas Pipeline
,”
Appl. Math. Model.
,
13
(2), pp.
79
85
.10.1016/0307-904X(89)90018-8
14.
,
A. J.
, and
Chaczykowski
,
M.
,
2001
, “
Comparison of Isothermal and Non-Isothermal Pipeline Gas Flow Models
,”
Chem. Eng. J.
,
81
(1-3), pp.
41
51
.10.1016/S1385-8947(00)00194-7
15.
Helgaker
,
J. F.
, and
Ytrehus
,
T.
,
2012
, “
Coupling Between Continuity/Momentum and Energy Equation in 1D Gas Flow
,”
Energy Proc.
,
26
, pp.
82
89
.10.1016/j.egypro.2012.06.013
16.
Barley
,
J.
,
2012
, “
Thermal Decoupling: An Investigation
,”
Proceedings to 43rd PSIG Annual Meeting
, Santa Fe, New Mexico, May 15–18.
17.
Modisette
,
J.
,
2012
, “
Instability and Other Numerical Problems in Finite Difference Pipeline Models
,”
Proceedings to 43rd PSIG Annual Meeting
, Santa Fe, New Mexico, May 15–18.
18.
Abbaspour
,
M.
,
Chapman
,
K. S.
, and
Glasgow
,
L. A.
,
2010
, “
Transient Modeling of Non-Isothermal, Dispersed Two-Phase Flow in Natural Gas Pipelines
,”
Appl. Math. Model.
,
34
(2), pp.
495
507
.10.1016/j.apm.2009.06.023
19.
Chaczykowski
,
M.
,
2009
, “
Sensitivity of Pipeline Gas Flow Model to the Selection of the Equation of State
,”
Chem. Eng. Res. Design
,
87
(12), pp.
1596
1603
.10.1016/j.cherd.2009.06.008
20.
Starling
,
K. E.
,
1973
,
Fluid Thermodynamic Properties for Light Petroleum Systems
, Gulf Publishing Company Houston, TX.
21.
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
(4), pp.
133
156
.10.1680/ijoti.1939.13150
22.
Helgaker
,
J. F.
,
2013
, “
Modeling Transient Flow in Long Distance Offshore Natural Gas Pipelines
,” Ph.D. thesis, Norwegian University of Science and Technology, Trondheim, Norway.