Mathematical Modeling of Non-Isothermal Multi-Component Fluid Flow in Pipes Applying to Rapid Gas Decompression in Rich and Base Gases

The paper presents a one-dimensional transient mathematical model of compressible non-isothermal multicomponent fluid mixture flow in a pipe. The set of the mass, momentum and enthalpy conservation equations for gas phase is solved in the model. Thermo-physical properties of multi-component gas mixture are calculated by solving the Equation of State (EOS) model. The Soave-Redlich-Kwong (SRK-EOS) model is chosen. Gas mixture viscosity is calculated on the basis of the Lee-Gonzales- Eakin (LGE) correlation. Numerical analysis of rapid gas decompression process in rich and base natural gases is made on the basis of the proposed mathematical model. The model is successfully validated on the experimental data [1]. The proposed mathematical model shows a very good agreement with the experimental data [1] in a wide range of pressure values and predicts the decompression in rich and base gas mixtures much better than analytical and mathematical models, which are available from the open source literature.

Authors:

• Evgeniy Burlutskiy

Keywords:

• Mathematical model
• Multi-Component gas mixture flow
• Rapid Gas Decompression

References:

[1] K.K. Botros, W. Studzinski, J. Geerligs, A. Glover, "Measurement of
decompression wave speed in rich gas mixtures using a decompression
tube," American Gas Association Proceedings -(AGA-2003), 2003.
[2] R.J. Eiber, T.A. Bubenik, W.A. Maxey, "Fracture control for natural gas
pipelines," PRCI Report Number L51691, 1993.
[3] R.J. Eiber, L. Carlson, B. Leis, "Fracture control requirements for gas
transmission pipelines," Proceedings of the Fourth International
Conference on Pipeline Technology, p. 437, 2004.
[4] K.K. Botros, W. Studzinski, J. Geerligs, A. Glover, "Determination of
decompression wave speed in rich gas mixtures," The Canadian Journal
of Chemical Engineering, vol. 82, pp. 880-891, 2004.
[5] K.K. Botros, J. Geerligs, J. Zhou, A. Glover, "Measurements of flow
parameters and decompression wave speed follow rapture of rich gas
pipelines, and comparison with GASDECOM," International Journal of
Pressure Vessels and Piping, vol. 84, pp. 358-367, 2007.
[6] K.K. Botros, J. Geerligs, R.J. Eiber, "Measurement of decompression
wave speed in rich gas mixtures at high pressures (370 bars) using a
specialized rupture tube," Journal of Pressure Vessel Technology, vol.
132, 051303-15, 2010.
[7] K.K. Botros, J. Geerligs, B. Rothwell, L. Carlson, L. Fletcher, P.
Venton, "Transferability of decompression wave speed measured by a
small-diameter shock tube to full size pipelines and implications for
determining required fracture propagation resistance," International
Journal of Pressure Vessels and Piping, vol. 87, pp. 681-695, 2010.
[8] GADECOM, Computer code for the calculation of gas decompression
speed that is included in "Fracture Control Technology for Natural Gas
Pipelines", by R.J. Eiber, T.A. Bubenik, W.A. Maxey, NG-18 report
208, AGA Catalog N L51691, 1993.
[9] K.H. Bendiksen, D. Maines, R. Moe, S. Nuland, "The Dynamic Two-
Fluid Model OLGA: Theory and Application," SPE Production Eng.,
vol. 6, N 2, pp. 171-180, 1991.
[10] G.B. Wallis, "One-dimensional two-phase flows," McGraw Hill, New
York, 1969.
[11] P.R.H. Blasius, "Das Aehnlichkeitsgesetz bei Reibungsvorgangen in
Fluessigkeiten," Forschungsheft, vol. 131, pp. 1-41, 1913.
[12] G. Soave, "Equilibrium constants from a modified Redlich-Kwong
equation of state," Chemical Engineering Science, vol. 27, pp. 1197-
1203, 1979.
[13] A.L. Lee, M.N. Gonzales, B.E. Eakin, "The viscosity of natural gases,"
Journal of Petroleum Technology, pp. 997-1000, 2010.
[14] S. Patankar, "Numerical heat transfer and fluid flow," Hemisphere
Publishing, New York, 1980.