UFR 3-34 Test Case: Difference between revisions
m (Dave.Ellacott moved page Lib:UFR 3-34 Test Case to UFR 3-34 Test Case over redirect) |
|||
(72 intermediate revisions by 3 users not shown) | |||
Line 8: | Line 8: | ||
= Test Case = | = Test Case = | ||
==Brief Description of the Study Test Case== | ==Brief Description of the Study Test Case== | ||
A 3D sketch of the experimental setup is shown in Fig. 1. The configuration presents a | A 3D sketch of the experimental setup is shown in [[UFR_3-34_Test_Case#figure1|Fig. 1]]. | ||
The configuration presents a Glauert-Goldschmied type body consisting of a relatively long fore body and a relatively short concave | |||
ramp comprising the aft part of the model mounted between two glass endplate frames with both | ramp comprising the aft part of the model mounted between two glass endplate frames with both | ||
leading edge and trailing edges faired smoothly with a wind tunnel splitter plate. | leading edge and trailing edges faired smoothly with a wind tunnel splitter plate. | ||
Line 19: | Line 20: | ||
|} | |} | ||
Major geometric and flow parameters of the TC are presented in Fig.2 and summarized in | Major geometric and flow parameters of the TC are presented in [[UFR_3-34_Test_Case#figure2|Fig. 2]] | ||
Table 1 (note that in the “baseline” experiment considered here the slot shown in Fig.2 was | and summarized in [[UFR_3-34_Test_Case#table1|Table 1]] | ||
closed). | (note that in the “baseline” experiment considered here the slot shown in | ||
[[UFR_3-34_Test_Case#figure2|Fig. 2]] was closed). | |||
<div id="figure2"></div> | <div id="figure2"></div> | ||
Line 27: | Line 29: | ||
|[[Image:UFR3-34_Fig2.png|613px]] | |[[Image:UFR3-34_Fig2.png|613px]] | ||
|- | |- | ||
|align="center"|'''Figure 2:''' Schematic of TC geometry | |align="center"|'''Figure 2:''' Schematic of TC geometry and flow parameters | ||
|} | |} | ||
Line 52: | Line 54: | ||
(Case 3), at https://turbmodels.larc.nasa.gov/nasahump_val.html, and also enter the ERCOFTAC | (Case 3), at https://turbmodels.larc.nasa.gov/nasahump_val.html, and also enter the ERCOFTAC | ||
Classic Collection http://cfd.mace.manchester.ac.uk/ercoftac under number C.83. The data set | Classic Collection http://cfd.mace.manchester.ac.uk/ercoftac under number C.83. The data set | ||
includes: streamwise distributions of the surface pressure and skin-friction coefficients, | includes: streamwise distributions of the surface pressure and skin-friction coefficients, | ||
<math>{C_P}</math> and <math>{C_f}</math>, | |||
the region 0.63 < x/c < 1.39. | and mean velocity and Reynolds stresses fields in the tunnel center-plane, roughly covering | ||
the region 0.63 < ''x''/''c'' < 1.39. | |||
== Test Case Experiments == | == Test Case Experiments == | ||
A photo of the experimental rig is presented in Fig.3. | A photo of the experimental rig is presented in [[UFR_3-34_Test_Case#figure3|Fig. 3]]. | ||
<div id="figure3"></div> | <div id="figure3"></div> | ||
Line 68: | Line 71: | ||
The experiments were performed in the NASA Langley 20′′×28′′ shear flow tunnel. The flow was | The experiments were performed in the NASA Langley 20′′×28′′ shear flow tunnel. The flow was | ||
nominally 2D, although with side-wall effects (3D flow) expected near the endplates. The | nominally 2D, although with side-wall effects (3D flow) expected near the endplates. The | ||
reference “chord” length of the model, c, is defined as the length of the hump on the wall and is | reference “chord” length of the model, ''c'', is defined as the length of the hump on the wall and is | ||
equal to 420mm. The maximum thickness of the hump, h, is equal to 53.7mm. | equal to 420mm. The maximum thickness of the hump, ''h'', is equal to 53.7mm. | ||
The model was equipped with 153 centre-span static pressure ports and 20 dynamic pressure | The model was equipped with 153 centre-span static pressure ports and 20 dynamic pressure | ||
ports in the vicinity of the separated flow region. Sixteen spanwise pressure ports were located | ports in the vicinity of the separated flow region. Sixteen spanwise pressure ports were located | ||
on the fore body (x/c = 0.19) and on the ramp at (x/c = 0.86). Two-dimensional PIV data were | on the fore body (''x''/''c'' = 0.19) and on the ramp at (''x''/''c'' = 0.86). | ||
Two-dimensional PIV data were | |||
acquired in a plane, along the model centreline and normal to the surface, starting from right | acquired in a plane, along the model centreline and normal to the surface, starting from right | ||
upstream of the slot and ending well beyond the reattachment location at x/c ≈ 1.4. Stereoscopic | upstream of the slot and ending well beyond the reattachment location at ''x''/''c'' ≈ 1.4. | ||
PIV (3D) data were acquired in planes perpendicular to the flow direction, arranged to intersect | Stereoscopic PIV (3D) data were acquired in planes perpendicular to the flow direction, arranged to intersect | ||
the 2D plane from x/c = 0.7 to 1.3 in steps of approximately 0.1. Oil-film interferometry was | the 2D plane from ''x''/''c'' = 0.7 to 1.3 in steps of approximately 0.1. | ||
Oil-film interferometry was | |||
used to quantify the skin friction over the entire model, from the region upstream of the hump to | used to quantify the skin friction over the entire model, from the region upstream of the hump to | ||
beyond the reattachment location. | beyond the reattachment location. | ||
Two-dimensionality of the flow in the separated and reattachment | |||
Two-dimensionality of the flow in the separated and reattachment regions was thoroughly | |||
assessed via three methods: by considering the spanwise pressures on the ramp in the separated | assessed via three methods: by considering the spanwise pressures on the ramp in the separated | ||
region, performing 3D PIV measurements in planes perpendicular to the flow direction, and by | region, performing 3D PIV measurements in planes perpendicular to the flow direction, and by | ||
Line 85: | Line 91: | ||
was also measured on the fore body of the model. It was shown that the spanwise variations of | was also measured on the fore body of the model. It was shown that the spanwise variations of | ||
the flow parameters are small (e.g., at the test condition, the pressure variation over the central | the flow parameters are small (e.g., at the test condition, the pressure variation over the central | ||
half of the model (–0.25 ≤ z/c ≤ 0.25) is | half of the model (–0.25 ≤ ''z''/''c'' ≤ 0.25) is ∆C<sub>P</sub> = ±0.005) | ||
At the inflow location (x/c = –2.14), pitot-probe and hot-wire anemometer data were compared | and that departures from twodimensionality are observed mainly near the wall. | ||
At the inflow location (''x''/''c'' = –2.14), pitot-probe and hot-wire anemometer data were compared | |||
with 2D and 3D PIV. Inflow skin friction was also documented using oil-film interferometry. | with 2D and 3D PIV. Inflow skin friction was also documented using oil-film interferometry. | ||
Resulting inflow velocity profile which was used in the numerical studies as a benchmark for | Resulting inflow velocity profile which was used in the numerical studies as a benchmark for | ||
imposing boundary conditions at the inflow of the computational domain (see the next section) is | imposing boundary conditions at the inflow of the computational domain (see the next section) is | ||
shown in Fig.4. | shown in [[UFR_3-34_Test_Case#figure4|Fig. 4]]. | ||
<div id="figure4"></div> | <div id="figure4"></div> | ||
Line 102: | Line 110: | ||
Experimental uncertainties reported in the original experimental papers are | Experimental uncertainties reported in the original experimental papers are | ||
summarized in | summarized in | ||
Table 2. Note that there is also an operational uncertainty associated with the endplates blockage | [[UFR_3-34_Test_Case#table2|Table 2]]. | ||
effect, which should be compensated in the quasi-2D simulations (see next sub-section). | Note that there is also an operational uncertainty associated with the endplates blockage | ||
effect, which should be compensated in the quasi-2D simulations | |||
(see [[UFR_3-34_Test_Case#CFD_Codes_and_Methods|next sub-section]]). | |||
Line 116: | Line 126: | ||
|Velocity||Less than 3% | |Velocity||Less than 3% | ||
|- | |- | ||
|Reynolds shear stress|Up to 20% | |Reynolds shear stress||Up to 20% | ||
|- | |- | ||
|Root mean square of streamwise and wall-normal velocity fluctuations (normal Reynolds stresses) | |Root mean square of streamwise and wall-normal velocity fluctuations (normal Reynolds stresses) | ||
Line 125: | Line 135: | ||
As mentioned in the Introduction, the present document focuses on the three groups of CFD | As mentioned in the Introduction, the present document focuses on the three groups of CFD | ||
studies of the considered test case, namely, on its predictions obtained with the use of different | studies of the considered test case, namely, on its predictions obtained with the use of different | ||
RANS models [7], hybrid RANS-LES models [8], [9] and WRLES [6] (see Table 3). | RANS models [‌[[UFR_3-34_References#7|7]]], | ||
hybrid RANS-LES models [‌[[UFR_3-34_References#8|8]]], [‌[[UFR_3-34_References#9|9]]] | |||
and WRLES [‌[[UFR_3-34_References#6|6]]] | |||
(see [[UFR_3-34_Test_Case#table3|Table 3]]). | |||
<div id="table3"></div> | |||
{|align="center" border=1 cellpadding=5px width=650px | |||
|+ style="caption-side:top" align="center"|'''Table 3:''' Summary of considered CFD studies | |||
|- | |||
!Source!!Turbulence Modelling approach | |||
|- | |||
|NASA Turbulence Modeling Resource Portal [‌[[UFR_3-34_References#7|7]]]||RANS | |||
|- | |||
|Go4Hybrid EU Project [‌[[UFR_3-34_References#8|8]]], [‌[[UFR_3-34_References#9|9]]] | |||
|Non-zonal and zonal RANS-LES hybrids | |||
|- | |||
|A. Uzun & M.R. Malik [‌[[UFR_3-34_References#6|6]]]||Wall Resolved LES | |||
|} | |||
Naturally, corresponding numerical procedures and codes are quite different. | Naturally, corresponding numerical procedures and codes are quite different. | ||
Line 135: | Line 164: | ||
and unstructured (FUN3D code https://fun3d.larc.nasa.gov/) grids. The grids used in the | and unstructured (FUN3D code https://fun3d.larc.nasa.gov/) grids. The grids used in the | ||
computations are built according to the well-know guidelines for solving RANS equations (the | computations are built according to the well-know guidelines for solving RANS equations (the | ||
grids are available on the portal [7]) and both codes return virtually identical results, thus | grids are available on the portal [‌[[UFR_3-34_References#7|7]]]) | ||
and both codes return virtually identical results, thus | |||
suggesting code- and (indirectly) grid-independence of the obtained solutions. | suggesting code- and (indirectly) grid-independence of the obtained solutions. | ||
The RANS-LES computations carried out in [8], [9] employed two enhanced hybrid models | The RANS-LES computations carried out | ||
in [‌[[UFR_3-34_References#8|8]]], [‌[[UFR_3-34_References#9|9]]] | |||
employed two enhanced hybrid models | |||
equipped with special tools for the “Grey-Area” Mitigation (GAM): a non-zonal method Delayed | equipped with special tools for the “Grey-Area” Mitigation (GAM): a non-zonal method Delayed | ||
DES (DDES) [10] with Shear-Layer Adapted (SLA) subgrid length scale [11] and a zonal | DES (DDES) [‌[[UFR_3-34_References#10|10]]] with Shear-Layer Adapted (SLA) subgrid length scale | ||
Improved DDES (IDDES) [12] combined with Synthetic Turbulence Generator (STG) [13] or | Δ<sub>SLA</sub> [‌[[UFR_3-34_References#11|11]]] | ||
Synthetic Eddy Method (SEM) [14], [15] for imposing boundary conditions at the RANS-IDDES | and a zonal | ||
interface. The shear-layer adapted subgrid length scale [11] replaces the length scale of the | Improved DDES (IDDES) [‌[[UFR_3-34_References#12|12]]] | ||
original DDES equal to the maximum grid spacing . It accounts for | combined with Synthetic Turbulence Generator (STG) [‌[[UFR_3-34_References#13|13]]] or | ||
Synthetic Eddy Method (SEM) [‌[[UFR_3-34_References#14|14]]], [‌[[UFR_3-34_References#15|15]]] | |||
for imposing unsteady boundary conditions at the RANS-IDDES | |||
interface. The shear-layer adapted subgrid length scale [‌[[UFR_3-34_References#11|11]]] | |||
replaces the length scale of the | |||
original DDES equal to the maximum grid spacing Δ<sub>max</sub> = max{Δ<sub>x</sub>, Δ<sub>y</sub>, Δ<sub>z</sub>}. | |||
It accounts for | |||
strong anisotropy typical of the grids in the initial part of separated shear layers and the quasi-2D | strong anisotropy typical of the grids in the initial part of separated shear layers and the quasi-2D | ||
character of the flow in this region. This substitution results in a considerable reduction of the | character of the flow in this region. This substitution results in a considerable reduction of the | ||
Line 150: | Line 188: | ||
A matrix of the performed simulations and appropriate references to original publications which | A matrix of the performed simulations and appropriate references to original publications which | ||
provide a detailed outline of the approaches used are presented in Table 4. | provide a detailed outline of the approaches used are presented in [[UFR_3-34_Test_Case#table4|Table 4]]. | ||
<div id="table4"></div> | |||
{|align="center" border=1 cellpadding=5px width=650px | |||
|+ style="caption-side:top" align="center"|'''Table 4:''' Matrix of hybrid RANS-LES computations [‌[[UFR_3-34_References#6|6]]], [‌[[UFR_3-34_References#8|8]]], [‌[[UFR_3-34_References#9|9]]] | |||
|- | |||
!Hybrid Method!!GAM tool!!Location of RANS-LES interface!!CFD Code | |||
|- | |||
|Zonal SST-based IDDES [12]||SEM [14], [15]||x/c = 0.5<br/>x/c = -1.0||Compressible TAU code with hybrid LD2 scheme [16] | |||
|- | |||
|Zonal SST-based IDDES [12]||STG [13]||x/c = 0.5<br/> x/c = -1.0|| Incompressible branch of NTS code [17] | |||
|- | |||
|Non-zonal SST-based DDES [10]||Shear-Layer Adapted length scale, Δ<sub>SLA</sub> [11] | |||
|No interface||Incompressible branch of NTS code [17] | |||
|} | |||
The non-zonal simulation was conducted with the use of the well-established in-house finite- | The non-zonal simulation was conducted with the use of the well-established in-house finite- | ||
Line 156: | Line 210: | ||
unstructured compressible finite-volume code TAU of DLR was employed. | unstructured compressible finite-volume code TAU of DLR was employed. | ||
The incompressible branch of the NTS code used in the simulations is based on the flux- | The incompressible branch of the NTS code used in the simulations is based on the flux-difference | ||
splitting method of Rogers and Kwak [‌[[UFR_3-34_References#18|18]]]. | |||
The approximation of the inviscid fluxes | |||
in the code depends on the turbulence representation approach: in the zonal RANS-IDDES | in the code depends on the turbulence representation approach: in the zonal RANS-IDDES | ||
computations, it is a | computations, it is a 3<sup>rd</sup>-order upwind-biased scheme in the RANS zone and a 4<sup>th</sup>-order | ||
scheme in the WMLES zone, whereas for the non-zonal DDES the hybrid (weighted | central scheme in the WMLES zone, whereas for the non-zonal DDES the hybrid (weighted 3<sup>rd</sup>-order | ||
upwind-biased / | upwind-biased / 4<sup>th</sup>-order central) scheme [‌[[UFR_3-34_References#19|19]]] is used. | ||
the | The viscous fluxes are approximated with | ||
scheme with sub-iterations is applied. | the 2<sup>nd</sup>-order central scheme. | ||
For the time integration, an implicit 2<sup>nd</sup>-order backward Euler scheme with sub-iterations is applied. | |||
The TAU code employs a low-dissipation low-dispersion (LD2) scheme [16] which is based on | The TAU code employs a low-dissipation low-dispersion (LD2) scheme [‌[[UFR_3-34_References#16|16]]] | ||
the | which is based on the 2<sup>nd</sup>-order energy-conserving skew-symmetric convection operator | ||
level of | combined with a minimal level of 4<sup>th</sup>-order artificial matrix dissipation for stabilization. | ||
The central flux terms employ an | |||
additional gradient extrapolation that effectively increases the discretization stencil and is used to | additional gradient extrapolation that effectively increases the discretization stencil and is used to | ||
reduce the dispersion error of the scheme. Both ingredients are essential in terms of turbulence | reduce the dispersion error of the scheme. Both ingredients are essential in terms of turbulence | ||
Line 178: | Line 235: | ||
grid (“mandatory” grid in the Go4Hybrid project). | grid (“mandatory” grid in the Go4Hybrid project). | ||
The computational domain in the XY-plane is shown in Fig. 5. Its size in the spanwise direction | The computational domain in the ''XY''-plane is shown in [[UFR_3-34_Test_Case#figure5|Fig. 5]]. | ||
is equal to 0. | Its size in the spanwise direction | ||
the ATAAC project [4]). As recommended in [7], the contour of the upper (slip) wall of the | is equal to 0.4''c'' (this has been proven sufficient to arrive at a span-independent solution within | ||
the ATAAC project [‌[[UFR_3-34_References#4|4]]]). | |||
As recommended in [‌[[UFR_3-34_References#7|7]]], the contour of the upper (slip) wall of the | |||
domain was modified (moved somewhat downwards in the area above the hump) in order to | domain was modified (moved somewhat downwards in the area above the hump) in order to | ||
compensate the blockage effect of the endplates (see Fig. 1 above). | compensate the blockage effect of the endplates (see [[UFR_3-34_Test_Case#figure1|Fig. 1]] above). | ||
Line 195: | Line 254: | ||
The inflow and outflow conditions in the simulations performed with the use of the TAU and | The inflow and outflow conditions in the simulations performed with the use of the TAU and | ||
NTS codes were somewhat different since the former code is compressible and in the latter one | NTS codes were somewhat different since the former code is compressible and in the latter one | ||
its incompressible branch was employed. In both cases velocity and turbulence quantities (k and | its incompressible branch was employed. In both cases velocity and turbulence quantities (''k'' and | ||
''ω'') at the inflow boundary (''x''/''c'' = -2.14) were specified from a precursor SST RANS of the zero | |||
pressure gradient boundary layer carried out to match the experimental value of the momentum | pressure gradient boundary layer carried out to match the experimental value of the momentum | ||
thickness based Re-number (see Fig. 4 above). In the | thickness based Re-number | ||
<math>{\text{Re}_{\theta_\text{inflow}}=\theta_\text{inflow} U_\infty/\nu=7200}</math> | |||
(see [[UFR_3-34_Test_Case#figure4|Fig. 4]] above). In the | |||
incompressible simulations (NTS code), the inflow static pressure was defined by the linear | incompressible simulations (NTS code), the inflow static pressure was defined by the linear | ||
extrapolation of the values in the nearby interior cells of the domain, whereas in the compressible | extrapolation of the values in the nearby interior cells of the domain, whereas in the compressible | ||
Line 213: | Line 274: | ||
|[[Image:UFR3-34_Fig6.png|650px]] | |[[Image:UFR3-34_Fig6.png|650px]] | ||
|- | |- | ||
|align="center"|'''Figure 6:''' Entire grid in XY-plane used in hybrid RANS-LES simulations and its zoomed-in view [8],[9] | |align="center"|'''Figure 6:''' Entire grid in ''XY''-plane used in hybrid RANS-LES simulations and its zoomed-in view [‌[[UFR_3-34_References#8|8]]], [‌[[UFR_3-34_References#9|9]]] | ||
|} | |} | ||
The grid in the XY-plane used in all the simulations is shown in Fig. 6. It has 511×127 cells in | The grid in the ''XY''-plane used in all the simulations is shown in [[UFR_3-34_Test_Case#figure6|Fig. 6]]. | ||
the x- and y-directions, respectively. Near the bottom wall, the grid is wall-orthogonal and yields | It has 511×127 cells in | ||
a first cell height in wall units () less than 1.0 in the whole domain. For -1.0 < x/c < 0.65, the | the ''x''- and ''y''-directions, respectively. Near the bottom wall, the grid is wall-orthogonal and yields | ||
grid is also refined in the x-direction in order to satisfy demands to WMLES grids in the outer | a first cell height in wall units (Δy<sub>1</sub><sup>+</sup>) less than 1.0 in the whole domain. | ||
flow units (), which allows placing the RANS-LES interface at any cross-section | For -1.0 < ''x''/''c'' < 0.65, the | ||
within this region. In the separation (“focus”) region (0.67 < x/c < 1.5) both x- and maximum y- | grid is also refined in the ''x''-direction in order to satisfy demands to WMLES grids in the outer | ||
steps are equal to 0. | flow units (Δ''x'' / ''δ'' < 0.1), which allows placing the RANS-LES interface at any cross-section | ||
step size equal to 0. | within this region. In the separation (“focus”) region (0.67 < ''x'' / ''c'' < 1.5) | ||
both ''x''- and maximum ''y''-steps are equal to 0.005''c''. | |||
In the spanwise direction, the grid has 80 cells and is uniform with a | |||
step size equal to 0.005''c'' (Δ''z'' / ''δ'' < 0.05), | |||
thus resulting in a nearly cubic grid in the focus region. | |||
As a result, the total grid has around 5.2 million cells, that is, rather modest by today's standards. | As a result, the total grid has around 5.2 million cells, that is, rather modest by today's standards. | ||
The physical time step is set to | The physical time step is set to Δ''t'' = 0.002''c'' / ''U''<sub>∞</sub> | ||
which results in a convective CFL number of | |||
well below 0.5 in the focus region. The total simulation time was around 30 Convective Time | well below 0.5 in the focus region. The total simulation time was around 30 Convective Time | ||
Units, CTU (1 CTU = c/ | Units, CTU (1 CTU = ''c'' / ''U''<sub>∞</sub>), of which the last 20 CTUs were used for obtaining time-averaged | ||
flow data and turbulence statistics. | flow data and turbulence statistics. | ||
The last numerical database for the 2DWMH flow discussed in the present document is that | The last numerical database for the 2DWMH flow discussed in the present document is that | ||
accumulated in the course of the study [6] where compressible Wall Resolved LES of this flow | accumulated in the course of the study [‌[[UFR_3-34_References#6|6]]] | ||
has been performed (see Table 3). The code used in this study has been previously successfully | where compressible Wall Resolved LES of this flow | ||
applied to LES of different wall-bounded turbulent flows (see, e.g., [20] | has been performed (see [[UFR_3-34_Test_Case#table3|Table 3]]). | ||
The code used in this study has been previously successfully | |||
applied to LES of different wall-bounded turbulent flows | |||
(see, e.g., [‌[[UFR_3-34_References#20|20]]] – [‌[[UFR_3-34_References#22|22]]]. | |||
It solves the | |||
unsteady fully three-dimensional compressible Navier-Stokes equations discretized on multi- | unsteady fully three-dimensional compressible Navier-Stokes equations discretized on multi- | ||
block structured overset grids either with no SGS model (Implicit LES or ILES) or with the static | block structured overset grids either with no SGS model (Implicit LES or ILES) or with the static | ||
SGS model of Vreman [23]. | SGS model of Vreman [‌[[UFR_3-34_References#23|23]]]. | ||
All the spatial derivatives in the code were computed employing an optimized pre-factored | All the spatial derivatives in the code were computed employing an optimized pre-factored | ||
fourth-order compact finite-difference scheme [24], which offers improved dispersion | fourth-order compact finite-difference scheme [‌[[UFR_3-34_References#24|24]]], | ||
which offers improved dispersion | |||
characteristics. To eliminate spurious high-frequency oscillations that may arise from several | characteristics. To eliminate spurious high-frequency oscillations that may arise from several | ||
sources (grid stretching, unresolved fluctuations, approximation of physical boundary conditions) | sources (grid stretching, unresolved fluctuations, approximation of physical boundary conditions) | ||
and to ensure numerical stability, a sixth-order compact filtering scheme [25], [26] was | and to ensure numerical stability, a sixth-order compact filtering | ||
scheme [‌[[UFR_3-34_References#25|25]]], [‌[[UFR_3-34_References#26|26]]] was | |||
employed. Finally, in order to maintain high-order accuracy throughout the entire computational | employed. Finally, in order to maintain high-order accuracy throughout the entire computational | ||
domain, a sixth-order accurate explicit Lagrangian interpolation [27] was applied whenever | domain, a sixth-order accurate explicit Lagrangian | ||
interpolation [‌[[UFR_3-34_References#27|27]]] was applied whenever | |||
overset grids were used. A Beam-Warming type approximately factorized implicit scheme was | overset grids were used. A Beam-Warming type approximately factorized implicit scheme was | ||
used for the time advancement (more details of the simulation methodology can be found in [20], | used for the time advancement (more details of the simulation methodology can be found | ||
[21], [28], [29]). | in [‌[[UFR_3-34_References#20|20]]], [‌[[UFR_3-34_References#21|21]]], [‌[[UFR_3-34_References#28|28]]], [‌[[UFR_3-34_References#29|29]]]). | ||
The simulations were performed in the same domain in XY-plane as used in the hybrid RANS- | The simulations were performed in the same domain in ''XY''-plane as used in the | ||
LES computations (see Fig. 5) with two different widths, one being the same as that in the hybrid | hybrid RANS-LES computations (see [[UFR_3-34_Test_Case#figure5|Fig. 5]]) | ||
simulations and the other two times less, i.e., 0. | with two different widths, one being the same as that in the hybrid | ||
spanwise, and wall-normal directions in the wall units | simulations and the other two times less, i.e., 0.2''c''. The grid-spacings in the streamwise, | ||
the boundary layer (up to about y+ | spanwise, and wall-normal directions in the wall units Δ''x''<sup>+</sup> , | ||
Δ''z''<sup>+</sup>, and Δ''y''<sup>+</sup>in the near wall part of | |||
the boundary layer (up to about ''y''<sup>+</sup> ≈ 200) | |||
were equal to about 25, 12.5 and 0.8, respectively, and | |||
then the grid was coarsened by a factor of two in both the streamwise and spanwise directions. | then the grid was coarsened by a factor of two in both the streamwise and spanwise directions. | ||
This was possible thanks to the solver overset-grid capability (a fragment of the two-level | This was possible thanks to the solver overset-grid capability (a fragment of the two-level | ||
overset-grid in YZ-plane is shown in Fig.7). As a result, the total grid had about 825 million | overset-grid in ''YZ''-plane is shown in [[UFR_3-34_Test_Case#figure7|Fig. 7]]). | ||
As a result, the total grid had about 825 million | |||
points, which ensured resolving the boundary layer approaching the hump by 90 to 100 points | points, which ensured resolving the boundary layer approaching the hump by 90 to 100 points | ||
and the thinnest part of the boundary layer over the hump by about 30 points in the wall-normal | and the thinnest part of the boundary layer over the hump by about 30 points in the wall-normal | ||
Line 270: | Line 347: | ||
|} | |} | ||
The time step of all simulations was 2.5×10- | The time step of all simulations was 2.5×10<sup>-4</sup>''c''/''a''<sub>ref</sub> | ||
2.5×10- | (''a''<sub>ref</sub> is the reference speed of sound) or | ||
2.5×10<sup>-3</sup>''c''/''U''<sub>ref</sub>. Statistical results were averaged over about 10 CTU. | |||
Similar to the hybrid simulations outlined above, the top boundary of the computational domain | Similar to the hybrid simulations outlined above, the top boundary of the computational domain | ||
was treated as an inviscid wall while viscous adiabatic boundary conditions where imposed on | was treated as an inviscid wall while viscous adiabatic boundary conditions where imposed on | ||
the lower (hump) wall. In order to eliminate reflections from the outlet boundary, a sponge zone | the lower (hump) wall. In order to eliminate reflections from the outlet boundary, a sponge zone | ||
was placed downstream of the physical region of interest (x/c > 1.6). At the outflow boundary of | was placed downstream of the physical region of interest (''x''/''c'' > 1.6). | ||
the domain (x/c = 4) the back pressure was specified (set slightly below the upstream pressure) | At the outflow boundary of | ||
the domain (''x''/''c'' = 4) the back pressure was specified (set slightly below the upstream pressure) | |||
and characteristic boundary conditions were applied, whereas the inflow conditions were based | and characteristic boundary conditions were applied, whereas the inflow conditions were based | ||
on the characteristic relaxation boundary conditions [30] that inject turbulent fluctuations | on the characteristic relaxation boundary conditions [‌[[UFR_3-34_References#30|30]]] | ||
that inject turbulent fluctuations | |||
(generated by the inflow generation technique) in the boundary layer while allowing upstream- | (generated by the inflow generation technique) in the boundary layer while allowing upstream- | ||
travelling waves to exit the domain. | travelling waves to exit the domain. |
Latest revision as of 11:09, 15 March 2018
Semi-Confined Flows
Underlying Flow Regime 3-34
Test Case
Brief Description of the Study Test Case
A 3D sketch of the experimental setup is shown in Fig. 1. The configuration presents a Glauert-Goldschmied type body consisting of a relatively long fore body and a relatively short concave ramp comprising the aft part of the model mounted between two glass endplate frames with both leading edge and trailing edges faired smoothly with a wind tunnel splitter plate.
Figure 1: 3D sketch of experimental setup [1], [2] |
Major geometric and flow parameters of the TC are presented in Fig. 2 and summarized in Table 1 (note that in the “baseline” experiment considered here the slot shown in Fig. 2 was closed).
Figure 2: Schematic of TC geometry and flow parameters |
Parameter | Notation | Value |
---|---|---|
Free stream velocity | 34.6 m/s | |
Hump chord | 0.42 m | |
Crest height | 0.0537 m | |
Reynolds number | 936,000 | |
Mach number | 0.1 |
The experimental data are available at http://cfdval2004.larc.nasa.gov/case3expdata.html
(Case 3), at https://turbmodels.larc.nasa.gov/nasahump_val.html, and also enter the ERCOFTAC
Classic Collection http://cfd.mace.manchester.ac.uk/ercoftac under number C.83. The data set
includes: streamwise distributions of the surface pressure and skin-friction coefficients,
and ,
and mean velocity and Reynolds stresses fields in the tunnel center-plane, roughly covering
the region 0.63 < x/c < 1.39.
Test Case Experiments
A photo of the experimental rig is presented in Fig. 3.
Figure 3: Photo of the experimental rig |
The experiments were performed in the NASA Langley 20′′×28′′ shear flow tunnel. The flow was nominally 2D, although with side-wall effects (3D flow) expected near the endplates. The reference “chord” length of the model, c, is defined as the length of the hump on the wall and is equal to 420mm. The maximum thickness of the hump, h, is equal to 53.7mm. The model was equipped with 153 centre-span static pressure ports and 20 dynamic pressure ports in the vicinity of the separated flow region. Sixteen spanwise pressure ports were located on the fore body (x/c = 0.19) and on the ramp at (x/c = 0.86). Two-dimensional PIV data were acquired in a plane, along the model centreline and normal to the surface, starting from right upstream of the slot and ending well beyond the reattachment location at x/c ≈ 1.4. Stereoscopic PIV (3D) data were acquired in planes perpendicular to the flow direction, arranged to intersect the 2D plane from x/c = 0.7 to 1.3 in steps of approximately 0.1. Oil-film interferometry was used to quantify the skin friction over the entire model, from the region upstream of the hump to beyond the reattachment location.
Two-dimensionality of the flow in the separated and reattachment regions was thoroughly assessed via three methods: by considering the spanwise pressures on the ramp in the separated region, performing 3D PIV measurements in planes perpendicular to the flow direction, and by means of the surface oil-film flow visualization. Spanwise distribution of the surface pressure was also measured on the fore body of the model. It was shown that the spanwise variations of the flow parameters are small (e.g., at the test condition, the pressure variation over the central half of the model (–0.25 ≤ z/c ≤ 0.25) is ∆CP = ±0.005) and that departures from twodimensionality are observed mainly near the wall.
At the inflow location (x/c = –2.14), pitot-probe and hot-wire anemometer data were compared with 2D and 3D PIV. Inflow skin friction was also documented using oil-film interferometry. Resulting inflow velocity profile which was used in the numerical studies as a benchmark for imposing boundary conditions at the inflow of the computational domain (see the next section) is shown in Fig. 4.
Figure 4: Experimental profile of streamwise velocity at x/c = -2.14 (momentum thickness Reθ=7200) |
Experimental uncertainties reported in the original experimental papers are
summarized in
Table 2.
Note that there is also an operational uncertainty associated with the endplates blockage
effect, which should be compensated in the quasi-2D simulations
(see next sub-section).
Quantity | Experimental Uncertainty |
---|---|
Pressure coefficient, | ± 0.001 |
Velocity | Less than 3% |
Reynolds shear stress | Up to 20% |
Root mean square of streamwise and wall-normal velocity fluctuations (normal Reynolds stresses) | Up to 14% |
CFD Codes and Methods
As mentioned in the Introduction, the present document focuses on the three groups of CFD studies of the considered test case, namely, on its predictions obtained with the use of different RANS models [7], hybrid RANS-LES models [8], [9] and WRLES [6] (see Table 3).
Source | Turbulence Modelling approach |
---|---|
NASA Turbulence Modeling Resource Portal [7] | RANS |
Go4Hybrid EU Project [8], [9] | Non-zonal and zonal RANS-LES hybrids |
A. Uzun & M.R. Malik [6] | Wall Resolved LES |
Naturally, corresponding numerical procedures and codes are quite different.
In particular, the RANS based computations, results of which are presented on the NASA Turbulence Modeling Resource portal were carried out with the use of two long-standing compressible Navier-Stokes CFD codes developed at NASA Langley Research Center for solving 2D and 3D compressible flows on structured (CFL3D code https://cfl3d.larc.nasa.gov/) and unstructured (FUN3D code https://fun3d.larc.nasa.gov/) grids. The grids used in the computations are built according to the well-know guidelines for solving RANS equations (the grids are available on the portal [7]) and both codes return virtually identical results, thus suggesting code- and (indirectly) grid-independence of the obtained solutions.
The RANS-LES computations carried out in [8], [9] employed two enhanced hybrid models equipped with special tools for the “Grey-Area” Mitigation (GAM): a non-zonal method Delayed DES (DDES) [10] with Shear-Layer Adapted (SLA) subgrid length scale ΔSLA [11] and a zonal Improved DDES (IDDES) [12] combined with Synthetic Turbulence Generator (STG) [13] or Synthetic Eddy Method (SEM) [14], [15] for imposing unsteady boundary conditions at the RANS-IDDES interface. The shear-layer adapted subgrid length scale [11] replaces the length scale of the original DDES equal to the maximum grid spacing Δmax = max{Δx, Δy, Δz}. It accounts for strong anisotropy typical of the grids in the initial part of separated shear layers and the quasi-2D character of the flow in this region. This substitution results in a considerable reduction of the subgrid eddy-viscosity, thus ensuring a rapid development of 3D turbulence structures.
A matrix of the performed simulations and appropriate references to original publications which provide a detailed outline of the approaches used are presented in Table 4.
Hybrid Method | GAM tool | Location of RANS-LES interface | CFD Code |
---|---|---|---|
Zonal SST-based IDDES [12] | SEM [14], [15] | x/c = 0.5 x/c = -1.0 |
Compressible TAU code with hybrid LD2 scheme [16] |
Zonal SST-based IDDES [12] | STG [13] | x/c = 0.5 x/c = -1.0 |
Incompressible branch of NTS code [17] |
Non-zonal SST-based DDES [10] | Shear-Layer Adapted length scale, ΔSLA [11] | No interface | Incompressible branch of NTS code [17] |
The non-zonal simulation was conducted with the use of the well-established in-house finite-
volume structured multi-block code NTS, whereas for the zonal ones, along with this code, an
unstructured compressible finite-volume code TAU of DLR was employed.
The incompressible branch of the NTS code used in the simulations is based on the flux-difference splitting method of Rogers and Kwak [18]. The approximation of the inviscid fluxes in the code depends on the turbulence representation approach: in the zonal RANS-IDDES computations, it is a 3rd-order upwind-biased scheme in the RANS zone and a 4th-order central scheme in the WMLES zone, whereas for the non-zonal DDES the hybrid (weighted 3rd-order upwind-biased / 4th-order central) scheme [19] is used. The viscous fluxes are approximated with the 2nd-order central scheme. For the time integration, an implicit 2nd-order backward Euler scheme with sub-iterations is applied.
The TAU code employs a low-dissipation low-dispersion (LD2) scheme [16] which is based on the 2nd-order energy-conserving skew-symmetric convection operator combined with a minimal level of 4th-order artificial matrix dissipation for stabilization. The central flux terms employ an additional gradient extrapolation that effectively increases the discretization stencil and is used to reduce the dispersion error of the scheme. Both ingredients are essential in terms of turbulence resolution with the unstructured TAU code. Note that in the present zonal RANS-WMLES computations, the LD2 scheme is only active in the respective IDDES region downstream of the interface. The temporal discretization is based on an implicit dual-time stepping scheme which is also of 2nd-order accuracy.
All the hybrid simulations were performed in the same computational domain and on the same grid (“mandatory” grid in the Go4Hybrid project).
The computational domain in the XY-plane is shown in Fig. 5. Its size in the spanwise direction is equal to 0.4c (this has been proven sufficient to arrive at a span-independent solution within the ATAAC project [4]). As recommended in [7], the contour of the upper (slip) wall of the domain was modified (moved somewhat downwards in the area above the hump) in order to compensate the blockage effect of the endplates (see Fig. 1 above).
Figure 5: Schematic of computational domain in XY-plane and wall boundary conditions [8], [9] |
The inflow and outflow conditions in the simulations performed with the use of the TAU and
NTS codes were somewhat different since the former code is compressible and in the latter one
its incompressible branch was employed. In both cases velocity and turbulence quantities (k and
ω) at the inflow boundary (x/c = -2.14) were specified from a precursor SST RANS of the zero
pressure gradient boundary layer carried out to match the experimental value of the momentum
thickness based Re-number
(see Fig. 4 above). In the
incompressible simulations (NTS code), the inflow static pressure was defined by the linear
extrapolation of the values in the nearby interior cells of the domain, whereas in the compressible
simulations (TAU code) flux boundary conditions were used together with an extrapolation of
the pressure from the interior. At the outflow boundary, in the NTS code a constant static
pressure is specified and all the other flow variables are defined by linear extrapolation from the
interior cells of the domain, which is similar to TAU’s flux conditions with an imposed fixed
exit pressure. Finally, in all simulations the boundary conditions in the spanwise direction were
periodic.
Figure 6: Entire grid in XY-plane used in hybrid RANS-LES simulations and its zoomed-in view [8], [9] |
The grid in the XY-plane used in all the simulations is shown in Fig. 6.
It has 511×127 cells in
the x- and y-directions, respectively. Near the bottom wall, the grid is wall-orthogonal and yields
a first cell height in wall units (Δy1+) less than 1.0 in the whole domain.
For -1.0 < x/c < 0.65, the
grid is also refined in the x-direction in order to satisfy demands to WMLES grids in the outer
flow units (Δx / δ < 0.1), which allows placing the RANS-LES interface at any cross-section
within this region. In the separation (“focus”) region (0.67 < x / c < 1.5)
both x- and maximum y-steps are equal to 0.005c.
In the spanwise direction, the grid has 80 cells and is uniform with a
step size equal to 0.005c (Δz / δ < 0.05),
thus resulting in a nearly cubic grid in the focus region.
As a result, the total grid has around 5.2 million cells, that is, rather modest by today's standards.
The physical time step is set to Δt = 0.002c / U∞
which results in a convective CFL number of
well below 0.5 in the focus region. The total simulation time was around 30 Convective Time
Units, CTU (1 CTU = c / U∞), of which the last 20 CTUs were used for obtaining time-averaged
flow data and turbulence statistics.
The last numerical database for the 2DWMH flow discussed in the present document is that accumulated in the course of the study [6] where compressible Wall Resolved LES of this flow has been performed (see Table 3). The code used in this study has been previously successfully applied to LES of different wall-bounded turbulent flows (see, e.g., [20] – [22]. It solves the unsteady fully three-dimensional compressible Navier-Stokes equations discretized on multi- block structured overset grids either with no SGS model (Implicit LES or ILES) or with the static SGS model of Vreman [23].
All the spatial derivatives in the code were computed employing an optimized pre-factored fourth-order compact finite-difference scheme [24], which offers improved dispersion characteristics. To eliminate spurious high-frequency oscillations that may arise from several sources (grid stretching, unresolved fluctuations, approximation of physical boundary conditions) and to ensure numerical stability, a sixth-order compact filtering scheme [25], [26] was employed. Finally, in order to maintain high-order accuracy throughout the entire computational domain, a sixth-order accurate explicit Lagrangian interpolation [27] was applied whenever overset grids were used. A Beam-Warming type approximately factorized implicit scheme was used for the time advancement (more details of the simulation methodology can be found in [20], [21], [28], [29]).
The simulations were performed in the same domain in XY-plane as used in the hybrid RANS-LES computations (see Fig. 5) with two different widths, one being the same as that in the hybrid simulations and the other two times less, i.e., 0.2c. The grid-spacings in the streamwise, spanwise, and wall-normal directions in the wall units Δx+ , Δz+, and Δy+in the near wall part of the boundary layer (up to about y+ ≈ 200) were equal to about 25, 12.5 and 0.8, respectively, and then the grid was coarsened by a factor of two in both the streamwise and spanwise directions. This was possible thanks to the solver overset-grid capability (a fragment of the two-level overset-grid in YZ-plane is shown in Fig. 7). As a result, the total grid had about 825 million points, which ensured resolving the boundary layer approaching the hump by 90 to 100 points and the thinnest part of the boundary layer over the hump by about 30 points in the wall-normal direction.
Figure 7: Two-level overset-grid system at x/c = -2 [6] |
The time step of all simulations was 2.5×10-4c/aref (aref is the reference speed of sound) or 2.5×10-3c/Uref. Statistical results were averaged over about 10 CTU.
Similar to the hybrid simulations outlined above, the top boundary of the computational domain was treated as an inviscid wall while viscous adiabatic boundary conditions where imposed on the lower (hump) wall. In order to eliminate reflections from the outlet boundary, a sponge zone was placed downstream of the physical region of interest (x/c > 1.6). At the outflow boundary of the domain (x/c = 4) the back pressure was specified (set slightly below the upstream pressure) and characteristic boundary conditions were applied, whereas the inflow conditions were based on the characteristic relaxation boundary conditions [30] that inject turbulent fluctuations (generated by the inflow generation technique) in the boundary layer while allowing upstream- travelling waves to exit the domain.
Note finally that in the framework of WRLES several sensitivity tests were carried out aimed at
finding out effects of the grid-resolution and span size of the computational domain and of the
operational uncertainty associated with the blockage effect of end plates.
Contributed by: E. Guseva, M. Strelets — Peter the Great St. Petersburg Polytechnic University (SPbPU)
© copyright ERCOFTAC 2024