07.02.2013 Views

Optimization and Computational Fluid Dynamics - Department of ...

Optimization and Computational Fluid Dynamics - Department of ...

Optimization and Computational Fluid Dynamics - Department of ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

200 Laurent Dumas<br />

7.3.1.1 Genetic Algorithms (GA)<br />

GAs are global optimization methods directly inspired from the Darwinian<br />

theory <strong>of</strong> evolution <strong>of</strong> species [11]. They require following the evolution <strong>of</strong><br />

a certain number Np <strong>of</strong> possible solutions, also called population. A fitness<br />

value is associated to each element (or individual) xi ∈O<strong>of</strong> the population<br />

that is inversely proportional to J(xi) in case <strong>of</strong> a minimization problem. The<br />

population is regenerated Ng times by using three stochastic principles called<br />

selection, crossover <strong>and</strong> mutation, that mimic the biological law <strong>of</strong> “survival<br />

<strong>of</strong> the fittest”.<br />

The GA that will be used in the drag reduction problem in Sect. 7.4 acts<br />

in the following way: at each generation, Np/2 couples are selected by using<br />

a roulette wheel process with respective parts based on the fitness rank <strong>of</strong><br />

each individual in the population. To each selected couple, the crossover <strong>and</strong><br />

mutation principles are then successively applied with a respective probability<br />

pc <strong>and</strong> pm. The crossover <strong>of</strong> two elements consists in creating two new<br />

elements by doing a barycentric combination <strong>of</strong> them with r<strong>and</strong>om <strong>and</strong> independent<br />

coefficients in each coordinate. The mutation principle consists <strong>of</strong><br />

replacing a member <strong>of</strong> the population by a new one r<strong>and</strong>omly chosen in its<br />

neighborhood. A one-elitism principle is added in order to be sure to keep in<br />

the population the best element <strong>of</strong> the previous generation.<br />

Thus, the algorithm can be written as:<br />

• Choice <strong>of</strong> an initial population P1 = {x1 i ∈O, 1 ≤ i ≤ Np}<br />

• ng =1.Repeatuntilng = Ng<br />

• Evaluate {J(x ng<br />

i ), 1 ≤ i ≤ Np} <strong>and</strong> m =min{J(x ng<br />

i ), 1 ≤ i ≤ Np}<br />

• 1-elitism: if ng ≥ 2 & J(Xng−1)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!