Optimal Design of a Hydrolysis Sugar Membrane Purification System Using a Superstructure-Based Approach

As an alternative to gasoline, bioethanol can be produced from lignocellulosic biomass through hydrolysis using an ionic solution containing zinc chloride (ZnCl2). This method allows for a high yield of glucose from lignocellulose, but entails the removal of ZnCl2 from the hydrolysate using multiple nanofiltration membranes before the fermentation of glucose. This paper presents a mathematical technique for designing such a multistage membrane separation system. The optimization model for the synthesis of membrane networks is based on a superstructure with all feasible interconnections between the membrane units, and consists of mass balances, logical constraints and product specifications. A case study of the separation of a bagasse hydrolysis solution is used to demonstrate the application of the proposed model. Results show that using both types of nanofiltration membranes allows higher ZnCl2 removal ratios at each membrane unit, hence a decrease in the number of membrane units required and a reduction of about 35% in capital cost compared to the cases in which only one membrane type is used. Further analysis is performed to examine the effect of membrane performance on the economics of the separation system.


Introduction
With declining petroleum resources and the pursuit of clean energy, researchers have paid more attention to developing efficient technology for alternative fuels. Biomassderived sugars are considered a crucial raw material for liquid biofuel production through biological and chemical conversion. To that end, conventional sugar crops may not be preferred due to the competition with food uses. Alternatively, lignocellulosic biomass (e.g., energy crops, agricultural waste and forest residues) can be a good source of sugars. Lignocellulose is one of the most abundant plant-based materials on earth, amounting to approximately 200 billion metric tons annually [1]. Technological advancements in converting lignocellulosic biomass to sugars for biofuel production are expected to provide a promising solution to energy autonomy.
Lignocellulosic biomass consists of complicated carbohydrates with three basic polymers: hemicellulose, cellulose and lignin. Researchers at the National Renewable Energy Laboratory (NREL) of the U.S. Department of Energy have established five conversion technology platforms for base chemicals [2]. Among those, the sugar platform is the most promising process for the production of alternative fuels (e.g., ethanol) and other chemicals. Traditional lignocellulose hydrolysis technology can be divided into acid hydrolysis and Processes 2021, 9,168 2 of 16 enzyme hydrolysis. The latter is time-consuming and requires suitably adjusted conditions for bacteria to survive. Acid hydrolysis can be further divided into concentrated acid hydrolysis and diluted acid hydrolysis. The former usually uses sulfuric acid with a concentration of 30-70%. At lower temperatures, the cellulose can be completely dissolved and converted to oligosaccharides, which can then be hydrolyzed to glucose with a yield of up to 90%. However, the main drawbacks of the hydrolysis of cellulose using concentrated sulfuric acid are the high cost of anticorrosive devices and the necessity of recycling used acid to avoid corrosion and environmental pollution. These factors limit the application of this method.
The Purdue Research Foundation patented the hydrolysis of cellulose to glucose using zinc chloride (ZnCl 2 ) [3]. More recently, the Industrial Technology Research Institute (ITRI) developed an ionic ZnCl 2 solution to improve the yield of converting cellulose to glucose. This ionic solution is much cheaper than conventional ionic liquids, and can achieve a more than 90% high yield of glucose from lignocellulosic biomass [4]. However, the resulting glucose product stream has a high concentration of ZnCl 2 , which needs to be removed before the fermentation of glucose. The separation of glucose and ZnCl 2 thus became an important issue.
Purification of reducing sugars after ionic-solution hydrolysis is considered a critical step of the entire bioethanol production process, because the separation of the salt from glucose is often expensive and may encounter contamination from solvent extraction. The ITRI has studied the separation of high-concentration ZnCl 2 from the hydrolysis sugar product by nanofiltration (NF) and proven the capability of this approach [5]. The separation performance of NF membranes is between that of ultrafiltration (UF) and reverse osmosis (RO) membranes; however, NF membranes have a wider range of applications than conventional membranes [6]. So far, the commercial application of NF membranes is mainly in desalination, especially in food processes. Some previous studies [7][8][9][10] have focused on the effects of solute charge, pH and temperature on the retention effect of NF membranes.
With many results, NF has proven itself to be a good candidate for sugar/salt separation [5]. However, the required product specifications cannot be met with a single membrane; a number of membrane units are needed to achieve the desired effect. This presents the problem of designing an NF separation system featuring a network of membrane units, with some arranged in series and some in parallel. Moreover, although simulation is useful in assessing the performance of a given network configuration, optimization is needed to find the optimal network configuration and determine the performance targets in terms of capital cost, water consumption and the number of membrane units. This work aims to develop an optimization approach to the design of NF separation systems.
Superstructure optimization is a powerful tool for process design, enabling systematic and integrated analysis of alternative structures [11]. Superstructure-based approaches have been applied to the synthesis of various membrane separation networks. El-Halwagi [12] introduced the problem of synthesizing RO networks and developed a mixed integer nonlinear programming (MINLP) model for desalination and dephenolization applications. Voros et al. [13,14] proposed an alternative superstructure resulting in an easier-to-tackle nonlinear programming (NLP) model. Lu et al. [15] presented a simplified superstructure and a pressure vessel model for the design of multistage RO systems for seawater desalination under different feed concentrations and product specifications. Lu et al. [16] studied an RO-based seawater desalination system for producing freshwater of multiple concentration levels. A numerical model describing the membrane transport of RO in spiral-wound modules was introduced by Du et al. [17]. Optimization approaches to the design of RO networks under boron constraints were also proposed [18,19]. Alnouri and Linke [20] proposed compact superstructure representations and a structured synthesis strategy to facilitate the identification of optimal membrane networks. Abejón et al. [21] optimized an RO network for the ultrapurification of hydrogen peroxide for electronic applications. Chauhan et al. [22] developed a superstructure-based optimization approach for hybrid seawater RO-NF membrane desalination and salt production. More recently, Parra et al. [23] proposed the use of metamodels to speed up the solution of the nonlinear RO network optimization problem.
Apart from desalination, membrane networks are also designed for gas separation. Qi and Henson [24] developed an MINLP model for spiral-wound membrane systems for the separation of CO 2 /CH 4 mixtures in natural gas treatment and enhanced oil recovery. Kookos [25] considered countercurrent hollow-fiber membrane modules and presented a targeting approach to optimize the membrane material for the production of nitrogen or oxygen enriched air. Various flow patterns (cross, countercurrent and concurrent) and all possible membrane network configurations are embedded in the superstructure of Uppaluri et al. [26]. These authors later presented a targeting-design approach to optimize the pressure and network configuration of gas separation membranes [27]. Arias et al. [28] presented an NLP model to determine the optimal number of membrane stages and operating conditions for CO 2 capture in coal-fired power plants. Shafiee et al. [29] presented an automated flowsheet synthesis method to optimize a hybrid membrane/cryogenic CO 2 capture process. Ramírez-Santos et al. [30] presented a global optimization approach for multistage membrane separation systems for capturing CO 2 from blast furnace gas. Lee et al. [31] developed and applied an optimization framework to the design of multistage membrane processes for CO 2 capture from flue gas.
Superstructure-based approaches were also developed for the synthesis of thermal membrane distillation networks [32] and bio-separation networks [33]. More recently, Baaqeel and El-Halwagi [34] considered the integration of multiple effect distillation and membrane distillation for multi-period capacity planning in desalination networks. Oke et al. [35] considered the use of thermally driven membrane distillation to purify flowback wastewater for fracturing water reuse. Mores et al. [36] carried out an optimization study of membrane-based separation systems with up to four stages for CO 2 capture from flue gas of power plants. Tao et al. [37] carried out an integrated design of multi-stage membrane separation systems with uncertain concentrations and flowrates, proposing three landfill gas purification strategies. Martín-Hernández et al. [38] developed a superstructure model integrating amine absorption, pressure swing adsorption and membrane separation to determine the optimal process and operating conditions for biomethane production. However, as yet there was no such approach developed for the separation of hydrolysis solutions in bioethanol production. Chemical hydrolysis using a ZnCl 2 solution has been recognized as a promising method for converting lignocellulosic biomass to glucose [4]. The subsequent removal of ZnCl 2 from the hydrolysate using NF membranes (for the hydrolysis sugar solution to be fermented) is critical to the economic viability of bioethanol production through chemical hydrolysis [5]. In this paper, a mathematical model is developed for the design and optimization of NF membrane networks for sugar/salt separation. The optimization model is based on a generic superstructure that embeds all possible interconnections between the membrane units, and is useful in determining the optimal network configuration of the membrane separation system. The rest of the paper is organized as follows. A formal problem statement is given in Section 2. The model formulation is presented in Section 3. A case study is then presented in Section 4 to illustrate the proposed approach. Finally, conclusions and prospects for future work are given in Section 5.

