Open access peer-reviewed chapter

Modeling and Simulation for Steady State and Transient Pipe Flow of Condensate Gas

By Li Changjun, Jia Wenlong and Wu Xia

Submitted: November 29th 2010Reviewed: April 5th 2011Published: September 22nd 2011

DOI: 10.5772/23621

Downloaded: 4566

1. Introduction

Condensate gas is mainly demonstrated by methane. However, it also contains a lot of heavier contents like C5 or C5+ and some non-hydrocarbon mixture as well (Mokhatab et al, 2006). After recovering from gas wells, condensate gas needs liquid separation, gas purification and condensate stabilization treatment in the processing plant to meet the quality requirements. Processing plants far away from the gas well with long distances of two-phase flow in one condensate gas pipeline will take less investment than adjacent process plant with two single phase pipelines which are dry gas pipeline and liquid phase pipeline (Li, 2008).

If the operation temperature somewhere in the condensate gas pipeline is lower than the gas dew point, liquid condensation would occur, subjecting the pipeline to two phase flow (Potocnik, 2010). While gas and its condensate flow simultaneously, mass transfer takes place continuously due to the change in pressure and temperature conditions. This leads to compositional changes and associated fluid property changes and also makes the hydraulic and thermal calculations of condensate gas more complex than normal gas. The condensate gas pipeline model which is established and solved based on the principle of fluid mechanics can simulate hydraulic and thermal parameters under various operation conditions. By means of technical support, this model is of great importance in the pipeline design and safety operation aspects (Mokhatab, 2009).


2. Thermodynamic model

The purpose of the thermodynamic model is three-fold. First, it defines the transition between single phase/two phase conditions (point of condensate inception in the pipeline or gas dew point). Second, it is used for the prediction of properties for the flowing fluids (gas and its condensate). And lastly, it derives the mass exchange between the flowing phases (Adewumi et al, 1990 ; Estela-Uribe et al, 2003). This work uses the BWRS equation of state (EOS) to implement the thermodynamic model as it has proven reliable for gas condensate system (McCain et al, 1990). Most property predictions are derived from the equation of state (i.e., densities values, densities values and their derivatives with respect to pressure and temperature, departure enthalpies, heat capacities, and Joule-Thompson coefficients). Additionally, phase equalibria are calculated on the basis of flash calculation method. Expressions for such parameters as fugacity are elaborated in standard textbooks, where the theory and relevant procedures for flash calculation are well documented (API, 2005).

3. General model of condensate gas pipeline

3.1. Basic assumptions

One of the most fundamental approaches used to model two phase flow is the two-fluid model (Ayala et al, 2003). It consists of separate mass, momentum and energy conservation equations written for each of the phases. This results in up to six differential equations. Furthermore, the model is more complicate as relating parameters in the equations vary with the fluid flowing. Thus, simple and reasonable assumptions could help to reduce the unknowns in the model (Hasan et al, 1992). The hydrodynamic model has three major inherent assumptions:

  1. Gas and liquid average flow velocities are calculated according to the section area occupied by each phase respectively.

  2. Mass transfer takes place between gas and liquid phase. However, two phases are assumed to be at thermodynamic equilibrium at every point within the pipe.

  3. In the transient process, the pipeline assumed to be isothermal. The liquid hold up within pipeline, the wall friction force of each phase, the drag force between two phases are assumed to be the same as steady state.

3.2. Description of general model

The general model for two phase flow is built on the basis of mass, momentum and energy conservation equations of each phase (Li et al, 2009; Ayala et al, 2003).

Gas phase continuity equation:


Liquid phase continuity equation:


Where, Δm˙glis mass rate of phase change from gas to liquid, (kg/s.m); Δm˙lgis mass rate of phase change from liquid to gas, (kg/s.m). HLis liquid hold up; φis gas hold up; ..ρ density of the gas phase, kg/m3; ρlis density of the liquid phase, kg/m3; Ais pipe cross sectional area, m2; tis time, s; x is length along the pipe length, m; wgis velocity of the gas phase, m/s; wlis velocity of the liquid phase, m/s.

Δm˙glis defined as (3).

Where, Pis pressure, Pa; Tis temperature, K; mgis gas mass fraction in two-phase fluid system; mlis liquid mass fraction in two-phase fluid system; Ysis mass fraction of gas.

Gas phase momentum equation.


