# Numerical Simulation of an Industrial Absorber for Dehydration of Natural Gas Using Triethylene Glycol

Authors: Kenneth Kekpugile Dagde; Jackson Gunorubon Akpa
Journal: Journal of Engineering
 (2014)

---

## Abstract

Models of an absorber for dehydration of natural gas using triethylene glycol are presented. The models were developed by applying the law of conservation of mass and energy to predict the variation of water content of gas and the temperature of the gas and liquid with time along the packing height. The models were integrated numerically using the finite divided difference scheme and incorporated into the MATLAB code. The results obtained agreed reasonably well with industrial plant data obtained from an SPDC TEG unit in Niger-Delta, Nigeria. Model prediction showed a percentage deviation of 8.65% for gas water content and 3.41% and 9.18% for exit temperature of gas and liquid, respectively.


---

## Body

## 1. Introduction

Natural gas needs to be dried before pipeline transport, because the water molecules present in the gas in both vapour and liquid state form hydrates which cause flow restrictions and pressure drops and lower the heating value of gas and corrode pipelines and other equipment. Other problems associated with the presence of water molecules are foaming, degradation, puking, corrosion, low pH, oxidation, thermal decomposition, inadequate absorber design for flow conditions, and salt contamination. Extensive literature is available on common gas dehydration systems including solid and liquid desiccant and refrigeration-based systems [1, 2]. There are several methods of dehydrating natural gas. The most common of these are liquid desiccant (glycol) dehydration and solid desiccant dehydration [3, 4]. Among these gas dehydration processes, absorption is the most common technique, where the water vapor in the gas stream becomes absorbed in a liquid solvent stream. Glycols are the most widely used absorption liquids as they approximate the properties that meet commercial application criteria [5, 6]. Several glycols have been found suitable for commercial application. Triethylene glycol (TEG) is by far the most common liquid desiccant used in natural gas dehydration as it exhibits most of the desirable criteria of commercial suitability [2]. The glycol absorber (contactor) contains trays that provide an adequate intimate contact area between the gas and the glycol. One other option to the tray TEG contactor is the use of structured packing. Structured packing was developed as an alternative to random packing to improve mass transfer control by use of a fixed orientation of the transfer surface. The combination of high gas capacity and reduced height of an equilibrium stage, compared with tray contactors, makes the application of structured packing desirable for both new contactor designs and existing tray contactor capacity upgrades. Hence, the structured packing may offer potential cost savings over trays [1].Optimization of glycol dehydration unit of a natural gas plant is generally aimed at developing a suitable mathematical model which, when tested with plant data, will aid in deciding the best operating conditions required to reduce natural gas water content to the standard pipeline specification of less than 7 lb H2O/MMSCF of gas [7, 8]. Triethylene glycol (TEG) would be used as the absorbent for this process and would be regenerated in a glycol dehydration unit to 99% purity. This, however, is not the case in most of these units. Jaćimović et al. [9] simulated a reactive absorption system for the absorption of CO2 in a packed column using methyl diethanolamine (MDEA) as the solvent. Steady state conditions and plug flow were assumed for the gas phase, leading to a set of ordinary differential equations. In Richardson et al. [10], a mathematical model for the wet scrubbing of CO2 using chilled ammonia was studied. Diffusion and conduction terms were included in the development of the unsteady state models. These models predict the variation of the concentration of the reactants and products with time across the packed height, as well as the variation of the temperature of the system with time across the packed height. The partial differential equations developed were solved using the numerical technique of MATLAB by applying the Robin, Neumann, and Dirichlet boundary conditions (BC) [11]. A similar study on CO2 absorption was carried out by Ahmed et al. [12] using a highly concentrated monoethanolamine (MEA). Most studies on gas dehydration using TEG were simulated using special packages like HYSYS used by Peyghambarzadeh and Jafarpour [13] and the parameters used in their models cannot be easily obtained without extensive experimental studies; thus the model cannot be adapted for simulation of industrial absorber unit. In this paper, models for a functional industrial absorber are presented. The results from the models are compared with data obtained from functional full-scale industrial absorber plant.


## 2. Model Development

The most common method for dehydration in the natural gas industry is the use of a liquid desiccant contactor (absorber) process. In this process, the wet gas is contacted with lean solvent (triethylene glycol) as the absorbent. The water in the gas is absorbed in the lean solvent, producing a rich solvent stream and a dry gas. The dehydrated gas leaves at the top of the column while the glycol leaves at the bottom. Figure1 depicts the hypothetical representation of the dehydrator.Figure 1
Schematic of the dehydrator.The entering wet gas enters the bottom of the absorber and flows up counter currently against the lean triethylene glycol, which enters at the top of the absorber. The triethylene glycol absorbs water vapour from the wet gas as it flows down the column and leaves the bottom of the column rich in water, whereas dry gas leaves from the top of the dehydrator. Therefore, the mass diffusion principles governing this operation will be used in developing the mathematical models for the dehydrator. The models would be developed using the principle of conservation of mass and energy to predict the variation of water content in the gas and the variation of temperature of the gas and triethylene glycol across the height of the dehydrator.


## 3. Model Formulation/Assumptions

The following assumptions were made to develop the model.(i)
Since the column requirement is a diameter ≤ 0.65 m, a packing height of ≤6 m and the fluid are corrosive coupled with a minimum pressure drop across the column, and packed column is preferred to plate column [14, 15].(ii)
The absorber is well lagged; hence, the heat losses are negligible.(iii)
Since the water vapour in the wet gas is the only diffusing component, no diffusing term would be considered for the liquid phase.(iv)
The effect of change in total molar flow rate is ignored, and an average value is assumed constant [16, 17].(v)
Vapour-liquid equilibrium relationship is described using Raoult’s law and Antoine’s equation used for calculation of vapour pressure [18].