Problem Statement
The problem addressed in this paper can be formally stated as follows.
• Given a feed mixture containing a set of components i ∈ I. The flowrate of the feed and the concentrations of the components are known. • A set of membrane units j ∈ J are used to separate the feed mixture. The performance data for the membranes (i.e., cut fractions and recoveries/removal ratios) are specified.
• In this work, NF membranes are used to separate ZnCl 2 and glucose in a bagasse hydrolysis solution. The objective is to synthesize an optimal NF membrane network that achieves the required product specifications at the minimum capital cost. Figure 1 shows a schematic diagram of an NF membrane unit, where the inlet stream is separated into a retentate and a permeate. Equations (1) and (2) describe the flowrate and component balances, respectively.

Model Formulation
where f in j , f R j and f P j are the flowrates of the inlet stream, retentate and permeate of unit j, respectively; c in ij , c R ij and c P ij are the respective concentrations of component i. The cut fraction (θ j ) and the recovery (r ij ) are defined as in Equations (3) and (4), respectively. The recovery of the components depends on the membrane type.
Processes 2021, 9, x FOR PEER REVIEW 4 of 16  A set of membrane units ∈ J are used to separate the feed mixture. The performance data for the membranes (i.e., cut fractions and recoveries/removal ratios) are specified.  In this work, NF membranes are used to separate ZnCl2 and glucose in a bagasse hydrolysis solution. The objective is to synthesize an optimal NF membrane network that achieves the required product specifications at the minimum capital cost. Figure 1 shows a schematic diagram of an NF membrane unit, where the inlet stream is separated into a retentate and a permeate. Equations (1) and (2) describe the flowrate and component balances, respectively.

