Next Article in Journal
Neural-Network-Based Building Energy Consumption Prediction with Training Data Generation
Next Article in Special Issue
Data-Driven Robust Optimization for Steam Systems in Ethylene Plants under Uncertainty
Previous Article in Journal
Model Prediction and Optimization of Waste Lube Oil Treated with Natural Clay
Previous Article in Special Issue
Utility Paths Combination in HEN for Energy Saving and CO2 Emission Reduction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of the Economy and Controllability of Pressure Swing Distillation with Two Energy-Saving Modes for Separating a Binary Azeotrope Containing Lower Alcohols

1
School of Chemistry and Chemical Engineering, Yangtze Normal University, Fuling 408100, China
2
Key Laboratory of Gas Process Engineering, School of Chemistry and Chemical Engineering, Southwest Petroleum University, Chengdu 610500, China
3
China Petroleum Engineering & Construction Corporation Southwest Company, Chengdu 610041, China
*
Author to whom correspondence should be addressed.
Submission received: 30 July 2019 / Revised: 6 October 2019 / Accepted: 8 October 2019 / Published: 12 October 2019
(This article belongs to the Special Issue Process Optimization and Control)

Abstract

:
The pressure swing distillation (PSD) with two different energy-saving modes are put forward to separate a binary azeotrope containing lower alcohols: benzene/methanol. A comparison of the economy and controllability for the partial and fully heat integrated pressure swing distillation (HIPSD) is made by detailed simulation analysis. The optimal operating parameters of partial and fully HIPSD processes are obtained by minimizing total annual cost (TAC). These results show that the fully HIPSD mode saves 5.88% TAC compared with the partial HIPSD mode. Meanwhile, this paper proposes that the composition slope profile can help to select the temperature control stage (TCS), when the temperature profile in the column is rising rapidly near the bottom and the maximum of temperature slope value occurs in the bottom of the column. Several control structures are developed to check the rationality of the selection of the TCS and evaluate the industrial application. These results illustrate the composition/temperature cascade control structure for the PSD with two energy-saving modes can both get good control performances, and the purities of benzene and methanol can be brought close back to the initial value. However, the fully HIPSD mode can only handle much smaller composition disturbances (<10%) compared with the partial HIPSD mode. Therefore, the selection of energy-saving modes for the separation process should weigh economy against controllability.

1. Introduction

Several binary azeotropic systems containing low carbon alcohols, such as benzene and methanol, are widely used in the chemical and petrol chemical industries [1,2,3]. To separate and recycle the azeotropic mixture, several enhanced distillations including pressure swing distillation (PSD) [4,5,6], azeotropic distillation (AD) [7,8], and extractive distillation (ED) [9,10,11] have been applied to separate the azeotrope successfully. For example, Shen et al. put emphasis on the separation of azeotropic systems for the ED method [12]. Shirsat and Wang et al. optimized the ED and PSD method of separating the isobutyl alcohol/isobutyl acetate and tetrahydrofuran/ethanol azeotropic system, and indicated that the cost and product contamination should be concerned when selecting the separation method [13,14]. Compared with AD and ED methods, PSD has drawn more and more attention for the benefit of no third component added into the azeotropic system. Fulgueras and Maloney [15] pointed out that PSD method was economically feasible when the mass fraction of azeotropic composition changed more than 5% (10% is better) in the range of pressure 10 bar. The PSD process has two columns which includes low pressure column (LPC) and high pressure column (HPC). For the minimum boiling azeotropes, the high-purity component can be obtained at the bottom of LPC and HPC. For the maximum boiling azeotropes, the opposite is true.
To reduce heating and cooling loads of thermal utilities and improve heat recovery capacity of the system, some methods such as heat integration are proposed. In general, two modes of heat integration are usually applied in the internal system [6,16,17]. One mode is the heat integration between condenser and reboiler, which includes partial and fully heat integration. The heat integration between the rectifying and stripping section is another mode. In the past years, the economy of PSD processes was compared with different heat integrated modes [1,18,19]. Luyben [1] explored the separation of tetrahydrofuran/water azeotropic system by PSD method with two different heat integrated modes. Wang et al. [19] made a comparison of partial and fully heat integrated pressure swing distillation (HIPSD) modes for the separation of methanol/tetrahydrofuran azeotropic system. These results show that the economy of the PSD process depends on the different azeotropic systems and heat-integrated modes.
Not only the economy but also the controllability of the separation process should be considered in the determination of a suitable separation method. Many researchers pointed out that the economy of the fully HIPSD mode was better but the controllability was weaker than that of the partial HIPSD mode [1,17,19]. The main reason is that the distillate stream of the HPC is used to provide the heat source of the reboiler in the LPC and no auxiliary reboiler is needed, leading to a degree of freedom missed in the full HIPSD process. Luyben [1] studied the energy consumption and controllability of partial and fully HIPSD modes for the separating water/tetrahydrofuran azeotropic system. It was revealed that the energy consumption for the fully HIPSD mode was lowest among the three heat integration modes, but its dynamic controllability was worse compared with the other two heat integration modes. As a result, it is necessary to weigh the energy consumption and controllability for different heat integrated modes. Zhang et al. explored the separation of a methyl acetate/methanol azeotropic system by different heat-integrated modes. It can be seen that the fully HIPSD mode has a little economic benefit compared with partial HIPSD mode, but its controllable performance is also poor [6]. Zhu et al. contrasted the economics and controllability for separating of ethanol and toluene azeotropic system by PSD process with three modes of heat integration [3]. The results indicated that the total annual cost (TAC) for the fully heat integration is lowest, but this kind of heat integration mode can handle only small disturbances very well. Similar studies were done for separating an n-butyl alcohol and toluene system [20], a methanol and tetrahydrofuran system [2], and an ethanol and chloroform system [21].
The selection of the temperature control stage (TCS) is very important to achieve the excellent control performance because the signal of TCS is detected by the temperature controller and transmitted to the adjusted variable. In general, five criteria including product fluctuation minimum criterion, constant temperature criterion, singular value decomposition, sensitivity criterion, and slope criterion are usually used to select the suitable TCS. Zhang et al. used singular value decomposition and sensitivity criterion to select the TCS [6]. Hosgor et al. explored the separation of the chloroform/methanol azeotropic system by fully HIPSD process and used slope criterion to select the TCS [22]. Zhang et al. presented a derivative slope criterion to choose the TCS when the temperature distribution is quite close to linear [23].
So far, there has been no article published that compares PSD with different heat integrated modes for separating benzene/methanol azeotrope through economy and controllability angles. In this work, the feasibility of the PSD process with two energy-saving modes for the separation of the benzene/methanol azeotrope were proved, and a comparison of the economy and controllability for the partial and fully HIPSD process were made in detail. Firstly, the partial and fully HIPSD processes were optimized based on the minimization of the TAC. Secondly, a composition slope criterion was presented to assist the selection of the suitable TCS for the HPC, when the temperature distribution in the column is rising rapidly near the bottom and the maximum of temperature slope value occurs in the bottom of the column. The different degrees of the disturbance in the feed flowrate and composition were introduced to verify the rationality of the selection of TCS and evaluate the industrial application of two processes.

