Thermodynamic Modeling and Assessment of Ionic Liquid-Based CO2

Jun 17, 2014 - Dawei Shang , Xinyan Liu , Lu Bai , Shaojuan Zeng , Qiuxia Xu ... Jun Yin , Jie Gao , Ming Tong , Xin Chen , Wanzhong Kang , Yanbo Zhou...
0 downloads 0 Views 1003KB Size
Subscriber access provided by Maastricht University Library

Article

Thermodynamic Modeling and Assessment of Ionic Liquid-based CO2 Capture Processes Ying Huang, Xiangping Zhang, Xin Zhang, Haifeng Dong, and Suojiang Zhang Ind. Eng. Chem. Res., Just Accepted Manuscript • Publication Date (Web): 17 Jun 2014 Downloaded from http://pubs.acs.org on June 18, 2014

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 free 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 accessible to all readers and 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.

Industrial & Engineering Chemistry Research 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 43

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

Thermodynamic Modeling and Assessment of Ionic Liquid-based CO2 Capture Processes Ying Huang1,2, Xiangping Zhang1,*, Xin Zhang1,2, Haifeng Dong1,2, Suojiang Zhang1,* 1

Beijing Key Laboratory of Ionic Liquids Clean Process, State Key Laboratory of Multiphase

Complex Systems, Key Laboratory of Green Process and Engineering, Institute of Process Engineering, Chinese Academy of Sciences, Beijing 100190, China 2

School of Chemistry and Chemical Engineering, University of Chinese Academy of Sciences,

Beijing 100049, China ABSTRACT Ionic liquid (IL)-amine hybrid solvents has been experimentally proved to be effective for CO2 capture. This paper provided rigorous thermodynamic models, process simulation and cost estimation of a potential design of IL-based CO2 capture processes. Three ILs ([Bmim][BF4], [Bmim][DCA] and [Bpy][BF4]) were investigated to blend with MEA aqueous solution. The physicochemical properties of the ILs were predicted by several temperature-dependent correlations. Phase equilibria were modeled based on Henry’s law and NRTL equation, and the calculated values were in good agreement with the experimental data. The simulation results show that the [Bpy][BF4]-MEA process can save about 15% regeneration heat duty compared to the conventional MEA process, which is attributed to the reduction of sensible and latent heat. Moreover, a modified [Bpy][BF4]-MEA process via adding inter-cooling and lean vapor recompression presents 12% and 13.5% reduction in overall equivalent energy penalty and capture cost compared with the conventional MEA process, respectively.

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

Page 2 of 43

1. INTRODUCTION Carbon capture and storage (CCS) is recognized as an option to reduce the greenhouse gas emissions1, including three typical technologies, i.e., pre-combustion capture, post-combustion capture and the oxyfuel process2. Among which post combustion capture via amine absorption such as monoethanolamine (MEA) has been commercially used for decades3-5 with several disadvantages, including the loss of solvent, large stripping energy consumption, unfavourable degradation6 and corrosion problems. Thus it is highly required to develop new solvents and processes. In recent years, ionic liquids (ILs) have been paid more attentions7-12 because of their nonvolatility, thermal stability and tunable chemistry, which are perceived to be helpful in reducing energy consumption and environmental pollution. Conventional ILs such as 1-butyl-3-methylimidazolium

tetrafluoroborate

([Bmim][BF4])13

and

1-butyl-3-methylimidazolium bis((trifluoromethyl)sulfonyl)imide ([Bmim][Tf2N])14 present relatively well CO2 solubilities under high pressure, thus they are not suitable for CO2 capture from flue gas under atmospheric pressure. Task-specific ILs with amine groups can chemically absorb up to 1 mole CO2 per mole IL15. However, it may be difficult to employ the task-specific ILs in industrial applications because of complex synthesis and purification steps, high cost16 and increased viscosity after the absorption reaction17. Alternatively, IL-amine hybrid solvents seem to be effective for CO2 capture, such as IL-MEA18 and IL-MDEA19-21 aqueous solutions. They are prior to current amine absorption in terms of reducing solvent volatility, energy consumption and the corrosion rate22, 23. Recently, Yang et al.24 conducted a CO2 capture experiment in an absorption−desorption loop system using an IL-MEA hybrid 2

ACS Paragon Plus Environment

Page 3 of 43

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

solvent (30 wt % MEA + 40 wt % [Bmim][BF4] + 30 wt % H2O). They observed that the MEA loss for the mixed solution was much lower than that for the MEA solution. Besides, no ionic liquid loss was detected. These results proved the possibility and feasibility of IL-amine solvents for CO2 capture, including the CO2 removal efficiency and the solvents stability. LaFrate et al.

25

detected the degradation of two IL-MEA solvents. The results revealed that

[Emim][OTf]

(1-ethyl-3-methylimidazolium

trifluoromethanesulfonate)

make

the

degradation of the solvent faster, which was induced by amine alkylation by the anion([OTf]). Moreover, Janiczek et al.26 set up a technical-scale-plant for CO2 absorption from natural gas streams utilizing an imidazolium-based IL, which indicated the potential of industrial application of ILs. Process simulation and assessment is essential for a new process to be developed in industry. However, such research related to IL-based processes was rarely reported probably due to the shortage of rigorous thermodynamic models for complex IL-containing systems27. Shiflett et al.28 simulated a capture process using pure [Bmim][Ac] and claimed a 16% reduction in energy consumption compared to the commercial MEA process. Basha et al.29 developed a conceptual CO2 capture process applying three flash drums as regenerators. They considered [Hmim][Tf2N] could be used as a physical solvent for high efficient CO2 capture from warm shifted fuel gas streams with relatively high CO2 content (23.87 mol.%). The above pioneering works provide ideas and methods for the simulation of CO2 capture process using pure IL, while they may be not available for flue gas with lower CO2 concentration. Thus it is imperative to simulate and assess the capture processes using IL-amine hybrid solvents. Besides developing new solvents, efforts were also concentrated on the process 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

Page 4 of 43

modification to improve the process performance. Numerous modifications can be made on the basis of a standard process such as split flow, rich split, vapor recompression, and absorber inter-stage cooling30. Karimi et al.31 found that vapor recompression was the best configuration with the lowest increase in total capture cost and the plant complexity. In general, in order to achieve the energy benefits, the capital cost and the complexity of the plant will probably increase32. So the satisfactory process should the one with the best balance between additional cost and decreased energy penalty33. Taken together, it is meaningful to employ an IL-amine hybrid solvent in a modified process configuration to develop an efficient capture process. The aim of this work is to provide thermodynamic models and system assessment method for new ionic liquid-based CO2 capture processes. Ionic liquids were created as new (pseudo) components in Aspen Plus (version 8.4), requiring the import of the basic properties and parameters such as critical properties and heat capacities, and the parameters of thermodynamic models which were regressed via experimental data from literature and our study. The rigorous thermodynamic models of the MEA-H2O-IL systems were verified to ensure the reliabilities of the simulation results. Based on the results of process simulation, the energetic analysis and cost estimation of IL-based processes were performed, taking the conventional MEA-based process as the reference case. 2. SELECTION OF ILS The anion of an ionic liquid has a greater influence on the CO2 solubility than the cation34. In this work, three hydrophilic ionic liquids ([Bmim][BF4] , [Bmim][DCA] and [Bpy][BF4]) were selected as component of the IL-amine hybrid solvents based on the following reasons. [Bmim][BF4] has been studied frequently for mixing with MEA. Taib et al.18 measured CO2 4

ACS Paragon Plus Environment

Page 5 of 43

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

