CHEMICAL ENGINEERING TRANSACTIONS VOL. 61, 2017 A publication of The Italian Association of Chemical Engineering Online at www.aidic.it/cet Guest Editors: Petar S Varbanov, Rongxin Su, Hon Loong Lam, Xia Liu, Jiří J Klemeš Copyright © 2017, AIDIC Servizi S.r.l. ISBN 978-88-95608-51-8; ISSN 2283-9216 Estimation of the Number of Distillation Sequences with Dividing Wall Column for Multi-component Separation Yunlu Zhanga, Guangyue Hanb, Wei Suna* aBeijing Key Lab of Membrane Science and Technology, College of Chemical Engineering, Beijing University of Chemical Technology, 100029 Beijing, China bDepartment of Mathmatics, The University of Hongkong, Hongkong, China sunwei@mail.buct.edu.cn As an important separation unit, distillation column is widely applied in petrochemical and other process industry. For separating multicomponent mixtures, distillation is conducted sequentially in industry. Both individual column and distillation sequence optimization are efficient ways for saving energy consumption. Distillation sequence is usually evaluated by number of distillation subproblems, subgroups and distillation sequences. Distillation sequence has been well studied based on simple column assumption. Dividing wall column (DWC), which is an atypical distillation column for separating a multicomponent feed mixture into three output streams, as a thermally coupled distillation column, has been proposed and applied in distillation sequence. Usually sharp split is also assumed in most literatures on DWCs. A distillation sequence with DWC will give more number of feasible sequences. It is important to estimate the total available number of distillation sequences theoretically. In this work, distillation sequences with both simple column and DWC are considered. Inferential deduction method has been used to explore the number of distillation sequences for multi-component sharp splits. The three general term formulas are obtained with the assumption of sharp split. Under different assumptions, the corresponding numbers of distillation sequences are also discussed. 1. Introduction Distillation is the primary separation technique widely adopted in process industry, by which final products can be obtained, while most other techniques, such as absorption and extraction, need further processing (Mustafa et al., 2014). Traditional distillation column is assumed to separate one feed stream into two output streams. For the separation of multicomponent mixture, a sequence of distillation columns will be needed. Thus, a large number of columns are applied in process industry. It is reported that there are over 40,000 distillation columns all over the world (Kiss et al., 2013). Despite its flexibility and wide application, one important concern on it is its considerable energy consumption, which can account for more than 50 % of plant operating cost (Kiss et al., 2007). Therefore, many research efforts have been concentrated on energy-efficient distillations in terms of both individual column and distillation sequence. Distillation sequence synthesis is one of a significant way to save energy in distillation process. The estimation of the number of all possible sequences is very necessary information in searching an optimal distillation sequence in process synthesis. Distillation sequence has been well studied based on simple column assumption (Muhammad et al., 2015), i.e. in each column one feed is separated into two streams without component mixing between two output streams, which is also named as sharp split. The general term formula of the simple column distillation sequence number has been achieved (Thompson and King, 1972). Non-sharp split has also been studied in recent years. Systematic synthesis of functionally distinct distillation systems which including non- sharp splits for five-component separations is presented by a step-by-step enumeration method in the work (Rong et al., 2005). It was conducted for a specific industrial application, and there are more than one solutions for a single case. However, the discussion based on sharp split assumption has theoretical significance to general separation problems. Dividing wall column (DWC) was first proposed by Wright (1949). It is an atypical distillation column with an internal, vertical partition wall, which effectively accommodates two conventional distillation columns into one DOI: 10.3303/CET1761055 Please cite this article as: Zhang Y., Han G., Sun W., 2017, Estimation of the number of distillation sequences with dividing wall column for multi-component separation, Chemical Engineering Transactions, 61, 343-348 DOI:10.3303/CET1761055 343 shell for separating a multicomponent feed mixture into three output streams (Mohamad et al., 2015). As a thermally coupled distillation column, it is reported that about 30 % energy consumption and equipment investment cost can be saved by the application of DWC (Triantafyllou and Smith, 1992). When including DWC, the number of columns in each sequence can be reduced. However, at the meantime, the number of feasible sequences will increase rapidly, almost beyond the speed of exponential growth. Process synthesis in industrial practice showed that the economic performance of the distillation sequence highly depends on the completeness of the sequence alternatives under consideration. The sequence with DWC only at the last separation step was discussed for a five component sequence (Du et al., 2016) by enumeration. More general study on the number of distillation sequences with DWC hasn’t been reported yet. Thus, it is important to theoretically estimate the total available number of distillation sequences including both simple column and DWC. In this work, distillation sequences with both simple column and dividing wall column (DWC) are considered. An inferential deduction to calculate the number of distillation sequences, distillation subproblems, and distillation subgroups for multi-component separation based on sharp split has been presented. Three general term formulas are obtained and discussed under different assumptions. 2. Problem description Distillation sequence is usually described by the number of distillation sequences, distillation subproblems, and distillation subgroups. Following the expression used in Analysis and Synthesis of Chemical Process (Zhang et al., 2011), the number of the distillation sequences is defined as 𝑆𝑅 which refers to the possible number of sequences for separating 𝑅 components into 𝑅 pure products by simple column. The number of distillation subproblems is defined as 𝑈 which corresponds to a possible number of separation units. Separation sequences are different combinations of subproblems. The number of distillation subgroups is defined as 𝐺 which is number of streams with adjacent components in a multi-component separation, as the feed or final product of each separator, or subproblem. In this work, only sharp split is considered in both simple column and DWC. For the convenience of understanding, several new variables are defined in the following discuss. 2.1 Distillation sequences The recursive formula and general term formula of the number of distillation sequence with simple column only were given in (Seider et al., 1999) as follows: 𝑆𝑅 = ∑ 𝑆𝑗 𝑅−1 𝑗=1 𝑆𝑅−𝑗 (1) 𝑆𝑅 = [2(𝑅−1)]! 𝑅!(𝑅−1)! (2) where 𝑅 is the number of components (𝑅 = 1,2,3, … , 𝑛). Considering both simple column and DWC, 𝑎𝑅 is used to replace the 𝑆𝑅 in Eq(1), which represents the number of separation sequence for 𝑅 component separation if a simple column is picked at the current step, and 𝑏𝑅 as the number of separation sequence for 𝑅 component separation if a DWC is picked at the current step. Choosing DWC at the current step, there are (𝑅 − 2) separation choices for the first column. Let the number of components appearing at the top product of the column be 𝑗, the number of components appearing at the side product of the column be 𝑘, then the number of components appearing at the bottom product of the column will be (𝑅 − 𝑗 − 𝑘). Therefore, the recursive formula will be three items multiplication as in Eq(4). The calculation is valid from 𝑅 ≥ 2. 𝑎𝑅 = ∑ 𝑆𝑗 𝑅−1 𝑗=1 𝑆𝑅−𝑗 (3) 𝑏𝑅 = ∑ 𝑆𝑗 ∑ 𝑆𝑘 𝑅−𝑗−1 𝑘=1 𝑅−2 𝑗=1 𝑆𝑅−𝑗−𝑘 (for 𝑅 ≥ 2 ) (4) The total number of distillation sequences 𝑆𝑅 can be calculated as follows: 𝑆𝑅 = 𝑎𝑅 + 𝑏𝑅 (5) For calculation convenience, it is assumed that 𝑎1 = 1, while there is actually no column needed if only one component exits, 𝑏1 = 0, 𝑏2 = 0, as one or two component separation won’t be able to use a DWC. Thus, 𝑆1 = 𝑎1 + 𝑏1 = 1 (6) Obviously, 𝑎2 = 1, so that 344 𝑆2 = 𝑎2 + 𝑏2 = 1 (7) With these initial values, recursive formulas can be executed by R programming. The result of the number of the distillation sequence is shown in Table 1. Table 1: The result of the number of the distillation sequences with DWC 𝑅 2 3 4 5 6 7 8 9 10 11 … 𝑆𝑅 1 3 10 38 154 654 2871 12925 59345 276835 … Fortunately, the integer sequence of the result has been studied by Hanna (2005), and the general term formula (10) is achieved. For the expression simplicity, it is assumed that: 𝑅 = 𝑄 + 1 (8) 𝑆𝑅 = 𝑆𝑄+1 (9) 𝑆𝑄+1 = ∑ 𝐶2𝑄−𝑘 𝑄+𝑘 𝐶𝑄+𝑘 𝑘 𝑄+1 [ 𝑄 2 ] 𝑘=0 (10) 2.2 Distillation subproblems As the four component separation example shown in Figure 1, elements in the first column of Figure 1(a) and 1(b) represent a four component subgroup with different separation choices, i.e. subproblem. Elements in the second column of each subplot represent the outcome of subproblem in its left column; same logic applies in the third column too. (a) Simple column subproblems (b) DWC subproblems Figure 1: Four component separation subproblems It can be counted out that the number of subproblems for four component separation with simple column is 10, which can be calculated by the general term formula (Zhang et al., 2011) written as follows: 𝑈 = ∑ 𝑗(𝑅 − 𝑗)𝑅−1𝑗=1 = 𝑅(𝑅−1)(𝑅+1) 6 (11) Considering DWC into the problem, 𝑈𝑎 is used to replace the 𝑈 in Eq(11), which represents the subproblem number with simple column, and 𝑈𝑏 as the subproblem number with DWC. So that the 𝑈 can be calculated as follows: 𝑈 = 𝑈𝑎 + 𝑈𝑏 (12) Obviously, at the step of 𝑅 component separation, there are 𝐶𝑅−1 2 ways which is a kind of combination number to apply DWC. Thus, 𝑢𝑅 is defined at current separation step. 𝑢𝑅 = 𝐶𝑅−1 2 (13) When 𝑛 < 3, 𝑢𝑏1 = 0, 𝑢𝑏2 = 0 ,𝑈𝑏 = 0, 𝑈 = 𝑈𝑎 = 𝑅(𝑅−1)(𝑅+1) 6 (14) 345 When 𝑛 ≥ 3, 𝑈𝑎 is remained as Eq(11), while 𝑈𝑏 = 1𝑢𝑅 + 2𝑢𝑅−1 + 3𝑢𝑅−2 + 4𝑢𝑅−3 + ⋯ + (𝑛 − 2)𝑢3 = ∑ 𝑗𝑢𝑅+1−𝑗 𝑅−2 𝑗=1 (15) So that 𝑈 can be expressed as the general term formula followed: 𝑈 = 𝑈𝑎 + 𝑈𝑏 = (𝑛−1)𝑛(𝑛+1)(𝑛+2) 24 (16) 2.3 Distillation subgroups In terms of the definition, the general term formula of 𝐺 (Zhang et al., 2011) is constant no matter which type of column is included. 𝐺 = ∑ 𝑖𝑅𝑗=1 = 𝑅(𝑅+1) 2 (17) Take the four component separation as an example, all possibilities are shown in Figure 2. Table 2: The number of distillation sequences, distillation subproblems, and distillation subgroups 𝑅 𝑆𝑅 𝑈 𝐺 2 1 1 3 3 3 5 6 4 10 15 10 5 38 35 15 6 154 70 21 7 654 126 28 8 2871 210 36 9 12925 330 45 10 59345 495 55 11 276835 715 66 Figure 2: Four component separation subgroups The number of distillation sequences, distillation subproblems, and distillation subgroups mentioned above can be calculated by the general term formulas. The results are shown in the Table 2. 3. Extension Three general term formulas are obtained for separation sequence with simple column and DWC. This analysis can be applied if more complicated column structure is considered, e.g. a column with one feed stream and more number of output streams. 3.1 Case for more number of output streams from a column A special column structure with one feed stream and four output streams is considered as an example here. Choosing a column with one feed stream and four output streams at the current step, there are (𝑅 − 3) separation choices for this type of column. Let the number of components appearing at the top product of the column be 𝑗, the number of components appearing at the first side product of the column be 𝑘, the number of components appearing at the second side product of the column be 𝑙, then the number of components appearing at the bottom product of the column will be (𝑅 − 𝑗 − 𝑘 − 𝑙). Therefore, the recursive formula will be four items 346 multiplication as Eq(18). The 𝑆𝑅 can be calculated as Eq(19). The 𝑎𝑅 and 𝑏𝑅 are defined in Eq(3) and Eq(4) respectively. 𝑐𝑅 = ∑ 𝑆𝑗 ∑ 𝑆𝑘 𝑅−𝑗−2 𝑘=1 ∑ 𝑆𝑙 𝑅−𝑗−𝑘−1 𝑙=1 𝑅−3 𝑗=1 𝑆𝑅−𝑗−𝑘−𝑙 (18) 𝑆𝑅 = 𝑎𝑅 + 𝑏𝑅 + 𝑐𝑅 (19) 𝑐1 = 0, 𝑐2 = 0, 𝑐3 = 0 as one, two or three component separation won’t be able to use this type of column. With the initial value 𝑎1 = 1, 𝑎2 = 1, 𝑏1 = 0, 𝑏2 = 0, the recursive formulas can be executed by R programming. The result of the number of the distillation sequence is shown in Table 3. Table 3: The result of the number of the distillation sequences with four output streams column 𝑅 2 3 4 5 6 7 8 9 10 11 … 𝑆𝑅 1 3 11 44 189 850 3951 18832 91542 452075 … With the same method mentioned in the section 2.1, the general term formula (Vladimir, 2011) can be obtained as Eq(20). The general binomial distribution definition is available at 𝑆𝑅 = 1 𝑅 ∑ [∑ 𝑏𝑖𝑛𝑜𝑚𝑖𝑎𝑙(𝑗, 𝑅 − 3𝑘 + 2𝑗 − 1)𝑘𝑗=0 ∙ 𝑏𝑖𝑛𝑜𝑚𝑖𝑎𝑙(𝑘, 𝑗) 𝑅−1 𝑘=1 ] ∙ 𝑏𝑖𝑛𝑜𝑚𝑖𝑎𝑙(𝑅 + 𝑘 − 1, 𝑅 − 1) (for 𝑅 ≥ 2 ) (20) From the discussion above, these series of recursive formulas have a possibility to extend to a separator with any number of feed streams and any number of output streams, which possibly exists in practical separation sequence. 3.2 Case for DWC only included at the last separation step In certain industrial process, pure product is more attractive. Due to the operating challenge, DWC is only used to separate three components into three pure products at the last separation step in a sequence. In this case, 𝐺 is constant with the Eq(17) as always. 𝑈𝑏 in Eq(15) only left with the last item, shown as Eq(21), with the definition of 𝑢𝑅 in Eq(13), 𝑈𝑏 can be calculated with the result 𝑛 − 2. The 𝑈 can be obtained with the Eq(22). 𝑈𝑏 = (𝑛 − 2)𝑢3 (21) 𝑈 = 𝑅(𝑅−1)(𝑅+1) 6 + (𝑛 − 2) (22) DWC is used only at the last separation step, in other words, DWC is used only in three components separation. Thus, 𝑆𝑅 defined in Eq(5) is suitable only when 𝑅 = 3, i.e. the sequence with DWC is only considered when 𝑅 = 3 . 𝑆𝑅 will be calculated according to Eq(1), as no DWC will be applied at other separation step. Based on Eq(3) and Eq(4), 𝑆3 can be expressed as Eq(23). With the initial value of 𝑆1 and 𝑆2, 𝑆3 can be calculated as follows: 𝑆3 = 𝑆1𝑆2 + 𝑆2𝑆1 + 𝑆1𝑆1𝑆1 = 1×1 + 1×1 + 1×1×1 = 3 (23) Set 𝑆3 with the new initial value 3, the recursive formula given in Eq(1) can be executed to get new results by R programming, which is shown in Table 4. Table 4: The result of the number of the distillation sequences with DWC used at the last separation step 𝑅 2 3 4 5 6 7 8 9 10 11 … 𝑆𝑅 1 3 7 20 63 208 711 2,496 8,944 32,578 … With the same method mentioned in the section 2.1, the general term formula (Vladimir, 2014) of the number of the distillation sequences can be obtained as Eq(26). The general binomial distribution definition is also available at WolframMathWorld (2017). For the expression simplicity, it is assumed that: 𝑅 = 𝑄 + 1 (24) 𝑆𝑅 = 𝑆𝑄+1 (25) 𝑆𝑄 = ∑ 𝑏𝑖𝑛𝑜𝑚𝑖𝑎𝑙(𝑄−2𝑚+1,𝑚)∙𝑏𝑖𝑛𝑜𝑚𝑖𝑎𝑙(2𝑄−4𝑚,𝑄−2𝑚) 𝑄+1−2𝑚 𝑄 2 𝑚=0 (for 𝑄 ≥ 1 ) (26) 4. Conclusions The number of distillation sequences, distillation subproblems, and distillation subgroups with both simple column and DWC are discussed under the only consumption of sharp split. Corresponding general term formulas are obtained. The development procedure of recursive formulas can be extended to a separator with 347 any number of feed streams and any number of output streams, while the number of distillation subgroups remains constant no matter how complicated the column is included. Under the assumption of DWC used only at the last separation step, the corresponding formulas are obtained as well. If non-sharp split is included in a process, the related work for multi-component separation will be even more complicated, which can be further discussed for a given separation requirement. The impact of thermodynamic properties of components in a separation system hasn’t been investigated yet, which could be helpful in reducing the number of feasible separation sequences. Acknowledgments The authors gratefully acknowledge the following institutions for their support: the National Natural Science Foundation of China (Grant No.21576015); the National Basic Research Program of China, 973 program (Grant No. 2013CB733600). References Du Zengzhi, Wang Jianhong, Zhai Chi, Zhang Yunlu, Zhu Zhexi, Sun Wei, 2016, < aiche.org/conferences/aiche- annual-meeting/2016/proceeding/paper/449j-study-on-distillation-sequence-dividing-wall-column>, accessed 13.02.2017. Kiss, A.A., 2013, Design, 2. Control and Economics of Distillation, Advanced Distillation Technologies: Design, Control and Applications, John Wiley & Sons, Ltd, Hoboken, US, 37-65. Kiss A.A., Pragt H., Strien C.V., 2007, Overcoming Equilibrium Limitations in Reactive Dividing-Wall Columns, Computer Aided Chemical Engineering, 24(07), 467-472. Mustafa M.F., Samad N.A.F.A., Ibrahim K.A., Hamid M.K.A., Mustafa M.F., 2014, Methodology Development for Designing Energy Efficient Distillation Column Systems, Energy Procedia, 61, 2550-2553. Othman M.R., Amran U.I., Rangaiah G.P., 2015, Process Optimization of DWC for Fatty Acid Fractionation using Taguchi Methods of Experimental Design, Chemical Engineering Transactions, 45,925-930, doi: 10.3303/CET1545155. Paul D. Hanna, 2005, , accessed 15.02.2017. Rong B.G., Kraslawski A., Turunen I., Systematic Synthesis of Functionally Distinct New Distillation Systems for Five-Component Separations, 2005, Comput. Chem. Eng., 20, 823–828. Thompson, R.W., C.J. King, 1972, AIChE J. 18, 941. Seider W.D., Seader J.D., Lewin D.R., 1999, Process Design Principles Synthesis, Analysis, and Evaluation, New York: JOHN WILEY& SONS, INC. Thompson R.W., King C.J., AIChE J. 1972, 18, 941. Triantafyllou C., Smith R., 1992, The design and optimisation of fully thermally coupled distillation columns: Process design, Chemical Engineering Research & Design, 70, 118-132. Vladimir Kruchinin, 2011, , accessed 24.02.2017. Vladimir Kruchinin, 2014, , accessed 26.02.2017. WolframMathWorld, accessed 08.06.2017 Wright R O., 1949, Fractionation apparatus, United States patent, US 2471134. Zhang Weidong, Sun Wei, Liu Junteng, 2011, Analysis and Synthesis of Chemical Process, Beijing: Chemical Industry Press. Zaine M.Z., Mustafa M.F., Ibrahim K.A., Ibrahim N., Hamid M.K.A., 2015, Sustainable Energy Efficient Distillation Columns Sequence, Chemical Engineering Transactions, 45, 1207-1212, doi: 10.3303/CET1545202. 348