Experiments, Modeling, and Simulation of CO2 Dehydration by Ionic

Aug 5, 2019 - Experiments, Modeling, and Simulation of CO2 Dehydration by Ionic ..... of the hydrogen bonding interaction (HHB) to excess enthalpy is ...
0 downloads 0 Views 512KB Size
Subscriber access provided by RUTGERS UNIVERSITY

Separations

Experiments, Modeling, and Simulation of CO2 Dehydration by Ionic Liquid, Triethylene Glycol, and Their Binary Mixtures Yifan Jiang, Mohsen Taheri, Gangqiang Yu, Jiqin Zhu, and Zhigang Lei Ind. Eng. Chem. Res., Just Accepted Manuscript • DOI: 10.1021/acs.iecr.9b02540 • Publication Date (Web): 05 Aug 2019 Downloaded from pubs.acs.org on August 6, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Experiments, Modeling, and Simulation of CO2 Dehydration by Ionic Liquid, Triethylene Glycol, and Their Binary Mixtures

Yifan Jiang, Mohsen Taheri, Gangqiang Yu, Jiqin Zhu,* and Zhigang Lei* State Key Laboratory of Chemical Resource Engineering, Beijing University of Chemical Technology, Box 266, Beijing 100029, China

ABSTRACT: Ionic liquids (ILs) as a new class of gas dehydration solvents, in combination with the traditional triethylene glycol (TEG) solvent were first proposed for carbon dioxide (CO2) dehydration. Among 210 preliminary ILs, the hydrophilic [BMIM]+[BF4]- was selected based on the COSMO-RS model-involved IL screening methodology. Solubility of CO2 in pure TEG, pure [BMIM]+[BF4]-, binary mixture of TEG + [BMIM]+[BF4]-, and ternary mixture of TEG + [BMIM]+[BF4]- + H2O were measured experimentally. Two new binary interaction parameters (BIPs) were introduced by correlating a series of experimental data with the UNIFAC-Lei model. The COSMO-RS model along with the reduced density gradient (RDG) method was applied to interpret the nature of interaction between molecules. The CO2 dehydration experiment was conducted in a laboratory-scale absorption tower. The process simulation indicates that, in comparison with pure TEG, the use of IL purely or mixed with TEG improves both separation performance and process energy penalty.

1

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1. INTRODUCTION With the rapid development of industrialization, emissions of greenhouse gases have caused severe environmental problems. Among the primary greenhouse gases, carbon dioxide (CO2) has the greatest impact on the global warming owning to the large volume of its constant anthropogenic emissions.1 The first solution for global warming is to control and reduce CO2 emissions. However, the current scientific and technological endeavours are not sufficient to change the inherent energy demand structure of the globe. Thus, separation and recovery of CO2 from industrial processes become the key point to the alleviation of destructive effect of this long-lived greenhouse gas. At present, the recovered CO2 is mainly utilized in two types of processes: one is to use it as an industrial raw material for processing and synthesis;2 the other is to seal it in the underground geological formation.3 In both applications, CO2 needs to be transported through pipelines. Pipeline transportation has certain requirement for moisture content. To avoid the corrosion and blockage of pipelines, dehydration must be carried out prior to the gas transportation. The triethylene glycol (TEG) dehydration method is widely applied in industry for CO2 dehydration, since TEG has the advantages of high moisture absorption and easy regeneration.4 However, there are some problems associated with the TEG dehydration process: (1) the reboiler temperature needs to be strictly controlled, because an excessive temperature will lead to the solvent decomposition;5 and (2) CO2 promotes the degradation process of TEG, thus reducing the solvent replacement time and increasing the operating cost.6 Over the last several years, the applications of ionic liquids (ILs) in separation and purification processes,7-10 organic synthesis and catalysis,11 electrochemistry,12 and material

2

ACS Paragon Plus Environment

Page 2 of 36

Page 3 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

chemistry13 have been studied by many researchers. In separation processes, the high stability of ILs makes them easy to regenerate, the low vapor pressure of ILs reduces atmospheric pollution, and the designability of ILs brings about the possibility of selection of a suitable IL for a given process.14 Therefore, in this work the usage of IL-TEG hybrid solvents in CO2 dehydration processes was investigated. IL and TEG mixed as the hybrid solvents can benefit from the advantages of each other, which can subsequently increase the separation efficiency of CO2 dehydration process. The aim of this work is divided into the following sections: (1) the most suitable IL for CO2 dehydration experiment was selected from 210 preliminary IL candidates, using a COSMO-RS-based IL solvent screening method and taking into account the Henry’s constants (Hi) of CO2 and H2O, the selectivity of CO2 to H2O (S = HCO2/HH2O), and other physical properties such as IL stability; (2) experimental measurements were performed to obtain the solubility of CO2 in pure TEG, pure IL, binary mixture of TEG + IL, and ternary mixture of TEG + IL + H2O, and the predicted values by UNIFAC-Lei model were compared with experimental data; (3) the CO2 dehydration experiment using IL-TEG as absorbent in an absorption tower was carried out at the laboratory-scale; (4) an in-depth analysis on CO2 dehydration mechanism by the reduced density gradient (RDG) method and COSMO-RS model was made; and (5) the CO2 gas dehydration process simulation was conducted within the ASPEN PLUS framework. 2. EXPERIMENTAL SECTION 2.1. Materials The chemical materials used in this study including CO2, TEG, and IL were purchased

3

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

from chemical markets. Detailed information on these chemicals can be found in Supporting Information (SI) Table S1. 2.2. Apparatus and Procedure 2.2.1. CO2 Solubility Measurement The solubility data of CO2 in pure IL, pure TEG, binary mixture of IL + TEG, and ternary mixture of IL + TEG + H2O were measured over a wide temperature (from 273.15 to 353.15 K) and pressure (from 0.5 to 3.5 MPa) range using a gas-liquid equilibrium (GLE) apparatus. Ethanol and silicone oil as the medium for cooling and heating, respectively, were used in the experimental apparatus to achieve the required low and high temperatures. Details on the structure and operating procedure of experimental apparatus have been given elsewhere.14 Fluctuations of the pressure gauge and temperature sensor used in the experiments were ± 0.001 MPa and ± 0.1 K, respectively. 2.2.2. Gas drying experiment The CO2 drying experiment was carried out in an absorption tower (Φ30 mm × 1 m). A scheme of gas drying experimental setup is presented in Figure S1. The gas from a CO2 cylinder passes through a regulator, a gas flow regulator, and a buffer tank containing H2O. The saturated H2O gas enters from the bottom of absorption tower. The solvent passes through the advection pump from the top of absorption tower so that the two streams are in countercurrent contact. The H2O content in the feeding CO2 gas was measured on-line by a dew point meter (the type RHD-601), around 20000 ppm (molar fraction basis) at 298.2 K and atmospheric pressure. The flow rate of the feeding CO2 gas was kept at 800 mL·min-1.

