Abstract Unsteady-state or transient two-phase flow, caused by any change in rates, pressures or temperature at any location in a two-phase flow line, may last from a few seconds to several hours. In general, these changes are an order of magnitude longer than the transient encountered during single-phase flow. The primary reason for this phenomenon is that the velocity of wave propagation in a two-phase mixture is significantly slower. Interfacial transfer of mass, momentum and energy further complicate the problem. It is primarily due to the numerical difficulties anticipated in accurately modeling transient two-phase flow that the state of the art in this important area is restricted to a handful of studies with direct applicability to petroleum and gas engineering. A limited amount of information on the subject of two-phase transport phenomena is available in the petroleum engineering literature. Most of the publications for two-phase flow of gas assume that temperature is constant over the entire length of the pipeline. This study is the first effort to simulate the non-isothermal, one-dimensional, transient homogenous two-phase flow gas pipeline system using two-fluid conservation equations. The modified Peng–Robinson equation of state is used to calculate the vapor–liquid equilibrium in multi-component natural gas to find the vapor and liquid compressibility factors. Mass transfer between the gas and the liquid phases is treated rigorously through flash calculation, making the algorithm capable of handling retrograde condensation. The liquid droplets are assumed to be spheres of uniform size, evenly dispersed throughout the gas phase. The method of solution is the fully implicit finite difference method. This method is stable for gas pipeline simulations when using a large time step and therefore minimizes the computation time. The algorithm used to solve the non-linear finite difference thermo-fluid equations for two-phase flow through a pipe is based on the Newton–Raphson method. The results show that the liquid condensate holdup is a strong function of temperature, pressure, mass flow rate, and mixture composition. Also, the fully implicit method has advantages, such as the guaranteed stability for large time step, which is very useful for simulating long-term transients in natural gas pipeline systems.
[1]
Michael A. Adewumi,et al.
Transients in Gas-Condensate Natural Gas Pipelines
,
1998
.
[2]
Barry D Ganapol.
Multiphase Transport Fundamentals, Reactor Safety Applications
,
1982
.
[3]
Michael A. Adewumi,et al.
Engineering design of gas-condensate pipelines with a compositional hydrodynamic model
,
1990
.
[4]
A. Dukler,et al.
Transient gas‐liquid flow in horizontal pipes: Modeling the flow pattern transitions
,
1978
.
[5]
R. Clift,et al.
Bubbles, Drops, and Particles
,
1978
.
[6]
D. Peng,et al.
A New Two-Constant Equation of State
,
1976
.
[7]
Haroun Mahgerefteh,et al.
Fast numerical simulation for full bore rupture of pressurized pipelines
,
1999
.
[8]
N. Chen.
An Explicit Equation for Friction Factor in Pipe
,
1979
.
[9]
Haroun Mahgerefteh,et al.
Modeling fluid phase transition effects on dynamic behavior of ESDV
,
2000
.
[10]
Ben Asante.
Two Phase Flow: Accounting for the Presence of Liquids in Gas Pipeline Design
,
2004
.
[11]
J. Hinze,et al.
Critical speeds and sizes of liquid globules
,
1949
.