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.

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

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

Saved successfully!

Ooh no, something went wrong!