Economics and controllability of conventional and ... - ACS Publications

comparison with a conventional extractive distillation (CED) configuration in ... extended-SEDC compared to CED can reduce TAC by 6.32%, 14.39% and ...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF DURHAM

Process Systems Engineering

Economics and controllability of conventional and intensified extractive distillation configurations for acetonitrile/methanol/benzene mixtures Chao Wang, Chen Wang, Yue Cui, Chao Guang, and Zhishan Zhang Ind. Eng. Chem. Res., Just Accepted Manuscript • DOI: 10.1021/acs.iecr.8b01875 • Publication Date (Web): 16 Jul 2018 Downloaded from http://pubs.acs.org on July 19, 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.

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 33 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

Economics and controllability of conventional and intensified extractive distillation configurations for acetonitrile/methanol/benzene mixtures Chao Wang, Chen Wang, Yue Cui, Chao Guang, Zhishan Zhang* College of Chemical and Environmental Engineering, Shandong University of Science and Technology, Qingdao 266590, China ABSTRACT: For the industrial separation problem of acetonitrile/methanol/benzene involving three azeotropes, this article proposes two novel intensified configurations with chlorobenzene as solvent by extending extractive dividing-wall column (EDWC) and side-stream extractive distillation column (SEDC) from single azeotropic to multi-azeotropic mixtures, and makes a comparison with a conventional extractive distillation (CED) configuration in aspects of economic and controllability. Steady-state designs are optimized via using the sequential iterations search based on minimum total annual cost (TAC). The results show that extended-EDWC and extended-SEDC compared to CED can reduce TAC by 6.32%, 14.39% and energy consumption by 4.26%, 20.59% respectively. As far as controllability is concerned, all proposed control structures can hold products at high purities with acceptable deviations and short settle time after introducing feed flow and composition disturbances. Overall, apparent economic benefits and high energy-efficiency provided by the extended-SEDC configuration can be achieved without a deterioration of control behavior for some multi-azeotropic mixtures. 1. INTRODUCTION In recent years, the separation problem involving two or more azeotropes in the chemical and pharmaceutical industries has arouse extensive interest in the study,1,2 and it is fairly challenging in two aspects of optimal design and dynamic control due to the system complexity. A few special distillation processes have been widely used for separating azeotropic mixtures, such as azeotropic distillation,3-7 extractive distillation,8-13 pressure-swing distillation.14-16 Pressure-swing distillation is very efficient only for pressure-sensitive mixtures and azeotropic distillation is likely to bring some challenges like multiple steady states and larger energy demand. By contrast, extractive distillation is the most promising method via introducing a solvent which can break all azeotropes and eliminate distillation boundaries for this separation issue, and also a 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

number of solvent selection methods have been proposed by some researchers,17,18 but at the same time it is also the energy intensive process. Generally, the remixing of internal streams with different compositions that occurs at feed points and along the column is an intrinsic source of thermodynamic inefficiency of this separation process, which exhibits concentration peaks of middle boiling components either above or below the feed stage.19, 20 To reduce or even eliminate the remixing effect for improving energy efficiency, extractive dividing-wall column (EDWC) and side-stream extractive distillation column (SEDC) were studied,21-27 and it has been proven that these thermally coupled configurations are hopeful alternative energy solutions. Certainly, any process intensification will inevitably bring the complicated dynamic behavior and control structure. Studies on the dynamic and control properties of EDWC have been reported by many literatures,28-33 and but there have been few reports of SEDC. Moreover, up to now most researches on extractive distillation concentrated on the separation of binary mixtures only containing one azeotrope. The focus of this work is to illustrate the possible applications of these intensified alternative configurations in the separation of multi-azeotropic mixtures with as an example of acetonitrile/methanol/benzene. To the best of our knowledge, only a few studies on pressure swing distillation has been reported for this system containing three binary azeotropes.34-38 In this work, two novel thermal coupled configurations, namely extended-EDWC and extended-SEDC, were developed for this specific separation with chlorobenzene as a solvent, and compared with the configuration of conventional extractive distillation (CED) from two aspects of economic and controllability. The optimization of sequential iterative search was carried out to minimize total annual cost (TAC). The steady state study showed promising results in economic benefits for two proposed configurations. The control structures of three configurations were established, and they represented satisfactory results in the feed flow and composition disturbances rejection capability. In addition, a comparison among them was presented to further evaluate operational aspects of extended EDWC and extended-SEDC implementations. 2. DESIGN BASIS 2.1 Thermodynamics. The wilson thermodynamic model is chosen to describe the vapor-liquid equilibrium in the system being studied because of the agreement between the predicted and experimental data in the literatures34,35, 39-41. All binary interaction parameters built in the Aspen 2

ACS Paragon Plus Environment

Page 2 of 33

Page 3 of 33

Plus (V8.4) are represented as a supplementary material in Supporting Information. From the ternary composition diagram shown in Figure 1, there are three binary homogenous azeotropes in the system of acetonitrile/methanol/benzene.

Figure 1. The ternary composition diagram for the acetonitrile/methanol/benzene system. Chlorobenzene is selected as a solvent for this separation process of extractive distillation. Figure 2 describes vapor-liquid equilibrium behaviors in connected with chlorobenzene at 1 atm. It can be seen that all azeotropes are broken at the solvent-to-feed ratio (S:F) of 1.0, and the separation in a solvent recovery column is straightforward, and the light-heavy component order becomes methanol-acetonitrile-benzene.

(a) 1

(b) 1 S:F=0 S:F=1

0.8

0.6

0.6

S:F=0 S:F=1

ymethanol

0.8

ymethanol

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

Industrial & Engineering Chemistry Research

0.4

0.4

0.2

0.2

0

0 0

0.2

0.4 0.6 xmethanol

0.8

0

1

0.2

3

ACS Paragon Plus Environment

0.4

0.6 xmethanol

0.8

1

Industrial & Engineering Chemistry Research

(c) 1

(d) 1 methanol

S:F=0 S:F=1

0.8

acetonitrile

0.8 ychlorobenzene

benzene

ybenzene

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 33

0.6 0.4 0.2

0.6 0.4 0.2

0

0 0

0.2

0.4 0.6 xbenzene

0.8

1

0

0.2

0.4 0.6 xchlorobenzene

0.8

1