### 3.1. Model Development

Material balance (gas phase): Figure2 shows the elemental packed volume and its flow.Figure 2
Elemental packed volume.Consider a homogeneous medium consisting of wet gas (A) and nondiffusive triethylene glycol (B). Let the packed bed be stationary (i.e., the molar average velocity of the mixture is zero), and the mass transfer may occur only by diffusion.Now consider a differential control volumedxdydz.


#### 3.1.1. Mass Balance

A general equation can be derived for a binary mixture of wet gas and nondiffusive triethylene glycol for diffusion and convection that also includes terms for unsteady-state diffusion and chemical reaction. Making the material balance on the wet gas on an element ofdx, dy, and dz fixed in space and shown in Figure 2,
(1)-[∂∂x(NA,x)+∂∂y(NA,y)+∂∂z(NA,z)]=∂CA∂t.For a packed column that is, stationary media, applying Fick’s law (1) reduces to
(2)-[∂∂x(DAB∂CA∂x)+∂∂y(DAB∂CA∂y)+∂∂z(DAB∂CA∂z)]=∂CA∂t.IfDAB is constant, (2) becomes
(3)-[∂2CA∂x2+∂2CA∂y2+∂2CA∂z2]=1DAB∂CA∂t.Since the absorber is in vertical position,(4)-d2CA∂x2=∂2CA∂y2=0.Equation (3) now becomes
(5)∂2CA∂z2=1DAB∂CA∂t.But(6)CA=CAO(1-yA).Differentiating (6),(7a)dCA=-CAOdyA,(7b)d2CA=-CAOd2yA.Substituting (7a) and (7b) into (5) gives
(8)∂yA∂t=DAB∂2yA∂z2.The model equation (8) can be used to predict the variation of water content of gas along the column height at different residence times.


#### 3.1.2. Energy Balance

The energy balance will be carried out using the principle of conservation of energy for both the gas and the liquid triethylene glycol. The glycol enters the column at a higher temperature, transferring some amount of heat to the gas, and hence gas phase energy balance is included.


#### 3.1.3. Energy Balance for the Gas Phase

Figure3 depicts the hypothetical representation of the differential element for energy balance of the gas phase within the packing height, where Tog and Tg are the inlet and outlet temperature of the gas, qz and dqz are the inlet quantity of heat and outlet quantity of heat from the packing space (dz), and dz is the incremental height of the packing space.Figure 3
Hypothetical representation of the energy balance within the differential packing height.Taking cognizance of the conduction of heat axially up the column due to molecular diffusion, the energy balance of the differential element applying the conservation principle gives(9)∂Tg∂t=m˙ACAg∂Tg∂z+Kg∂2TgCAgCpg∂z2+QAdzCAgCpg,
where Cpg is the specific capacity of water vapour in the gas stream, qz is the heat flux in the z-direction due to molecular conduction by Fourier’s law, A is the area of the packing space, and Q is the amount of heat transferred from the lost glycol to the gas steam. The heat transfer at constant pressure is given by Vuthaluru and Bahadori [19] as
(10)Q=LCPLdTL,
where L is the molar flow rate of triethylene glycol in mol/s and CPL and dTL are the heat capacity and temperature difference of the liquid glycol.Recall from dimensional analysis that(11)kgCAgCPg=αg,
where αg is the thermal diffusivity of water vapour in m2/S.Substituting (10) and (11) into (9) results into
(12)∂Tg∂t=-m˙ACAg∂Tg∂z+α∂2Tg∂z2+LCPLACAgCPg∂TL∂z.Letγ=-m˙/ACAg and β=(LCPL/ACAgCPg)(∂TL/∂z), giving
(13)∂Tg∂t=γ∂Tg∂z+α∂2Tg∂z2+β.


#### 3.1.4. Energy Balance for the Liquid Phase

Figure4 depicts the hypothetical representation of the inlet and outlet flow into and out from the differential packing bed in the column.Figure 4
Schematic energy balance for liquid phase.Similarly the energy balance for the liquid phase is made using the principles of conservation of energy taking cognizance that the triethylene glycol flows from the top to the bottom of the column to obtain(14)∂TL∂t=-m˙LACAL∂TL∂z+αL∂TL2∂z2,
where αL=kL/CALCpL is the thermal diffusivity in the liquid phase (triethylene glycol) in m2/s.Equations (8), (13), and (14) constitute the mass balance for the water content in feed gas and the energy balances for gas temperature and TEG temperature variations, respectively, in the absorber.


### 3.2. Operational Parameter and Solution Techniques

#### 3.2.1. Operational Parameters

The input and output operating conditions and the physical properties of the wet gas and glycol (density, molecular weight, and molar volume mass and thermal diffusivity) were estimated from an industrial plant [3, 20] and are presented in Tables 1, 2, and 3.Table 1
Inlet conditions [20].
Components
Input streams
Gas stream
Glycol stream
Weight %
Mol %
Weight %
Mol %
TeG
—
—
99.51
96.054
H2O
0.17
0.187
0.49
3.946
Gas
99.83
99.813
—
—
Total
100.00
100.00
100.00
100.00
Temperature °C
50
55Table 2
Outlet conditions [20].
Components
Output streams
Gas stream
Glycol stream
Weight %
Mol %
Weight %
Mol %
TeG
—
—
95.36
71.127
H2O
0.01
0.011
4.04
28.873
Gas
99.99
99.989
—
—
Total
100.00
100.00
100.00
100.00
Temperature °C
51.3
51Table 3
Physical properties of components [3].
Properties
TEG
H2O
GAS
Molar mass
150.17
18.02
19.83
Molar volume, m3/Kmol
0.01813
—
Mass diffusivity, m2/S
3.80
×
10
-
10
Thermal diffusivity, m2/S
2.338
×
10
-
5
Density, Kg/m3
1125
1000Empirical Evaluation of Mass Diffusivity. The mass diffusivity of water vapour in triethylene glycol (TEG) is evaluated using the formula [3, 21]
(15)D12=1.1728×10-16T(x2M2)1/2μ2V10.6,
where subscript 1 represents the water vapour in the gas, and subscript 2 represents triethylene glycol, where T=50°C = 323.15 K and at T=50°C, μ2=0.01355515 Pa·sec.Solvent Association Parameters.
X
2
=
1 for (TEG), V1=0.0183 m3/Kmol, and M2=150.17 Kg/Kmol.Substitution of these values into (15) gives
(16)D12=3.80×10-10m2/s.Inlet Gas and Glycol Water Content. The inlet gas and glycol water content (in weight %) were obtained from plant operating data and were analytically converted to mol% (assuming binary mixture) using the relations. Consider
(17)y1=x1/m1x1/m1+x2/m2,y2=1-y1,
where x1 and x2 are concentrations of gas and glycol in wt.%, respectively, y1 and y2 are their respective mol%, and m1 and m2 are their molecular weight.


