Research Article  Open Access
Optimal Solutions of Multiproduct Batch Chemical Process Using Multiobjective Genetic Algorithm with Expert Decision System
Abstract
Optimal design problem are widely known by their multiple performance measures that are often competing with each other. In this paper, an optimal multiproduct batch chemical plant design is presented. The design is firstly formulated as a multiobjective optimization problem, to be solved using the well suited non dominating sorting genetic algorithm (NSGAII). The NSGAII have capability to achieve fine tuning of variables in determining a set of non dominating solutions distributed along the Pareto front in a single run of the algorithm. The NSGAII ability to identify a set of optimal solutions provides the decisionmaker DM with a complete picture of the optimal solution space to gain better and appropriate choices. Then an outranking with PROMETHEE II helps the decisionmaker to finalize the selection of a best compromise. The effectiveness of NSGAII method with multiojective optimization problem is illustrated through two carefully referenced examples.
1. Introduction
Batch processes are used in production of many lowvolume but highvalueadded products (such as speciality chemicals, health care, food, agrochemicals,โฆetc.) because of operation flexibility in todayโs marketdriven environment. Manufactory of these products generally involves multi step synthesis [1]. In addition, if two or more products require similar processing steps, the same set of equipment is considered for at least economical reason. A batch plant producing multiple products is categorized as either a multiproduct plant or a multipurpose plant. Multiproduct plants produce multiple products following a sequential similar recipe. In such a plant, all the products follow the same path through the process and only one product is manufactured at a time. Each step is carried out on single equipment or on several parallel equipment units. Processing of other products is carried out using the same equipment in successive production runs or campaigns. In a multipurpose plant, each product follows one or more distinct processing paths; so more than one product may be produced simultaneously in such plants. The present work is directed toward the optimal design problems of multiproduct batch plants.
In conventional optimal design of a multiproduct plant, production requirements of each product and a total production time for all products are available and specified. The number, the required volume, and size of parallel equipment units in each stage are then determined to minimize the investment. It should be emphasized that batch plantsdesign has been for long identified as a key problem in chemical engineering as reported in literature [2โ9]. Formulation of batch plant design generally involves mathematical programming methods, such as linear programming (LP), nonlinear programming (NLP), mixedinteger linear programming (MILP) or mixedinteger nonlinear programming (MINLP). Mathematical programming or different optimization techniques, such as branch and bound, heuristics, genetic algorithm, simulated annealing, are thoroughly used to derive optimal solutions.
However, in reality the multiproduct design problem can be formulated as a multiobjective design optimization problem in which one seeks to minimize investment, operation cost, and total production time, and, simultaneously, to maximize the revenue. Recall that not much work has been reported in the literature on the multiobjective optimal design of a multiproduct batch plant. Huang and Wang [10] introduced a fuzzy decisionmaking approach for multiobjective optimal design problem of a multiproduct batch plant. A monotonic increasing or decreasing membership function is used to define the degree of satisfaction for each objective function and the problem is then represented as an augmented minmax problem formulated as MINLP models. To obtain a unique solution, the MINLP problem is solved using a hybrid differential evolution technique. Dedieu et al. [11] presented the development of a twostage methodology for multiobjective batch plant design and retrofit according to multiple criteria. The authors used a multiobjective genetic algorithm based on the combination of a singleobjective genetic algorithm and a Pareto sort procedure for proposing several plant structures and a discrete event simulator for evaluating the technical feasibility of the proposed configurations.
In the case of multiple objectives, an optimum solution with respect to all objectives may not exist. In most cases, the objective functions are in conflict, because in order to decrease any of the objective functions, we need to increase other objective functions. Recently, Solimanpur et al. [12] developed a sophisticated multiobjective integer programming model where the objectives considered were the maximization of total similarity between parts, the minimization of the total processing cost, the minimization of the total processing time and the minimization of the total investment needed for the acquisition of machines [13].
The presence of multiple objectives in a problem usually gives rise to a family of nondominated solutions, largely known as Paretooptimal solutions, where each objective component of any solution along the Pareto front can only be improved by degrading at least one of its other objective components. Since none of the solutions in the nondominated set is absolutely better than any other, any one of them is then an acceptable solution. As it is difficult to choose any particular solution for a multiobjective optimization problem without iterative interaction with the decision maker (DM) [14] one general approach is to establish first the entire set of Paretooptimal solutions, where an external Decision Maker (DM) direct intervention gives interactive information in the multiobjective optimization loop [15]. So, a satisfactory solution of the problem is found as soon as the knowledge is acquired [16]. Promethee II (Preference Ranking Organisation METHod for Enrichment Evaluationsโ2nd version) is a popular decision method that has been successfully applied in the selection of the final solution of multiobjective optimization problems. It generates a ranking of available points, according to the DM preferences, and the best ranked one is considered the favourite final solution. It is based on the concept of outranking relation, which is a binary relation defined between every pair of alternatives, in such way that, if is preferred to (according to the DM interests), then it is said that outranks . When these relations are defined between all pairs of alternatives, they are exploited according to some rules in order to rank all solutions from the best to the worst.
The first GA proposed for multiobjective optimization was VEGA [17]. This is a nonPareto based approach based on the selection of several relevant groups of individuals, each group being associated to a given objective. It is reported that the method tends to crowd results at extremes of the solution space, often yielding to poor convergence of the Pareto front. A more recent algorithm, based on scalarization with a weighted sum function, is proposed in Ishibuchi and Murata [18] where the weights are randomly chosen. Many successful evolutionary multiobjective optimization algorithms were developed based on the two ideas suggested by Goldberg [19]: Pareto dominance and niching. Pareto dominance is used to exploit the search space in the direction of the Pareto front and niching technique explores the search space along the front to keep diversity. The wellknown algorithms in this category include Multiobjective Genetic Algorithm: (MOGA) [16], Niched Pareto Genetic Algorithm: (NPGA) [20], Strength Pareto Evolutionary Algorithm: (SPEA) [21], Multiobjective Evolutionary Algorithm: (MOEA) [22], the Nondominated Sorting Genetic Algorithm (NSGA) proposed by Srinivas and Deb [23] was one of the first evolutionary algorithm for solving multiobjective optimization problems. Although NSGA has been successfully applied, the main criticisms of this approach has been its high computational complexity of nondominated sorting, lack of elitism, and need for specifying a tuneable parameter called sharing parameter. Recently, Deb et al. [24] reported an improved version of NSGA, which they called NSGAII, to address all the above issues.
The purpose of this study is to extend this methodology for solution of multiobjective optimal control problems under the framework of NSGAII. The efficiency of the proposed method is illustrated by solving multiobjective optimization problem.
2. Formulation of the Multiobjective Problem
The problem of multiproduct batch plant covered in this paper can be defined by assuming that the plant consists of a sequence of batch processing stages that are used to manufacture different products. At each stage j there are identical units in parallel operating out of phase, each with a size . Each product follows the same general processing sequence.
Batches are transferred from one stage to the next without any delay, that is, we consider a zerowait operating policy.
In the conventional design of a multiproduct batch plant, one seeks to minimize the investment cost by determining the optimal number, required volume and size of parallel equipment units in each stage for a specified production requirement of each product and the total production time. However, in reality the designer considers not only minimizing the investment but also minimizing the operation cost and total production time while maximizing the revenue, simultaneously
So, the multiobjective problem consists of determining the following parameters:
(i) the number of parallel units in stage (ii) the required volume of a unit in stage ,(iii) size of the batch of product at the end of the stages,(iv) the cycle time for product , (v) the production requirement of product and,(vi) the total production time,while satisfying certain constraints such as volume, time, and so forth.
The constraints are expressed as follows:
(1) Volume constraints. Volume has to be able to process all the products :
(2) Time constraint. The summation of available production time for all products is not more than the net total time for production
(3) The limiting cycle time for product i:
(4) Dimension constraints. Every unit has restricted allowable range
3. Elitist Nondominated Sorting Genetic Algorithm (NSGAII)
The NSGA II Pareto ranking algorithm is an elitist Deb et al. [24] system and maintains an external archive of the Pareto solutions. In contrast to the simple genetic algorithms that look for the unique solution, the multiobjective genetic algorithm tries to find as many elements of the Pareto set as possible. For the case of the NSGAII, this one is provided with operators who allow it to know the level of nondominance of every solution as well as the grade of closeness with other solutions; which allows it to explore widely inside the feasible region.
In a brief form, the functioning of the multiobjective genetic algorithm NSGAII can be described through the following steps.
Fast Nondominated Sort
A very efficient procedure, is used to arrange the solutions in fronts (nondominated arranging), in accordance with their aptitude values. This is achieved, creating two entities for each of the solutions. A domination count the number of solutions which dominates the solution , and a set (), that contains the solutions that are dominated for . The solutions of the first front have the higher status of nondominance in the Pareto sense.
Diversity Preservation
This is achieved, by means of the calculation of the crowding degree or closeness for each of the solutions inside the population. This quantity is obtained, by calculating the average distance of two points on either side of a particular solution along each of the objectives. This quantity serves as an estimate of the cuboid perimeter, formed by using the nearest neighbours as the vertices. There is also, an operator called CrowdedComparison (), which guides to the genetic algorithm, towards the Pareto optimal front, in accordance with the following criterion:
In accordance with the previous criterion, between two nondominated solutions, we prefer the solution with the better rank. Otherwise, if both solutions belong to the same front, then, we prefer the solution that is located in a lesser crowded region.
Initial Loop
Initially, a random parent population () of size is created. Later this one is ordained, using the procedure of nondominated arranging. Then the usual binary tournament selection, recombination and mutation operators are used to create a new population (), of size .
Main Loop
The NSGAII procedure can be explained, by describing the th generation just as it is showed in Figure 1. The procedure begins with the combination of and forming a new population called , then the population is sorted using the nondomination criterion. Since all previous and current population members are included in , elitism is ensure. The population has a size of , later, the different fronts of nondominated solutions are created, being the front that contains the better rank solutions. Figure 4 shows that, during the process of forming the new population , the algorithm takes all members of the fronts and , and some elements of the front ; this is, because N solutions are needed exactly for the new population to find them exactly solutions, the last front is ordained, which for this description is the number 3, arranging the solutions in descending order by means of the crowded comparison (), and selecting the best solutions needed to fill all population slots. After having the population , the genetic operators of selection, crossing and mutation, are used to create the new population of size . Finally it is mentioned that the selection process, the crowded comparison operator is used.
4. Description of the Process
To demonstrate the effectiveness of NSGAII on batch plant processes two examples are given here. The first example is about a batch plant consisting of three processing stages (mixer, reactor, and centrifuge) to manufacture two products, A and B. The second example treats four processing stages (mixer, reactor, extraction and centrifuge) to manufacture three products A, B and C. The data for examples 1 and 2 are illustrated, respectively, in Tables 1 and 2 (the processing times, size factor for the units and cost for each product).