Figure 2. The y-x plots of different components: (a) methanol-benzene, (b) methanol-acetonitrile, (c) benzene-acetonitrile, and (d) chlorobenzene-different component. 2.2 Economics. In this article, TAC is used as the objective function to be minimized by adopting the global sequential iterative optimization methods to screen the optimum design, including annual capital costs and operating costs. Major pieces of equipment are column vessels (including column internals) and heat exchangers (condenser and reboiler). Small items such as reflux drums, pumps, valves, and pipes are usually not considered at the conceptual design stage because of their lower costs compared with column vessels and heat exchangers. The costs of cooling water and solvent makeup are not included in operating costs because they are much lower than the costs of heat duties. The specific calculation procedure for the economic evaluation42-44 is given in Supporting Information. 3. STEADY-STATE DESIGN The commercial software Aspen Plus V8.4 is applied to simulate these separation processes of CED, extended-EDWC and extended-SEDC. A case study presented in the reference35 is used as the basis for our study. In that reference, a fresh feed is 1000 kg/h (28 kmol/h) with composition 70 wt% (78 mol%) methanol, 20 wt% (17.4 mol%) acetonitrile and 10 wt% (4.6 mol%) benzene. The product purities are specified as 99.5 mol% methanol, 99.5 mol% acetonitrile, 99.8 mol% benzene, and 99.995 mol% solvent, respectively. In the economic optimization, the number of stages (NT), feed location (NF), solvent location (NFS), side stream location (NSD) and solvent flow (SF) are chosen as optimal variables. The key design variables including reflux ratios, distillate flowrates and side stream flowrates are adjusted to meet the specified purity and recovery 4

ACS Paragon Plus Environment

Page 5 of 33 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

of key components. Additionally, the temperature of the recycle solvent and the operation pressure of distillation column are often fixed to a constant and not optimized in most of the literatures to date.45-48 Therefore, the temperature of the recycle solvent is fixed at 320 K and the operation pressure is set at 1 atm for simplification in this work, which is line in the work of Luyben.49,50 3.1. Design of CED. Figure 3 gives the flowsheet of the CED process with two extractive distillation columns and a solvent recovery column. The fresh feed is fed to the first extractive distillation column (EDC1) from the lower tray, and one recycle stream of solvent is fed to EDC1 from the upper tray. The distillate of EDC1 is the methanol product. The bottom from EDC1 is fed to the second extractive distillation column (EDC2) from the lower tray, and meanwhile the other recycle stream of solvent is fed to EDC2 from the upper tray. The distillate of EDC2 is the acetonitrile product. The bottom from EDC2 enters the solvent recovery column (SRC). The benzene product is obtained at the top of SRC. The high-purity solvent from the base of SRC after being cooled is split and recycled back to EDC1 and EDC2. A very small solvent makeup stream is required because of very small solvent losses in three product streams.

Figure 3. The optimal flowsheet of the CED configuration. The sequential iterative optimization procedure for the CED configuration is shown in Supporting Information. The key design variables including reflux ratios (RR1, RR2 and RR3) and distillate flowrates (D1, D2 and D3) of EDC1, EDC2 and SRC are adjusted to reach the separation requirements. There are still ten optimization variables in this configuration as follows: the number of stages (NT1, NT2 and NT3) and the fresh feed locations (NF1, NF2 and NF3) of EDC1, 5

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

EDC2 and SRC, the solvent locations (NFS1 and NFS2) and the solvent flowrates (SF1 and SF2) of EDC1 and EDC2. The vapor and liquid composition profiles of three columns in the CED configuration are shown in Figure 4. It can be clearly seen from Figure 4a and Figure 4b that the purities of acetonitrile and benzene achieve the maximum value on a certain tray near the bottom of EDC1 or EDC2, after that tray they decrease toward the bottom of EDC1 or EDC2. The effect of remixing with the heaviest solvent increases the energy consumption of the following separation process in order to repurify the acetonitrile and benzene. From Figure 4c, the distillate is the high purity benzene and the bottom is the high purity solvent in SRC.

(a) 1

1

Vapor mole fraction

0.6

Liquid mole frction

Methanol Benzene Acetonitrile Solvent

0.8

0.4 0.2 0

Methanol Benzene Acetonitrile Solvent

0.8 0.6 0.4 0.2 0

0

10

20

30 Stage

40

50

0

(b) 1

10

20

30 Stage

40

50

1

0.6

Liquid mole fraction

Methanol Benzene Acetonitrile Solvent

0.8

Vapor mole fraction

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 33

0.4 0.2 0 0

Methanol Benzene Acetonitrile Solvent

0.8 0.6 0.4 0.2 0

5 10 15 20 25 30 35 40 45 Stage

0

5

10 15 20 25 30 35 40 45 Stage

6

ACS Paragon Plus Environment

Page 7 of 33

1

(c) 1 Methanol Benzene Acetonitrile Solvent

0.8 0.6

Liquid mole fraction

Vapor mole fraction

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

Industrial & Engineering Chemistry Research

0.4 0.2 0

Methanol Benzene Acetonitrile Solvent

0.8 0.6 0.4 0.2 0

0

5

10

15 20 Stage

25

30

35

0

5

10

15

20 Stage

25

30

35

Figure 4. Vapor and liquid composition profiles of the CED configuration: (a) EDC1, (b) EDC2, and (c) SRC. 3.2. Design of extended-EDWC. Figure 5 shows the flowsheet of the newly proposed extended-EDWC configuration. This complex configuration features an extended-EDWC in which extractive distillation is carried out in both main column (MC) and side column (SC). The distillates of MC and SC are the methanol product and the acetonitrile product, respectively. The bottom from MC, which is mostly benzene and solvent, is fed to the solvent recovery column (SRC). The distillate of SRC is the benzene product. The high-purity solvent from the base of SRC after being cooled is split and recycled back to MC and SC.

Figure 5. The optimal flowsheet of the extended-EDWC configuration. The sequential iterative optimization procedure for the extended-EDWC configuration is given in Supporting Information. The key design variables including reflux ratios (RR1, RR2 and RR3) of 7

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

MC, SC and SRC, distillate flowrates (D1 and D3) of MC and SRC, and side vapor flowrates (VR) of MC are adjusted to reach the separation requirements. There are still ten optimization variables in this configuration as follows: the number of stages (NT1, NT2 and NT3) of MC, SC and SRC, the feed locations (NF1 and NF3) in MC and SRC, the side vapor withdrawn location (NSD) of MC, the solvent flowrates (SF1 and SF2) and the solvent locations (NFS1 and NFS2) of MC and SC. Figure 6 describes the vapor and liquid composition profiles of the extended-EDWC configuration. As shown in Figure 6a and Figure 6b, the high purity methanol can be obtained at the top of MC and the acetonitrile product can be obtained at the top of the SC. Moreover, the remixing phenomenon is alleviated, which results in the reduction of energy consumption. From Figure 6c, the high purity benzene and recycle solvent are obtained in the distillate and the bottom of SRC, respectively.

