UFR 3-36 Best Practice Advice: Difference between revisions

From KBwiki
Jump to navigation Jump to search
No edit summary
No edit summary
 
(25 intermediate revisions by 3 users not shown)
Line 1: Line 1:
=HiFi-TURB-DLR rounded step=
=HiFi-TURB-DLR rounded step=


{{UFRHeaderLib
{{UFRHeader
|area=3
|area=3
|number=36
|number=36
Line 13: Line 13:
== Key Physics ==
== Key Physics ==
The key flow physics to be accurately captured in this UFR are turbulent boundary layers subjected to an adverse pressure gradient over a curved surface with and without separation and reattachment.  
The key flow physics to be accurately captured in this UFR are turbulent boundary layers subjected to an adverse pressure gradient over a curved surface with and without separation and reattachment.  
== Numerical Modelling ==
== Numerical Modelling and Boundary Conditions ==
With the FV-Code TAU a second order discretization for convective fluxes was used including the turbulence model. For grid details please refer to the previous section ("[[Lib:UFR 3-36 Test Case#CFD Methods|CFD Methods]]"). Otherwise, standard rules of grid generation apply.
With the DLR-TAU Code a second order discretization for convective fluxes was used including the turbulence model. For grid details please refer to the previous section ("[[UFR 3-36 Test Case#CFD Methods|CFD Methods]]"). Otherwise, standard rules of grid generation apply.
To ensure the upper boundary will not influence the flow over the curved wall, a study was performed to find the optimal position for the far-field boundary condition. It was found that a distance to the curved wall less than <math>{179 H}</math> shows an influence on the computed pressure distributions. With a distance larger than <math>{179 H}</math> the influence vanishes. These results were generated during the initial design stage of the testcase when a cheaper negative Spalart-Allmaras turbulence model (SA-neg) was applied. Results of this study are given in [[Lib:UFR_3-36_Test_Case#figure7|Fig. 7]]. Very minor deviations (not visible) only vanished for a distance larger than <math>{179 H}</math>.
 
To ensure that the upper boundary will not influence the flow over the curved wall, a study was performed to find the optimal position for the far-field boundary. It was found that a distance to the curved wall less than <math>{179 H}</math> shows an influence on the computed pressure distributions. For distances larger than <math>{179 H}</math> the influence vanishes. These results were generated during the initial design stage of the testcase when a cheaper Spalart-Allmaras turbulence model ( the SA-neg variant) was applied. Results of this study are given in [[UFR_3-36_Test_Case#figure7|Fig. 7]]. Minor deviations (not visible) only vanished for a distance larger than <math>{179 H}</math>.


<div id="figure7"></div>
<div id="figure7"></div>
Line 23: Line 24:
!align="center" colspan="2"|Figure 7: Pressure distributions for varying distances of the far-field boundary condition  
!align="center" colspan="2"|Figure 7: Pressure distributions for varying distances of the far-field boundary condition  
|}
|}
Besides studying the far-field position, three simulations were performed to study the sensitivity of the inlet position: the upstream section of the curved step was extended, a Dirichlet boundary condition was applied and laminar-turbulent transition was considered applying a RANS-based transition transport model. The results showed no significant effect and did not change the incipient character of the flow, however, an influence of the predicted skin friction coefficient distribution of the SSG/LRR-<math>\omega</math> model in the recovery region was observed. As the RANS computations were performed to design the testcase, no DNS data were initially available to evaluate this influence.


== Physical Modelling ==
== Physical Modelling ==
The test case for this numerical experiment needs to be well-defined and all boundary conditions as well as the flow conditions at suitable reference positions need to be specified. The reason behind this is to produce a reliable basis for subsequent DNS computations of the test case. Therefore, a study of the grid space was conducted. Three key features were investigated: the far-field position and the inlet position and laminar-turbulent transition. Each of these features was varied to ensure they do not affect the flow behavior.
In the early design phase of the testcase, five different RANS turbulence models already released in the DLR-TAU code and one additional model modification were used for assessment of the flow.
 
Furthermore, five RANS turbulence models already released in TAU and one additional model modification were used for the simulations:


* Four variations of the one-equation Spalart-Allmaras RANS turbulence model:  
* Four variations of the one-equation Spalart-Allmaras RANS turbulence model:  
** SA-neg [&#8204;[[Lib:UFR_3-36_References#9|9]]]
** SA-neg [&#8204;[[UFR_3-36_References#9|9]]]
** Non-linear SA-neg model version with rotation and curvature correction (RC) and quadratic constitutive relation (QCR) [&#8204;[[Lib:UFR_3-36_References#10|10]]][&#8204;[[Lib:UFR_3-36_References#11|11]]]: SA-RC-QCR  
** Non-linear SA-neg model version with rotation and curvature correction (RC) and quadratic constitutive relation (QCR) [&#8204;[[UFR_3-36_References#10|10]]][&#8204;[[UFR_3-36_References#11|11]]]: SA-RC-QCR  
** SA model version with Low-Re-modifications [&#8204;[[Lib:UFR_3-36_References#12|12]]][&#8204;[[Lib:UFR_3-36_References#13|13]]]: SA-LRe  
** SA model version with Low-Re-modifications [&#8204;[[UFR_3-36_References#12|12]]][&#8204;[[UFR_3-36_References#13|13]]]: SA-LRe  
** Non-linear SA-neg model version with rotation and curvature correction (RC) and quadratic constitutive relation (QCR) as well as Low-Re-modifications: SA-RC-QCR-LRe  
** Non-linear SA-neg model version with rotation and curvature correction (RC) and quadratic constitutive relation (QCR) as well as Low-Re-modifications: SA-RC-QCR-LRe  
*Two variations of the two-equation Menter SST RANS turbulence model:  
*Two variations of the two-equation Menter SST RANS turbulence model:  
**SST2003 [&#8204;[[Lib:UFR_3-36_References#14|14]]]
**SST2003 [&#8204;[[UFR_3-36_References#14|14]]]
**Non-linear model version SST2003-RC-QCR [&#8204;[[Lib:UFR_3-36_References#15|15]]][&#8204;[[Lib:UFR_3-36_References#11|11]]]
**Non-linear model version SST2003-RC-QCR [&#8204;[[UFR_3-36_References#15|15]]][&#8204;[[UFR_3-36_References#11|11]]]
*The seven-equation omega-based Differential Reynolds stress turbulence model SSG/LRR-<math>\omega</math> [&#8204;[[Lib:UFR_3-36_References#16|16]]] including the length scale correction [&#8204;[[Lib:UFR_3-36_References#17|17]]]
*The seven-equation omega-based Differential Reynolds stress turbulence model SSG/LRR-<math>\omega</math> [&#8204;[[UFR_3-36_References#16|16]]] including the length scale correction [&#8204;[[UFR_3-36_References#17|17]]]
 
For the detailed analysis of the flow only the SSG/LRR-<math>\omega</math> model was evaluated, results for the other RANS turbulence models are not given in detail. An indication of the performance of these models is displayed in [[UFR_3-36_Test_Case#figure8|Fig. 8]].


<div id="figure8"></div>
<div id="figure8"></div>
Line 43: Line 46:
|[[Image:UFR3-36_fig8a.png|370px]]||[[Image:UFR3-36_fig8b.png|370px]]
|[[Image:UFR3-36_fig8a.png|370px]]||[[Image:UFR3-36_fig8b.png|370px]]
|-
|-
!align="center" colspan="2"|Figure 8: Skin friction and pressure coefficient distributions with different RANS turbulence models  
!align="center" colspan="2"|Figure 8: Skin friction (left) and pressure coefficient (right) distributions for different RANS turbulence models  
|}
|}


== Application Uncertainties ==
Considering the skin friction coefficient distribution, the SA-neg turbulence model shows very good agreement with the DNS results, while the corrections RC + QCR outlined above move the flow closer to separation. Adding the LowRe modifiction yields results in between the original SA-neg and the SA-neg-RC-QCR model. Predictions of the Menter-SST turbulence model are close to separation, adding the corrections RC + QCR yields separated flow. The SSG/LRR-<math>\omega</math> model shows a deviation in the recovery region downstream of the step which might be improved by extending the extent of the upstream section of the curved step.
The Reynolds numbers used for the test case were determined, on the one hand, low enough to be affordable for DNS computations and on the other hand, high enough to avoid re-laminarization in the accelerated region of the boundary layer upstream the APG region. The latter was ensured by applying two different criteria for possible re-laminarization, the acceleration parameter
 
In the ("[[UFR 3-36 Test Case#Evaluation|Evaluation]]") section the results of the DNS are compared in detail to the results of the SSG/LRR-<math>\omega</math> model of DLR and the Wilcox k-<math>\omega</math> model of UniBG. From the pressure coefficient distribution it is visible that the Wilcox k-<math>\omega</math> model overpredicts the suction and pressure peak compared to the DNS while the SSG/LRR-<math>\omega</math> model shows reverse behaviour but is closer to the DNS. The skin friction coefficient reveals a stronger tendency of the SSG/LRR-<math>\omega</math> model to separation compared with the DNS and the Wilcox k-<math>\omega</math> model. Downstream of the step both models deviate in a similar manner from the DNS, however, at x/H = 5 the SSG/LRR-<math>\omega</math> model is in better agreement with the DNS data. The streamwise velocity profiles reflect the results of the skin friction coefficient: the SSG/LRR-<math>\omega</math> model predicts a smaller velocity in the APG region between 0 < x/H < 5 compared to the DNS and here the Wilcox k-<math>\omega</math> model is closer to DNS, but the former model is in better agreement with DNS at 5 < x/H < 9. For the streamwise normal stresses the SSG/LRR-<math>\omega</math> model is not able to capture the characteristics visible in the DNS data at <math>0.001 \leq y_W/H < 0.01</math> upstream of the APG region, while these stresses are overpredicted close to the boundary-layer edge downstream of the APG region. For the shear-stresses, the SSG/LRR-<math>\omega</math> model better captures the DNS data compared to the Wilcox k-<math>\omega</math> model. Finally, for the turbulent kinetic energy, both RANS models show very similar distributions compared to the DNS.
 
Based on these results there is no clear indication as to which RANS turbulence model performs better for the testcase. Both models do not capture the flow in the APG-region correctly compared to the DNS data. From a design point of view, the results of the SSG/LRR-<math>\omega</math> model are more conservative as the predicted flow is closer to separation compared to the DNS data. Again, it needs to be emphasized that the SSG/LRR-<math>\omega</math> model was used to design the testcase and the results were generated in a "blind" manner i.e. the target DNS results were not known. For future computations of this testcase it is recommended to increase the extent of the upstream section of the flow domain when applying the SSG/LRR-<math>\omega</math> model.
 
== Application uncertainties ==
The Reynolds numbers for the testcase were determined, on the one hand, low enough to be affordable for DNS computations and on the other hand, high enough to avoid re-laminarization in the accelerated region of the turbulent boundary layer upstream of the APG region. The latter was ensured by applying two different criteria for possible re-laminarization, the acceleration parameter


<math>
<math>
Line 59: Line 68:
</math>  
</math>  


and compared to limiting values for re-laminarization according to the literature [&#8204;[[Lib:UFR_3-36_References#24|24]]], [&#8204;[[Lib:UFR_3-36_References#25|25]]] and [&#8204;[[Lib:UFR_3-36_References#26|26]]]. Here, <math> U_\tau </math> is the friction velocity. [[Lib:UFR_3-36_Test_Case#figure8|Fig. 8]] shows the distribution of <math> K_{acc} </math> and <math> \Delta_p </math> for the Reynolds number <math>Re_H = 78,490 </math> together with the threshold values for both parameters, which this Reynolds number does not exceed.
and compared to limiting values for re-laminarization according to the literature [&#8204;[[UFR_3-36_References#24|24]]], [&#8204;[[UFR_3-36_References#25|25]]] and [&#8204;[[UFR_3-36_References#26|26]]]. Here, <math> U_\tau </math> is the friction velocity. [[UFR_3-36_Test_Case#figure8|Fig. 9]] shows the distribution of <math> K_{acc} </math> and <math> \Delta_p </math> for the Reynolds number <math>Re_H = 78,490 </math> together with the threshold values for both parameters, not exceeded for this Reynolds number.


<div id="figure8"></div>
<div id="figure9"></div>
{|align="center" border="0" width="500"
{|align="center" border="0" width="500"
|[[Image:UFR3-36_fig8.png|500px]]
|[[Image:UFR3-36_fig8.png|500px]]
|-
|-
!align="center"|Figure 8: Applied criteria for re-laminarization   
!align="center"|Figure 9: Assessment of possible re-laminarization   
|}
|}


== Recommendations for Future Work ==
== Recommendations for Future Work ==
This test case was designed to be part of a family of test cases with varying pressure gradient and Reynolds number to provide a comprehensive data base for data-driven turbulence modeling. As soon as DNS data for the other configurations are available, these cases should be added to the test case description and evaluation.
This testcase was designed to be part of a family of testcases with varying pressure gradient and Reynolds number to provide a comprehensive data base for data-driven turbulence modeling. As soon as DNS data for the other configurations are available, these cases should be added to the test case description and RANS computations for these testcases should be added. Additionally, RANS turbulence models augmented by a Machine Learning method that were generated and tested on this testcase during the HiFi-TURB project should be described, their computational results added and evaluated in comparison to the traditional models.
<br/>
<br/>
----
----
Line 76: Line 85:
|organisation=Deutsches Luft-und Raumfahrt Zentrum (DLR)
|organisation=Deutsches Luft-und Raumfahrt Zentrum (DLR)
}}
}}
{{UFRHeaderLib
{{UFRHeader
|area=3
|area=3
|number=36
|number=36

Latest revision as of 15:55, 17 February 2023

HiFi-TURB-DLR rounded step

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References

Semi-confined flows

Underlying Flow Regime 3-36

Best Practice Advice

Key Physics

The key flow physics to be accurately captured in this UFR are turbulent boundary layers subjected to an adverse pressure gradient over a curved surface with and without separation and reattachment.

Numerical Modelling and Boundary Conditions

With the DLR-TAU Code a second order discretization for convective fluxes was used including the turbulence model. For grid details please refer to the previous section ("CFD Methods"). Otherwise, standard rules of grid generation apply.

To ensure that the upper boundary will not influence the flow over the curved wall, a study was performed to find the optimal position for the far-field boundary. It was found that a distance to the curved wall less than shows an influence on the computed pressure distributions. For distances larger than the influence vanishes. These results were generated during the initial design stage of the testcase when a cheaper Spalart-Allmaras turbulence model ( the SA-neg variant) was applied. Results of this study are given in Fig. 7. Minor deviations (not visible) only vanished for a distance larger than .

UFR3-36 fig7a.png UFR3-36 fig7b.png
Figure 7: Pressure distributions for varying distances of the far-field boundary condition

Besides studying the far-field position, three simulations were performed to study the sensitivity of the inlet position: the upstream section of the curved step was extended, a Dirichlet boundary condition was applied and laminar-turbulent transition was considered applying a RANS-based transition transport model. The results showed no significant effect and did not change the incipient character of the flow, however, an influence of the predicted skin friction coefficient distribution of the SSG/LRR- model in the recovery region was observed. As the RANS computations were performed to design the testcase, no DNS data were initially available to evaluate this influence.

Physical Modelling

In the early design phase of the testcase, five different RANS turbulence models already released in the DLR-TAU code and one additional model modification were used for assessment of the flow.

  • Four variations of the one-equation Spalart-Allmaras RANS turbulence model:
    • SA-neg [‌9]
    • Non-linear SA-neg model version with rotation and curvature correction (RC) and quadratic constitutive relation (QCR) [‌10][‌11]: SA-RC-QCR
    • SA model version with Low-Re-modifications [‌12][‌13]: SA-LRe
    • Non-linear SA-neg model version with rotation and curvature correction (RC) and quadratic constitutive relation (QCR) as well as Low-Re-modifications: SA-RC-QCR-LRe
  • Two variations of the two-equation Menter SST RANS turbulence model:
    • SST2003 [‌14]
    • Non-linear model version SST2003-RC-QCR [‌15][‌11]
  • The seven-equation omega-based Differential Reynolds stress turbulence model SSG/LRR- [‌16] including the length scale correction [‌17]

For the detailed analysis of the flow only the SSG/LRR- model was evaluated, results for the other RANS turbulence models are not given in detail. An indication of the performance of these models is displayed in Fig. 8.

UFR3-36 fig8a.png UFR3-36 fig8b.png
Figure 8: Skin friction (left) and pressure coefficient (right) distributions for different RANS turbulence models

Considering the skin friction coefficient distribution, the SA-neg turbulence model shows very good agreement with the DNS results, while the corrections RC + QCR outlined above move the flow closer to separation. Adding the LowRe modifiction yields results in between the original SA-neg and the SA-neg-RC-QCR model. Predictions of the Menter-SST turbulence model are close to separation, adding the corrections RC + QCR yields separated flow. The SSG/LRR- model shows a deviation in the recovery region downstream of the step which might be improved by extending the extent of the upstream section of the curved step.

In the ("Evaluation") section the results of the DNS are compared in detail to the results of the SSG/LRR- model of DLR and the Wilcox k- model of UniBG. From the pressure coefficient distribution it is visible that the Wilcox k- model overpredicts the suction and pressure peak compared to the DNS while the SSG/LRR- model shows reverse behaviour but is closer to the DNS. The skin friction coefficient reveals a stronger tendency of the SSG/LRR- model to separation compared with the DNS and the Wilcox k- model. Downstream of the step both models deviate in a similar manner from the DNS, however, at x/H = 5 the SSG/LRR- model is in better agreement with the DNS data. The streamwise velocity profiles reflect the results of the skin friction coefficient: the SSG/LRR- model predicts a smaller velocity in the APG region between 0 < x/H < 5 compared to the DNS and here the Wilcox k- model is closer to DNS, but the former model is in better agreement with DNS at 5 < x/H < 9. For the streamwise normal stresses the SSG/LRR- model is not able to capture the characteristics visible in the DNS data at upstream of the APG region, while these stresses are overpredicted close to the boundary-layer edge downstream of the APG region. For the shear-stresses, the SSG/LRR- model better captures the DNS data compared to the Wilcox k- model. Finally, for the turbulent kinetic energy, both RANS models show very similar distributions compared to the DNS.

Based on these results there is no clear indication as to which RANS turbulence model performs better for the testcase. Both models do not capture the flow in the APG-region correctly compared to the DNS data. From a design point of view, the results of the SSG/LRR- model are more conservative as the predicted flow is closer to separation compared to the DNS data. Again, it needs to be emphasized that the SSG/LRR- model was used to design the testcase and the results were generated in a "blind" manner i.e. the target DNS results were not known. For future computations of this testcase it is recommended to increase the extent of the upstream section of the flow domain when applying the SSG/LRR- model.

Application uncertainties

The Reynolds numbers for the testcase were determined, on the one hand, low enough to be affordable for DNS computations and on the other hand, high enough to avoid re-laminarization in the accelerated region of the turbulent boundary layer upstream of the APG region. The latter was ensured by applying two different criteria for possible re-laminarization, the acceleration parameter

as well as the pressure-gradient parameter

and compared to limiting values for re-laminarization according to the literature [‌24], [‌25] and [‌26]. Here, is the friction velocity. Fig. 9 shows the distribution of and for the Reynolds number together with the threshold values for both parameters, not exceeded for this Reynolds number.

UFR3-36 fig8.png
Figure 9: Assessment of possible re-laminarization

Recommendations for Future Work

This testcase was designed to be part of a family of testcases with varying pressure gradient and Reynolds number to provide a comprehensive data base for data-driven turbulence modeling. As soon as DNS data for the other configurations are available, these cases should be added to the test case description and RANS computations for these testcases should be added. Additionally, RANS turbulence models augmented by a Machine Learning method that were generated and tested on this testcase during the HiFi-TURB project should be described, their computational results added and evaluated in comparison to the traditional models.



Contributed by: Erij Alaya and Cornelia Grabe — Deutsches Luft-und Raumfahrt Zentrum (DLR)

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References


© copyright ERCOFTAC 2024