4

ACS Paragon Plus Environment

Page 4 of 36

Page 5 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

3. MODELING SECTION 3.1. UNIFAC-Lei Model The UNIFAC-Lei model, an efficient activity coefficient-based group contribution predictive model, has been extensively applied by many researchers.7,15-17 It extends the original UNIFAC model to IL systems, and is written as ln  i  ln  iC  ln  iR

(1)

where ln  Ci is directly related to the shape and size of groups, which contains the group area and volume parameters Rk and Qk (see SI Table S2 ); and ln  iR is only related to the areas and interactions of groups, which contains a pair of binary interaction parameters (BIPs)

amn and anm . Information on BIPs correlation is available in our previous work.18,19 In this work, six UNIFAC-Lei structural groups (CH2, CO2, OH, CH2O, H2O, and

 MIMBF4  )

are concerned, which are listed in Table 1.18,20,21 The BIPs ( amn and anm )

between CO2 and CH2O, as well as those between CH2O and

 MIMBF4  ,

are unknown.

Thus, they were acquired using the CO2 solubility data measured in this work as well as those collected from the literature.22 The objective function (OF) was used to acquire the unknown

amn and anm :  1 OF  min   N

N

 1

xcal  xexp   xexp 

(2)

where x is the solubility of CO2 (in mole fraction) in liquid phase; and N is the number of data points. For the BIPs acquired in this work, a three-fold cross validation method was used to verify the predictive capability of UNIFAC-Lei model. The solubility data of CO2 in TEG and in the mixture of TEG and IL are randomly divided into three groups of the same type. The BIPs were correlated with the two sets of data, and verified with the left one. Therefore, 5

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 36

three models (Model 1/3, Model 2/3, and Model 3/3) were set up in this procedure, and each model had 1/3 of the data. The average relative deviation (ARD) was applied to select the best values of BIPs. Table 2 lists the BIPs between CO2 and CH2O and between [MIM][BF4] and CH2O ( amn and anm ), along with the ARD values. It can be seen that ARD is less than 0.08 in most of the cases, indicating the validity of UNIFAC-Lei model. Gas-liquid equilibrium (GLE) for the CO2 + solvent systems can be written as

i yi P   i xi Pi s

(3)

where xi and yi are the mole fractions of CO2 in the liquid and gas phases, respectively; P and Pi S

are the system pressure and the saturated vapor pressure of CO2 at a certain

temperature;23 i is the gas phase fugacity coefficient of CO2;24 and γi is the CO2 activity coefficient calculated by the UNIFAC-Lei model. 3.2. COSMO-RS Model Screening a suitable IL by means of predictive thermodynamic models can save time and cost of the experiments. Herein, among 210 preliminary ILs, the statistical thermodynamic-based COSMO-RS model was utilized to explore a matching IL for gas dehydration experiments. The two-dimensional structures, names, and abbreviations of ILs used in this study are given in SI Table S3. Figure S2a shows the influence of cation and anion types on the COSMO-RS predicted HCO2 in ILs. Moreover, for most of the anions (e.g., [ClO4]-, [AC]-, and [TF2N]-), the longer the cationic alkyl chain length in ILs, the higher the solubility of CO2 in ILs. Therefore, [BMIM]+ with a shorter alkyl chain is selected as the cationic moiety. As shown in Figures S2b and S2c, anions play the dominant roles in determining the HH2O and selectivity of CO2 to H2O (S). Although S is the largest in the ILs

6

ACS Paragon Plus Environment

Page 7 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

with the anion [AC]-, the [AC]- is not a thermally stable anion.25 Moreover, it exhibits chemical absorption for CO2, which is clearly not conducive to the dehydration process.26 In addition, H2O has a higher solubility in ILs in which the anion is halogen; however, their melting points are high, that is, they are solid at room temperature.27,28 Therefore, the hydrophilic-IL [BMIM]+[BF4]- was selected as a suitable CO2 dehydration IL. In addition to the low melting temperature and good chemical stability, [BMIM]+[BF4]- is completely miscible with H2O at room temperature.25,29,30 Although the hydrolysis of BF4-based ILs may occur under certain conditions, the H2O content in feeding gas is relatively small and thus hydrolysis effect is not considered in this work. The COSMO-RS model is a powerful and efficient tool which allows to estimate the excess enthalpies for the binary or ternary mixture systems including ILs.31-33 In the calculation, the excess enthalpy is the sum of three specific interactions, i.e., electrostatic-misfit interaction (HMF), hydrogen bonding interaction (HHB), and van der Waals interaction (HvdW):

H m  H MF  H HB  H vdw

(4)

3.3. RDG Analysis RDG analysis is used for the interpretation of intermolecular and intramolecular weak interactions.34 Weak interactions mainly include hydrogen bonding (HB), van der Waals (vdW), and steric interactions. Herein, Multiwfn35 software was used for RDG calculation, and VMD36software was used for the display of weak interactions. Compared to the COSMO-RS model, RDG analysis further sheds light on the interaction mechanism between CO2/H2O molecule and the solvent molecule.

7

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

4. PROCESS SIMULATION The CO2 dehydration process simulation was conducted within the ASPEN PLUS framework, the rigorous equilibrium stage (EQ-stage) model being adopted. The employed IL ([BMIM]+[BF4]-) was added manually into the ASPEN program, while TEG, H2O, and CO2 were added from the ASPEN PLUS (V7.2) database. The newly obtained parameters along with the previously available UNIFAC-Lei BIPs were imported into the UNIFAC property model so that the processes can be simulated accurately based on the experimentally-derived parameters. 5. RESULTS AND DISCUSSION 5.1. Solubility of CO2 in Solvents Figures 1 and 2 show the solubility of CO2 (1) in pure TEG (2), [BMIM]+[BF4]- (3), the binary mixture of TEG (2) + [BMIM]+[BF4]- (3) at different contents (x2 = 0.2, 0.5, 0.8), and the ternary mixture of TEG (2) + [BMIM]+[BF4]- (3) + H2O (4) at different contents (x2 = 0.14, 0.35, 0.56), along with the UNIFAC-Lei predicted values. It can be observed that the results follows the general trend. Under the same temperature and pressure, the solubility of CO2 in mixed solvents increases with the increase of [BMIM]+[BF4]- content in the mixture. Furthermore, the UNIFAC-Lei model shows the reliable predictive power for the solubility of CO2 in either single (TEG or [BMIM]+[BF4]-) or mixed solvents. Comparison of the UNIFAC-Lei predicted values versus experimental data is given in SI Table S4. As a whole, the ARDs are less than 10%, thus, UNIFAC-Lei model is suitable for predicting the CO2 solubility in the hybrid TEG-IL solvents. As shown in Figure 2 and Figure S3, the presence of H2O in solvents reduces the CO2 solubility, which facilitates the acquisition of more CO2