Model Formulation
where in , R and P are the flowrates of the inlet stream, retentate and permeate of unit j, respectively; in , R and P are the respective concentrations of component i. The cut fraction ( ) and the recovery ( ) are defined as in Equations (3) and (4), respectively. The recovery of the components depends on the membrane type.  The modeling of a multistage NF system is based on the superstructure in Figure 2, where the membrane units are divided into three sections: the unit handling the feed mixture (j = MF), units handling retentates, or MR units (j ∈ J MR = {MR1, MR2, · · ·}), and units handling permeates, or MP units (j ∈ J MP = {MP1, MP2, · · ·}). Because the feed mixture is the only input to unit MF, the inlet flowrate and concentration of unit MF are the feed flowrate (F F ) and concentration (C F i ), as stated in Equations (5) and (6).
Equations (7) and (8) describe the inlet flowrate and component balances for unit j ∈ J MR , respectively. The inlet stream of unit j ∈ J MR may consist of the retentate from the previous unit (j − 1), retentates from units j ∈ J MP , and freshwater, which is added to control the concentration and viscosity of the solution.
where f j j is the flowrate from unit j to unit j; f FW j is the amount of freshwater added to unit j; C FW i is the concentration of component i in freshwater. Equations (9) and (10) describe the inlet flowrate and component balances for unit j ∈ J MP , respectively. The inlet stream of unit j ∈ J MP may consist of the permeate from the previous unit and permeates from units j ∈ J MR .
The outlet flowrate balances for unit MF are given by Equations (11) and (12). The retentate of unit MF may be sent to unit MR1 or collected as part of the retentate product. Similarly, the permeate of unit MF may be sent to unit MP1 or collected as part of the permeate product.
where f RP j is the flowrate from unit j to the retentate product, and f PP j is the flowrate from unit j to the permeate product.
Equations (13) and (14) describe the outlet flowrate balances for units j ∈ J MR . The retentates of these units may be sent to the next unit (j + 1) or collected as part of the retentate product, except the last unit (j = MRn), of which the retentate is part of the product. The permeates of units j ∈ J MR may be sent to units j ∈ J MP as inputs, collected as part of the permeate product, or purged as waste.
where f W j is the flowrate from unit j to waste. Equations (15) and (16) describe the outlet flowrate balances for units j ∈ J MP . The retentates of these units may be sent to units j ∈ J MR as inputs, collected as part of the retentate product, or purged as waste. The permeates of units j ∈ J MP may be sent to the next unit or collected as part of the permeate product, except the last unit (j = MPm), of which the permeate is part of the product.
Equations (17) and (18) describe the flowrate and component balances for the retentate product, respectively. The flowrate and component balances for the permeate product are given by Equations (19) and (20), respectively.
Processes 2021, 9,168 6 of 16 where f R,prod is the retentate product flowrate; f P,prod is the permeate product flowrate; c R,prod i and c P,prod i are the respective concentrations of component i. These variables are constrained by the product specifications.
In addition, the flowrate and component balances for the waste are given by Equations (21) and (22), respectively.
where f waste and c waste i are the waste flowrate and concentration (of component i), respectively. Equation (23) imposes the lower (F L j ) and upper limits (F U j ) to the inlet flowrate when unit j is used (y j = 1), and forces the inlet flowrate to zero when unit j is not used (y j = 0).
where y j is a binary variable denoting the existence of unit j.
Equations (24) and (25) state that unit j + 1 will not be used if unit j does not exist. These constraints ensure proper sequencing of the membrane units.
The objective function is to minimize the total capital cost of membrane units, as given in Equation (26). (26) where z represents the objective value; β is the base cost for a membrane unit with an inlet flowrate of F * ; γ is a constant.
Equations (1)-(26) constitute an MINLP model, for which global optimality may not always be guaranteed. The notation used in the formulation is given in Appendix A. In the next section, a case study is presented to illustrate the proposed approach. The MINLP model is implemented and solved in GAMS on a Core i7-9700K, 3.60 GHz processor, using BARON as the solver. Solutions were found with reasonable processing time.
where represents the objective value; is the base cost for a membrane unit with an inlet flowrate of * ; is a constant.
Equations (1)-(26) constitute an MINLP model, for which global optimality may not always be guaranteed. The notation used in the formulation is given in Appendix A. In the next section, a case study is presented to illustrate the proposed approach. The MINLP model is implemented and solved in GAMS on a Core i7-9700K, 3.60 GHz processor, using BARON as the solver. Solutions were found with reasonable processing time.