Where, Agis pipe cross sectional area occupied by gas phase, m2;Fgwis wall shear force of gas phase, N/m3; Fgiis interfacial drag force on gas phase, N/m3; gis acceleration of gravity, m/s2; θis pipeline slope, rad.

Liquid phase momentum equation.


Where, Alis pipe cross sectional area occupied by liquid phase, m2; Flwis wall shear force of liquid phase, N/m3; Fliis interfacial drag force on liquid phase, N/m3; wais transition velocity between the gas phase and liquid phase, m/s.

Gas-liquid phase mixture energy equation


Where, K is overall heat transfer co-efficiency, W/(m2 K); D is external diameter, m; T0is environmental temperature, K; hlis enthalpy of liquid phase in pipeline, kJ/kg; hgis enthalpy of the gas phase in pipeline, kJ/kg; zis pipeline elevation, m.

Add (7) and (8), obtain:


Equation (10) is transformed into (11) through eliminating the pressure terms.


In view of the slow transient behaviour in condensate gas pipeline, we can obtain equation (12) by ignoring velocity variation terms in equation (11) (Li et al, 1998):


Equations (1), (2), (9), (10), (12) construct the basic model for condensate gas pipe flow simulation.

4. Constitutive equations

The condensate gas flow model is one dimensional two-fluid multiphase hydrodynamic model which adapts to different flow patterns in pipeline. According to Cindric and Shoham, the flow patterns in horizontal pipeline are stratified flow, intermittent flow, annular flow, dispersed flow and these in vertical pipeline are bubble flow, slug flow, churn flow, annular flow (Mokhatab et al, 2006). Because of the constitutive equations is dependent on the flow pattern, one of the greatest difficulties in the analysis of two-phase flow in pipeline is defining appropriate constitutive equations for relating relevant forces-such as the steady drag force and interfacial force.

Considering low liquid hold up, the flow pattern in the condensate gas pipeline is stratified flow which has explicit interface between the liquid and gas phase, as depicted in Fig.1. Then, we can obtain the calculation methods of unknowns which dependent on the constitutive equations (Taitel et al, 1995 ; Chen et al, 1997 ).

Figure 1.

Stratified flow in condensate gas pipeline

Wall shear force of each phase is expressed as follow.


Where, k=gwhen the equation is applied for the gas phase; k=lwhen the equation is applied for the liquid phase; Skis defined as follow:

Sk=totalwallareawettedbyphasekTotoal volumeE14

τkwis defined as:


In which,λis Fanning factor which is calculated by Colebrook & White empirical correlation.

If Reynolds Number Re≤ 2000


If Re≥ 2000


Where, εis absolute roughness of pipeline wall, m.

Interfacial force between phases is defined as follow:




The interfacial friction factor λiis calculated with Hanrrity correlation.If wsgwsgt


If wsgwsgt




Where, wsgis reduced velocity of the gas phase, m/s; wsg.tis reduced velocity for indentifying the transition from stratified flow pattern to smooth stratified flow pattern, m/s; QGis flow rate of the gas phase, m3/s.

5. Steady state analysis of condensate gas pipeline

5.1. Basic equations

While steady operation, the variation of each parameter in equations (1), (2), (9), (10) with time can be ignored. Expand the equations above and the following equations used for steady state simulation can be obtained (Li et al, 2009):

Gas phase continuity equation:


Liquid phase continuity equation:


Gas-liquid phase mixture momentum equation:


Gas-liquid phase mixture energy equation:


The system of simultaneous differential equations composed of (25)-(28) can be written in their non-conservative form.


Where A=[a11a12a13a14a21a22a23a24a31a32a33a34a41a42a43a44],D=[D1D2D3D4],U=[PTwgwl],a11=Awgφ(ρgP)T,a12=Awgφ(ρgT)P,a13=Aρgϕ,a14=0,a21=AwlHL(ρlP)T,a22=AwlHL(ρlT)P,a23=0,a24=AρlHL,a31=A,a32=0,a33=ρgwgφA,a34=ρlwlHLA,a41=ρgwgφA(hgP)T+ρlwlHLA(hlP)Ta42=ρgwgφA(hgT)P+ρlwlHLA(hlT)P,a43=ρgwg2φA,a44=ρlwl2HLAD1=Δm˙glAρgwgdφdx,D2=Δm˙lgAρlwldHLdx,D3=FgwFlw(ρgφ+ρlHL)gAsinθD4=(ρggwgφA+ρlgwlHLA)dzdxKπD(TT0)Δm˙gl(hghl+wg2wl22)

5.2. Model solving

