UFR 4-06 Evaluation: Difference between revisions

From KBwiki
Jump to navigation Jump to search
m (Dave.Ellacott moved page Silver:UFR 4-06 Evaluation to UFR 4-06 Evaluation over redirect)
 
(13 intermediate revisions by 4 users not shown)
Line 1: Line 1:
{{UFR|front=UFR 4-06|description=UFR 4-06 Description|references=UFR 4-06 References|testcase=UFR 4-06 Test Case|evaluation=UFR 4-06 Evaluation|qualityreview=UFR 4-06 Quality Review|bestpractice=UFR 4-06 Best Practice Advice|relatedACs=UFR 4-06 Related ACs}}
{{UFR|front=UFR 4-06|description=UFR 4-06 Description|references=UFR 4-06 References|testcase=UFR 4-06 Test Case|evaluation=UFR 4-06 Evaluation|qualityreview=UFR 4-06 Quality Review|bestpractice=UFR 4-06 Best Practice Advice|relatedACs=UFR 4-06 Related ACs}}


Line 17: Line 16:
== Comparison of CFD calculations with Experiments ==
== Comparison of CFD calculations with Experiments ==


The experimental results of Clausen et al [[U4-06ref.htm#ref10 10]] are included on the [http://tmdb.ws.tn.tudelft.nl/database/test60/test60.html 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]).
The experimental results of Clausen et al [[UFR 4-06 References#ref10|[10]]] are included on the [http://tmdb.ws.tn.tudelft.nl/database/test60/test60.html ERCOFTAC database] for each measurement station identified in [[UFR 4-06 Test Case#Test Case Experiments|section 4]]. Armfield et al [[UFR 4-06 References#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 [[#figure5|5]], [[#figure6|6]], [[#figure7|7]], [[#figure8|8]] and [[#figure9|9]]. Measurement stations 3, 5 and 7 are located 100mm, 250mm and 405mm downstream of the diffuser inlet respectively (see [[UFR 4-06 Test Case#fig4|Figure 4]]).


<span id="figure5"></span>
<center>[[Image:U4-06d32_files_image017.jpg]]</center>
<center>[[Image:U4-06d32_files_image017.jpg]]</center>


<center>Figure 5 Radial Variation of Normalised Axial Velocity</center>
<center>Figure 5 Radial Variation of Normalised Axial Velocity</center>


<center>(<sub>[[Image:U4-06d32_files_image019.gif]]</sub>) ASM, (----) k-ε [[U4-06ref.htm#ref4 ]</center>
<center>([[Image:U4-06d32_files_image019.gif]]) ASM, (----) k-&epsilon; [[UFR 4-06 References#ref4|[4]]]</center>


<span id="figure6"></span>
<center>[[Image:U4-06d32_files_image021.jpg]]</center>
<center>[[Image:U4-06d32_files_image021.jpg]]</center>


<center>Figure 6 Radial Variation of Normalised Swirl Velocity</center>
<center>Figure 6 Radial Variation of Normalised Swirl Velocity</center>


<center>(<sub>[[Image:U4-06d32_files_image022.gif]]</sub>) ASM, (----) k-ε [[U4-06ref.htm#ref4 ]</center>
<center>([[Image:U4-06d32_files_image022.gif]]) ASM, (----) k-&epsilon; [[UFR 4-06 References#ref4|[4]]]</center>


<span id="figure7"></span>
<center>[[Image:U4-06d32_files_image024.jpg]]</center>
<center>[[Image:U4-06d32_files_image024.jpg]]</center>


<center>Figure 7 Radial Variation of Normalised Turbulent Kinetic Energy</center>
<center>Figure 7 Radial Variation of Normalised Turbulent Kinetic Energy</center>


<center>(<sub>[[Image:U4-06d32_files_image025.gif]]</sub>) ASM, (----) k-ε [[U4-06ref.htm#ref4 ]</center>
<center>([[Image:U4-06d32_files_image025.gif]]) ASM, (----) k-&epsilon; [[UFR 4-06 References#ref4|[4]]]</center>


<span id="figure8"></span>
<center>[[Image:U4-06d32_files_image027.jpg]]</center>
<center>[[Image:U4-06d32_files_image027.jpg]]</center>


<center>Figure 8 Radial Variation of Normalised Reynolds Stress <sub>[[Image:U4-06d32_files_image029.gif]] </sub></center>
<center>Figure 8 Radial Variation of Normalised Reynolds Stress [[Image:U4-06d32_files_image029.gif]] </center>


<center>(<sub>[[Image:U4-06d32_files_image030.gif]]</sub>) ASM, (----) k-ε [[U4-06ref.htm#ref4 ]</center>
<center>([[Image:U4-06d32_files_image030.gif]]) ASM, (----) k-&epsilon; [[UFR 4-06 References#ref4|[4]]]</center>


<span id="figure9"></span>
<center>[[Image:U4-06d32_files_image032.jpg]]</center>
<center>[[Image:U4-06d32_files_image032.jpg]]</center>


<center>Figure 9 Radial Variation of Normalised Reynolds Stress <sub>[[Image:U4-06d32_files_image034.gif]] </sub></center>
<center>Figure 9 Radial Variation of Normalised Reynolds Stress [[Image:U4-06d32_files_image034.gif]] </center>


<center>(<sub>[[Image:U4-06d32_files_image035.gif]]</sub>) ASM, (----) k-ε (the dotted line is the ASM without wall echo terms) [[U4-06ref.htm#ref4 ]</center>
<center>([[Image:U4-06d32_files_image035.gif]]) ASM, (----) k-&epsilon; (the dotted line is the ASM without wall echo terms) [[UFR 4-06 References#ref4|[4]]]</center>


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]].
Based on the findings of Armfield et al [[UFR 4-06 References#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 [[UFR 4-06 References#ref17|[17]]], Cho and Fletcher [[UFR 4-06 References#ref12|[12]]] and Page et al [[UFR 4-06 References#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.
The summary report of EDF [[UFR 4-06 References#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 <sub>[[Image:U4-06d32_files_image037.gif]]levels also in the initial diffuser region. The low Re k-</sub>ε model calculations also develop excessive levels of k and of the shear stresses [[Image:U4-06d32_files_image038.gif]] [[Image:U4-06d32_files_image040.gif]] [[Image:U4-06d32_files_image041.gif]] 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 [[Image:U4-06d32_files_image043.gif]]stress. For the other turbulence quantities the standard k-ε model gives results equally close to the measurements.''
''&bull; 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-&omega; model calculation was far too high, causing excessive k and [[Image:U4-06d32_files_image037.gif]] levels also in the initial diffuser region. The low Re k-&epsilon; model calculations also develop excessive levels of k and of the shear stresses [[Image:U4-06d32_files_image038.gif]] and [[Image:U4-06d32_files_image040.gif]], which persist throughout the whole diffuser. In particular, the large value of [[Image:U4-06d32_files_image041.gif]] 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 [[UFR 4-06 References#ref12|[12]]]] stay high all the way through the diffuser while the RSM calculations of EC Lyon [see reference [[UFR 4-06 References#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 [[UFR 4-06 References#ref12|[12]]]] with the standard k-&epsilon; model and two non-linear k-&epsilon; 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 [[Image:U4-06d32_files_image043.gif]] stress. For the other turbulence quantities the standard k-&epsilon; 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.''
''&bull; 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-&epsilon; 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 [[UFR 4-06 References#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-&epsilon; 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-&epsilon; model calculations also produced reasonable distributions of the axial and circumferential 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.”''
''&bull; Altogether it seems that the standard k-&epsilon; 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.''


<font size="-2" color="#888888">© copyright ERCOFTAC 2004</font><br />
<font size="-2" color="#888888">© copyright ERCOFTAC 2004</font><br />
Line 67: Line 71:


{{UFR|front=UFR 4-06|description=UFR 4-06 Description|references=UFR 4-06 References|testcase=UFR 4-06 Test Case|evaluation=UFR 4-06 Evaluation|qualityreview=UFR 4-06 Quality Review|bestpractice=UFR 4-06 Best Practice Advice|relatedACs=UFR 4-06 Related ACs}}
{{UFR|front=UFR 4-06|description=UFR 4-06 Description|references=UFR 4-06 References|testcase=UFR 4-06 Test Case|evaluation=UFR 4-06 Evaluation|qualityreview=UFR 4-06 Quality Review|bestpractice=UFR 4-06 Best Practice Advice|relatedACs=UFR 4-06 Related ACs}}
[[Category:Underlying Flow Regime]]

Latest revision as of 14:14, 12 February 2017

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References




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 [10] are included on the ERCOFTAC database for each measurement station identified in section 4. Armfield et al [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 5, 6, 7, 8 and 9. Measurement stations 3, 5 and 7 are located 100mm, 250mm and 405mm downstream of the diffuser inlet respectively (see Figure 4).

U4-06d32 files image017.jpg
Figure 5 Radial Variation of Normalised Axial Velocity
(U4-06d32 files image019.gif) ASM, (----) k-ε [4]

U4-06d32 files image021.jpg
Figure 6 Radial Variation of Normalised Swirl Velocity
(U4-06d32 files image022.gif) ASM, (----) k-ε [4]

U4-06d32 files image024.jpg
Figure 7 Radial Variation of Normalised Turbulent Kinetic Energy
(U4-06d32 files image025.gif) ASM, (----) k-ε [4]

U4-06d32 files image027.jpg
Figure 8 Radial Variation of Normalised Reynolds Stress U4-06d32 files image029.gif
(U4-06d32 files image030.gif) ASM, (----) k-ε [4]

U4-06d32 files image032.jpg
Figure 9 Radial Variation of Normalised Reynolds Stress U4-06d32 files image034.gif
(U4-06d32 files image035.gif) ASM, (----) k-ε (the dotted line is the ASM without wall echo terms) [4]

Based on the findings of Armfield et al [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 [17], Cho and Fletcher [12] and Page et al [9].

The summary report of EDF [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 U4-06d32 files image037.gif levels also in the initial diffuser region. The low Re k-ε model calculations also develop excessive levels of k and of the shear stresses U4-06d32 files image038.gif and U4-06d32 files image040.gif, which persist throughout the whole diffuser. In particular, the large value of U4-06d32 files image041.gif 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 [12]] stay high all the way through the diffuser while the RSM calculations of EC Lyon [see reference [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 [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 U4-06d32 files image043.gif 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 [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 circumferential 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


Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References