#### 3.2.2. Solution Techniques

A numerical solution based on the finite divided difference scheme was developed and keyed into MATLAB program to solve the condensed models for gas water content, gas temperature, and TEG temperature variations given in (8), (13), and (14), respectively.The developed finite divided difference schemes yield finite grids and computational stencils representing theyA, Tg, and TL, from which Boundary conditions were specified according to “Dirichlet BC.” These boundary conditions and initial conditions are given below.For the gas water content model,(18)yA(zo,t)=yAo,that  is,forz=zo=0,0≤t≤tm,yA(zn,t)=yAf,that  is,forz=zn=H,0≤t≤tm,
where yAo and yAf are initial and final water content in gas stream, respectively.The above boundary conditions explain that the initial gas water content is fixed at the inlet point of absorber column (z=0) and change with varying values of residence time ranging from 0 to tm. More so, the final gas water content is established at the outlet point of the absorber column (z=H) for changing residence time values ranging from 0 to tm.The initial condition is(19)yA(z,to)=0,that  is,fort=to=0,0≤z≤H.This implies that the gas water content is established at only zero residence time for varying absorber column height ranging from bottom to top of column.For the gas temperature model,(20)Tg(zo,t)=Tgi,that  is,forz=zo=O,0≤t≤tm,Tg(zn,t)=Tgo,that  is,forz=zn=H,0≤t≤tm,
where Tgi and Tgo are inlet and outlet gas temperatures, respectively.The above boundary conditions explain that the initial gas temperature is fixed at the inlet point of absorber column (z=0) and changes with varying values of residence time ranging from 0 to tm. More so, the final gas temperature is established at the outlet point of the absorber column (z=H) for changing residence time values ranging from 0 to tm.The initial condition is(21)Tg(z,to)=0,that  is,fort=to=0,0≤z≤H.This implies that the gas temperature is established at only zero residence time for varying absorber column height ranging from bottom to top of column.For the TEG temperature model,(22)TL(zo,t)=TLi,that  is,forz=zo=0,0≤t≤tm,TL(zn,t)=TLo,that  is,forz=zn=H,0≤t≤tm,
where TLi and TLo are inlet and outlet TEG temperatures, respectively.The above boundary conditions explain that the initial TEG temperature is fixed at the inlet point of absorber column (z=0) and changes with varying values of residence time ranging from 0 to tm. More so, the final TEG temperature is established at the outlet point of the absorber column (z=H) for changing residence time values ranging from 0 to tm.The initial condition is(23)TL(z,to)=0,that  is,fort=to=0,0≤z≤H.This implies that the TEG temperature is established at only zero residence time for varying absorber column height ranging from bottom to the top of column.


## 3.1. Model Development

Material balance (gas phase): Figure2 shows the elemental packed volume and its flow.Figure 2
Elemental packed volume.Consider a homogeneous medium consisting of wet gas (A) and nondiffusive triethylene glycol (B). Let the packed bed be stationary (i.e., the molar average velocity of the mixture is zero), and the mass transfer may occur only by diffusion.Now consider a differential control volumedxdydz.


### 3.1.1. Mass Balance

A general equation can be derived for a binary mixture of wet gas and nondiffusive triethylene glycol for diffusion and convection that also includes terms for unsteady-state diffusion and chemical reaction. Making the material balance on the wet gas on an element ofdx, dy, and dz fixed in space and shown in Figure 2,
(1)-[∂∂x(NA,x)+∂∂y(NA,y)+∂∂z(NA,z)]=∂CA∂t.For a packed column that is, stationary media, applying Fick’s law (1) reduces to
(2)-[∂∂x(DAB∂CA∂x)+∂∂y(DAB∂CA∂y)+∂∂z(DAB∂CA∂z)]=∂CA∂t.IfDAB is constant, (2) becomes
(3)-[∂2CA∂x2+∂2CA∂y2+∂2CA∂z2]=1DAB∂CA∂t.Since the absorber is in vertical position,(4)-d2CA∂x2=∂2CA∂y2=0.Equation (3) now becomes
(5)∂2CA∂z2=1DAB∂CA∂t.But(6)CA=CAO(1-yA).Differentiating (6),(7a)dCA=-CAOdyA,(7b)d2CA=-CAOd2yA.Substituting (7a) and (7b) into (5) gives
(8)∂yA∂t=DAB∂2yA∂z2.The model equation (8) can be used to predict the variation of water content of gas along the column height at different residence times.


### 3.1.2. Energy Balance

The energy balance will be carried out using the principle of conservation of energy for both the gas and the liquid triethylene glycol. The glycol enters the column at a higher temperature, transferring some amount of heat to the gas, and hence gas phase energy balance is included.


### 3.1.3. Energy Balance for the Gas Phase