solubilities in [Bmim][BF4]+MEA solution at pressure from 100 to 1600 kPa and observed that the solubilities were higher than conventional MEA solution. Based on a kinetic experiment, Lu et al.35 found that the values of the enhancement factor and the second-order reaction rate constant for the [Bmim][BF4]+MEA solution were higher than those in pure MEA aqueous. [Bmim][DCA] has very low viscosity (16.8mPa·s at 40℃ 36) which is favorable for the absorption process. Experiments results including the CO2 solubilities, density and viscosity of the [Bmim][DCA]-MEA hybrid solvent have been reported in our previous work37. [Bpy][BF4] is perceived to be more potential for large scale utility with lower cost38, toxicity39 and higher biodegradability40 than imidazolium-based ILs. 3. METHODOLOGY A four-step procedure of this study was conducted: 1) The critical properties of the ILs were calculated based on the FC-CS method41 proposed in our previous work. The temperature-dependent properties, including densities, viscosities, surface tensions, heat capacities and thermal conductivities were predicted and compared with the experimental data from literature. 2) Gas-liquid equilibrium (GLE) of CO2-IL was correlated using the experimental data from literature. Vapor-liquid equilibria (VLE) of the IL-MEA-H2O systems were measured for the thermodynamic study of the whole absorption system. 3) Four IL-based CO2 capture processes and the reference MEA process were simulated, including an IL-based process with configuration modification. 4) The MEA and IL-based processes were assessed and compared on the aspects of energy consumption and capture cost. 3.1 Experimental section 3.1.1 Materials.

All the above ILs were supplied by Linzhou Keneng Materials 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 43

Technology Co., Ltd. China and dried under vacuum for 48 hours at 333 K before used. The water contents of the ILs after drying were ensured <200 ppm by Karl-Fisher titration. The residual chloride in the ILs were < 0.002mol·L-1, which were measured by the ion chromatography (Dionex ICS-900). MEA (purity>99.8%) was purchased from Beijng Chemical Works and used without further purification. The deionized water was supplied by the local laboratory. 3.1.2 Apparatus and Operation.

The experimental apparatus was based on the

quasi-static ebulliometry method which was the similar to some researchers 42,43-45, as shown in Figure 1. It was composed by two ebulliometers (a working ebulliometer and a reference ebulliometer) connected to a buffering vessel, two condensers to minimize the composition variation, two temperature control and measurement systems to make the liquid boiling and record the temperature, two magnetic stirrers, and a pressure transducer system. The total variation of the liquid phase composition caused by the holdup amount of volatile component in the condenser was within 0.2%42. First, the hybrid solvents were obtained by mixing definite weight of the corresponding components weighted by an electronic analytical balance (BS124S, Sartorius Scientific Instrument Co. Ltd., China) with precision of ± 0.0001 g. Second, the sample was added into the working ebulliometer and the deionized water was added into the reference ebulliometer. Third, the system was evacuated to an absolute pressure of approximately 7 kPa, and then the solution was heated and stirred. When VLE was reached, the temperature of two ebulliometers was unchanged and recorded. A series of VLE data were obtained by increasing the system pressure. 6

ACS Paragon Plus Environment

Page 7 of 43

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.2 Prediction of Physicochemical Properties The physicochemical properties, phase equilibrium and chemical reaction of the MEA-H2O-CO2 system were modeled based on the electrolyte nonrandom two-liquid model. This property method has been commonly adopted for simulating the MEA scrubbing process27, 46-48

. For ionic liquids, the property prediction and phase equilibrium modeling were described

in detail in the subsequent sections. 3.2.1 Scalar Properties.

Critical properties are used in many corresponding states

correlations for volumetric, thermodynamic, and transport properties of gases and liquids49. They are also essential for a process simulation. Critical properties of the ILs were obtained from the FC-CS method reported in our previous work41. Other properties were obtained from literature50,51,52, 53,54. 3.2.2 Temperature-dependent Properties.

The temperature-dependent properties such as

heat capacity, density and viscosity were calculated by the equations (Eq.1-Eq.6) listed in Table 1. Ionic liquid is well-known as non-volatile liquid whose vapor pressure is difficult to be observed. Liquid heat capacity is a basic thermodynamic property used for specifying the amount of heat required to change the temperature of a liquid by a given amount, which is important in analyzing energy demand of the whole process. All the equation coefficients were obtained by minimizing the objective function (Eq. 7) based on available experimental data. The properties were calculated by the FC-CS method41 if there was no available experimental data. NP

O.F . = ∑ ( X i cal − X i exp )

2

(7)

i =1

The average absolute relative deviation (AARD) was defined as follows: 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

Page 8 of 43

NP

AARD(%) = 100 × ∑ X i cal / X i exp − 1.0 / N P

(8)

i =1

where X represents the temperature-dependent property, subscript ‘i’ denotes different investigated IL, Np denotes the total number of data points, and the superscripts ‘cal’ and ‘exp’ are the calculated values and the experimental data, respectively. 3.3 Phase Equilibria Modeling The phase equilibrium and kinetic reactions among the MEA-H2O-CO2 system were modeled using a so-called KEMEA thermodynamic package, which is available within the Aspen Plus modeling system based on the electrolyte nonrandom two-liquid (ELECNRTL) model. This property method has been commonly adopted for simulating the MEA scrubbing CO2 process by many researchers46-48. IL was created as a new component whose interaction parameters with other components were required and determined as the following described. Solubilities of N2 and O2 in ILs were very low55 at the magnitude of 10-4. Thus N2 and O2 were assumed as no solubility in ILs in the simulation. Moreover, their solubilities in MEA and H2O were also very low which have been included in the ‘KEMEA’ property package. 3.3.1 GLE of the CO2-IL Binary System.

The phase equilibrium relationship for

dissolved gases was modeled as following56:

ϕiV yi p = xiγ i* H iA

(9)

γ i* = (γ i / γ i∞ )

(10)

where γ i∞ is the infinite dilution activity coefficient of component i in the mixture. The Henry’s law constant of gas component i was defined as

fi ϕiV yi p H iA = lim = lim xi → 0 x xi →0 xi i

(11)

8

ACS Paragon Plus Environment

Page 9 of 43

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

where p, yi, and xi are total pressure, mole fraction of component i in vapor phase and liquid phase, respectively. ϕ iV is the fugacity coefficient in vapor phase. HiA is the Henry constants in Pa of component i in solvent A. The temperature dependence of the Henry constants was described by the following equation56:

ln H iA = aiA + biA / T + ciA ln T + d iAT

(12)

where aiA, biA, ciA and diA are the equation parameters. Due to the negligible vapor pressure of ILs, the gaseous phase was assumed to be pure CO2. The fugacity coefficient in vapor phase ϕ iV can be obtained by Redlich–Kwong (R–K) equation of state57, given in the following form.

ln ϕ = Z − 1 − ln( Z − bp / RT ) − (

a / R 2T 2.5 ) ln(1 + bp / ZRT ) b / RT

Z = pV / RT

p=

RT a − 0.5 V − b T V (V + b) a = 0.42748

(13) (14) (15)

R 2Tc2.5 pc

(16)

RTc pc

(17)

b = 0.08664

where ϕ is the fugacity coefficient, Z is the compressibility factor, R is the gas constant, a and b are the EOS constants, Pc and Tc are the critical pressure and critical temperature, respectively, and V is the molar volume. * Based on the calculated H CO2 and γ CO , the activity coefficient γ CO2 in the mixture were 2

modeled by the NRTL model58: 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

δ

lnγ i =

∑ xτ j =1

j

ji

G ji

δ

∑xG k =1

k

ki

 δ  x j Gij + ∑ δ j =1   ∑ xk Gkj  k =1

  τ ij −   

δ

Page 10 of 43

 Gmj     xk Gkj   

∑xτ

m =1

δ

∑ k =1

m mj

(18)

Gij = exp( −α ijτ ij ) ; τ ij = ( gij − g jj ) / RT = aij + bij / T ; α ij = α ji (i ≠ j ) where δ is the number of components; x represents the mole fraction; R is the gas constant; T is the absolute temperature; G ij is a dimensionless interaction parameter depending on the energy interaction parameter (gij) and the nonrandomness factor (αij). α ij was fixed at 0.3 used in most polar systems. The parameters of Henry constants and NRTL binary interaction parameters between CO2 and IL were generated based on the above models using experimental solubility data from literature.

3.3.2 VLE of the MEA-H2O-IL System.

The following equation (modified Raoult’s law)

was used for VLE calculation at the pressures less than 2 bars.

yi p = xiγ i pis