Case Study
A bagasse hydrolysis solution is to be processed for fermentation in a bioethanol production process. As shown in Figure 3, bagasse is hydrolyzed first using the ITRI's ZnCl2 solution. After the subsequent distillation and multistage filtration, the resulting hydrolysis solution containing 2 wt% glucose, 6 wt% ZnCl2 and 92 wt% water is further treated by NF to remove ZnCl2. Two cases are considered in this case study. Case 1 uses membrane performance data obtained from experiments [5], while case 2 assumes hypothetical parameters to analyze the impact of membrane performance on the optimal network configuration. In either case, the objective is to synthesize the minimum-cost NF membrane network that achieves the desired separation.

Case Study
A bagasse hydrolysis solution is to be processed for fermentation in a bioethanol production process. As shown in Figure 3, bagasse is hydrolyzed first using the ITRI's ZnCl 2 solution. After the subsequent distillation and multistage filtration, the resulting hydrolysis solution containing 2 wt% glucose, 6 wt% ZnCl 2 and 92 wt% water is further treated by NF to remove ZnCl 2 . Two cases are considered in this case study. Case 1 uses membrane performance data obtained from experiments [5], while case 2 assumes hypothetical parameters to analyze the impact of membrane performance on the optimal network configuration. In either case, the objective is to synthesize the minimum-cost NF membrane network that achieves the desired separation.