Figure3 depicts the hypothetical representation of the differential element for energy balance of the gas phase within the packing height, where Tog and Tg are the inlet and outlet temperature of the gas, qz and dqz are the inlet quantity of heat and outlet quantity of heat from the packing space (dz), and dz is the incremental height of the packing space.Figure 3
Hypothetical representation of the energy balance within the differential packing height.Taking cognizance of the conduction of heat axially up the column due to molecular diffusion, the energy balance of the differential element applying the conservation principle gives(9)∂Tg∂t=m˙ACAg∂Tg∂z+Kg∂2TgCAgCpg∂z2+QAdzCAgCpg,
where Cpg is the specific capacity of water vapour in the gas stream, qz is the heat flux in the z-direction due to molecular conduction by Fourier’s law, A is the area of the packing space, and Q is the amount of heat transferred from the lost glycol to the gas steam. The heat transfer at constant pressure is given by Vuthaluru and Bahadori [19] as
(10)Q=LCPLdTL,
where L is the molar flow rate of triethylene glycol in mol/s and CPL and dTL are the heat capacity and temperature difference of the liquid glycol.Recall from dimensional analysis that(11)kgCAgCPg=αg,
where αg is the thermal diffusivity of water vapour in m2/S.Substituting (10) and (11) into (9) results into
(12)∂Tg∂t=-m˙ACAg∂Tg∂z+α∂2Tg∂z2+LCPLACAgCPg∂TL∂z.Letγ=-m˙/ACAg and β=(LCPL/ACAgCPg)(∂TL/∂z), giving
(13)∂Tg∂t=γ∂Tg∂z+α∂2Tg∂z2+β.


### 3.1.4. Energy Balance for the Liquid Phase

Figure4 depicts the hypothetical representation of the inlet and outlet flow into and out from the differential packing bed in the column.Figure 4
Schematic energy balance for liquid phase.Similarly the energy balance for the liquid phase is made using the principles of conservation of energy taking cognizance that the triethylene glycol flows from the top to the bottom of the column to obtain(14)∂TL∂t=-m˙LACAL∂TL∂z+αL∂TL2∂z2,
where αL=kL/CALCpL is the thermal diffusivity in the liquid phase (triethylene glycol) in m2/s.Equations (8), (13), and (14) constitute the mass balance for the water content in feed gas and the energy balances for gas temperature and TEG temperature variations, respectively, in the absorber.


## 3.1.1. Mass Balance

A general equation can be derived for a binary mixture of wet gas and nondiffusive triethylene glycol for diffusion and convection that also includes terms for unsteady-state diffusion and chemical reaction. Making the material balance on the wet gas on an element ofdx, dy, and dz fixed in space and shown in Figure 2,
(1)-[∂∂x(NA,x)+∂∂y(NA,y)+∂∂z(NA,z)]=∂CA∂t.For a packed column that is, stationary media, applying Fick’s law (1) reduces to
(2)-[∂∂x(DAB∂CA∂x)+∂∂y(DAB∂CA∂y)+∂∂z(DAB∂CA∂z)]=∂CA∂t.IfDAB is constant, (2) becomes
(3)-[∂2CA∂x2+∂2CA∂y2+∂2CA∂z2]=1DAB∂CA∂t.Since the absorber is in vertical position,(4)-d2CA∂x2=∂2CA∂y2=0.Equation (3) now becomes
(5)∂2CA∂z2=1DAB∂CA∂t.But(6)CA=CAO(1-yA).Differentiating (6),(7a)dCA=-CAOdyA,(7b)d2CA=-CAOd2yA.Substituting (7a) and (7b) into (5) gives
(8)∂yA∂t=DAB∂2yA∂z2.The model equation (8) can be used to predict the variation of water content of gas along the column height at different residence times.


## 3.1.2. Energy Balance

The energy balance will be carried out using the principle of conservation of energy for both the gas and the liquid triethylene glycol. The glycol enters the column at a higher temperature, transferring some amount of heat to the gas, and hence gas phase energy balance is included.


## 3.1.3. Energy Balance for the Gas Phase

Figure3 depicts the hypothetical representation of the differential element for energy balance of the gas phase within the packing height, where Tog and Tg are the inlet and outlet temperature of the gas, qz and dqz are the inlet quantity of heat and outlet quantity of heat from the packing space (dz), and dz is the incremental height of the packing space.Figure 3
Hypothetical representation of the energy balance within the differential packing height.Taking cognizance of the conduction of heat axially up the column due to molecular diffusion, the energy balance of the differential element applying the conservation principle gives(9)∂Tg∂t=m˙ACAg∂Tg∂z+Kg∂2TgCAgCpg∂z2+QAdzCAgCpg,
where Cpg is the specific capacity of water vapour in the gas stream, qz is the heat flux in the z-direction due to molecular conduction by Fourier’s law, A is the area of the packing space, and Q is the amount of heat transferred from the lost glycol to the gas steam. The heat transfer at constant pressure is given by Vuthaluru and Bahadori [19] as
(10)Q=LCPLdTL,
where L is the molar flow rate of triethylene glycol in mol/s and CPL and dTL are the heat capacity and temperature difference of the liquid glycol.Recall from dimensional analysis that(11)kgCAgCPg=αg,
where αg is the thermal diffusivity of water vapour in m2/S.Substituting (10) and (11) into (9) results into
(12)∂Tg∂t=-m˙ACAg∂Tg∂z+α∂2Tg∂z2+LCPLACAgCPg∂TL∂z.Letγ=-m˙/ACAg and β=(LCPL/ACAgCPg)(∂TL/∂z), giving
(13)∂Tg∂t=γ∂Tg∂z+α∂2Tg∂z2+β.


## 3.1.4. Energy Balance for the Liquid Phase

