UFR 3-03 Best Practice Advice: Difference between revisions

From KBwiki
Jump to navigation Jump to search
(New page: {{UFR|front=UFR 3-03|description=UFR 3-03 Description|references=UFR 3-03 References|testcase=UFR 3-03 Test Case|evaluation=UFR 3-03 Evaluation|qualityreview=UFR 3-03 Quality Review|bestp...)
 
No edit summary
Line 13: Line 13:




= Evaluation =
= Best Practice Advice =


== Comparison of CFD calculations with Experiments ==
== Best Practice Advice for the UFR ==


Several model comparison studies for this testcase have been published by Menter (1993, 1994). Additional results for this case can be found in Bardina et al. (1997) and a comprehensive validation of EARSM models in comparison with conventional two-equation is given by Hellsten and Laine (2000) and Hsu et al. (2003). As not all relevant solutions are available as data sets, the comparison is based on pictures scanned from existing articles and one set of new simulation using CFX-5. The following comparisons will be presented:
This CS0 diffuser testcase is relatively simple to compute concerning the set-up of the flow and the boundary conditions. A grid of 61 x 121 nodes is recommended with equal spacing in the streamwise direction and a stretching in wall normal direction, which is consistent with the wall formulation of the CFD method. The use of a low-Reynolds number grid is recommended if allowed by the code and the wall treatment. The set-up in terms of the boundary conditions is described in chapter 6. The inlet profiles have to be taken from the experimental data, both for the velocity and the turbulence. Bulk profiles are not acceptable, as they alter the boundary layer development and the displacement effects. The outer boundary is to be treated as a slip wall using the contour obtained from the integration of the experimental velocity profiles.


'''<u>Set 1:</u>'''
The main issue for the accurate prediction of this case is the selection of the turbulence model. The best solutions, known to the author are from the SST model, the Spalart-Allmaras model and from ω-equation based Reynolds Stress models. Comparable results to the SST model have been obtained with <span lang="EN-US"><font face="Symbol">w-</font></span>equation based Explicit Algebraic Reynolds Stress Models (Hellsten and Laine 2000).


•        SST model
Turbulence modeling is an ever-developing area and it is likely that more and more models will prove in the future to be able to handle this testcase. This will contribute significantly to the reliability of CFD methods for flows with pressure-induced separation.


•        BSL model.
The CS0 case does not allow judging the ability of turbulence models to properly predict re-attachment of the flow and flow-recovery downstream of a separation zone. It is known from other studies that all RANS turbulence models seem to have significant problems in that area. Particularly models, which predict the correct separation onset, are prone to give larger separation zones and slower flow recovery than indicated by the data. Note that this does not mean that the other models are superior in that respect. It is only due to a cancellation of errors that models, which fail to predict the separation onset, have little problems in the downstream region. The user is advised to keep this deficiency in mind when judging CFD simulations of diffuser flows.


•        k-ω model
== Recommendations for Future Work ==


•        k-<span lang="DE"><font face="Symbol">e </font></span>model
Develop model enhancements or new models, which can predict both, flow separation and flow recovery downstream of reattachment. Note that a highly accurate separation prediction capability is essential to make the second step. The improvement cannot be judged with the CS0 case, but the Buice and Eaton (1997) experiments can serve as a guideline.
 
'''<u>Set 2:</u>'''
 
•        SST model
 
•        RNG k-ε model.
 
•        SA model
 
•        k-<span lang="DE"><font face="Symbol">e </font></span>model
 
'''<u>Set 3:</u>'''
 
•        k-<span lang="DE"><font face="Symbol">e </font></span>model
 
•        SSG model (ε-equation)
 
'''<u>Set 4:</u>'''
 
•        Comparison from Hsu et al. (2003)
 
'''<u>Set 5:</u>'''
 
•        Comparison by Hellsten and Laine (2000)
 
All simulations have been performed on low-RE grids using an integration to the surface.
 
The main concern with this flow is the correct prediction of the flow reversal near the wall and the proper displacement effect (reflected in the wall pressure distribution).
 
== 1.1        Results from the Present Author ==
 
Figure 4 shows a comparison of velocity profiles from the k-ε model (Jones and Launder 1972), the k-ω (Wilcox, 1993) model, the BSL and SST model (Menter, 1993, 1994) (Set 1). The simulations have been performed with the INS2D code and have been taken from Menter (1993). The velocity profiles show clearly that the k-ε is not suitable for flows with adverse pressure gradients and pressure induced separation. The model shows little sensitivity to the flow deceleration and predicts firmly attached velocity profiles. The k-ω model and the BSL model are somewhat closer to the data, but still in poor agreement with the experimental profiles. The improvement relative to the k-εmodelis a result of the improved near wall behavior of the ω-equation based models. It was shown by Wilcox (1993) that the ω-equation is superior to the ε-equation in its near wall prediction of the turbulent length-scale. The ε-equation overpredicts the near wall length-scale under adverse pressure gradient conditions, resulting in too high wall-shear stresses (Rodi and Scheuerer 1986). However, none of the three models (k-<span lang="EN-US"><font face="Symbol">e,</font></span> k-ω and BSL) accounts for the transport of the turbulent shear-stress, which is the main driving force in the momentum exchange. It was demonstrated by Johnson and King (1984) that this is important in non-equilibrium adverse pressure gradient flows. This effect is taken into account with the SST model, which results in a significantly improved representation of the experimental data.
 
[[Image:U3-03d32_files_image017.gif]]
 
<center>Figure 4<nowiki>: Comparison of computed velocity profiles obtained from turbulence models of Set 1 with the experimental data.</nowiki></center>
 
<span lang="EN-US"><font face="&quot;Times New Roman&quot;"><span><font face="Times">Figure 5</font></span> shows the wall shear-stress distribution for the same turbulence models of Set 1. The results are consistent with the velocity profile prediction. However, the wall shear-stress for the k-</font></span>ω and the BLS model is in better agreement with the experimental data than expected from the velocity profiles. This is a result of the proper inner scaling of these models. The only model in the comparison to reproduce the negative wall shear-stress in the separation bubble is the SST model, as expected from the velocity profiles.
 
<center>[[Image:U3-03d32_files_image019.jpg]]</center>
 
<center>Figure 5<nowiki>: Comparison of computed wall shear-stress distributions obtained from turbulence models of Set 1 with the experimental data.</nowiki></center>
 
<span lang="EN-US"><font face="&quot;Times New Roman&quot;"><span><font face="Times">Figure 6</font></span> shows the wall pressure distribution. The flattening of the experimental pressure distribution is a result of the acceleration of the outer flow due to the separation. Clearly, the models, which do not capture the velocity profiles also fail to predict this characteristic of the pressure distribution. </font></span>
 
<center>[[Image:U3-03d32_files_image021.jpg]]</center>
 
<center>Figure 6<nowiki>: Comparison of computed wall pressure distributions obtained from turbulence models of Set 1 with the experimental data.</nowiki></center>
 
<span lang="EN-US"><font face="&quot;Times New Roman&quot;"><span><font face="Times">Figure 7</font></span> shows a comparison of the principal turbulent shear stress <sub>[[Image:U3-03d32_files_image023.gif]]. The two sets of experimental data represent the shear-stress in the coordinate system of the wind tunnel (Carth) and in a streamline coordinate system. These systems are not identical, due to the velocity component normal to the wall. The most important profile is the first profile at (x/D=-0.091 m) as it is just upstream of the experimental separation point. The differences in the velocity profiles are a result of the differences in the turbulent shear stress. From the comparison of the results it is apparent that the tunnel-oriented shear-stress is the relevant quantity. The SST model reproduces this distribution and thereby is able to predict the correct shape of the velocity profiles. The larger the maximum of the predicted shear-stress is, the smaller is the tendency to predict the correct velocity profiles. </sub></font></span>
 
<center>[[Image:U3-03d32_files_image025.gif]]</center>
 
<center>Figure 7 Comparison of computed turbulent shear-stress profiles obtained from turbulence models of Set 1 with the experimental data.</center>
 
<span lang="EN-US"><font face="&quot;Times New Roman&quot;"><span><font face="Times">Figure 8</font></span> shows the maximum of the non-dimensional eddy-viscosity <sub>[[Image:U3-03d32_files_image027.gif]][[Image:U3-03d32_files_image029.gif]]is the maximum value of the eddy-viscosity profile at a given streamwise location, <sub>[[Image:U3-03d32_files_image031.gif]]is the velocity at the edge of the boundary layer and <sub>[[Image:U3-03d32_files_image033.gif]]is the displacement thickness. For equilibrium flows this value should be constant and of the order of 0.017. Figure 8 shows the reduction of this parameter due to the non-equilibrium flow predicted by the models. </sub></sub></sub></font></span>
 
<center>[[Image:U3-03d32_files_image035.gif]]</center>
 
<center>Figure 8 Comparison of maximum non-dimensional eddy-viscosities between Set 1 of turbulence models and experimental data (dotted: Clauser correlation)</center>
 
Figure 9 shows a comparison of the velocity profiles for Set 2 of the turbulence models also computed with INS2D. It shows solutions for the SST , the RNG k-<span lang="EN-US"><font face="Symbol">e (</font></span>Yakhotet al.<span lang="EN-US"><font face="Symbol">1992)</font></span>, the Spalart-Allmaras model and the standard k-ε model. The first observation is that the k-<span lang="EN-US"><font face="Symbol">e </font></span>RNG does not result in a significant improvement over the standard k-ε model. The second observation is that the SA model gives better results than the k-ε models and is close to the experimental data and the SST model. This explains the popularity of this model in aerodynamic solvers. Note however, that the SA model was not designed as a general-purpose model and should be used with caution in non-aerodynamic applications.
 
While for most industrial applications, the differences seen between the SA and the SST model for the CS0 case can be considered marginal, this is not the case for transonic flows with shock-boundary-layer interaction. These flows proved to be highly sensitive to small differences in model performance and even the marginal discrepancies seen in Figure 9 can lead to significant deviations in the pressure distribution between the two models.
 
<center>[[Image:U3-03d32_files_image037.jpg]]</center>
 
<center>Figure 9 Comparison of computed velocity profiles obtained from turbulence models of Set 2 with the experimental data.</center>
 
Figure 10 shows a comparison of the turbulent stress profiles for the same turbulence models of Set 2. Here the differences between the models are more visible and are consistent with the velocity profiles.
 
<center>[[Image:U3-03d32_files_image039.jpg]]</center>
 
<center>Figure 10 Comparison of computed turbulent shear-stress profiles obtained from turbulence models of Set 2 with the experimental data.</center>
 
Figure 11and Figure 12 show the wall shear-stress and the wall pressure distribution for Set 2 of turbulence models. Again, the plots are in agreement with the conclusions from the comparison of the velocity profiles.
 
<center>[[Image:U3-03d32_files_image041.jpg]]</center>
 
<center>Figure 11 Comparison of computed wall shear-stress distributions obtained from turbulence models of Set 2 with the experimental data.</center>
 
<center>[[Image:U3-03d32_files_image043.jpg]]</center>
 
<center>Figure 12 Comparison of computed wall pressure distributions obtained from turbulence models of Set 2 with the experimental data.</center>
 
A final comparison shows the velocity profiles from the k<span lang="EN-US"><font face="Symbol">-e</font></span> and the Spezial-Sarkar-Gatski model (Speziale et al 1991) computed with CFX-5. The last model is a Reynolds Stress model, which solves a separate transport equation for each element of the Reynolds stress tensor. The SSG model clearly gives a better representation of the displacement effect of the boundary layer and the shape of the velocity profiles than the k-ε model. Nevertheless, near the surface, the SSG model deviates from the experimental data and produces a positive wall shear-stress, as can be seen in Figure 13.
 
<center>[[Image:U3-03d32_files_image045.jpg]]</center>
 
<center>Figure 13 Comparison of computed velocity profiles obtained from turbulence models of Set 3 with the experimental data.</center>
 
<center> </center>
 
Despite the significant differences in the outer part of the boundary layer, the SSG model produces almost the same wall shear-stress as the k-ε model.
 
<center>[[Image:U3-03d32_files_image047.jpg]]</center>
 
<center>Figure 14 Comparison of computed wall shear-stress distributions obtained from turbulence models of Set 3 with the experimental data.</center>
 
<center> </center>
 
This is an example where an inconsistency is observed between the displacement effect seen in Figure 9 and the wall shear-stress distribution. The good agreement in the wall pressure distribution indicates that the displacement of the boundary layer has been computed properly, whereas the wall shear-stress distribution gives a fully attached flow field.
 
<center>[[Image:U3-03d32_files_image049.jpg]]</center>
 
<center>Figure 15 Comparison of computed wall pressure distributions obtained from turbulence models of Set 3 with the experimental data.</center>
 
One of the reasons for the unusual behavior of the SSG model near the surface lies most likely in the use of the <span lang="EN-US"><font face="Symbol">e-</font></span>equation in this model. While the SSG model does account properly for the transport of the shear-stresses, the ε-equation overpredicts the turbulent length-scale near the wall. This leads to a momentum exchange near the surface, which locally suppresses the separation bubble.
 
== 1.2              Results from Other Authors ==
 
The following results have been taken from Hellsten and Laine (2000) and Hsu et al. (2003) with permission from AIAA.
 
Figure 16 shows a comparison of velocity profiles computed with a variety of mostly non-linear stress strain models by Hsu et al. (2003). Due to the large number of curves it is difficult to see, but the agreement of the model with the data is significantly better than for the standard k-<span lang="EN-US"><font face="Symbol">e </font></span>model.
 
{| align="left"
| width="12" height="17" |
|-
|
|
[[Image:U3-03d32_files_image051.jpg]]
|}
 
<br clear="ALL" />
 
<center>Figure 16 Velocity profiles computed by Hsu et al (2003)</center>
 
Nevertheless, none of the models gives the correct size of the separation bubble. This can also be seen in Figure 17 and Figure 18 showing the wall shear stress and the wall pressure distributions.
 
{| align="left"
| width="96" height="1" |
|-
|
|
[[Image:U3-03d32_files_image053.jpg]]
|}
 
<br clear="ALL" />
 
<center>Figure 17 Wall shear stress distributions computed by Hsu et al. (2003)</center>
 
{| align="left"
| width="128" height="12" |
|-
|
|
[[Image:U3-03d32_files_image055.jpg]]
|}
 
<br clear="ALL" />
 
<center>Figure 18 Wall pressure distributions computed by Hsu et al. (2003)</center>
 
Only the ω-equation based EASM approaches the zero wall shear stress limit. The other models predict no separation.
 
<span lang="EN-US">Similar conclusions can be drawn from the results of Hellsten and Laine (2000) shown in <span lang="EN-US"><span lang="EN-US">Figure </span></span><span lang="EN-US"><span lang="EN-US">19</span></span> – <span lang="EN-US"><span lang="EN-US">Figure </span></span><span lang="EN-US"><span lang="EN-US">22</span></span>.</span>
 
<center>[[Image:U3-03d32_files_image057.jpg]]</center>
 
<center>Figure 19 Wall pressure distributions computed by Hellsten and Laine (2000)</center>
 
<center> </center>
 
<center>[[Image:U3-03d32_files_image059.jpg]]</center>
 
<center>Figure 20 Wall shear stress distributions computed by Hellsten and Laine (2000)</center>
 
<center>[[Image:U3-03d32_files_image061.jpg]]</center>
 
<center>Figure 21 |Velocity profiles computed by Hellsten and Laine (2000)</center>
 
<center> </center>
 
<center>[[Image:U3-03d32_files_image063.jpg]]</center>
 
<center>Figure 22 Turbulent shear stress profiles computed by Hellsten and Laine (2000)</center>
 
Again, the SST model is giving the best agreement with the data, closely followed by an ω-equation based EASM model.


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

Revision as of 07:51, 8 March 2009


Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References




2D Boundary layers with pressure gradients (A)

Underlying Flow Regime 3-03               © copyright ERCOFTAC 2004


Best Practice Advice

Best Practice Advice for the UFR

This CS0 diffuser testcase is relatively simple to compute concerning the set-up of the flow and the boundary conditions. A grid of 61 x 121 nodes is recommended with equal spacing in the streamwise direction and a stretching in wall normal direction, which is consistent with the wall formulation of the CFD method. The use of a low-Reynolds number grid is recommended if allowed by the code and the wall treatment. The set-up in terms of the boundary conditions is described in chapter 6. The inlet profiles have to be taken from the experimental data, both for the velocity and the turbulence. Bulk profiles are not acceptable, as they alter the boundary layer development and the displacement effects. The outer boundary is to be treated as a slip wall using the contour obtained from the integration of the experimental velocity profiles.

The main issue for the accurate prediction of this case is the selection of the turbulence model. The best solutions, known to the author are from the SST model, the Spalart-Allmaras model and from ω-equation based Reynolds Stress models. Comparable results to the SST model have been obtained with w-equation based Explicit Algebraic Reynolds Stress Models (Hellsten and Laine 2000).

Turbulence modeling is an ever-developing area and it is likely that more and more models will prove in the future to be able to handle this testcase. This will contribute significantly to the reliability of CFD methods for flows with pressure-induced separation.

The CS0 case does not allow judging the ability of turbulence models to properly predict re-attachment of the flow and flow-recovery downstream of a separation zone. It is known from other studies that all RANS turbulence models seem to have significant problems in that area. Particularly models, which predict the correct separation onset, are prone to give larger separation zones and slower flow recovery than indicated by the data. Note that this does not mean that the other models are superior in that respect. It is only due to a cancellation of errors that models, which fail to predict the separation onset, have little problems in the downstream region. The user is advised to keep this deficiency in mind when judging CFD simulations of diffuser flows.

Recommendations for Future Work

Develop model enhancements or new models, which can predict both, flow separation and flow recovery downstream of reattachment. Note that a highly accurate separation prediction capability is essential to make the second step. The improvement cannot be judged with the CS0 case, but the Buice and Eaton (1997) experiments can serve as a guideline.

© copyright ERCOFTAC 2004



Contributors: Florian Menter - AEA Technology


Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References