(19)

where p and pis are total pressure of the system and vapor pressure of the pure component i, respectively. yi and xi represent mole fraction of component i in the vapor phase and liquid phase, respectively. Due to the vapor mole fraction of IL was assumed as zero, Eq. 19 can be simplified to Eq. 20 and Eq. 21 for binary system and ternary system, respectively.

p = x1γ 1 p1s

(20)

p = x1γ 1 p1s + x3γ 3 p3s

(21)

where subscript ‘1’ denotes H2O, ‘3’ denotes MEA.

10

ACS Paragon Plus Environment

Page 11 of 43

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

Vapor pressure of pure H2O and MEA can be calculated by the Antoine Equation. Based on the vapor pressure of liquid mixture of different composition, the phase equilibrium criterion (Eq. 20 and Eq.21) combining with the NRTL model (Eq. 18) were used to calculate binary parameters of IL-H2O and IL-MEA systems.

3.4 Process Assessment 3.4.1 Process simulation.

The MEA and IL-based processes were all simulated in Aspen

Plus, among which the columns were modeled by RadFrac (an equilibrium stage model). The flue gas specification from the Shanxi Coke Plant in China was as following, temperature: 200 ℃ ; pressure: 1.01bar; mole flow rate: 24000kmol/h; mole composition: 6.37% CO2, 69.46% N2, 3.66% O2, 20.51% H2O. The plant produces methanol from CH4 reforming, thus the captured CO2 can be utilized as an additive to adjust the H/C ratio to an appropriate value for methanol production, which was perceived to be beneficial both in economy and environment, similar to the idea of Taghdisian et al.59. The flue gas has been desulphurized before it is sent to the CO2 capture unit. Figure 2 illustrates the conventional capture process, including the following main parts: (1) Gas pretreatment. The flue gas is pressurized to 1.1 bar using a blower and cooled down to approximately 308K after passing the washing column. (2) CO2 removal step. The cooled flue gas then contacts the solvent in a packed counter-current absorber. The purified gas is vented to the atmosphere after a MEA recovery section at the top of the absorber. (3) Heat exchanger section. The rich solvent from the absorber bottom is preheated by the lean solvent from the stripper bottom via a lean/rich heat exchanger and then pumped into the stripper. (4) Solvent regeneration step. The preheated rich solvent is regenerated in the stripper to release CO2. CO2 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

Page 12 of 43

and steam are cooled in a condensate reflux section at the top of the column, with condensed water being recycled to the top of the stripping column. The lean solvent from the bottom is cooled to 313K by passing the lean/rich heat exchanger and lean solvent cooler before recycling to the absorber. (5) Compression step. CO2 is dried and compressed to 2.5 MPa to enter into the methanol production plant. A modified capture process adding an absorber inter-cooling (AIC) and a lean vapor recompression (LVR) is depicted in Figure 3. Such modifications have been proved to be efficient in reducing energy consumption both by experiment and simulation31,60. The AIC section draws a side stream from the absorber and feeds it to the intercooler to cool down to a specific temperature. The cooled liquid is returned to the column at a location just beneath where the liquid is taken. Lower temperature of the solvent may be benefit in increasing absorption capacity thus decrease the solvent demand. In LVR scheme, semi-lean solvent is extracted from the stripper and flashed. The flashed lean vapor is recompressed to the stripper to substitute part of steam. The stage numbers of absorber and stripper were determined by sensitivity analysis depicted in Supporting Information. Murphree efficiency of the absorber was set as 25% (typically Murphree efficiency can be set as 10-30%)61. The flow rate of the inlet washing water of the MEA recovery section was varied to maintain the water balance of the system, and to lower the MEA slip to an acceptable value of 1 ± 0.05 ppmV. The equipment specifications of the modified process were set referring to literature31,62 as shown in Table S11 of Supporting Information.

12

ACS Paragon Plus Environment

Page 13 of 43

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.4.2 Energetic analysis.

The specific regeneration heat duty is often an indicator of

process performance, but it is not sufficient to rank different capture process. Thus the overall energy consumption was expressed by the equivalent energy penalty considering both electrical power and thermal energy. As adopted by other researchers31,

63

, the turbine

efficiency to produce electrical power from steam was assumed to be 75%, the temperature of steam was assumed to be 10 K higher than the reboiler temperature (Treb) and that steam condensed at 313 K in the turbine.

 313  Weq = 0.75Qreb 1 −  + Wcomp + Wpump  Treb + 10  3.4.3 Cost estimation methodology.

(22)

There are several metrics to measure the cost of CCS

system, including the cost of CO2 avoided, cost of CO2 captured, cost of CO2 abated (or reduced), and the increased cost of electricity64. In this study, cost of CO2 captured was selected. The total capture cost (TCC) was calculated from the sum of annual capital cost (ACC) and total operating cost (TOC) based on Table 2 which were common used in the literature33, 65. The purchased equipment cost (PEC) for columns, heat exchangers and pumps were obtained according to the NETL report66. Cost for those equipment out of the capacity scale given in the report were computed using the six tenths rule67 as depicted in Eq. 23. All the equipment costs were updated to the year 2008 with the Chemical Engineering Cost Index. Then TCI was computed on the basis of Table S1 in Supporting Information, referring to Abu-Zahra et al.65 and Schach et al.33

S  CB = C A  B   SA 

0.6

(23)

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

Page 14 of 43

where CB is the cost of equipment having size SB, CA is the known cost of equipment having corresponding size SA (same units as SB), and SB/SA is the ratio known as the size factor, dimensionless. Annual capital cost (ACC) was computed by Eq. 24, which converted TCI to a constant series of payments for every year of the project life.

ACC =

TCI ((1 + ir ) − 1) / ir (1 + ir )n n

(24)

Where ir denotes the interest rate and n denotes project life time. The total operating cost was divided into two parts. One was variable operating cost from MEA make up and utility consumption such as steam, cooling water and electricity etc. Another was fixed operating cost, mainly composed by operating labor, plant overheads, research and development (R&D) cost etc. Detailed computation basis was shown in Table S2 in Supporting Information. The utility price and calculation percentage for fixed operating cost (FOC) were taken from Mores et al.68, Peters et al.69 and Karimi et al.31. The price of MEA was 1250 US$/t68, while The price of IL was estimated at 6600$/t under the assumption of industrial production in the future based on the estimation by Linzhou Keneng Materials Technology Co., Ltd. China.

4. RESULTS AND DISCUSSION 4.1Validation of Thermodynamic Property Models Scalar properties from literature were listed in Table S4 of Supporting Information. The equation coefficients for vapor pressures, heat capacities, densities, viscosities, surface tensions and thermal conductivities were listed in Table S3 of Supporting Information. As shown in Figure 4, the predicted value for these properties display high accordance with the 14

ACS Paragon Plus Environment

Page 15 of 43

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

experimental data with the AARD less than 0.5%.

4.2 Phase Equilibrium Calculation Results Calculated Henry constants were listed in Table 3. NRTL binary interaction parameters of CO2, H2O and MEA with different IL were obtained as shown in Table 4. As shown in Figure 5, the calculated CO2 solubilities agree well with the experimental solubilities with an AARD of 0.9%. Experimental and calculated vapor pressure data can be found in Tables S5−S10 of Supporting Information with an overall AARD of 0.8%. The good accordance between experimental and calculated values in Figures 6−8 validate the applicability of the binary interaction parameters for all systems investigated. All the experimental results indicate that ILs decrease the vapor pressures of solvents due to diluting effect and affinity to H2O. So IL can improve the boiling temperature of the solvent thus increase the reboiler temperature of the IL-based processes.

4.3 Comparison of the MEA and IL-based Processes The following processes were simulated: (1) reference case of 30wt % conventional MEA process; (2) three conventional process configuration using IL-MEA solvents ([Bmim][BF4]-MEA, [Bmim][DCA]-MEA and [Bpy][BF4]-MEA). After that the IL-based solvent with lowest regeneration heat duty was used in the modified process. The main simulation results to assess a CO2 capture process were summarized in Table 5. The following discussion compared the [Bpy][BF4]-based processes with the MEA process in term of lean loading, solvent specific demand, cooling water demand, energy consumption and capture cost.

