Startseite Optimization of Pressure-Swing Distillation by Evolutionary Techniques: Separation of Ethanol-Water and Acetonitrile-Water Mixtures
Artikel
Lizenziert
Nicht lizenziert Erfordert eine Authentifizierung

Optimization of Pressure-Swing Distillation by Evolutionary Techniques: Separation of Ethanol-Water and Acetonitrile-Water Mixtures

  • Adjay Sagar S und Imran Rahman EMAIL logo
Veröffentlicht/Copyright: 13. Oktober 2017
Veröffentlichen auch Sie bei De Gruyter Brill

Abstract

Complete separation can be achieved in selective homogeneous azeotropic mixtures by exploiting the pressure sensitive nature of the system. In the present work the optimal number of trays, feed location and reflux ratio for sequential column systems encountered in continuous pressure swing distillation (PSD) have been determined by use of two evolutionary techniques. Two industrially relevant systems: ethanol-water and acetonitrile-water have been considered. The Napthali-Sandholm model is solved to obtain the concentration and temperature profiles. The objective is to minimize the total cost using Genetic Algorithm (GA) and Differential Evolution (DE) for the two azeotropic systems. The techniques offer attractive features like applicability to discontinuous and non-differentiable search spaces.

Appendix

A

A.1 Costing and design of the Column

TotalCost=(EnergyCost,C1)+(Depreciation+Interest+Maintenance)×(FixedCost,C2)

The energy costs were calculated according to:

C1=Qr×Cs×ND×24λsteam

Qr Is the reboiler duty, λsteam is the latent heat of evaporation of steam, Cs refers to the steam cost, ND denotes the number of yearly working days.

Cpack=Ac×N×HETP×Cpack0

N is the number of stages, Ac Represents the column area, for area calculations it is necessary to know the velocity.

We choose the velocity to be 80% of the flooding velocity.

uf=K1ρlρvρv
FLV=LwVwρvρl

LW=Liquid mass flow rate, kg/s

VW=Vapour mass flow rate, kg/s

uf=Flooding vapour velocity, m/s

K1 Can be found from K1 v/s FLV diagram in Coulson Richardson Chemical engineering, Vol. 6, Page 585 [24].

Ccol=π×1.4×dc×Nst×HETP×Ws×ρs×Csteel

dc, Ws ,ρs and Csteel. Refer to the column diameter, column thickness, density and cost of steel respectively.

Cvacuumpump=1.5Lakhs

A heuristic was taken to consider the fixed cost of the pump, energy calculations were done considering the power necessity of the pump to be 8 units of electricity and at a nominal cost of Rs. 6/unit of Electricity.

Assuming depreciation, interest, and maintenance costs of 18%, 15% and 2%, respectively.

CT=C1+0.35×C2

A.2 Napthali-Sandholm Method

The Napthali-Sandholm method is one of the most robust techniques and is superior in convergence compared to the successive substitution methods but constraints have to put on the iterated variables for smooth convergence. The simple distillation flow diagram is shown in Figure 14.

Figure 14: Schematic of simple distillation column.
Figure 14:

Schematic of simple distillation column.

Material Balance Equations:

F(1,j)+L2×x(2,j)+V0×x(1,j)L1×x(1,j)V1×y(1,j)=m(1,j)j=1,2....c
F(i,j)+Li+1×x(i+1,j)+Vi1×y(i1,j)Li×x(i,j)Vi×y(i,j)=m(i,j)

                                     j=1,2....ci=1,2....n1

F(n,j)+Ln+1×y(n,j)+Vn1×x(n1,j)Ln×x(n,j)Vn×y(n,j)=m(n,j)

                                     j=1,2....c

Summation Equations:

j=1c(K(i,j)1)×x(i,j)=0(RachfordRiceEquation)

y(i,j)=K(i,j)×x(i,j)i=1,2...nj=1,2...c

K values are predicted using the UNIQUAC model.

The Jacobian of the Method forms a Block-Tridiagonal Matrix of the form:

J=B1C1000A2B2C2000A3000Cn1000An1Bn

Where Ai=M(i)X(i1), X(i1)=[x(i,1),x(i1,2)x(i1,c),T(i1)]

M(i)=[m(i,1),m(i1,2)m(i1,c),m(i1,c+1)]

WhereAi=m(i,1)x(i1,1)m(i,1)x(i1,2)m(i,1)T(i1))m(i,2)x(i1,1)m(i,2)T(i1))m(i,c)x(i1,1)m(i,c)T(i1))

Bi=M(i)X(i)Ci=M(i)X(i+1)

JΔX=m

ΔX=XnewXold