Steady state condensate gas model is formed by 5 equations which are (29) and (12). There are five unknowns, liquid holdup (HL), pressure (P), temperature (T), gas and liquid velocity (wgand wl), in the model. Thus, the closure of the model is satisfied.

To solve the model, the liquid hold up is obtained by solving (12) firstly. And then, pressure (P), temperature (T), velocity of the gas phase (wg), and the velocity of the liquid phase (wl) are obtained by solving (29). The procedures for solving (29) are presented in details as follow:

  1. Suppose the pipeline is composed of a lot of pipes with different slope. Divide each pipe into small blocks with the step length of △x and input the start point data.

  2. Establish steady equation (29) on each block section. Input the boundary conditions at the initial point of pipeline (pressure, temperature, gas velocity, and liquid velocity). According to the thermodynamic model, calculate the thermophysic parameters such as density of the gas and liquid phase, gas fraction. Because there is no slip between the two phases at initial point, the liquid hold up can be gained by its relationship with mass flow rate of the gas phase and liquid phase.

  3. Set dU/dxas unknowns, and simplify (29) with Gaussian elimination method, then we can obtain more explicit form of (29).

  4. Work out pressure (Pi), temperature (Ti), gas and liquid velocity (wgiand wli) by four-order Runge - Kutta Method.

  5. Figure out liquid holdup (HLi)by equation (12).

  6. Resolve equations (29) by Adams predictor-corrector formula until the reasonable unknowns of this grid section are all gotten.

  7. Repeat the second step to the sixth step until reach the last block section which is also the end of this pipeline.

In order to make the numerical calculation converges more quickly, the Adams predictor-corrector and Runge - Kutta Method should be applied simultaneous. As the two methods have four-order accuracy, the desired accuracy also can be improved. The flow chart of the whole solving procedures is depicted in Fig.3.

Figure 3.

Solving procedures of steady state model


6. Transient analysis of condensate gas pipeline

6.1. Basic equations of transient analysis

Opposite to the steady state simulation, the parameters in the general model are dependent on time. Expand (1), (2) and (10), and the following equations can be obtained (Masella et al, 1998).

Gas phase continuity equation:


Liquid phase continuity equation:


Momentum equation


The transient flow model can be represented by (31) ~ (33) and (12). The unknowns are the pressure P, flow velocity of the gas phasewg, flow velocity of the liquid phase wland liquid holdupHL. Notice that (31) ~ (33) are a set of partial differential equations so that they can be recast to the following matrix form.


Where, A=[a11a12a13a21a22a23a31a32a33],B=[b11b12b13b21b22b23b31b32b33],D=[D1D2D3],U=[Pwgwl]a11=Awgφ(ρgP)T,a12=Aρgϕ,a13=0,a14=Aρgwg,a21=AwlHL(ρlP)T,a22=0,a23=AρlHL,a24=Aρlwl,a31=A,a32=ρgwgφA,a33=ρlwlHLA,.a34=0.,b11=Aϕ(ρgP)T,b12=0,b13=0,b14=Aρg,b21=AHL(ρlP)T,b22=0,b23=0,b24=Aρl,b31=0,b32=ρgϕA,b33=ρlHLA,b34=0D1=Δm˙glb14ϕta14ϕx,D2=Δm˙lgb24φta24HLx,D3=FgwFlw(ρgϕ+ρlHL)gAsinθThe characteristic determinant of (34) is as follow:


Where, λ0and λ1are the eigenvalues of Aand B.If λ=λ0/λestablished, one solution of λis 0 and the other two non-zero solutions should be expressed as follow:


The condition for having real solutions of (36) is:


In (37), (Pρl)Tand (Pρg)Trepresent square of liquid and gas isothermal wave velocities separately. Thus the right side of (37) is greater than the wave velocity, which means (37) can be established for raw gas pipelines.

According to the analysis above, it can be concluded that under the giving range of operation conditions, all the eigenvaluse of (35) are real different number. Therefore, (34) is strictly hyperbolic and this type will not change as its well posedness.

6.2. Boundary conditions

The boundary conditions include the input and output conditions as well as conditions, such as valves shut off, compressors shut down and etc., which will lead to operating changes. For this transient model, the boundary conditions are considered as: the starting pressure and flow rate over time, ending pressure and flow rate over time.


6.3. Solutions of the transient model