5. Results and Discussion
5.1. Example 1
A fourobjective optimization problem is considered and expressed in (1)โ(4). The set of decision variables consists of the batch size, the total production time, the number of parallel units at each stage, the cycle time for each product, and the required volume of a unit in each stage. Since the number of parallel units at each stage is an integer decision variable, we code this variable as a binary variable. All other decision variables are coded as real numbers. Thus, there are 3 integer variables and 10 real variables. In addition to the constraints expressed by (5)โ(8), we consider bounds on objective functions as additional constraints to generate feasible nondominated solutions in the range desired by the decisionmaker, to have 19 constraints in all
Then NSGAII is employed to solve the optimization problem with the following parameters: maximum number of generation up to 200, population size 500, probability of crossover 0.85, probability of mutation 0.05, distribution index for the simulated crossover operation 10 and distribution index for the simulated mutation operation 20.
The Paretooptimal solutions for example 1 are presented in Figure 2. The revenue () increases with the increase in operation cost (), while the investment cost () decreases. When all the four objective functions are considered simultaneously, solutions obtained in the present study show improvement as by Huang and Wang [10] results for the same problem. For example, let us consider the solution presented by Huang and Wang [10] with unit reference membership level for all objectives: , , . The solution (1) presented in Table 3 of the present study improves the above solution while is comparable.
 