If, the vapour flow rates are assumed to be the constant, the energy equations need not be solved for, The MES or the n×(c+1) form of equations is used for simulations. This system of Non-Linear equations is solved using Thomas Algorithm for Block Tridiagonal matrices. In, general thomas algorithm for tridiagonal matrices has a time complexity O(n) compared to Gauss Elimination with a time complexity of O(n3). Iterations are carried out until tolerance103.

This iterative scheme is sensitive to the initial Guesses and fails to converge if the initial values are too far from the feasible solution. All simulations have been carried out on MATLAB 2015b. The method discussed before is applicable to simple columns. If the columns are interlinked, then the system no longer remains Tridiagonal [15].

A.3 Summary of fitness functions

Table 7:

Summary of fitness functions tested.

S.NoPGw1w2w3w4r1r2Fitness functionCostDistillateBottomConverged fitness value
115501106107.5w1(1w4(1+ecostw2)1)w3|dist(1)0.99|+w4|bot(2)0.99|r1+r23.681070.8970.8131.65
230601010810-w1(1(1+ecostw2)1)w3(1+|dist(1)0.99|+|bot(2)0.99|)1r1+r23.731070.8040.73112.87
3306010-7.5-costw1|dist(1)0.99|+w3|bot(2)0.99|r1*r24.811070.9860.8276.04
43060----cost|dist(1)0.99||bot(2)0.99|r1*r27.801070.9890.8660
55015101087.5-ecostw2w1|dist(1)0.99|+w3|bot(2)0.99|r1*r28.101070.9890.8660
66030101087.5-ecostw2w1|dist(1)0.99|+w3|bot(2)0.99|r1*r28.191070.9890.8660
76030101087.5-ecostw2w1|dist(1)0.99|+w3|bot(2)0.99|r1+10*r23.721070.9870.7311.68
83060101087.5-ecostw2w1|dist(1)0.99|+w3|bot(2)0.99|r1*10*r27.841070.9890.8670
91530101087.5-ecostw2w1|dist(1)0.99|+w3|bot(2)0.99|r1*10*r27.931070.9890.8670
103045101087.5-ecostw2w1|dist(1)0.99|+w3|bot(2)0.99|r1*10*r27.411070.9890.8650
114560101087.5-ecostw2w1|dist(1)0.99|+w3|bot(2)0.99|r1*10*r27.521070.9890.8650
  1. P = Population Size, G = Generation Size

References

[1] Haelssig JB,Considerations for Various Recovery Schemes in Ethanol Production by Fermentation. Industrial & Engineering&Chemistry Research 2008:6185–6191. DOI: 10.1021/ie0715005 .Suche in Google Scholar

[2] Patwardhan AM, Shirsat SP, SD Dawande. Pressure Swing Distillation Simulation and Optimization of Azeotropic IBA + IBAc Mixture. International Journal of Computer Applications :9–13. Vol 66(10), 2013.Suche in Google Scholar

[3] Lewis WK.Dehydrating Alcohol and the like. In: WK Lewis, editors. United States Patent , 1928:676–700. 3.Suche in Google Scholar

[4] Onkar D, Sathe V, Yogesh M. THF-Water mixture separation using Extractive and Azeotropic distillation: Review, Modeling and Simulation in Aspen plus. International Journal of Science and Engineering. 2016;3:10–16.Suche in Google Scholar

[5] Iqbal A, Ahmad SA. Pressure swing distillation of azeotropic mixture – A simulation study. Perspectives in Science. 2016 9;8:4–6. DOI: 10.1016/j.pisc.2016.01.001.Suche in Google Scholar

[6] Yu M, Yang J, Zheng Z. Simulation and optimization for separation of acetonitrile and water by pressure swing distillation and extractive distillation. Journal Of Central South University (Science And Technology) , 2014:2966–2971. Vol 45(9).Suche in Google Scholar

[7] 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. Computers & Chemical Engineering. 2016 12;95:97–107. DOI: 10.1016/j.compchemeng.2016.09.014.Suche in Google Scholar

[8] Modla G. Energy saving methods for the separation of a minimum boiling point azeotrope using an intermediate entrainer. Energy IJSE, 2013:103–109. DOI: 10.1016/j.energy.2012.12.011 .Suche in Google Scholar

[9] Tabari A, Ahmad A. A semicontinuous approach for heterogeneous azeotropic distillation processes. Computers & Chemical Engineering 2015:183–190. DOI: 10.1016/j.compchemeng.2014.12.005 .Suche in Google Scholar

[10] Cortez-Gonzalez J, Segovia-Hernandez JG , Hernandez S, Gutierrez-Antonio C , Briones-Ramírez A ,Rong BG.Genetic Algorithms in the Design of Configurations for Distillation of Quaternary Mixtures using Less than N-1 columns with Thermally Coupling. Computer Aided Chemical Engineering. 2012;30:677–681. DOI: 10.1016/j.compchemeng.2014.12.005.Suche in Google Scholar

