Evaluation AC1-09: Difference between revisions

From KBwiki
Jump to navigation Jump to search
m (Dave.Ellacott moved page SilverP:Evaluation AC1-09 to Evaluation AC1-09)
 
(11 intermediate revisions by the same user not shown)
Line 73: Line 73:
<div id="figure9"></div>
<div id="figure9"></div>
{|align="center" border="0" width="700"
{|align="center" border="0" width="700"
!Experiment
|[[Image:AC1-09_fig9a.png|230px]]||[[Image:AC1-09_fig9b.png|230px]]||[[Image:AC1-09_fig9c.png|230px]]
|[[Image:AC1-09_fig9a.png|230px]]||[[Image:AC1-09_fig9b.png|230px]]||[[Image:AC1-09_fig9c.png|230px]]
|-
|-
!NLR HPF DX-LES
|[[Image:AC1-09_fig9d.png|230px]]||[[Image:AC1-09_fig9e.png|230px]]||[[Image:AC1-09_fig9f.png|230px]]
|[[Image:AC1-09_fig9d.png|230px]]||[[Image:AC1-09_fig9e.png|230px]]||[[Image:AC1-09_fig9f.png|230px]]
|-
|-
!CASS SST-DES
|[[Image:AC1-09_fig9g.png|230px]]||[[Image:AC1-09_fig9h.png|230px]]||[[Image:AC1-09_fig9i.png|230px]]
|[[Image:AC1-09_fig9g.png|230px]]||[[Image:AC1-09_fig9h.png|230px]]||[[Image:AC1-09_fig9i.png|230px]]
|-
|-
!DLR SA-DDES
|[[Image:AC1-09_fig9j.png|230px]]||[[Image:AC1-09_fig9k.png|230px]]||[[Image:AC1-09_fig9l.png|230px]]
|[[Image:AC1-09_fig9j.png|230px]]||[[Image:AC1-09_fig9k.png|230px]]||[[Image:AC1-09_fig9l.png|230px]]
|-
|-
!CFSE SST-IDDES
|[[Image:AC1-09_fig9m.png|230px]]||[[Image:AC1-09_fig9n.png|230px]]||[[Image:AC1-09_fig9o.png|230px]]
|[[Image:AC1-09_fig9m.png|230px]]||[[Image:AC1-09_fig9n.png|230px]]||[[Image:AC1-09_fig9o.png|230px]]
|-
|-
!NUM SA-DDES
|[[Image:AC1-09_fig9p.png|230px]]||[[Image:AC1-09_fig9q.png|230px]]||[[Image:AC1-09_fig9r.png|230px]]
|[[Image:AC1-09_fig9p.png|230px]]||[[Image:AC1-09_fig9q.png|230px]]||[[Image:AC1-09_fig9r.png|230px]]
|-
|-
!NLR SST-DDES
|[[Image:AC1-09_fig9s.png|230px]]||[[Image:AC1-09_fig9t.png|230px]]||[[Image:AC1-09_fig9u.png|230px]]
|[[Image:AC1-09_fig9s.png|230px]]||[[Image:AC1-09_fig9t.png|230px]]||[[Image:AC1-09_fig9u.png|230px]]
|-
|-
!align="center" colspan="3"|Figure 9: Resolved turbulent kinetic energy at three cross planes (note the logarithmic scale)
!align="center" colspan="4"|Figure 9: Resolved turbulent kinetic energy at three cross planes (note the logarithmic scale)
|}
|}


Line 97: Line 104:
the resolved turbulent kinetic energy, but still the details of the distribution
the resolved turbulent kinetic energy, but still the details of the distribution
differ from the experiment. There are also a few unexpected results: the DLR
differ from the experiment. There are also a few unexpected results: the DLR
SA-DDES computation shows significant RMS levels near the symmetry plane ($\eta
SA-DDES computation shows significant RMS levels near the symmetry plane <math>{\left.\eta=0\right.}</math>
= 0$) and on the lower side of the wing. The CFSE SST-IDDES shows the same
and on the lower side of the wing. The CFSE SST-IDDES shows the same
behaviour for the first two stations and also has very high RMS levels below
behaviour for the first two stations and also has very high RMS levels below
the vortex at <math>{\left.x/c_r = 0.2\right.}</math>, which then drop to very low levels at <math>{\left.x/c_r = 0.4\right.}</math>.
the vortex at <math>{\left.x/c_r = 0.2\right.}</math>, which then drop to very low levels at <math>{\left.x/c_r = 0.4\right.}</math>.
Line 115: Line 122:


