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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

8 Multi-objective <strong>Optimization</strong> in Convective Heat Transfer 233<br />

Table 8.1 Variables defining the linear piecewise channel<br />

Variable Symbols Range<br />

module length L [0.8;2.0]<br />

corrugation height h [0.0;0.5]<br />

forward edge angle ϕin [10 ◦ ;60 ◦ ]<br />

backward edge angle ϕout [10 ◦ ;60 ◦ ]<br />

translation <strong>of</strong> upper wall transl [−0.5L;0.5L]<br />

Table 8.2 Control points for the NURBS-pr<strong>of</strong>ile <strong>and</strong> parameters required<br />

Point x y DOFs<br />

1 0 0 0<br />

2 (L − Lwave)/2 0 0<br />

3 x2 + ∆x 23 0 1<br />

4 x5 + ρ 54 cos ϑ54 y5 + ρ 54 sinϑ54 2<br />

5 x1 + ∆x5 y1 + ∆y5 2<br />

6 x5 + ρ 56 cos ϑ56 y5 + ρ 56 sinϑ56 2<br />

7 x8 − ∆x 87 0 1<br />

8 x2 + Lwave 0 1<br />

9 L 0 1<br />

8.5.1.2 2D NURBS Parametrization<br />

In order to generate the NURBS channel during the optimization process,<br />

let’s start by defining first the wall pr<strong>of</strong>ile. As a good compromise between the<br />

number <strong>of</strong> DOFs <strong>and</strong> the geometrical complexity, a 9 control-point periodic<br />

cubic NURBS has been chosen, to ensure the periodicity <strong>of</strong> the channel itself.<br />

A large number <strong>of</strong> DOFs allows to describe minutely the pr<strong>of</strong>ile, but if this<br />

number is excessive, it would make the optimization process quite difficult<br />

<strong>and</strong> expensive. For this reason, it has been also decided to fix to a unitary<br />

value the curve’s weights. This, together with the uniform knots distribution<br />

we adopted, makes our NURBS curve practically equivalent to a B-Spline<br />

curve [18, 38].<br />

As depicted in Fig. 8.5(d), both the first <strong>and</strong> the last three aligned control<br />

points are needed to maintain the entrance <strong>and</strong> the exit <strong>of</strong> the pr<strong>of</strong>ile parallel<br />

to the x direction. The remaining ones give freedom to the wavy section. The<br />

parameters required for the definition <strong>of</strong> the lower pr<strong>of</strong>ile are explained in<br />

Table 8.2 <strong>and</strong> numbered as in Fig. 8.5(d). The symbol ∆ means the difference<br />

between the coordinates <strong>of</strong> two points while ρ ij is the module <strong>and</strong> ϑij the<br />

phase <strong>of</strong> a polar-coordinate system centered on the point i. The phase is<br />

positive counterclockwise with respect to the positive direction <strong>of</strong> the x axis.<br />

The pr<strong>of</strong>ile is again constructed, for convenience, in order to have the wavy

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

Saved successfully!

Ooh no, something went wrong!