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.

234 Marco Manzan, Enrico Nobile, Stefano Pieri <strong>and</strong> Francesco Pinto<br />

part centered on the wall. Therefore, the value <strong>of</strong> parameter x2 is equal to<br />

(L − Lwave)/2. The number <strong>of</strong> DOFs is now 10.<br />

8.5.1.3 Channel Construction<br />

In Figs. 8.5(a) <strong>and</strong> 8.5(b) both types <strong>of</strong> channel are depicted, i.e., linear piecewise<br />

channel <strong>and</strong> NURBS channel, <strong>and</strong> they are constructed in the same way.<br />

Once the lower wall pr<strong>of</strong>ile has been obtained, the upper one, as illustrated<br />

in Fig. 8.1, is made first by a simple translation in y-direction <strong>of</strong> the former<br />

in order to obtain the height <strong>of</strong> the channel, followed by a translation in<br />

x-direction. This guarantees the realistic desire to construct the channels <strong>of</strong><br />

e.g., a finned heat exchanger, by simple juxtaposition <strong>of</strong> identical wavy plates.<br />

This action introduces in both cases another DOF which defines the x translation<br />

<strong>of</strong> the upper pr<strong>of</strong>ile, variable transl in Tables 8.1 <strong>and</strong> 8.4. In order<br />

to avoid linear-dependent channels, the translation range is bounded to one<br />

half-period in both the positive <strong>and</strong> negative x-direction. The y-translation is<br />

instead fixed: in this way the average height <strong>of</strong> the channel is set to 0.5, that<br />

is half <strong>of</strong> the non-dimensional hydraulic diameter (Dh). Finally the periodic<br />

duct module is cut by two straight lines, representing the inlet <strong>and</strong> the outlet<br />

section, as sketched in Fig. 8.1. Overall, we are left with 5 DOFs for the linear<br />

piecewise model, <strong>and</strong> 11 DOFs for the NURBS channel.<br />

8.5.1.4 Three-dimensional Extension<br />

Once the two-dimensional results have been obtained, it has been tried to<br />

encourage further mixing in the flow, <strong>and</strong> therefore an augmentation <strong>of</strong> the<br />

heat transfer rate by forcing a non-zero value <strong>of</strong> the z-component <strong>of</strong> the velocity<br />

vector [47]. Three-dimensional analysis has been performed by simple<br />

extrusion from two-dimensional modules. For this purpose, the best twodimensional<br />

optimized channels in terms <strong>of</strong> low friction factor <strong>and</strong> high Nusselt<br />

number have been selected first, <strong>and</strong> one new additional variable was<br />

introduced such as the extrusion angle. The constructing procedure is described<br />

in Fig. 8.2. Due to the exploratory nature <strong>of</strong> this study, the extrusion<br />

length is fixed to the same value as the channel height.<br />

8.5.2 CC Module<br />

The CAD package CATIA [11] has been utilized to describe the heat transfer<br />

surface geometry, thanks to its key feature <strong>of</strong> allowing the generation <strong>of</strong><br />

parametric drawings. In Fig. 8.6, the parameters required to describe the geometry<br />

<strong>of</strong> a CC surface are presented. Points P1, P2, P3 <strong>and</strong> P4 define the

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

Saved successfully!

Ooh no, something went wrong!