Figure4 depicts the hypothetical representation of the inlet and outlet flow into and out from the differential packing bed in the column.Figure 4
Schematic energy balance for liquid phase.Similarly the energy balance for the liquid phase is made using the principles of conservation of energy taking cognizance that the triethylene glycol flows from the top to the bottom of the column to obtain(14)∂TL∂t=-m˙LACAL∂TL∂z+αL∂TL2∂z2,
where αL=kL/CALCpL is the thermal diffusivity in the liquid phase (triethylene glycol) in m2/s.Equations (8), (13), and (14) constitute the mass balance for the water content in feed gas and the energy balances for gas temperature and TEG temperature variations, respectively, in the absorber.


## 3.2. Operational Parameter and Solution Techniques

### 3.2.1. Operational Parameters

The input and output operating conditions and the physical properties of the wet gas and glycol (density, molecular weight, and molar volume mass and thermal diffusivity) were estimated from an industrial plant [3, 20] and are presented in Tables 1, 2, and 3.Table 1
Inlet conditions [20].
Components
Input streams
Gas stream
Glycol stream
Weight %
Mol %
Weight %
Mol %
TeG
—
—
99.51
96.054
H2O
0.17
0.187
0.49
3.946
Gas
99.83
99.813
—
—
Total
100.00
100.00
100.00
100.00
Temperature °C
50
55Table 2
Outlet conditions [20].
Components
Output streams
Gas stream
Glycol stream
Weight %
Mol %
Weight %
Mol %
TeG
—
—
95.36
71.127
H2O
0.01
0.011
4.04
28.873
Gas
99.99
99.989
—
—
Total
100.00
100.00
100.00
100.00
Temperature °C
51.3
51Table 3
Physical properties of components [3].
Properties
TEG
H2O
GAS
Molar mass
150.17
18.02
19.83
Molar volume, m3/Kmol
0.01813
—
Mass diffusivity, m2/S
3.80
×
10
-
10
Thermal diffusivity, m2/S
2.338
×
10
-
5
Density, Kg/m3
1125
1000Empirical Evaluation of Mass Diffusivity. The mass diffusivity of water vapour in triethylene glycol (TEG) is evaluated using the formula [3, 21]
(15)D12=1.1728×10-16T(x2M2)1/2μ2V10.6,
where subscript 1 represents the water vapour in the gas, and subscript 2 represents triethylene glycol, where T=50°C = 323.15 K and at T=50°C, μ2=0.01355515 Pa·sec.Solvent Association Parameters.
X
2
=
1 for (TEG), V1=0.0183 m3/Kmol, and M2=150.17 Kg/Kmol.Substitution of these values into (15) gives
(16)D12=3.80×10-10m2/s.Inlet Gas and Glycol Water Content. The inlet gas and glycol water content (in weight %) were obtained from plant operating data and were analytically converted to mol% (assuming binary mixture) using the relations. Consider
(17)y1=x1/m1x1/m1+x2/m2,y2=1-y1,
where x1 and x2 are concentrations of gas and glycol in wt.%, respectively, y1 and y2 are their respective mol%, and m1 and m2 are their molecular weight.


### 3.2.2. Solution Techniques

A numerical solution based on the finite divided difference scheme was developed and keyed into MATLAB program to solve the condensed models for gas water content, gas temperature, and TEG temperature variations given in (8), (13), and (14), respectively.The developed finite divided difference schemes yield finite grids and computational stencils representing theyA, Tg, and TL, from which Boundary conditions were specified according to “Dirichlet BC.” These boundary conditions and initial conditions are given below.For the gas water content model,(18)yA(zo,t)=yAo,that  is,forz=zo=0,0≤t≤tm,yA(zn,t)=yAf,that  is,forz=zn=H,0≤t≤tm,
where yAo and yAf are initial and final water content in gas stream, respectively.The above boundary conditions explain that the initial gas water content is fixed at the inlet point of absorber column (z=0) and change with varying values of residence time ranging from 0 to tm. More so, the final gas water content is established at the outlet point of the absorber column (z=H) for changing residence time values ranging from 0 to tm.The initial condition is(19)yA(z,to)=0,that  is,fort=to=0,0≤z≤H.This implies that the gas water content is established at only zero residence time for varying absorber column height ranging from bottom to top of column.For the gas temperature model,(20)Tg(zo,t)=Tgi,that  is,forz=zo=O,0≤t≤tm,Tg(zn,t)=Tgo,that  is,forz=zn=H,0≤t≤tm,
where Tgi and Tgo are inlet and outlet gas temperatures, respectively.The above boundary conditions explain that the initial gas temperature is fixed at the inlet point of absorber column (z=0) and changes with varying values of residence time ranging from 0 to tm. More so, the final gas temperature is established at the outlet point of the absorber column (z=H) for changing residence time values ranging from 0 to tm.The initial condition is(21)Tg(z,to)=0,that  is,fort=to=0,0≤z≤H.This implies that the gas temperature is established at only zero residence time for varying absorber column height ranging from bottom to top of column.For the TEG temperature model,(22)TL(zo,t)=TLi,that  is,forz=zo=0,0≤t≤tm,TL(zn,t)=TLo,that  is,forz=zn=H,0≤t≤tm,
where TLi and TLo are inlet and outlet TEG temperatures, respectively.The above boundary conditions explain that the initial TEG temperature is fixed at the inlet point of absorber column (z=0) and changes with varying values of residence time ranging from 0 to tm. More so, the final TEG temperature is established at the outlet point of the absorber column (z=H) for changing residence time values ranging from 0 to tm.The initial condition is(23)TL(z,to)=0,that  is,fort=to=0,0≤z≤H.This implies that the TEG temperature is established at only zero residence time for varying absorber column height ranging from bottom to the top of column.


## 3.2.1. Operational Parameters