Equation (34) is a set of non-liner partial differential equations and the method of characteristic (MOC) and implicit difference method are always adopted to gain numerical solutions. By both methods have advantages and disadvantages. There is strict restriction with the time length of each step (Fig.4), thus the stability of the implicit difference method is better. But at each time layer, differential equations of all the grids should be solved simultaneously, which causes the computing time is large. Different with the implicit difference method, for MOC, the time length of each step is restricted in a relatively short range by its stability conditions. The advantage of MOC is the unknowns of each grid at each time layer can be solved dependently, and the amount of equations solved simultaneously is reduced. Considering the transient feature of low-liquid loading multiphase flow in pipeline is slow transient behaviour which means the transient process lasts long time, the former method is adopted.

Mesh the length variable xand time variable tinto grids as shown in Fig. 4. In each block, the differential equations of (31) ~ (33) can be obtained by the implicit different method (Li et al, 2011).

Figure 4.

Differential blocks of pipelines

Gas phase continuity equation:


In which ξis the weight coefficient. Setλ=2ΔtΔx, then above equation reduces to:


Similarly, liquid phase continuity equation is as follow:


Momentum equation:


Based on the above three differential equations, the pressure, flow velocity of the gas phase as well as that of the liquid phase can be obtained. The solution procedures are as follow:

  1. Divide pipeline system into series of pipeline segments with different angles. Mesh each segment into blocks with the step length of △x;

  2. Get the operation parameters on each grid node by steady state solutions;

  3. Put the parameters on 0 time layer into (43) ~ (45) and get the solutions of operation parameters (the pressure P, gas flow rate wgand liquid flow rate wl) at 1 time layer;

  4. Get the solution of liquid holdupHLat 1 time layer by (12);

  5. Repeat step 3 ~ step 4 until reach the giving time layer;

The flow chart of the solution procedures is shown in Fig. 5.

7. Application

7.1. Steady state simulation

Table 1 is the basic data of a condensate gas pipeline. The pressure at starting point is 11.0MPa and the temperature is 330K. The mass of gas flow rate is 29.0kg/s and liquid flow is 0kg/s. The components of the condensate gas are shown in Table 2. Calculate the operation parameters of the pipeline by steady state model and determine the location where phase change occurs.

Absolute roughness
Total diathermanous factor (W/m2.K)Ambient
Step length

Table 1.

Basic data of the condensate gas pipeline

Figure 5.

Solving procedures of transient model


Table 2.

Basic data of the condensate gas pipeline

The solutions can be illustrated as Fig. 6. According to the solutions above, the liquid begins to condensate from the section of 12.3km because the hydrocarbons enter the anti-condensate region.

Figure 6.

The steady state pressure variations of the condensate gas pipeline

Figure 7.

The steady state temperature variations of the condensate gas pipeline

Figure 8.

The steady state gas velocity variations of the condensate gas pipeline

Figure 9.

The steady state liquid velocity variations of the condensate gas pipeline

Figure 10.

The steady state liquid holdup variations of the condensate gas pipeline

The feature of condensate gas pipelines is phase change may occur during operating. This leads to a lot of new phenomena as follow:

  1. It can be seen from Fig.6 that the pressure drop curve of two phase flow is significantly different from of gas flow even the liquid holdup is quite low. The pressure drop of gas flow is non-linear while the appearance of liquid causes a nearly linear curve of the pressure drop. This phenomenon is expressed that the relatively low pressure in the pipeline tends to increase of the gas volume flow; the appearance of condensate liquid and the temperature drop reduce the gas volume flow.

  2. It can be seen from Fig. 7 that the temperature drop curve of two phase flow is similar to single phase flow. The temperature drop gradient of the first half is greater than the last half because of larger temperature difference between the fluid and ambient.

  3. It can be seen from Fig. 8 and Fig.9 that the appearance of two phase flow lead to a reduction of gas flow velocity as well as an increase of liquid flow velocity. The phenomenon also contributes to the nearly linear drop of pressure along the pipeline.

  4. The sharp change of liquid flow velocity as shown in Fig. 9 is caused by phase change. The initial flow velocity of liquid is obtained by flash calculation which makes no consideration of drag force between the phases. Therefore, an abrupt change of the flow rate before and after the phase change occurs as the error made by the flash calculation cannot be ignored. The two-fluid model which has fully considerate of the effect of time is adopted to solve the flow velocity after phase change and the solutions are closer to realistic. It is still a difficulty to improve the accuracy of the initial liquid flow rate at present. The multiple boundaries method is adapted to solve the steady state model. But the astringency and steady state need more improve while this method is applied to non-linear equations.

  5. As shown in Fig.10, the liquid hold up increases behind the phase transition point (two-phase region). Due to the increasing of the liquid hold up is mainly constraint by the phase envelope of the fluid, increasing amount is limited.