4.3.1 Lean loading (αlean).

Lean loadings of the MEA process and the IL-based 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

Page 16 of 43

processes were studied to achieve the minimal regeneration heat duty. At the same time, the lean solvent circulation rate was varied to achieve the same CO2 removal capacity. As shown in Figure 9, regeneration heats of different processes present the similar variation with the increase of αlean. The minimum regeneration heat duty can be reached at αlean of 0.25 for the MEA process, 0.22 for the [Bmim][DCA]-MEA process, 0.2 for the [Bmim][BF4]-MEA and the [Bpy][BF4]-MEA process, which were defined as the optimum lean solvent loadings. The reboiler heat duty of 3.7 GJth/t CO2 for the MEA process from this work is close to the value reported by many researches62, 70, 71, which indicates the simulation measures and results are valid. The [Bpy][BF4]-MEA hybrid solvent was used in the modified IL-based process due to the lowest heat duty. The main reasons for such changing trend are as follows. At low values of αlean, the latent energy of generating steam is dominant in the overall heat consumption. With the increase of

αlean, the latent heat decreases till constant, while the sensible heat increases and begins to play a key role in the total regeneration heat consumption due to higher solvent flow rate required for a fixed CO2 removal ratio. According to the experimental results by Kim et al.72, heat of absorption almost keeps constant within the αlean range of 0.16 to 0.5. Thus the minimum heat duty can be reached at the coordinate of sensible heat and latent heat.

4.3.2 Solvent Specific Demand.

Higher value of ∆CO2 loading results in a significant

decrease of solvent specific demand. The solvent circulation rates of the [Bpy][BF4]-MEA process and the modified [Bpy][BF4]-MEA decrease by 22% and 28% compared to the MEA process, which result in lower regeneration sensible heat. Besides, the powers of solvent pumps reduce thus saving a small amount of electrical power. 16

ACS Paragon Plus Environment

Page 17 of 43

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

4.3.3 Cooling Water Demand.

The cooling water demand of the [Bpy][BF4]-MEA

process is less than that of the MEA process. The reduction of circulating solvent flow rate contributes to the reduction of cooling water demand of MEA recovery section and lean cooler. The reduction of lean solvent temperature out of lean/rich exchanger also cut down the cooling water demand of lean cooler. The reduction of vaporized water causes the decrease of cooling water demand of stripper condenser. Ultimately, the cooling water specific demand of the [Bpy][BF4]-MEA process decreases by about 20% compared to the MEA process. Although the modified [Bpy][BF4]-MEA process require more cooling water than the [Bpy][BF4]-MEA process due to the additional inter-cooler, the amount is still less than the reference MEA process.

4.3.4 Energy Consumption.

Regeneration heat duty and total equivalent energy penalty

of the MEA and IL-based processes are described in Figure 10. It is observed that the IL-based processes require less regeneration heat duty than the MEA process, which is attributed to the decrease of sensible and latent heat during regeneration. The sensible heat is decreased because of two factors: (1) The heat capacity of the IL-MEA hybrid solvent is lower than that of 30% MEA solution because the heat capacity of IL is obviously lower than water; (2) The solvent flow rate of the IL-MEA hybrid solvent is lower than MEA solution, which is resulted from higher value of ∆CO2 loading determined by the lean loading optimization. The latent heat decreases due to the reduction of the amount of water vaporized, which is achieved by using IL to reduce vapor pressure of the solvent or increase the reboiler temperature. The modified [Bpy][BF4]-MEA process cut down thermal energy by 31% than the conventional MEA process, because LVR section can provide additional stripping steam to reduce the latent heat, 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

Page 18 of 43

and AIC section decrease the solvent demand to reduce the sensible heat. Besides the thermal energy consumption, the electrical power required for the [Bpy][BF4]-MEA process is a little lower than the MEA process because of lower pump power. However, power consumption of the modified [Bpy][BF4]-MEA process is higher than the MEA process because of the additional compressor. All in all, the [Bpy][BF4]-MEA process and the modified [Bpy][BF4]-MEA process saved the overall 5.8% and 12% of the equivalent energy compared to the reference MEA process, which indicate the IL-based processes are energy-saving compared to the MEA process.

4.3.5 Cost analysis.

The differences of calculation criterion, process configuration,

disposal capacity and CO2 concentration should be considered when comparing the capture cost with other researches. These will lead to different capture cost results, for instance, 55$ by Hassan et al.73, 74€ by Raynal et al.71 and 112$ by Mores et al68. The capture cost of 70$ per ton CO2 by this work is reasonably among the above values. The total capital investment (TCI) is depicted in Figure 11. It is observed that the fixed capital investment (reflecting equipment cost) of the [Bpy][BF4]-MEA process is slightly lower than the MEA process, which is resulted from the reduced cost of pumps and heat exchangers. Even so, TCI of the [Bpy][BF4]-MEA process is a little higher than the MEA process due to increased solvent investment. For the modified [Bpy][BF4]-MEA process, each item of the capital investment is the highest among the three processes, which can be attributed to solvent investment and high purchased equipment cost due to additional flash tank and compressor. Figure 12 presents the distribution of total capture cost (TCC) per ton CO2 captured. It can be firstly observed that total operating cost occupy more than 70% of the total CO2 capture cost. 18

ACS Paragon Plus Environment

Page 19 of 43

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

The main sinks for electrical power cost are the compressor, flue gas inlet blower and pumps. Steam is consumed as a heat source during regeneration. The above costs are the dominating contributors (above 80%) to the total operational cost (TOC). This illustrates the importance of reducing regeneration energy. With comparison to the reference MEA process, the [Bpy][BF4]-MEA process and the modified [Bpy][BF4]-MEA process offer steam cost reduction of 21.8% and 36.1%, respectively. However, the total capital investments (TCI) of IL-based processes are larger than MEA process mainly due to the high cost of IL. Even so, the overall capture cost of IL-based and modified IL-based process are reduced by 11% and 13.5%, respectively compared to the MEA process. This can be illustrated by the low percentage of solvent cost in the total capture cost (1.3% and 1.2% for the IL-based and modified IL-based process, respectively). In general, the cost analysis results indicate the IL-based processes are cost-efficient.

5. CONCLUSIONS The present work assessed the techno-economic performance of IL-based CO2 capture processes based on the rigorous thermodynamic modeling and process simulation results, taking a conventional MEA-based process as the reference case. The thermodynamic models related to the IL-based capture system, including the physicochemical properties and phase equilibria were developed. Phase equilibria of the CO2-IL binary systems were modeled by the R-K equation of state, Henry's law and NRTL model. Experiments were conducted to obtain the VLE data of the H2O-IL-MEA systems. The calculated results displayed good accordance with the experimental data, which indicated the models and parameters were reliable to be used in process simulation. 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 43

The [Bpy][BF4]-MEA process display 15% and 11% reduction in regeneration heat duty and capture cost respectively compared to the reference MEA process. Similarly, the modified [Bpy][BF4]- MEA process with inter-cooling and lean vapor recompression can accomplish 31% and 13.5% reduction in regeneration heat consumption and capture cost, respectively. Although ionic liquid has high price, total capture cost of the IL-based processes are still reduced because solvent investment merely occupy less than 1.5% of the total capture cost. Detailed calculation of energy distribution explains that the reduction in regeneration heat duty is resulted from the contribution of both sensible and latent heat. Besides, the above two IL-based processes are better than the MEA process in term of circulated solvent, cooling water demand and total equivalent energy penalty. All in all, the IL-based solvent combining with the process modification could realize energy-saving and cost-efficient carbon capture, which provides a perspective capture technology in the future.

ASSOCIATED CONTENT Supporting Information Supporting Information provides the following tables and discussion: the calculation method and basis of capital and operation cost, obtained parameters of correlations for physicochemical properties prediction, experimental and calculated vapor pressures of three IL-MEA hybrid solvents, determination of stage numbers of the absorber and stripper, column profiles of the absorber and stripper of the MEA and IL-based processes, mass and energy balance information of the process streams. This information is available free of charge via the internet at http://pubs.acs.org. 20