∈J
where represents the objective value; is the base cost for a membrane unit with an inlet flowrate of * ; is a constant.
Equations (1)-(26) constitute an MINLP model, for which global optimality may not always be guaranteed. The notation used in the formulation is given in Appendix A. In the next section, a case study is presented to illustrate the proposed approach. The MINLP model is implemented and solved in GAMS on a Core i7-9700K, 3.60 GHz processor, using BARON as the solver. Solutions were found with reasonable processing time.

Case Study
A bagasse hydrolysis solution is to be processed for fermentation in a bioethanol production process. As shown in Figure 3, bagasse is hydrolyzed first using the ITRI's ZnCl2 solution. After the subsequent distillation and multistage filtration, the resulting hydrolysis solution containing 2 wt% glucose, 6 wt% ZnCl2 and 92 wt% water is further treated by NF to remove ZnCl2. Two cases are considered in this case study. Case 1 uses membrane performance data obtained from experiments [5], while case 2 assumes hypothetical parameters to analyze the impact of membrane performance on the optimal network configuration. In either case, the objective is to synthesize the minimum-cost NF membrane network that achieves the desired separation.  Figure 3. Bagasse hydrolysis and subsequent separation processes.  Table 1 shows the data for membrane units. These performance data are obtained experimentally using commercially available NF membranes of Dow Chemicals [5]. It can be seen that the removal ratio of ZnCl 2 varies according to the membrane type and the inlet concentration. Therefore, additional constraints in Equations (27)-(31) are needed to assign removal ratios to membrane units:

Case 1: Membrane Network Synthesis Using Experimental Data
Processes 2021, 9, 168 8 of 16 ∑ k∈K ∑ l∈L y jkl = y j ∀j ∈ J (30) where y jkl is a binary variable denoting if membrane unit j is of type k and its inlet ZnCl 2 concentration lies in interval l; U is an arbitrary large number; ε is a small value; RR kl is the ZnCl 2 removal ratio of membrane type k when the inlet ZnCl 2 concentration lies in interval l. Equation (27) is active when the inlet ZnCl 2 concentration lies in interval I, in which case one membrane type is selected for unit j (∑ k∈K y jk,I = 1). Similarly, Equations (28) and (29) are active when the inlet ZnCl 2 concentration lies in intervals II and III, respectively. Equation (30) states that if unit j is used, then its inlet ZnCl 2 concentration lies in one of the concentration intervals and a membrane type is selected. The retention of ZnCl 2 in the retentate of unit j is calculated from the assigned ZnCl 2 removal ratio in Equation (31). With two membrane types, the objective function in Equation (26) is rewritten as in Equation (32).
where β k is the base cost for a membrane unit of type k with an inlet flowrate of F * . Three scenarios are analyzed in this case. Scenario 1 considers the use of only NF1 membranes, scenario 2 only NF2 membranes, while scenario 3 both NF1 and NF2 membranes. Assuming a feed flowrate of 3 t/h (to unit MF), the inlet flowrate of MR units is fixed at the base flowrate of 3 t/h. For each scenario, the NF membrane network is designed to meet the specifications that the retentate product has a minimum glucose/ZnCl 2 ratio of 10 and a minimum glucose recovery of 90%, as given in Equations (33) and (34). No specifications for the permeate product are set in this case.
Solving the MINLP model (Equation (32) subject to Equations (1)-(25), (27)-(31), (33) and (34)) for scenarios 1-3 yields the results in Table 2. It can be seen that the specified glucose/ZnCl 2 ratio of 10 in the retentate product is achieved in all scenarios, while the specified glucose recovery of 90% in the retentate product can be achieved only in scenarios 2 and 3. With a maximum number of 20 for MR units, the maximum achievable glucose recovery (in the retentate product) in scenario 1 is 78.4%, which is determined by maximizing the glucose recovery in Equation (33). The reason that the glucose recovery specification is not met in scenario 1 is that the ZnCl 2 removal ratio of NF1 is only 10% when the inlet ZnCl 2 concentration ≤ 1.5 wt%. Therefore, the glucose recovery specification cannot be met in scenario 1 unless the number of MR units is up to 30. Figure 4 shows the inlet concentrations and the corresponding retention of glucose and ZnCl 2 in MF and MR units. Obviously, most of the inlet ZnCl 2 concentrations of MR units lie in concentration interval I (≤1.5 wt%), in which the ZnCl 2 removal ratio of NF2 is twice that of NF1 (see Table 1). This explains why a better solution is found in scenario 2, where both specifications for the retentate product are achieved with fewer membrane units and a lower capital cost. When both membrane types are used, NF2 can be used for membrane units with inlet ZnCl 2 concentrations in interval I and NF1 in intervals II and III to give higher ZnCl 2 removal ratios, as shown in Table 3. However, to minimize the capital cost, NF1 may be selected instead of NF2 for units with an inlet ZnCl 2 concentration in interval I (MR9 in Table 3), if the product specifications can still be met. Accordingly, scenario 3 has the best solution, with significant cost reductions of 36.9% and 34.8% compared to scenarios 1 and 2, respectively.  Table 3. Selection of membrane types for scenario 3. NF1  3  MR1  NF1  3  MR2  NF1  3  MR3  NF2  3  MR4  NF2  3  MR5  NF2  3  MR6  NF2  3  MR7  NF2  3  MR8  NF2  3  MR9  NF1  3  MR10  NF2  3  MP1  NF1  7.433 Processes 2021, 9, x FOR PEER REVIEW 10 o  Table 3. Selection of membrane types for scenario 3. MF  NF1  3  MR1  NF1  3  MR2  NF1  3  MR3  NF2  3  MR4  NF2  3  MR5  NF2  3  MR6  NF2  3  MR7  NF2  3  MR8  NF2  3  MR9  NF1  3  MR10  NF2  3  MP1  NF1 7.433  Figure 5 shows the retentate concentrations and glucose/ZnCl 2 ratios in MF and MR units. Other details of the NF membrane networks for scenarios 1-3 are given in Tables S1-S3 in the Supplementary Materials, including the flowrates of streams between membrane units and streams to products.