8

ACS Paragon Plus Environment

Page 8 of 36

Page 9 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

product in the gas dehydration processes. The detailed CO2 solubility data are listed in SI Tables S5 and S6. The Henry’s constant is also an important thermodynamic property, which is calculated by37

H 1 (T,P)  lim

x1  0

f1L P (T , P, y1 )  lim 1 x1 x1 0 x1

(5)

where H 1 (T,P) is the Henry’s constant of CO2 (1) in TEG, [BMIM]+[BF4]-, or TEG + [BMIM]+[BF4]- (expressed in mole fraction); f1L is the fugacity of CO2 in liquid phase; and

1 (T , P, y1 ) is the fugacity coefficient of pure CO2 in gas phase.24 The Henry’s constants of CO2 in pure and mixed solvents by linear regression at temperatures (273.15, 293.15, 313.15, 333.15, and 353.15 K) are given in Table 3, together with the calculated values by UNIFAC-Lei model. Obviously, [BMIM]+[BF4]- has the lower Henry's constants than other solvents, that is, CO2 is more soluble in [BMIM]+[BF4]- than in other solvents investigated in this work. Moreover, the Henry’s constants predicted by UNIFAC-Lei model were compared with those obtained experimentally. It is clear from Table 3 that both agree well. 5.2. CO2 Dehydration Process Experiment Using TEG-IL as Absorbent The solvent used in the CO2 dehydration process experiment at temperature of 298.2 K and atmospheric pressure was TEG (x2 = 0.5) + [BMIM]+[BF4]- (x3 = 0.5). The main reason for choosing this proportion of mixed solvents is based on the dehydration effect of mixed hydrophilic solvents at the same molar ratio. Evidently, the solvent flow rate has a significant effect on separation performance. Figure 3a shows that the H2O content in outlet gas is related to solvent flow rate. As the solvent flow rate increases, the H2O content in outlet gas first drops rapidly, and then becomes stable. Meanwhile, the H2O content in solvent itself 9

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

also greatly affects the H2O content in outlet gas. As shown in Figure 3b, at a given flow rate, the H2O content in outlet gas has a linear relationship with the H2O content in solvent. Thus, in practice these two factors should be taken into account together. The detailed gas dehydration experimental data are given in SI Table S7. 5.3. Excess Enthalpy for CO2-Liquid Systems The excess enthalpies of one ternary mixture (CO2 + TEG (x2=0.5) + [BMIM]+[BF4](x3=0.5)) and two binary mixtures (CO2 + TEG and CO2 + [BMIM]+[BF4]-) at 298.15 K are presented in Figure 4. At low CO2 concentration region, the negative excess enthalpy of CO2-TEG system indicates that the van der Waals interaction (HvdW) is dominant, bringing out the exothermicity as depicted in Figure 4a. However, when the mole fraction of CO2 exceeds 0.8, the contribution of hydrogen bonding interaction (HHB) to excess enthalpy is more than HvdW, while the electrostatic-misfit interaction (HMF) can almost be ignored. For the CO2-[BMIM]+[BF4]- system, Figure 4b presents the predicted excess enthalpies. Unlike the CO2-TEG system, the excess enthalpy of CO2 + [BMIM]+[BF4]- is negative in the whole concentration range, indicating the exothermic process. Moreover, the absolute value of excess enthalpy of CO2 + [BMIM]+[BF4]- is much greater than that of CO2 + TGE, indicating a stronger interaction between CO2 and [BMIM]+[BF4]-, which makes the solubility of CO2 in ILs much greater than in TEG. To understand the solubility mechanism of CO2 in binary mixtures, the excess enthalpies of CO2 + TEG (x2=0.5) + [BMIM]+[BF4]- (x3=0.5) are shown in Figure 4c. The Hm, HHB, HvdW, and HMF in ternary system seem to lie between those in binary systems. More specially, the solubility of CO2 in binary mixtures is also between those in pure TEG and IL, which means that excess enthalpy also reflects the change in solubility

10

ACS Paragon Plus Environment

Page 10 of 36

Page 11 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

in some degree. Similarly, the trend of excess enthalpies of CO2 + TEG (x2=0.2) + [BMIM]+[BF4]- (x3=0.8) and CO2 + TEG (x2=0.8) + [BMIM]+[BF4]- (x3=0.2) is also related to the solubility. Their excess enthalpies are shown in SI Figure S4. The excess enthalpies (Hm, HHB, HvdW, and HMF) of binary and ternary mixtures at 298.15 K are depicted in Figure 5. It was found that both Hm and HvdW become more negative with the increase of [BMIM]+[BF4]- content in solution, exhibiting the strong exothermicity. However, the strong exothermicity of the mixture is favorable for improving the solubility of CO2 in IL.38 Moreover, HvdW makes the largest contribution to excess enthalpy, indicating that the difference in CO2 interactions with TEG and [BMIM]+[BF4]- mainly results from van der Waals force. As for HHB and HMF, with the increase of [BMIM]+[BF4]- content, HHB becomes smaller, whereas HMF becomes higher, but still lower than HHB. Thus, HvdW and HHB have a more significant contribution to excess enthalpy. Figure 6a shows the Hm of equimolar CO2-solvent mixtures calculated by COSMO-RS model versus the Henry’s constants calculated by UNIFAC-Lei model at 298.15 K. As the solubility of CO2 in solvent increases (corresponding to the small Henry’s constants), exothermic effect of the mixture becomes stronger (corresponding to the reduced Hm ). Figure 6b further analyzes the effect of HHB, HvdW, and HMF on the solubility of CO2 in solvent. As the solubility of CO2 in solvent increases, HvdW and HMF also increase, while HHB shows an opposite trend. The detailed Henry’s constants and excess enthalpies are given in SI Table S8. Figure 7 shows the comparison between the excess enthalpy of H2O + TEG + [BMIM]+[BF4]- and that of CO2 + TEG + [BMIM]+[BF4]-. Under the same molar

11

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