[11] Gutiérrez-Antonio C. Multiobjective Stochastic Optimization of Dividing-wall Distillation Columns Using a Surrogate Model Based on Neural Networks. Chemical and Biochemical Engineering Quarterly. 2016 1 3;29(4):491–504. DOI: 10.15255/CABEQ.2014.2132.Suche in Google Scholar

[12] Franke MB. MINLP optimization of a heterogeneous azeotropic distillation process: Separation of ethanol and water with cyclohexane as an entrainer. Computers & Chemical Engineering. 2016 6;89:204–221. DOI: 10.1016/j.compchemeng.2016.03.027.Suche in Google Scholar

[13] Ramanathan SP. Mukherjee S. Dahule RK. Ghosh S. Rahman I. Tambe SS, Ravetkar DD Kulkarni BD Optimization of Continuous Distillation Columns Using Stochastic Optimization Approaches. Chemical Engineering Research and Design. 2001 4;79(3):310–322. DOI: 10.1205/026387601750281671.Suche in Google Scholar

[14] Rahman I, Ahmad A, Kumar P, Kulkarni BD.Optimization of a Continuous Process for the Recovery of Lactic Acid Using Differential Evolution Algorithm. Chemical Product and Process Modeling. 2008 2 3;3(1). DOI: 10.2202/1934-2659.1111.Suche in Google Scholar

[15] Hofeling BS, Seader JD.A modified Naphtali-Sandholm method for general systems of interlinked, multistaged separators. AIChE Journal. 1978 11;24(6):1131–1134. DOI: 10.1002/aic.690240631.Suche in Google Scholar

[16] Beebe AH, Coulter KE, Lindsay RA, Baker EM.Equilibria in Ethanol-Water System at Pressures Less Than Atmospheric. Industrial & Engineering Chemistry 1942:1501–1504.10.1021/ie50396a019 .10.1021/ie50396a019Suche in Google Scholar

[17] Voutsas EC, Pamouktsis C, Argyris D, Pappa GD.Measurements and thermodynamic modeling of the ethanol–water system with emphasis to the azeotropic region. Fluid Phase Equilibria 2011:135–141.10.1016/j.fluid.2011.06.009 66.10.1016/j.fluid.2011.06.009Suche in Google Scholar

[18] Zhang Z, Ming Lv, Huang D, Jia P, Sun D, Li W. Isobaric Vapor–Liquid Equilibrium for the Extractive Distillation of Acetonitrile + Water Mixtures Using Dimethyl Sulfoxide at 101.3 kPa. Journal of Chemical & Engineering Data 2013:3364–3369.10.1021/je400531a 89.10.1021/je400531aSuche in Google Scholar

[19] Babu BV, Angira R, Nilekar A.Computers and Chemical Engineering. 2005;29:1041–1045.10.1016/j.compchemeng.2004.11.010Suche in Google Scholar

[20] Hegerty B , Hung C , Kasprak K.A Comparative Study on Differential Evolution and Genetic Algorithms for Some Combinatorial Problems. proceedings of 8th Mexican International Conference on Artificial Intelligence 2009.Suche in Google Scholar

[21] Skiborowski M, Rautenberg M, Marquardt W., A Hybrid Evolutionary–Deterministic Optimization Approach for Conceptual Design. Industrial & Engineering Chemistry Research. 2015 10 21;54(41):10054–10072.10.1021/acs.iecr.5b01995.10.1021/acs.iecr.5b01995Suche in Google Scholar

[22] Tun LK, Matsumoto H., Application Methods for Genetic Algorithms for the Search of Feed Positions in the Design of a Reactive Distillation Process. Procedia Computer Science. 2013;22:623–632.10.1016/j.procs.2013.09.143.10.1016/j.procs.2013.09.143Suche in Google Scholar

[23] Shahandeh H, Ivakpour J, Kasiri N, Internal and external HIDiCs (heat-integrated distillation columns) optimization by genetic algorithm. Energy. 2014 1;64:875–886.10.1016/j.energy.2013.10.042.10.1016/j.energy.2013.10.042Suche in Google Scholar

[24] Sinnott RK.Coulson and Richardson's Chemical Engineering., 1993;6:567 4th edition Oxford: Pergamon.Suche in Google Scholar

Received: 2017-3-15
Revised: 2017-7-25
Accepted: 2017-7-25
Published Online: 2017-10-13

© 2018 Walter de Gruyter GmbH, Berlin/Boston

Heruntergeladen am 30.11.2025 von https://www.degruyterbrill.com/document/doi/10.1515/cppm-2017-0007/pdf?lang=de
Button zum nach oben scrollen