Bounds for objective function: ,, , . 
(a)
(b)
Figure 3 presents the relationships between some chosen decision variables. The large set of multiple optimal solutions provides the decision maker with immediate information about the relationship among the several objective criteria and a set of feasible solutions. Thus, it helps the decisionmaker to select a highly confident choice of solution. The fixed optimal plant structure as 221 corresponds to a two mixers, two reactors, and one centrifuge design. The optimal solution is shown in Figure 4.
(a)
(b)
(c)
5.2. Example 2
The batch plant in this case consists of four processing stages to manufacture three products A, B, and C, with fourobjective optimization problem as expressed in (1)โ(4). The set of decision variables remains the same as that in example 1. But we deal with 4 integer variables, 14 real variables, and 31 constraints which includes bounds on objective functions.
The same model equations of example 1 are used here except the processing time, in (7). The time required to process one batch of product i in stage j is expressed as:
where , and are constants and is the batch size for product .
Thus, the processing time is not a constant, but depends on the decision parameters of the batch size. Table 2 presents the necessary data for the problem.
The constrained multiobjective MINLP problem is solved by NSGAII with the same set of NSGAII parameters as used in example 1.
As in example 1, the revenue () increases as operation cost () increases, while the investment cost () decreases following operation cost (). The Paretooptimal solutions for example 2 are presented in Figure 5. The relationships of the various decision variables are shown in Figure 6. Let mention that when all four objective functions are considered simultaneously, the solutions obtained in the present study improve significantly the results presented by Huang and Wang [10] for the same problem. For example, the solution presented by Huang and Wang [10] with unit reference membership level for all objectives the solution (1) presented in Table 4 of the present study improves the above solution while is comparable.
 
