UFR 4-19 Evaluation

From KBwiki
Jump to navigation Jump to search

Converging-diverging transonic diffuser

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References

Confined flows

Underlying Flow Regime 4-19

Evaluation

Comparison of CFD Calculations with Experiments

For both "weak" and "strong" Mach number cases,for comparison with the calculations, experimental data for four measurements stations after the shock-wave (h1, h2, h3, h4) are selected for the longitudinal velocity as well as pressure distributions along the bottom and the top wall of the diffuser. A representative figure of the shock-wave positions and the measurement stations is shown in fig.9.The x-location of the axial measurement stations are non-dimensionalized with the diffuser throat height,resulting in the following four stations: h1 = 2.822, h2 = 4.611, h3 = 6.340 and h4 = 7.493.

UFR4-19 Fig9.png
Figure 9: Shock-waves positions and experimental measurements

The weak Mach number case

The modelling results of the three adopted turbulence models are compared with the available experimental data for the "weak" Mach number case in figs.10 and 11.

All three turbulent models are able to capture the maximum velocity value. The main differences are observed in the boundary layer in the wall regions. The EVMs have similar behavior as they underpredict the axial velocity,especially in the top wall region of the diverging part of the diffuser. On the other hand,the RSM overpredicts the axial velocity on both walls providing thinner boundary layers. The pressure coefficient distributions along both walls of the diffuser are plotted in fig.11.

All the turbulence models capture the pressure coefficient distribution in the converging and the diverging part of the diffuser with a very good accuracy.However,in the region near the shock-wave (the region where the minimum pressure value is computed),the good level of accuracy is not maintained. The EVMs compute a less intense shock-wave,as it can be concluded from fig.11. Regarding the RSM,it computes a shock-wave with a greater value for the maximum Mach number,giving a lower pressure coefficient value.

The static pressure (Pa) and the Mach number contours for all the turbulence models are presented in figs.12 and 13. The flow field representation is almost the same for all the models with the RSM providing the largest Mach number value in comparison with the EVMs.

UFR4-19 Fig12.jpg
Figure 12: Weak shock wave case: Contours of static pressure
UFR4-19 Fig13.jpg
Figure 13: Weak shock wave case: Contours of Mach number

The strong Mach number case

The computed velocity distributions for the EVMs and RSM in comparison with the experimental data are shown in fig.14.

UFR4-19 Fig14a.png|UFR4-19 Fig14b.png
Figure 14: Strong shock-wave case: Velocity distributions

From the first two measurements stations it is evident that the experimental data show a recirculation region formed after the ?strong? shock-wave on the top wall of the diverging part of the diffuser. From fig.14 it can be seen that the EVMs are able to capture the boundary layer separation but with a smaller negative velocity resulting in a smaller recirculation region. Additionally,the LEVM computes also a small recirculation zone on the bottom wall, which has an impact on the development of the attached boundary layer further downstream. As a result,the attached boundary layers of the LEVM are thicker in comparison with the computed boundary layers of the other two turbulence models and in relation to the experiment,fig.14. For the first station at x/hthroat = 2.882 only the NLEVM is able to capture the velocity maximum in the core of the flow. However, it gives a lower value of the maximum Mach number as can be concluded from the pressure-coefficient distributions shown in fig.15. Additionally, the NLEVM is the model that presents the greatest accuracy in the calculation of the boundary layer on the bottom wall of the diffuser. Finally, the RSM is not able to capture the flow recirculation on the upper wall, but its results are closer to the experimental data on the bottom wall and in comparison with the LEVM. It must be mentioned that the RSM’s inability to model the boundary layer separation is probably related to the minor modification that has been introduced. Although this modification to the RSM provided numerical stability and results more consistent with the flow physics regarding the intensity of the Mach number, it provided non-physical results regarding the recirculation region since it was not able to predict it as shown in fig.18. Hence, the RSM computation for the “strong” Mach number case with the model in its current form is not appropriate for reaching to general conclusions regarding the performance of RSM models for this flow. Further studies on the modelling of the pressure-strain correlation term are needed.

The pressure coefficient distributions for the bottom and the top diverging diffuser wall are presented in fig.15. The NLEVM predicts a higher value of the coefficient and a lower value at the Mach number that was measured. On the other hand, the RSM and the LEVM compute pressure coefficient values slightly lower than the measured ones for the top wall and almost exact for the bottom wall in comparison with the experiments.


The contours of the static pressure (Pa) distribution and the Mach number for the strong shock-wave case are shown in figs.16 and 17. The NLEVM predicts smaller Mach number values and has a similar behavior as in the weak shock-wave case, in comparison to the other two turbulence models.


UFR4-19 Fig16.jpg
Figure 16: Strong shock wave case: Contours of static pressure
UFR4-19 Fig17.jpg
Figure 17: Strong shock wave case: Contours of Mach number

Table 4 presents the measured and the modelled separation and reattachment point on the diverging part of the diffuser. Additionally, streamlines focusing on the recirculation region are depicted in fig.18, which also show the separation and reattachment points


From table 4 it can be concluded that the EVMs underpredict the total length of the recirculation region. The LEVM computes a smaller recirculation region since the separation point is located further downstream than the measurements show. The NLEVM predictions are much closer to the experiment for the locations of both separation and reattachment points. The RSM, as adopted in the current study, is not able to produce any negative values for the axial velocity and thus computes the flow to be fully attached on the top wall as can be seen from the streamlines in fig.16. As a concluding remark for all the adopted models, the results obtained with the NLEVM for the predicted onset and end of the recirculation region are closer to the experimental data.

UFR4-19 Fig18.png
Figure 18: Streamlines for the strong Mach number case

Turbulent kinetic energy and Reynolds-stresses distribution

For the turbulence quantities calculated by the turbulence models, e.g. Reynolds-stresses and turbulent kinetic energy distributions, there are no experimental data available. Hence, in order to assess the more advanced turbulence model behavior regarding the turbulence development in the diffuser transonic flow, a comparison of the normal Reynolds-stress distributions along the same selected locations is presented only for the computed turbulence quantities. The Reynolds-stress values are normalized with the maximum value of the longitudinal stress computed by the RSM, in order to have a qualitative relative comparison since the RSM computes higher longitudinal Reynolds-stress than the EVMs. The normalized normal Reynolds stress distributions i.e. (u'u') ̅, (v'v') ̅and (w'w') ̅, for the three adopted turbulence models are presented in fig.19 for the “weak” Mach number case. For the “strong” Mach number case they are given in VY. As a general observation, the RSM provides larger values for the normal Reynolds-stresses. Furthermore, the Reynolds-stress anisotropy is clearly better represented by the RSM. The NLEVM can also capture partly the anisotropy (higher values for the longitudinal stress) but yields (v'v') ̅≈(w'w') ̅ as discussed further below. Finally, as is well known, LEVM cannot predict the anisotropy as it basically gives a value of 2/3k for each normal stress component and should not be used for calculating the individual components.





Contributed by: Z. Vlahostergios, K. Yakinthos — Aristotle University of Thessaloniki, Greece

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References


© copyright ERCOFTAC 2024