2. Steady State Design and Optimization

2.1. Property Method

The accuracy of the steady-state simulation depends on the selection of the thermodynamic model. Since the benzene/methanol azeotrope is a strongly non-ideal system, several thermodynamic models including Wilson, NRTL and UNIQUAC equations are used for calculating the vapor-liquid equilibrium (VLE) data at 1 bar. The average relative error between calculated and experimental data is listed in the Table 1. As seen from this table, the average relative error of the Wilson equation is minimum. Figure 1 shows the T-x-y diagram of benzene/methanol azeotropic system with experimental data and calculated data at 1 bar. From Figure 1, the calculated data and experimental data are in good agreement. Therefore, the Wilson equation was selected as the suitable property method.

2.2. Optimization Method

TAC including both equipment cost and annual operating cost is a common measure for the economic evaluation in chemical and pharmaceutical process [18]. Equations (1)–(3) show the detailed calculation. Note that the payback period time is set as 3 years and the operation time is set at 8000 h/year. The equipment cost for the partial and fully HIPSD processes mainly includes two distillation columns and heat exchangers. Compared with heat exchangers and distillation columns, the costs for the pumps, reflux drum, pipes and valves are very low, which can even be ignored in the conceptual design phase [18,24]. The diameters of HPC and LPC are computed via the “tray sizing” function in the Aspen plus. The annual operating cost mainly includes the annual steam cost and cooling water cost. The related calculating formulae and parameters are listed in the Table 2 according to previous work of Luyben and Douglas [19,24].
TAC   =   total   equipment   cos t payback   period + annual   operating   cos t   ,
Total equipment cost = distillation columns + heat exchangers cost,
Annual operating cost = annual steam cost + annual cooling water cost.

2.3. Pressure Determination

Figure 2 displays the azeotropic composition and temperature in the pressure range of 0.1–10 bar. From Figure 2, the variation of the azeotropic composition is about 50% when the pressure changes from 0.1 to 10 bar which sufficiently indicates that the PSD method for separating benzene/methanol azeotropic system is economic and feasible.
Furthermore, the key of PSD is the selection of the pressures of two columns. The difference of azeotropic composition increases with the increase of pressure difference, which leads to decrease of the material circulation and energy consumption. However, it increases the cost of the equipment. For the LPC, the complexity and cost of the refrigeration system increase when the operating pressure is too low. In general, the pressure of the LPC is set as 1 bar, and the top vapor stream (about 57 °C) of the LPC can be condensed by the cooling water (32 °C). For the HPC, the selection principle of the pressure is that the top stream can be used to heat the bottom liquid stream of LPC [6,25]. Table 3 gives the TAC of several cases with different pressures in HPC. Note that the temperature difference between the bottom liquid stream (about 67 °C at 1 bar) and the top vapor stream (about 116 °C at 6 bar) is more than 25 °C. Thus, the heat integration can be implemented in the internal system. From Table 3, the TAC is minimum for the HPC at 6 bar while the pressure of HPC at range of 7–10 bar. The main reason is that the equipment cost and operating cost have a significant increase when the operation pressure of HPC is greater than or equal to 7 bar. Therefore, the optimum operation pressure of HPC is selected as 6 bar.

2.4. Optimization

2.4.1. Optimization of Distillation Sequence

For the partial and fully HIPSD processes, the distillation sequence including LP-HP and HP-LP should be first determined via the minimization of the TAC. For the LP-HP distillation sequence, the raw material which contains 70 wt% methanol and 30 wt% benzene is first fed into the LPC. The methanol product is obtained in the bottom of LPC and the azeotropic mixture of benzene and methanol is fed to the HPC. The bottom product of the HPC is the high-purity benzene. The top azeotrope in HPC is used as the heat sources of LPC reboiler and then recycled to the LPC. For the HP-LP distillation sequence, the first column is HPC and the second column is LPC. When the raw material with same composition is fed in the HPC, the high-purity methanol product is first obtained at the bottom of the first column and the high-purity benzene is produced at the bottom of the second column. Table 4 lists the economic comparison for partial and fully HIPSD processes with different distillation sequences. The results indicate that the former one is more economical than the latter one for the partial and fully HIPSD processes.

2.4.2. Optimization of Operating Parameters

