Improvement of Multipopulation Genetic Algorithms Convergence Time Maria

Improvement of Multi-population Genetic Algorithms Convergence Time Maria Angelova, Tania Pencheva maria. angelova@clbme. bas. bg, tania. pencheva@clbme. bas. bg

Fermentation processes (FP) are widely used in different branches of industry – in the production of pharmaceuticals, chemicals and enzymes, yeast, foods and beverages. n n Fermentation processes are: characterized as complex, dynamic systems with interdependent and time-varying process variables; described by non-linear models with a very complex structure. An important step for adequate modeling of non-linear models of FP is the choice of a certain optimization procedure for model parameter identification.

Aims of the investigation n n The influence of five of the main genetic algorithm parameters to be investigated for six modifications of multipopulation genetic algorithms (Mp. GA) towards convergence time: generation gap - GGAP crossover rate - XOVR mutation rate - MUTR insertion rate - INSR migration rate - MIGR

Genetic algorithms (GA) : - are a direct random search technique for finding global optimal solution in complex multidimensional search space; - are based on mechanics of natural selection and natural genetics; - have advantages such as hard problems solving, noise tolerance, easy to interface and hybridize; - are proved to be very suitable for the optimization of highly non-linear problems, - are applied in the area of biotechnology, especially for parameter identification of fermentation process models.

Multi-population genetic algorithms n Simple genetic algorithm (SGA) works with a population of coded parameters set called “chromosomes”. Each of these artificial chromosomes is composed of binary strings (or genes) of certain length (number of binary digits). Each gene contains information for the corresponding parameter. n Multi-population genetic algorithm (Mp. GA) is a single population genetic algorithm, in which many populations, called subpopulations, evolve independently from each other for a certain number of generations. After a certain number of generations (isolation time), a number of individuals are distributed between the subpopulations.

Mp. GA modifications n n n Six kinds of Mp. GA are investigated towards improvement of algorithms convergence time. Mp. GA differ from each other in the sequence of execution of main genetic operators’ selection, crossover and mutation: Mp. GA-SCM (coming from sequence selection, crossover, mutation); Mp. GA-CMS (crossover, mutation, selection); Mp. GA-SMC (selection, mutation, crossover); Mp. GA-MCS (mutation, crossover, selection); Mp. GA-SC (selection, crossover); Mp. GA-CS (crossover, selection) is newly developed here, provoked by the promising results obtained when selection operator is processed after crossover in SGA.

Mp. GA-CS The main idea of this modification is that the individuals are reproduced processing only crossover and avoiding mutation. In the beginning, Mp. GA-CS generates a random population of n chromosomes, i. e. suitable solutions for the problem. In order to prevent the loss of reached good solution by crossover, selection has been processed after crossover. Parents’ genes combine to form a whole new chromosome during the crossover. After the reproduction, the Mp. GA-CS calculates the objective function for the offspring and the best fitted individuals from the offspring are selected to replace the parents, according to their objective values. When a certain number of generations is fulfilled, the Mp. GA-CS is terminated.

Range of investigated genetic algorithm parameters Very big generation gap value does not improve performance of GA, especially regarding how fast the solution will be found. Mutation is randomly applied with low probability, typically in the range 0. 01 and 0. 1. A higher crossover rate introduces new strings more quickly into the population. A low crossover rate may cause stagnation due to the lower exploration rate. Insertion rate is a general measure how many of the individuals produced at each population are inserted into the new generation. Migration rate characterized the number of exchanged individuals.

Mathematical model of S. cerevisiae fed-batch cultivation where X, S, E, O 2 and O 2* are concentrations of biomass, substrate (glucose), ethanol, [g. l-1], oxygen and dissolved oxygen saturation, [%]; F – feeding rate, [l. h-1]; V – volume of bioreactor, [l]; – volumetric oxygen transfer coefficient, [h-1]; Sin – glucose concentration in the feeding solution, [g. l-1]; , q. S, q. E and are respectively specific rates of growth, substrate utilization, ethanol production and dissolved oxygen consumption, [h-1].