(a) 1

1

Vapor mole fraction

0.8 0.6

Liquid mole fraction

Methanol Benzene Acetonitrile Solvent

0.4 0.2

Methanol Benzene Acetonitrile Solvent

0.8 0.6 0.4 0.2

0

0 10

20

30 40 Stage

(b) 1

50

60

Methanol Benzene Acetonitrile Solvent

0.8 0.6 0.4 0.2

0

10

20

30 40 Stage

50

60

1 Liquid mole fraction

0

Vapor mole fraction

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 33

Methanol Benzene Acetonitrile Solvent

0.8 0.6 0.4 0.2 0

0 0

10

20Stage 30

40

0

50

10

8

ACS Paragon Plus Environment

20

Stage

30

40

50

Page 9 of 33

(c) 1

1 Methanol Benzene Acetonitrile Solvent

Methanol

0.8

Benzene

Liquid mole fraction

Vapor mole fraction

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

Acetonitrile

0.6

Solvent

0.4 0.2 0 0

5

10

15 Stage

20

25

0.8 0.6 0.4 0.2 0

30

0

5

10

15 Stage

20

25

30

Figure 6. Vapor and liquid composition profiles of the extended EDWC configuration: (a) MC, (b) SC, and (c) SRC. 3.3. Design of extended-SEDC. Figure 7 gives the flowsheet of the newly proposed extended-SEDC process. This complex configuration features two side-stream extractive distillation columns. The methanol product is removed from the top of the first side-stream extractive distillation column (SEDC1). A liquid side stream (LS1) is withdrawn near the bottom of SEDC1 and fed to the second side-stream extractive distillation column (SEDC2), which contains all benzene and acetonitrile and partial solvent. The distillate of SEDC2 is the acetonitrile product. A liquid side stream (LS2) is withdrawn near the bottom of SEDC2 and fed to the solvent recovery column (SRC), which contains all benzene and partial solvent. The benzene product is obtained at the top of SRC. The high-purity solvent from the base of three columns after being cooled is split and recycled back to SEDC1 and SEDC2.

9

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

Figure 7. The optimal flowsheet of the extended-SEDC configuration. The sequential iterative optimization procedure for the extended-SEDC configuration is shown in Supporting Information. The key design variables including the reflux ratios (RR1, RR2 and RR3) and the distillate flowrates (D1, D2 and D3) of SEDC1, SEDC2 and SRC, and the side liquid flowrates (LS1 and LS2) of SEDC1 and SEDC2 are adjusted to reach the separation requirements. There are still twelve optimization variables in this configuration as follows: the number of stages (NT1, NT2 and NT3) and the feed locations (NF1, NF2 and NF3) of SEDC1, SEDC2 and SRC, the solvent locations (NFS1 and NFS2) and the solvent flowrates (SF1 and SF2) of SEDC1 and SEDC2 and the side liquid withdrawn locations (NSD1 and NSD2) of SEDC1 and SEDC2. The vapor and liquid composition profiles of the extended-SEDC configuration are shown in Figure 8. As shown in Figure 8a and Figure 8b, three high purity products can be obtained at the top of the SEDC1, SEDC2 and SRC, respectively. Moreover, the remixing effect is apparently alleviated, which results in the dominant reduction of energy consumption. From Figure 8c, the high purity benzene and solvent can be obtained in the distillate and bottom of SRC, respectively.

(a) 1

1

Vapor mole fraction

0.8 0.6

Liquid mole fraction

Methanol Benzene Acetonitrile Solvent

0.4 0.2

Methanol Benzene Acetonitrile Solvent

0.8 0.6 0.4 0.2 0

0 0

10

20

30 40 Stage

50

0

60

(b) 1

10

20

30 40 Stage

50

60

1 Methanol Benzene Acetonitrile Solvent

0.6

0.8

Liquid mole fraction

0.8 Vapor mole fraction

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 10 of 33

Methanol Benzene Acetonitrile Solvent

0.6 0.4

0.4

0.2

0.2

0

0 0

10

20

30 Stage

40

50

0

60

10

10

ACS Paragon Plus Environment

20

30 Stage

40

50

60

Page 11 of 33

(c) 1 Methanol Benzene Acetonitrile Solvent

0.8 0.6 0.4 0.2

Liquid mole fraction

1

Vapor mole fraction

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

Industrial & Engineering Chemistry Research

0 0

5

10

15 20 Stage

25

Methanol Benzene Acetonitrile Solvent

0.8 0.6 0.4 0.2 0

30

0

5

10

15 20 Stage

25

30

Figure 8. Vapor and liquid composition profiles for the extended-SEDC configuration: (a) SEDC1 (b) SEDC2, and (c) SRC. Table 1 presents the head-to-head comparison of three configurations in terms of optimized results and economic costs. It can be seen that two proposed thermal coupled configurations have obvious advantages in economic benefits and energy consumption compared with the CED configuration. This high energy efficiency in two intensified configurations reflects the reduction or elimination of the remixing effect, especial for the extended-SEDC configuration. The reason for this difference is that two side liquid streams (LS1 and LS2) in SEDC1 and SEDC2 can reduce the remixing effect in the compositions of acetonitrile and benzene, but one vapor split (VR) in MC can reduce partly the remixing effect in the compositions of acetonitrile and not work for the component of benzene. Table 1. Optimal Results for Three Separation Configurations CED

extended-EDWC

extended-SEDC

EDC1

EDC2

SRC

MC

SC

SRC

SEDC1

SEDC2

SRC

NT

50

44

35

56

48

33

58

57

31

NFS

29

21



27

21



29

24



NF

43

33

12

41

49

10

43

35

11

NSD







49





50

43



ID(m)

0.82

0.58

0.67

0.664

0.949

0.541

0.236

RR

2.09

2.43

43.94

2.12

1.86

43.83

2.11

2.11

5.3

Side stream (kmol/h)







18.43





25.53

15.13



SF (kmol/h)

100

40



105

33



100

53



1.13

11

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

Capital cost (105$/a)

1.87

1.00

1.21

2.57

1.18

2.22

1.09