For the partial HIPSD mode, the extra energy needs to be supplied (removed) by an auxiliary reboiler (condenser) since the heat duty of reboiler for LPC (QR1) is not equal to the heat duty of condenser for HPC (QC2). When the material composition and the operating pressures are fixed, there are still several variables including the number of stages in LPC and HPC (NT1 and NT2), and the feed positions (NF1, NF2, NR) and the reflux ratios of the LPC and HPC (RR1 and RR2) need to be optimized via set the minimum TAC as the objective function. The purities of the benzene and methanol as the constraint variables are set as 99.9 wt% via the “design specification/vary” function whose manipulating variables is the distillate to feed ratios of two columns. Figure 3 shows the sequential iterative optimization procedure for partial HIPSD process. From this figure, RR1 and RR2 are the innermost loops, NF1, NF2, and NR are the secondary loops, and NT1 and NT2 are the outermost loops. The optimized operating parameters are acquired through multiple iterations of the internal cycle to external cycle parameter.
Figure 4 shows the effect of the seven variables on TAC in the partial HIPSD process. As can be seen in these figures, the effects of each variable on TAC was first decrease then increase within limits. The detailed operating parameters and equipment parameters for the partial HIPSD process is displayed in Figure 5. From this figure, about 331.283 kW energy to the base of the LPC is supplied by the top vapor stream of the HPC and the extra 39.871 kW energy is supplied by the low-pressure steam to remain the energy requirement. The TAC of the partial HIPSD process is $198,891. Compared with the no heat integrated PSD process, the equipment cost and annual energy cost of this process have saved $48,503 and $77,607, respectively.
Compared with the partial HIPSD mode, there is no auxiliary reboiler needed in the fully HIPSD mode. The primary reason is that the top vapor stream of the HPC just can provide the required energy of the reboiler of LPC. Only one heat exchanger serves as the condenser and reboiler of two columns. In order to make the −QC1 equal to QR2, the RR2 is varied to implement the fully heat integration. Meanwhile, the distillate to feed ratios of HPC and LPC are adjusted to maintain the specific purities. The number of the variables in this process is less than the partial HIPSD process. The six variables including NT1, NT2, NF1, NF2, NR, and RR1 should be optimized to obtain the best configuration by minimizing TAC. The optimization procedure and detailed optimizing flowsheet for the fully HIPSD mode are displayed in the Figure 6 and Figure 7. The results indicate that the fully HIPSD process with a TAC of $187,189 has more economic benefit than the partial HIPSD process.

3. Control Structure Design

3.1. Selection of TCS

The selection of TCS is very important to keep the product purity at the specific value. Figure 8 shows the temperature distribution and the temperature slope distribution for the LPC and HPC. From the temperature distributions of two columns, the temperature increases quickly in the rectification section of two columns because the benzene and methanol build up near the bottom. According to the slope criterion by Luyben [24], there is one peak at stage 13 in the temperature slope distribution for LPC, and stage 13 is selected as the suitable TCS. For the HPC, the maximum value of the temperature slope is stage 21 which is the bottom of this column. In theory, the product purity can be achieved by keeping the temperature constant in the top or bottom of a distillation column. However, in reality, if the 21st stage is chosen as TCS, it will be very difficult to maintain the temperature as a constant value. The main reason is that the product stream of benzene is removed from the bottom of HPC where the temperature changes very quickly.
The principle of selecting TCS based on slope criterion is that the temperature difference between two adjacent plates and the composition variation of key component correspond, one by one [24]. Therefore, this paper proposes that the composition and its slope profile can be used to assist the selection for the TCS. Figure 9 displays the composition distribution and its slope distribution of the HPC. From this figure, only one fairly sharp composition break appears on stage 20. Meanwhile, as can be seen in Figure 8b, the temperature slope value of stage 20 is the maxima, when removing the value of the temperature slope on stage 21. In summary, stage 20 in the HPC is selected as the temperature sensitive plate. In the next section, the rationality of selecting temperature sensitive plate will be confirmed by introducing the different degree disturbances with flowrate and composition in the control structures.

3.2. Control Strategies for Partial HIPSD Process

3.2.1. Control Structure CS1

The dynamic characteristics of the partial HIPSD process was studied after optimizing the steady-state process. Under the heuristic method [24], the size of the reflux tank and bottom sump should be computed and entered into the “Dynamics/ Reflux drum and sump”. The volumes of the reflux tank and bottom sump are calculated by Equation (4), the qv is the volume flow which is obtained by the “Profile/Hydraulics/volume flow liquid from stage 1”, and the time of liquid holdup is 10 min to make them 100% full. In general, the ratio of length (or height) to diameter is set as 2. Table 5 lists the size of devices in partial HIPSD process. Besides, the pressure drop of the control valve is set to 3 bar when the valve opening is 50%. Note that the outlet pressure of all control valves must be equal to the pressure of the feed position.
V = 10 × 60 × qv.
The basic control structure of the partial HIPSD process is built after checking the pressure of the steady-state flowsheet. The corresponding flowsheet, control panels, and the flowsheet equation are shown in the Figure 10. It can be seen that four kinds of controller including flow controller, pressure controller, level controller, and temperature controller are used in the flowsheet. For the LPC, the feed flowrate is controlled by flow controller. The levels of the reflux rank and the sump are kept by adjusting the stream flowrate of the top and sump. The pressure is controlled by adjusting the heat removal rate of the condenser. The temperature of the 13th stage is adjusted by the heat duty of the auxiliary reboiler in the LPC. However, for the HPC, the pressure controller has only a display function because the pressure is floating as the different operating condition, and the temperature controller TC2 is used to control the temperature of the 20th stage by adjusting the heat duty of the reboiler. In addition, two lag time modules with 1 min are embedded in the temperature control loops. According to the Tyreus–Luyben tuning by Luyben [24], the tuning parameters of TC1 and TC2 are determined by the relay-feedback rules. The relevant tuning results are given in Table 6.
In order to implement the partial HIPSD, some flowsheet equations used for calculating the heat duty of the composable condenser/reboiler (QC/R) and QAux are shown in Figure 10c. The QC/R is calculated by the first equation which involves the heat transfer area (AC/R), total heat transfer coefficient (KC/R), and the temperature difference between condenser of HPC and reboiler of LPC (∆tC/R). According to the steady-state simulation result, the AC/R calculated is 11.803 m2 and the ∆tC/R is 49.48 °C. The total KC/R is assumed as 0.0020448 GJ/(h m2 K). The QC2 and QR1 are 331.283 and 371.154 kW, respectively. Therefore, the QAux calculated by another equation is 39.871 kW.
The dynamic performance of control structure CS1 for the partial HIPSD process was tested by introducing the feed flowrate and composition disturbances with different degrees at t = 2 h. Figure 11 displays the response of this basic control structure to 5%, 10%, and 20% step changes at the specified value of the FC controller. The black, green, and red solid lines are shown the positive disturbances, while the corresponding dashed lines show the negative disturbances. From the stability of the temperatures on two sensitive trays, when the different degrees of disturbances are introduced in this system, the temperatures can be brought back to the set values by adjusting QAux and QR2 in about 2 h. From the view of the transient deviation, the larger the disturbance in the system, the larger the transient deviation in the product purities. When the disturbances are not more than 10%, the two product purities return close to their set value in about 2 h. However, the purity of benzene in the HPC (B2) is away from the desired purity when the system encounters +20% feed flowrate disturbance (99.2 mol% for the new steady state) and the −20% composition disturbance (97.2 mol% for the new steady state).
The results illustrate that the selection of the temperature sensitive plate in the HPC is reasonable. However, the benzene purity cannot turn back to the specified purity, which is mainly because the pressure of the HPC is not fixed as the operating condition. An improved control structure should be developed to overcome this problem.