The steady state model can simulate the variation of parameters at steady state operation. Actually, there is not absolute steady state condition of the pipeline. If more details of the parameters should be analyzed, following transient simulation method is adopted.

7.2. Transient simulation

Take the previous pipeline as an example, and take the steady state steady parameters as the initial condition of the transient simulation. The boundary condition is set as the pressure at the inlet of pipeline drops to 10.5MPa abruptly at the time of 300s after steady state. The simulation results are shown in Fig11-Fig.15.

Figure 11.

Pressure variation along the pipeline

Figure 12.

Temperature variation along the pipeline

Figure 13.

Velocity of the gas phase variation along the pipeline

Figure 14.

Velocity of the liquid phase variation along the pipeline

Figure 15.

Liquid hold up variation along the pipeline

Compared with steady state, the following features present.

  1. Fig.11 depicts the pressure along the pipeline drops continuously with time elapsing after the inlet pressure drops to 10.5MPa at the time of 300s as the changing of boundary condition.

  2. Fig.12 shows the temperature variation tendency is nearly the same as steady state. The phenomenon can be explained by the reason that the energy equation is ignored in order to simplify the transient model. The approximate method is reasonable because the temperature responses slower than the other parameters.

  3. As depicted in Fig.13, there are abrupt changes of the gas phase velocity at the time of 300s. The opposite direction flow occurs because the pressure at the inlet is lower than the other sections in the pipeline. However, with the rebuilding of the new steady state, the velocity tends to reach a new steady state.

  4. Fig.14 shows the velocity variation along the pipeline. Due to the loss of pressure energy at the inlet, the liquid velocity also drops simultaneously at the time of 300s. Similar to gas velocity, after 300s, the liquid velocity increases gradually and tends to reach new steady state with time elapsing.

  5. Due to the same liquid hold up equation is adopted in the steady state and transient model, the liquid hold up simulated by the transient model and steady state mode has almost the same tendency (Fig.15). However, the liquid hold up increases because of the temperature along the pipeline after 300s is lower than that of initial condition.

Sum up, the more details of the results and transient process can be simulated by transient model. There are still some deficiencies in the model, which should be improved in further work.

8. Conclusions

In this work, a general model for condensate gas pipeline simulation is built on the basis of BWRS EOS, continuity equation, momentum equation, energy equation of the gas and liquid phase. The stratified flow pattern and corresponding constitutive equation are adopted to simplify the model.

By ignoring the parameters variation with time, the steady state simulation model is obtained. To solve the model, the four-order Runge - Kutta method and Gaussian elimination method are used simultaneously. Opposite to steady state model, the transient model is built with consideration of the parameters variation with time, and the model is solved by finite difference method. Solving procedures of steady-state and transient models are presented in detail.

Finally, this work simulated the steady-state and transient operation of a condensate gas pipeline. The pressures, temperatures, velocity of the gas and liquid phase, liquid hold up are calculated. The differences between the steady-state and transient state are discussed. The results show the model and solving method proposed in this work are feasible to simulate the steady state and transient flow in condensate gas pipeline. Nevertheless, in order to expand the adaptive range the models, more improvements should be implemented in future work (Pecenko et al, 2011).


This paper is a project supported by sub-project of National science and technology major project of China (No.2008ZX05054) and China National Petroleum Corporation (CNPC) tackling key subject: Research and Application of Ground Key Technical for CO2 flooding, JW10-W18-J2-11-20.

© 2011 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike-3.0 License, which permits use, distribution and reproduction for non-commercial purposes, provided the original is properly cited and derivative works building on this content are distributed under the same license.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Li Changjun, Jia Wenlong and Wu Xia (September 22nd 2011). Modeling and Simulation for Steady State and Transient Pipe Flow of Condensate Gas, Thermodynamics - Kinetics of Dynamic Systems, Juan Carlos Moreno Piraján, IntechOpen, DOI: 10.5772/23621. Available from:

chapter statistics

4566total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Extended Irreversible Thermodynamics in the Presence of Strong Gravity

By Hiromi Saida

Related Book

First chapter

Thermodynamics of Ligand-Protein Interactions: Implications for Molecular Design

By Agnieszka K. Bronowska

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us