0.34

Energy cost (105$/a)

1.48

0.52

0.73

1.91

0.72

1.55

0.53

0.10

TAC (105$/a)

3.35

1.52

1.94

4.48

1.90

3.77

1.62

0.44

Total (105$/a)

6.81

6.38

5.83

4. DYNAMIC CONTROL The proposed configurations in this paper have the advantages of more energy efficiency and economic benefits than other distillation alternatives in separating the multi-azeotropic system, but control complexity and less robust dynamic behaviors are undesirable side effects since more control variables are introduced. Therefore, to examine operational aspects of three configurations, their control structures are developed and the dynamic performances are tested when disturbance occur in feed flowrate and composition. For the control implementation, the steady state design is exported to a pressure-driven simulation in Aspen Dynamics. Reflux drum and base volume are specified to provide 5 min of holdup for 50% liquid level. Pumps and valves are sized to give proper pressure drops to handle changes in flowrate.42 4.1 Control of CED. In order to select properly a temperature control location, the open loop sensitivity analysis is carried out with ±0.1% changes in the reboiler duty (QR1, QR2 and QR3) of EDC1, EDC2 and SRC and with ±0.1% changes in the reflux ratio (RR3) of SRC, respectively. As shown in Figure 9, the temperature control point of each column is selected as the following: the 46th tray for EDC1 (T1,46), the 36th tray for EDC2 (T2,36), and the 7th tray for SRC (T3,7). Note that only one temperature control tray that occurs in the rectifying can be determined by using the sensitivity criteria of reflux ratio (RR3) and reboiler duty (QR3), as shown in Figure 9c and Figure 9d, thus one composition controller has to be installed to control the impurity in the bottom.

(b) 12

(a) 60

+0.1%QR2

+0.1%QR1 -0.1%QR1

40

8

-0.1%QR2

4 Gain

20

Gain

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 33

0

0

-4

-20

-8

-40

-12 0

10

20

30 Stage

40

50

60

0

10

12

ACS Paragon Plus Environment

20

30 Stage

40

50

Page 13 of 33

(c) 60

(d) 60 +0.1%RR3

40

-0.1%RR3

20

+0.1%QR3

40

-0.1%QR3

20 Gain

Gain

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

Industrial & Engineering Chemistry Research

0

0

-20

-20

-40

-40

-60

-60 0

5

10

15

20 25 Stage

30

35

40

0

5

10

15 20 Stage

25

30

35

Figure 9. The open-loop sensitivity analysis of the CED configuration: (a) EDC1, (b) EDC2 and (c, d) SRC. Figure 10 shows the control structure developed for the CED configuration. Flow controllers are set as Proportional-Integral (PI) controller with a gain (KC) of 0.5 and an integral time (τ1) of 0.3min. All pressure controllers are tightly tuned with KC =20 and τ1 = 12 min. All level controller are proportional with KC = 2 and τ1 = 9999 min. All temperature and composition controllers have deadtimes of 1 min and 3 min, respectively, and are tuned using relay-feedback testing and Tyreus-Luyben tuning rules.42 All control loops are described as follows: (1) Fresh feed is flow controlled. (2) The column pressure is controlled by manipulating the condenser duty. (3) The base level of EDC1 and EDC2 are held by manipulating the bottom flowrates, and that of SRC is held by manipulating the small solvent makeup. (4) Reflux drums levels of EDC1 and EDC2 are held by manipulating the distillate flowrates respectively, and that of SRC is held by manipulating the reciprocal of reflux flowrate because of the high reflux ratio. (5) All reflux ratios are kept constant. (6) Two proportional controllers that solvent flowrates are ratioed to feed flowrates are introduced because the recycle solvent needs to split properly for two extraction distillation processes. (7) The control point T1,46 is controlled by manipulating heat input to the EDC1 reboiler. (8) The control point T2,36 is controlled by manipulating heat input to the EDC2 reboiler. (9) The control point T3,7 is controlled by manipulating the reciprocal of the reflux ratio of SRC. 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

(10) The benzene composition in the bottoms of SRC is controlled by manipulating heat input to the SRC reboiler. (11) The temperature of the recycle solvent is controlled by manipulating the cooling water flowrate to the solvent cooler.

Figure 10. The overall control structure of the CED configuration. In order to demonstrate the effectiveness of the developed control structure, a small magnitude of 5% disturbances in the feed flow is introduced since valves saturation is easy to occur in the control of the complex extraction distillation process. The magnitude of the disturbances in the feed composition is similar with the reported literatures,51, 52 that is, increasing from 20 % to 22 wt % acetonitrile, 10 % to 11 wt % benzene, and decreasing from 20 % to 18wt % acetonitrile, 10 % to 9 wt % benzene. Figure 11a gives the close loop responses for feed flow disturbances. All temperature controllers responded relatively fast and the product compositions are maintained quite close to the specified values. The time required to reach the steady state again is approximately 5 h. 14

ACS Paragon Plus Environment

Page 14 of 33

Page 15 of 33

Figure 11b gives the close loop response results for feed composition disturbances. All temperature controllers handled the disturbance reasonably well by bringing the temperatures back to their setpoints. The purities of acetonitrile and benzene after 5h have a small deviation from their desired values. Overall, the proposed control structure shows good dynamic behavior when facing feed flow and composition disturbances.

99

133

95

131

T2,36(℃)

T1,46(℃)

(a)

129

91

127

87 0

5

10 Time/h

15

0

20

xD1,methanol(mol%)

102

10 Time/h

15

20

0.9955

101 T3,7(℃)

5

0.996

103

100 99

0.995

0.9945

98

0.994

97 0

5

10 Time/h

15

0

20 1

0.9955

0.999

xD3,benzene(mol%)

0.996 xD2,acetonitrile(mol%)

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

Industrial & Engineering Chemistry Research

0.995

0.9945 0.994

5

10 Time/h

15

20

15

20

0.998 0.997 0.996 0.995

0

5

10 Time/h

15

20

0

(b)

15

ACS Paragon Plus Environment

5

10 Time/h

134 T2,36(℃)

132 130 128 126 124

0

5

10 Time/h

15

20

0

0.996

101

0.9955

xD1,methanol(mol%)

102

T3,7(℃)

Page 16 of 33

136

99 97 95 93 91 89 87 85

100 99

5

10 Time/h

15

20

0.995

0.9945

98 0

5

10 Time/h

15

0.994

20

0

5

10 Time/h

15

20

0

5

10 Time/h

