these simulation numerique et modelisation de l'ecoulement autour ...

these simulation numerique et modelisation de l'ecoulement autour ... these simulation numerique et modelisation de l'ecoulement autour ...

21.04.2014 Views

A. Analysis of the small-scale reference LES results Small-scale simulations were performed in order to learn more about the fine structure of the flow around and inside a perforated plate in an isothermal configuration corresponding to the LARA experiment. A complete description of the methodology can be found in Mendez et al. 30 In FCFC experiments, the flow is known to be different depending on the number of rows of perforations considered. This is very difficult to handle from a modeling point of view. Furthermore, as already said, this dependance on the row number is not obvious to transpose in full burners. It has been decided to deal with this difficulty by considering the case where the perforated plate is infinite. The computational domain can then be reduced to a small box containing only one perforation, using periodic boundary conditions in order to reproduce the periodicity of the staggered pattern (see Fig. 3). This is also consistent with the construction of a local model. CALCULATION DOMAIN a b PRIMARY FLOW SECONDARY FLOW INFINITE PERFORATED WALL 24 d 2 d 10 d y x z 11.68 d 6.74 d Figure 3. From the infinite plate to the ‘bi-periodic’ calculation domain. (a): Geometry of the infinite perforated wall. (b): Calculation domain centered on a perforation; the bold arrows correspond to the periodic directions. The dimensions of the computational domain are provided. The calculation grid contains 1,500,000 tetrahedral cells: fifteen points describe the diameter of the hole and on the average the first off-wall point is located 5 wall units apart from the wall. Typically the cells along the wall to be cooled and in the hole are sized to a height of 0.3mm (recall that the aperture diameter is 5 mm). As natural mechanisms that normally drive the flow in a non-periodic domain are absent in such a configuration, an appropriate method has to be set up to generate and sustain the flow: the main tangential flow at both sides of the plate is enforced by a constant source term added to the momentum equation, as is usually done in channel flow simulations. The effusion flow through the hole is sustained by imposing a uniform vertical mass flow rate at the bottom boundary of the domain. The source terms and the uniform vertical mass flow rate were tuned in order to reproduce the operating conditions given in section IIB. A 9 of 26

pressure drop of 41 Pa is effectively imposed in the simulation. The resulting bulk velocity in the hole is V j ≈ 5.67 m s −1 and the mass flow rate through the hole is q = 0.126 g s −1 . Note that the mass density in the hole is approximately ρ j ≈ 1.13 kg m −3 . These values of mass density and velocity, ρ j and V j , are used as reference for non-dimensional quantities. The simulations in such a periodic configuration have proved to provide results that reproduce very well the global structure of the flow observed in the LARA experiment and comparisons with experimental profiles at row 9 show good agreement. 30 The numerical fields have been averaged over 20 flow through times (FTT). This time-averaged solution of the flow is analyzed here in terms of modeling. Time-averaged quantities are denoted by the ¯. operator. The model has to reproduce correctly the momentum and energy fluxes at the perforated plate (injection and suction sides), at a given mass flow rate and a given geometry. Indeed, the mass flow rate through the plate is supposed to be known: wether it is imposed by the user, or it can be calculated thanks to a relation between the pressure drop and the mass flow rate thanks to a discharge coefficient, in cases where both sides of the plate are computed. The analysis of small-scale data can support the modeling effort by answering to two main questions: • Among the terms contributing to the wall fluxes, which are the dominant ones? • The mass flow rate being known, is it possible to model the dominant fluxes? Note that because the flow computed is isothermal, only information about the momentum fluxes at the perforated plate are relevant. Momentum fluxes are calculated over two planes located just above (for the injection side) and just below (for the suction side) the perforated plate (see Fig. 4). The perforated plate is considered as a boundary condition made of two parts: the hole surface and the solid surface. The normal to the total surface, taken in the outward direction from the fluid point of view, is noted ⃗n. INJECTION SIDE SUCTION SIDE ⃗n S h S h S s S s S W = S s + S h ⃗n S W = S s + S h Figure 4. Schematic of the planes where the fluxes are assessed. From the momentum conservation equation, both the viscous (τ ik , k = 1, 2, 3) and the inviscid (Pδ ik + ρV i V k , k = 1, 2, 3) terms contribute to the flux associated to the momentum 10 of 26

A. Analysis of the small-scale reference LES results<br />

Small-scale <strong>simulation</strong>s were performed in or<strong>de</strong>r to learn more about the fine structure of the<br />

flow around and insi<strong>de</strong> a perforated plate in an isothermal configuration corresponding to the<br />

LARA experiment. A compl<strong>et</strong>e <strong>de</strong>scription of the m<strong>et</strong>hodology can be found in Men<strong>de</strong>z <strong>et</strong><br />

al. 30 In FCFC experiments, the flow is known to be different <strong>de</strong>pending on the number of<br />

rows of perforations consi<strong>de</strong>red. This is very difficult to handle from a mo<strong>de</strong>ling point of<br />

view. Furthermore, as already said, this <strong>de</strong>pendance on the row number is not obvious to<br />

transpose in full burners. It has been <strong>de</strong>ci<strong>de</strong>d to <strong>de</strong>al with this difficulty by consi<strong>de</strong>ring the<br />

case where the perforated plate is infinite. The computational domain can then be reduced<br />

to a small box containing only one perforation, using periodic boundary conditions in or<strong>de</strong>r<br />

to reproduce the periodicity of the staggered pattern (see Fig. 3). This is also consistent<br />

with the construction of a local mo<strong>de</strong>l.<br />

CALCULATION<br />

DOMAIN<br />

a<br />

b<br />

PRIMARY FLOW<br />

SECONDARY FLOW<br />

INFINITE<br />

PERFORATED<br />

WALL<br />

24 d<br />

2 d<br />

10 d<br />

y<br />

x<br />

z<br />

11.68 d 6.74 d<br />

Figure 3. From the infinite plate to the ‘bi-periodic’ calculation domain. (a): Geom<strong>et</strong>ry<br />

of the infinite perforated wall. (b): Calculation domain centered on a perforation; the bold<br />

arrows correspond to the periodic directions. The dimensions of the computational domain<br />

are provi<strong>de</strong>d.<br />

The calculation grid contains 1,500,000 t<strong>et</strong>rahedral cells: fifteen points <strong>de</strong>scribe the diam<strong>et</strong>er<br />

of the hole and on the average the first off-wall point is located 5 wall units apart<br />

from the wall. Typically the cells along the wall to be cooled and in the hole are sized to a<br />

height of 0.3mm (recall that the aperture diam<strong>et</strong>er is 5 mm).<br />

As natural mechanisms that normally drive the flow in a non-periodic domain are absent<br />

in such a configuration, an appropriate m<strong>et</strong>hod has to be s<strong>et</strong> up to generate and sustain the<br />

flow: the main tangential flow at both si<strong>de</strong>s of the plate is enforced by a constant source<br />

term ad<strong>de</strong>d to the momentum equation, as is usually done in channel flow <strong>simulation</strong>s. The<br />

effusion flow through the hole is sustained by imposing a uniform vertical mass flow rate<br />

at the bottom boundary of the domain. The source terms and the uniform vertical mass<br />

flow rate were tuned in or<strong>de</strong>r to reproduce the operating conditions given in section IIB. A<br />

9 of 26

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

Saved successfully!

Ooh no, something went wrong!