Unit Membrane Type Inlet Flowrate (t/h)
In addition, without specifications for the permeate product, quite some ZnCl 2 is found in the waste stream (see Table 2). In fact, ZnCl 2 is more expensive than glucose. Therefore, specifications for the permeate product will be set in case 2 for ZnCl 2 recovery and waste reduction. Figure 5 shows the retentate concentrations and glucose/ZnCl2 ratios in MF an units. Other details of the NF membrane networks for scenarios 1-3 are given in S1-S3 in the Supplementary Materials, including the flowrates of streams between brane units and streams to products. In addition, without specifications for the permeate product, quite some Z found in the waste stream (see Table 2). In fact, ZnCl2 is more expensive than g Therefore, specifications for the permeate product will be set in case 2 for ZnCl2 re and waste reduction.

Case 2: Membrane Network Synthesis Using Hypothetical Parameters
Although the use of both NF1 and NF2 membranes in case 1 resulted in a significant reduction in capital cost, there were still a large number of membrane units (=12) in the network, due to low ZnCl 2 removal ratios (<50%) in most of the membrane units. In this case, an alternative membrane type with hypothetical performance is used. The parameters are shown in Table 4. Compared to NF1 and NF2, the hypothetical membrane (NF3) has lower glucose retention, a constant ZnCl 2 removal ratio of 70% (r ZnCl 2 ,j = 30%) and a higher capital cost for the same unit size. As in case 1, the flowrate of the feed mixture to unit MF is assumed to be 3 t/h, and the inlet flowrate of MR units is fixed at the base flowrate of 3 t/h. The NF membrane network is designed to meet the specifications for both retentate and permeate products. In addition to a minimum glucose/ZnCl 2 ratio of 10 and a minimum glucose recovery of 90% in the retentate product, in this case the permeate product is required to have a minimum ZnCl 2 /glucose ratio of 10 and a minimum ZnCl 2 recovery of 90%, as given in Equations (35) and (36).
Solving the MINLP model of Equation (26) subject to Equations (1)- (25) and (33)-(36) yields the results in Table 5. It can be seen that with the specifications for the permeate product, there is no longer a waste stream. Also, with a much higher ZnCl 2 removal ratio, the product specifications are met with only six membrane units, hence a considerable cost reduction of 30.2% compared to the solution for scenario 3 of case 1. While the retentate product just meets the glucose/ZnCl 2 ratio and glucose recovery specifications, the permeate product overshoots and achieves a ZnCl 2 /glucose ratio of 29.1 and a ZnCl 2 recovery of 97%. This is because not much glucose remains in the permeate product when the specifications for the retentate product are met. Figure 6 shows the corresponding NF membrane network, with some flowrates between membrane units given in Table 6.