15

20

1 0.999

xD2,acetonitrile(mol%)

0.999

xD3,benzene(mol%)

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

T1,46(℃)

Industrial & Engineering Chemistry Research

0.998

0.996

0.997

0.993

0.996

0.99

0.995 0

5

10 Time/h

15

20

Figure 11. The dynamic responses for disturbances of the CED configuration: (a) feed flow (b) feed composition. 4.2 Control of extended-EDWC. As shown in Figure 12, temperature control locations of MC and SRC are selected also by the open loop sensitivity analysis, which is carried out with ±0.1% changes in the reboiler duty (QR1 and QR2) of MC and SRC, and with ±0.1% changes in the reflux ratio (RR2 and RR3) of SC and SRC, respectively. MC have two temperature control points of the 43th and 50th tray (T1,43, T1,50), and the temperature control point of SRC is the 6th tray (T3,6). However, it is hard to identify the temperature control point for SC based on Figure 12b, showing the open loop sensitivity analysis with ±0.1% changes in the reflux ratio (RR2) of SC. Therefore, the slope criterion shown in Figure 12e is used for selecting a temperature control 16

ACS Paragon Plus Environment

Page 17 of 33

location of SC, that is, the 36th tray (T2,21).

(a) 30

(b) 0.2 +0.1%QR1

20

0

0

Gain

Gain

0.1

-0.1%QR1

10

-0.1

-10

+0.1%RR2

-0.2

-20 -30 10

20

30 40 Stage

50

0

60

(c) 80

10

20

30 40 Stage

50

60

(d) 60

60

+0.1%RR3

40

40

-0.1%RR3

20

20

Gain

Gain

-0.1%RR2

-0.3 0

0

+0.1%QR2 -0.1%QR2

0

-20

-20

-40

-40 -60

-60 0

5

10

15 20 Stage

25

30

35

0

5

10

15 Stage

20

25

30

(e) 120 110 T(℃)

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

100 90 80 70 0

10

20

30 Stage

40

50

Figure 12. The temperature sensitive tray of the extended-EDWC configuration: (a) open-loop sensitivity analysis for MC, (b) open-loop sensitivity analysis for SC, (c, d) open-loop sensitivity analysis for SRC and (e) temperature profile of SC. Figure 13 shows the control structure developed for the extended-EDWC configuration. There are several differences except for basic control loops compared to that of traditional extractive distillation as follows: (1) A cascade combination of composition and temperature control is used in SC, which achieves both fast control and the maintenance of acetonitrile product purity. (2) The temperature control point of T1,43 is controlled by manipulating the ratio of the side vapor 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

flow (VR) to the reboiler heat duty (QR1), with the reflux ratio (RR1) being held constant.

Figure 13. The overall control structure of the extended-EDWC configuration. In order to demonstrate the effectiveness of the developed control structure, the same feed disturbances as the CED configuration are introduced. Figure 14a gives the close loop responses for feed flow disturbances. All temperature control points return quite effectively to their set points and all the product compositions are maintained quite close to the specified values. The time required to reach the steady state again is approximately 5 h. Figure 14b gives the close loop responses for feed composition disturbances. All temperature controllers also handle the disturbance reasonably well by bringing the temperatures back to their setpoints. The purity of acetonitrile after 5h has a small deviation from its desired value. These results demonstrate fully the proposed control structure can achieve the controllable operating of the extend-EDWC configuration. 18

ACS Paragon Plus Environment

Page 18 of 33

Page 19 of 33

(a) 113

92

112

T1,43(℃)

T1,50(℃)

90

111

88

110

86

109 108

84 0

5

10 Time/h

15

20

0

111

104.8

110

T2,21(℃)

104.4

15

20

109

104

108

103.6 103.2

107 0

5

10 Time/h

15

20

0

0.998

5

10 Time/h

15

20

15

20

0.998 xD2,acetonitrile(mol%)

xD1,methanol(mol%)

10 Time/h

T3,6(℃)

105.2

5

0.996

0.994

0.992

0.996

0.994

0.992 0

5

0

5

10 Time/h

15

20

0

1 xD3,benzene(mol%)

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

Industrial & Engineering Chemistry Research

0.999 0.998 0.997 0.996 10 Time/h

15

20

19

ACS Paragon Plus Environment

5

10 Time/h

Industrial & Engineering Chemistry Research

(b) 110.6

94

110.5

92

T1,50(℃)

T1,43(℃)

110.4

90

110.3

88

110.2

86

110.1

84

110 5

10 Time/h

15

0

20

105

111

104.6

110

104.2

10 Time/h

15

20

109

103.8

108

103.4 103

107 0

5

10 Time/h

15

20

0

0.998

5

10 Time/h

15

20

1 xD2,acetonitrile(mol%)

xD1,methanol(mol%)

5

T3,6(℃)

T2,21(℃)

0

0.996

0.994

0.998 0.996 0.994 0.992 0.99

0.992 0

5

10 Time/h

15

0

20

5

10 Time/h

15

20

1

xD3,benzene(mol%)

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 33

0.999 0.998 0.997 0.996 0

5

10 Time/h

15

20

Figure 14. The dynamic responses for disturbances of the extended-EDWC control structure: (a) feed flowrate and (b) feed composition. 20

ACS Paragon Plus Environment

Page 21 of 33

4.3 Control of extended-SEDC. Similarly, the open loop sensitivity analysis is carried out with ±0.1% changes in the reboiler duty (QR1, QR2 and QR3) of EDC1, EDC2 and SRC and with ±0.1% changes in the reflux ratio (RR3) of SRC, respectively, for the temperature control tray. It can be seen from Figure 15 that the temperature control tray of SEDC1 is the 47th tray (T1,47), and that of SEDC2 is the 37th tray (T1,37) and that of SRC is the 7th tray (T3,7). Figure 16 shows the overall control structure. Special instructions except for basic control loops are required as follows: (1) A cascade combination of composition-flow control (CC3-FC2) is installed to keep the methanol composition in the side stream LS1 constant, and as well a cascade combination of composition-flow control (CC4-FC3) is settled to keep the acetonitrile composition in the side stream LS2 constant. (2) The CB impurity in the distillate of SEDC1 is controlled by the reflux ratio (RR1), and that of SEDC2 is controlled by the reflux ratio (RR2). (3) The solvent flow (SF2) is ratioed to the initial feed flow (F1).

(a) 300

(b) 40 +0.1%QR1 -0.1%QR1

200

