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 ...
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
- Page 95 and 96: 4.3 Sensibilité des simulations au
- Page 97 and 98: Chapitre 5 Simulations numériques
- Page 99 and 100: 2 S. Mendez and F. Nicoud COMBUSTIO
- Page 101 and 102: 4 S. Mendez and F. Nicoud (c) The i
- Page 103 and 104: 6 S. Mendez and F. Nicoud CALCULATI
- Page 105 and 106: 8 S. Mendez and F. Nicoud Name Numb
- Page 107 and 108: 10 S. Mendez and F. Nicoud x/d z/d
- Page 109 and 110: 12 S. Mendez and F. Nicoud main, it
- Page 111 and 112: 14 S. Mendez and F. Nicoud 1.0 a 1.
- Page 113 and 114: 16 S. Mendez and F. Nicoud ROWS 1 3
- Page 115 and 116: 18 S. Mendez and F. Nicoud Figure 1
- Page 117 and 118: 20 S. Mendez and F. Nicoud 30 ◦ 3
- Page 119 and 120: 22 S. Mendez and F. Nicoud & Mahesh
- Page 121 and 122: 24 S. Mendez and F. Nicoud Figure 1
- Page 123 and 124: 26 S. Mendez and F. Nicoud Figure 1
- Page 125 and 126: 28 S. Mendez and F. Nicoud film, th
- Page 127 and 128: 30 S. Mendez and F. Nicoud Region t
- Page 129 and 130: 32 S. Mendez and F. Nicoud Expressi
- Page 131 and 132: 34 S. Mendez and F. Nicoud geometri
- Page 133 and 134: 36 S. Mendez and F. Nicoud (f) Two
- Page 135 and 136: 38 S. Mendez and F. Nicoud Mendez,
- Page 137 and 138: Chapitre 6 Un modèle adiabatique p
- Page 139 and 140: ρ Mass density, kg/m 3 P T V i τ
- Page 141 and 142: cence of the jets. This is particul
- Page 143 and 144: Simulations. The small-scale LES re
- Page 145: chosen to compare with numerical re
- Page 149 and 150: Region total plate hole solid wall
- Page 151 and 152: This equality is almost verified at
- Page 153 and 154: V U σ V inj jet cotan(α′ ) V U
- Page 155 and 156: PERIODIC Z INFLOW 1 WALL 24 OUTFLOW
- Page 157 and 158: part of the plate. Downstream of th
- Page 159 and 160: streamwise momentum flux. As a cons
- Page 161 and 162: 8 Fric, T. and Roshko, A., “Vorti
- Page 163 and 164: Discussion sur la modélisation de
- Page 165 and 166: 20 15 10 y/d 5 0 -5 -10 0.0 0.2 0.4
- Page 167 and 168: Chapitre 7 Simulations numériques
- Page 169 and 170: Primary flow (burnt gases) Secondar
- Page 171 and 172: Figure 7: Nusselt number over the s
- Page 173 and 174: the perforated plate. A Large-Eddy
- Page 175 and 176: Les tableaux de flux permettent de
- Page 177 and 178: Conclusion générale Dans cette é
- Page 179 and 180: Conclusion générale En ce qui con
- Page 181 and 182: Bibliographie AKSELVOLL, K. & MOIN,
- Page 183 and 184: BIBLIOGRAPHIE CORON, X. 2001 Étude
- Page 185 and 186: BIBLIOGRAPHIE HALE, C. A., PLESNIAK
- Page 187 and 188: BIBLIOGRAPHIE LIEUWEN, T. & YANG, V
- Page 189 and 190: BIBLIOGRAPHIE MOST, A. & BRUEL, P.
- Page 191 and 192: BIBLIOGRAPHIE SCHILDKNECHT, M., MIL
- Page 193: Annexes
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