composition, the excess enthalpies of CO2 and solvent are negative, in contrast to the positive excess enthalpies of H2O and solvent (see Figure 7a). In addition, Figure 7b further reveals that under the same molar composition, HHB plays the key role in the excess enthalpy of H2O and solvent, while HvdW does in the excess enthalpy of CO2 and solvent. This indicates that the main interaction between H2O and solvent is hydrogen bond formed to remove the moisture from gas mixture. 5.4. RDG Analysis for the CO2/H2O + TEG + IL Systems Figure 8a shows a large green iso-surface between TEG and [BMIM]+[BF4]-, indicating that the two compounds are bound together by vdW interaction. Furthermore, inside the [BMIM]+[BF4]-, the hydrogen atom on the imidazole ring forms a strong HB with the fluorine atom. The red region inside the imidazole ring indicates a strong steric effect. For the TEG + [BMIM]+[BF4]- + CO2 system (see Figure 8b), the green region between CO2 and TEG indicates vdW interaction, while the yellow-green region represents a weak spatial steric effect. For the TEG + [BMIM]+[BF4]- + H2O system (see Figure 8c), the dark blue region indicates that the hydrogen atoms in H2O molecule form a strong HB with the fluorine atoms in [BF4]-. In addition, there is a green region between the oxygen atoms in H2O molecule and the hydrogen atoms in the methylene group of TEG, indicating the formation of a weak HB. It can be inferred that the formation of HB plays the key role in the removal of H2O molecules from gas mixture, which is consistent with the excess enthalpy analysis aforementioned. The molecular structures optimized by Gaussian 09 software are shown in SI Figure S5. 5.5. Process Simulation on CO2 Gas Dehydration with TEG, IL, and the Mixed TEG +

12

ACS Paragon Plus Environment

Page 12 of 36

Page 13 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

IL as Absorbents Figure 3 shows that the simulated values by the EQ-stage model with the UNIFAC property model agree well with the experimental data, thus verifying the reliability of UNIFAC-Lei model for CO2 dehydration using the mixed TEG + IL solvent as absorbent. Figure 9 shows the CO2 dehydration process flowsheet using three different absorbents (pure IL, pure TEG, and 50 wt% TEG + 50 wt% IL). As shown in Figure 9a, the pure IL contacts countercurrent with the moisture-containing CO2 gas in absorption column, in which the product gas is discharged at the top, while the IL rich in moisture passes through a heat exchanger and then enters the flash tank. The recycled absorbent passes through the condenser and enters the absorber. Since TEG is more volatile than IL, for the dehydration process using TEG or TEG + IL as absorbent, a distillation column is used instead of flash tank for solvent recovery (see Figure 9b). The following specifications are considered in process simulation: (1) maintaining the same specification of gas feed and solvent, and (2) the molar fraction of H2O in product gas doesn’t exceed 600 ppm. Optimized specifications of the three processes are given in SI Tables S9. As given in Table 4, the pure IL process achieves the best separation performance and the lowest energy consumption under the same gas feed (1000 kg/h, 2% mole fraction in water content) and absorbent doses (3500 kg/h). However, for pure TEG process, both separation performance and energy consumption are the worst. In addition, compared to the pure TEG process, both separation performance and energy consumption of the TEG + IL process are improved due to the addition of IL. Therefore, the application of TEG + IL process can not only reduce the cost of pure IL in some degree, but also achieve a better separation performance.

13

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

6. CONCLUSIONS In this work, the mixed IL and TEG was first proposed as absorbent for gas dehydration. The hydrophilic IL ([BMIM]+[BF4]-) was selected after considering such factors as solubility and chemical stability. Moreover, the COSMO-RS model and reduced density gradient (RDG) method were used together to reveal the different types of interactions. It was proven that the strong hydrogen bond (HB) formed by H2O molecule and the anion ([BF4]-) plays a key role in gas dehydration. The solubility data of CO2 in pure TEG, pure [BMIM]+[BF4]-, binary mixture of TEG + [BMIM]+[BF4]-, and ternary mixture of TEG + [BMIM]+[BF4]- + H2O were measured. The UNIFAC-Lei predicted values show a good consistency with experimental data. The binary mixture of TEG and [BMIM]+[BF4]- as absorbent exhibits an excellent CO2 dehydration performance. Furthermore, the process simulation on CO2 dehydration was conducted within the ASPEN PLUS framework. In comparison with the benchmark pure TEG solvent, pure IL or a combination of IL-TEG solvent can improve both separation performance and energy consumption. In the end, it is noting that no extra streams and equipments are needed to add based on the original TEG process, when using the TEG + IL solvent as absorbent.

14

ACS Paragon Plus Environment

Page 14 of 36

Page 15 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research



ASSOCIATED CONTENT

Supporting Information The Supporting Information is available free of charge on the ACS Publications website. Detailed experimental solubility data and predicted values by the UNIFAC-Lei model, the flowsheet of gas dehydration experiment, and the optimized process parameters can be found in the online version (xls).  AUTHOR INFORMATION Corresponding Author *Phone:

+86-10-64433695; E-mail: zhujq@mail.buct.edu.cn (J. Zhu).

E-mail: leizhg@mail.buct.edu.cn (Z. Lei). ORCID Jiqin Zhu: 0000-0003-1187-0522 Zhigang Lei: 0000-0001-7838-7207 Yifan Jiang: 0000-0003-4855-9340 Gangqiang Yu: 0000-0002-3595-6972 Notes The authors declare no competing financial interest.  ACKNOWLEDGMENTS This work is financially supported by the National Key R&D Plan of China (No. 2018YFB0604702) and the National Natural Science Foundation of China under Grant (No. U1862103).

15

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

 REFERENCES (1) Work, S. N. Kinetics of core material dissolution in the presence of inhibitors for application in geologic carbon sequestration. Rice University, 2010. (2) Ramachandriya, K. D.; Kundiyana, D. K.; Wilkins, M. R.; Terrill, J. B.; Atiyeh, H. K.; Huhnke, R. L. Carbon dioxide conversion to fuels and chemicals using a hybrid green process. Appl. Energy 2013, 112, 289-299. (3) Murai, S.; Fujioka, Y. Challenges to the Carbon Dioxide Capture and Storage (CCS) Technology. IEEJ Trans. Electr. Electron. Eng. 2010, 3, 37-42. (4) Kemper, J.; Sutherland, L.; Watt, J.; Santos, S. Evaluation and Analysis of the Performance of Dehydration Units for CO2 Capture. Energy Procedia 2014, 63, 7568-7584. (5) Netusil, M.; Ditl, P. Comparison of three methods for natural gas dehydration. J. Nat. Gas Chem. 2011, 20, 471-476. (6) Bahadori, A.; Vuthaluru, H. B.; Mokhatab, S. Analyzing solubility of acid gas and light alkanes in triethylene glycol. J. Nat. Gas Chem. 2008, 17, 51-58. (7) Taheri, M.; Dai, C.; Lei, Z. CO2 capture by methanol, ionic liquid, and their binary mixtures: Experiments, modeling, and process simulation. AIChE J. 2018, 64, 2168-2180. (8) Zhu, Z.; Hu, J.; Geng, X.; Qin, B.; Ma, K.; Wang, Y.; Gao, J. Process design of carbon dioxide and ethane separation using ionic liquid by extractive distillation. J. Chem. Technol. Biotechnol. 2018, 93, 887-896. (9) Hu, Y.; Su, Y.; Jin, S.; Chien, I. L.; Shen, W. Systematic approach for screening organic and ionic liquid solvents in homogeneous extractive distillation exemplified by the tert-butanol dehydration. Sep. Purif. Technol. 2019, 211, 723-737. (10) Zhu, Z.; Ri, Y.; Li, M.; Jia, H.; Wang, Y.; Wang, Y. Extractive distillation for ethanol dehydration using imidazolium-based ionic liquids as solvents. Chem. Eng. Process. 2016, 109, 190-198. (11) Weishi, M.; Tak Hang, C. Ionic-liquid-supported synthesis: a novel liquid-phase strategy for organic synthesis. Acc. Chem. Res. 2006, 39, 897-908. (12) Shen, Y.; Zhang, Y.; Qiu, X.; Guo, H.; Li, N.; Ivaska, A. Polyelectrolyte-functionalized ionic liquid for electrochemistry in supporting electrolyte-free aqueous solutions and application in amperometric flow injection analysis. Green Chem. 2007, 9, 746-753.