20 0

0

Gain

Gain

100 -100 -200 -300

-20

+0.1%QR2

-40

-0.1%QR2

-60 0

10

20

30 40 Stage

50

60

0

10

20

30 40 Stage

50

60

(d) 300

(c) 300 200

+0.1%RR3

200

+0.1%QR3

100

-0.1%RR3

100

-0.1%QR3

Gain

Gain

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

Industrial & Engineering Chemistry Research

0

0

-100

-100

-200

-200

-300 0

10

20 Stage

30

-300

40

0

10

20 Stage

30

40

Figure 15. The open-loop sensitivity analysis of the extended-SEDC: (a) SEDC1, (b) SEDC2, and (c, d) SRC.

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

Figure 16. The overall control structure of the extended-SEDC configuration. The effectiveness of the developed control structure is demonstrated by introducing the same feed disturbances as the CED configuration as well. Figure 17a gives the close loop responses for feed flow disturbances. All temperature control point return quite effectively and all the product compositions are maintained quite close to the specified values. The time required to reach the steady state again is approximately 15 h. Figure 17b gives the close loop responses for feed composition disturbances. All temperature controllers handle the disturbance reasonably well by bringing the temperatures back to their setpoints. The purity of acetonitrile after 10h has an acceptable deviation from its desired values. Likewise, these results prove that the extend-SEDC configuration can implement with the good controllability based on the proposed control structure.

(a)

22

ACS Paragon Plus Environment

Page 22 of 33

117.48

98.4

117.04 T2,37(℃)

98.6

98.2 98 97.8 5

10 Time/h

15

116.6 116.16 115.72

20

0

93.6

0.996

93.4

0.9955

xD1,methanol(mol%)

T3,7(℃)

0

93.2 93 92.8 0

5

10 Time/h

15

20

15

20

5

10 Time/h

15

20

15

20

0.994 0

0.9985

xD3,benzene(mol%)

0.9955

0.9945

10 Time/h

0.9945

0.999

0.995

5

0.995

0.996 xD2,acetonitrile(mol%)

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

T1,47(℃)

Page 23 of 33

0.998 0.9975 0.997

0.994 0

5

10 Time/h

15

0

20

(b)

23

ACS Paragon Plus Environment

5

10 Time/h

Industrial & Engineering Chemistry Research

116.8

98.24

116.7 T2,37(℃)

T1,47(℃)

98.2 98.16 98.12

116.6 116.5 116.4

98.08

116.3 0

5

10 Time/h

15

20

0

0.996

93.6

0.9955

T3,7(℃)

xD1,methanol(mol%)

94

93.2 92.8

5

10 Time/h

15

20

0.995 0.9945 0.994

92.4 0

5

10 Time/h

15

20

0.997

0.999

0.996

0.9985

xD3,benzene(mol%)

xD2,acetonitrile(mol%)

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 33

0.995 0.994

0

5

10 Time/h

15

20

0

5

10 Time/h

15

20

0.998

0.9975

0.993

0.997 0

5

10 Time/h

15

20

Figure 17. The dynamic responses for disturbances of the extended-SEDC configuration: (a) feed flowrate and (b) feed composition. The proposed control structures for three complex configurations are effective in managing these operations to achieve the separation objectives. But there are some differences between them that need to be highlighted. First, the control structure for the CED configuration uses only one composition controllers that is slow and expensive, thus it has the advantages of fast responses and economy in the industrial application. Second, the control structures for both intensified configurations have more degrees of freedom, and adopt more composition controllers to maintain the product purities and complex cascade controls for improving the control quality. At last, the control structure for the extended-SEDC configuration is a bit less responsive to feed disturbances 24

ACS Paragon Plus Environment

Page 25 of 33 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

than other two control structures, and this is not to be expected and needs to be further improved. The data for all temperature and composition controllers of three control structures, such as transmitter ranges, tuning parameters and controller output ranges can be found in Supporting Information. 5. CONCLUSION

This article demonstrates three newly proposed extractive distillation configurations including CED, extended-EDWC and extended-SEDC for the separation of acetonitrile/methanol/benzene with chlorobenzene as solvent. The optimization results of the steady-state design show that two intensified configurations can provide great savings in TAC and energy consumption compared to the CED configuration, especially for the extended-SEDC. Dynamic control studies represent that they are also comparable to CED in terms of controllability under the same feed disturbances, although requiring more complex control structures. To sum up, all of the proposed configurations in this article are the essential alternatives without any potential dynamic control problems for separating the similar multi-azeotropic mixtures with acetonitrile/methanol/benzene, and particularly attractive is the extended-SEDC flowsheet from the energetic and economic standpoints. ■ ASSOCIATED CONTENT S

Supporting Information

The Supporting Information is available free of charge on the ACS publications website at DOI: (1) Correlation parameters for Wilson model for acetonitrile/methanol/benzene. (2) The necessary formula and parameter for economic evaluation. (3) Temperature-composition controllers tuning parameters. (4) The sequential iterative optimization procedure for the CED configuration. (5) The sequential iterative optimization procedure for the extended-EDWC configuration. (6) The sequential iterative optimizaton procedure for the extended-SEDC configuration. ■ AUTHOR INFORMATION *E-mail: [email protected]. ORCID Zhishan Zhang: 0000-0001-8291-4727 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

Notes The authors declare no competing financial interests. ■ NOMENCLATURE AC = heat transfer area of condenser (m2) AR = heat transfer area of reboiler (m2) CED = conventional extractive distillation CSCinst = installed cost of column shell ($) EDC = extractive distillation column EDWC = extractive dividing-wall column extended-EDWC = the configuration of extended extractive dividing-wall column extend-SEDC = the configuration of extended side-stream extractive distillation column F = feed flow (kmol/h) Hcol = column height (m) Htray = column height between top and bottom tray (m) ID = column diameter (m) KC = gain coefficient KU = ultimate gains LS = side liquid flow(kmol/h) NF = fresh feed location NFS = solvent location NSD = side stream location NT = the number of stages PU = periods QC = condenser duty (kW) QR = reboiler duty (kW) RR = reflux ratio SEDC = side-stream extractive distillation column S:F = solvent to feed ratio SF = solvent flow (kmol/h)

26

ACS Paragon Plus Environment

Page 26 of 33

Page 27 of 33 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

SRC = solvent recovery column T = stage temperature (℃) TAC = total annual cost ($/a) TCinst = installed cost of trays ($) U = heat transfer coefficient (kW/(K·m2)) VR = side vapor flow (kmol/h) x = liquid mole fraction τI = integral time constant (min)

