Subscriber access provided by Queen Mary, University of London
Article
Reversible reaction-assisted intensification process for separating the azeotropic mixture of ethanediol and 1,2butanediol: Vapor-liquid equilibrium and Economic evaluation Hong Li, Zhenyu Zhao, Jie Qin, Rui Wang, Xingang Li, and Xin Gao Ind. Eng. Chem. Res., Just Accepted Manuscript • DOI: 10.1021/acs.iecr.7b04921 • Publication Date (Web): 21 Mar 2018 Downloaded from http://pubs.acs.org on March 22, 2018
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.
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 32 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
1
Reversible reaction-assisted intensification process for
2
separating the azeotropic mixture of ethanediol and
3
1,2-butanediol: Vapor-liquid equilibrium and Economic
4
evaluation
5
Hong Li, Zhenyu Zhao, Jie Qin, Rui Wang, Xingang Li, Xin Gao∗
6
School of Chemical Engineering and Technology, National Engineering Research Center of Distillation
7
Technology, Collaborative Innovation Center of Chemical Science and Engineering (Tianjin), Tianjin University,
8
Tianjin 300072, China
9
Abstract: Vapor-liquid equilibrium data of six binary systems composed of ethylene
10
glycol
11
2,3-butanediol (2,3-BD) were measured and regressed by NRTL model. With the
12
binary interaction parameters regressed by experimental data, four separation
13
processes to obtain polymer grade EG from the synthesis products were simulated,
14
including conventional distillation route, azeotropic distillation route, distillation
15
coupled with liquid-liquid extraction route and reaction-assisted distillation route. The
16
total annual costs and profits of the four routes were estimated and compared, which
17
indicates that the annual profit was mainly determined by the yield of EG. The results
18
also demonstrate the reaction-assisted distillation route was proved to be the best
19
strategy currently among those separation methods.
20
Keywords: Ethylene glycol, Vapor-liquid equilibrium, Process simulation, Economic
21
evaluation, Reactive separation ∗
(EG),
1,2-propylene
glycol
(1,2-PG),
1,2-butanediol
Corresponding author: Xin Gao E-mail:
[email protected]. Tel: +86-022-27404701. 1
ACS Paragon Plus Environment
(1,2-BD)
and
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
1. Introduction
2
Ethylene glycol (EG) is widely used in petrochemical industry and other fields, such
3
as polymer ester manufacturing and anti-freezer. It is mainly produced by ethylene
4
epoxidation and ethylene epoxide hydrolysis. However, the worldwide oil shortage
5
and the increasing demand of EG call for alternative non-oil routes to produce
6
glycols.1 Among those routes, C1 chemical process based on syngas mainly from coal
7
is proved to be promising.2 The liquid product of indirect syngas-to-EG reaction
8
contains 48.0wt%-51.0wt% methanol, 39.0wt%-43.0wt% EG, 1.5wt-2.5wt% water,
9
2.6wt-3.5wt% ethanol (EtOH), 2.0wt%-2.5wt% glycols including 1,2-butanediol
10
(1,2-BD), 2,3-butanediol (2,3-BD), 1,2-propylene glycol (1,2-PG), 1.5wt%-2.1wt%
11
unreacted methyl nitrite (MN), dimethyl oxatate (DMO), methyl glycolate (MG),
12
dimethyl carbonate (DMC) and 0.3wt-0.5wt% heavy glycols, including diethylene
13
glycol (DEG), triethylene glycol (TEG) and others.3
14
Currently, the traditional mode of separating the mixture consumes excessive amount
15
of energy because of the application of distillation column with large reflux ratio and
16
tray number to remove other glycols from EG, resulting in poor competitive
17
advantage of coal-based glycol production. Therefore, several methods have been
18
proposed to enhance the separation process, for example, azeotropic distillation.
19
Much work4,
20
reports on industrial application have rarely been seen. Besides, extraction is expected
21
to be used in separation process but the recovery and consumption of extracting
22
solvent should be taken into account.6 Moreover, reaction-assisted separation for the
5
has been devoted to selecting proper entrainer for the system but
2
ACS Paragon Plus Environment
Page 2 of 32
Page 3 of 32 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
1
glycols developed in recent years was proved to be of great potential.7 In addition,
2
other novel separation ideas have been reported, such as extractive distillation2, 8 and
3
adsorption9, but they didn’t quite come off due to high toxicity or large wasting and
4
were dismissed in this paper.
5
Those ideas are all great but it is still not clear which route is more applicable to the
6
glycol system. Therefore, in this paper, four separation processes were designed and
7
evaluated based on process simulation by Aspen Plus. To ensure the accuracy of
8
process simulation, the vapor-liquid equilibrium data of glycol systems should be
9
firstly measured to retrieve the parameters of thermodynamic model.
10
2. Vapor-liquid equilibrium data
11
2.1 Measurement of thermodynamic data
12
The isobaric vapor-liquid equilibrium (VLE) data of the glycols has already been
13
measured by many researchers but it is still not systematic, especially the data of
14
binary systems composed of 1,2-BD, 2,3-BD, 1,2-PG and EG.10-15 Therefore, the
15
double circulating vapor-liquid equilibrium still was used to take measurements. The
16
temperature was measured with a precise digital thermometer (1552A-12-DL model,
17
provided by Fluke Calibration) with the accuracy ±0.01K. The pressure was regulated
18
by a pressure control system, including one vacuum pump, one digital manometer
19
(AOB-20 model) with the accuracy ±0.1 kPa and an aneroid barometer (DYM3 model)
20
with the accuracy ±0.2 kPa.
21
The VLE measurements for the systems (EG +1,2-BD), (EG +1,2-PG), (EG +2,3-BD),
22
(1,2-BD +1,2-PG), (1,2-BD +2,3-BD), (1.2-PG +2,3-BD) at pressure of 10.0kPa, 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 32
1
30.0kPa and 101.3 kPa were performed in the still with 110 cm3 capacity. Liquid
2
phase and condensed vapor phase were continuously circulated to ensure rapid
3
establishment of vapor-liquid equilibrium. The system was thoroughly cleaned by the
4
standard sample before experiment to eliminate contaminants. During the experiments,
5
the equilibrium was assumed to be reached when the pressure and temperature
6
maintained for 45 minutes.
7
The composition of the sample was analyzed by gas chromatograph (GC7890Ⅱ,
8
Techcomp (China) Ltd. ) equipped with a hydrogen flame ionization detector (FID)
9
and a DB-WAX capillary column (30m × 0.25 mm × 0.25 µm, Agilent). The
10
temperature of injector, oven and detector were kept at 523.15K, 393.15K and
11
543.15K, respectively. The carrier gas was high purity nitrogen (mass fraction 0.9999)
12
and the flow rate was 1 mL•min−1. Each VLE data was measured at least three times
13
to ensure accuracy.
14
results by using the external standard curve to convert the peak areas into mole
15
fractions of each sample. Taken the measurement error into consideration, the
16
uncertainties of the vapor and liquid mole fractions were both calculated to be within
17
0.001. The isobaric VLE data (T, xi, yi) for the six binary mixtures were measured and
18
shown in Table S1-S6 in Supporting Information. The activity coefficients γi of the
19
liquid phase could be calculated from the experiment data by Eq.(1), as the vapor
20
phase was assumed as ideal gas under low experimental pressures.
The external standard was chosen to obtain quantitative analysis
(1) 21
where
was the saturated vapor pressure of pure substance i at the system 4
ACS Paragon Plus Environment
Page 5 of 32 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
1
temperature T, which was calculated by Extended Antoine equation based on Aspen
2
Plus database. yi is the vapor phase mole fraction, xi is the liquid phase mole fraction,
3
and P is the total pressure of the system. The VLE data of EG +1,2-PG and 1,2-PG
4
+1,2-BD at 10.0kPa and EG +1,2-BD at 101.3kPa has already been measured by other
5
researchers10,13 and was not measured repeatedly.
6
2.2 Consistency test and data regression
7
To verify the measured quality of VLE values, all experimental results were checked
8
with Herington test and van Ness test proposed by Kang et.al16. As Herington17
9
suggested, if
, the pertinent experimental isobaric VLE data can be
10
considered thermodynamically consistent. In the process, D and J were calculated as
11
Eq.(2) and Eq.(3).
(2)
(3)
12
where Tmax and Tmin were the highest and lowest temperature in the system,
13
respectively. The γ values of the systems (2,3-BD +1,2-PG, 2,3-BD +1,2-BD, EG
14
+1,2-PG, 1,2-PG +1,2-BD) were all between 0.95 and 1.10 at 10.0kPa, 30.0kPa and
15
101.3kPa, hence the area test was not employed. The values of D -Jfor the system
16
{EG +1,2-BD} were calculated to be 1.48, 0.09, 1.05 at 10.0kPa, 30.0kPa, 101.3kPa,
17
respectively, while the values of the binary system {2,3-BD +EG} were 3.71, 2.27
18
and 1.88 at 10.0kPa, 30.0kPa and 101.3kPa, which meant the experimental data
19
passed the Herington test. 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 32
1
Van Ness point test described by Fredenslund18 was applied to confirm the accuracy
2
of the experimental data. The index average absolute deviation (AAD) was defined by
3
following equation: (4)
4
where Ci is the independent variable (bubble point temperature, vapor phase mole
5
fraction) of the pure component i, N is the number of experimental data points, the
6
superscript exp denotes data from experiments and cal denotes data calculated with
7
thermodynamic models.
8
The experimental data were correlated by minimizing object function Q, which was
9
defined by the sum of squared differences between calculated and experimental data
10
as Eq.(5): N
Q = ∑[( i =1
Te ,i - Tm ,i σ T ,i
)2 + (
Pe,i - Pm ,i σ P ,i
)2 + (
xe ,i - xm ,i σ x ,i
)2 + (
y e ,i - y m ,i σ y ,i
)2 ]
(5)
11
where y, x represent mole fraction of pure component in the vapor phase and in the
12
liquid phase. T and P are temperature and pressure. i represents each pure component
13
in the binary mixtures, N is the number of experimental data points, and σ represents
14
the estimated standard deviation of measured experimental values. e, m represent
15
estimated value and measured value, respectively.
16
Correlated interaction parameters, together with the average absolute deviation (AAD)
17
of equilibrium temperature, vapor phase mole fraction for the binary systems between
18
experimental values and calculated values (obtained from Wilson, NRTL and
19
UNIQUAC models) were shown in Table 1.
20
Table 1 Correlated Interaction Parameters and the Average Absolute Deviations (AAD) between 6
ACS Paragon Plus Environment
Page 7 of 32 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 and Calculated Values
1
AAD (T)/K
AAD kPa
11.2114 1025.1681 -481.7063
1,2-PG(1) +EG(2) -890.5198 0.0059 -164.0382 0.3 0.0060 102.8876 0.0061
0.0285 0.0402 0.0404
6.7907e-06 9.3288e-06 9.3784e-06
30.0
Wilson NRTL UNIQUAC
-303.8596 379.5311 -149.2076
-435.0841 407.0284 -160.0364
0.3
0.0086 0.0087 0.0087
0.0290 0.0367 0.0363
2.2166e-05 2.9864e-05 2.9381e-05
101.3
Wilson NRTL UNIQUAC
-303.8596 306.2324 -59.9569
-502.8533 598.3744 -301.4157
0.3
0.0068 0.0081 0.0080
0.0898 0.0942 0.0932
2.9896e-04 3.0985-04 3.0700e-04
10.0
Wilson NRTL UNIQUAC
-365.3776 -274.0664 76.2143
1,2-BD(1) +EG(2) 177.9303 0.0074 471.9168 0.3 0.0074 -137.4402 0.0074
0.0055 0.0094 0.0057
1.3376e-06 2.2593e-06 1.3666e-06
30.0
Wilson NRTL UNIQUAC
-553.1914 -344.9806 137.4240
218.5276 688.6855 -253.4244
0.3
0.0039 0.0046 0.0040
0.0137 0.0225 0.0132
1.1115e-05 1.8333e-05 1.0805e-05
101.3
Wilson NRTL UNIQUAC
-357.4847 -234.7664 12.3188
129.9490 476.8356 -87.5923
0.3
0.0030 0.0032 0.0030
0.0099 0.0110 0.0098
3.3322e-05 3.6732e-05 3.2597e-05
10.0
Wilson NRTL UNIQUAC
-466.6218 -343.8058 154.1807
2,3-BD(1) +EG(2) 251.7104 0.0085 550.6957 0.3 0.0086 -233.8679 0.0084
0.0468 0.0313 0.0450
1.0411e-05 6.8939e-06 1.0131e-05
30.0
Wilson NRTL UNIQUAC
-345.4576 -293.5719 72.6219
182.9973 477.0140 -126.3332
0.3
0.0039 0.0038 0.0039
0.0152 0.0206 0.0149
1.2047e-05 1.6138e-05 1.2037e-05
101.3
Wilson NRTL UNIQUAC
-174.4196 -94.9450 -91.3732
44.6094 227.5011 31.7564
0.3
0.0082 0.0082 0.0082
0.0118 0.0094 0.0128
4.0027e-05 3.3662e-05 4.2203e-05
10.0
Wilson NRTL UNIQUAC
214.5182 385.0055 -192.5574
1,2-PG(1) +1,2-BD(2) -309.1564 0.0039 -290.5894 0.3 0.0040 150.5187 0.0039
0.0494 0.0510 0.0492
1.1516e-05 1.1953e-05 1.1472e-05
30.0
Wilson NRTL UNIQUAC
-104.6842 -139.0940 98.9830
82.7958 167.2831 -119.0317
0.3
0.0066 0.0066 0.0066
0.0197 0.0193 0.0186
1.6738e-05 1.6403e-05 1.5832e-05
101.3
Wilson NRTL UNIQUAC
242.3635 408.5374 -225.1837
-369.2433 -300.8118 171.5125
0.3
0.0047 0.0050 0.0047
0.0122 0.0130 0.0122
3.8772e-05 4.0436e-05 3.8509e-05
model
B12/K
B21/K
10.0
Wilson NRTL UNIQUAC
p/kPa
α
AAD (y)
7
ACS Paragon Plus Environment
(P)/
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 32
10.0
Wilson NRTL UNIQUAC
-304.6439 -296.5335 149.0168
2,3-BD(1) +1,2-PG(2) 215.6571 0.0043 389.0352 0.3 0.0042 -188.3749 0.0043
30.0
Wilson NRTL UNIQUAC
-346.6179 -316.4019 165.1422
240.3541 416.5423 -210.8901
101.3
Wilson NRTL UNIQUAC
219.9860 333.0377 -29.0645
-299.5595 -268.0018 24.5504
0.0093 0.0092 0.0087
0.0030 0.0042 0.0226
1.1982e-05 1.716e-05 7.0655e-05
Wilson NRTL UNIQUAC
277.1605 418.4868 -259.2443
2,3-BD(1) +1,2-BD(2) -454.6703 0.0054 -310.2117 0.3 0.0061 192.6849 0.0056
0.0931 0.0792 0.0864
2.0524e-05 1.7419e-05 1.9031e-05
Wilson
361.2147
-817.1036
0.0053
0.0081
6.3589e-06
NRTL
799.0103
-464.4162
0.0056
0.0427
3.3472e-05
UNIQUAC
-457.6764
282.4736
0.0052
0.0251
1.9698e-05
Wilson NRTL UNIQUAC
-412.5052 -366.2347 198.3000
290.0799 469.8231 -249.2973
0.0041 0.0042 0.0041
0.0547 0.0548 0.0561
1.7443e-04 1.74244-04 1.7866e-04
10.0
30.0
101.3
0.0210 0.0246 0.0217
4.692e-06 5.4823e-06 4.8438e-06
0.3
0.0058 0.0056 0.0058
0.0032 0.0053 0.0032
3.1805e-06 4.3514e-06 3.1468e-06
0.3
0.3
0.3
1
The AAD of bubble point temperature, pressure and vapor phase mole fraction are
2
less than 0.1K, 1Pa and 0.01 respectively in all of cases, which indicate that the three
3
models agree well with experimental data and can be used in industrial application.
4
Among those models, non-random two liquid (NRTL) model, as a common tool for
5
fitting in many literatures19-21, has been reported for VLE and LLE in multicomponent
6
systems consisting alcohols, water, aromatic hydrocarbons, e.g. toluene and glycols
7
that are all involved in the simulation while the other two models have rarely been
8
reported. Therefore, NRTL thermodynamic model was applied to simulate the
9
separation process. The activity coefficient calculated in Aspen Plus is shown as
10
Eq.(6). (6)
8
ACS Paragon Plus Environment
Page 9 of 32 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
1
where
2 3 4
In order to retireve the binary interaction parameters correspongding the models, the
5
data regression module of Aspen Plus was employed based on the least-squares
6
method subject to the maximum likehood principle following Britt-Luecke
7
algorithm22.The model parameters regressed by experimental data were shown in
8
Table 2. Table 2 Model parameters fitted by experimental data
9 component i
1,2-PG
1,2-BD
2,3-BD
1,2-PG
2,3-BD
2,3-BD
component j
EG
EG
EG
1,2-BD
1,2-PG
1,2-BD
0 0 306.2324 598.3744 0.3
0 0 -234.7664 476.8356 0.3
0 0 -94.9450 227.5011 0.3
0 0 408.5374 -300.8118 0.3
0 0 333.0377 -268.0018 0.3
0 0 -366.2347 469.8231 0.3
(K) (K)
10
The NRTL parameters of the binary systems composed of methanol, ethanol, glycols,
11
DMC, DMO and MG were adjusted based on the reported literatures1, 23-25.
12
According to VLE results correlated by NRTL model, the key to obtain EG with high
13
quality is removal of other glycols, including 1,2-PG, 2,3-BD as these components
14
form close boiling point mixtures at the entire range of their compositions, especially
15
1,2-BD which forms azeotrope with EG shown in Figure 1. The difficulty to separate
16
the mixture of glycols only by conventional distillation requires application of other
17
special distillation technology.
18
3. Processes descriptions and simulations
19
To obtain polymer grade EG, three separation processes were proposed and simulated 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
1
by Aspen Plus to calculate equipment sizes and energy cost, including azeotropic
2
distillation route, distillation coupled with liquid-liquid extraction and distillation
3
coupled with acetalization reaction. Here separation process currently applied in
4
conventional industry was also simulated as an investment for comparison. The
5
distillaton columns were firstly designed by DSTWU modules to determine tray
6
number, reflux ratio and distillate rate. The columns were then rated by RadFrac
7
modules. The operational conditions of the equipments were optimized by sensitive
8
analysis with the object to decrease energy cost, e.g. reduce reflux ratio and reboiler
9
duty under the prerequisite of product quality.
10
3.1 Conventional distillation route (CD route)
11
The four-column distillation process widely used in coal-based ethylene glycol
12
industry was shown in Figure 2. The operational conditions of columns and main
13
information of streams were also demonstrated in the figure while some details were
14
omitted due to space reasons.
15
Note that methanol could be served as raw materials to produce MN and DMO in
16
previous CO coupling reaction.3 Therefore, the majority of methanol was recovered
17
by Col-11 with 57 stages, where the reflux ratio was set at 1.9 to ensure the methanol
18
purity according to simulation results. The recovery of methanol was determined
19
according to sensitive analysis of Col-11. From Figure 3, when the recovery of
20
methanol is too high a sharp increase of reboiler duty occurs due to increasing
21
difficulty of separation according to the tendency of minimum reflux ratio Rmin. Thus
22
the methanol recovery was reasonably set at 95.0%. B11 from the bottom stream of 10
ACS Paragon Plus Environment
Page 10 of 32
Page 11 of 32 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
1
Col-11 was separated by Col-12 to remove the remaining methanol, by-product
2
ethanol and other light components. D12 was assumed to be used as fuel and was
3
unnecessary to be purified. Then unreacted MG and DMO was recovered by Col-13.
4
For separation of glycols, it is widely used as expediency in industry to distillate a
5
part of EG carrying almost all other glycols to the top of column, of which the
6
concentration of EG should be up to 95wt% so that it could be still sold as qualified
7
product at relative lower price. Then EG (> 99.9wt%) could be gained from the
8
sidestream while other heavy glycols (DEG, TEG) were removed from the system at
9
the bottom of Col-14. Obviously, it is necessary to apply large tray number and large
10
reflux ratio to Col-14 to ensure the purity of EG at the sidestream. The concentration
11
profile of EG in Col-14 (at RR=24) was shown in Figure 4. As it will keep up to
12
99.9wt% below the 15th plate in vapor phase, the drew-position was accordingly set at
13
the 16th plate considering tower performance. Col-14 was operated at 0.05bar to
14
decrease reboiler temperature so that cheaper low pressure steam could be used.
15
3.2 Azeotropic distillation route(AD route)
16
The form of close boiling point mixture and even azeotropic mixture discounts EG
17
quality and yield in CD route and hence causes large diseconomy in manufacture.
18
Therefore, many researches have been devoted to development of azeotropic
19
distillation technology for glycols system.
20
The application of azeotropic distillation to separating glycols system was firstly
21
proposed by Berg26 in 1990. He listed about 30 potential entrainers to enhance the
22
separation via heterogeneous azeotropic distillation. The past two decades has 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
1
witnessed developing performance of azeotropic distillation technology in separation
2
of EG and 1,2-BD. Song5 has used homotopy-Newton united method to predict
3
azeotropic component and proved ethylene benzene was the best entrainer by fuzz
4
mathematics. Wang27 simulated separation process for EG and 1,2-BD through
5
azeotropic distillation and compared it with conventional distillation. The result
6
showed that azeotropic distillation is superior to conventional distillation on the
7
operating cost and investment. Li28 conducted several experiments to separate EG and
8
1,2-BD and the purity of EG can be up to 99.7wt% and recovery reached 80.32%.
9
Furthermore, the feasibility of application of azeotropic distillation to separating
10
glycols has been studied by experiments.4 Based on those researches, a feasible
11
azeotropic distillation process with the optimized results was shown in Figure 5.
12
Similar to CD route, methanol and DMO should be recovered firstly and meanwhile
13
light impurities, including ethanol and water were removed from the system by
14
Col-21, Col-22 and Col-23 which were similar with Col-11, Col-12 and Col-13 and
15
not shown repeatedly in Figure 5. Then the mixture of glycols was pumped into
16
Col-24, where the mixture of EG, 1,2-BD, 1,2-PG, 2,3-BD was distillated out. The
17
majority of EG mixed with heavy glycols from the bottom was separated by Col-25.
18
1,2-BD in the mixture from the top of Col-24 was expected be removed by azeotropic
19
distillation using ethylbenzene (EB) as entrainer. The model parameters of binary
20
mixture of EG and ethylbenzene were adjusted according to Yang’s work29 to ensure
21
the accuracy of simulation. According to the result of property analysis in Aspen Plus,
22
1,2-PG is bad for the process of azeotropic distillation. Therefore, the mixture should 12
ACS Paragon Plus Environment
Page 12 of 32
Page 13 of 32 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
1
be pretreated to remove a part of 1,2-PG from the system by Col-26 with 69 trays to
2
ensure quality of the final product. The bottom product of Col-26 was pumped into
3
Col-27, together with EB, where 1,2-BD was carried by EB to the top and thus EG
4
(>99.9wt%) could be obtained at the bottom. Then ethylbenzene-rich phase from the
5
decanter flowed back to Col-27 while glycol-rich phase was treated by Col-28 to
6
recover EB.
7
The mass ratio of entrainer to feed was 0.5-0.8 according to the result of sensive
8
analysis shown in Figure 6. Under the goal that the recovery of 1,2-BD in Col-27 was
9
set at 0.995, the addition of entrainer makes it easier to separate the two glycols,
10
therefore the minimum tray number and input energy decrease with the increase of EB.
11
However, too much entrainer will cause the increase of column capacity and thus the
12
increase of energy cost.
13
3.3 Distillation coupled with liquid-liquid extraction route(ED route)
14
Apart from azeotropic distillation, there are some researches on application of
15
liquid-liquid extraction on separating EG and 1,2-BD. It has been proved that toluene
16
was the best agent among aromatics, ethers, alkanes as well as ketones.6, 9, 28 Sun23
17
simulated several extraction processes for separating EG and 1,2-BD using toluene as
18
extracting solvent. Based on those researches, a detailed conceptual design and
19
operation of separation process was shown in Figure 7. Col-31-Col-33 were similar
20
with Col-11-Col-13 in Figure 2 and not shown.
21
After the recovery of methanol, MG and DMO, the stream B33 was pumped to the
22
extraction section. The extraction column was operated at 313K and 1 bar based on 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 32
1
the experimental results in the literature28, in which the feed is counter-currently
2
brought into contact with toluene i.e. extracting solvent. The extract stream leaving
3
the extraction column contains 1,2-BD, coextracted EG and extracting solution.
4
The process designs included extraction and distillation, which were simulated by
5
using the Extract and RadFrac modules of Aspen Plus, respectively while the
6
equilibrium stage model was used, of which the distribution coefficient of toluene was
7
estimated by UNIFAC group-contribution method with the help of LLE data20. To
8
ensure the quality of EG product, the mass fraction of 1,2-BD in raffinate phase
9
should be below 0.001, then the minimum solvent to feed mole ratio was calculated to
10
be 0.1 using Eq. (7). (7)
11
where S is flow rate of the solvent, F is flow rate of the feed, Xin is mass fraction of
12
the solute i.e. 1,2-BD in the feed, Xout is mass fraction of the solute in the raffinate
13
phase, Yin is the mass fraction of the solute in the entering solvent, and Yout is the mass
14
fraction of the solute in the extractant phase, which is in equilibrium with Xin. As
15
suggested by King30, an appropriate estimation of the optimal solvent to feed ratio
16
should be 1.15 ~2 times of (S/F)min to ensure product quality and in this work it was
17
set to be 1.5. Therefore, the flow rate of solvent was calcultated to be about 5000kg/h.
18
The extraction phase (D34) including the majority of toluene and 1,2-BD flowing out
19
from the top of extraction tower was separated by Col-36 to recover the solvent while
20
the raffinate phase (B34) from the bottom was separated by Col-35, where the mixture
21
of 1,2-PG and toluene was distillated from the system and treated as fuel. The bottom 14
ACS Paragon Plus Environment
Page 15 of 32 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
1
liquid of Col-35 was separated by Col-38, where EG (>99.9wt%) was obtained at the
2
top. Col-36-Col-38 were operated under lower pressure to reduce operational
3
temperatures to keep products/solvent from spoilling and save energy.
4
3.4 Reaction-assisted distillation route (RD route)
5
The novel idea has been developed in recent years that chemical reactions could be
6
used to enhance separation process. Dhale’s work31 has shown excellent performance
7
in the recovery of EG, 1,2-PG and 1,2-BD from aqueous solution through
8
acetalization reaction. The acetals formed by acetaldehyde (ALK) and glycols were
9
easier to be separated and then additional reactive distillation columns were needed
10
for hydrolysis of acetals back to EG and 1,2-PG, which showed great advantage for
11
EG with low concentration. As for manufacturing of EG, Huang et al.32 studied
12
reaction kinetics and developed a reactive distillation process for separation of EG
13
and 1,2-BD. Based on those previous researches, a conceptual design was established
14
shown in Figure 8.
15
Similar to AD route in Figure 5, light components were removed by Col-41-Col-43
16
and the majority of EG was obtained after distillation by Col-44 and Col-45. Here
17
they were omitted as they were similar with AD route. The brief introduction of
18
Reaction-assisted distillation (RD) route begins with the stream D44 from Col-44.
19
Unlike AD route, the distillates from Col-44 (similar with Col-24) was pumped,
20
together with acetaldehyde into the reactor R41, which was operated under 1.5 bar,
21
333.15 K. A 2:1 mole ratio of acetaldehyde to glycols was set to ensure the conversion
22
of glycols based on the previous study32. Besides, the conversions of EG and 1,2-BD 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 32
1
were fixed to 0.61 and 0.70, respectively, based on the previously published
2
experimental results7. The conversions of 1,2-PG and 2,3-BD have not been reported
3
and they were assumed to be 0.70, due to their chemical properties were similar with
4
1,2-BD and their concentrations were too little to have effect on energy cost. The
5
acetalization reactions of glycols following Eq(8)-(11) occurred in the reactor. OH H
O
-H2O
+ O
OH
(ALK)
(EG)
(8) O
(2-MD)
OH H
O
-H2O
+
(9) O
OH
O
(ALK)
(1,2-PG)
(2,4-DMD)
OH H
O
(10)
-H2O
+ O
OH
O
(ALK)
(1,2-BD)
(4-EMD)
OH H
O
-H2O
+ OH
(2,3-BD)
O
(11) O
(2,4,5-TMD)
(ALK)
6
The distillation tower Col-46 was to recover the unreacted acetaldehyde and glylcols.
7
The distillation tower Col-47 and Col-48 were aimed to separate the acetals.Then EG
8
could be obtained by acetals hydrolysis in the continuous reactive distillation column
9
R42, where the conversions of hydrolysis reaction were set at 100% according to the
10
experimental results in a pilot scale column as mentioned by Atul D. Dhale33. The
11
hydrolysis of 2,4-DMD, 4-EMD and 2,4,5-TMD occurred in R43, where acetaldehyde
12
was recovered. Then water in S42 and S43 was removed by Col-48 and Col-49,
13
respectively. 16
ACS Paragon Plus Environment
Page 17 of 32 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
1
The process included reactor and distillation tower, which were simulated using the
2
Rstoic and RadFrac modules of Aspen Plus, respectively, while the reactive
3
distillation columns R42-R44 were simulated by using RadFrac modules. The binary
4
interaction parameters of acetals and water were obtained from regression of the VLE
5
data measured by Chopade et al.34,
6
equilibrium (LLE) data of the binary system 2MD +4EMD were measured in the
7
preliminary work7.
8
4. Economic estimation
9
After the initial steady state design of each section, optimization was performed to
10
find a set of optimal design and operating variables with minimizing total annual cost
11
(TAC), including adjusting tray number and reflux ratio, choosing the optimal
12
position to feed and draw-out36. According to Luyben’s work37, the four processes
13
were evaluated by annual cost and profit. Total annual cost (TAC) contains capital
14
investigation and operational cost.
35
. The VLE data as well as liquid-liquid
(12) 15
Payback period was usually set at 3 years in industry. The capital investigation
16
included columns and heat exchangers and the purchase costs of pumps and other
17
equipments were too small to be taken in account38.
18
Stainless steel was specified as the material in coolers, reboilers, tanks, and internals
19
to avoid iron contamination in the final product. Therefore, the cost of columns was
20
calculated by Eq.(13) based on the height of column L and the diameter of column D
21
esitimated by Tray Sizing module in Aspen Plus following Eq.(14). 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 32
(13) (14) 1
where NT is tray number of each column.
2
The estimation of reactor cost followed Eq. (13), either. Under the assumption that the
3
reactor could be considered as CSTR(continuous stirred tank reactor) and the ratio of
4
D and L was generally set at 3. D and L were estimated by reactor volume, which was
5
calculated by Eq. (15). (15)
6
where Q is the flow rate of input stream, t is space time depending on optimal time of
7
reaction at different operating conditions. α is the ratio of effective volume to actual
8
volume, which is always set at 1.5~2 in industrial design to take the cost of additional
9
volumes and internals into account. The cost of catalyst is too low to be considered.
10
The cost of heat exchangers was calculated following Eq.16 based on the heat transfer
11
area of each heat exchanger. The latter was estimated by Aspen Exchanger Design and
12
Rating following Eq.(16). (16)
13
Energy cost was calcuated based on the heat duty of heat exchangers. For reboilers
14
and heaters, the price of high pressure (41barg, 251℃) steam was 9.87$/GJ while the
15
prices of middle pressure (10barg, 184℃) and low pressure (5barg, 160℃) steam
16
were 8.22$/GJ and 7.78GJ$, respectively. In the four separation processes, the
17
condensers were all cooled by cooling water, whose price was 0.354$/GJ.
18
Apart from input energy, those separations need additional materials for supplement, 18
ACS Paragon Plus Environment
Page 19 of 32 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
1
including ethylbenzene in AD route, toluene in ED route as well as acetaldehyde in
2
RD route. The output value was determined by production of EG. The prices of
3
chemicals were inquired in Molbase database. The ecomomic estimation results of the
4
four separation processses were shown in Table 3. Table 3 Economic estimation of separation processes
5 Price
6
Column&Reactor/10 $ heat exchanger/106$ energy cost/106$/a additional material//106$ TAC/106$ output/(106$/a) profit/(106$/a)
CD
AD
ED
RD
9.72 4.36 44.00 0.00 48.69 296.28 247.59
8.99 6.25 17.94 0.40 23.42 317.75 294.33
7.88 5.63 25.10 3.74 33.34 318.42 285.08
9.19 6.08 23.76 0.01 28.86 336.76 307.89
6
Table 3 reveals that TAC is mainly determined by energy cost and the annual profit
7
depends mainly on the yield of EG (output). In CD route, the existence of 1,2-PG,
8
1,2-BD and 2,3-BD will result in decrease of EG output and increase of imported
9
energy and hence cause reduction of annual profit. Compared with the CD route, the
10
application of azeotropic distillation technology improves the yield of EG and
11
reduced the energy cost of separation, which demonstrated great potential. However,
12
ethylbenzene, as is proved to be the best entrainer though, has danger of serious
13
damage to health by prolonged exposure through inhalation and if swallowed.
14
Therefore, the application of azeotropic distillation technology on separation of
15
glycols system in industry requires development of non-toxic entrainer, e.g.
16
ionic-liquid. The application of liquid-liquid extraction on separating EG and 1,2-BD
17
is able to increase the productivity of EG but the process to remove 1,2-PG and
18
recover extracting agent still consumes large quantity of energy. In addition, owing to
19
great quantity and high toxicity of extracting solvent, ED route shows no advantage 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
1
whether from economic angle or from the perspective of the environment. By
2
comparison, reaction-assisted separation with acetaldehyde is an effectively way to
3
solve the puzzle of separating glycols mixture due to low energy consumption and
4
large EG output while additional products such as 1,2-PG and 1,2-BD could be
5
obtained. Thus, the RD route was calculated to be the best strategy for separation of
6
liquid product of ethylene glycol synthesis from syngas. Furthermore, as an idea
7
developed in recent years, reverse reacion-assisted separation needs more researches,
8
such as, developing new suitable reverse reactions e.g. reaction with butanone and
9
propionaldehyde39, application of reactive distillation to promote reaction conversion
10
and save energy, etc.
11
5. Conclusion
12
In this work, the rigorous steady state design and optimization work of four separation
13
processes for liquid products of EG synthesis from syngas were investigated. The
14
thermodynamic data of glycol system was measured to ensure the accuracy of process
15
simulation by Aspen Plus.
16
The result of economic analysis indicated that azeotropic distillation owned great
17
potential due to lower energy cost compared to CD route but it still needed promotion,
18
for example, development of entrainer. Among those separation routes, the distillation
19
process coupling with acetalization reaction was currently proved to be the best
20
strategy for the glycols system in industry due to high yield of EG and low energy
21
cost.
22
Supporting Information 20
ACS Paragon Plus Environment
Page 20 of 32
Page 21 of 32 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
1
Tables of vapor-liquid equilibrium data measured for the binary mixture composed of
2
EG, 1,2-PG, 1,2-BD and 2,3-BD are available in the Supporting Information.
3
Acknowledgements
4
The authors are grateful for the financial support from National Natural Science
5
Foundation of China (Nos. 21690084, 21336007).
6
Corresponding Author
7
*E-mail:
[email protected] (Xin Gao) Tel: +86-022-27404701
8
Notes
9
The authors declare no competing financial interest.
10
References
11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34
(1) Sun, L.; Fu, J.; Li, W.; Zhan, X. L.; Xiao, W. D. Binary vapor-liquid equilibrium of methyl glycolate and ethylene glycol. Fluid Phase Equilib. 2006, 250, 33-36. (2) Chen, Y.; Hung, S.; Lee, H.; Chien, I. Energy-Saving Designs for Separation of a Close-Boiling 1,2-Propanediol and Ethylene Glycol Mixture. Ind. Eng. Chem. Res. 2015, 54, 3828-3843. (3) Yu, B.; Chien, I. Chemical Engineering Research and Design Design and optimization of dimethyl oxalate ( DMO ) hydrogenation process to produce ethylene glycol ( EG ). Chem. Eng. Res. Des. 2017, 121, 173-190. (4) Niu, Y.; Liu, Z.; Qiao, K.; Chen, M. Study on Separation of Glycol and 1,2-Propanediol with Heterogeneous Azeotropic Distillation. Contemp. Chem. Ind. (China) 2011, 40, 560-561. (5) Song, G. Study on the Separation of Glycol-Propanediol-Butanediol System Using Azeotropic Distillation. Thesis, Tianjin University, Tianjin, 2006. (6) Dai, C.; Pan, G.; Liu, Y.; Li, Q. A novel method to separate ethylene and 1,2-butanediol with azeotropic distillation coupled with extraction technology. China Patent CN 103396290-A, 2013. (7) Li, H.; Huang, W.; Li, X.; Gao, X. Application of the Aldolization Reaction in Separating the Mixture of Ethylene Glycol and 1,2-Butanediol: Thermodynamics and New Separation Process. Ind. Eng. Chem. Res. 2016, 55, 9994-10003. (8) Figueirêdo, M. F.; Brito, K. D.; Ramos, W. B.; Vasconcelos, L. G. S.; Brito, R. P.; Brito, K. D.; Ramos, W. B.; Vasconcelos, L. G. S.; Brito, R. P., Optimization of the Design and Operation of Extractive Distillation Processes Optimization of the Design and Operation of Extractive Distillation Processes. Sep. Sci. Technol. 2015, 50, 2238-2247. (9) Xiao, J.; Zhong, L.; Guo, Y. Separation method for glycol, propylene glycol and butylene glycol. China Patent CN 102372600 A, 2012. (10) Zhu, L. Yan, J.; Xiao, W. Determination and correlation of VLE data for ethylene glycol and 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
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
1,2-butanedioi system. Chem. Eng. (China) 2012, 40, 34-37. (11) Chiavone-filho, O.; Proustj, P., Vapor-Liquid Equilibria for Glycol Ether + Water Systems. J. Chem. Eng. Data 1993, 38, 128-131. (12) Yue, H.; Zhao, Y.; Ma, X.; Gong, J. Ethylene glycol: properties, synthesis, and applications. Chem. Soc. Rev. 2012, 41, 4218. (13) Yang, Q.; Jin, F.; Feng, X.; Zhi, J.; Yang, C., Isobaric Vapor–Liquid Equilibrium for Two Binary Systems ( n -Butanol + 1,4-Butanediol and γ-Butyrolactone + 1,4-Butanediol) at p = (30.0, 50.0, and 70.0) kPa. J. Chem. Eng. Data 2016, 61, 3034-3040. (14) Zhong, Y.; Wu, Y.; Zhu, J.; Chen, K.; Wu, B.; Ji, L. Thermodynamics in separation for the ternary system 1,2-ethanediol + 1,2-propanediol + 2,3-butanediol. Ind. Eng. Chem. Res. 2014, 53, 12143-12148. (15) Yang, C.; Feng, X.; Sun, Y.; Yang, Q.; Zhi, J. Isobaric Vapor–Liquid Equilibrium for Two Binary Systems{Propane-1,2-diol + Ethane-1,2-diol and Propane-1,2-diol + Butane-1,2-diol} at p = (10.0, 20.0, and 40.0) kPa. J. Chem. Eng. Data 2015, 60, 1126-1133. (16) J.W., Kang; V., Diky; R.D. Chirico; J.W. Magee; C.D. Muzny; I. Abdulagatov; A.F. Kazakov; Michael. F. Quality assessment algorithm for vapour-liquid equilibrium data. J. Chem. Eng. Data 2010, 55, 3631-3640. (17) Herington, E. F. G. Tests for the consistency of experimental isobaric vapour-liquid equilibrium data. J. Pet. Ins. 1951, 37, 457-470. (18) Hack C. W. Winkle M. V. Vapor-Liquid Equilibria of the Diacetone Alcohol–Water System at Subatmospheric Pressures. Ind. Eng. Chem. Res., 1954, 46, 2392-2395. (19) Esteban J.; Ladero M.; García-Ochoa F. Liquid–liquid equilibria for the systems ethylene carbonate+ethylene glycol+glycerol; ethylene carbonate+glycerol carbonate+glycerol and ethylene carbonate+ethylene glycol+glycerol carbonate+glycerol at catalytic reacting temperatures. Chem. Eng.g Res. Des., 2014, 94, 440-448. (20) Mohsen-Nia M.; Doulabi F. S. M.; Manousiouthakis V. I. (Liquid + liquid) equilibria for ternary mixtures of (ethylene glycol + toluene + n-octane). J. Chem. Thermodyn. 2008, 40, 1269-1273. (21) Cháfer, A.; Burguet, M. C.; Monton, J. B.; Lladosa, E., Liquid–liquid equilibria of the systems dipropyl ether+n-propanol+water and dipropyl ether+n-propanol+ethylene glycol at different temperatures. Fluid Phase Equilib. 2007, 262, 76-81. (22) Britt, H. I.; Luecke, R. H., The Estimation of Parameters in Nonlinear, Implicit Models. Technometrics 1973, 15, 233-247. (23) Sun, L. Binary Vapor-liquid Equilibrium of Methyl Glycolate and Ethylene glycol and Separation Process for Hydrogenation of Dimethyl Oxalate to Ethylene Glycol. Thesis, East China University of Science and Technology, Shanghai, 2007. (24) Matsuda, H.; Takahara, H.; Fujino, S.; Constantinescu, D.; Kurihara, K.; Tochigi, K.; Ochi, K.; Gmehling, J., Fluid Phase Equilibria Selection of entrainers for the separation of the binary azeotropic system methanol + dimethyl carbonate by extractive distillation. Fluid Phase Equilib. 2011, 310, 166-181. (25) Ma, X.; Liu, X.; Li, Z.; Xu, G. Vapor–liquid equilibria for the ternary system methanol + dimethyl carbonate + dimethyl oxalate and constituent binary systems at different temperatures. Fluid Phase Equilib. 2004, 221, 51-56. (26) Berg, L. Recovery of ethylene glycol from butanediol isomers by azeotropic distillation. U.S. 22
ACS Paragon Plus Environment
Page 22 of 32
Page 23 of 32 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
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
Patent, US4966658. 1990. (27) Wang, Z. Separation Process for the Product of Ethylene Glycol Synthesis from Syngas. Thesis, East China University of Science and Technology, Shanghai, 2011. (28) Li, Q. Process Separation of ethylene glycol from the polyols mixture. Thesis, Jilin University, Jilin, 2014. (29) Yang, Y., Huang, M., Ma, Xiaoxun. Isobaric vapor-liquid equilibrium(VLE) data of the binary system 1,2-Ethanediol+Ethylbenzene at 3.55kPa. J. Northwest University(China), 2017. (30) King, C. J. Separaton Processes. McGraw-Hill: New York, 1971. (31) Dhale A. D.; Myrant L K; Chopade S. P. Propylene glycol and ethylene glycol recovery from aqueous solution via reactive distillation. Chem. Eng. Sci. 2004, 59, 2881-2890. (32) Huang, W.; Li, H.; Wang, R.; Li, X.; Gao, X. Chemical Engineering & Processing : Process Intensification Application of the aldolization reaction in separating the mixture of ethylene glycol and 1,2-butanediol: Kinetics and reactive distillation. Chem. Eng. Process.: Process Intensification 2017, 120, 173-183. (33) Dhale, A. D.; Myrant, L. K.; Chopade, S. P.; Jackson, J. E.; Miller, D. J. Propylene glycol and ethylene glycol recovery from aqueous solution via reactive distillation. Chem. Eng. Sci. 2004, 59, 2881-2890. (34) Chopade, S. P.; Dhale, A. D.; Clark, A. M.; Kiesling, C. W.; Myrant, L. K.; Jackson, J. E.; Miller, D. J. Vapor-Liquid-Liquid Equilibrium ( VLLE ) and Vapor Pressure Data for the Systems 2-Methyl-1,3-dioxolane ( 2MD ) + Water and 2,4-Dimethyl-1,3-dioxolane ( 24DMD )+ Water. J. Chem. Eng. Data 2003, 46, 44-47. (35) Shubham P. Chopade; Atul D. Dhale; Angela M. Clark; Chris W. Kiesling; Laurie K. Myrant; James E. Jackson.; Dennis J. Miller. Vapor−Liquid−Liquid Equilibrium (VLLE) and Vapor Pressure Data for the Systems 2-Methyl-1,3-dioxolane (2MD) + Water and 2,4 -Dimethyl-1,3-dioxolane (24DMD) + Water. J. Chem. Eng. Data 2003, 1, 44-47. (36) Chien, I Lung, B. Y. Yu, and Z. J. Ai. Design of Azeotropic Distillation Systems. Chemical Engineering Process Simulation. 2017. (37) Luyben, W. L., Comparison of extractive distillation and pressure-swing distillation for acetone/chloroform separation. Comput. Chem. Eng. 2013, 50, 1-7. (38) Dye, R. F. Ethylene glycols technology. Korean J. Chem. Eng. 2001, 18, 571-579. (39) Li, X.; Wang, R.; Na, J.; Li, H.; Gao, X., Reversible Reaction-Assisted Intensification Process for Separating the Azeotropic Mixture of Ethanediol and 1,2-Butanediol: Reactants Screening. Ind. Eng. Chem. Res. 2018, 57, 710-717.
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
470.5
Page 24 of 32
liquid phase vapor phase
470.0 469.5 469.0 468.5 468.0 467.5 467.0 0.0
0.2
0.4
0.6
0.8
1.0
Mass fraction of 1,2-BD Figure 1. Binary phase diagram of the binary system {EG +1,2-BD}.
ACS Paragon Plus Environment
Page 25 of 32 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. Flowsheet and simulation results of CD route.
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 32
1000
80
Rmin
70
Reboiler Duty 800
60 50
600
40 400
30 20
200 10 0 0.90
0.92
0.94
0.96
0.98
0
Recovery of Methanol Figure 3. Minimum reflux ratio and reboiler heat duty of Col-11 at different methanol recovery.
ACS Paragon Plus Environment
Page 27 of 32 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
1.00 0.98 0.96 0.94 0.92 0.90
Vapor phase Liquid phase
0.88 0.86
0
5
10
15
20
76
77
78
79
80
81
82
Stage number Figure 4. The concentration profile of EG in vapor and liquid phase of Col-14.
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
)LJXUH
)ORZVKHHW DQG VLPXODWLRQ UHVXOWV RI $' URXWH
ACS Paragon Plus Environment
Page 28 of 32
Page 29 of 32 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
12
70
NT,min Reboiler duty
60
50
10
8
40 6 30 4 20 2 10
0.2
0.4
0.6
0.8
1.0
EB/Feed Figure 6. Minimum stage number and reboiler heat duty of Col-27 at different ratio of EB to feed.
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 7. Flowsheet and simulation results of ED route.
ACS Paragon Plus Environment
Page 30 of 32
Page 31 of 32 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 8. Flowsheet and simulation results of RD route.
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 Graph
ACS Paragon Plus Environment
Page 32 of 32