CFD Simulations AC1-09: Difference between revisions

From KBwiki
Jump to navigation Jump to search
 
(75 intermediate revisions by the same user not shown)
Line 1: Line 1:
{{ACHeader
{{ACHeader
|area=1
|area=1
Line 26: Line 25:


Computations have been performed with different DES-type methods as listed in
Computations have been performed with different DES-type methods as listed in
Table \ref{ref-models}. For the underlying RANS model, the SA model (SA-DES;
[[CFD_Simulations_AC1-09#table1|Table 1]]. For the underlying RANS model, the SA model (SA-DES;
Spalart \emph{et al.}, 1997), the SST
[[CFD_Simulations_AC1-09#8|Spalart ''et al.'', 1997]]), the SST
<math>{k}</math> &ndash; <math>{\omega}</math> model (SST-DES; Travin
<math>{k}</math>&nbsp;&ndash;&nbsp;<math>{\omega}</math> model
\emph{et al.}, 2002), and the TNT
(SST-DES; [[CFD_Simulations_AC1-09#11|Travin&nbsp;''et&nbsp;al.'', 2002]]), and the TNT
<math>{k}</math> &ndash; <math>{\omega}</math> model (X-LES; Kok \emph{et al.}
<math>{k}</math>&nbsp;&ndash;&nbsp;<math>{\omega}</math> model
2004) have been employed. The main difference between SST-DES and X-LES is that
(X-LES; [[CFD_Simulations_AC1-09#2|Kok&nbsp;''et&nbsp;al.'', 2004]])
have been employed. The main difference between SST-DES and X-LES is that
in SST-DES the blended turbulent length scale is only used to define the
in SST-DES the blended turbulent length scale is only used to define the
dissipation term in the <math>{k}</math>-equation, whereas in X-LES it is used to define the
dissipation term in the <math>{k}</math>-equation, whereas in X-LES it is used to define the
dissipation term as well as the eddy-viscosity coefficient. The delayed approach
dissipation term as well as the eddy-viscosity coefficient. The delayed approach
of Spalart \emph{et al.} (2006), shielding attached boundary layers against
of [[CFD_Simulations_AC1-09#9|Spalart&nbsp;''et&nbsp;al.'' (2006)]], shielding
attached boundary layers against
inadvertently switching to LES (so-called shear-stress depletion), is used in
inadvertently switching to LES (so-called shear-stress depletion), is used in
all DES and X-LES computations (denoted as DDES and DX-LES) but one. Finally,
all DES and X-LES computations (denoted as DDES and DX-LES) but one. Finally,
CFSE has employed the improved variant IDDES (Shur \emph{et al.}, 2008). Table
CFSE has employed the improved variant IDDES
\ref{ref-models} also lists the type of solvers and grids.
([[CFD_Simulations_AC1-09#7|Shur&nbsp;''et&nbsp;al.'', 2008]]). [[CFD_Simulations_AC1-09#table1|Table&nbsp;1]] also lists the type of solvers and grids.
 
 
==Flow and Inlet Conditions==
In the present experiments, two flow  conditions  with  different  flow
rates  at  the  particle-laden  primary  inlet  were  considered.  The
resulting maximum gas velocities in the primary jet for the  two  cases
were 12.5 and 7.4&nbsp;m/s, respectively. The flow rate in the annular inlet
was adjusted to give a maximum velocity of about 18&nbsp;m/s.  The  maximum
tangential velocities for both cases were about 13&nbsp;m/s,  corresponding
to a swirl vane angle of 30°. The resulting swirl number was about  0.5
in both cases.
 
The associated mass flow rates for the gas and the particles, the  flow
Reynolds number, the swirl number and other experimental conditions are
listed in [[Test_Data_AC3-12#table1|Table&nbsp;1]].
The mass flow rates of  the  primary  and  secondary
annular jets were calculated from the pressure drops across the orifice
flow meters. The flow Reynolds  number  was  obtained  with  the  total
volume flow rate at the inlet and the outer diameter of the annulus (''D<sub>3</sub>''&nbsp;=&nbsp;64 mm).
The swirl number was calculated as the  ratio  of  the  axial
flux of angular momentum to the axial flux of  linear  momentum, which
was obtained by integration across both the primary and annular inlets.
 
Furthermore, the particle mass flow rates and  the  properties  of  the
glass beads are given in [[Test_Data_AC3-12#table1|Table&nbsp;1]].
The particles have a  smooth  surface
and are spherical in shape. Only less than about 2% of particles  were
non-spherical or fragments, which resulted in small  errors  in  sizing
the beads by the PDA. Such a particle material is ideal for PDA studies
in particulate two-phase-flows. The particle size distribution obtained
by a PDA measurement (18,000 samples)  is  given  in
[[Test_Data_AC3-12#figure3|Figure&nbsp;3]].
Since during the experiment some of the smaller particles were not  collected
in the cyclone separator but were collected in the  paper  filter,  the
particle material was frequently renewed in order to guarantee that the
particles always have the same size distribution. This was  ensured  by
measuring the particle size distribution at  the  inlet  from  time  to
time. The effects of particle damage  could  not  be  observed  in  the
present measurements.


<div id="table1"></div>
<div id="table1"></div>
{| width="600px"
{|border="0" align="center" cellpadding="10"
! style="text-align:left;"|Partner
! style="text-align:left;" width="250px"|Partner
! Model
! style="text-align:left;" width="150px"|Model
! Solver
! style="text-align:left;" | Solver
! Grid
! style="text-align:left;" | Grid
|-
|-
| National Aerospace Laboratory, Amsterdam '''NLR'''
| National Aerospace Laboratory, Amsterdam '''NLR'''
Line 111: Line 73:
| FINE/Open (unstructured)
| FINE/Open (unstructured)
| Common
| Common
|}
<!--
<div id="table1"></div>
{|align="center" border="1" cellpadding="5" width="600px"
|colspan="3" align="center"|'''Table 1. Flow conditions and particle properties'''
|-
|-
! !!Case 1!!Case 2
|align="center" colspan="4"|'''Table 1:''' Models, solvers and grids used by different partners
|-
|colspan="3" style="background-color:lightgray;"|'''''Air flow'''''
|-
|Gas temperature [K]||align="center"|300||align="center"|300
|-
|Gas density [kg/m3]||align="center"|1.18||align="center"|1.18
|-
|Dynamic viscosity [kg/(m s)]||align="center"|18.4&times;10<sup>-6</sup>||align="center"|18.4&times;10<sup>-6</sup>
|-
|Mass flow rate of the primary jet [g/s]||align="center"|9.9||align="center"|6.0
|-
|Mass flow rate of the secondary jet [g/s]||align="center"|38.3||align="center"|44.6
|-
|Ratio of mass flow rates secondary/primary||align="center"|3.87||align="center"|7.43
|-
|Inlet Reynolds number (with D<sub>3</sub> = 64 mm||align="center"|52400||align="center"|54500
|-
|Swirl number||align="center"|0.47||align="center"|0.49
|-
|colspan="3" style="background-color:lightgray;"|'''''particle phase'''''
|-
|Particle mass flow rate [g/s]||align="center"|0.34||align="center"|1.0
|-
|Particle loading in the primary jet||align="center"|0.034||align="center"|0.17
|-
|colspan="3" style="background-color:lightgray;"|'''''particle properties'''''
|-
|Particle number mean diameter [µm]||colspan="2" align="center"|45
|-
|Particle material density [kg/m<sup>3</sup>]||colspan="2" align="center"|2500
|-
|Stokesian particle response time [ms]||colspan="2" align="center"|15.2
|-                      |
|Refractive index||colspan="2" align="center"|1.52
|}
|}
-->




Line 165: Line 83:
the leading edges in the present case, may be significantly delayed. Two
the leading edges in the present case, may be significantly delayed. Two
approaches that improve the prediction of free shear layers have been
approaches that improve the prediction of free shear layers have been
considered for X-LES: a stochastic SGS model (Kok and Van der Ven, 2009) and a
considered for X-LES: a stochastic SGS model
high-pass-filtered (HPF) SGS model (Kok and Van der Ven, 2012). For the present
([[CFD_Simulations_AC1-09#4|Kok and Van der Ven, 2002]]) and a
case, the effect of the HPF model was significant, but the effect of the
high-pass-filtered (HPF) SGS model ([[CFD_Simulations_AC1-09#5|Kok and Van der Ven, 2012]]).
For the present case, the effect of the HPF model was significant, but the effect of the
stochastic model was limited when used in combination with the HPF model.
stochastic model was limited when used in combination with the HPF model.
Therefore, results are mainly presented with only the HPF model. However, some
Therefore, results are mainly presented with only the HPF model. However, some
Line 173: Line 92:
both approaches can be combined with other DES variants as well.
both approaches can be combined with other DES variants as well.


Some numerical details of the different computations are summarized in Table
Some numerical details of the different computations are summarized in
\ref{ref-methods}. It gives the spatial and temporal order of accuracy of the
[[CFD_Simulations_AC1-09#table2|Table 2]]. It gives the spatial and temporal order of accuracy of the
different numerical methods as well as the chosen time-step size and the time
different numerical methods as well as the chosen time-step size and the time
sample used to compute statistical data. The latter two are expressed in
sample used to compute statistical data. The latter two are expressed in
Line 180: Line 99:
free-stream velocity and the mean aerodynamic chord. NLR used a fourth-order
free-stream velocity and the mean aerodynamic chord. NLR used a fourth-order
symmetry-preserving low-dispersion finite-volume scheme with sixth-order
symmetry-preserving low-dispersion finite-volume scheme with sixth-order
artificial dissipation (Kok, 2009). DLR, CASS, and NUM used a standard
artificial dissipation ([[CFD_Simulations_AC1-09#3|Kok, 2009]]).
DLR, CASS, and NUM used a standard
second-order JST-type finite-volume scheme with fourth-order artificial
second-order JST-type finite-volume scheme with fourth-order artificial
dissipation (Jameson \emph{et al.}, 1981), but with a strongly reduced level
dissipation ([[CFD_Simulations_AC1-09#1|Jameson&nbsp;''et&nbsp;al.'', 1981]]),
but with a strongly reduced level
of artificial dissipation.  CFSE used a fourth-order central finite-volume
of artificial dissipation.  CFSE used a fourth-order central finite-volume
scheme with JST-type artificial dissipation.
scheme with JST-type artificial dissipation.
Line 192: Line 113:
the free-stream velocity. The experiment reports a dominant frequency of
the free-stream velocity. The experiment reports a dominant frequency of
<math>{St = f c_\text{mac}/u_{\infty} = 2}</math>
<math>{St = f c_\text{mac}/u_{\infty} = 2}</math>
implying 1,333 time steps per period. Schiavetta
implying 1,333 time steps per period.
\emph{et al.} (2007) also report higher relevant frequencies, e.g., for shear
[[CFD_Simulations_AC1-09#6|Schiavetta&nbsp;''et&nbsp;al'' (2007)]] also report higher relevant frequencies, e.g., for shear
layer instabilities <math>{St = 8}</math> &ndash; <math>{10}</math>.
layer instabilities <math>{St = 8}</math> &ndash; <math>{10}</math>.
These frequencies are also well resolved
These frequencies are also well resolved
by the time step. Note that a time sample of 10 CTU implies that twenty periods
by the time step. Note that a time sample of 10 CTU implies that twenty periods
of the dominant frequency have been captured.
of the dominant frequency have been captured.
<div id="table2"></div>
{|border="0" align="center" border="0" cellpadding="5"
! style="text-align:left;" |Partner
! style="text-align:left;" |Model
! style="text-align:left;" |Space
! style="text-align:left;" |Time
! style="text-align:left;" |Grid
! style="text-align:left;" |Time step (CTU)
! style="text-align:left;" |Time sample (CTU)
|-
| '''NLR'''
| HPF DX-LES
| 4<sup>th</sup>
| 2<sup>nd</sup>
| Common
| 3.75&times;10<sup>-4</sup>
| 11.25
|-
| '''NLR'''
| HPF DX-LES
| 4<sup>th</sup>
| 2<sup>nd</sup>
| Fine
| 2.50&times;10<sup>-4</sup>
| 8.00
|-
| '''NLR'''
| SST-DDES
| 4<sup>th</sup>
| 2<sup>nd</sup>
| Common
| 3.75&times;10<sup>-4</sup>
| 10.50
|-
| '''CASS'''
| SST-DES
| 2<sup>nd</sup>
| 2<sup>nd</sup>
| Own
| 3.65&times;10<sup>-4</sup>
| 14.70
|-
| '''DLR'''
| SA-DDES
| 2<sup>nd</sup>
| 2<sup>nd</sup>
| Modified
| 4.00&times;10<sup>-4</sup>
| 10.00
|-
| '''CFSE'''
| SST-IDDES
| 4<sup>th</sup>
| 2<sup>nd</sup>
| Common
| 3.75&times;10<sup>-4</sup>
| 3.75
|-
| '''NUM'''
| SA-DDES
| 2<sup>nd</sup>
| 2<sup>nd</sup>
| Common
| 9.20&times;10<sup>-4</sup>
| 3.12
|-
|align="center" colspan="7"|'''Table 2:''' Numerical methods used by different partners
|}


==Computational Domain==
==Computational Domain==


The geometries of the wing and the sting are given analytically by Chu and
The geometries of the wing and the sting are given analytically by
Luckring (1996). In the computational domain, the sting is simplified: the shape
[[Description_AC1-09#3|Chu and Luckring (1996)]].
In the computational domain, the sting is simplified: the shape
of the sting upstream of the wing trailing edge conforms to the analytic
of the sting upstream of the wing trailing edge conforms to the analytic
definition, but it is extended as a straight cylinder from the trailing edge of
definition, but it is extended as a straight cylinder from the trailing edge of
Line 231: Line 223:


<div id="figure2"></div>
<div id="figure2"></div>
{|align="center" border="0" width="700"
{|align="center" border="0" width="700" cellspacing="10"
|[[Image:AC1-09_fig2a.png|345px]]||[[Image:AC1-09_fig2b.png|345px]]
|[[Image:AC1-09_fig2a.png|345px]]||[[Image:AC1-09_fig2b.png|345px]]
|-
|align="center"|a) Surface grid|| align="center"|b) Grid in symmetry plane ''y''&nbsp;=&nbsp;0
|-
|-
!align="center" colspan="2"|Figure 2: Impression of common multi-block structured grid around the delta wing (red lines are block edges)
!align="center" colspan="2"|Figure 2: Impression of common multi-block structured grid around the delta wing (red lines are block edges)
Line 271: Line 265:


<div id="figure3"></div>
<div id="figure3"></div>
{|align="center" border="0" width="700"
{|align="center" border="0" width="700" cellspacing="10"
|[[Image:AC1-09_fig3a.png|345px]]||[[Image:AC1-09_fig3b.png|345px]]
|[[Image:AC1-09_fig3a.png|345px]]||[[Image:AC1-09_fig3b.png|345px]]
|-
|align="center"|a) NLR SST-DDES||align="center"|b) NLR HPF DX-LES
|-
|-
!align="center" colspan="2"|Figure 3: Instantaneous isosurface of Q = &Omega;<sup>2</sup> - S<sup>2</sup> (Q = 365 u<sup>2</sup><sub>&infin;</sub>/&nbsp;c<sub>r</sub><sup>2</sup>), coloured with vorticity magnitude &Omega;
!align="center" colspan="2"|Figure 3: Instantaneous isosurface of Q = &Omega;<sup>2</sup> - S<sup>2</sup> (Q = 365 u<sup>2</sup><sub>&infin;</sub>/&nbsp;c<sub>r</sub><sup>2</sup>), coloured with vorticity magnitude &Omega;
Line 291: Line 287:


<div id="figure4"></div>
<div id="figure4"></div>
{|align="center" border="0" width="700"
{|align="center" border="0" width="700" cellspacing="10"
|[[Image:AC1-09_fig4a.png|345px]]||[[Image:AC1-09_fig4b.png|345px]]
|[[Image:AC1-09_fig4a.png|345px]]||[[Image:AC1-09_fig4b.png|345px]]
|-
|align="center"|a) Mean pressure coefficient||align="center"|b) RMS of pressure coefficient
|-
|-
|[[Image:AC1-09_fig4c.png|345px]]||[[Image:AC1-09_fig4d.png|345px]]
|[[Image:AC1-09_fig4c.png|345px]]||[[Image:AC1-09_fig4d.png|345px]]
|-
|-
!align="center" colspan="2"|Figure 4: Solutions on common grid and fine grid of HPF DX-LES at station x/c<sub>r</sub> = 0.4 (Note: subfigures a) and b) contain the solutions on both the left side and the right side of the wing)
|align="center"|c) Resolved turbulent kinetic energy <math>{k/u_\infty^2}</math>on common grid
|align="center"|d) Resolved turbulent kinetic energy <math>{k/u_\infty^2}</math>on fine grid
|-
!align="center" colspan="2"|Figure 4: Solutions on common grid and fine grid of HPF DX-LES at station x/c<sub>r</sub>&nbsp;=&nbsp;0.4 (Note: subfigures a) and b) contain the solutions on both the left side and the right side of the wing)
|}
|}




<div id="figure5"></div>
<div id="figure5"></div>
{|align="center" border="0" width="700"
{|align="center" border="0" width="700" cellspacing="10"
|[[Image:AC1-09_fig5a.png|345px]]||[[Image:AC1-09_fig5b.png|345px]]
|[[Image:AC1-09_fig5a.png|345px]]||[[Image:AC1-09_fig5b.png|345px]]
|-
|align="center"|a) Mean pressure coefficient
|align="center"|b) RMS of pressure coefficient
|-
|-
|[[Image:AC1-09_fig5c.png|345px]]||[[Image:AC1-09_fig5d.png|345px]]
|[[Image:AC1-09_fig5c.png|345px]]||[[Image:AC1-09_fig5d.png|345px]]
|-
|-
!align="center" colspan="2"|Figure 5: Solutions on common grid and non-isotropic grid for Stochastic HPF X-LES at station x/c<sub>r</sub> = 0.2 (Note: subfigures a) and b) contain the solutions on both the left side and the right side of the wing)
|align="center"|c) Resolved turbulent kinetic energy <math>{k/u_\infty^2}</math> with far field at 3c<sub>r</sub>
|align="center"|d) Resolved turbulent kinetic energy <math>{k/u_\infty^2}</math> with far field at 15c<sub>r</sub>
|-
!align="center" colspan="2"|Figure 5: Solutions on common grid with different locations of far field for Stochastic HPF X-LES at station x/c<sub>r</sub>&nbsp;=&nbsp;0.4 (Note: subfigures a) and b) contain the solutions on both the left side and the right side of the wing)
|}
|}


Line 316: Line 323:
vortex breakdown location. On this grid, a significantly lower level of resolved
vortex breakdown location. On this grid, a significantly lower level of resolved
turbulence was obtained in the initial vortex, leading to a significantly
turbulence was obtained in the initial vortex, leading to a significantly
stronger vortex and influencing the solution further downstream (Figure
stronger vortex and influencing the solution further downstream ([[CFD_Simulations_AC1-09#figure6|Figure 6]]).
\ref{fig-griddist-dependence}).




<div id="figure6"></div>
<div id="figure6"></div>
{|align="center" border="0" width="700"
{|align="center" border="0" width="700" cellspacing="10"
|[[Image:AC1-09_fig6a.png|345px]]||[[Image:AC1-09_fig6b.png|345px]]
|[[Image:AC1-09_fig6a.png|345px]]||[[Image:AC1-09_fig6b.png|345px]]
|-
|align="center"|a) Mean pressure coefficient
|align="center"|b) RMS of pressure coefficient
|-
|-
|[[Image:AC1-09_fig6c.png|345px]]||[[Image:AC1-09_fig6d.png|345px]]
|[[Image:AC1-09_fig6c.png|345px]]||[[Image:AC1-09_fig6d.png|345px]]
|-
|align="center"|c) Resolved turbulent kinetic energy <math>{k/u_\infty^2}</math> on common grid
|align="center"|d) Resolved turbulent kinetic energy <math>{k/u_\infty^2}</math> on non-isotropic grid
|-
|-
!align="center" colspan="2"|Figure 6: Solutions on common grid and non-isotropic grid for Stochastic HPF X-LES at station x/c<sub>r</sub> = 0.2 (Note: subfigures a) and b) contain the solutions on both the left side and the right side of the wing)
!align="center" colspan="2"|Figure 6: Solutions on common grid and non-isotropic grid for Stochastic HPF X-LES at station x/c<sub>r</sub> = 0.2 (Note: subfigures a) and b) contain the solutions on both the left side and the right side of the wing)
Line 334: Line 346:
Another question is how long the time sample should be to obtain statistical
Another question is how long the time sample should be to obtain statistical
convergence. Using the stochastic X-LES method, a range of time samples has been
convergence. Using the stochastic X-LES method, a range of time samples has been
considered, from 6.6 CTU to 21 CTU, see Figure \ref{fig-samples-dependence}.
considered, from 6.6 CTU to 21 CTU, see [[CFD_Simulations_AC1-09#figure7|Figure 7]].
(For the smallest time sample, 13.2 periods of the dominant frequency are still
(For the smallest time sample, 13.2 periods of the dominant frequency are still
captured.) For the mean flow quantities, statistical convergence is obtained
captured.) For the mean flow quantities, statistical convergence is obtained
Line 349: Line 361:


<div id="figure7"></div>
<div id="figure7"></div>
{|align="center" border="0" width="700"
{|align="center" border="0" width="700" cellspacing="10"
|[[Image:AC1-09_fig7a.png|345px]]||[[Image:AC1-09_fig7b.png|345px]]
|[[Image:AC1-09_fig7a.png|345px]]||[[Image:AC1-09_fig7b.png|345px]]
|-
|align="center"|a) Mean pressure coefficient
|align="center"|b) RMS of pressure coefficient
|-
|-
!align="center" colspan="2"|Figure 7: Solutions with varying time samples for Stochastic HPF X-LES at station x/c<sub>r</sub> = 0.4 (Note: the figures contain the solutions on both the left side and the right side of the wing)
!align="center" colspan="2"|Figure 7: Solutions with varying time samples for Stochastic HPF X-LES at station x/c<sub>r</sub> = 0.4 (Note: the figures contain the solutions on both the left side and the right side of the wing)

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

Solution Strategy

Detached Eddy Simulations (DES) have been performed by a number of partners in the EU-project ATAAC. DES (Spalart, 2009) is a hybrid RANS–LES approach that was originally based on the Spalart–Allmaras (SA) model, but that has also been extended to models. In DES switching between LES and RANS is effectively achieved by defining the turbulence length scale employed in the turbulence model as the minimum of the LES and RANS length scales. The LES length scale is the filter width, which is defined as the maximum of the mesh size in all three computational directions at each grid point. The RANS length scale depends on the RANS model employed: essentially the wall distance for the SA model and for the model. Note that for the SA model the RANS length scale is static and therefore the RANS–LES interface is fixed, whereas for the model the RANS–LES interface is dynamic.

Computations have been performed with different DES-type methods as listed in Table 1. For the underlying RANS model, the SA model (SA-DES; Spalart et al., 1997), the SST  –  model (SST-DES; Travin et al., 2002), and the TNT  –  model (X-LES; Kok et al., 2004) have been employed. The main difference between SST-DES and X-LES is that in SST-DES the blended turbulent length scale is only used to define the dissipation term in the -equation, whereas in X-LES it is used to define the dissipation term as well as the eddy-viscosity coefficient. The delayed approach of Spalart et al. (2006), shielding attached boundary layers against inadvertently switching to LES (so-called shear-stress depletion), is used in all DES and X-LES computations (denoted as DDES and DX-LES) but one. Finally, CFSE has employed the improved variant IDDES (Shur et al., 2008). Table 1 also lists the type of solvers and grids.

Partner Model Solver Grid
National Aerospace Laboratory, Amsterdam NLR HPF DX-LES & SST-DDES ENSOLV (structured) Common & fine
Airbus Defence and Space, München CASS SST-DES TAU (unstructured) Own
German Aerospace Center, Göttingen DLR SA-DDES TAU (unstructured) Modified
Computational Fluids and Structures Engineering, Lausanne CFSE SST-IDDES NSMB (structured) Common
NUMECA International, Brussels NUM SA-DDES FINE/Open (unstructured) Common
Table 1: Models, solvers and grids used by different partners


One of the difficulties of DES consists of the development of resolved turbulence as the method switches from RANS to LES. In particular, the development of instabilities in free shear layers, such as those separating from the leading edges in the present case, may be significantly delayed. Two approaches that improve the prediction of free shear layers have been considered for X-LES: a stochastic SGS model (Kok and Van der Ven, 2002) and a high-pass-filtered (HPF) SGS model (Kok and Van der Ven, 2012). For the present case, the effect of the HPF model was significant, but the effect of the stochastic model was limited when used in combination with the HPF model. Therefore, results are mainly presented with only the HPF model. However, some of the sensitivity studies were performed with the stochastic model. Note that both approaches can be combined with other DES variants as well.

Some numerical details of the different computations are summarized in Table 2. It gives the spatial and temporal order of accuracy of the different numerical methods as well as the chosen time-step size and the time sample used to compute statistical data. The latter two are expressed in convective time units (CTU), i.e., they are made dimensionless with the free-stream velocity and the mean aerodynamic chord. NLR used a fourth-order symmetry-preserving low-dispersion finite-volume scheme with sixth-order artificial dissipation (Kok, 2009). DLR, CASS, and NUM used a standard second-order JST-type finite-volume scheme with fourth-order artificial dissipation (Jameson et al., 1981), but with a strongly reduced level of artificial dissipation. CFSE used a fourth-order central finite-volume scheme with JST-type artificial dissipation.

For time integration, all partners used the second-order backward implicit scheme. A time step of CTU and a mesh width of (see grid characteristics below) imply a convective CFL number of approximately based on the free-stream velocity. The experiment reports a dominant frequency of implying 1,333 time steps per period. Schiavetta et al (2007) also report higher relevant frequencies, e.g., for shear layer instabilities . These frequencies are also well resolved by the time step. Note that a time sample of 10 CTU implies that twenty periods of the dominant frequency have been captured.


Partner Model Space Time Grid Time step (CTU) Time sample (CTU)
NLR HPF DX-LES 4th 2nd Common 3.75×10-4 11.25
NLR HPF DX-LES 4th 2nd Fine 2.50×10-4 8.00
NLR SST-DDES 4th 2nd Common 3.75×10-4 10.50
CASS SST-DES 2nd 2nd Own 3.65×10-4 14.70
DLR SA-DDES 2nd 2nd Modified 4.00×10-4 10.00
CFSE SST-IDDES 4th 2nd Common 3.75×10-4 3.75
NUM SA-DDES 2nd 2nd Common 9.20×10-4 3.12
Table 2: Numerical methods used by different partners

Computational Domain

The geometries of the wing and the sting are given analytically by Chu and Luckring (1996). In the computational domain, the sting is simplified: the shape of the sting upstream of the wing trailing edge conforms to the analytic definition, but it is extended as a straight cylinder from the trailing edge of the wing down to the far-field boundary (whereas the analytical definition of the sting diameter increases at some distance from the trailing edge).

NLR has generated a common multi-block structured grid, consisting of 22 blocks and 6.3 million grid cells. The grid has a conical structure over a large part of the wing: the grid covering the main vortex is essentially isotropic at each chord-wise station (outside the boundary layer) and the mesh width grows in all directions (including the stream-wise direction) together with the main vortex, going from approximately to . In other words, the grid resolution relative to the main vortex is kept constant. Only in a small region near the apex, the conical structure is not fully maintained, avoiding a grid singularity. The first grid cell at the wing surface has a height of . The far-field boundary is located at three root chord lengths from the wing in all directions, i.e., given that the apex is at the origin, the extent of the computational domain is given by , , and . To study grid sensitivity, also a finer grid with the mesh width reduced by a factor 2/3 in all directions (21.4 million grid cells) as well as a modified grid with the far field located at 15 root chord lengths (8.7 million grid cells) have been generated. Additionally, CASS has generated a fine, unstructured grid (17.4 million grid cells, far field at 15 root chord lengths) that also essentially has a conical structure. An impression of the common grid is given in Figure 2.


AC1-09 fig2a.png AC1-09 fig2b.png
a) Surface grid b) Grid in symmetry plane y = 0
Figure 2: Impression of common multi-block structured grid around the delta wing (red lines are block edges)

Boundary Conditions

At the wing and sting surfaces, a no-slip boundary condition is applied for the velocity and an adiabatic boundary condition is applied for the temperature. At the far field, a uniform flow is applied according the given free-stream conditions. Two far field locations have been considered, as indicated above.

Application of Physical Models

As indicated above, DES computations are applied with either the SA model, the SST model, or the TNT model as underlying (RANS) turbulence model. No laminar-turbulent transition is applied, i.e., the boundary layers are assumed to be fully turbulent.

An impression of the flow is given in Figure 3, showing the instantaneous vortical structures using the Q-criterion. The main vortex above the wing is essentially captured with LES. The RANS part of the computation concerns the attached boundary layers and the near-wall region below the main vortex.

To illustrate the difficulty of capturing the instabilities developing in the free shear layers, Figure 3 shows the flow structures for a standard DES computation (SST-DDES) and for a computation with the high-pass filtered SGS model (HPF DX-LES). In the standard DES computation, the leading-edge shear layers are stable over approximately the first half of the wing (the helical structures that are visible are steady). Only towards the trailing edge, the shear layer becomes unstable and some large-scale unsteady structures are formed. In contrast, including the HPF SGS model allows the shear layer to become unstable close to the apex of the wing. Subsequently, much finer turbulent structures are formed in the shear layers and the main vortex.


AC1-09 fig3a.png AC1-09 fig3b.png
a) NLR SST-DDES b) NLR HPF DX-LES
Figure 3: Instantaneous isosurface of Q = Ω2 - S2 (Q = 365 u2/ cr2), coloured with vorticity magnitude Ω

Numerical Accuracy

For the X-LES computations of NLR, several sensitivity studies have been performed. For LES computations in general, grid convergence is a troublesome matter; the best one can do is to show the sensitivity of the result to the grid resolution. The results on the fine grid (mesh width and time step both reduced by a factor 2/3) were found to lie close to the results on the common grid, as illustrated in Figure 4 for station (other stations show similar results). Furthermore, the sensitivity to the location of the far field, placed at three or fifteen root chord lengths from the wing, was found to be weak, as illustrated in Figure 5. Therefore, most computations have been performed with the nearest far field (smallest grid size).


AC1-09 fig4a.png AC1-09 fig4b.png
a) Mean pressure coefficient b) RMS of pressure coefficient
AC1-09 fig4c.png AC1-09 fig4d.png
c) Resolved turbulent kinetic energy on common grid d) Resolved turbulent kinetic energy on fine grid
Figure 4: Solutions on common grid and fine grid of HPF DX-LES at station x/cr = 0.4 (Note: subfigures a) and b) contain the solutions on both the left side and the right side of the wing)


AC1-09 fig5a.png AC1-09 fig5b.png
a) Mean pressure coefficient b) RMS of pressure coefficient
AC1-09 fig5c.png AC1-09 fig5d.png
c) Resolved turbulent kinetic energy with far field at 3cr d) Resolved turbulent kinetic energy with far field at 15cr
Figure 5: Solutions on common grid with different locations of far field for Stochastic HPF X-LES at station x/cr = 0.4 (Note: subfigures a) and b) contain the solutions on both the left side and the right side of the wing)


A strong sensitivity was found, however, to the grid point distribution. Initially, a grid had been generated with a uniform mesh width in -direction (equal to ), while it was conical in the other two directions. This grid was not isotropic, except in the vicinity of the expected vortex breakdown location. On this grid, a significantly lower level of resolved turbulence was obtained in the initial vortex, leading to a significantly stronger vortex and influencing the solution further downstream (Figure 6).


AC1-09 fig6a.png AC1-09 fig6b.png
a) Mean pressure coefficient b) RMS of pressure coefficient
AC1-09 fig6c.png AC1-09 fig6d.png
c) Resolved turbulent kinetic energy on common grid d) Resolved turbulent kinetic energy on non-isotropic grid
Figure 6: Solutions on common grid and non-isotropic grid for Stochastic HPF X-LES at station x/cr = 0.2 (Note: subfigures a) and b) contain the solutions on both the left side and the right side of the wing)



Another question is how long the time sample should be to obtain statistical convergence. Using the stochastic X-LES method, a range of time samples has been considered, from 6.6 CTU to 21 CTU, see Figure 7. (For the smallest time sample, 13.2 periods of the dominant frequency are still captured.) For the mean flow quantities, statistical convergence is obtained with the smallest time sample (up to plotting accuracy). For the fluctuation quantities (in particular the RMS of ), a larger sample of about 15 CTU is needed. The smallest time sample, however, is also acceptable, because the differences with the statistically converged results are small relative to the differences found between the partners (as shown below) and the differences with the experiment. Note that the time sample should not include the transient in which the flow is still developing. The length of this transient has not been extensively investigated, but should be checked for each computation separately as it will depend on the initial solution.


AC1-09 fig7a.png AC1-09 fig7b.png
a) Mean pressure coefficient b) RMS of pressure coefficient
Figure 7: Solutions with varying time samples for Stochastic HPF X-LES at station x/cr = 0.4 (Note: the figures contain the solutions on both the left side and the right side of the wing)

CFD Results

The CFD results consist of the following data:

  • Time-averaged and RMS surface pressure coefficient at five chord stations:
    , , , , and .
  • Time-averaged velocity field and resolved turbulent kinetic energy field at five cross planes:
    , , , , and .

References

  1. A. Jameson, W. Schmidt, and E. Turkel (1981) Numerical Solutions of the Euler Equations by Finite Volume Methods Using Runge--Kutta Time-Stepping Schemes, AIAA Paper 1981-1259.
  2. J. C. Kok, H. S. Dol, B. Oskam, and H. van der Ven (2004) Extra-Large Eddy Simulation of Massively Separated Flows, AIAA paper 2004-264.
  3. J. C. Kok (2009) A high-order low-dispersion symmetry-preserving finite-volume method for compressible flow on curvilinear grids, Journal of Computational Physics, 228:6811–6832. NLR-TP-2008-775.
  4. J. C. Kok and H. van der Ven (2009) Destabilizing free shear layers in X-LES using a stochastic subgrid-scale model. In S. H. Peng and P. Doerffer and W. Haase, editors, Progress in Hybrid RANS–LES Modelling, volume 111 of Notes on Numerical Fluid Mechanics and Multidisciplinary Design, pages 179--189. Springer. NLR-TP-2009-327.
  5. J. C. Kok and H. van der Ven (2012) Capturing free shear layers in hybrid RANS--LES simulations of separated flow, In Third Symposium `Simulation of Wing and Nacelle Stall', 21–22 June Braunschweig, Germany. NLR-TP-2012-333.
  6. L. A. Schiavetta, K. J. Badcock, R. M. and Cummings (2007) Comparison of DES and URANS for unsteady vortical flows over delta wings, AIAA Paper 2007-1085.
  7. M. L. Shur, P. R. Spalart, M. Kh. Strelets, and A. K. Travin (2008) A hybrid RANS—LES approach with delayed-DES and wall-modelled LES capabilities, International Journal of Heat and Fluid Flow, 29:1638--1649.
  8. P. R. Spalart, W.-H. Jou, M. Strelets, and S. R. Allmaras (1997) Comments on the feasibility of LES for wings, and on a hybrid RANS/LES approach. In C. Liu and Z. Liu, editors, Advances in DNS/LES, Proc.\ 1st AFOSR Int.\ Conf.\ on DNS/LES, 1997, Ruston (LA), USA.<?div>
  9. P. R. Spalart, S. Deck, M. L. Shur, K. D. Squires, M. Kh. Strelets, and A. Travin (2006) A new version of detached-eddy simulation, resistant to ambiguous grid densities. Theor. Comp. Fluid Dyn. 20:181–95.
  10. P. R. Spalart (2009) Detached-Eddy Simulation, Annual Review of Fluid Mechanics, 41:181–202.
  11. A. Travin, M. Shur, M. Strelets, and P. R. Spalart (2002) Physical and numerical upgrades in the detached-eddy simulation of complex turbulent flows. In R. Friedrich and W. Rodi, editors, Advances in LES of Complex Flows, pages 239–354. Kluwer Academic Publishers. Proc.\ of the Euromech Colloquium 412, 2000, Munich, Germany.
  12. A. K. Travin, M. L. Shur, P. R. Spalart, and M. Kh. Strelets (2006) Improvement of delayed detached-eddy simulation for LES with wall modelling. In P. Wesseling and E. Oñate and J. Périaux, editors, ECCOMAS CFD 2006, 5–8 September, Egmond aan Zee, The Netherlands.




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