ACS Paragon Plus Environment

Page 21 of 43

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

AUTHOR INFORMATION Corresponding Author * E-mail: [email protected] and [email protected] Tel./fax: +86-010-82544875.

ACKNOWLEDGEMENTS This research was supported by the Key Program of National Natural Science Foundation of China (21036007), the National Basic Research Program of China (973 Program) (No. 2013CB733506) and the National Natural Science Foundation of China (51274183, 51306194). We also appreciate the support from the Science and Technology Innovation Team of Cross and Cooperation of Chinese Academy of Sciences.

NOTATIONS List of symbols pis

P T Cp Z i*,RA

Z Ai Bi Ci Vi Tr M Np Pc Tb Tbr Tc Vc fi HiA

vapor pressure of pure component i pressure, Pa Temperature, K liquid heat capacity, J·mol-1·K-1 parameters of racket liquid molar volume model compressibility factor the first coefficient of Andrade liquid viscosity equation the second coefficient of Andrade liquid viscosity equation the third coefficient of Andrade liquid viscosity equation liquid molar volume, cm3·mol-1 reduced temperature molar mass, g·mol-1 total number of data points critical pressure, bar normal boiling temperature, K reduced temperature at the normal boiling point critical temperature, K critical volume, cm3·mol-1 fugacity of component i Henry’s law constant of component i in solvent A, Pa 21

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

αlean αrich

the first coefficient of Henry’s law constant the second coefficient of Henry’s law constant the third coefficient of Henry’s law constant the fourth coefficient of Henry’s law constant composition of component i in liquid composition of component i in vapor activity coefficient the infinite dilution activity coefficient of component i total equivalent work, kJ ·kg-1 CO2 the reboiler temperature, K heat duty of reboiler, GJ ·ton-1 CO2 work compressor and blower, kJ ·kg-1 CO2 work of pumps, kJ ·kg-1 CO2 CO2 loading in lean solvent CO2 loading in rich solvent

Abbreviations IL MEA [Bmim][BF4] [Bmim][DCA] [Bpy][BF4] [Bmim][Tf2N] [Bmim][Ac] [Emim][OTf] GLE VLE FC-CS O.F. AIC LVR AARD Nabs Ns TCC TCI TOC VOC FOC

ionic liquid monoethanolamine 1-butyl-3-methylimidazolium tetrafluoroborate 1-butyl-3-methylimidazolium dicyanamide 1-butylpyridinium tetrafluoroborate 1-butyl-3-methylimidazolium bis((trifluoromethyl)sulfonyl)imide 1-butyl-3-methylimidazolium Acetate 1-ethyl-3-methylimidazolium trifluoromethanesulfonate gas-liquid phase equilibrium vapor-liquid phase equilibrium the fragment contribution-corresponding states method objective function absorber inter-cooling lean vapor recompression average absolute relative deviation stage number of absorber stage numbers of stripper total capture cost total capital investment total operating cost variable operating cost fixed operating cost

Greek letters ρ η σ

density, g·cm-3 viscosity, mPa·s surface tension, mN·m-1

aiA, biA ciA diA xi yi

γ

γ i∞

Weq Treb Qreb Wcomp Wpump

Page 22 of 43

22

ACS Paragon Plus Environment

Page 23 of 43

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

λ ω

liquid thermal conductivity, W·m-1·K-1 acentric factor

Superscripts cal exp

calculated property experimental property

REFERENCES 1.

Agency, I. E., World Energy Outlook, 2012. OECD/IEA: 2012.

2.

Kanniche, M.; Gros-Bonnivard, R.; Jaud, P.; Valle-Marcos, J.; Amann, J.-M.; Bouallou, C. Pre-combustion,

post-combustion and oxy-combustion in thermal power plant for CO2 capture. Appl. Therm. Eng. 2010, 30, (1), 53-62. 3.

Artanto, Y.; Jansen, J.; Pearson, P.; Thong, D.; Cottrell, A.; Meuleman, E.; Feron, P. Performance of MEA

and amine-blends in the CSIRO PCC pilot plant at Loy Yang Power in Australia. Fuel 2012, 101, 264-275. 4.

Huang, B.; Xu, S.; Gao, S.; Liu, L.; Tao, J.; Niu, H.; Cai, M.; Cheng, J. Industrial test and techno-economic

analysis of CO2 capture in Huaneng Beijing coal-fired power station. Appl. Energ. 2010, 87, (11), 3347-3354. 5.

Mariz, C. Carbon dioxode recovery: Large scale design trends. J. Can. Pet. Technol. 1998, 37, (7).

6.

Gao, J.; Wang, S.; Zhao, B.; Qi, G.; Chen, C. Pilot-Scale Experimental Study on the CO2 Capture Process

with Existing of SO2: Degradation, Reaction Rate, and Mass Transfer. Energ. Fuel 2011, 25, (12), 5802-5809. 7.

Zhang, X. P.; Zhang, X. C.; Dong, H. F.; Zhao, Z. J.; Zhang, S. J.; Huang, Y. Carbon capture with ionic

liquids: overview and progress. Energ Environ Sci 2012, 5, (5), 6668-6681. 8.

Liu, X.; Zhang, X.; Zhou, G.; Yao, X.; Zhang, S. All-atom and united-atom simulations of guanidinium-based

ionic liquids. Sci. China Chem. 2012, 55, (8), 1573-1579. 9.

Gimeno, M. P.; Mayoral, M. C.; Andrés, J. M. Influence of Temperature on CO2 Absorption Rate and

Capacity in Ionic Liquids. Energ. Fuel 2013, 27, (7), 3928-3935. 10. Zhang, S. J.; Zhang, X. P.; Zhao, Y. S.; Zhao, G. Y.; Yao, X. Q.; Yao, H. W. A novel ionic liquids-based scrubbing process for efficient CO2 capture. Science China-Chemistry 2010, 53, (7), 1549-1553. 11. Zhang, S. J.; Yuan, X. L.; Chen, Y. H.; Zhang, X. P. Solubilities of CO2 in 1-butyl-3-methylimidazolium hexafluorophosphate and 1,1,3,3-tetramethylguanidium lactate at elevated pressures. J. Chem. Eng. Data 2005, 50, (5), 1582-1585. 12. Zhang, S. J.; Chen, Y. H.; Ren, R. X. F.; Zhang, Y. Q.; Zhang, J. M.; Zhan, X. P. Solubility of CO2 in sulfonate ionic liquids at high pressure. J. Chem. Eng. Data 2005, 50, (1), 230-233. 13. Chen, Y. H.; Zhang, S. J.; Yuan, X. L.; Zhang, Y. Q.; Zhang, X. P.; Dai, W. B.; Mori, R. Solubility of CO2 in imidazolium-based tetrafluoroborate ionic liquids. Thermochim. Acta 2006, 441, (1), 42-44. 14. Aki, S. N. V. K.; Mellein, B. R.; Saurer, E. M.; Brennecke, J. F. High-Pressure Phase Behavior of Carbon Dioxide with Imidazolium-Based Ionic Liquids. J. Phys. Chem. B 2004, 108, (52), 20355-20365. 15. Goodrich, B. F.; de la Fuente, J. C.; Gurkan, B. E.; Lopez, Z. K.; Price, E. A.; Huang, Y.; Brennecke, J. F. Effect of Water and Temperature on Absorption of CO2 by Amine-Functionalized Anion-Tethered Ionic Liquids. J. Phys. Chem. B 2011, 115, (29), 9140-9150. 16. Ramdin, M.; de Loos, T. W.; Vlugt, T. J. H. State-of-the-Art of CO2 Capture with Ionic Liquids. Ind. Eng. Chem. Res. 2012, 51, (24), 8149-8177. 17. Brennecke, J. F.; Gurkan, B. E. Ionic Liquids for CO2 Capture and Emission Reduction. J. Phys. Chem. Lett. 2010, 1, (24), 3459-3464. 18. Taib, M. M.; Murugesan, T. Solubilities of CO2 in aqueous solutions of ionic liquids (ILs) and monoethanolamine (MEA) at pressures from 100 to 1600kPa. Chem. Eng. J. 2012, 181–182, (0), 56-62. 23

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 24 of 43