16

ACS Paragon Plus Environment

Page 16 of 36

Page 17 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

(13) Brandt, A.; Hallett, J. P.; Leak, D. J.; Murphy, R. J.; Welton, T. The effect of the ionic liquid anion in the pretreatment of pine wood chips. Green Chem. 2010, 12, 672-679. (14) Lei, Z.; Chen, B.; Koo, Y. M.; MacFarlane, D. R. Introduction: Ionic Liquids. Chem. Rev. 2017, 117, 6633-6635. (15) Kamgar, A.; Esmaeilzadeh, F. Prediction of H2S solubility in [hmim][Pf6], [hmim][Bf4] and [hmim][Tf N] using UNIQUAC, NRTL and COSMO-RS. J. Mol. Liq. 2016, 220, 631-634. (16) Akbari, A.; Rahimpour, M. R. Prediction of the solubility of carbon dioxide in imidazolium based ionic liquids using the modified scaled particle theory. J. Mol. Liq. 2018, 255, 135-147. (17) Yu, G.; Dai, C.; Lei, Z. Modified UNIFAC-Lei Model for Ionic Liquid–CH4 Systems. Ind. Eng. Chem. Res. 2018, 57, 7064-7076. (18) Lei, Z.; Zhang, J.; Li, Q.; Chen, B. UNIFAC model for ionic liquids. Ind. Eng. Chem. Res. 2014, 48, 2697-2704. (19) Lei, Z.; Dai, C.; Liu, X.; Xiao, L.; Chen, B. Extension of the UNIFAC Model for Ionic Liquids. Ind. Eng. Chem. Res. 2012, 51, 12135-12144. (20) Skjoldjorgensen, S.; Kolbe, B.; Gmehling, J.; Rasmussen, P. Vapor-Liquid Equilibria by UNIFAC Group Contribution. Revision and Extension. Ind. Eng. Chem. Process Des. Dev. 1979, 18, 2352-2355. (21) Dai, C.; Wei, W.; Lei, Z.; Li, C.; Chen, B. Absorption of CO2 with methanol and ionic liquid mixture at low temperatures. Fluid Phase Equilib. 2015, 391, 9-17. (22) Wise, M.; Chapoy, A. Carbon dioxide solubility in triethylene glycol and aqueous solutions. Fluid Phase Equilib. 2016, 419, 39-49. (23) Shiflett, M. B.; Yokozeki, A. Solubility and diffusivity of hydrofluorocarbons in room‐temperature ionic liquids. AIChE J. 2010, 52, 1205-1219. (24) Span, R.; Wagner, W. A New Equation of State for Carbon Dioxide Covering the Fluid Region from the Triple㏄oint Temperature to 1100 K at Pressures up to 800 MPa. J. Phys. Chem. Ref. Data 1996, 25, 1509-1596. (25) Cao, Y.; Mu, T. Comprehensive Investigation on the Thermal Stability of 66 Ionic Liquids by Thermogravimetric Analysis. Ind. Eng. Chem. Res. 2014, 53, 8651–8664. (26) Gomez-Coma, L.; Garea, A.; Irabien, A. Carbon dioxide capture by [emim][Ac] ionic liquid in a polysulfone hollow fiber membrane contactor. Int. J. Greenhouse Gas Control 2016, 52, 401-409. (27) Ngo, H. L.; Lecompte, K.; Hargens, L.; Mcewen, A. B. Thermal properties of imidazolium ionic 17

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

liquids. Thermochim. Acta 2000, 357, 97-102. (28) Brünig, T.; Krekić, K.; Bruhn, C.; Pietschnig, R. Calorimetric Studies and Structural Aspects of Ionic Liquids in Designing Sorption Materials for Thermal Energy Storage. Chem-Eur J. 2016, 22, 16200-16212. (29) Suarez, P. A. Z.; Einloft, S.; Dullius, J. E. L.; Souza, R. F. D.; Dupont, J. Synthesis and physical-chemical properties of ionic liquids based on 1-n-butyl-3-methylimidazolium cation. J. Chim. Phys. Phys.-Chim. Biol. 1998, 95, 1626-1639. (30) Maiti, A.; Kumar, A.; Rogers, R. D. Water-clustering in hygroscopic ionic liquids-an implicit solvent analysis. Phys. Chem. Chem. Phys. 2012, 14, 5139-5146. (31) Kurnia, K. A.; Coutinho, J. A. P. Overview of the Excess Enthalpies of the Binary Mixtures Composed of Molecular Solvents and Ionic Liquids and Their Modeling Using COSMO-RS. Ind. Eng. Chem. Res. 2013, 52, 13862-13874. (32) Klamt, A. The COSMO and COSMO-RS solvation models. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2011, 1, 699-709. (33) Mokrushina, L.; Buggert, M.; Smirnova, I.; Arlt, W.; Schomäcker, R. COSMO-RS and UNIFAC in Prediction of Micelle/Water Partition Coefficients. Ind. Eng. Chem. Res. 2007, 46, 6501-6509. (34) Johnson, E. R.; Keinan, S.; Morisánchez, P.; Contrerasgarcía, J.; Cohen, A. J.; Yang, W. Revealing noncovalent interactions. J. Am. Chem. Soc. 2010, 132, 6498. (35) Lu, T.; Chen, F. Multiwfn: a multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580-592. (36) Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graphics 1996, 14, 33-38. (37) Lee, B.-C.; Outcalt, S. L., Solubilities of gases in the ionic liquid 1-n-butyl-3-methylimidazolium bis (trifluoromethylsulfonyl) imide. J. Chem. Eng. Data 2006, 51, 892-897. (38) Palomar, J.; Gonzalezmiquel, M.; Polo, A.; Rodriguez, F. Understanding the Physical Absorption of CO2 in Ionic Liquids Using the COSMO-RS Method. Ind. Eng. Chem. Res. 2011, 50, 3452-3463.