■ REFERENCES (1) Zhao, Y.; Zhao, T.; Jia, H.; Li, X.; Zhu, Z.; Wang, Y. Optimization of the composition of mixed entrainer for economic extractive distillation process in view of the separation of tetrahydrofuran/ethanol/water ternary azeotrope. J. Chem. Technol. Biotechnol. 2017, 92, 2433-2444. (2) Wang, Y.; Zhang, X.; Liu, X.; Bai, W.; Zhu, Z.; Wang, Y.; Gao, J. Control of extractive distillation process for separating heterogeneous ternary azeotropic mixture via adjusting the solvent content. Sep. Purif. Technol. 2018, 191, 8-26. (3) Abu-Eishah, S.I.; Luyben, W.L. Design and control of a two column azeotropic distillation system. Ind. Eng. Chem. Process Des. Dev. 1985, 24, 132-140. (4) Thong, D.Y-C.; Jobson, M. Multicomponent homogeneous azeotropic distillation 1. Assessing product feasibility. Chem. Eng. Sci. 2001a, 56, 4369-4391. (5) Thong, D.Y-C.; Jobson, M. Multicomponent homogeneous azeotropic distillation 2. Column design. Chem. Eng. Sci. 2001b, 56, 4393-4416. (6) Thong, D.Y-C.; Jobson, M. Multicomponent homogeneous azeotropic distillation 3. Column sequence synthesis. Chem. Eng. Sci. 2001c, 56, 4417-4432. (7) Shi, P.; Xu, D.; Ding, J.; Wu, J.; Ma, Y.; Gao, J.; Wang, Y. Separation of azeotrope (2,2,3,3-tetrafluoro-1-propanol

+

water)

via

heterogeneous

azeotropic

distillation

by

energy-saving dividing-wall column: Process design and control strategies. Chem. Eng. Res. Des. 2018, 135, 52-66. (8) Yi, C-C.; Huang, W-C.; Chien, I-L. Energy-efficient heterogeneous extractive distillation system for the separation of close-boiling cyclohexane/cyclohexene mixture. J. Taiwan. Inst. 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

Chem. E. 2018, 87, 26-35. (9) Liang, K.; Li, W.; Luo, H.; Xia, M.; Xu, C. Energy-efficient extractive distillation process by combining preconcentration column and entrainer recovery column. Ind. Eng. Chem. Res. 2014, 53, 7121-7131. (10) Liao, B.; Lei, Z.; Xu, Z.; Zhou, R.; Duan, Z. New process for separating propylene and propane by extractive distillation with aqueous acetonitrile. Chem. Eng. J. 2001, 84, 581-586. (11) Lei, Z.; Zhou, R.; Duan, Z. Process improvement on separation C4 by extractive distillation. Chem. Eng. J. 2002, 85, 379-386. (12) Lei, Z.; Li, C.; Chen, B. Extractive distillation: a review. Sep. Purif. Rev. 2003, 32, 121-213. (13) Kossack, S.; Kraemer, K.; Gani, R.; Marquardt, W. A systematic synthesis framework for extractive distillation processes. Chem. Eng. Res. Des. 2008, 86, 781-792. (14) Zhang, Z.; Zhang, Q.; Li, G.; Liu, M.; Gao, J. Design and control of methyl acetate-methanol separation via heat-integrated pressure-swing distillation. Chinese. J. Chem. Eng. 2016, 24, 1584-1599 (15) Wang, Y.; Bu, G.; Wang, Y.; Zhao, T.; Zhang, Z.; Zhu, Z. Application of a simulated annealing algorithm to design and optimize a pressure-swing distillation process. Comput. Chem. Eng. 2016, 95, 97-107. (16) Liang, S.; Cao, Y.; Liu. X.; Li, X.; Zhao, Y.; Wang, Y.; Wang, Y. Insight into pressure-swing distillation from azeotropic phenomenon to dynamic control. Chem. Eng. Res. Des. 2017, 117, 318-335. (17) Foucher, E.R.; Doherty, M.F.; Malone, M.F. Automatic screening of entrainers for homogeneous azeotropic distillation. Ind. Eng. Chem. Res. 1991, 30, 760-772. (18) Gmehling, J.; Möllmann, C. Synthesis of distillation process using thermodynamic models and the Dortmund Data Bank. Ind. Eng. Chem. Res. 1998, 37, 3112-3123. (19) Dejanović, I.; Matijašević, L.; Olujić, Ž. Dividing wall column—A breakthrough towards sustainable distilling. Chem. Eng. Process. 2010, 49, 559-580. (20) Kiss, A.A. Advanced distillation technologies design, control and applications; John Wiley & Sons: West Sussex, 2013. (21) Bravo-Bravo, C.; Segovia-Hernández, J.G.; Gutiérrez-Antonio, C.; Durán, A.L.; Bonilla-Petriciolet, A.; Briones- Ramírez, A. Extractive dividing wall column: design and 28

ACS Paragon Plus Environment

Page 28 of 33

Page 29 of 33 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