19. Ahmady, A.; Hashim, M. A.; Aroua, M. K. Density, viscosity, physical solubility and diffusivity of CO2 in aqueous MDEA+[bmim][BF4] solutions from 303 to 333K. Chem. Eng. J. 2011, 172, (2–3), 763-770. 20. Ahmady, A.; Hashim, M. A.; Aroua, M. K. Experimental Investigation on the Solubility and Initial Rate of Absorption

of

CO2

in

Aqueous

Mixtures

of

Methyldiethanolamine

with

the

Ionic

Liquid

1-Butyl-3-methylimidazolium Tetrafluoroborate. J. Chem. Eng. Data 2010, 55, (12), 5733-5738. 21. Zhao, Y. S.; Zhang, X. P.; Zeng, S. J.; Zhou, Q.; Dong, H. F.; Tian, X. A.; Zhang, S. J. Density, Viscosity, and Performances of Carbon Dioxide Capture in 16 Absorbents of Amine plus Ionic Liquid + H2O, Ionic Liquid + H2O, and Amine + H2O Systems. J. Chem. Eng. Data 2010, 55, (9), 3513-3519. 22. Camper, D.; Bara, J. E.; Gin, D. L.; Noble, R. D. Room-Temperature Ionic Liquid−Amine Solutions: Tunable Solvents for Efficient and Reversible Capture of CO2. Ind. Eng. Chem. Res. 2008, 47, (21), 8496-8498. 23. Hasib-ur-Rahman, M.; Bouteldja, H.; Fongarland, P.; Siaj, M.; Larachi, F. ç. Corrosion Behavior of Carbon Steel in Alkanolamine/Room-Temperature Ionic Liquid Based CO2 Capture Systems. Ind. Eng. Chem. Res. 2012, 51, (26), 8711-8718. 24. Yang, J.; Yu, X.; Yan, J.; Tu, S.-t. CO2 capture using amine solution mixed with ionic liquid. Ind. Eng. Chem. Res. 2014, 53, 2790-2799. 25. LaFrate, A. L.; Huffman, M. C.; Brown, N.; Shannon, M. S.; Belmore, K.; Bara, J. E.; Brown, A. E. Accelerated aging and qualitative degradation pathway analysis of CO2 capture solvents containing ionic liquids. Energ. Fuel 2012, 26, (8), 5345-5349. 26. Janiczek, P.; Kalb, R. S.; Thonhauser, G.; Gamse, T. Carbon dioxide absorption in a technical-scale-plant utilizing an imidazolium based ionic liquid. Sep. Purif. Technol. 2012, 97, 20-25. 27. Kumar, S.; Cho, J. H.; Moon, I. Ionic liquid-amine blends and CO2BOLs: Prospective solvents for natural gas sweetening and CO2 capture technology—A review. Int. J. Greenh. Gas Con. 2014, 20, (0), 87-116. 28. Shiflett, M. B.; Drew, D. W.; Cantini, R. A.; Yokozeki, A. Carbon Dioxide Capture Using Ionic Liquid 1-Butyl-3-methylimidazolium Acetate. Energ. Fuel 2010, 24, 5781-5789. 29. Basha, O. M.; Keller, M. J.; Luebke, D. R.; Resnik, K. P.; Morsi, B. I. Development of a Conceptual Process for Selective CO2 Capture from Fuel Gas Streams Using [hmim][Tf2N] Ionic Liquid as a Physical Solvent. Energ. Fuel 2013, 27, (7), 3905-3917. 30. Cousins, A.; Wardhaugh, L. T.; Feron, P. H. M. Preliminary analysis of process flow sheet modifications for energy efficient CO2 capture from flue gases using chemical absorption. Chem. Eng. Res. Des. 2011, 89, (8), 1237-1251. 31. Karimi, M.; Hillestad, M.; Svendsen, H. F. Capital costs and energy considerations of different alternative stripper configurations for post combustion CO2 capture. Chem. Eng. Res. Des. 2011, 89, (8), 1229-1236. 32. Cousins, A.; Wardhaugh, L. T.; Feron, P. H. M. A survey of process flow sheet modifications for energy efficient CO2 capture from flue gases using chemical absorption. Int. J. Greenh. Gas Con. 2011, 5, (4), 605-619. 33. Schach, M.-O.; Schneider, R. d.; Schramm, H.; Repke, J.-U. Techno-Economic Analysis of Postcombustion Processes for the Capture of Carbon Dioxide from Power Plant Flue Gas. Ind. Eng. Chem. Res. 2010, 49, (5), 2363-2370. 34. Cadena, C.; Anthony, J. L.; Shah, J. K.; Morrow, T. I.; Brennecke, J. F.; Maginn, E. J. Why Is CO2 So Soluble in Imidazolium-Based Ionic Liquids? J. Am. Chem. Soc. 2004, 126, (16), 5300-5308. 35. Lu, B.; Wang, X.; Xia, Y.; Liu, N.; Li, S.; Li, W. Kinetics of Carbon Dioxide Absorption into Mixed Aqueous Solutions of MEA + [Bmim]BF4 Using a Double Stirred Cell. Energ. Fuel 2013, 27, (10), 6002-6009. 36. Sánchez, L. G. n.; Espel, J. R.; Onink, F.; Meindersma, G. W.; Haan, A. B. d. Density, Viscosity, and Surface Tension of Synthesis Grade Imidazolium, Pyridinium, and Pyrrolidinium Based Room Temperature Ionic Liquids. J. Chem. Eng. Data 2009, 54, (10), 2803-2812. 24

ACS Paragon Plus Environment

Page 25 of 43

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