Bounds for objective function: , ,, 
(a)
(b)
(a)
(b)
(c)
(d)
In this example, the plant structure evolved as optimal is: two mixers, two reactors, two extractors, and one centrifuge as presented in Figure 7.
Implementation of a tradeoff analysis is dependent upon the availability of the decisionmakerโs preferences.
6. Conclusion
A multiobjective decision in a batch plant process design is considered and a non dominating sorting genetic algorithm (NSGAII) is developed to get an optimal zone containing solutions under the concept of Pareto set. NSGAII capability has been proved in evolving the entire set of nondominating solutions along the Pareto front in a single run of the algorithm. Thus, the Decision Maker (DM) is provided with the best tradeoff operating zone. Furthermore, a better confident choice of design among several compromises of the decision maker can be achieved if the decision variables effects on the objective functions are known.
Finally, the large set of solutions presents a useful base for further alternative approaches to fulfil the DM targets.
The inherent dynamic nature of batch processes allows for their ability to handle variations in feedstock and product specifications and provides the flexibility required for multiproduct or multipurpose facilities. They are thus best suited for the manufacture of lowvolume, highvalue products, such as specialty chemicals, pharmaceuticals, agricultural, food, and consumer products, and most recently the constantly growing spectrum of biotechnologyenabled products. Reduced time to market, lower production costs, and improved flexibility are all critical success factors for batch processes.
Nomenclature
Size of the batch of product at the end of the stages (kg)  
Operation cost in stage ($)  
Operation cost of product to be produced ($/kg)  
Price of product ($/kg)  
Total production time ()  
Number of stages in the batch process  
Number of products to be produced  
Number of parallel units in stage  
Production requirement of product (kg)  
Size factor of product in stage (L/kg)  
Cycle time for product ()  
Required volume of a unit in stage (L)  
โ  Greek Symbols 
Cost coefficient for unit  
Cost exponent for unit  
Processing time of product in stage (). 
References
 A. Lamghabbar, S. Yacout, and M. S. Ouali, โConcurrent optimization of the design and manufacturing stages of product development,โ International Journal of Production Research, vol. 42, no. 21, pp. 4495โ4512, 2004. View at: Publisher Site  Google Scholar
 D.M. Cao and X.G. Yuan, โOptimal design of batch plants with uncertain demands considering switch over of operating modes of parallel units,โ Industrial and Engineering Chemistry Research, vol. 41, no. 18, pp. 4616โ4625, 2002. View at: Publisher Site  Google Scholar
 W. Chunfeng and Z. Xin, โAnts foraging mechanism in the design of multiproduct batch chemical process,โ Industrial and Engineering Chemistry Research, vol. 41, no. 26, pp. 6678โ6686, 2002. View at: Publisher Site  Google Scholar
 H. D. Goel, M. P. C. Weijnen, and J. Grievink, โOptimal reliable retrofit design of multiproduct batch plants,โ Industrial and Engineering Chemistry Research, vol. 43, no. 14, pp. 3799โ3811, 2004. View at: Publisher Site  Google Scholar
 S.K. Heo, K.H. Lee, H.K. Lee, I.B. Lee, and J. H. Park, โA new algorithm for cyclic scheduling and design of multipurpose batch plants,โ Industrial and Engineering Chemistry Research, vol. 42, no. 4, pp. 836โ846, 2003. View at: Publisher Site  Google Scholar
 J. M. Montagna, A. R. Vecchietti, O. A. Iribarren, J. M. Pinto, and J. A. Asenjo, โOptimal design of protein production plants with time and size factor process models,โ Biotechnology Progress, vol. 16, no. 2, pp. 228โ237, 2000. View at: Publisher Site  Google Scholar
 A. Chakraborty, A. Malcolm, R. D. Colberg, and A. A. Linninger, โOptimal waste reduction and investment planning under uncertainty,โ Computers & Chemical Engineering, vol. 28, no. 67, pp. 1145โ1156, 2004. View at: Publisher Site  Google Scholar
 L. Cavin, U. Fischer, F. Glover, and K. Hungerbühler, โMultiobjective process design in multipurpose batch plants using a Tabu Search optimization algorithm,โ Computers & Chemical Engineering, vol. 28, no. 4, pp. 459โ478, 2004. View at: Publisher Site  Google Scholar
 T. Pinto, A. P. F. D. BarbosaPóvoa, and A. Q. Novais, โOptimal design and retrofit of batch plants with a periodic mode of operation,โ Computers & Chemical Engineering, vol. 29, no. 6, pp. 1293โ1303, 2005. View at: Publisher Site  Google Scholar
 H.J. Huang and F.S. Wang, โFuzzy decisionmaking design of chemical plant using mixedinteger hybrid differential evolution,โ Computers & Chemical Engineering, vol. 26, no. 12, pp. 1649โ1660, 2002. View at: Publisher Site  Google Scholar
 S. Dedieu, L. Pibouleau, C. AzzaroPantel, and S. Domenech, โDesign and retrofit of multiobjective batch plants via a multicriteria genetic algorithm,โ Computers & Chemical Engineering, vol. 27, no. 12, pp. 1723โ1740, 2003. View at: Publisher Site  Google Scholar
 M. Solimanpur, P. Vrat, and R. Shankar, โA multiobjective genetic algorithm approach to the design of cellular manufacturing systems,โ International Journal of Production Research, vol. 42, no. 7, pp. 1419โ1441, 2004. View at: Publisher Site  Google Scholar
 C. Dimopoulos, โMultiobjective optimization of manufacturing cell design,โ International Journal of Production Research, vol. 44, no. 22, pp. 4855โ4875, 2006. View at: Publisher Site  Google Scholar
 J.P. Brans and B. Mareschal, โThe PROMCALC and GAIA decision support system for multicriteria decision aid,โ Decision Support Systems, vol. 12, no. 45, pp. 297โ310, 1994. View at: Publisher Site  Google Scholar
 V. Bhaskar, S. K. Gupta, and A. K. Ray, โApplications of multiobjective optimization in chemical engineering,โ Reviews in Chemical Engineering, vol. 16, no. 1, pp. 1โ54, 2000. View at: Google Scholar
 C. Fonseca and P. Fleming, โAn overview of evolutionary algorithms in multiobjective optimization,โ Evolutionary Computation, vol. 3, no. 1, pp. 1โ16, 1993. View at: Publisher Site  Google Scholar
 J. Schaffer, โMultiple objective optimization with vector evaluated genetic algorithms,โ in Proceedings of the 1st International Conference on Genetic Algorithms and Their Applications (ICGA '85), pp. 93โ100, Pittsburgh, Pa, USA, July 1985. View at: Google Scholar
 H. Ishibuchi and T. Murata, โA multiobjective genetic local search algorithm and its application to flowshop scheduling,โ IEEE Transactions on Systems, Man and Cybernetics, Part C, vol. 28, no. 3, pp. 392โ403, 1998. View at: Publisher Site  Google Scholar
 D. E. Goldberg, Genetic Algorithms in Search, Optimization and Machine Learning, AddisonWesley, Reading, Mass, USA, 1989.
 J. Horn, N. Nafpliotis, and D. E. Goldberg, โA niched Pareto genetic algorithm for multiobjective optimization,โ in Proceedings of the 1st IEEE Conference on Evolutionary Computation, IEEE World Congress on Computational Intelligence (ICEC '94), vol. 1, pp. 82โ87, Orlando, Fla, USA, June 1994. View at: Publisher Site  Google Scholar
 E. Zitzler and L. Thiele, โMultiobjective evolutionary algorithms: a comparative case study and the strength Pareto approach,โ IEEE Transactions on Evolutionary Computation, vol. 3, no. 4, pp. 257โ271, 1999. View at: Publisher Site  Google Scholar
 K. C. Tan, T. H. Lee, and E. F. Khor, โEvolutionary algorithms with dynamic population size and local exploration for multiobjective optimization,โ IEEE Transactions on Evolutionary Computation, vol. 5, no. 6, pp. 565โ588, 2001. View at: Publisher Site  Google Scholar
 N. Srinivas and K. Deb, โMuiltiobjective optimization using nondominated sorting in genetic algorithms,โ Evolutionary Computation, vol. 2, no. 3, pp. 221โ248, 1994. View at: Publisher Site  Google Scholar
 K. Deb, A. Pratap, S. Agarwal, and T. Meyarivan, โA fast and elitist multiobjective genetic algorithm: NSGAII,โ IEEE Transactions on Evolutionary Computation, vol. 6, no. 2, pp. 182โ197, 2002. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2009 Diab Mokeddem and Abdelhafid Khellaf. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.