UFR 4-06 Evaluation
Swirling diffuser flow
Underlying Flow Regime 4-06 © copyright ERCOFTAC 2004
Evaluation
Comparison of CFD calculations with Experiments
The experimental results of Clausen et al U4-06ref.htm#ref10 10 are included on the ERCOFTAC database for each measurement station identified in [U4-06exp.htm#section4 section 4]. Armfield et al U4-06ref.htm#ref4 4 presented a comparison of the CFD results with the experimental data including radial variation of axial velocity, radial variation of swirl velocity, radial variation of turbulent kinetic energy and radial variation of Reynolds stresses. Examples of these comparisons are shown in Figures [U4-06eva.htm#figure5 5], [U4-06eva.htm#figure6 6], [U4-06eva.htm#figure7 7], [U4-06eva.htm#figure8 8] and [U4-06eva.htm#figure9 9]. Measurement stations 3, 5 and 7 are located 100mm, 250mm and 405mm downstream of the diffuser inlet respectively (see [U4-06exp.htm#fig4 Figure 4]).
) ASM, (----) k-ε [[U4-06ref.htm#ref4 ]
) ASM, (----) k-ε [[U4-06ref.htm#ref4 ]
) ASM, (----) k-ε [[U4-06ref.htm#ref4 ]
) ASM, (----) k-ε [[U4-06ref.htm#ref4 ]
) ASM, (----) k-ε (the dotted line is the ASM without wall echo terms) [[U4-06ref.htm#ref4 ]Based on the findings of Armfield et al U4-06ref.htm#ref4 4, overall, good prediction of the mean and turbulence quantities was obtained using both turbulence models with a two layer wall function based approach. However, the ASM gave slightly better overall predictions of the turbulence quantities. The use of a two layer, rather than a single layer, wall function was found to be necessary to predict accurately the level, location and the axial variation of the near wall peak in turbulence quantities. These findings have been further supported by the work of Armfield and Fletcher U4-06ref.htm#ref17 17, Cho and Fletcher U4-06ref.htm#ref12 12 and Page et al U4-06ref.htm#ref9 9.
The summary report of EDF U4-06ref.htm#ref11 11 discussed the results of eleven simulations of the test case flow, as described in section 5, and the findings are reproduced below.
In some of the calculations, the measured k-distribution at the inflow boundary was not used as the initial profile in the boundary layer, but different (mainly lower) distributions were used. The initial k-level in the k-ω model calculation was far too high, causing excessive k and 
levels also in the initial diffuser region. The low Re k-ε model calculations also develop excessive levels of k and of the shear stresses 
and 
, which persist throughout the whole diffuser. In particular, the large value of 
near the wall is strange and entirely unrealistic. Also, the swirl velocity in these calculations is somewhat higher than in the experiments at most stations, indicating that the swirl number was not exactly matched to the experiments. The ASM and Reynolds stress model calculations of k and shear stresses are also initially all too high. The ASM and the RSM calculations from TH Aachen [see reference [U4-06ref.htm#ref12 12]] stay high all the way through the diffuser while the RSM calculations of EC Lyon [see reference [U4-06ref.htm#ref12 12]] fall below the experimental levels towards the outlet of the diffuser. There are hence considerable differences in the calculated turbulence quantities between the calculations using nominally the same LRR Reynolds stress model. The best argument concerning the turbulence quantities was obtained by U. Brussels [see reference [U4-06ref.htm#ref12 12]] with the standard k-ε model and two non-linear k-ε model versions. The Hirsch-Khodak version produces the better results and is in fact the only model that can simulate reasonably well the distribution of the 
stress. For the other turbulence quantities the standard k-ε model gives results equally close to the measurements.
The reduction in U velocity in the center portion towards the outlet of the diffuser is under-predicted by most models. Only the low Re k-ε model calculation produces nearly the correct reduction. This may have to do with the fact that, compared with the other calculations, this model yields higher swirl velocity and shear stresses which are, however, also higher than the experimentally observed values. The LRR calculation from Lyon [see reference [U4-06ref.htm#ref12 12]] also yields lower U-velocity in the center portion in better agreement with the measurements, but strangely this is associated with smaller shear stresses (significantly below the experimental levels) than in the other models. With these two exceptions, there is generally good agreement among the various predictions about the development of the U and W velocity profiles; only at the station closest to the outlet there is considerable spread, the k-ε model yielding too high U velocity at the axis. The streamwise distribution of pressure in the diffuser is calculated well by all models. The various k-ε model calculations also produced reasonable distributions of the axial and circumfrential wall shear stress, while some of the other models yield too high values.
Altogether it seems that the standard k-ε model with wall functions gives results in this case which, apart from an insufficient reduction in the center line velocity, are reasonable for all quantities that are of practical interest. More complex models do not seem to yield a marked improvement. It appears that the swirl is not strong enough in this case to show the superiority of Reynolds stress models reported in the literature for swirling flows.
© copyright ERCOFTAC 2004
Contributors: Chris Carey - Fluent Europe Ltd