Specific rates where – maximum growth rates of substrate and ethanol, [h-1]; k. S, k. E – saturation constants of substrate and ethanol, [g. l-1]; Yij – yield coefficients, [g. g-1]. Optimization criterion: where Y is the experimental data, Y* – model predicted data, Y = [X, S, E, O 2].

Influence of GGAP in Mp. GA with three genetic operators Influence of GGAP has been investigated towards model accuracy and convergence time.

Influence of GGAP in Mp. GA with two genetic operators Influence of GGAP has been again investigated towards model accuracy and convergence time.

Comparison of Mp. GA results n The optimization criterion values obtained with six kinds of Mp. GA are very similar - there is no loss of accuracy. The obtained results can be grouped: Mp. GA-SCM with Mp. GA-SMC and Mp. GA-CMS with Mp. GA-MCS, but the convergence time in second group is much bigger than the first group. n Two algorithms without mutation execution, Mp. GA-SC and Mp. GA-CS, can be grouped together too. In cases when algorithms are implemented only with two operators the calculation time is much less but for the expenses of model accuracy. n Proceeding selection operator before crossover and mutation (no matter their order) needs much less computational time at GGAP, XOVR, MUTR, MIGR and INSR.

Results concerning considered GA parameters n The GGAP is the most sensitive from five investigated parameters concerning the convergence time. Up to 40% (in case of Mp. GASCM, ) can be saved using GGAP = 0. 5 instead of 0. 9 without loss of accuracy. n Exploring different values of crossover rate no such time saving is realized but it should be pointed that values of 0. 85 for XOVR can be assumed as more appropriate. n Exploring MUTR values of 0. 02 can be assumed as more appropriate. n In INSR and MIGR no tendency of influence can be drawn.

Optimal GA parameter values GGAP = 0. 5, XOVR = 0. 85, MUTR = 0. 02, INSR = 0. 9 and MIGR = 0. 1. Because of the similarity of the results obtained with all six kinds of algorithms the results obtained by the developed here Mp. GA-CS, are presented. As a result of parameter identification, the values of model parameters are respectively: S = 0. 98 [h-1], E = 0. 13 [h-1], k. S = 0. 13 [g·l-1], k. E = 0. 84 [g·l-1], YSX = 0. 42 [g·g-1], YEX = 1. 67 [g·g-1], = 96. 2329 [h-1], YOS = 766. 7862 [g·g-1], YOE = 125. 5165 [g·g-1], while CPU time was 288. 6720 s and J = 0. 0221. Presented results from Mp. GA-CS application for parameter identification of S. cerevisiae fed-batch cultivation show the effectiveness of GA for solving complex nonlinear problems.

Experimental and model data for biomas and substrate concentration

Experimental and model data for ethanol and dissolved oxygen concentration

Analysis and conclusions Altogether six kinds of multi-population genetic algorithms have been examined: - Four of them are with exchanged operators’ sequence of selection, crossover and mutation operators; - Two modifications are without performing of mutation operator. The influence of some of genetic algorithm parameters, namely GGAP, XOVR, MUTR, INSR and MIGR, has been examined for all six kinds of genetic algorithms and the most sensitive - GGAP has been distinguished aiming to improve the convergence time. As “favorite” among the considered here algorithms Mp. GA-SCM has been marked as the fastest one. Up to almost 40% from calculation time can be saved in the case of Mp. GA-SCM application using GGAP = 0. 5 instead of 0. 9 without loss of model accuracy. All modifications of Mp. GA show the effectiveness of genetic algorithms for solving complex nonlinear problems.

IMACS’ 11 Improvement of Multi-population Genetic Algorithms Convergence Time ACKNOWLEDGEMENTS This work is partially supported by the European Social Fund and Bulgarian Ministry of Education, Youth and Science under Operative Program “Human Resources Development”, grant BG 051 PO 001 -3. 3. 04/40 and National Science Fund of Bulgaria, grant DID 02 -29 “Modeling Processes with Fixed Development Rules”. Thank you for your attention!