3.2.2. Control Structure CS2

As shown in Figure 11, a large transient deviation appeared in the benzene product purity and a long time is needed to restore a new steady state when this system faces a great disturbance. The pressure-compensated temperature control structure, which has been described in detail and confirmed in some literature [1,17,25], are applied to improve this problem according to the previous control scheme.
The first step to build this control scheme is to compute the VLE data of the liquid mixture on stage 20 according to the steady-state results. Equation (5), which is obtained by the bubble point temperature–pressure diagram (its slope is 5.94765), is used for calculating the corresponding pressure-compensated temperature TPC, where the temperature T is for the temperature on stage 20 in the HPC, the pressure P is for the pressure in the top of the HPC. To implement this control scheme in Aspen Dynamic, the third equation in the “flowsheet equation” is used, which is shown in Figure 12. The TPC calculated is the input signal of the deadtime block “dead2” and then is fed to TC2. The improved control structure is shown in Figure 12.
T P C = T ( P 6 ) × 5.94765 .
The effectiveness of this improved control structure is estimated in Figure 13 for three different degrees of disturbance in feed flowrate and composition. It is noteworthy that the degree of the disturbance and time introduced are the same as the previous control scheme. From Figure 13a, the pressure in the HPC are floated in the range of 5–7 bar when this system faces 5%, 10%, and 20% disturbances. The largest deviation of purity in the bottom stream of HPC reduces to −0.0005 wt% and 0.0005 wt% instead of −0.007 wt% and 0.0008 wt% after 20% and −20% variance in the feed flowrate. The improvement in the methanol product purity is not obvious when the six kinds of disturbance take place in the time equal to 2 h.
Figure 13b demonstrates the responses for 5%, 10%, and 20% composition disturbances. As can be seen in this figure, the pressure in the HPC is also not fixed, but it changes with the disturbances. From the benzene purity in the bottom stream B2, the purity of benzene shows the value of 99.97 wt% and 99.83 wt% for 20% and −20% disturbances in the composition. Compared with basic control structure, the transient deviation decreases from 2.8 to 0.17 wt%. However, the methanol purity in the six cases is similar to those found in the previous control structure.
These results illustrate that the pressure-compensated temperature control structure has a significant improvement in the transient deviation in the xB2 (benzene), but the effect in improving the purity of methanol is not obvious when the same disturbances occur at 2 h.

3.2.3. Control Structure CS3

The purities of the methanol and benzene cannot be driven to the initial setpoint when the system encounters the disturbances with different degrees in the feed flowrate and composition even though the control structure CS2 can obtain a robust control. A composition/temperature cascade control structure, which integrates the advantages of both composition control and temperature control, was put forward to improve the control structure.
To set up the composition/temperature control structure, two composition controllers (CC1 and CC2) with 3 min deadtime, which detect the content of methanol and benzene in the product streams B1 and B2, were inserted in the loops of TC1 and TC2. In the composition/temperature cascade loop, the composition and temperature controllers are the major and assistant controllers, respectively. It is interesting to note that the operating variables of the temperature controllers are the same as the previous control structure, but the input signal which is the output signal of the composition controller is not a constant. The tuning parameters of the major controllers CC1 and CC2 need to be determined by using the Tyreus–Luyben tuning rules and rerunning the relay-feedback tests. The ultimate gain and period of for CC1 are 14.96 and 25.08 min, and for CC2 are 15.40 and 18.48 min. Figure 14 shows control structure CS3 and control-panels.
Figure 15 displays the dynamic responses of the control structure encountering the different degrees of disturbance. The solid lines with black, green, and red colors signify the feed flowrate and the methanol composition increases of 5%, 10%, and 20%, and the corresponding dashed lines signify the feed flowrate and the methanol composition decreases of 5%, 10%, and 20% in Figure 15a,b. As can be seen from Figure 15, the temperatures on the two sensitive plates are not maintained at the initial temperatures of 63.31 and 123.29 °C, but change with the different degrees of disturbance. The larger the variation in the flowrate or composition, the larger the deviation in the initial value. The purities of two products have small transient deviation and can be brought back to the setting point when this control structure faces the ±5%, ±10%, and ±20% disturbances in the feed flowrate and composition. The results illustrate that the control structure CS3 can provide a good dynamic control performance.

3.3. Control Strategies of Fully HIPSD Process

3.3.1. Control Structure CS1

As mentioned earlier, the fully HIPSD process for the separation of the benzene and methanol system has significant economic superiority. Therefore, the dynamic performance of this process was also explored by the Aspen Dynamic. According to the temperature slope and composition slope criterions, the stages 13 and 20 in the LPC and HPC were chosen as the TCS. Figure 16 demonstrates the first control structure, control-panel, and flowsheet equations for the fully HIPSD process. From Figure 16a, the kind and number of the controller are similar to the partial HIPSD process. All of the controllers are proportional control, and the flowrate controller is set as the conventional tuning parameters with τI of 0.3 min and KC of 0.5. Four level controllers have the τI of 9999 min and KC of 2. Note that the temperature T13th is controlled by changing RR1. The corresponding τI and KC are calculated by running the relay-feedback rules. The gain and integral time of TC1 and TC2 are 9.97 and 3.39 min, 17.16 and 5.28 min, respectively.
The top of the HPC and bottom of the LPC share only one heat exchanger, the heat duty of the heat exchanger is equal to the product of KC/R (assumption for 0.0020448 GJ/ (h m2 K)), AC/R, and ∆tC/R. The “flowsheet equation” function in the Aspen dynamic aids in implementing the fully heat integration. The QR1 and QC2 are respectively calculated by the first and second equations. The performance of the control structure will be estimated in the next paragraph.
The disturbances with different degrees were utilized to test the control performance for the fully HIPSD process. Figure 17a gives the dynamic responses when 5%, 10%, 20% step change in the feed flowrate occur at the time equal to 2 h. The temperatures on two TCS are maintained to the set point (63.5 and 124.8 °C) in approximately 2 h after the flowrate disturbances occurred at 2 h. This is caused by the controllers “TC1”and “TC2” that can readjust the RR1 and QR2 for achieving a new steady-state condition. The product purities of five cases (±5%, ±10%, and −20% flowrate disturbances) are held quite close to the specified purity, but the transient deviation in the purity of benzene is 0.9 wt% when the flowrate changes from 1000 to 1200 kg/h. The main culprit is that the increase of the feed flowrate leads to a sharp increase in the pressure of HPC and results in the concentration of methanol on the sensitive stage, but the temperature on stage 20 is controlled at a constant.
Figure 17b illustrates the dynamic responses of the control structure for the feed composition disturbance with three different degrees. The solid and dashed lines are for positive and negative change, respectively. The black, green, and red lines demonstrate the 3%, 5%, and 10% step change in the feed composition. The results of the response for the six cases are similar to the six cases in Figure 17a. The larger the variance in the composition, the larger the deviation in the benzene purity. The deviation (0.8 wt%) occurred when the methanol concentration decreases from 70 to 63 wt%.
When the methanol concentration changes more than 7%, the flowrate of the recycling stream increases to the point which caused the liquid level of the reflux tank in the LPC to exceed 100% and the system was forced to shut down. This is caused by more benzene in the feed stream making the product in the bottom of HPC increase and the pressure rise rapidly.
These results illustrate that the selection of the temperature sensitive plate is suitable, and this control structure can implement the robust control. However, it cannot maintain the purity of the product at the specification after encountering a large disturbance in the composition. Therefore, it is necessarily to develop several improved structures to overcome this problem.