The input and output operating conditions and the physical properties of the wet gas and glycol (density, molecular weight, and molar volume mass and thermal diffusivity) were estimated from an industrial plant [3, 20] and are presented in Tables 1, 2, and 3.Table 1
Inlet conditions [20].
Components
Input streams
Gas stream
Glycol stream
Weight %
Mol %
Weight %
Mol %
TeG
—
—
99.51
96.054
H2O
0.17
0.187
0.49
3.946
Gas
99.83
99.813
—
—
Total
100.00
100.00
100.00
100.00
Temperature °C
50
55Table 2
Outlet conditions [20].
Components
Output streams
Gas stream
Glycol stream
Weight %
Mol %
Weight %
Mol %
TeG
—
—
95.36
71.127
H2O
0.01
0.011
4.04
28.873
Gas
99.99
99.989
—
—
Total
100.00
100.00
100.00
100.00
Temperature °C
51.3
51Table 3
Physical properties of components [3].
Properties
TEG
H2O
GAS
Molar mass
150.17
18.02
19.83
Molar volume, m3/Kmol
0.01813
—
Mass diffusivity, m2/S
3.80
×
10
-
10
Thermal diffusivity, m2/S
2.338
×
10
-
5
Density, Kg/m3
1125
1000Empirical Evaluation of Mass Diffusivity. The mass diffusivity of water vapour in triethylene glycol (TEG) is evaluated using the formula [3, 21]
(15)D12=1.1728×10-16T(x2M2)1/2μ2V10.6,
where subscript 1 represents the water vapour in the gas, and subscript 2 represents triethylene glycol, where T=50°C = 323.15 K and at T=50°C, μ2=0.01355515 Pa·sec.Solvent Association Parameters.
X
2
=
1 for (TEG), V1=0.0183 m3/Kmol, and M2=150.17 Kg/Kmol.Substitution of these values into (15) gives
(16)D12=3.80×10-10m2/s.Inlet Gas and Glycol Water Content. The inlet gas and glycol water content (in weight %) were obtained from plant operating data and were analytically converted to mol% (assuming binary mixture) using the relations. Consider
(17)y1=x1/m1x1/m1+x2/m2,y2=1-y1,
where x1 and x2 are concentrations of gas and glycol in wt.%, respectively, y1 and y2 are their respective mol%, and m1 and m2 are their molecular weight.


## 3.2.2. Solution Techniques

A numerical solution based on the finite divided difference scheme was developed and keyed into MATLAB program to solve the condensed models for gas water content, gas temperature, and TEG temperature variations given in (8), (13), and (14), respectively.The developed finite divided difference schemes yield finite grids and computational stencils representing theyA, Tg, and TL, from which Boundary conditions were specified according to “Dirichlet BC.” These boundary conditions and initial conditions are given below.For the gas water content model,(18)yA(zo,t)=yAo,that  is,forz=zo=0,0≤t≤tm,yA(zn,t)=yAf,that  is,forz=zn=H,0≤t≤tm,
where yAo and yAf are initial and final water content in gas stream, respectively.The above boundary conditions explain that the initial gas water content is fixed at the inlet point of absorber column (z=0) and change with varying values of residence time ranging from 0 to tm. More so, the final gas water content is established at the outlet point of the absorber column (z=H) for changing residence time values ranging from 0 to tm.The initial condition is(19)yA(z,to)=0,that  is,fort=to=0,0≤z≤H.This implies that the gas water content is established at only zero residence time for varying absorber column height ranging from bottom to top of column.For the gas temperature model,(20)Tg(zo,t)=Tgi,that  is,forz=zo=O,0≤t≤tm,Tg(zn,t)=Tgo,that  is,forz=zn=H,0≤t≤tm,
where Tgi and Tgo are inlet and outlet gas temperatures, respectively.The above boundary conditions explain that the initial gas temperature is fixed at the inlet point of absorber column (z=0) and changes with varying values of residence time ranging from 0 to tm. More so, the final gas temperature is established at the outlet point of the absorber column (z=H) for changing residence time values ranging from 0 to tm.The initial condition is(21)Tg(z,to)=0,that  is,fort=to=0,0≤z≤H.This implies that the gas temperature is established at only zero residence time for varying absorber column height ranging from bottom to top of column.For the TEG temperature model,(22)TL(zo,t)=TLi,that  is,forz=zo=0,0≤t≤tm,TL(zn,t)=TLo,that  is,forz=zn=H,0≤t≤tm,
where TLi and TLo are inlet and outlet TEG temperatures, respectively.The above boundary conditions explain that the initial TEG temperature is fixed at the inlet point of absorber column (z=0) and changes with varying values of residence time ranging from 0 to tm. More so, the final TEG temperature is established at the outlet point of the absorber column (z=H) for changing residence time values ranging from 0 to tm.The initial condition is(23)TL(z,to)=0,that  is,fort=to=0,0≤z≤H.This implies that the TEG temperature is established at only zero residence time for varying absorber column height ranging from bottom to the top of column.


## 4. Results and Discussion

Table4 shows the comparison between plant data and predictions from model (see (8), (13), and (14)), indicating that the predicted results agree reasonably well with the plant data. These results show a deviation ranging from 3.41 to 9.18 percent.Table 4
Comparison between plant data and model predictions.
Process parameter
Model prediction
Plant data
% deviation
Final gas water content
7.93
×
10
-
7
7.24
×
10
-
7
8.65
Gas outlet temperature (°C)
44.52
43
3.41
TEG outlet temperature (°C)
48.45
44
9.18Profiles presented and discussed here will subsequently reveal the following: variations of gas water content with time and axial height of packing in the column, variation of temperature of triethylene glycol (TEG) with column height at different thermal diffusivities, variation of temperature of gas with column height at different residence times, variation of gas water content across column height at different mass diffusivities, and variation of temperature of triethylene glycol across the column height at different residence times.