optimization. Ind. Eng. Chem. Res. 2010, 49, 3672-3688. (22) Gutiérrez-Guerra, R.; Segovia-Hernández, J.G.; Hernández, S. Reducing energy consumption and CO2 emissions in extractive distillation. Chem. Eng. Res. Des. 2009, 87, 145-152. (23) Kiss, A.A.; Ignat, R.M. Innovative single step bioethanol dehydration in an extractive dividing-wall column. Sep. Purif. Technol. 2012, 98, 290-297. (24) Kiss, A.A.; Suszwalak, D.J-.P.C. Enhanced bioethanol dehydration by extractive and azeotropic distillation in dividing-wall columns. Sep. Purif. Technol. 2012, 86, 70-78. (25) Timoshenko, A.V.; Anokhina, E.A.; Morgunov, A.V.; Rudakov, D.G. Application of the partially thermally coupled distillation flowsheets for the extractive distillation of ternary azeotropic mixtures. Chem. Eng. Res. Des. 2015, 104, 139-155. (26) Anokhina, E.; Timoshenko, A. Criterion of the energy effectiveness of extractive distillation in the partially thermally coupled columns. Chem. Eng. Res. Des. 2015, 99, 165-175. (27) Tututi-Avila, S.; Medina-Herrera, N.; Hahn, J.; Jiménez-Gutiérrez, A. Design of an energy efficient side-stream extractive distillation system. Comput. Chem. Eng. 2017, 102, 17-25. (28) Tututi-Avila, S.; Jiménez-Gutiérrez, A.; Hahn, J. Control analysis of an extractive dividing-wall column used for ethanol dehydration. Chem. Eng. Process. 2014, 82, 88-100. (29) Luyben, W.L. Control of comparison of conventional and thermally coupled ternary extractive distillation processes. Chem. Eng. Res. Des. 2016, 106, 253-262. (30) Wang, X.; Xie, L.; Tian, P.; Tian, G. Design and control of extractive dividing wall column and pressure-swing distillation for separating azeotropic mixtures of acetonitrile/N-propanol. Chem. Eng. Process. 2016, 110, 172-187. (31) Li, L.; Guo. L.; Tu, Y.; Yu, N.; Sun. L.; Tian. Y.; Li, Q. Comparison of different extractive distillation processes for 2-methoxyethanol/toluene separation: Design and control. Comput. Chem. Eng. 2017, 99, 117-134. (32) Modla, G. Energy saving methods for the separation of a minimum boiling point azeotrope using an intermediate entriner. Energy. 2013, 50, 103-109. (33) Zhao, Y.; Ma, K.; Bai, W.; Du, D.; Zhu, Z.; Wang, Y.; Gao, J. Energy-saving thermally coupled ternary extractive distillation process by combining with mixed entrainer for separating ternary mixture containing bioethanol. Energy. 2018, 148, 296-308. 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 33

(34) Zong, L. Study on the separation process of methanol+acetonitrile+benzene system Thesis. Hebei University of Technology. 2013. (35) Zhu, Z.; Xu, D.; Liu, X.; Zhang, Z.; Wang, Y. Separation of acetonitrile/methanol/benzene ternary azeotrope via triple column pressure-swing distillation. Sep. Purif. Technol. 2016, 169, 66-77. (36) Zhu, Z.; Xu, D.; Jia, H.; Zhao, Y.; Wang, Y. Heat integration and control of a triple-column pressure-swing distillation process. Ind. Eng. Chem. Res. 2017, 56, 2150-2167. (37) Luyben, W.L. Control of a triple-column pressure-swing distillation process. Sep. Purif. Technol. 2017, 174, 232-244. (38) Zhu, Z.; Xu, D.; Wang, Y.; Geng, X.; Wang, Y. Effect of multi-recycle streams on the triple-column pressure-swing distillation optimization. Chem. Eng. Res. Des. 2017, 127, 215-222. (39) Maher, P.J.; Smith, B.D. Vapor- liquid equilibrium data for binary systems of chlorobenzene with acetone, acetonitrile, ethyl acetate, ethylbenzene, methanol, and 1-pentene. J. Chem. Eng. Data. 1979, 24, 363- 377. (40)

Ohta,

T.;

Nagata,

I.

Vapor-liquid

equilibria

for

the

ternary

systems

acetonitrile-2-butanone-benzene and acetonitrile-methanol-benzene at 328.15 K. J. Chem. Eng. Data. 1983, 28, 398- 402. (41) Rolemberg, M.P.; Krähenbühl, M.A. Vapor-liquid equilibria of binary and ternary mixtures of benzene, cyclohexane, and chlorobenzene at 40 kPa and 101.3 kPa. J. Chem. Eng. Data. 2001. 46, 256- 260. (42) Luyben, W.L. Distillation Design and Control Using AspenTM Simulation; John Wiley & Sons: Hoboken, New Jersey, 2013. (43) Olujić, Ž.; Sun, L.; de Rijke, A.; Jansens, P.J. Conceptual design of an internally heat integrated propylene-propane splitter. Energy. 2006, 31, 3083-3096. (44) Zhang, Z.; Liu, M.; Zhang, Q.; Li, G.; Gao, J. Simulation and control of DME/methanol/water liquid sidestream distillation column. J Shandong Univ Sci Techno: Nat Sci. 2017, 36, 63-71. (45) Li, Y.; Xu, C. Pressure-swing distillation for separating pressure-insensitive minimum boiling azeotrope methanol/toluene via introducing a light entrainer : design and control. Ind. Eng. Chem. Res. 2017, 56, 4017-4037. 30

ACS Paragon Plus Environment

Page 31 of 33 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

(46) Ghuge, P.D.; Mali, N.A., Joshi, S.S. Comparative analysis of extractive and pressure swing distillation for separation of THF-water separation. Comput. Chem. Eng. 2017, 103, 188-200. (47) Gu, J.; You, X.; Tao, C.; Li, J.; Shen, W.; Li, J. Improved design and optimization for separating tetrahydrofuran-water azeotrope through extractive distillation with and without heat integration by varying pressure. Chem. Eng. Res. Des. 2018, 133, 303-313. (48) Wang, C.; Wang, C.; Guang, C.; Zhang, Z. Comparison of extractive distillation separation sequences for acetonitrile/methanol/benzene multi-azeotropic mixtures. J. Chem. Technol. Biotechnol. 2018, 10.1002/jctb.5693 (49) Luyben, W.L.; Chien, I-L. Design and control of distillation systems for separation azeotropes; John Wiley & Sons: Hoboken, New Jersey, 2010. (50) Luyben, W.L. Comparison of extractive distillation and pressure-swing distillation for acetone-methanol separation. Ind. Eng. Chem. Res. 2008. 47, 2696-2707. (51) Luyben, W.L. Control composition of conventional extractive distillation with a new split-feed configuration. Chem. Eng. Process. 2016. 107, 29-41. (52) Luyben, W.L. Control of heat-integrated extractive distillation processes. Comput. Chem. Eng. 2018, 111,267-277.

ABSTRACT GRAPHIC

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

TABLE OF CONTENTS

ABSTRACT 1. INTRODUCTION 32

ACS Paragon Plus Environment

Page 32 of 33

Page 33 of 33 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

2. DESIGN BASIS 2.1 Thermodynamics 2.2 Economics 3. STEADY-STATE DESIGN 3.1. Design of CED 3.2. Design of extended-EDWC 3.3 Design of extended-SEDC 4. DYNAMIC CONTROL 4.1 Control of CED. 4.2 Control of extended-EDWC 4.3 Control of extended-SEDC 5. CONCLUSION ASSOCIATED CONTENT Supporting information

AUTHOR INFORMATION NOMENCLATURE REFERENCES ABSTRACT GRAPHIC

33

ACS Paragon Plus Environment