3.3.2. Control Structure CS2

As described above, the high purity benzene product cannot be obtained when the large disturbance is introduced in the previous control structure. The pressure-compensated temperature can improve the problem that the changing pressure with operating condition results in the large transient deviation in the benzene product in the HPC.
The VLE data also needs to be recalculated before establishing the new control structure. The major reason is that the operating parameters for the fully HIPSD process are different from that for the partial HIPSD process which will result in the difference in the composition on stage 20 in the HPC. Equation (6) shows the relationship of the temperature on stage 20 in the HPC, the pressure in the top of the HPC, and the pressure-compensated temperature TPC. In Aspen Dynamic, TPC can be calculated by the third equation in the “flowsheet equation” function. As shown in Figure 18, the input signal of the “dead 2” block is the TPC. Note that the output signal of the temperature controller TC2, which is the ratio of reboiler heat duty 1.49297 GJ/h to the flowrate 1000 kg/h, is 0.0015 and not connected the reboiler, but the multiplier block. The effectiveness of this improved structure will be illustrated in the next section.
T P C = T ( P 6 ) × 7.5722 .
The disturbances in the flowrate and composition with several different degrees were used to test the dynamic performance. The obtained results are demonstrated in Figure 19. Figure 19a shows the response of the flowrate with 5%, 10%, and 20% step change, and Figure 19b shows the response of the methanol component with 3%, 6%, and 10% step change.
From the variation of the pressure in the HPC, the pressure is floating from 5 to 7 bar when the system encounters the 12 kinds of disturbances. From the benzene purity in the bottom stream B2, the product purity can return to the specified value after facing the 5%, 10%, and 20% flowrate disturbances. Meanwhile, the product purity also can be held close to the specification when the methanol composition in the feed is varied 3%, 6%, and 10%. The most purity deviation decreases more sharply than control structure CS1. The time to reach a new steady state was reduced from ~3 to 1 h after experiencing the same disturbance in the flowrate and composition. However, the methanol product purity and the transient deviation in the purity are similar to the response for the control structure CS1.
These results illustrate that the pressure-compensated temperature control structure can improve the purity of the benzene, decrease the transient deviation in the purity and the time to reach a new steady state, but play a small role in the increase of methanol purity.

3.3.3. Control Structure CS3

To keep the methanol purity at the set point, the third control structure which is the composition/temperature cascade control structure was set up to implement the aim.
As for Figure 20, the composition controller CC1 with 3 min is added to the temperature control loop. The purity of the methanol product is detected and input in the CC1 composition controller, and the output signal is the setpoint of the TC1 temperature controller. Note that the setpoint is not fixed, but changes with the disturbances in the flowrate and composition. The CC1 composition controller, which is cascade with TC1, is tuned by running the relay-feedback tests and using the Tyreus–Luyben tuning rules. The corresponding tuning parameters of CC1 are KC = 35.82, t = 31.68 min.
In the next part, we will test the dynamic performances by introducing the same disturbances in flowrate and composition at 2 h to verify the effectiveness of the improved structure.
The responses for the different degree of disturbances are displayed in Figure 21. The different color lines stand for the same variation in the feed flowrate and composition in the basic control structure. From Figure 21, the product purity of methanol can be brought back to the specified value after encountering feed flowrate changes of 50, 100, and 200 kg/h and the composition changes of 2.1, 4.2, and 7 wt% methanol. The larger the disturbances in the flowrate and composition, the larger the transient deviation in the initial purity. However, the transient deviation is smaller than those in the previous structure control. Meanwhile, the time to reach a new steady state reduces from 4 to 1 h.
Compared with previous control structure, the benzene purity is almost unchanged when the system faces the disturbances in the feed flowrate and composition. The methanol purity deviation in the new steady state is maintained within acceptable range when there is a 3%, 6%, and 10% step change in the methanol composition. This reveals that the control structure CS3 can provide a good dynamic control performance.

4. Conclusions

The PSD with two energy-saving modes for the separation of the benzene/methanol azeotropic system was studied by Aspen Plus and Aspen Dynamics. We contrasted the economy and controllability of the PSD with two heat integration modes. The optimization results show that the fully HIPSD process with a TAC of $187,189 has some economic benefit compared with the partial HIPSD process with a TAC of $198,891.
The temperature profile in the HPC rose rapidly near the bottom of this column and the maximum of temperature slope value occurred in the last stage, which is the bottom of the HPC. The composition slope profile was proposed to assist with the selection of the temperature sensitive plate. According to the two criterions, stages 13 and stage 20 in the LPC and HPC were chosen as the TCS. Several control structures were developed to evaluate the industrial application and check the rationality of the selection of the temperature sensitive plate. The different degrees of disturbances in the feed flowrate and composition were introduced in two processes at 2 h. The results illustrate that the selection of the TCS is reasonable and the composition/temperature control structure can provide rust control performance for the PSD with two heat integrated modes. Compared with the partial HIPSD mode, the fully HIPSD mode can handle much smaller composition disturbance despite little economic benefit.
These works can provide a reference for selecting TCS, especially for azeotropic systems whose maximum value of temperature slope profile occurs in the bottom of the column. Meanwhile, the findings can also help engineers to weigh the economy against controllability for selection of the suitable separating process.