37. Xu, F.; Gao, H.; Dong, H.; Wang, Z.; Zhang, X.; Ren, B.; Zhang, S. Solubility of CO2 in aqueous mixtures of monoethanolamine and dicyanamide-based ionic liquids. Fluid Phase Equilib. 2014, 365, (0), 80-87. 38. Yunus, N. M.; Mutalib, M. I. A.; Man, Z.; Bustam, M. A.; Murugesan, T. Solubility of CO2 in pyridinium based ionic liquids. Chem. Eng. J. 2012, 189–190, (0), 94-100. 39. Couling, D. J.; Bernot, R. J.; Docherty, K. M.; Dixon, J. K.; Maginn, E. J. Assessing the factors responsible for ionic liquid toxicity to aquatic organisms via quantitative structure–property relationship modeling. Green Chem. 2006, 8, (1), 82-90. 40. Docherty, K. M.; Dixon, J. K.; Kulpa Jr, C. F. Biodegradability of imidazolium and pyridinium ionic liquids by an activated sludge microbial community. Biodegradation 2007, 18, (4), 481-493. 41. Huang, Y.; Dong, H.; Zhang, X.; Li, C.; Zhang, S. A new fragment contribution-corresponding states method for physicochemical properties prediction of ionic liquids. AIChE J. 2013, 59, (4), 1348-1359. 42. Zhao, J.; Jiang, X.-C.; Li, C.-X.; Wang, Z.-H. Vapor pressure measurement for binary and ternary systems containing a phosphoric ionic liquid. Fluid Phase Equilib. 2006, 247, (1–2), 190-198. 43. Shen, C.; Lu, Y.-z.; Li, X.-m.; Li, C.-x. Effect of the Ionic Liquid Triethylmethylammonium Dimethylphosphate on the Vapor Pressure of Water, Methanol, Ethanol, and Their Binary Mixtures. J. Chem. Eng. Data 2011, 56, (5), 1933-1940. 44. Li, G.; Li, Z.; Asselin, E. Determination and Modeling of Vapor–Liquid Equilibria for the Sulfuric Acid + Water + Butyl Acetate + Ethanol System. Ind. Eng. Chem. Res. 2013, 52, (9), 3481-3489. 45. Wang, J. F.; Li, C. X.; Wang, Z. H. Measurement and Prediction of Vapor Pressure of Binary and Ternary Systems Containing 1-Ethyl-3-methylimidazolium Ethyl Sulfate. J. Chem. Eng. Data 2007, 52, (4), 1307-1312. 46. Duan, L.; Zhao, M.; Yang, Y. Integration and optimization study on the coal-fired power plant with CO2 capture using MEA. Energy 2012, 45, (1), 107-116. 47. Oexmann, J.; Kather, A. Minimising the regeneration heat duty of post-combustion CO2 capture by wet chemical absorption: The misguided focus on low heat of absorption solvents. Int. J. Greenh. Gas Con. 2010, 4, (1), 36-43. 48. Pellegrini, L. A.; Moioli, S.; Gamba, S. Energy saving in a CO2 capture plant by MEA scrubbing. Chem. Eng. Res. Des. 2011, 89, (9), 1676-1683. 49. Reid, R. C.; Prausnitz, J. M.; Poling, B. E. The properties of gases and liquids. McGraw Hill,New York, 1987. 50. Olivier-Bourbigou, H.; Magna, L. Ionic liquids: perspectives for organic and catalytic reactions. J. Mol. Catal. A: Chem. 2002, 182–183, (0), 419-437. 51. Kabo, G. J.; Paulechka, Y. U.; Kabo, A. G.; Blokhin, A. V. Experimental determination of enthalpy of 1-butyl-3-methylimidazolium iodide synthesis and prediction of enthalpies of formation for imidazolium ionic liquids. J. Chem. Thermodyn. 2010, 42, (10), 1292-1297. 52. Emel'yanenko, V. N.; Verevkin, S. P.; Heintz, A. The Gaseous Enthalpy of Formation of the Ionic Liquid 1-Butyl-3-methylimidazolium Dicyanamide from Combustion Calorimetry, Vapor Pressure Measurements, and Ab Initio Calculations. J. Am. Chem. Soc. 2007, 129, (13), 3930-3937. 53. Verevkin, S. P.; Emel’yanenko, V. N.; Zaitsau, D. H.; Heintz, A.; Muzny, C. D.; Frenkel, M. Thermochemistry of imidazolium-based ionic liquids: experiment and first-principles calculations. PCCP 2010, 12, (45), 14994-15000. 54. Fredlake, C. P.; Crosthwaite, J. M.; Hert, D. G.; Aki, S. N. V. K.; Brennecke, J. F. Thermophysical Properties of Imidazolium-Based Ionic Liquids. J. Chem. Eng. Data 2004, 49, (4), 954-964. 55. Jacquemin, J.; Gomes, M. F. C.; Husson, P.; Majer, V. Solubility of carbon dioxide, ethane, methane, oxygen, nitrogen, hydrogen, argon, and carbon monoxide in 1-butyl-3-methylimidazolium tetrafluoroborate between temperatures 283 K and 343 K and at pressures close to atmospheric. J. Chem. Thermodyn. 2006, 38, (4), 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

Page 26 of 43

490-502. 56. Kolbe, B.; Kleiber, M.; Rarey, J., Chemical thermodynamics for process simulation. Wiley. com: 2012. 57. Redlich, O.; Kwong, J. On the Thermodynamics of Solutions. V. An Equation of State. Fugacities of Gaseous Solutions. Chem. Rev. 1949, 44, (1), 233-244. 58. Renon, H.; Prausnitz, J. M. Local compositions in thermodynamic excess functions for liquid mixtures. AIChE J. 1968, 14, (1), 135-144. 59. Taghdisian, H.; Farhadi, F.; Pishvaie, M. R. An optimization-oriented green design for methanol plants. J. Chem. Technol. Biotechnol. 2012, 87, (8), 1111-1120. 60. Knudsen, J. N.; Andersen, J.; Jensen, J. N.; Biede, O. Evaluation of process upgrades and novel solvents for the post combustion CO2 capture process in pilot-scale. Energy Procedia 2011, 4, (0), 1558-1565. 61. Husebye, J.; Brunsvold, A. L.; Roussanaly, S.; Zhang, X. Techno Economic Evaluation of Amine based CO2 Capture: Impact of CO2 Concentration and Steam Supply. Energy Procedia 2012, 23, (0), 381-390. 62. Sipöcz, N.; Tobiesen, F. A. Natural gas combined cycle power plants with CO2 capture – Opportunities to reduce cost. Int. J. Greenh. Gas Con. 2012, 7, (0), 98-106. 63. Oyenekan, B. A.; Gary, T. Energy performance of stripper configurations for CO2 capture by aqueous amines. Ind. Eng. Chem. Res. 2006, 45, (8), 2457-2464. 64. Metz, B.; Davidson, O.; de Coninck, H.; Loos, M.; Meyer, L. IPCC special report on carbon dioxide capture and storage; Intergovernmental Panel on Climate Change, Geneva (Switzerland). Working Group III: 2005. 65. Abu-Zahra, M. R. M.; Niederer, J. P. M.; Feron, P. H. M.; Versteeg, G. F. CO2 capture from power plants: Part II. A parametric study of the economical performance based on mono-ethanolamine. Int. J. Greenh. Gas Con. 2007, 1, (2), 135-142. 66. Loh, H.; Lyons, J.; White, C. Process Equipment Cost Estimation, Final Report. National Energy Technoloy Center, DOE/NETL-2002/1169 2002. 67. Sweeting, J., Project Cost Estimating: Principles and Practice. IChemE: 1997. 68. Mores, P.; Rodriguez, N.; Scenna, N.; Mussati, S. CO2 capture in power plants: Minimization of the investment and operating cost of the post-combustion process using MEA aqueous solution. Int. J. Greenh. Gas Con. 2012, 10, 148-163. 69. Peters, L.; Hussain, A.; Follmann, M.; Melin, T.; Hägg, M. B. CO2 removal from natural gas by employing amine absorption and membrane technology—A technical and economical analysis. Chem. Eng. J. 2011, 172, (2–3), 952-960. 70. Desideri, U.; Paolucci, A. Performance modelling of a carbon dioxide removal system for power plants. Energy Convers. Manage. 1999, 40, (18), 1899-1915. 71. Raynal, L.; Bouillon, P.-A.; Gomez, A.; Broutin, P. From MEA to demixing solvents and future steps, a roadmap for lowering the cost of post-combustion carbon capture. Chem. Eng. J. 2011, 171, (3), 742-752. 72. Kim, I.; Svendsen, H. F. Heat of Absorption of Carbon Dioxide (CO2) in Monoethanolamine (MEA) and 2-(Aminoethyl)ethanolamine (AEEA) Solutions. Ind. Eng. Chem. Res. 2007, 46, (17), 5803-5809. 73. Hassan, S. M. N.; Douglas, P. L.; Croiset, E. Techno-Economic Study of CO2 Capture from an Existing Cement Plant Using MEA Scrubbing. Int. J. Green Energy 2007, 4, (2), 197-220.

26

ACS Paragon Plus Environment

Page 27 of 43

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 CAPTIONS Figure 1. The experimental apparatus for saturated vapor pressure measurement. Figure 2. Conventional process flow diagram. Figure 3. Modified process flow diagram. Figure 4. Experimental and calculated heat capacities(a), densities(b), viscosities(c), surface tensions(d) and thermal conductivities(e) of the ILs: dots and lines denote experimental data from literature and calculated value, respectively. Figure 5. Total pressures of the systems of CO2-[Bmim][BF4](a), CO2-[Bmim][DCA](b) and CO2-[Bpy][BF4](c): dots and lines denote experimental data from literature and calculated value. Figure 6. Experimental and predicted vapor pressure data for the binary system [H2O(1)-[Bmim][BF4] (2)] and ternary system [H2O(1)-[Bmim][BF4](2)-MEA(3)]: dots and lines denote experimental data reported here and calculated value, respectively. Figure 7. Experimental and predicted vapor pressure data for the binary system [H2O(1)-[Bmim][DCA] (2)] and ternary system [H2O(1)-[Bmim][DCA](2)-MEA(3)]: dots and lines denote experimental data reported here and calculated value, respectively. Figure 8. Experimental and predicted vapor pressure data for the binary system [H2O(1)-[Bpy][BF4](2)] and ternary system [H2O(1)-[Bpy][BF4](2)-MEA(3)]: dots and lines denote experimental data reported here and calculated value, respectively. Figure 9. Effect of αlean on the regeneration thermal energy requirement of MEA and IL-based processes. Figure 10. Regeneration heat duty and total equivalent energy penalty of three CO2 capture processes. Figure 11. Distribution of TCI of the three CO2 capture processes. Figure 12. Distribution of TCC of the three CO2 capture processes.