### 4.1. Variation of Water Content of Gas with Column Height at Different Residence Times

It can be deduced from Figure5 that the water content of the gas reduces as the gas moves from the bottom of the column to the top. It can also be deduced that the higher the residence time of the gas in the column, the higher the rate of removal of the water vapour from the gas. This holds true since a relatively smaller time is needed to establish equilibrium between the water vapour in the gas and that in the liquid phase [10, 22]. This means that as the residence time increases, say, to 200 seconds, the water vapour returns to the vapour phase again implying that the water content in gas increases. It can also be deduced from Figure 5 that, at a height of approximately 7 m and above, the gas water content variation becomes steady.Figure 5
Variation of gas water content (mole fraction) from bottom of column.In addition, the solutions of the model will be represented as a three-dimensional surface plot in Figure6. The purpose of the surface plots is to visualize the propagation of the gas water content in time and space and to make conclusions based on the subsequent trends. The surface plots are not intended to give the exact numerical values but for visualization.Figure 6
Surface plot showing gas water content propagation along the column.The natural gas propagates from the base of the absorber and initially holds a water concentration of 0.187 mol%. The low resistance in the gas bulk will cause the gas and liquid bulk phases to reach chemical equilibrium virtually instantaneous. The steep transient observed at the lower part of the column confirms the trend illustrated by Figure6; it is also in agreement with plant data. Also, as operation proceeds half way up the column, the absorption of water from natural gas becomes numerically insignificant and remains practically constant.


### 4.2. Variation of Temperature of Triethylene Glycol with the Column Height at Different Residence Times

In Figure7, the temperature of the absorbing solvent, triethylene glycol (TEG), reduces gradually as it travels from the top of the column to the bottom. Initially, triethylene glycol enters the column at a temperature of 50°C and leaves at 46.5–48°C depending on the residence time. It is observed that the temperature change becomes smaller as the residence time increases, resulting into a very steep slope at time = 200 seconds. This is obviously because more water vapour has been absorbed at higher residence time.Figure 7
Variation of temperature of triethylene glycol down the column.The surface plot in Figure8 further visualizes the TEG temperature reduction at varying residence times. The largest TEG temperature change is achieved at a residence time of 40 to 80 seconds, while a little change is achieved at 160 to 200 seconds.Figure 8
Surface plot showing triethylene glycol temperature propagation along the column.


### 4.3. Variation of Temperature of Gas with Column Height at Different Residence Times

The gas enters the column at a temperature of 42.5°C and increases very slowly until it leaves the column at a slightly higher temperature of 46.5°C. However, at a residence timet=40 seconds and at column height z=12 m, the outlet temperature of gas is approximately 44.5°C. It can be noticed in Figure 9 that the temperature becomes lower at higher residence time. This implies that if the residence time is reduced further, the required outlet gas temperature would be achieved at approximately 120 seconds. It should be noted that the increase in residence time results from the transfer of heat from the liquid stream to the gas stream.Figure 9
Variation of the temperature of the gas stream with column height from column bottom.The surface plot in Figure10 reveals that the gas temperature variation is not widely distributed. It can be further deduced that the lowest trend in gas temperature is attained at residence time of 120 seconds, while the largest change is obtained at residence time of approximately 40 to 80 seconds.Figure 10
Surface plot showing gas temperature variation along the column.


### 4.4. Variation of Water Content of Gas with Column Height at Different Mass Diffusivities

The mass diffusivity is the property of a material that determines the rate at which a given component is transferred across a concentration gradient. This property is a vital parameter in this work. From Figure11, it is evident that, given a fixed time of 40 seconds, at higher mass diffusivities, the rate of transfer of water vapour from gas to the liquid stream decreases slightly as up the column. Also, as the mass diffusivity reduces, the rate of transfer of water vapour from gas to the triethylene glycol stream increases sharply. This implies that the mass diffusion coefficient of the gas through the TEG should be as low as ≤3.80 × 10−10 m2/s for optimal absorption.Figure 11
Variation of water content of gas at different mass diffusivities att=40 seconds.The surface plot in Figure12 shows that the gas water content variation widely spreads across the column height at different mass diffusivities. Similar final gas water content values are obtained at mass diffusivities of 3.80 × 10−9 m2/s and 3.80 × 10−10 m2/s.Figure 12
Surface plot showing gas water content variation at different mass diffusivities when time = 40 seconds.


### 4.5. Variation of Temperature of Triethylene Glycol (TEG) with Column Height at Different Thermal Diffusivities

Thermal diffusivity is the property of a material which describes the rate at which heat flows through the material. Water vapour, being a better heat transfer agent, has a value of 2.336 × 10−5 m2/s, whereas liquid water has a value of 1.4 × 10−5 m2/s. From Figure 13, it is observed generally that the temperature of the solvent triethylene glycol decreases sharply down the column as the thermal diffusivity decreases. This decrease is faster as the thermal diffusivity decreases, resulting in a steep slope at the lowest thermal diffusivity of 2.338 × 10−5 m2/s. It is imperative to note that these plots were taken at a residence time of t=40 seconds.Figure 13
Variation of temperature of TEG with thermal diffusivities att=40 seconds.In addition, the surface plot in Figure14 clearly visualizes the TEG temperature decrease along the column and at varying thermal diffusivities. A relatively large TEG temperature change exists between thermal diffusivities of 2.338 × 10−4 and 2.338 × 10−5 m2/s; conversely, there is a insignificant temperature change at 2.338 × 10−1 m2/s.Figure 14
Surface plot of TEG temperature at different thermal diffusivities att=40 seconds.


## 4.1. Variation of Water Content of Gas with Column Height at Different Residence Times