Author Contributions

Conception and design: Z.Z., L.L.; Process simulation and optimization: L.L., Z.Z.; Establishment of control structure: H.L., H.H.; Analysis and interpretation: L.L., H.H.; Data collection: H.L., H.H.; Writing the article: L.L.; Critical revision of the article: L.L., Z.Z., H.L., H.H.; Final approval of the article: L.L., Z.Z., H.L., H.H.

Funding

This research was funded by the project of science and technology research program of Chongqing Education Commission of China (No. KJQN201901429).

Acknowledgments

The authors gratefully acknowledge the anonymous reviewers for the insightful comments and suggestions, which improved the quality of this paper to a large extent. Additionally, the authors also gratefully acknowledge the financial support by the project of science and technology research program of Chongqing Education Commission of China (No. KJQN201901429).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Luyben, W.L. Design and Control of a Fully Heat-Integrated Pressure-Swing Azeotropic Distillation System. Ind. Eng. Chem. Res. 2008, 47, 2681–2695. [Google Scholar] [CrossRef]
  2. Wang, Y.; Zhang, Z.; Zhang, H.; Zhang, Q. Control of Heat Integrated Pressure-Swing-Distillation Process for Separating Azeotropic Mixture of Tetrahydrofuran and Methanol. Ind. Eng. Chem. Res. 2015, 54, 1646–1655. [Google Scholar] [CrossRef]
  3. Zhu, Z.; Wang, L.; Ma, Y.; Wang, W.; Wang, Y. Separating an Azeotropic Mixture of Toluene and Ethanol Via Heat Integration Pressure Swing Distillation. Comput. Chem. Eng. 2015, 76, 137–149. [Google Scholar] [CrossRef]
  4. Lv, L.Z.L.; Liu, H.; Li, H.; Sun, S. Comparison of Continuous Homogenous Azeotropic and Pressure-Swing Distillation for a Minimum Azeotropic System Ethyl Acetate N-Hexane Separation. Chin. J. Chem. Eng. 2018, 26, 2023–2033. [Google Scholar]
  5. Li, X.; Zhao, Y.; Qin, B.; Zhang, X.; Wang, Y.; Zhu, Z. Optimization of Pressure-Swing Batch Distillation with and without Heat Integration for Separating Dichloromethane/Methanol Azeotrope Based on Minimum Total Annual Cost. Ind. Eng. Chem. Res. 2017, 56, 4104–4112. [Google Scholar] [CrossRef]
  6. Zhang, Z.; Zhang, Q.; Li, G.; Liu, M.; Gao, J. Design and Control of Methyl Acetate-Methanol Separation Via Heat-Integrated Pressure-Swing Distillation. Chin. J. Chem. Eng. 2016, 24, 1584–1599. [Google Scholar] [CrossRef]
  7. Zhao, L.; Lyu, X.; Wang, W.; Shan, J.; Qiu, T. Comparison of Heterogeneous Azeotropic Distillation and Extractive Distillation Methods for Ternary Azeotrope Ethanol/Toluene/Water Separation. Comput. Chem. Eng. 2017, 100, 27–37. [Google Scholar] [CrossRef]
  8. Li, W.; Zhong, L.; He, Y.; Meng, J.; Yao, F.; Guo, Y.; Xu, C. Multiple Steady-States Analysis and Unstable Operating Point Stabilization in Homogeneous Azeotropic Distillation with Intermediate Entrainer. Ind. Eng. Chem. Res. 2015, 54, 7668–7686. [Google Scholar] [CrossRef]
  9. Lladosa, E.; Monton, J.B.; Burguet, M. Separation of di-N-Propyl ether and n-Propyl Alcohol by Extractive Distillation and Pressure-Swing Distillation: Computer Simulation and Economic Optimization. Chem. Eng. Process. Process Intensif. 2011, 50, 1266–1274. [Google Scholar] [CrossRef]
  10. Shen, W.; Benyounes, H.; Gerbaud, V. Extension of Thermodynamic Insights on Batch Extractive Distillation to Continuous Operation. 1. Azeotropic Mixtures with a Heavy Entrainer. Ind. Eng. Chem. Res. 2013, 52, 4606–4622. [Google Scholar] [CrossRef] [Green Version]
  11. Shen, W.; Dong, L.; Wei, S.A.; Li, J.; Benyounes, H.; You, X.; Gerbaud, V. Systematic Design of an Extractive Distillation for Maximum-Boiling Azeotropes with Heavy Entrainers. Aiche J. 2015, 61, 3898–3910. [Google Scholar] [CrossRef]
  12. Shen, W.; Benyounes, H.; Gerbaud, V. Extractive Distillation: Recent Advances in Operation Strategies. Rev. Chem. Eng. 2015, 31, 13–26. [Google Scholar] [CrossRef]
  13. Shirsat, S.P. Separation of Isobutyl Alcohol and Isobutyl Acetate by Extractive Distillation and Pressure-Swing Distillation: Simulation and Optimization. Sep. Purif. Technol. 2015, 154, 366. [Google Scholar] [CrossRef]
  14. Wang, Y.; Cui, P.; Ma, Y.; Zhang, Z. Extractive Distillation and Pressure-Swing Distillation for THF/Ethanol Separation. J. Chem. Technol. Biotechnol. 2015, 90, 1463–1472. [Google Scholar] [CrossRef]
  15. Fulgueras, A.M.; Poudel, J.; Kim, D.S.; Cho, J. Optimization Study of Pressure-Swing Distillation for the Separation Process of a Maximum-Boiling Azeotropic System of Water-Ethylenediamine. Korean J. Chem. Eng. 2015, 33, 46–56. [Google Scholar] [CrossRef]
  16. Suo, X.; Ye, Q.; Li, R.; Dai, X.; Yu, H. The Partial Heat-Integrated Pressure-Swing Reactive Distillation Process for Transesterification of Methyl Acetate with Isopropanol. Chem Eng. Process. Process Intensif. 2016, 107, 42–57. [Google Scholar] [CrossRef]
  17. Luo, B.; Feng, H.; Sun, D.; Zhong, X. Control of Fully Heat-Integrated Pressure Swing Distillation for Separating Isobutyl Alcohol and Isobutyl Acetate. Chem. Eng. Process. Process Intensif. 2016, 110, 9–20. [Google Scholar] [CrossRef]
  18. Luyben, W.L. Methanol/Trimethoxysilane Azeotrope Separation Using Pressure-Swing Distillation. Ind. Eng. Chem. Res. 2014, 53, 5590–5597. [Google Scholar] [CrossRef]
  19. Wang, Y.; Cui, P.; Zhang, Z. Heat-Integrated Pressure-Swing-Distillation Process for Separation of Tetrahydrofuran/Methanol with Different Feed Compositions. Ind. Eng. Chem. Res. 2014, 53, 7186–7194. [Google Scholar] [CrossRef]
  20. Zhang, H. Heat Integration Process for Separating System of Toluene and N-Butyl Alcohol with Pressure Swing Distillation; University of Science and Technology: Qingdao, China, 2015. [Google Scholar]
  21. Zhang, Q. Dynamics and Control of Pressure Swing Distillation Process for Azeotropic Systems of Methanol/Ethanol-Tetrahydrofuran Qingdao; University of Science and Technology: Qingdao, China, 2015. [Google Scholar]
  22. Hosgor, E.; Kucuk, T.; Oksal, I.N.; Kaymak, D.B. Design and Control of Distillation Processes for Methanol–Chloroform Separation. Comput. Chem. Eng. 2014, 67, 166–177. [Google Scholar] [CrossRef]
  23. Zhang, Q.; Peng, J.; Zhang, K. Separation of an Azeotropic Mixture of Dimethyl Carbonate and Methanol Via Partial Heat Integration Pressure Swing Distillation. Asia Pac. J. Chem. Eng. 2017, 12, 50–64. [Google Scholar] [CrossRef]
  24. Luyben, W.L. Design and Control of Distillation Systems for Separating Azeotropes; John Wiley & Son: Hoboken, NJ, USA, 2008. [Google Scholar]
  25. Cao, Y.; Li, M.; Wang, Y.; Zhao, T.; Li, X.; Zhu, Z.; Wang, Y. Effect of Feed Temperature on Economics and Controllability of Pressure-Swing Distillation for Separating Binary Azeotrope. Chem. Eng. Process. Process Intensif. 2016, 110, 160–171. [Google Scholar] [CrossRef]