18

ACS Paragon Plus Environment

Page 18 of 36

Page 19 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Table Captions Table 1. Group Binary Interaction Parameters (amn and anm) in the UNIFAC-Lei Model

Table 2. Three-fold Cross Validation Results

Table 3. Henry’s Constants of CO2 in Solvents Obtained Experimentally (Hexp) and Predicted by UNIFAC-Lei Model (Hpred) at Different Temperatures

Table 4. Comparison of the Simulation Results among the Three CO2 Dehydration Processes

19

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 36

Table 1. Group Binary Interaction Parameters (amn and anm) in the UNIFAC-Lei Model

a

m

n

amn

anm

CH2 CH2 CH2 CH2 CH2 CO2 CO2 CO2 CO2 [MIM][BF4] [MIM][BF4] [MIM][BF4] OH OH CH2O

CO2 [MIM][BF4] OH CH2O H2O [MIM][BF4] OH CH2O H2O OH CH2O H2O CH2O H2O H2O

107.7 1108.51 986.5 251.5 1318 -14.4413 794.9 137.1718a 497 131.24 -42.2773a -408 28.06 353.5 -314.7

6339 588.74 156.4 83.36 300 430.7991 65.65 748.9029a 386.91 -13.77 -152.6024a 242.88 237.7 -229.1 540.5

Group binary interaction parameters obtained in this work; others from refs.18,20,21

20

ACS Paragon Plus Environment

Page 21 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46

Industrial & Engineering Chemistry Research

Table 2. Three-fold Cross Validation Results

m

CO2 CO2 CO2 [MIM][BF4] [MIM][BF4] [MIM][BF4]

n

Models

No. of data points in training set

No. of data points in validation set

CH2O Model 1/3 44 CH2O Model 2/3 45 CH2O Model 3/3 45 best values for interaction parameters CH2O Model 1/3 78 CH2O Model 2/3 78 CH2O Model 3/3 78 best values for interaction parameters

23 22 22 39 39 39

amn

anm

195.6120 137.1718 137.0440 137.1718 -42.2773 -36.7135 -38.8464 -42.2773

569.1142 748.9029 749.5391 748.9029 -152.6024 -153.6099 -142.1957 -152.6024

21

ACS Paragon Plus Environment

ARDs for validation set

ARDs for all data

0.0779 0.0732 0.0824

0.0772 0.0768 0.0769

0.0389 0.0461 0.0437

0.0427 0.0427 0.0427

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 36

Table 3. Henry’s Constants of CO2 in Solvents Obtained Experimentally (Hexp) and Predicted by UNIFAC-Lei Model (Hpred) at Different Temperatures Solvents

T (K)

Hexp(MPa)

Hpred(MPa)

RDs

TEG

273.15 293.15 313.15 333.15 353.15

7.61 11.13 17.05 23.37 32.61

7.00 10.97 16.55 24.12 34.13

0.0804 0.0148 0.0293 0.0321 0.0468

[BMIM]+[BF4]-

273.15

4.05

3.63

0.1023

293.15 313.15 333.15 353.15 273.15 293.15 313.15 333.15 353.15 273.15 293.15 313.15 333.15 353.15 273.15 293.15 313.15 333.15 353.15

6.67 9.66 13.49 16.82 4.84 7.03 10.12 14.74 20.02 5.65 8.51 12.20 16.37 22.00 6.48 10.48 15.93 18.89 27.63

5.64 8.46 12.28 17.34 4.19 6.49 9.71 14.07 19.82 5.06 7.86 11.76 17.03 23.99 6.09 9.50 14.28 20.75 29.28

0.1538 0.1249 0.0901 0.0304 0.1335 0.0768 0.0405 0.0455 0.0099 0.1035 0.0764 0.0361 0.0406 0.0901 0.0606 0.0929 0.1036 0.0985 0.0595

TEG(x2=0.2) + ILs(x3=0.8)

TEG(x2=0.5) + ILs(x3=0.5)

TEG(x2=0.8) + ILs(x3=0.2)

22

ACS Paragon Plus Environment

Page 23 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Table 4. Comparison of the Simulation Results among the Three CO2 Dehydration Processes Contents

Absorbents

Product stream

Streams

Solvent

Pure IL

Pure TEG

TEG + IL

Temperature (℃)

25

25

25

Mass flowrate (kg·h-1)

974

980

978

H2O content of the CO2 product (molar fraction)

101 ppm

563 ppm

393 ppm

CO2 loss (kg·h-1)

17.72

11.72

13.72

Mass flowrate (kg·h-1)

3500

3500

3500

CO2 content in the recycled solvent (molar fraction)

384ppm

2000ppm

2282ppm

0.147

-0.008 0.16 0.16

-0.0035 0.152 0.152

1.201

1.362

1.385

Flash drum Heat duty (GJ/h) Condenser (GJ/h) Desorption column Reboiler (GJ/h) Total heating duty (GJ/h) Heat Heat duty (GJ/h) exchanger

Heat dutya

0.147

Cooler1

Heat duty (GJ/h)

-0.146

-0.195

-0.161

Cooler2

Heat duty (GJ/h)

-0.024

-0.0008

-0.0001

-0.17 0.184 3.56

-0.2038 0.214

-0.1646 0.209

Total cooling duty (GJ/h) PSEP (GJ/kg)b Energy consumption of vacuum pump (kW) aCO

2

dehydration process with IL: total heating duty = energy consumption of flash drum, total cooling

duty = energy consumption of cooler 1 + energy consumption of cooler 2; CO2 dehydration process with pure TEG or TEG + IL: total heating duty = energy consumption of reboiler, total cooling duty = energy consumption of cooler 1 + energy consumption of cooler 2 + energy consumption of condenser. bPSEP

=

energy

demand

/

flow

rate

of

23

ACS Paragon Plus Environment

captured

H2O

(GJ/kg).

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure Captions Figure 1. Solubility of CO2 (1) in TEG, [BMIM]+[BF4]-, and their mixtures at 273.15 K (a), 293.15 K (b), 313.15 K (c), 333.15 K (d), and 353.15 K (e). Lines, results predicted by the UNIFAC-Lei model; scattered points, experimental data. (●) and (—), TEG; (□) and (– –), TEG (x2=0.8) + [BMIM]+[BF4]- (x3=0.2); (△) and (– –), TEG (x2=0.5) + [BMIM]+[BF4](x3=0.5); (○) and (– –), TEG (x2=0.2) + [BMIM]+[BF4]- (x3=0.8); (▲) and (—), [BMIM]+[BF4]-.

