UFR 4-18 Evaluation: Difference between revisions

From KBwiki
Jump to navigation Jump to search
m (Dave.Ellacott moved page SilverP:UFR 4-18 Evaluation to UFR 4-18 Evaluation)
 
(95 intermediate revisions by 3 users not shown)
Line 1: Line 1:
= Flow and heat transfer in a pin-fin array =
= Flow and heat transfer in a pin-fin array =
{{UFRHeader
{{UFRHeader
Line 11: Line 10:
<!--{{LoremIpsum}}-->
<!--{{LoremIpsum}}-->
== Comparison of CFD Calculations with Experiments ==
== Comparison of CFD Calculations with Experiments ==
{{Demo_UFR_Eval}}
<!--{{Demo_UFR_Eval}}-->
<br/>
<br/>
=== Global comparisons ===
=== Global comparisons ===
Line 17: Line 16:
==== Pressure loss coefficient ====
==== Pressure loss coefficient ====


The following table gives the numerical values of the presure loss coefficient in Ames et al. experiment as well as those obtained in the main computations.
The following table gives the numerical values of the presure loss coefficient in Ames et al. experiment as well as those obtained in the main computations. Note that the experimental data shown here have been sent by the authors of the experiment and are different from the ones found in Ames et al. (2005).  


Figure ?? shows these values and compare them to existing correlations from the literature:
Figure 25 shows these values and compares them to existing correlations from the literature (note that Ames et al. experimental values are very close to Metzger et al. correlation):


* Metzger et al. (1982) :
* Metzger et al. (1982) :
Line 29: Line 28:
<math>f = (0.25 + 0.1175/(P/D-1)^{1.08})  Re^{-0.16}</math>
<math>f = (0.25 + 0.1175/(P/D-1)^{1.08})  Re^{-0.16}</math>


Note that the experimental results have an error of 2.5% at the higher Reynolds numbers.
Table 5 compares various comptations to the experiment and gives the relative errors.
Large Eddy Simulation shows a very good agreement with the correlations and Ames et al. experimental results. The relative error is equal to 1% end 3% for the two Reynolds numbers, respectively. for both computed Reynolds numbers.  The EB-RSM shows a very satisfactory pressure loss coefficient at <math>Re_D=10000</math> (1% error) but a worse pressure loss coefficient at <math>Re_D=30000</math> (10% error). The k&omega;-SST shows globally good results (7% and 6% errors) and the &phi;-model poorer results with an error of 17% and 12% for the two Reynolds numbers, respectively.  
Large Eddy Simulation shows a very good agreement with the correlations and Ames et al. experimental results. The relative error is equal to 1% end 3% for the two computed Reynolds numbers, respectively.  The EB-RSM model shows a very satisfactory pressure loss coefficient at <math>Re_D=10000</math> (1% error) but the result is worse at <math>Re_D=30000</math> (10% error). The k&omega;-SST model shows globally good results (7% and 6% errors) and the &phi;-model poorer results with an error of 17% and 12% for the two Reynolds numbers, respectively.  


Wall resolved LES is then the best candidate to obtain almost perfect pressure loss coefficients in the present configuration. The EB-RSM five very good results at the lowest Reynolds number but a substantial error is observed at <math>Re_D=30000</math>. Although the k&omega;-SST exhibits pretty good results, the physics of the present flow is not well captured and one can't rely on this model. The &phi;-model has to be avoided for the present configuration in order to compute the pressure loss coefficient.  
Wall resolved LES is therefore a good candidate to obtain satisfactory pressure loss coefficients in the present configuration. The EB-RSM model gives very good results at the lowest Reynolds number but a substantial error is observed at <math>Re_D=30000</math>. Although the k&omega;-SST model exhibits pretty good results, the physics of the present flow is not well captured as it will be shown later and one can't rely on this model. The &phi;-model gives the worst results and is not reliable.  


{| border="1"
{| border="1"
|+Pressure loss coefficient values for the different computations
|+ '''Table 5:''' Pressure loss coefficient values for the different computations compared to Ames et al. experimental values.
!  
!  
! Ames et al. <math>Re_D=3000</math>
! Ames et al. <math>Re_D=3000</math>
Line 75: Line 74:
| 6
| 6
|}
|}




[[Image:Friction_Final_Results.jpg]]<br/>
[[Image:Friction_Final_Results.jpg]]<br/>
'''Figure ??:''' Pressure loss coefficient
'''Figure 25:''' Pressure loss coefficient


==== Average Nusselt Number ====
==== Average Nusselt Number ====


The following table gives the numerical values of the average Nusselt numbers in Ames et al. experiment as well as those obtained in the main computations.
Table 6 gives the numerical values of the average Nusselt numbers on the bottom wall in Ames et al. experiment as well as those obtained in the main computations (go to "Test Case Studies" section to find the definition of the Nusselt number and its average on the bottom wall).  


Figure ?? shows these values and compare them to existing correlations from the literature:
Figure 26 shows these values and compares them to existing correlations from the literature (Note again that the Ames et al. experimental results are close to the Metzger correlation):


* Metzger et al. (1982) :
* Metzger et al. (1982) :
Line 96: Line 93:
<math>Nus = 0.153*Re^{0.686}</math>
<math>Nus = 0.153*Re^{0.686}</math>


Note that the experimental results have an error of +/- 12%, +/-11.4%, and +/-10.5% for the 3 000, 10 000, and 30 000 Reynolds numbers, respectively, in the endwall regions adjacent to the pins and +/- 9% away from the pin.
Note that the local experimental results have an error of +/- 12%, +/-11.4%, and +/-10.5% for the 3 000, 10 000, and 30 000 Reynolds numbers, respectively, in the endwall regions adjacent to the pins and +/- 9% away from the pin.
 
Table 6 compares the computations with the experiment and gives the relative errors. One can see that the two models using an eddy viscosity (&phi;-model and k&omega;-SST) are too far from the experimental results, except the &phi;-model at <math>Re_D=30000</math>. However, as it will be shown later, this model exhibits wrong physics. LES and EB-RSM yield very satisfactory values in the error interval. The superiority of these two models will be discussed later when dealing with the local normalized Nusselt number on the bottom wall. Ames and Dvorak (2006a) used standard, Realizable k-&epsilon; and RNG k-&epsilon; models. They showed a substantial underestimation of the average Nusselt number when the pins were heated.
Delibra et al. (2009) used the URANS &zeta;-f model (very close to the &phi;-model utilized in the present work) at the 10000 and 30000 Reynolds numbers. They obtained average Nusselt numbers of 46.2 and 122.3, respectively. These values are further from the experimental ones. The temperature was imposed at the wall in Delibra et al. (2009). Imposing the temperature at the wall in the present work with the &phi;-model led to average Nusselt numbers of 41 and 100, respectively. This shows that imposing the temperature gives higher Nusselt numbers than the ones obtained by imposing the heat flux with this model. The effect of the temperature boundary condition has also been tested with the LES at <math>Re_D=10000</math> and the EB-RSM model at <math>Re_D=30000</math>. The average Nusselt numbers were equal to 46 and 115.2, respectively. No clear conclusion can be drawn for these two models which give a more realistic physics than the ones using an eddy viscosity concerning the effect of imposing the temperature instead of the heat flux at the bottom wall. However, one can state that there is a difference between the two methods.
Delibra et al. (2008) also reported LES computations. The one carried out for <math>Re_D=30000</math> was too coarse with a short time integration and will not be considered here. LES at <math>Re_D=10000</math> gave an average Nusselt number equal to 44.3 which is lower than the one obtained in the present work even with a fixed temperature at the wall. The error compared to the experimental value is equal to 18% which is too high. This discrepancy can be attributed to two factors:
 
* the mesh was probably still too coarse, not in the wall normal directions but in the stream-wise and span-wise ones (the mesh in the present study contains 76 million cells and their mesh contained 5.5 million cells),
 
* the boundary condition using an imposed temperature is less representative of the experiment.
 
Finally, Delibra et al. (2010) switched to use a hybrid RANS/LES approach to overcome the understimation of the Nusselt number but the results where not much better than the ones obtained with the URANS approach based on the &zeta;-f model.


{| border="1"
{| border="1"
|+Average Nusselt number on the bottom wall
|+'''Table 6:''' Average Nusselt number on the bottom wall
! Ames et al. 3000
!
! Ames et al. 10000
! Ames3000
! Ames et al. 30000
! LES3000
! LES3000
! Ames10000
! LES10000
! LES10000
! DLES10000(*)
! EBRSM10000
! EBRSM10000
! PHI10000
! PHI10000
! kwSST10000
! kwSST10000
! DURANS10000(*)
! Ames30000
! EBRSM30000
! EBRSM30000
! PHI30000
! PHI30000
! kwSST30000
! kwSST30000
! DURANS30000(*)
|-
|-
| Value
| 22.1
| 22.1
| 21.6
| 54.1
| 54.1
| 111.5
| 21.6
| 48.6
| 48.6
| 44.1
| 48.1
| 48.1
| 37.6
| 37.6
| 33.1
| 33.1
| 46.2
| 111.5
| 114.8
| 114.8
| 93.3
| 93.3
| 72
| 72
| 122.3
|-
| Error (%)
| -
| 2
| -
| 10
| 18
| 11
| 30
| 39
| 14
| -
| 3
| 16
| 35
| 9
|}
|}
 
(*) Delibra et al. results.


[[Image:NusAv_Final_Results.jpg]]<br/>
[[Image:NusAv_Final_Results.jpg]]<br/>
'''Figure ??:''' Averaged Nusselt number on the bottom wall
'''Figure 26:''' Average Nusselt number on the bottom wall


=== Local comparisons ===
=== Local comparisons ===


==== Mean pressure coefficient ====
==== Mean pressure coefficient ====
Figures 27 to 29 give the profiles of the pressure coefficient along the midline of pins 1, 2, 3 and 5 for the three Reynolds numbers.
Note that the experimental results have an uncertainty of +/-0.075 at <math>Re_D=3000</math> and of +/-0.025 at <math>Re_D=10000</math> and <math>Re_D=30000</math>.
LES results are in relatively good agreement with the experimental data. EB-RSM exhibits decent results as well. The two eddy viscosity models show poor results in particular the &phi;-model. The prediction of the <math>C_p</math> coefficient may be related to the wake prediction downstream of the pins.


[[Image:Cp_midplane_Re3000_rows1_2_3_5.jpg]]<br/>
[[Image:Cp_midplane_Re3000_rows1_2_3_5.jpg]]<br/>
'''Figure ??:''' Pressure coefficient along the midline of the pins - LES at <math>Re_D = 3000</math>
'''Figure 27:''' Pressure coefficient along the midline of the pins - LES at <math>Re_D = 3000</math>


[[Image:Cp_midplane_Re10000_rows1_2_3_5.jpg]]<br/>
[[Image:Cp_midplane_Re10000_rows1_2_3_5.jpg]]<br/>
'''Figure ??:''' Pressure coefficient along the midline of the pins - LES and URANS computations at <math>Re_D = 10000</math>
'''Figure 28:''' Pressure coefficient along the midline of the pins - LES and URANS computations at <math>Re_D = 10000</math>


[[Image:Cp_midplane_Re30000_rows1_2_3_5.jpg]]<br/>
[[Image:Cp_midplane_Re30000_rows1_2_3_5.jpg]]<br/>
'''Figure ??:''' Pressure coefficient along the midline of the pins - URANS computations at <math>Re_D = 30000</math>
'''Figure 29:''' Pressure coefficient along the midline of the pins - URANS computations at <math>Re_D = 30000</math>


==== Mean velocity profiles ====
==== Mean velocity profiles and contours ====
 
Figures 30 to 32 show the mean axial velocity component along line B for the different computations. The two eddy viscosity models exhibit a dramatic deficit of the mean axial velocity. On the other hand, LES and EB-RSM perform well compared to the experimental data with a slight overestimation of the mean stream-wise velocity peaks close to the walls in several locations, in particular departing from row 3. Note that the convergence study shown in the "Test Case Studies" section proved that the needed refinement to obtain a convergent solution around pin 2 is very important (with the first two levels of refinement, a deficit in the central velocity was observed).
 
Figure 33 shows the means axial velocity component along line A1 for the highest Reynolds number. One can see again that the linear eddy viscosity models are unable to predict the present flow. The results obtained with the EB-RSM model are much closer to the experimental data. However, the flow redistribution doesn't seem to be similar to the one obtained experimentally. This may come from the periodic boundary conditions used in the present work v.s. confined conditions in the experiment.
 
Figures 34 and 35 display the dimensional mean stream-wise velocity field in the midplane (Z=H). The wakes downstream of the cylinders are similar in LES and EB-RSM. The wakes are shorter downstream of the cylinders from row 4. Longer wakes are obtained for rows 1 and 2, then the lengths decrease. This is not observed at all for the &phi;-model and the k&omega;-SST model. One can state at this stage that eddy viscosity models are not able to predict the present flow. The values obtained with these two models are still plotted, but without analyzing deeply the results.


[[Image:U_Re3000_lineB_rows1_2_3_5.jpg]]<br/>
[[Image:U_Re3000_lineB_rows1_2_3_5.jpg]]<br/>
'''Figure ??:''' Mean stream-wse velocity component along line B - LES at <math>Re_D = 3000</math>
'''Figure 30:''' Mean stream-wise velocity component along line B - LES at <math>Re_D = 3000</math>


[[Image:U_Re10000_lineB_rows1_2_3_5.jpg]]<br/>
[[Image:U_Re10000_lineB_rows1_2_3_5.jpg]]<br/>
'''Figure ??:''' Mean stream-wse velocity component along line B - LES and URANS computations at <math>Re_D = 10000</math>
'''Figure 31:''' Mean stream-wise velocity component along line B - LES and URANS computations at <math>Re_D = 10000</math>


[[Image:U_Re30000_lineB_rows1_2_3_5.jpg]]<br/>
[[Image:U_Re30000_lineB_rows1_2_3_5.jpg]]<br/>
'''Figure ??:''' Mean stream-wse velocity component along line B - URANS computations at <math>Re_D = 30000</math>
'''Figure 32:''' Mean stream-wise velocity component along line B - URANS computations at <math>Re_D = 30000</math>


[[Image:U_Re30000_lineA1_rows2_3_5.jpg]]<br/>
[[Image:U_Re30000_lineA1_rows2_3_5.jpg]]<br/>
'''Figure ??:''' Mean stream-wse velocity component along line A1 - URANS computations at <math>Re_D = 30000</math>
'''Figure 33:''' Mean stream-wise velocity component along line A1 - URANS computations at <math>Re_D = 30000</math>
 
[[Image:Umean_z1_Re10000.jpg]]<br/>
'''Figure 34:''' Dimensional mean stream-wise velocity (m/s) in the midplane z=H at <math>Re_D = 10000</math>
 
[[Image:Umean_z1_Re30000.jpg]]<br/>
'''Figure 35:''' Dimensional mean stream-wise velocity (m/s)in the midplane z=H at <math>Re_D = 30000</math>


==== ''R.m.s.'' velocity profiles ====
==== ''R.m.s.'' velocity profiles ====
Figures 36 to 38 and figure 39 show the ''r.m.s.'' value of the stream-wise velocity along line B and line A1, respectively. The eddy viscosity models are very far from the experimental measurements (of course one cannot expect linear eddy viscosity models to yield realistic normal stresses, i.e. u', but only the sum k). On the other hand, LES and EB-RSM perform well compared to the experimental data with a certain superiority of LES. One can notice an overestimation of the ''r.m.s.'' value of the stream-wise velocity component with the EB-RSM model in the central part of row 5 (along both lines B and A1).
Figures 40 to 42 give the detailed profiles of the resolved and modelled stresses for the two Reynolds numbers at which the test case is performed with URANS approaches along lines B and A1, respectively. It is assumed here that the total Reynolds stress is given by the sum of the modelled part and the resolved one: <math>R^{tot}_{ij}=\overline{u'_iu'_j}+\overline{R}_{ij}</math>, where <math>u'_i</math> is the velocity fluctuation (resolved part) and <math>\overline{R}_{ij}</math> the time averaged value of the modelled Reynolds stress <math>{R}_{ij}</math>.
First, the EB-RSM model doesn't degenrate to a Large Eddy Simulation. The modelled part does't vanish and enhances the total level of turbulence. The linear eddy viscosity models fail naturally to predict the normal stresses. Figures 43 and 44 give the modelled and resolved stresses in the mid-plane in order to have a better overview of the behavior of this quantity (the whole computational domain is not represented to enhance clarity). The k&omega;-SST model is clearly steady upstream of row 5 for both Reynolds numbers. The &phi;-model shows some unsteadiness starting from row 3 at Re=10000 but is totally steady at Re=30000. Note that at this Reynolds number, the &phi;-model exhibited unsteadiness everywhere with a coarser mesh. This is an exemple of a computation which might exhibit unsteadiness at relatively coarse meshes and gives a steady solution after refinement (see Fadai-Ghotbi et al. (2008)).


[[Image:Urms_Re3000_lineB_rows1_2_3_5.jpg]]<br/>
[[Image:Urms_Re3000_lineB_rows1_2_3_5.jpg]]<br/>
'''Figure ??:''' ''R.m.s.'' of the stream-wise velocity component along line B - LES at <math>Re_D = 3000</math>
'''Figure 36:''' ''R.m.s.'' of the stream-wise velocity component along line B - LES at <math>Re_D = 3000</math>


[[Image:Urms_Re10000_lineB_rows1_2_3_5.jpg]]<br/>
[[Image:Urms_Re10000_lineB_rows1_2_3_5.jpg]]<br/>
'''Figure ??:''' ''R.m.s.'' of the stream-wise velocity component along line B - LES and URANS computations at <math>Re_D = 10000</math>
'''Figure 37:''' ''R.m.s.'' of the stream-wise velocity component along line B - LES and URANS computations at <math>Re_D = 10000</math>
 
[[Image:Urms_Re30000_lineB_rows1_2_3_5.jpg]]<br/>
'''Figure 38:''' ''R.m.s.'' of the stream-wise velocity component along line B - URANS computations at <math>Re_D = 30000</math>
 
[[Image:Urms_Re30000_lineA1_rows2_3_5.jpg]]<br/>
'''Figure 39:''' ''R.m.s.'' of the stream-wise velocity component along line A1 - URANS computations at <math>Re_D = 30000</math>


[[Image:Urms_Re10000_lineB_rows1_2_3_5_details.jpg]]<br/>
[[Image:Urms_Re10000_lineB_rows1_2_3_5_details.jpg]]<br/>
'''Figure ??:''' Resolved and modelled contributions to the ''r.m.s.'' of the stream-wise velocity component along line B - URANS computations at <math>Re_D = 10000</math>
'''Figure 40:''' Resolved and modelled contributions to the ''r.m.s.'' of the stream-wise velocity component along line B - URANS computations at <math>Re_D = 10000</math>


[[Image:Urms_Re30000_lineB_rows1_2_3_5.jpg]]<br/>
[[Image:Urms_Re30000_lineB_rows1_2_3_5_details.jpg]]<br/>
'''Figure ??:''' ''R.m.s.'' of the stream-wise velocity component along line B - URANS computations at <math>Re_D = 30000</math>
'''Figure 41:''' Resolved and modelled contributions to the ''r.m.s.'' of the stream-wise velocity component along line B - URANS computations at <math>Re_D = 30000</math>


[[Image:Urms_Re30000_lineB_rows1_2_3_5_details.jpg]]<br/>
[[Image:Urms_Re30000_lineA1_rows2_3_5_details.jpg]]<br/>
'''Figure ??:''' Resolved and modelled contributions to the ''r.m.s.'' of the stream-wise velocity component along line B - URANS computations at <math>Re_D = 30000</math>
'''Figure 42:''' Resolved and modelled contributions to the ''r.m.s.'' of the stream-wise velocity component along line A1 - URANS computations at <math>Re_D = 30000</math>


[[Image:Urms_Re30000_lineA1_rows2_3_5.jpg]]<br/>
[[Image:Mod_Res_Re10000.jpg]]<br/>
'''Figure ??:''' ''R.m.s.'' of the stream-wise velocity component along line A1 - URANS computations at <math>Re_D = 30000</math>
'''Figure 43:''' Resolved (left) and modelled (right) contributions to the dimensional ''r.m.s.'' of the stream-wise velocity component (m/s) in the midplane - LES and URANS computations at <math>Re_D = 10000</math>  


[[Image:Urms_Re30000_lineA1_rows2_3_5_details.jpg]]<br/>
[[Image:Mod_Res_Re30000.jpg]]<br/>
'''Figure ??:''' Resolved and modelled contributions to the ''r.m.s.'' of the stream-wise velocity component along line A1 - URANS computations at <math>Re_D = 30000</math>
'''Figure 44:''' Resolved (left) and modelled (right) contributions to the dimensional ''r.m.s.'' of the stream-wise velocity component (m/s) in the midplane - URANS computations at <math>Re_D = 30000</math>


==== Nusselt number on the bottom wall ====  
==== Nusselt number on the bottom wall ====  
Figures 45 to 47 show contours of the Nusselt number normalized by its averaged value over the bottom wall for the 8 rows.
LES computations exhibit an excellent behaviour compared to experimental values. The Nusselt number seems however overestimated at the locations of the horseshoe vortices and underestimated in the wake of the third row, in particular at the highest Reynolds number. Recall here that the error made on the average of the normalized Nusselt number is equal to 3% and 1% for <math>Re=3000</math> and <math>Re=10000</math>, respectively.
The &phi;-model and the k&omega;-SST model are inapt for predicting the physics and thus local quantities such as the Nusselt number.
The EB-RSM model exhibits a very satisfactory global behaviour in particular at the highest Reynolds number. The Nusselt number is however underestimated in the wake of the first cylinder as the unsteadiness is lower than predicted by LES. The prediction of the next rows is satisfactory although an overestimation of heat transfer in the wake of the cylinders for the last 4 rows is observed.
Figures 48 to 50 show detailed profiles of the normalized Nusselt number along particular lines (see "Test Case Studies" section). LES exhibits excellent results in particular after the second pin. The EB-RSM model shows some discrepencies for <math>Re=10000</math> and very good results at <math>Re=30000</math> in particular starting from row 5 (where the statistically periodic behavior is supposed to be reached).
[[Image:Image_Nusselt_Re3000.jpg]]<br/>
[[Image:Image_Nusselt_Re3000.jpg]]<br/>
'''Figure ??:''' Local Nusselt number on the bottom wall - LES computation at <math>Re_D = 3000</math> vs. experimental results from Ames et al.
'''Figure 45:''' Local Nusselt number normalized by its average (<math>Nus/Nus_{av}</math>) on the bottom wall - LES computation at <math>Re_D = 3000</math> vs. experimental results from Ames et al.


[[Image:Image_Nusselt_Re10000.jpg]]<br/>
[[Image:Image_Nusselt_Re10000.jpg]]<br/>
'''Figure ??:''' Local Nusselt number on the bottom wall - LES and URANS computations at <math>Re_D = 10000</math> vs. experimental results from Ames et al.
'''Figure 46:''' Local Nusselt number normalized by its average (<math>Nus/Nus_{av}</math>) on the bottom wall - LES and URANS computations at <math>Re_D = 10000</math> vs. experimental results from Ames et al.


[[Image:Image_Nusselt_Re30000.jpg]]<br/>
[[Image:Image_Nusselt_Re30000.jpg]]<br/>
'''Figure ??:''' Local Nusselt number on the bottom wall - URANS computations at <math>Re_D = 30000</math> vs. experimental results from Ames et al.
'''Figure 47:''' Local Nusselt number normalized by its average (<math>Nus/Nus_{av}</math>) on the bottom wall - URANS computations at <math>Re_D = 30000</math> vs. experimental results from Ames et al.


[[Image:Nus_o_NusAv_different_y_Re3000_LES_y0_y1.25_improved.jpg]]<br/>
[[Image:Nus_o_NusAv_different_y_Re3000_LES_y0_y1.25_improved.jpg]]<br/>
'''Figure ??:''' Local Nusselt number on the bottom wall - <math>Re_D = 3000</math> - top: Y=0, bottom: Y=1.25 - symb.: Exp. - black: LES
'''Figure 48:''' Local Nusselt number normalized by its average (<math>Nus/Nus_{av}</math>) - <math>Re_D = 3000</math> - top: Y=0, bottom: Y=1.25 - symb.: Exp. - black: LES


[[Image:Nus_o_NusAv_different_y_Re10000_all_y0_y1.25_improved.jpg]]<br/>
[[Image:Nus_o_NusAv_different_y_Re10000_all_y0_y1.25_improved.jpg]]<br/>
'''Figure ??:''' Local Nusselt number on the bottom wall - <math>Re_D = 10000</math> - top: Y=0, bottom: Y=1.25 - symb.: Exp. - black: LES - red: EBRSM - blue: Phi-model - green: kw-SST
'''Figure 49:''' Local Nusselt number normalized by its average (<math>Nus/Nus_{av}</math>) on the bottom wall - <math>Re_D = 10000</math> - top: Y=0, bottom: Y=1.25 - symb.: Exp. - black: LES - red: EBRSM - blue: Phi-model - green: kw-SST


[[Image:Nus_o_NusAv_different_y_Re30000_y0_y1.25_improved.jpg]]<br/>
[[Image:Nus_o_NusAv_different_y_Re30000_y0_y1.25_improved.jpg]]<br/>
'''Figure ??:''' Local Nusselt number on the bottom wall - <math>Re_D = 30000</math> - top: Y=0, bottom: Y=1.25 - symb.: Exp. - red: EBRSM - blue: Phi-model - green: kw-SST  
'''Figure 50:''' Local Nusselt number normalized by its average (<math>Nus/Nus_{av}</math>) on the bottom wall - <math>Re_D = 30000</math> - top: Y=0, bottom: Y=1.25 - symb.: Exp. - red: EBRSM - blue: Phi-model - green: kw-SST  





Latest revision as of 14:43, 12 February 2017

Flow and heat transfer in a pin-fin array

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References

Confined Flows

Underlying Flow Regime 4-18

Evaluation

Comparison of CFD Calculations with Experiments


Global comparisons

Pressure loss coefficient

The following table gives the numerical values of the presure loss coefficient in Ames et al. experiment as well as those obtained in the main computations. Note that the experimental data shown here have been sent by the authors of the experiment and are different from the ones found in Ames et al. (2005).

Figure 25 shows these values and compares them to existing correlations from the literature (note that Ames et al. experimental values are very close to Metzger et al. correlation):

  • Metzger et al. (1982) :

for , Failed to parse (syntax error): {\displaystyle f = 0.317 {Re_D}^{−0.132} } , for , Failed to parse (syntax error): {\displaystyle f = 1.76 {Re_D}^{−0.318}}

  • Jacob (1938) :

Table 5 compares various comptations to the experiment and gives the relative errors. Large Eddy Simulation shows a very good agreement with the correlations and Ames et al. experimental results. The relative error is equal to 1% end 3% for the two computed Reynolds numbers, respectively. The EB-RSM model shows a very satisfactory pressure loss coefficient at (1% error) but the result is worse at (10% error). The kω-SST model shows globally good results (7% and 6% errors) and the φ-model poorer results with an error of 17% and 12% for the two Reynolds numbers, respectively.

Wall resolved LES is therefore a good candidate to obtain satisfactory pressure loss coefficients in the present configuration. The EB-RSM model gives very good results at the lowest Reynolds number but a substantial error is observed at . Although the kω-SST model exhibits pretty good results, the physics of the present flow is not well captured as it will be shown later and one can't rely on this model. The φ-model gives the worst results and is not reliable.

Table 5: Pressure loss coefficient values for the different computations compared to Ames et al. experimental values.
Ames et al. LES3000 Ames et al. LES10000 EBRSM10000 PHI10000 kwSST10000 Ames et al. EBRSM30000 PHI30000 kwSST30000
Value 0.111 0.1076 0.095 0.0939 0.0963 0.0789 0.0883 0.07 0.0771 0.0615 0.0744
Error (%) - 3 - 1 1 17 7 - 10 12 6


Friction Final Results.jpg
Figure 25: Pressure loss coefficient

Average Nusselt Number

Table 6 gives the numerical values of the average Nusselt numbers on the bottom wall in Ames et al. experiment as well as those obtained in the main computations (go to "Test Case Studies" section to find the definition of the Nusselt number and its average on the bottom wall).

Figure 26 shows these values and compares them to existing correlations from the literature (Note again that the Ames et al. experimental results are close to the Metzger correlation):

  • Metzger et al. (1982) :

  • Van Fossen (1982) :

Note that the local experimental results have an error of +/- 12%, +/-11.4%, and +/-10.5% for the 3 000, 10 000, and 30 000 Reynolds numbers, respectively, in the endwall regions adjacent to the pins and +/- 9% away from the pin.

Table 6 compares the computations with the experiment and gives the relative errors. One can see that the two models using an eddy viscosity (φ-model and kω-SST) are too far from the experimental results, except the φ-model at . However, as it will be shown later, this model exhibits wrong physics. LES and EB-RSM yield very satisfactory values in the error interval. The superiority of these two models will be discussed later when dealing with the local normalized Nusselt number on the bottom wall. Ames and Dvorak (2006a) used standard, Realizable k-ε and RNG k-ε models. They showed a substantial underestimation of the average Nusselt number when the pins were heated. Delibra et al. (2009) used the URANS ζ-f model (very close to the φ-model utilized in the present work) at the 10000 and 30000 Reynolds numbers. They obtained average Nusselt numbers of 46.2 and 122.3, respectively. These values are further from the experimental ones. The temperature was imposed at the wall in Delibra et al. (2009). Imposing the temperature at the wall in the present work with the φ-model led to average Nusselt numbers of 41 and 100, respectively. This shows that imposing the temperature gives higher Nusselt numbers than the ones obtained by imposing the heat flux with this model. The effect of the temperature boundary condition has also been tested with the LES at and the EB-RSM model at . The average Nusselt numbers were equal to 46 and 115.2, respectively. No clear conclusion can be drawn for these two models which give a more realistic physics than the ones using an eddy viscosity concerning the effect of imposing the temperature instead of the heat flux at the bottom wall. However, one can state that there is a difference between the two methods. Delibra et al. (2008) also reported LES computations. The one carried out for was too coarse with a short time integration and will not be considered here. LES at gave an average Nusselt number equal to 44.3 which is lower than the one obtained in the present work even with a fixed temperature at the wall. The error compared to the experimental value is equal to 18% which is too high. This discrepancy can be attributed to two factors:

  • the mesh was probably still too coarse, not in the wall normal directions but in the stream-wise and span-wise ones (the mesh in the present study contains 76 million cells and their mesh contained 5.5 million cells),
  • the boundary condition using an imposed temperature is less representative of the experiment.

Finally, Delibra et al. (2010) switched to use a hybrid RANS/LES approach to overcome the understimation of the Nusselt number but the results where not much better than the ones obtained with the URANS approach based on the ζ-f model.

Table 6: Average Nusselt number on the bottom wall
Ames3000 LES3000 Ames10000 LES10000 DLES10000(*) EBRSM10000 PHI10000 kwSST10000 DURANS10000(*) Ames30000 EBRSM30000 PHI30000 kwSST30000 DURANS30000(*)
Value 22.1 21.6 54.1 48.6 44.1 48.1 37.6 33.1 46.2 111.5 114.8 93.3 72 122.3
Error (%) - 2 - 10 18 11 30 39 14 - 3 16 35 9

(*) Delibra et al. results.

NusAv Final Results.jpg
Figure 26: Average Nusselt number on the bottom wall

Local comparisons

Mean pressure coefficient

Figures 27 to 29 give the profiles of the pressure coefficient along the midline of pins 1, 2, 3 and 5 for the three Reynolds numbers.

Note that the experimental results have an uncertainty of +/-0.075 at and of +/-0.025 at and .

LES results are in relatively good agreement with the experimental data. EB-RSM exhibits decent results as well. The two eddy viscosity models show poor results in particular the φ-model. The prediction of the coefficient may be related to the wake prediction downstream of the pins.


Cp midplane Re3000 rows1 2 3 5.jpg
Figure 27: Pressure coefficient along the midline of the pins - LES at

Cp midplane Re10000 rows1 2 3 5.jpg
Figure 28: Pressure coefficient along the midline of the pins - LES and URANS computations at

Cp midplane Re30000 rows1 2 3 5.jpg
Figure 29: Pressure coefficient along the midline of the pins - URANS computations at

Mean velocity profiles and contours

Figures 30 to 32 show the mean axial velocity component along line B for the different computations. The two eddy viscosity models exhibit a dramatic deficit of the mean axial velocity. On the other hand, LES and EB-RSM perform well compared to the experimental data with a slight overestimation of the mean stream-wise velocity peaks close to the walls in several locations, in particular departing from row 3. Note that the convergence study shown in the "Test Case Studies" section proved that the needed refinement to obtain a convergent solution around pin 2 is very important (with the first two levels of refinement, a deficit in the central velocity was observed).

Figure 33 shows the means axial velocity component along line A1 for the highest Reynolds number. One can see again that the linear eddy viscosity models are unable to predict the present flow. The results obtained with the EB-RSM model are much closer to the experimental data. However, the flow redistribution doesn't seem to be similar to the one obtained experimentally. This may come from the periodic boundary conditions used in the present work v.s. confined conditions in the experiment.

Figures 34 and 35 display the dimensional mean stream-wise velocity field in the midplane (Z=H). The wakes downstream of the cylinders are similar in LES and EB-RSM. The wakes are shorter downstream of the cylinders from row 4. Longer wakes are obtained for rows 1 and 2, then the lengths decrease. This is not observed at all for the φ-model and the kω-SST model. One can state at this stage that eddy viscosity models are not able to predict the present flow. The values obtained with these two models are still plotted, but without analyzing deeply the results.

U Re3000 lineB rows1 2 3 5.jpg
Figure 30: Mean stream-wise velocity component along line B - LES at

U Re10000 lineB rows1 2 3 5.jpg
Figure 31: Mean stream-wise velocity component along line B - LES and URANS computations at

U Re30000 lineB rows1 2 3 5.jpg
Figure 32: Mean stream-wise velocity component along line B - URANS computations at

U Re30000 lineA1 rows2 3 5.jpg
Figure 33: Mean stream-wise velocity component along line A1 - URANS computations at

Umean z1 Re10000.jpg
Figure 34: Dimensional mean stream-wise velocity (m/s) in the midplane z=H at

Umean z1 Re30000.jpg
Figure 35: Dimensional mean stream-wise velocity (m/s)in the midplane z=H at

R.m.s. velocity profiles

Figures 36 to 38 and figure 39 show the r.m.s. value of the stream-wise velocity along line B and line A1, respectively. The eddy viscosity models are very far from the experimental measurements (of course one cannot expect linear eddy viscosity models to yield realistic normal stresses, i.e. u', but only the sum k). On the other hand, LES and EB-RSM perform well compared to the experimental data with a certain superiority of LES. One can notice an overestimation of the r.m.s. value of the stream-wise velocity component with the EB-RSM model in the central part of row 5 (along both lines B and A1).

Figures 40 to 42 give the detailed profiles of the resolved and modelled stresses for the two Reynolds numbers at which the test case is performed with URANS approaches along lines B and A1, respectively. It is assumed here that the total Reynolds stress is given by the sum of the modelled part and the resolved one: , where is the velocity fluctuation (resolved part) and the time averaged value of the modelled Reynolds stress . First, the EB-RSM model doesn't degenrate to a Large Eddy Simulation. The modelled part does't vanish and enhances the total level of turbulence. The linear eddy viscosity models fail naturally to predict the normal stresses. Figures 43 and 44 give the modelled and resolved stresses in the mid-plane in order to have a better overview of the behavior of this quantity (the whole computational domain is not represented to enhance clarity). The kω-SST model is clearly steady upstream of row 5 for both Reynolds numbers. The φ-model shows some unsteadiness starting from row 3 at Re=10000 but is totally steady at Re=30000. Note that at this Reynolds number, the φ-model exhibited unsteadiness everywhere with a coarser mesh. This is an exemple of a computation which might exhibit unsteadiness at relatively coarse meshes and gives a steady solution after refinement (see Fadai-Ghotbi et al. (2008)).


Urms Re3000 lineB rows1 2 3 5.jpg
Figure 36: R.m.s. of the stream-wise velocity component along line B - LES at

Urms Re10000 lineB rows1 2 3 5.jpg
Figure 37: R.m.s. of the stream-wise velocity component along line B - LES and URANS computations at

Urms Re30000 lineB rows1 2 3 5.jpg
Figure 38: R.m.s. of the stream-wise velocity component along line B - URANS computations at

Urms Re30000 lineA1 rows2 3 5.jpg
Figure 39: R.m.s. of the stream-wise velocity component along line A1 - URANS computations at

Urms Re10000 lineB rows1 2 3 5 details.jpg
Figure 40: Resolved and modelled contributions to the r.m.s. of the stream-wise velocity component along line B - URANS computations at

Urms Re30000 lineB rows1 2 3 5 details.jpg
Figure 41: Resolved and modelled contributions to the r.m.s. of the stream-wise velocity component along line B - URANS computations at

Urms Re30000 lineA1 rows2 3 5 details.jpg
Figure 42: Resolved and modelled contributions to the r.m.s. of the stream-wise velocity component along line A1 - URANS computations at

Mod Res Re10000.jpg
Figure 43: Resolved (left) and modelled (right) contributions to the dimensional r.m.s. of the stream-wise velocity component (m/s) in the midplane - LES and URANS computations at

Mod Res Re30000.jpg
Figure 44: Resolved (left) and modelled (right) contributions to the dimensional r.m.s. of the stream-wise velocity component (m/s) in the midplane - URANS computations at

Nusselt number on the bottom wall

Figures 45 to 47 show contours of the Nusselt number normalized by its averaged value over the bottom wall for the 8 rows.

LES computations exhibit an excellent behaviour compared to experimental values. The Nusselt number seems however overestimated at the locations of the horseshoe vortices and underestimated in the wake of the third row, in particular at the highest Reynolds number. Recall here that the error made on the average of the normalized Nusselt number is equal to 3% and 1% for and , respectively.

The φ-model and the kω-SST model are inapt for predicting the physics and thus local quantities such as the Nusselt number.

The EB-RSM model exhibits a very satisfactory global behaviour in particular at the highest Reynolds number. The Nusselt number is however underestimated in the wake of the first cylinder as the unsteadiness is lower than predicted by LES. The prediction of the next rows is satisfactory although an overestimation of heat transfer in the wake of the cylinders for the last 4 rows is observed.

Figures 48 to 50 show detailed profiles of the normalized Nusselt number along particular lines (see "Test Case Studies" section). LES exhibits excellent results in particular after the second pin. The EB-RSM model shows some discrepencies for and very good results at in particular starting from row 5 (where the statistically periodic behavior is supposed to be reached).

Image Nusselt Re3000.jpg
Figure 45: Local Nusselt number normalized by its average () on the bottom wall - LES computation at vs. experimental results from Ames et al.

Image Nusselt Re10000.jpg
Figure 46: Local Nusselt number normalized by its average () on the bottom wall - LES and URANS computations at vs. experimental results from Ames et al.

Image Nusselt Re30000.jpg
Figure 47: Local Nusselt number normalized by its average () on the bottom wall - URANS computations at vs. experimental results from Ames et al.

Nus o NusAv different y Re3000 LES y0 y1.25 improved.jpg
Figure 48: Local Nusselt number normalized by its average () - - top: Y=0, bottom: Y=1.25 - symb.: Exp. - black: LES

Nus o NusAv different y Re10000 all y0 y1.25 improved.jpg
Figure 49: Local Nusselt number normalized by its average () on the bottom wall - - top: Y=0, bottom: Y=1.25 - symb.: Exp. - black: LES - red: EBRSM - blue: Phi-model - green: kw-SST

Nus o NusAv different y Re30000 y0 y1.25 improved.jpg
Figure 50: Local Nusselt number normalized by its average () on the bottom wall - - top: Y=0, bottom: Y=1.25 - symb.: Exp. - red: EBRSM - blue: Phi-model - green: kw-SST




Contributed by: Sofiane Benhamadouche — EDF

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References


© copyright ERCOFTAC 2024