The occurrence of vortex breakdown can be seen in the time-averaged
The occurrence of vortex breakdown can be seen in the time-averaged
velocity field (Figure \ref{fig-vel}).
velocity field ([[Evaluation_AC1-09#figure11|Figure 11]]).


<div id="figure11"></div>
<div id="figure11"></div>
{|align="center" border="0" width="700"
{|align="center" border="0" width="700"
!Experiment
|[[Image:AC1-09_fig11a.png|230px]]||[[Image:AC1-09_fig11b.png|230px]]||[[Image:AC1-09_fig11c.png|230px]]
|[[Image:AC1-09_fig11a.png|230px]]||[[Image:AC1-09_fig11b.png|230px]]||[[Image:AC1-09_fig11c.png|230px]]
|-
|-
!NLR HPF DX-LES
|[[Image:AC1-09_fig11d.png|230px]]||[[Image:AC1-09_fig11e.png|230px]]||[[Image:AC1-09_fig11f.png|230px]]
|[[Image:AC1-09_fig11d.png|230px]]||[[Image:AC1-09_fig11e.png|230px]]||[[Image:AC1-09_fig11f.png|230px]]
|-
|-
!CASS SST-DES
|[[Image:AC1-09_fig11g.png|230px]]||[[Image:AC1-09_fig11h.png|230px]]||[[Image:AC1-09_fig11i.png|230px]]
|[[Image:AC1-09_fig11g.png|230px]]||[[Image:AC1-09_fig11h.png|230px]]||[[Image:AC1-09_fig11i.png|230px]]
|-
|-
!DLR SA-DDES
|[[Image:AC1-09_fig11j.png|230px]]||[[Image:AC1-09_fig11k.png|230px]]||[[Image:AC1-09_fig11l.png|230px]]
|[[Image:AC1-09_fig11j.png|230px]]||[[Image:AC1-09_fig11k.png|230px]]||[[Image:AC1-09_fig11l.png|230px]]
|-
|-
!CFSE SST-IDDES
|[[Image:AC1-09_fig11m.png|230px]]||[[Image:AC1-09_fig11n.png|230px]]||[[Image:AC1-09_fig11o.png|230px]]
|[[Image:AC1-09_fig11m.png|230px]]||[[Image:AC1-09_fig11n.png|230px]]||[[Image:AC1-09_fig11o.png|230px]]
|-
|-
!NUM SA-DDES
|[[Image:AC1-09_fig11p.png|230px]]||[[Image:AC1-09_fig11q.png|230px]]||[[Image:AC1-09_fig11r.png|230px]]
|[[Image:AC1-09_fig11p.png|230px]]||[[Image:AC1-09_fig11q.png|230px]]||[[Image:AC1-09_fig11r.png|230px]]
|-
|-
!NUM SST-DDES
|[[Image:AC1-09_fig11s.png|230px]]||[[Image:AC1-09_fig11t.png|230px]]||[[Image:AC1-09_fig11u.png|230px]]
|[[Image:AC1-09_fig11s.png|230px]]||[[Image:AC1-09_fig11t.png|230px]]||[[Image:AC1-09_fig11u.png|230px]]
|-
|-
!align="center" colspan="3"|Figure 11: Time-averaged velocity field at three cross planes (contour plot of velocity normal to plane; vector plot of velocity components in plane)
!align="center" colspan="4"|Figure 11: Time-averaged velocity field at three cross planes (contour plot of velocity normal to plane; vector plot of velocity components in plane)
|}
|}


Line 151: Line 165:


The occurrence of secondary separation can also, in principle, be seen in the
The occurrence of secondary separation can also, in principle, be seen in the
time-averaged velocity field (Figure \ref{fig-vel}). For all computations, the
time-averaged velocity field ([[Evaluation_AC1-09#figure11|Figure 11]]). For all computations, the
time-averaged velocity vectors are parallel to the surface below the main
time-averaged velocity vectors are parallel to the surface below the main
vortex. Following the velocity vectors towards the leading edge, they point
vortex. Following the velocity vectors towards the leading edge, they point
Line 160: Line 174:
Therefore, no conclusion can be drawn about secondary separation in the
Therefore, no conclusion can be drawn about secondary separation in the
experiment from the time-averaged velocity field. The distribution of the
experiment from the time-averaged velocity field. The distribution of the
turbulent kinetic energy (Figure \ref{fig-tke}), however, does suggest a shear
turbulent kinetic energy ([[Evaluation_AC1-09#figure9|Figure 9]]), however, does suggest a shear
layer separating from the surface between the main vortex and the leading edge,
layer separating from the surface between the main vortex and the leading edge,
indicating secondary separation in the experiment as well.
indicating secondary separation in the experiment as well.

Latest revision as of 15:19, 11 February 2017


Front Page

Description

Test Data

CFD Simulations

Evaluation

Best Practice Advice

Vortex breakdown above a delta wing with sharp leading edge

Application Challenge AC1-09   © copyright ERCOFTAC 2024

Comparison of Test Data and CFD

An overview is given of the computational results obtained by the different partners.

In terms of the mean pressure coefficient, the different computations are reasonably clustered in the vicinity of the experimental data (Figure 8). The largest difference with the experiment is found at the first station (), where all computations show a clear suction peak due to the main vortex (and some also a smaller peak due to secondary separation), whereas the experiment shows a plateau. In the computations, the main vortex starts to develop immediately from the apex. In the experiment, however, it seems that the vortex starts further downstream along the leading edge, more resembling the vortex development of a round leading edge. Possibly, this may be related to the relative bluntness of the leading-edge angle. There may also be a Reynolds number effect, because at a higher Reynolds number (), the experiment does show a clear suction peak (Furman and Breitsamter, 2009). At the other four stations, the computed pressure distributions compare reasonably well with the experiment in terms of the location of the main suction peak and the level of the pressure plateau outboard of the main peak, indicating no or only weak secondary separation. The level of the suction peak is underpredicted at the station where the experiment shows vortex breakdown, but the computations do not (see below). At the other stations, the suction peak in most computations is close to or just below the experiment.

AC1-09 fig8a.png AC1-09 fig8b.png
AC1-09 fig8c.png AC1-09 fig8d.png
AC1-09 fig8e.png
Figure 8: Time-averaged pressure coefficient at five chord-wise stations

The level of the main suction peak, which indicates the strength of the main vortex, depends on the level of resolved turbulence. If the turbulence is underresolved, then the vortex will be too strong; if the level of turbulence is too high, then the vortex will be too weak. Hence, the initial development of the vortex and turbulence in the vortex is vital. With a standard DES method, the shear layer emanating from the leading edge is stable over a significant part of the wing (Figure 3a), so that there is insufficient resolved turbulence, leading to a strong, compact, main vortex. Including the high-pass-filtered (HPF) SGS model, the level of the subgrid stresses is significantly reduced and the shear layer becomes unstable close to the apex. Much higher levels of resolved turbulence are then obtained, weakening the main vortex and more in-line with the experiment. Consequently, the resolved turbulent kinetic energy compares fairly well with the experiment, at all three cross planes (see Figure 9). Only the peak level is overpredicted (consistent with the underprediction of the suction peak). In particular, the experiment shows a clear local minimum at the core of the vortex from station and further downstream, whereas the computation still shows a clear peak at and a plateau or a slight local minimum further downstream. Other results with DDES and IDDES methods compare less favourably with the experiment and also show significant variation between the different computations. In particular for the upstream sections, the DDES results show low levels of resolved turbulence (some more than others, most likely depending on the numerical method and in particular the level of numerical dissipation), as explained above.

Experiment AC1-09 fig9a.png AC1-09 fig9b.png AC1-09 fig9c.png
NLR HPF DX-LES AC1-09 fig9d.png AC1-09 fig9e.png AC1-09 fig9f.png
CASS SST-DES AC1-09 fig9g.png AC1-09 fig9h.png AC1-09 fig9i.png
DLR SA-DDES AC1-09 fig9j.png AC1-09 fig9k.png AC1-09 fig9l.png
CFSE SST-IDDES AC1-09 fig9m.png AC1-09 fig9n.png AC1-09 fig9o.png
NUM SA-DDES AC1-09 fig9p.png AC1-09 fig9q.png AC1-09 fig9r.png
NLR SST-DDES AC1-09 fig9s.png AC1-09 fig9t.png AC1-09 fig9u.png
Figure 9: Resolved turbulent kinetic energy at three cross planes (note the logarithmic scale)

In terms of the RMS of the pressure coefficient, a large variation between the results is found (Figure 10). A number of results strongly overpredict the level, while some others initially show a very low level (correlating to a low level of resolved turbulence). Only the HPF DX-LES computation essentially has the right level at all stations, consistent with the resolved turbulent kinetic energy, but still the details of the distribution differ from the experiment. There are also a few unexpected results: the DLR SA-DDES computation shows significant RMS levels near the symmetry plane and on the lower side of the wing. The CFSE SST-IDDES shows the same behaviour for the first two stations and also has very high RMS levels below the vortex at , which then drop to very low levels at . This makes the reliability of these two computations questionable.

AC1-09 fig10a.png AC1-09 fig10b.png
AC1-09 fig10c.png AC1-09 fig10d.png
AC1-09 fig10e.png
Figure 10: RMS of pressure coefficient at five chord-wise stations

The occurrence of vortex breakdown can be seen in the time-averaged velocity field (Figure 11).

Experiment AC1-09 fig11a.png AC1-09 fig11b.png AC1-09 fig11c.png
NLR HPF DX-LES AC1-09 fig11d.png AC1-09 fig11e.png AC1-09 fig11f.png
CASS SST-DES AC1-09 fig11g.png AC1-09 fig11h.png AC1-09 fig11i.png
DLR SA-DDES AC1-09 fig11j.png AC1-09 fig11k.png AC1-09 fig11l.png
CFSE SST-IDDES AC1-09 fig11m.png AC1-09 fig11n.png AC1-09 fig11o.png
NUM SA-DDES AC1-09 fig11p.png AC1-09 fig11q.png AC1-09 fig11r.png
NUM SST-DDES AC1-09 fig11s.png AC1-09 fig11t.png AC1-09 fig11u.png
Figure 11: Time-averaged velocity field at three cross planes (contour plot of velocity normal to plane; vector plot of velocity components in plane)

In the experiment, a very low axial velocity is present in the vortex core at , indicating that the initial vortex breakdown has taken place well before this location. All computations clearly show vortex breakdown at , but not at , with exception of the DLR SA-DDES result. This result shows an initial sign of vortex breakdown at , but not as strongly as the experiment, indicating the initial vortex breakdown is also downstream of the experiment. Thus, it appears that none of the computations obtain the same vortex breakdown location as the experiment, even though there is a significant variation in for example resolved turbulence between the computations.

The occurrence of secondary separation can also, in principle, be seen in the time-averaged velocity field (Figure 11). For all computations, the time-averaged velocity vectors are parallel to the surface below the main vortex. Following the velocity vectors towards the leading edge, they point away from the surface, which indicates secondary separation. For the experiment, the flow pattern is less clear. Note in particular that the velocity field of the experiment is not tangential to the wing surface below the primary vortex, indicating some error in the measurements or postprocessing. Therefore, no conclusion can be drawn about secondary separation in the experiment from the time-averaged velocity field. The distribution of the turbulent kinetic energy (Figure 9), however, does suggest a shear layer separating from the surface between the main vortex and the leading edge, indicating secondary separation in the experiment as well.

Based on the available results, no definitive conclusion can be drawn about which basic DES-type method (SA-DES, SST-DES, or X-LES) is most suitable for this test case. All partners have used different models in combination with different solvers. Thus, the influence of the numerical methods on the results may still be significant. Also, statistical convergence has not been verified for all results. One can conclude, however, that the X-LES method with the HPF SGS model shows the best comparison with the experiment in terms of resolved turbulence and pressure fluctuations. This can be attributed not to the X-LES method itself, but to the HPF SGS model, which strongly improves the initial development of resolved turbulence in the shear layer. This modification of the SGS model could also be used in other DES-type methods. Using the stochastic variant of X-LES does not significantly alter the results once the HPF SGS model has been included. The location of vortex breakdown is closest to the experiment for the DLR SA-DDES computation; this computation, however, appears to be less reliable when looking at the resolved turbulent kinetic energy and in particular the RMS levels of pressure. Finally, the delayed modification of DES (or X-LES) does not appear to have a strong effect on the results for this case and on this grid.




Contributed by: J.C. Kok, H. van der Ven (National Aerospace Laboratory NLR Amsterdam, The Netherlands), E. Tangermann (Airbus Defence and Space München, Germany), S. Sanchi (Computational Fluids and Structures Engineering Lausanne, Switzerland), A. Probst and K.A. Weinman (German Aerospace Center DLR Göttingen, Germany), L. Temmerman (NUMECA International Brussels, Belgium) — '

Front Page

Description

Test Data

CFD Simulations

Evaluation

Best Practice Advice


© copyright ERCOFTAC 2024