Figure 2. Solubility of CO2 (1) in the binary mixtures (TEG + [BMIM]+[BF4]-) and in the ternary mixtures (TEG + [BMIM]+[BF4]- + H2O) at 293.15 K (a, d), 313.15 K (b), and 333.15 K (c). Lines, results predicted by the UNIFAC-Lei model; scattered points, experimental data. (●) and (—), TEG (x2=0.14) + [BMIM]+[BF4]- (x3=0.56) + H2O (x4=0.3); (▲) and (—), TEG (x2=0.35) + [BMIM]+[BF4](x3=0.35) + H2O (x4=0.3); (■) and (—), TEG (x2=0.56) + [BMIM]+[BF4]- (x3=0.14) + H2O (x4=0.3); (★) and (—), TEG (x2=0.5) + [BMIM]+[BF4]- (x3=0.5).

Figure 3. (a) The effect of solvent volume flowrate (VS) on the H2O content (mole fraction) in outlet gas (y1) when the WH2O in the feeding solvent is 500 ppm (●), 1000ppm (▲), 2000ppm (■) , respectively; and (b) the effect of the H2O content WH2O (mass fraction) in solvent on the H2O content (mole fraction) in outlet gas (y1) when VS is 5 ml/min. Lines, predicted values by the EQ stage model; scattered points, experimental data.

24

ACS Paragon Plus Environment

Page 24 of 36

Page 25 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Figure 4. Excess enthalpies calculated by the COSMO-RS model at 298.15 K. CO2 (1) + TEG (2) (a); CO2 (1) + [BMIM]+[BF4]- (2) (b); and CO2 (1) + TEG (x2=0.5) + [BMIM]+[BF4]- (x3=0.5) (c).

Figure 5. The intermolecular interactions Hm (a), HvdW (b), HHB (c), and HMF (d) for the excess enthalpies of binary and ternary mixtures at 298.15 K. (▲) [BMIM]+[BF4]-; (□) TEG (x2=0.2) + [BMIM]+[BF4]- (x3=0.8); ( △ ) TEG (x2=0.5) + [BMIM]+[BF4]- (x3=0.5); (○) TEG (x2=0.8) + [BMIM]+[BF4]- (x3=0.2); (●) TEG.

Figure 6. Excess enthalpies (Hm) of equimolar CO2-solvent mixtures calculated by COSMO-RS model versus the Henry’s constants (H) calculated by UNIFAC-Lei model at 298.15 K (a) and the effect of HvdW, HHB, and HMF on CO2 solubility at 298.15 K (b). (▲) [BMIM]+[BF4]-; (□) TEG (x2=0.2) + [BMIM]+[BF4]- (x3=0.8); ( △ ) TEG (x2=0.5) + [BMIM]+[BF4]- (x3=0.5); (○) TEG (x2=0.8) + [BMIM]+[BF4]- (x3=0.2); (●) TEG.

Figure 7. Excess enthalpies of the ternary mixtures: (a) solid lines, H2O + TEG + [BMIM]+[BF4]-; dotted lines, CO2 + TEG + [BMIM]+[BF4]-; black lines, TEG (x2=0.2) + [BMIM]+[BF4]- (x3=0.8); red lines, TEG (x2=0.5) + [BMIM]+[BF4]- (x3=0.5); blue lines, TEG (x2=0.8) + [BMIM]+[BF4]- (x3=0.2). (b) H2O (x1=0.3) + TEG (x2=0.35) + [BMIM]+[BF4](x3=0.35) and CO2 (x1=0.3) + TEG (x2=0.35) + [BMIM]+[BF4]- (x3=0.35).

25

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 8. Color-filled reduced density gradient (RDG) maps. Isovalue of RDG is set to 0.5, and the value of sign (λ2)ρ on the surfaces is represented by filling color ranging from -0.03 to 0.02 au. Blue means the strong attractive interactions, and red means the strong nonbonded overlap.

Figure 9. The CO2 dehydration processes with pure IL (a) and pure TEG or TEG + IL (b) as absorbents.

26

ACS Paragon Plus Environment

Page 26 of 36

2.5

3.0

2.0

2.5 2.0

P (MPa)

1.5 1.0

1.5 1.0

0.5

0.5

(a) 0.0 0.0

0.1

0.2

4.0

x1

0.3

0.4

0.5

0.0 0.0 4.0

3.5

3.5

3.0

3.0

2.5

2.5

2.0

0.1

0.2

x1

0.3

0.4

2.0

1.5

1.5

1.0

1.0

0.5 0.0 0.0

(b)

P (MPa)

P (MPa)

0.5

(c) 0.1

x1

0.2 4.0

0.0 0.0

0.3

(d) 0.1

x1

0.2

3.5 3.0 2.5

P (MPa)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

P (MPa)

Page 27 of 36

2.0 1.5 1.0 0.5

(e)

0.0 0.0

0.1 x1

0.2

Figure 1. Solubility of CO2 (1) in TEG, [BMIM]+[BF4]-, and their mixtures at 273.15 K (a), 293.15 K (b), 313.15 K (c), 333.15 K (d), and 353.15 K (e). Lines, results predicted by the UNIFAC-Lei model; scattered points, experimental data. (●) and (—), TEG; (□) and (– –), TEG (x2=0.8) + [BMIM]+[BF4]- (x3=0.2); (△) and (– –), TEG (x2=0.5) + [BMIM]+[BF4](x3=0.5); (○) and (– –), TEG (x2=0.2) + [BMIM]+[BF4]- (x3=0.8); (▲) and (—), [BMIM]+[BF4]-.

27

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

4.0

5.0

3.0 P (MPa)

P (MPa)

4.0 3.0 2.0

2.0 1.0

1.0

(b)

(a)

0.0 0.0

0.1

0.2 x1

0.3

0.0 0.0

0.4

5.0

5.0

4.0

4.0

3.0

3.0

P (MPa)

P (MPa)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 36

2.0 1.0 0.0 0.0

x1

0.2

2.0 1.0

(c)

0.1 x1

0.1

(d)

0.2

0.0 0.0

0.1

x1

0.2

0.3

Figure 2. Solubility of CO2 (1) in the binary mixtures (TEG + [BMIM]+[BF4]-) and in the ternary mixtures (TEG + [BMIM]+[BF4]- + H2O) at 293.15 K (a, d), 313.15 K (b), and 333.15 K (c). Lines, results predicted by the UNIFAC-Lei model; scattered points, experimental data. (●) and (—), TEG (x2=0.14) + [BMIM]+[BF4]- (x3=0.56) + H2O (x4=0.3); (▲) and (—), TEG (x2=0.35) + [BMIM]+[BF4](x3=0.35) + H2O (x4=0.3); (■) and (—), TEG (x2=0.56) + [BMIM]+[BF4]- (x3=0.14) + H2O (x4=0.3); (★) and (—), TEG (x2=0.5) + [BMIM]+[BF4]- (x3=0.5).