Comparison and Discussion
It is worth mentioning that Kuo et al. [5] presented a multistage NF process for the removal of ZnCl2 from a bagasse hydrolysis solution. Through three NF membrane units in series, the concentration of ZnCl2 decreases by 74.8% and the glucose/ZnCl2 ratio increases from 0.245 to 1.268. There is no recovery of ZnCl2 in their design. In this work, NF networks were synthesized to increase the glucose/ZnCl2 ratio from 0.333 in the feed mixture to 10 in the retentate product stream. The proposed superstructure optimization approach allows the optimal system configuration to be determined systematically according to product specifications, resulting in a sophisticated system achieving improved performance. Moreover, the waste stream can be eliminated by recovering ZnCl2 in the permeate product stream.

Conclusions
A superstructure-based approach has been developed in this paper for the separation of hydrolysis solutions in bioethanol production using NF membranes. The optimization model for NF membrane network synthesis consists mainly of mass balances for individual membrane units and their interconnections, with the objective of minimizing the total capital cost of membrane units. A case study of a bagasse hydrolysis solution for glucose/ZnCl2 separation was solved to illustrate the proposed approach. First, experimental data for the performance of two types of NF membranes were used. The results show that using both membrane types gives higher ZnCl2 removal ratios and thus a better solution

Comparison and Discussion
It is worth mentioning that Kuo et al. [5] presented a multistage NF process for the removal of ZnCl 2 from a bagasse hydrolysis solution. Through three NF membrane units in series, the concentration of ZnCl 2 decreases by 74.8% and the glucose/ZnCl 2 ratio increases from 0.245 to 1.268. There is no recovery of ZnCl 2 in their design. In this work, NF networks were synthesized to increase the glucose/ZnCl 2 ratio from 0.333 in the feed mixture to 10 in the retentate product stream. The proposed superstructure optimization approach allows the optimal system configuration to be determined systematically according to product specifications, resulting in a sophisticated system achieving improved performance. Moreover, the waste stream can be eliminated by recovering ZnCl 2 in the permeate product stream.

Conclusions
A superstructure-based approach has been developed in this paper for the separation of hydrolysis solutions in bioethanol production using NF membranes. The optimization model for NF membrane network synthesis consists mainly of mass balances for individual membrane units and their interconnections, with the objective of minimizing the total capital cost of membrane units. A case study of a bagasse hydrolysis solution for glucose/ZnCl 2 separation was solved to illustrate the proposed approach. First, experimental data for the performance of two types of NF membranes were used. The results show that using both membrane types gives higher ZnCl 2 removal ratios and thus a better solution with significant reductions in the number of membrane units and capital cost. Hypothetical parameters were then used to analyze the effect of membrane performance on the optimal network configuration. It was found that despite the lower glucose retention, the higher ZnCl 2 removal ratio allows the product specifications to be achieved with much fewer membrane units at a considerably reduced total capital cost. This indicates that the economic performance of the NF membrane separation system can be improved by increasing the ZnCl 2 removal ratio. Such analysis is thus useful for debottlenecking the system. The present study focuses on structural optimization using a shortcut model for membrane units. Future work will involve detailed membrane models to optimize the operating conditions (e.g., pressure). A more sophisticated objective function will also be used to refine the membrane network. In addition, although developed for NF membrane network synthesis, the proposed mathematical model is generic in nature and can be extended to other separation problems.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A. Nomenclature
Notation used in the MINLP model is given below.
Appendix A.1. Indices and Sets i ∈ I components j ∈ J membrane units j ∈ J MP membrane units handling permeates k ∈ K membrane units handling retentates k ∈ K membrane types l ∈ L inlet ZnCl 2 concentration intervals