Figure 1. T-x-y diagram of benzene/methanol system with experimental data and calculated data at 1 bar.
Figure 1. T-x-y diagram of benzene/methanol system with experimental data and calculated data at 1 bar.
Processes 07 00730 g001
Figure 2. Effect of pressure on azeotropic composition and temperature.
Figure 2. Effect of pressure on azeotropic composition and temperature.
Processes 07 00730 g002
Figure 3. Sequential iterative optimization procedure for the partial HIPSD process.
Figure 3. Sequential iterative optimization procedure for the partial HIPSD process.
Processes 07 00730 g003
Figure 4. The effect of the seven variables on total annual cost (TAC) in the partial HIPSD process: (a) RR1, (b) RR2, (c) NF1, (d) NF2, (e) NR, (f) NT1, (g) NT2.
Figure 4. The effect of the seven variables on total annual cost (TAC) in the partial HIPSD process: (a) RR1, (b) RR2, (c) NF1, (d) NF2, (e) NR, (f) NT1, (g) NT2.
Processes 07 00730 g004
Figure 5. Flowsheet of the best optimal partial HIPSD process.
Figure 5. Flowsheet of the best optimal partial HIPSD process.
Processes 07 00730 g005
Figure 6. Sequential iterative optimization procedure for the fully HIPSD process.
Figure 6. Sequential iterative optimization procedure for the fully HIPSD process.
Processes 07 00730 g006
Figure 7. Flowsheet of the best optimal fully HIPSD process.
Figure 7. Flowsheet of the best optimal fully HIPSD process.
Processes 07 00730 g007
Figure 8. Temperature distribution and their slope distribution: (a) high pressure column (HPC), (b) low pressure column (LPC).
Figure 8. Temperature distribution and their slope distribution: (a) high pressure column (HPC), (b) low pressure column (LPC).
Processes 07 00730 g008
Figure 9. Composition distribution and its slope distribution of HPC.
Figure 9. Composition distribution and its slope distribution of HPC.
Processes 07 00730 g009
Figure 10. Control structure CS1, control panels, and flowsheet equations for the partial HIPSD process.
Figure 10. Control structure CS1, control panels, and flowsheet equations for the partial HIPSD process.
Processes 07 00730 g010
Figure 11. The dynamic responses of the basic control structure.
Figure 11. The dynamic responses of the basic control structure.
Processes 07 00730 g011aProcesses 07 00730 g011b
Figure 12. Control structure CS2 and flowsheet equations.
Figure 12. Control structure CS2 and flowsheet equations.
Processes 07 00730 g012
Figure 13. The dynamic responses of control structure CS2.
Figure 13. The dynamic responses of control structure CS2.
Processes 07 00730 g013aProcesses 07 00730 g013b
Figure 14. Control structure CS3 and control-panel.
Figure 14. Control structure CS3 and control-panel.
Processes 07 00730 g014
Figure 15. The dynamic responses of composition/temperature cascade control structure.
Figure 15. The dynamic responses of composition/temperature cascade control structure.
Processes 07 00730 g015aProcesses 07 00730 g015b
Figure 16. Basic control structure, control-panel, and flowsheet equations for the fully HIPSD process.
Figure 16. Basic control structure, control-panel, and flowsheet equations for the fully HIPSD process.
Processes 07 00730 g016
Figure 17. The dynamic responses of the basic control structure.
Figure 17. The dynamic responses of the basic control structure.
Processes 07 00730 g017
Figure 18. Pressure-compensated temperature control structure of flowsheet, control panel, and flowsheet equations.
Figure 18. Pressure-compensated temperature control structure of flowsheet, control panel, and flowsheet equations.
Processes 07 00730 g018
Figure 19. The dynamic responses of pressure-compensated temperature control structure.
Figure 19. The dynamic responses of pressure-compensated temperature control structure.
Processes 07 00730 g019
Figure 20. Control structure CS3 of flowsheet and control panel.
Figure 20. Control structure CS3 of flowsheet and control panel.
Processes 07 00730 g020
Figure 21. The dynamic responses of control structure CS3.
Figure 21. The dynamic responses of control structure CS3.
Processes 07 00730 g021
Table 1. The error analysis of the thermodynamic models.
Table 1. The error analysis of the thermodynamic models.
Thermodynamic ModelsWilsonNRTLUNIQUAC
Average relative error/%0.92952.2471.8563
Table 2. Basis economics and equipment sizing.
Table 2. Basis economics and equipment sizing.
ParametersData
TAC TAC = total   equipment   cos t payback   period + annual   operating   cos t  
Payback period = 3 years
Total equipment costDistillation column costEquipment cost = 17,640 × (D, m)1.066 × (L, m)0.802
Length: H = 1.2 × 0.61 × (NT − 2)
Condensers or Condenser/ReboilerEquipment cost = 7296 × (Ac, m2)0.65
Heat-transfer coefficient (KC) = 0.852 kW/K·m2
Heat transfer area: Ac = Q ( K c × Δ t )
ReboilersEquipment cost = 7296 × (Ac, m2)0.65
Heat-transfer coefficient (KR) = 0.568 kW/K·m2
Heat transfer area: Ac = Q ( K R × Δ t )
Annual operating costAnnual stream costAnnual stream cost = Q × CS × 8000 × 3600
LP steam (6 bar, 160 °C): $7.78/GJ
MP steam (11 bar, 184 °C): $8.22/GJ
HP steam (42 bar, 254 °C): $9.88/GJ
Annual cooling water costAnnual cooling water cost = Q × Cc × 8000 × 3600
Cooling water: $0.354/GJ
Table 3. Case studies of pressure swing distillation (PSD) with partial heat integration.
Table 3. Case studies of pressure swing distillation (PSD) with partial heat integration.
VariablesCase1Case2Case3Case4Case5
P1 (bar)11111
P2 (bar)678910
NT1/NT220/2020/2020/2020/2020/20
NF1/NR/NF210/6/1410/10/1011/7/129/10/1510/10/10
RR1/RR20.7/11.5/21.5/41.3/24/2
ID1/ID20.437/0.3940.539/0.4610.529/0.5380.512/0.4240.623/0.260
Rrec (kg/h)987.597829.923770.849770.094450.316
QC1/QC2−404.579/0−532.573/0−524.103/−215.221−484.754/0−693.918/0
QR1/QR239.144/410.21869.3/510.740/789.31632.111/504.469439.747/307.238
Equipment cost (105 $/y)3.270843.884734.310633.650643.87671
Energy cost (105 $/y)1.048111.474381.922031.315631.78340
TAC (105 $/y)2.140792.769293.358912.532513.07564
Table 4. Economic comparison for partial and fully heat integrated pressure swing distillation (HIPSD) processes with different distillation sequences.
Table 4. Economic comparison for partial and fully heat integrated pressure swing distillation (HIPSD) processes with different distillation sequences.
VariablesPartial HIPSDFully HIPSD
LP-HPHP-LPLP-HPHP-LP
P1 (bar)1616
P2 (bar)6161
NT1/NT220/2020/2019/2120/20
NF1/NR/NF210/6/146/16/1012/9/146/16/10
RR1/RR20.7/12/1.50.61/1.0972/3.768
ID1/ID20.437/0.3940.858/0.6630.313/0.40.773/0.912
Rrec (kg/h)987.5971703.14971.0511664.69
QC1/QC2−404.579/0−662.357/−784.539−376.467/−348.977−1320.5/−1462.78
QR1/QR239.144/410.2181494.72/662.357348.977/421.2761509.54/1320.51
Equipment cost (105 $/y)3.270846.209012.537825.94179
Energy cost (105 $/y)1.048113.496640.982313.09342
TAC (105 $/y)2.140795.566311.828255.07401
Table 5. The size of devices for the partial HIPSD process.
Table 5. The size of devices for the partial HIPSD process.
EquipmentLow Pressure ColumnHigh Pressure Column
Reflux tankVolume flow (m3/s)0.0007229420.00073135
Volume (m3)0.4340.439
Diameter/length (m)0.65/1.300.65/1.31
SumpVolume flow (m3/s)0.0007200350.00138284
Volume (m3)0.4320.830
Diameter/length (m)0.65/1.300.81/1.62
Table 6. Tuning parameters of temperature controllers TC1 and TC2.
Table 6. Tuning parameters of temperature controllers TC1 and TC2.
ParameterTC1TC2
Control variableT1,13T2,20
Operational variableQAR1QR2
Control variable range/°C0–126.620–248.89
Operating variable range/GJ/h0–0.2850–2.91
Final gain/%/%15.3610.02
Final cycle/min3.63.6
Gain/%/%4.803.13
Integration time/min7.927.92

Share and Cite

MDPI and ACS Style

Lv, L.; Li, H.; Zhang, Z.; Huang, H. Comparison of the Economy and Controllability of Pressure Swing Distillation with Two Energy-Saving Modes for Separating a Binary Azeotrope Containing Lower Alcohols. Processes 2019, 7, 730. https://0-doi-org.brum.beds.ac.uk/10.3390/pr7100730

AMA Style

Lv L, Li H, Zhang Z, Huang H. Comparison of the Economy and Controllability of Pressure Swing Distillation with Two Energy-Saving Modes for Separating a Binary Azeotrope Containing Lower Alcohols. Processes. 2019; 7(10):730. https://0-doi-org.brum.beds.ac.uk/10.3390/pr7100730

Chicago/Turabian Style

Lv, Liping, Hang Li, Zheng Zhang, and Huisheng Huang. 2019. "Comparison of the Economy and Controllability of Pressure Swing Distillation with Two Energy-Saving Modes for Separating a Binary Azeotrope Containing Lower Alcohols" Processes 7, no. 10: 730. https://0-doi-org.brum.beds.ac.uk/10.3390/pr7100730

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop