Optimization and Computational Fluid Dynamics - Department of ...
Optimization and Computational Fluid Dynamics - Department of ...
Optimization and Computational Fluid Dynamics - Department of ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
218 Marco Manzan, Enrico Nobile, Stefano Pieri <strong>and</strong> Francesco Pinto<br />
Since there is no single optimum to be found, we use a multi-objective genetic<br />
algorithm <strong>and</strong> the so-called Pareto dominance concept.<br />
The results obtained are very encouraging, <strong>and</strong> the procedure described<br />
can be applied, in principle, to even more complex convective problems.<br />
8.1 Introduction<br />
The common approach when solving therm<strong>of</strong>luid problems numerically is to<br />
first prescribe the geometry, boundary conditions, <strong>and</strong> thermophysical properties<br />
<strong>and</strong> then solve the governing equations for velocity, pressure, temperature,<br />
turbulent kinetic energy, etc. Problems <strong>of</strong> this type are referred to<br />
here as analysis problems. In design practice, the engineer usually tries different<br />
geometries, chooses other materials with different properties, <strong>and</strong> so<br />
on until satisfactory performance is obtained. Such cut-<strong>and</strong>-try method relies<br />
on the experience <strong>and</strong> skill <strong>of</strong> the designer to obtain any improvement but<br />
optimal performance is rarely achieved. Furthermore, this simple approach<br />
becomes impractical when the number <strong>of</strong> design variables is large, when there<br />
is more than one objective, <strong>and</strong> when there are several constraints to be satisfied.<br />
Therefore, in such cases, it is convenient, if not m<strong>and</strong>atory, to adopt<br />
an optimization strategy. However, the integration <strong>of</strong> numerical optimization<br />
techniques as part <strong>of</strong> the design process is still not very common today,<br />
particularly for complex heat transfer problems. In order to proceed in a systematic<br />
way, a basic requirement in shape optimization is to define the shape<br />
<strong>of</strong> the system to be optimized in terms <strong>of</strong> (known) functions <strong>and</strong> (unknown)<br />
parameters. This task is typically accomplished by means <strong>of</strong> a parametric<br />
computer-aided design (CAD) system, <strong>and</strong> by making use <strong>of</strong> an optimization<br />
procedure, automatic variations <strong>of</strong> the parameters associated with the geometric<br />
model lead to the creation <strong>of</strong> a variety <strong>of</strong> feasible shapes, which are<br />
then subjected to numerical analysis.<br />
The aim <strong>of</strong> this article is to describe a general strategy for automatic,<br />
multi-objective shape optimization <strong>of</strong> heat exchanger modules. This represents<br />
a truly multi-objective optimization problem, since it is desired, from a<br />
design point <strong>of</strong> view, to maximize the heat transfer rate in order to, for example,<br />
reduce the volume <strong>of</strong> the equipment <strong>and</strong> to minimize the friction losses<br />
which are proportional to the pumping power required. These two goals are<br />
clearly conflicting, <strong>and</strong> therefore no single optimum can be found. For this<br />
reason we use a Multi-Objective Genetic Algorithm (MOGA), <strong>and</strong> the socalled<br />
Pareto dominance which allows us to obtain a design set rather than<br />
a single design.<br />
<strong>Optimization</strong> <strong>of</strong> two-dimensional wavy channels is obtained by means <strong>of</strong><br />
an unstructured finite-element (FE) solver, for a fluid <strong>of</strong> Pr<strong>and</strong>tl number<br />
Pr = 0.7, representative <strong>of</strong> air, assuming fully developed velocity <strong>and</strong> temperature<br />
fields, <strong>and</strong> steady laminar conditions. The geometry <strong>of</strong> the channels