Published Paper
Published Paper
net/publication/259996289
CITATIONS READS
27 1,268
3 authors:
S.M. Peyghambarzadeh
Islamic Azad University
91 PUBLICATIONS 2,469 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
Flare Gas Recovery project at the gas processing plant View project
All content following this page was uploaded by Mohammadreza Bohloul on 08 October 2017.
a r t i c l e i n f o a b s t r a c t
Article history: In this study, some new experimental data on the solubility of CO2 in N-methyl-2-pyrrolidone (NMP)
Received 16 October 2013 is reported at different temperatures of 293 K, 313 K and 333 K and different pressures 0.5–1.1 MPa. The
Received in revised form solubility of CO2 was evaluated by pressure decay method in a batch stirred absorption cell. The solubility
19 December 2013
data was also envisaged by using the thermodynamic modeling of the problem with 1-parameter and
Accepted 23 December 2013
2-parameter versions of Peng–Robinson equation of state (PR-EOS). Using the experimental data, the
Available online 8 January 2014
values of the binary interaction parameters (kij and lij ) were optimized. The mean deviation between
the measured solubility data and the predicted values by the thermodynamic model was calculated to
Keywords:
CO2 solubility
be less than 3.5%. The results from both versions of PR-EOS are significantly similar. Therefore, it can
N-methyl-2-pyrrolidone (NMP) be concluded that the application of 1-parameter model for NMP–CO2 system is simpler approach with
Vapor–liquid equilibrium (VLE) sufficient accuracy to represent the vapor–liquid equilibria.
Absorption © 2013 Elsevier B.V. All rights reserved.
PR-EOS
0378-3812/$ – see front matter © 2013 Elsevier B.V. All rights reserved.
http://dx.doi.org/10.1016/j.fluid.2013.12.019
M.R. Bohloul et al. / Fluid Phase Equilibria 365 (2014) 106–111 107
Table 1
The specifications and source of the materials used in this study.
Name of compound Source Final purity Purity measuring method Purification method
NMP Bandar Imam Higher than 99.5 wt% Purities were obtained from Materials were used without
Petrochemical suppliers any further purification
Company (BIPC), Iran
CO2 Balon Gas Company, 99.9 mol%
Iran
[13] presented the experimental data on the solubility of CO2 in 2.2. Apparatus and procedure
amine salt solution. Moreover, they illustrated the phase equilibria
of CO2 –H2 O–NaCl system using PR-EOS. The equipment used in this study is shown illustratively in
Georgia et al. [14] studied the solubility of CO2 in the aque- Fig. 1. This type of batch absorption cell was extensively utilized
ous phase. They investigated the solubility of H2 O in the CO2 rich in the previous investigations as well (e.g. Young et al. [17] and
phase through the PR-EOS coupled with the conventional van der Malyanah et al. [7]). The apparatus contains gas reservoir, reg-
Waals one fluid mixing rules, the PR-EOS coupled with the universal ulator, middle cell, absorption cell, vacuum pump, pressure and
mixing rules (UMR), and the CPA-PR EOS. temperature measuring instruments, and indicators. The middle
Rafel et al. [15] resulted in phase equilibrium data at cell and absorption cell were made up of 316 stainless steel with
high pressures for the binary and ternary systems formed by internal volume of approximately 500 cm3 and 320 cm3 , respec-
propane + NMP + methanol. Consequently, thermodynamic mod- tively. The middle cell was employed to set the gas temperature
eling was executed by means of experimental data along with for the gas absorption process. The temperature of the CO2 gas in
PR-EOS using the conventional quadratic van der Waals mix- the middle cell and absorption cell was measured with a Pt-100
ing rule (vdW2). Yizhuan and Chan [16] accomplished the thermo-resistance. The pressure of the gas in the middle cell and
thermodynamic modeling of CO2 solubility in aqueous solu- absorption cell was constantly measured using a pressure sensor
tion of NaCl and Na2 SO4 through the electrolyte NRTL activity (Sensys, Model: PSCH0025BCIJ). The middle cell and the absorp-
model. tion cell were situated in a water bath, and the temperature of the
In this paper, experimental data on the solubility of CO2 in water bath was sustained the steady temperature with the pre-
NMP is presented at different temperatures of 293 K, 313 K, and cision of ±0.2 K. A magnetic stirrer is used in the absorption cell
333 K and different pressures up to 1.1 MPa. It must be empha- to attain the equilibrium condition faster. The kinetic data highly
sized that scarce experimental data on the solubility of CO2 depends on the experimental conditions such as speed and size of
in NMP have been presented up to now. These data are then stirrer. The purpose of using a stirrer in all of the experiments is
correlated using both 1-parameter and 2-parameter versions of to increase the mass transfer rate and reduction of the equilibrium
PR-EOS. Furthermore, the values of binary interaction parameters time. In all of the experiments, stirring was performed continuously
(kij and lij ) are also optimized using the collected experimental in the absorption cell at a speed of 300 rpm, maintaining a smooth
data. interfacial area.
In this study, experimental data on the solubility of CO2 in NMP
was collected at temperatures of 293 K, 313 K, and 333 K. The pres-
sure of the gas from CO2 reservoir is attuned using a regulator.
2. Experimental Subsequently, the gas is sent into the middle cell to adjust the
temperature; ultimately, it enters the absorption cell and the gas
2.1. Materials absorption process begins. It is important to mention that the inlet
CO2 is heated prior to entering the middle cell in order to pre-
N-methyl-2-pyrrolidone (NMP) with the purity higher than vent the gas freezing due to sudden pressure drop in the regulator.
99.5 wt%, and carbon dioxide with the purity of 99.9 mol% were Moreover, 100 ml of NMP is kept in the absorption cell. Prior to per-
used in this research. Table 1 shows the specifications and source forming the experiment, the middle cell, piping and absorption cell
of these materials. are evacuated by a vacuum pump (JB Industries DV-200N, USA). The
PCO2 vs. t data is documented until the gas–liquid equilibrium state Table 2
Summary of the measurement errors.
attains. The solubility of CO2 is evaluated by measuring the differ-
ence between the initial and final pressure (equilibrium point) in Parameter Accuracy of the Contribution in the
the absorption cell. The major feature of this experimental proce- measurement uncertainty of x.
device
dure is that no analysis of the liquid phase is necessitated. The real
gas law is applied to determine the number of absorbed moles of (Pi –Pf ) ±0.01 bar 37.7%
CO2 : VG ±1 ml 15.3%
VL ±1 ml 27.9%
Pi .Vi Pf .Vf T ±0.2 ◦ C 19.1%
nCO2 = − (1)
R.T.Zi R.T.Zf
2.3. Uncertainty analysis Pressure decay during absorption of CO2 in NMP at three dif-
ferent temperatures was presented in Figs. 2–4. It is imperative
Uncertainty analysis is carried out by calculating the error of the to declare that the pressure reduction in the absorption cell would
measurements. The uncertainty of the equilibrium concentration present the amount of gas absorption, and this would continue until
comes from the errors in the measurement of gas and liquid vol- the solvent (NMP) becomes saturated with gas, and equilibrium is
umes, initial and equilibrium pressures and operating temperature established between the two phases.
as follows: The time needed to reach the equilibrium was less than 2 h in
n (P0 − Pf )VG /(ZRT ) all the completed tests. At equilibrium point, the pressure has not
x= = (3) changed for about 20 min. The number of moles of CO2 absorbed in
n + nL (P0 − Pf )VG /(ZRT ) + (L VL /M)
solvent was calculated using Eq. (1).
2 2 2
ıx ∂x ı(P0 − Pf ) ∂x ıT
= × + ×
x ∂(P0 − Pf ) P0 − Pf ∂T T
2 2
∂x ıVL ∂x ıVG
+ × + × (4)
∂VL VL ∂VG VG
where x is the CO2 solubility in NMP, P0 and Pf are initial and equilib-
rium pressures of CO2 , VG and VL are gas volume and liquid volume
in the absorption cell, Z is CO2 compressibility factor, T is operating
temperature, L and M are density and molecular weight of NMP.
According to the above uncertainty analysis described by Mof-
fat [18], the measurement errors of the equilibrium mole fraction
was ±3.5%. The contributions of the main parameters involving in
the evaluation of equilibrium mole fraction according to Eq. (4)
were summarized in Table 2. The analysis showed that pressure
measurement had the greatest effect on the uncertainty of the equi-
librium mole fraction. Furthermore, to check the reproducibility
of the experiments, some runs were repeated. The repetition of Fig. 3. Pressure decay during CO2 absorption in NMP at 313.15 K.
M.R. Bohloul et al. / Fluid Phase Equilibria 365 (2014) 106–111 109
Table 3
Values of the critical properties, acentric factors and molecular weights.
0.077796.R.Tc
b= (7)
Pc
2
˛(T ) = [1 + k(1 − Tr0.5 )] (8)
2
k = 0.37464 + 1.54226ω − 0.26992ω (9)
Fig. 4. Pressure decay during CO2 absorption in NMP at 333.15 K. 3.2.2. Mixing rules
EOS can be used for multicomponent systems by employing
the appropriate mixing rules to determine their parameters for
The equilibrium data for the CO2 –NMP binary system from this
mixtures. The parameters of EOS are considered to present the
study are shown in Fig. 5. In this figure, equilibrium data of the
attractive and repulsive forces between the molecules [23]. The
present work were compared with the equilibrium data of the same
attractive (amix ) and repulsive (bmix ) parameters of the fluid mixture
system which was published by Rajasingam et al. [2] and Sweeney
were calculated with mixing rules as follows [2]:
[19] in the literature and a good agreement can be observed. It
must be mentioned that Sweeney [19] presented his solubility data
N
N
at two different temperatures using Henry’s law constant. There- amix = xi .xj .(1 − kij )(ai .aj )1/2 (10)
fore, his data were shown as straight lines in Fig. 5. At a similar i j
pressure, the solubility of CO2 decreases with increasing tempera-
ture, and at a constant temperature, the solubility of CO2 increases
N
N
bi + bj
with increasing pressure. The experimental data confirmed that bmix = xi .xj .(1 − lij ) (11)
2
solubility of CO2 has a near-linear dependence on pressure. i j
Table 4
Optimized values of binary interaction parameters (kij and lij ) used in PR-EOS.
Table 5
The values of solubility calculated by 1-parameter (1-p) and 2-parameter (2-p) versions of PR-EOS at various temperatures for CO2 + NMP binary system.
P (MPa) x(exp) x(1-p) x(2-p) (1-p) AARD (%) (2-p) AARD (%)
T = 293.15 K
0.980 0.1494 ± 0.052 0.1445 0.1470 2.6 1.6
0.864 0.1295 ± 0.045 0.1287 0.1293 0.6 0.1
0.766 0.1124 ± 0.039 0.1144 0.1144 1.7 1.7
0.698 0.0981 ± 0.034 0.1044 0.1041 6.4 6.1
0.564 0.0847 ± 0.002 0.0847 0.0839 0.0 0.9
Average AARD 2.2 2.0
T = 313.15 K
1.120 0.1181 ± 0.004 0.1188 0.1181 0.5 0.0
0.955 0.1039 ± 0.003 0.1017 0.1013 2.1 2.5
0.889 0.0896 ± 0.003 0.0949 0.0946 5.9 5.5
0.718 0.0820 ± 0.002 0.0770 0.0770 6.0 6.0
0.658 0.0707 ± 0.002 0.0707 0.0707 0.0 0.0
Average AARD 2.9 2.8
T = 333.15 K
1.120 0.0974 ± 0.003 0.0908 0.0909 6.7 6.6
1.040 0.0846 ± 0.002 0.0846 0.0846 0.0 0.9
0.935 0.0750 ± 0.002 0.0761 0.0759 1.4 1.2
0.839 0.0718 ± 0.002 0.0685 0.0682 4.5 5.0
0.696 0.0556 ± 0.001 0.0570 0.0576 2.5 3.5
Average AARD 3.0 3.2
of the binary interaction parameters (kij and lij ) were obtained by effect of second parameter (lij ) on the accuracy of PR-EOS is neg-
minimizing the average absolute relative deviations (AARD) with ligible. Therefore, it is better to work with 1-parameter model for
respect to mole fraction: the sake of simplicity.
100 xcal − xexp
M
List of symbols
AARD (%) = xexp (12)
M CO2 carbon dioxide
i=1 i
NMP N-methyl-2-pyrrolidone
where M is the number of data points, xcal and xexp are the calcu- R gas constant
lated and the experimental values of the mole fraction, respectively. T temperature
Optimized values of binary interaction parameters (kij and lij ) are Vi initial volume of the gas
presented in Table 4. Vf final volume of the gas
Z compressibility factor
3.4. Solubility prediction using PR-EOS nCO2 the number of absorbed moles of CO2 in NMP
PR-EOS Peng–Robinson equation of state
1-Parameter and 2-parameter versions of PR-EOS along with the E volume expansion
optimized binary interaction parameters were applied to calculate X CO2 solubility
the solubility of CO2 in NMP. The resulted values of solubility are VL liquid volume
presented and compared with the experimental data in Table 5 for VG gas volume
three different temperatures. The values of solubility of CO2 with 2- A attractive parameter
parameter model were slightly closer to the experimental data. This b repulsive parameter
shows that the effect of second parameter (lij ) on the accuracy of PR- molar volume
EOS is negligible and it is better to work with 1-parameter model ω acentric factor
for the sake of simplicity. The difference between the experimental Tc critical temperature
data and the predicted values derived from PR-EOS were acquired Pc critical pressure
using the average absolute relative deviations (AARD). Tr reduced temperature
amix attractive parameter for the mixture
bmix repulsive parameter for the mixture
4. Conclusion kij , lij binary interaction parameters
MW molecular weight
The solubility of CO2 in NMP was examined at different temper- AARD average absolute relative deviation
atures and pressures. The experimental solubility data increases
with decreasing temperature and increasing pressure. In addition,
Acknowledgements
the results showed that NMP has weak ability to absorb CO2 at
pressures up to 1.1 MPa. The values of binary interaction param- This research was supported by the National Iranian South Oil
eters (kij and lij ) implemented in PR-EOS were optimized using Company (NISOC, Ahwaz, Iran) (www.nisoc.ir). Special thanks of
the experimental data. Then the solubility data were calculated authors to Dr. Lucien (The University of New South Wales, Sydney,
using PR-EOS with the optimized binary interaction parameters. Australia), Dr. Shima Azizi (Islamic Azad University of Mahshahr
The mean deviation between the experimental data and the pre- Branch, Mahshahr, Iran) and Transport Phenomena Research Cen-
dicted values by thermodynamic model was determined to be less ter of Islamic Azad University of Mahshahr Branch.
than 3.5%. It was shown that 1-parameter and 2-parameter versions
of PR-EOS have similar accuracy. However, the predicted values of References
solubility of CO2 with 2-parameter model were slightly closer to
the experimental data. Consequently, it can be concluded that the [1] A. Kohl, R. Nielsen, Gas Purification, fifth ed., Gulf Publishing, Houston, TX, 1997.
M.R. Bohloul et al. / Fluid Phase Equilibria 365 (2014) 106–111 111
[2] R. Rajasingam, L. Lioe, Q.T. Pham, F.P. Lucien, Solubility of carbon dioxide in di- [13] S. Portier, C. Rochelle, Modelling CO2 solubility in pure water and NaCl-type
methylsulfoxide and N-methyl-2-pyrrolidone at elevated pressure, J. Supercrit. waters from 0 to 300◦ C and from 1 to 300 bar: Application to the utsira forma-
Fluids 31 (2004) 227–234. tion at sleipner, Chem. Geol. 217 (2005) 187–199.
[3] I. Dalmolin, E. Skovroinski, A. Biasi, M.L. Corazza, C. Dariva, J. Vladimir Oliveira, [14] G.D. Pappaa, C. Perakisa, N. Tsimpanogiannisb, C. Voutsas, Thermodynamic
Solubility of carbon dioxide in binary and ternary mixtures with ethanol and modeling of the vapor–liquid equilibrium of the CO2 /H2 O mixture, Fluid Phase
water, Fluid Phase Equilib. 245 (2006) 193–200. Equilib. 284 (2009) 56–63.
[4] C. Wang, A.N. Soriano, Z. Yang, M. Li, Solubility of carbon dioxide in the sol- [15] R.M. Charin, M.L. Corazza, P.M. Ndiaye, M.A. Mazutti, J.V. Oliveira,
vent system (2-amino-2-methyl-1-propanol + sulfolane + water), Fluid Phase Phase equilibrium data and thermodynamic modeling of the system
Equilib. 291 (2010) 195–200. propane + NMP + methanol at high pressures, J. Supercrit. Fluids 55 (2010)
[5] M. Esther, A. Trejo, Gas solubility of CO2 in aqueous solutions of 23–31.
N-methyldiethanolamine and diethanolamine with 2-amino-2-methyl-1- [16] Y. Yan, C. Chen, Thermodynamic modeling of CO2 solubility in aqueous solu-
propanol, Fluid Phase Equilib. 218 (2004) 261–267. tions of NaCl and Na2 SO4 , J. Supercrit. Fluids 55 (2010) 623–634.
[6] J. Chen, T. Liu, W. Yuan, L. Zhao, Solubility and diffusivity of CO2 in [17] M.D. Bermejo, A. Martín, L.J. Florusse, C.J. Peters, M.J. Cocero, The influence of
polypropylene/micro-calcium carbonate composites, J. Supercrit. Fluids 77 Na2 SO4 on the CO2 solubility in water at high pressure, Fluid Phase Equilib. 238
(2013) 33–43. (2005) 220–228.
[7] M. Taib, T. Murugesan, Solubilities of CO2 in aqueous solutions of ionic liquids [18] R.J. Moffat, Describing the uncertainties in experimental results, Exp. Therm.
(ILs) and monoethanolamine (MEA) at pressures from 100 to 1600 kPa, Chem. Fluid Sci. 1 (1988) 3–17.
Eng. J. 181–182 (181) (2012) 56–60. [19] C.W. Sweeney, Solubilities and partial molar anthalpies of solution for polar
[8] E. Aronu, S. Gondal, E. Hessen, T. Warberg, A. Hartono, A. Hoff, F. Svendsen, gas–liquid systems determined by gas chromatography, Chromatographia 18
Solubility of CO2 in 15, 30, 45 and 60 mass% MEA from 40 to 120 ◦ C and model (1984) 663–667.
representation using the extended UNIQUAC framework, Chem. Eng. Sci. 66 [20] D. Peng, D.B. Robinson, A new two-constant equation of state, Ind. Eng. Chem.
(2011) 6393–6406. 15 (1976) 59–64.
[9] M. Vahidi, N. Matin, M. Goharrokhi, M. Hosseini Jenab, M. Abedinzadegan Abdi, [21] A.M. Abudour, S.A. Mohammad, A.M. Gasem, Modeling high-pressure phase
H. Najibi, Correlation of CO2 solubility in N-methyldiethanolamine + piperazine equilibria of coalbed gases/water mixtures with the Peng–Robinson equation
aqueous solutions using extended Debye–Hückel model, J. Chem. Thermodyn. of state, Fluid Phase Equilib. 319 (2012) 77–89.
41 (2009) 1272–1278. [22] K.S. Abdulkadirova, J. Peters, J.V. Sengersa, M.A. Anisimov, An isomor-
[10] H. Pahlavanzadeh, S. Nourani, M. Saber, Experimental analysis and modelling of phic Peng–Robinson equation for phase-equilibria properties of hydro-
CO2 solubility in AMP (2-amino-2-methyl-1-propanol) at low CO2 partial pres- carbon mixtures in the critical region, J. Supercrit. Fluids 55 (2010)
sure using the models of Deshmukh–Mather and the artificial neural network, 594–602.
J. Chem. Thermodyn. 43 (2011) 1775–1783. [23] A. Danesh, PVT and Phase Behavior of Petroleum Reservoir Fluids, first ed.,
[11] X. Yuan, S. Zhang, J. Liu, X. Lu, Solubilities of CO2 in hydroxyl ammonium ionic Elsevier, Amsterdam, Netherlands, 1998.
liquids at elevated pressures, Fluid Phase Equilib. 257 (2007) 195–200. [24] M.A. McHugh, V.J. Krukonis, Supercritical Fluid Extraction Principles and
[12] H.K. Balsora, M.K. Mondal, Solubility of CO2 in aqueous TSP, Fluid Phase Equilib. Practice Butterworth–Heinemann Series in Chemical Engineering, 2nd ed.,
328 (2012) 21–24. 1994.