27

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 28 of 43

Figure1. The experimental apparatus for saturated vapor pressure measurement. 1-a: working ebulliometer (Type CL-II); 1-b: reference ebulliometer (Type CL-II); 2: magnetic stirrer; 3: heating rod; 4: heating mantle; 5: thermometer; 6:condenser; 7:vacuum jacket; 8:inner casing; 9:buffer bottle; 10- manometer; 11-vacuum pump; 12-evacuation valve; 13-pressure adjusting valve

28

ACS Paragon Plus Environment

Page 29 of 43

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 2. Conventional process flow diagram.

Figure 3. Modified process flow diagram.

29

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 30 of 43

Figure 4(a)

Figure 4(b)

30

ACS Paragon Plus Environment

Page 31 of 43

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(c)

Figure 4(d)

31

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 32 of 43

Figure 4(e)

Figure 4. Experimental and calculated heat capacities(a), densities(b), viscosities(c), surface tensions(d) and thermal conductivities(e) of the ILs: dots and lines denote experimental data from literature and calculated value, respectively.

32

ACS Paragon Plus Environment

Page 33 of 43

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 5(a)

Figure 5(b)

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 43

Figure 5(c)

Figure 5. Total pressures of the systems of CO2-[Bmim][BF4](a), CO2-[Bmim][DCA](b) and CO2-[Bpy][BF4](c): dots and lines denote experimental data from literature11,70,71 and calculated value.

34

ACS Paragon Plus Environment

Page 35 of 43

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 6. Experimental and predicted vapor pressure data for the binary system [H2O(1) - [Bmim][BF4] (2)] and ternary system [H2O(1)-[Bmim][BF4](2)-MEA(3)]: dots and lines denote experimental data reported here and calculated value, respectively.

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

Page 36 of 43

Figure 7. Experimental and predicted vapor pressure data for the binary system [H2O(1) [Bmim][DCA] (2)] and ternary system [H2O(1) -[Bmim][DCA](2)-MEA(3)]: dots and lines denote experimental data reported here and calculated value, respectively.

Figure 8. Experimental and predicted vapor pressure data for the binary system [H2O(1) - [Bpy][BF4] (2)] and ternary system [H2O(1)-[Bpy][BF4](2)-MEA(3)]: dots and lines denote experimental data reported here and calculated value, respectively. 36

ACS Paragon Plus Environment

Page 37 of 43

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 9. Effect of αlean on the regeneration thermal energy requirement of MEA and IL-based processes.

Figure 10. Regeneration heat duty and total equivalent energy penalty of three CO2 capture processes.

37

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 38 of 43

Figure 11. Distribution of TCI of the three CO2 capture processes.

Figure 12. Distribution of TCC of the three CO2 capture processes.

38

ACS Paragon Plus Environment

Page 39 of 43

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. Prediction Equations for Some Temperature Dependent Properties Table 2. Calculation Measures of Total CO2 Capture Cost (TCC) Table 3. The Parameters of Temperature Dependent of Henry Constants Table 4. NRTL Binary Parameters of CO2, H2O and MEA with Different IL Table 5. Key Parameters for the MEA and IL-based Processes

39

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 40 of 43

Table 1. Prediction Equations for Some Temperature Dependent Properties property

specific equation

No.

vapor pressure

C2 i ln Pi s = C1i + T + C3i

(1)

heat capacity

C pi = C1'i + C 2' i T + C3' i T 2

(2)

notation description

ref.

Pi s is vapor pressure in Pa and 46

T is in K, C1i−C3i are equation coefficients. Cpi is liquid molar heat capacity in J·kmol-1·K-1, C ' − C 3' i are

72

1i

equation coefficients. -3 ρ i is density in g·cm , Tci and

density

ρi =

M i Pci RTci  Z

*, RA i

(1 + d i (1 − Tr ) 

viscosity

ln ηi = Ai + Bi / T + Ci ln T

surface tension

σ i = C '' (1 − T / Tci ) ( C

thermal conductivity

1+ (1−T )2/7  r  

'' '' '' 2 '' 3 2 i + C3 i Tri + C4 iTri + C5 iTri )

1i

(3)

Pci are critical temperature and

73

pressure, respectively. Z i*,RA and

(4)

di are equation coefficients. ηi is the liquid viscosity in mPa·s, Ai, Bi and Ci are equation coefficients. σi is the liquid surface tension in

(5)

mN·m-1, C1''i − C 5''i are equation

46

72

coefficients. λi is the liquid thermal conductivity in W·m-1·K-1,

λi = C ''' + C ''' T + C ''' T 2 1i

2i

3i

(6)

C1'''i − C3'''i are equation

72

coefficients.

40

ACS Paragon Plus Environment

Page 41 of 43

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 2. Calculation Measures of Total CO2 Capture Cost (TCC) cost item

calculation basis

TCC

annual capital cost (ACC), total operating cost (TOC)

ACC

total capital investment (TCI)

TCI

fixed capital investment (FCI),working capital, startup cost, initial solvent cost

FCI

direct cost (DC), indirect cost(IC)

DC, IC

purchased equipment cost (PEC) with a coefficient

TOC

variable operating cost (VOC), fixed operating cost (FOC)

41

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 42 of 43

Table 3. The Parameters of Temperature Dependent of Henry Constants parametersa ILs

aiA

biA

ciA

d iA

[Bmim][BF4]

43.962

-10864.735

7.708

-0.119

[Bmim][DCA]

150.2171

-9598.69

-15.9565

-0.0379

[Bpy][BF4]

21.259

7113

-11.1

0.1146

parametersa: i denotes CO2 and A denotes various ionic liquid.

Table 4. NRTL Binary Parameters of CO2, H2O and MEA with Different IL component i CO2 H2O MEA CO2 H2O MEA CO2 H2O MEA

parameters

component j

[Bmim][BF4]

[Bmim][DCA]

[Bpy][BF4]

aij

aji

bij

bji

6.713 -0.9689 -0.265 -3.9246 1.1886 -0.365 81.456 -0.8891 0.165

-16.928 1.114 -0.771 -126.974 -5.514 -20.1705 6.1485 -6.084 -20.471

-2577.397 85 3.495 1178.485 -199 -271.995 -23369.5 248 185.9951

6048.75 732 10.866 32278.18 11510 28964.8 -2021.2 4962.67 9884.8

42

ACS Paragon Plus Environment

Page 43 of 43

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 5. Key Parameters for the MEA and IL-based Processes process

MEA

a

solvent specific demand, m3·ton-1 CO2

21.88

17.07

15.85

CO2 lean loading(αlean)

0.25

0.20

0.20

CO2 rich loading (αrich)

0.469

0.454

0.478

∆ loading of rich and lean solvent

0.219

0.254

0.278

cooled lean solvent temperature, K

334.45

331.15

319.35

heated rich solvent temperature, K

376.85

379.85

366.85

reboiler temperature, K

392.75

399.95

398.85

regeneration heat duty, GJth·ton CO2

3.72

3.17

2.51

electrical power consumption, kWh·ton-1CO2

117.86

115.14

133.09

1047.01

985.89

920.41

143.46

114.53

132.08

-1

-1

total equivalent energy penalty, kJ·kg CO2 -1

cooling water specific requirement, ton·ton CO2

IL

b

M-IL

a

IL: the [Bpy]BF4-MEA capture process M-IL: the modified [Bpy]BF4-MEA capture process

b

43

ACS Paragon Plus Environment