It can be deduced from Figure5 that the water content of the gas reduces as the gas moves from the bottom of the column to the top. It can also be deduced that the higher the residence time of the gas in the column, the higher the rate of removal of the water vapour from the gas. This holds true since a relatively smaller time is needed to establish equilibrium between the water vapour in the gas and that in the liquid phase [10, 22]. This means that as the residence time increases, say, to 200 seconds, the water vapour returns to the vapour phase again implying that the water content in gas increases. It can also be deduced from Figure 5 that, at a height of approximately 7 m and above, the gas water content variation becomes steady.Figure 5
Variation of gas water content (mole fraction) from bottom of column.In addition, the solutions of the model will be represented as a three-dimensional surface plot in Figure6. The purpose of the surface plots is to visualize the propagation of the gas water content in time and space and to make conclusions based on the subsequent trends. The surface plots are not intended to give the exact numerical values but for visualization.Figure 6
Surface plot showing gas water content propagation along the column.The natural gas propagates from the base of the absorber and initially holds a water concentration of 0.187 mol%. The low resistance in the gas bulk will cause the gas and liquid bulk phases to reach chemical equilibrium virtually instantaneous. The steep transient observed at the lower part of the column confirms the trend illustrated by Figure6; it is also in agreement with plant data. Also, as operation proceeds half way up the column, the absorption of water from natural gas becomes numerically insignificant and remains practically constant.


## 4.2. Variation of Temperature of Triethylene Glycol with the Column Height at Different Residence Times

In Figure7, the temperature of the absorbing solvent, triethylene glycol (TEG), reduces gradually as it travels from the top of the column to the bottom. Initially, triethylene glycol enters the column at a temperature of 50°C and leaves at 46.5–48°C depending on the residence time. It is observed that the temperature change becomes smaller as the residence time increases, resulting into a very steep slope at time = 200 seconds. This is obviously because more water vapour has been absorbed at higher residence time.Figure 7
Variation of temperature of triethylene glycol down the column.The surface plot in Figure8 further visualizes the TEG temperature reduction at varying residence times. The largest TEG temperature change is achieved at a residence time of 40 to 80 seconds, while a little change is achieved at 160 to 200 seconds.Figure 8
Surface plot showing triethylene glycol temperature propagation along the column.


## 4.3. Variation of Temperature of Gas with Column Height at Different Residence Times

The gas enters the column at a temperature of 42.5°C and increases very slowly until it leaves the column at a slightly higher temperature of 46.5°C. However, at a residence timet=40 seconds and at column height z=12 m, the outlet temperature of gas is approximately 44.5°C. It can be noticed in Figure 9 that the temperature becomes lower at higher residence time. This implies that if the residence time is reduced further, the required outlet gas temperature would be achieved at approximately 120 seconds. It should be noted that the increase in residence time results from the transfer of heat from the liquid stream to the gas stream.Figure 9
Variation of the temperature of the gas stream with column height from column bottom.The surface plot in Figure10 reveals that the gas temperature variation is not widely distributed. It can be further deduced that the lowest trend in gas temperature is attained at residence time of 120 seconds, while the largest change is obtained at residence time of approximately 40 to 80 seconds.Figure 10
Surface plot showing gas temperature variation along the column.


## 4.4. Variation of Water Content of Gas with Column Height at Different Mass Diffusivities

The mass diffusivity is the property of a material that determines the rate at which a given component is transferred across a concentration gradient. This property is a vital parameter in this work. From Figure11, it is evident that, given a fixed time of 40 seconds, at higher mass diffusivities, the rate of transfer of water vapour from gas to the liquid stream decreases slightly as up the column. Also, as the mass diffusivity reduces, the rate of transfer of water vapour from gas to the triethylene glycol stream increases sharply. This implies that the mass diffusion coefficient of the gas through the TEG should be as low as ≤3.80 × 10−10 m2/s for optimal absorption.Figure 11
Variation of water content of gas at different mass diffusivities att=40 seconds.The surface plot in Figure12 shows that the gas water content variation widely spreads across the column height at different mass diffusivities. Similar final gas water content values are obtained at mass diffusivities of 3.80 × 10−9 m2/s and 3.80 × 10−10 m2/s.Figure 12
Surface plot showing gas water content variation at different mass diffusivities when time = 40 seconds.


## 4.5. Variation of Temperature of Triethylene Glycol (TEG) with Column Height at Different Thermal Diffusivities

Thermal diffusivity is the property of a material which describes the rate at which heat flows through the material. Water vapour, being a better heat transfer agent, has a value of 2.336 × 10−5 m2/s, whereas liquid water has a value of 1.4 × 10−5 m2/s. From Figure 13, it is observed generally that the temperature of the solvent triethylene glycol decreases sharply down the column as the thermal diffusivity decreases. This decrease is faster as the thermal diffusivity decreases, resulting in a steep slope at the lowest thermal diffusivity of 2.338 × 10−5 m2/s. It is imperative to note that these plots were taken at a residence time of t=40 seconds.Figure 13
Variation of temperature of TEG with thermal diffusivities att=40 seconds.In addition, the surface plot in Figure14 clearly visualizes the TEG temperature decrease along the column and at varying thermal diffusivities. A relatively large TEG temperature change exists between thermal diffusivities of 2.338 × 10−4 and 2.338 × 10−5 m2/s; conversely, there is a insignificant temperature change at 2.338 × 10−1 m2/s.Figure 14
Surface plot of TEG temperature at different thermal diffusivities att=40 seconds.


## 5. Conclusion

Mathematical models of the absorber of a glycol dehydration facility were developed using the principles of conservation of mass and energy. The models could predict the variation of the water content of gas in mole fraction and the gas and liquid (TEG) temperatures across the packing height. The models developed contain contributions from bulk and diffusion flows. The models were validated using the initial conditions from a functional industrial TEG unit in Nigeria to ascertain if the outlet conditions predicted by the models meet the industrial plant outlet values.



---

*Source: 693902-2014-07-20.xml*