28

ACS Paragon Plus Environment

Page 29 of 36

700

700 (a)

600

600

500

500 y1 (ppm)

y1 (ppm)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

400 300

400 300

200

200

100

100

0

0

5

10

VS (ml/min)

0

15

(b)

0

1000

2000

3000

4000

WH2O (ppm)

Figure 3. (a) Effect of solvent volume flowrate (VS) on the H2O content (mole fraction) in outlet gas (y1) when the WH2O in the feeding solvent is 500 ppm (●), 1000ppm (▲), 2000ppm (■) , respectively; and (b) effect of the H2O content WH2O (mass fraction) in solvent on the H2O content (mole fraction) in outlet gas (y1) when VS is 5 ml/min. Lines, predicted values by the EQ stage model; scattered points, experimental data.

29

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

1500

1000

HCO2-IL(J·mol-1)

HHB

1000 HCO2-TEG(J·mol-1)

500 HMF

0

-500

Hm

-1000

HvdW

-1500 0.0

0.2

0.4

0.6

x1

0.8

500

HHB

0

HMF

-500 Hm

-1000 -1500 -2000

HvdW

-2500

(a)

-3000 0.0

1.0

0.2

0.4

x1

0.6

(b) 0.8 1.0

1000 HCO2-IL+TEG (J·mol-1)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 36

HHB

500 0

HMF

-500 Hm

-1000 -1500 -2000

HvdW

-2500 0.0

0.2

0.4

x1

0.6

(c)

0.8

1.0

Figure 4. Excess enthalpies calculated by the COSMO-RS model at 298.15 K. CO2 (1) + TEG (2) (a); CO2 (1) + [BMIM]+[BF4]- (2) (b); and CO2 (1) + TEG (x2=0.5) + [BMIM]+[BF4]- (x3=0.5) (c).

30

ACS Paragon Plus Environment

Page 31 of 36

500

0 -500

Hm (J·mol-1)

HvdW (J·mol-1)

0

-1000

-500

-1500

-1000

-2000

(b)

(a) -1500 0.0

0.2

0.4

x1

0.6

0.8

-2500 0.0

1.0

1400 1200

0.4

0.2

0.4

x1

0.6

0.8

1.0

0.6

0.8

1.0

(d)

(c) 400

HMF (J·mol-1)

800 600 400

300 200 100

200 0 0.0

0.2

500

1000

HHB (J·mol-1)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

0.2

0.4

x1

0.6

0.8

1.0

0 0.0

x1

Figure 5. The intermolecular interactions Hm (a), HvdW (b), HHB (c), and HMF (d) for the excess enthalpies of binary and ternary mixtures at 298.15 K. (▲) [BMIM]+[BF4]-; (□) TEG (x2=0.2) + [BMIM]+[BF4]- (x3=0.8); ( △ ) TEG (x2=0.5) + [BMIM]+[BF4]- (x3=0.5); (○) TEG (x2=0.8) + [BMIM]+[BF4]- (x3=0.2); (●) TEG.

31

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

(a)

0

Hm (J·mol-1)

-200 -400 -600 -800

-1000 -1200

6

8

10

12

14

IL

0.2TEG+0.8IL

0.5TEG+0.5IL

1200 900 600 300 0 -300 -600 -900 -1200 -1500 -1800

0.8TEG+0.2IL

(b)

TEG

H (MPa)

Hm (J·mol-1)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 36

■ HMF ■ HHB ■ HvdW

12.20

10.91

8.72

7.20

6.26

H (MPa)

Figure 6. Excess enthalpies (Hm) of equimolar CO2-solvent mixtures calculated by COSMO-RS model versus the Henry’s constants (H) calculated by UNIFAC-Lei model at 298.15 K (a) and the effect of HvdW, HHB, and HMF on CO2 solubility at 298.15 K (b). (▲) [BMIM]+[BF4]-; (□) TEG (x2=0.2) + [BMIM]+[BF4]- (x3=0.8); ( △ ) TEG (x2=0.5) + [BMIM]+[BF4]- (x3=0.5); (○) TEG (x2=0.8) + [BMIM]+[BF4]- (x3=0.2); (●) TEG.

32

ACS Paragon Plus Environment

Page 33 of 36

2000

Hm (J·mol-1)

1500 (a) 1000 500 0 -500

-1000 -1500 -2000 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 xCO2 (xH2O)

1200 800 Hm (J·mol-1)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

■ HMF

(b)

■ HHB ■ HvdW

400 0

-400 -800

-1200 CO2

H2O

Figure 7. Excess enthalpies of the ternary mixtures: (a) solid lines, H2O + TEG + [BMIM]+[BF4]-; dotted lines, CO2 + TEG + [BMIM]+[BF4]-; black lines, TEG (x2=0.2) + [BMIM]+[BF4]- (x3=0.8); red lines, TEG (x2=0.5) + [BMIM]+[BF4]- (x3=0.5); blue lines, TEG (x2=0.8) + [BMIM]+[BF4]- (x3=0.2). (b) H2O (x1=0.3) + TEG (x2=0.35) + [BMIM]+[BF4](x3=0.35) and CO2 (x1=0.3) + TEG (x2=0.35) + [BMIM]+[BF4]- (x3=0.35).

33

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 36

(b) TEG + [BMIM]+[BF4]- + CO2

(a) TEG + [BMIM]+[BF4]-

(c) TEG + [BMIM]+[BF4]- + H2O -0.02

-0.03

H-bond

vdW

Steric

Figure 8. Color-filled reduced density gradient (RDG) maps. Isovalue of RDG is set to 0.5, and the value of sign (λ2)ρ on the surfaces is represented by filling color ranging from -0.03 to 0.02 au. Blue means the strong attractive interactions, and red means the strong nonbonded overlap. 34

ACS Paragon Plus Environment

Page 35 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

H2O CO2 product gas IL

Cooler 2 Cooler 1

Heat exchanger Flash

Feed

Pump

(a) Offgas TEG or TEG+IL H2O

CO2 product gas Cooler 2 Cooler 1

Desorption column

Heat exchanger Reboiler

Feed

Pump

(b) Figure 9. The CO2 dehydration processes with pure IL (a) and pure TEG or TEG + IL (b) as absorbents. 35

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table of Content (TOC) Graphic

36

ACS Paragon Plus Environment

Page 36 of 36