UFR 3-31 Test Case: Difference between revisions

From KBwiki
Jump to navigation Jump to search
m (Dave.Ellacott moved page SilverP:UFR 3-31 Test Case to UFR 3-31 Test Case)
 
(179 intermediate revisions by 2 users not shown)
Line 1: Line 1:
= Flow over curved backward-facing step =
= Flow over curved backward-facing step =
{{UFRHeader
{{UFRHeader
Line 13: Line 12:
<!--{{LoremIpsum}}-->
<!--{{LoremIpsum}}-->
== Brief Description of the Study Test Case ==
== Brief Description of the Study Test Case ==
The geometry under consideration is shown in Fig. \ref{fig:geometry}. The  
The geometry under consideration is shown in [[UFR_3-31_Test_Case#figure2|Fig. 2]]. The  
rounded ramp of height <math>{H}</math> is placed in a high-aspect-ratio duct with upstream  
rounded ramp of height <math>{H}</math> is placed in a high-aspect-ratio duct with upstream  
height of <math>{8.52H}</math>. In the simulations, the flow is assumed to be spanwise  
height of <math>{8.52H}</math>. In the simulations, the flow is assumed to be spanwise  
homogeneous, with the spanwise slab being <math>{3.7H}</math>. The assumption of homogeneity  
homogeneous, with the spanwise slab being <math>{3.7H}</math>. The assumption of homogeneity  
is justified by the fact that the experimental ratio of duct depth to the step  
is justified by the fact that the experimental ratio of duct depth to the step  
height was 38. In the experiment \cite{zhang2010experimental}, tripped boundary layers were allowed to  
height was 38. In the experiment
[&#8204;[[UFR_3-31_References#31|31]]],
tripped boundary layers were allowed to  
develop on both walls for a distance of about <math>{25H}</math>. The Reynolds number, based  
develop on both walls for a distance of about <math>{25H}</math>. The Reynolds number, based  
on <math>H</math> and the inlet free-stream velocity <math>{U_{in}}</math>, is 13,700. At <math>{x/H = -7.4}</math>,  
on <math>H</math> and the inlet free-stream velocity <math>{U_{in}}</math>, is 13,700. At <math>{x/H = -7.4}</math>,  
Line 24: Line 25:
<math>{Re_\theta = 1192}</math>, and the boundary-layer thickness is <math>{\delta_{99} = 0.83H}</math>.
<math>{Re_\theta = 1192}</math>, and the boundary-layer thickness is <math>{\delta_{99} = 0.83H}</math>.


The step geometry is based on that used originally by Song and Eaton \cite{song2004reynolds}.  
<div id="figure2"></div>
{|align="center" border="0" width="750"
|[[Image:UFR3-31_fig2.png|740px]]
|-
!align="center"|Figure 2: Flow domain and grid (every 8th point shown, fine mesh).
|}
 
 
The step geometry is based on that used originally by Song and Eaton
[&#8204;[[UFR_3-31_References#27|27]]].
In order to enlarge the separated region, the height of the step was increased  
In order to enlarge the separated region, the height of the step was increased  
by a factor of 1.5. This adaptation was undertaken interactively with a parallel  
by a factor of 1.5. This adaptation was undertaken interactively with a parallel  
experiment by Zhang and Zhong \cite{zhang2010experimental}. The step shape is  
experiment by Zhang and Zhong
[&#8204;[[UFR_3-31_References#31|31]]].
The step shape is  
described by the following three relations, with the origin <math>{x/H = 0}</math> being the  
described by the following three relations, with the origin <math>{x/H = 0}</math> being the  
upstream edge of the ramp:
upstream edge of the ramp:
Line 33: Line 45:
<!--
<!--
{{DisplayEquation|eqn=y_{wall} = (1 - R_1) + \sqrt{R_1^2 - x_2} \qquad \mathrm{for} \qquad 0 < x/H \leq 2.3}}
{{DisplayEquation|eqn=y_{wall} = (1 - R_1) + \sqrt{R_1^2 - x_2} \qquad \mathrm{for} \qquad 0 < x/H \leq 2.3}}
{{DisplayEquation|eqn=y_{wall} = y_2 - \sqrt{\frac{1}{4}R_1^2 - ( x_2 - x)^2} \qquad \mathrm{for} \qquad 2.3 < x/H \leq 2.835}}
{{DisplayEquation|no=1|eqn=y_{wall} = y_2 - \sqrt{\frac{1}{4}R_1^2 - ( x_2 - x)^2} \qquad \mathrm{for} \qquad 2.3 < x/H \leq 2.835}}
{{DisplayEquation|eqn=y_{wall} = R_2 - \sqrt{R_2^2 - (3 - x)^2} \qquad \mathrm{for} \qquad 2.835 < x/H \leq 2.937}}
{{DisplayEquation|eqn=y_{wall} = R_2 - \sqrt{R_2^2 - (3 - x)^2} \qquad \mathrm{for} \qquad 2.835 < x/H \leq 2.937}}
-->
-->
<div id="eqn1"></div>
<center><math>{
<center><math>{
\begin{align}
\begin{align}
y_{wall} &= (1 - R_1) + \sqrt{R_1^2 - x_2} \qquad \mathrm{for} \qquad 0 < x/H \leq 2.3 \\
y_{wall} &= (1 - R_1) + \sqrt{R_1^2 - x_2} \qquad \mathrm{for} \qquad 0 < x/H \leq 2.3 \\
y_{wall} &= y_2 - \sqrt{\frac{1}{4}R_1^2 - ( x_2 - x)^2} \qquad \mathrm{for} \qquad 2.3 < x/H \leq 2.835 \\
y_{wall} &= y_2 - \sqrt{\frac{1}{4}R_1^2 - ( x_2 - x)^2} \qquad \mathrm{for} \qquad 2.3 < x/H \leq 2.835
\qquad\qquad (1) \\
y_{wall} &= R_2 - \sqrt{R_2^2 - (3 - x)^2} \qquad \mathrm{for} \qquad 2.835 < x/H \leq 2.937 \\
y_{wall} &= R_2 - \sqrt{R_2^2 - (3 - x)^2} \qquad \mathrm{for} \qquad 2.835 < x/H \leq 2.937 \\
\end{align}
\end{align}
Line 76: Line 90:
and the skin-friction is deduced from the LDA velocity measurements, assuming the first point  
and the skin-friction is deduced from the LDA velocity measurements, assuming the first point  
being in the viscous sub-layer.
being in the viscous sub-layer.
The experimental study is described in detail in [&#8204;[[UFR_3-31_References#31|31]]].


== CFD Methods ==
== CFD Methods ==
Line 82: Line 97:
The implicitly filtered LES momentum and continuity equations for incompressible  
The implicitly filtered LES momentum and continuity equations for incompressible  
flow were solved over a general non-orthogonal, boundary-fitted, multi-block  
flow were solved over a general non-orthogonal, boundary-fitted, multi-block  
finite volume mesh, using LES-STREAM (e.g. \cite{fishpool2009stability} for  
finite volume mesh, using LES-STREAM (e.g.
more details and related references). In total, the default mesh covering the duct-flow domain  
[&#8204;[[UFR_3-31_References#7|7]]]
for more details and related references). In total, the default mesh covering the duct-flow domain  
contained around 24 million nodes. The sensitivity to the resolution is addressed  
contained around 24 million nodes. The sensitivity to the resolution is addressed  
below. The variables are stored in a co-located manner. The solution is based on  
below. The variables are stored in a co-located manner. The solution is based on  
a fractional-step time-marching method, with the time derivative approximated by  
a fractional-step time-marching method, with the time derivative approximated by  
a third-order Gear scheme \cite{fishpool2009stability}. The fluxes are approximated  
a third-order Gear scheme
[&#8204;[[UFR_3-31_References#7|7]]].
The fluxes are approximated  
by second-order centred approximations. Within the fractional-step algorithm, a  
by second-order centred approximations. Within the fractional-step algorithm, a  
provisional velocity field results from advancing the solution with the flux  
provisional velocity field results from advancing the solution with the flux  
Line 95: Line 113:
of a V-cycle multigrid scheme. In order to suppress unphysical oscillations,  
of a V-cycle multigrid scheme. In order to suppress unphysical oscillations,  
associated with pressure-velocity decoupling, a practice equivalent to that  
associated with pressure-velocity decoupling, a practice equivalent to that  
introduced in \cite{rhie1983numerical} is adopted. Fishpool \& Leschziner   
introduced in
\cite{fishpool2009stability} demonstrate that the loss of accuracy associated  
[&#8204;[[UFR_3-31_References#23|23]]]
is adopted. Fishpool &amp; Leschziner   
[&#8204;[[UFR_3-31_References#7|7]]]
demonstrate that the loss of accuracy associated  
with the smoothing introduced by this practice is minimal.  
with the smoothing introduced by this practice is minimal.  
The subgrid-scale stresses were approximated using either the dynamic
The subgrid-scale stresses were approximated using either the dynamic
Smagorinsky model of \cite{germanoetal91} and \cite{lilly92} or the mixed-time
Smagorinsky model of
scale model of \cite{inagakietal05}, the latter being local, involving
[&#8204;[[UFR_3-31_References#10|10]]] and [&#8204;[[UFR_3-31_References#16|16]]]
or the mixed-time scale model of
[&#8204;[[UFR_3-31_References#11|11]]],
the latter being local, involving
no averaging and hence being computationally cheaper.  
no averaging and hence being computationally cheaper.  
The mixed-time scale model was adopted in the additional two grid-sensitivity  
The mixed-time scale model was adopted in the additional two grid-sensitivity  
Line 110: Line 134:
to emerge with greater clarity.
to emerge with greater clarity.


<div id="table1"></div>
{|align="center" border="1" cellpadding="10"
{|align="center" border="1" cellpadding="10"
!Computation!!Grid <math>(nx,ny,nz)</math>!!<math>x/H</math> (separation)!!<math>x/H</math> (reattachment)
!Computation!!Grid <math>(nx,ny,nz)</math>!!<math>x/H</math> (separation)!!<math>x/H</math> (reattachment)
Line 121: Line 146:




To explore the sensitivity to grid resolution, two simulations were performed (Table 1). The separated region (from <math>{x/H=0}</math> to <math>{x/H=4}</math>) was resolved using  
To explore the sensitivity to grid resolution, two simulations were performed
([[UFR_3-31_Test_Case#table1|Table 1]]).
The separated region (from <math>{x/H=0}</math> to <math>{x/H=4}</math>) was resolved using  
130 points in the streamwise direction for the coarse mesh, and 280 points in  
130 points in the streamwise direction for the coarse mesh, and 280 points in  
the fine mesh. In both grids, the first grid-cell centre is located below <math>{y^+=1}</math> at  
the fine mesh. In both grids, the first grid-cell centre is located below <math>{y^+=1}</math> at  
Line 129: Line 156:
respectively, downstream of <math>{x/H=0}</math>, with values of order 5 within the separation  
respectively, downstream of <math>{x/H=0}</math>, with values of order 5 within the separation  
zone. The maximum spanwise cell dimension, <math>{z^+}</math>, is 10 for both the coarse and fine  
zone. The maximum spanwise cell dimension, <math>{z^+}</math>, is 10 for both the coarse and fine  
simulations (cf. \cite{bentaleb2012large} for more details).
simulations (cf.  
[&#8204;[[UFR_3-31_References#3|3]]]
for more details).


Downstream of separation, the ratio of maximum cell dimension to the Kolmogorov length  
Downstream of separation, the ratio of maximum cell dimension to the Kolmogorov length  
Line 143: Line 172:
Turbulent boundary layers were imposed at the domain inlet (<math>{x/H = -7.4}</math>) on both  
Turbulent boundary layers were imposed at the domain inlet (<math>{x/H = -7.4}</math>) on both  
the upper and lower walls. To this end, data were taken from a precursor simulation,  
the upper and lower walls. To this end, data were taken from a precursor simulation,  
performed with the aid of a recycling method \cite{lund1998generation}, for a stretch  
performed with the aid of a recycling method
[&#8204;[[UFR_3-31_References#17|17]]],
for a stretch  
of a boundary layer developing in the actual channel. At a streamwise position at which  
of a boundary layer developing in the actual channel. At a streamwise position at which  
<math>{Re_\theta = 1190}</math>, every realization of the velocity field was stored over a period  
<math>{Re_\theta = 1190}</math>, every realization of the velocity field was stored over a period  
of <math>{1120H/U_{in}}</math>.  This data set was then used to explicitly prescribe the inlet  
of <math>{1120H/U_{in}}</math>.  This data set was then used to explicitly prescribe the inlet  
condition. Comparisons between lower-wall boundary-layer data returned by the present  
condition. Comparisons between lower-wall boundary-layer data returned by the present  
precursor simulation and DNS results of Schlatter ''et&nbsp;al'' \cite{schlatter2010assessment},  
precursor simulation and DNS results of Schlatter ''et&nbsp;al''
for a similar Reynolds number, <math>{Re_\theta = 1400}</math>, are shown on figure
[&#8204;[[UFR_3-31_References#25|25]]],
\ref{fig:inflowValidation}. The present plots pertain to a location immediately  
for a similar Reynolds number, <math>{Re_\theta = 1400}</math>, are shown on [[UFR_3-31_Test_Case#figure3|Fig. 3]].
downstream of the domain inlet, while Fig. 3 in \cite{bentaleb2012large} shows  
The present plots pertain to a location immediately  
distribution derived from the precursor simulation. The slight difference  
downstream of the domain inlet, while Fig. 3 in
relative to flat-plate boundary-layer data are due to a non-zero streamwise  
[&#8204;[[UFR_3-31_References#3|3]]]
shows distribution derived from the precursor simulation. The slight difference  
relative to flat-plate boundary-layer data is due to a non-zero streamwise  
pressure extending from the curved step upstream towards the inlet.
pressure extending from the curved step upstream towards the inlet.
<div id="figure3"></div>
{|align="center" border="0" width="750"
|[[Image:UFR3-31_fig3a.png|370px]]||[[Image:UFR3-31_fig3b.png|370px]]
|-
!align="center" colspan="2"|Figure 3: Mean flow and normal Reynolds stress profiles at the inflow. DNS data from [&#8204;[[UFR_3-31_References#25|25]]], for <math>{Re_\theta=1400}</math>.
|}


The time step corresponds to a maximum cell-CFL value of around 0.5. Simulations were  
The time step corresponds to a maximum cell-CFL value of around 0.5. Simulations were  
Line 167: Line 209:
contains three types of files, all extracted from the fine grid Large-Eddy simulation:  
contains three types of files, all extracted from the fine grid Large-Eddy simulation:  


#[[Media:UFR3-31_wallQuantities.dat|wallQuantities.dat]] contains boundary-layer data computed on the bottom wall. The file contains the following columns:<center>x, y_wall, C_p, tau_w, delta&lowast;, theta, delta_99</center>where x<math>{\equiv x/H}</math>, y_wall is given by&nbsp;[[UFR_3-31_Test_Case#eqn1|Eq.&nbsp;1)]],
and
<br/><center><math>{
\begin{align}
\mathrm{C\_p} &= \frac{p-p_0}{\frac{1}{2} \rho U_{in}^2}  \\
\mathrm{tau\_w} &= \frac{1}{Re_{in}} \left.\frac{\partial U}{\partial n}\right|_{y=0}  \\
\mathrm{delta*}&= \int_{y_{\min}}^{\infty} \left( 1 - \frac{U- U_{\min}}{U_{in}-U_{\min}} \right) dy \\
\mathrm{theta}&= \int_{y_{\min}}^{\infty} \frac{U- U_{\min}}{U_{in}-U_{\min}} \left( 1 - \frac{U- U_{\min}}{U_{in}-U_{\min}} \right) dy
\end{align}
}</math></center><br/>
::where <math>p_0</math> is taken at the bottom left corner of the domain and delta_99 is the distance to the wall for which <math>{U=0.99 U_{in}}</math>. <math>{U_{\min}}</math> is computed at each streamwise location by <math>{U_{\min} = \min(U(y),0)}</math>,  for <math>{y_{wall}\leq y\leq y_\infty}</math>.
<ol><li value="2">[[Media:UFR3-31_profiles.dat|profiles.dat]] contains the  velocity and Reynolds-Stress profiles in the vertical direction:<center>x, y, U, V, W, u', v', w', u'v', u'w', v'w'</center> for 28 different streamwise locations (<math>{x/H}</math> = -7.31 (first cell center from the inflow), -7, -6, -5, -4, -3, -2, -1, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 6, 7, 8, 9, 10, 11, 12, 13 and 14).</li><br/>
<li>'''balance.#.dat''' (with '''#'''&nbsp;=&nbsp;[[Media:UFR3-31_balance.k.dat|k]], [[Media:UFR3-31_balance.uu.dat|uu]], [[Media:UFR3-31_balance.vv.dat|vv]], [[Media:UFR3-31_balance.ww.dat|ww]] or [[Media:UFR3-31_balance.uv.dat|uv]]) contains the  budgets for <math>\ k</math>, <math>{\overline{uu}}</math>, <math>{\overline{vv}}</math>, <math>{\overline{ww}}</math> and  <math>{\overline{uv}}</math> respectively at the same 28 streamwise locations as in [[Media:UFR3-31_profiles.dat|profiles.dat]]. The different terms in these files are given by:</li></ol>
{| align="center" border="1" cellpadding="10"
|col||Abbrev.||<math>k</math>-equation||<math>{\overline{u_iu_j}}</math> equations
|-
|3||convec.||<math>{U_i \frac{\partial k}{\partial x_i}}</math>
|<math>{U_k \frac{\partial \overline{u_iu_j}}{\partial x_k}}</math>
|-
|4||prod.||<math>{-\overline{u_i u_j}\frac{\partial U_i}{\partial x_j}}</math>
|<math>{-\overline{u_i u_k}\frac{\partial U_j}{\partial x_k}-\overline{u_j u_k}\frac{\partial U_i}{\partial x_k}}</math>
|-
|5||t-diff.||<math>{-\frac{1}{2} \frac{\partial \overline{u_i u_i u_k}}{\partial x_k}}</math>
|<math>{-\frac{\partial \overline{u_i u_j u_k}}{\partial x_k}}</math>
|-
|6||p-diff.||<math>{ -\frac{1}{2} \frac{\partial \overline{u_ip}}{\partial x_i}}</math>
|<math>{- \left(\frac{\partial \overline{u_ip}}{\partial x_j} +  \frac{\partial \overline{u_jp}}{\partial x_i}\right)}</math>
|-
|7||p-str.||<math>{\ 0}</math>||<math>{\overline{p\frac{\partial u_i}{\partial x_j}} + \overline{p\frac{\partial u_j}{\partial x_i}}}</math>
|-
|8||v-diff.||<math>{\frac{1}{Re_{in}} \frac{\partial^2 k}{\partial x_k \partial x_k}}</math>
|<math>{\frac{1}{Re_{in}} \frac{\partial^2 \overline{u_iu_j}}{\partial x_k \partial x_k} }</math>
|-
|9||diss.||<math>{\ -(3) + (4) + (5) + (6) + (8)}</math>
|<math>{\ -(3) + (4) + (5) + (6) + (7) + (8)}</math>
|}
<br/>
===Reynolds-Averaged Navier-Stokes computations===
The objective of the RANS computations presented in this document is to guide
CFD practitioners and modellers for computing the present step curved flow.
All RANS computations have been performed by the author at CD-adapco with a commercial CFD software (STAR-CCM+,
ver. 7.02), using the default models, and the results should be
easily reproduced using other CFD software. Note that better
results have been obtained by using specific modification of the models
(e.g. low-Re damping modification of the <math>\omega</math> equation, or curvature correction),
but the objective here is not to show the best possible results, rather to
provide a set of comparison for the most classical formulations.


# wallQuantities.dat contains boundary-layer data computed on the bottom wall. The file contains the following columns:
The four different turbulence models are
the standard low-Re <math>{k-\varepsilon}</math> model of Launder and Sharma
[&#8204;[[UFR_3-31_References#14|14]]]
but with an alternative damping function
[&#8204;[[UFR_3-31_References#15|15]]],
the realizable <math>{k-\varepsilon}</math> model
[&#8204;[[UFR_3-31_References#26|26]]]
used in a ''two-layer approach'' (originally introduced in
[&#8204;[[UFR_3-31_References#24|24]]]
as an alternative to low-Reynolds number approach) in which the epsilon equation
is not solved near the wall but the turbulent length-scale is specified as
a function of the wall-distance,
the standard <math>{k-\omega}</math> model
[&#8204;[[UFR_3-31_References#30|30]]],
and the SST <math>{k-\omega}</math> model
[&#8204;[[UFR_3-31_References#19|19]]].
Detailed formulation of each model can be found in the STAR-CCM+ documentation
and is not repeated here.  
Note for instance that a low-Re damping of <math>\omega</math> (available as an option),
while very much applicable to this flow configuration, has not been selected
here as it might not be available as standard in all CFD packages.
The purpose of this part is purely to draw attention to some
particulars when computing this type of spatially developping flow using standard statistical
methods.


The grid is the same structured grid as that used in the fine LES described in
the previous section, only in 2D because of the spanwise homogeneity. The
resolution is thus quite fine, and several tests with coarser grids have 
shown that the results are grid-independent. The near-wall flow is fully resolved ,
with <math>{y^+\approx1}</math> at the inlet. Note that due to the very thin separated
region, high-Reynolds number models are not applicable for this test case:
computations with high-Reynolds number models return very small or no separation at all.


x, y_wall, C_p, tau_w, delta&lowast;, theta, delta_99
Spatially developping flows are very sensitive to the definition of the boundary conditions
and most notably to the inflow condition. For all computations, the outflow
boundary is treated as a pressure outlet, with a negligible impact on the results
if the position of the outflow is located at <math>{x/H>14}</math>. For the top wall, two options have been
tested: a free-slip or a no-slip condition. The only condition which needs specific
attention is the inflow. Different specifications of the inflow length-scale have been
tested, only the first (and simplest) one has been used in this report, but details
on setting-up the second one are given
[[UFR_3-31_Best_Practice_Advice#Application_Uncertainties|here]].
The velocity and turbulence-quantity profiles at the inlet are given by:


#the velocity vector (only <math>{U}</math> and <math>{V}</math> are necessary) are extracted from the LES database, at <math>{x/H=-7.31}</math>,
#the turbulent kinetic energy profile at <math>{x/H=-7.31}</math> is also extracted from the same database,
#The length-scale quantity is more problematic and cannot be directly extracted from the  database (except for ''a&nbsp;priori'' testing. The turbulent eddy-viscosity is over-estimated, which in turn leads to a too high production of the turbulent kinetic energy immediately  downstream of the inlet. It was thus decided to impose a constant value for both <math>\varepsilon</math> and <math>\omega</math> at the inlet, and to allow the flow to develop downstream. The choice  of the inlet value for these parameters will thus be dictated by only one requirement: matching the streamwise evolution of the wall-shear stress (Cf. Fig.&nbsp;[[UFR_3-31_Evaluation#figure7|7]]) from <math>{x/H=-4}</math>, ie. around <math>{5H}</math> upstream of the separation point.


where x<math>{\equiv x/H}</math>, y_wall is given by Eq. \ref{eq:ywall}, and
The second method is similar to that used for the large-eddy simulation: it requires
a pre-computation of a spatially developing flow on a flat-plate, to determine all
the correct quantities at the inflow. A third method would be to extract an eddy-viscosity
profile directly from the database and use it at the inflow.


<center><math>{
\begin{align}
\mathrm{C\_p} &= \frac{p-p_0}{\frac{1}{2} \rho U_{in}^2}  \\
\mathrm{tau\_w} &= \frac{1}{Re_{in}} \frac{\partial U}{\partial n}|y=0  \\
\mathrm{delta*}&= \int_{y_{\min}}^{\infty} \left( 1 - \frac{U- U_{\min}}{U_{in}-U_{\min}} \right) dy \\
\mathrm{theta}&= \int_{y_{\min}}^{\infty} \frac{U- U_{\min}}{U_{in}-U_{\min}} \left( 1 - \frac{U- U_{\min}}{U_{in}-U_{\min}} \right) dy
\end{align}
}</math></center>


where $P_0$ is taken at the bottom left corner of the domain
and delta\_99 is the distance to the wall for which $U=0.99 U_{in}$.
$U_{\min}$ is computed at each streamwise location by $U_{\min} = \min(U(y),0)$,
for $y_{wall}\leq y\leq y_\infty$.
<br/>
<br/>
----
----
{{ACContribs
{{ACContribs
|authors=Sylvain Lardeau
|authors=Sylvain Lardeau
|organisation=CD-adapco
|organisation=CD-adapco, London, UK
}}
}}
{{UFRHeader
{{UFRHeader

Latest revision as of 13:41, 12 February 2017

Flow over curved backward-facing step

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References

Semi-confined flows

Underlying Flow Regime 3-31

Test Case Study

Brief Description of the Study Test Case

The geometry under consideration is shown in Fig. 2. The rounded ramp of height is placed in a high-aspect-ratio duct with upstream height of . In the simulations, the flow is assumed to be spanwise homogeneous, with the spanwise slab being . The assumption of homogeneity is justified by the fact that the experimental ratio of duct depth to the step height was 38. In the experiment [‌31], tripped boundary layers were allowed to develop on both walls for a distance of about . The Reynolds number, based on and the inlet free-stream velocity , is 13,700. At , the computational inlet, the momentum-thickness Reynolds number is , and the boundary-layer thickness is .

UFR3-31 fig2.png
Figure 2: Flow domain and grid (every 8th point shown, fine mesh).


The step geometry is based on that used originally by Song and Eaton [‌27]. In order to enlarge the separated region, the height of the step was increased by a factor of 1.5. This adaptation was undertaken interactively with a parallel experiment by Zhang and Zhong [‌31]. The step shape is described by the following three relations, with the origin being the upstream edge of the ramp:


with , , and and for , for .

Test Case Experiments

The experimental working section is 5.49m long and 1.2m wide. The ramp height is 31.5mm, with an upstream section 660m long. The incoming flow velocity is fixed at 6.5m/s. Thirty-six pressure tapings are made along the ramp between and 90mm with a spacing of 10mm or 20mm. To prevent an interference with the LDA measurement on the central streamwise plane of the ramp, the pressure tapings are located 25mm off this plane. Combined with the existing pressure tapings on the ceiling of the working section downstream of the ramp (150mm in spacing), the coverage of the pressure measurement is extended further up to 600mm. The velocity at the inlet of the test section is monitored using a Pitot tube. The pressure measurement system consists of a low-range pressure transducer and a scani-valve with forty-eight ports. The pressure sensor has a pressure range of -2.5mbar to 2.5mbar and an accuracy of 0.25% full scale span. The system is remotely controlled by a computer equipped with a data acquisition card. The pressure data are sampled at 1kHz and 10,000 data points are used to produce the time-averaged pressure at each tapping.

A three-component Dantec LDA system is used to measure the boundary layer development around the ramp-down section in both the streamwise and the spanwise direction. The laser beam generated by the 5W Argon ion laser is separated by a series of beam separators and colour filters inside the optic transmitter box to produce three pairs of beams with a wavelength of 512nm, 488nm and 476nm respectively. The averaged sampling rate in the freestream is about 1kHz and a lower rate of about 100Hz is obtained in the near-wall region. The number of realisations used to obtain the time-averaged velocity and turbulence fluctuations is set at 10K, and the time threshold for acquiring the data at each point is 30 seconds. The Dantec 3D traverse specially designed for this LDA system has a resolution of 6.25μm in the three directions. The uncertainty of measurements is about 4 ×10-4 m/s in the freestream and 1.6 ×10-3 m/s in the near-wall region. All velocity components and all Reynolds stresses have been obtained, and the skin-friction is deduced from the LDA velocity measurements, assuming the first point being in the viscous sub-layer. The experimental study is described in detail in [‌31].

CFD Methods

Large-Eddy Simulations

Numerical treatment

The implicitly filtered LES momentum and continuity equations for incompressible flow were solved over a general non-orthogonal, boundary-fitted, multi-block finite volume mesh, using LES-STREAM (e.g. [‌7] for more details and related references). In total, the default mesh covering the duct-flow domain contained around 24 million nodes. The sensitivity to the resolution is addressed below. The variables are stored in a co-located manner. The solution is based on a fractional-step time-marching method, with the time derivative approximated by a third-order Gear scheme [‌7]. The fluxes are approximated by second-order centred approximations. Within the fractional-step algorithm, a provisional velocity field results from advancing the solution with the flux operators. This is then corrected through the pressure field by projecting the provisional solution onto a divergence-free velocity field. To this end, the pressure is computed as a solution to the pressure-Poisson problem with the aid of a V-cycle multigrid scheme. In order to suppress unphysical oscillations, associated with pressure-velocity decoupling, a practice equivalent to that introduced in [‌23] is adopted. Fishpool & Leschziner [‌7] demonstrate that the loss of accuracy associated with the smoothing introduced by this practice is minimal. The subgrid-scale stresses were approximated using either the dynamic Smagorinsky model of [‌10] and [‌16] or the mixed-time scale model of [‌11], the latter being local, involving no averaging and hence being computationally cheaper. The mixed-time scale model was adopted in the additional two grid-sensitivity simulations for the baseline case, and it is the one used for the results presented in the database. The rationale of using the simpler mixed-time scale model for the grid-sensitivity study is rooted primarily in the observation that this model gave somewhat lower time-averaged SGS viscosity values than the dynamic model, and thus allowed any grid-sensitivity specifically associated with variations in resolution to emerge with greater clarity.

Computation Grid (separation) (reattachment)
Fine (768,160,192) 0.83 4.36
Coarse (448,128,192) 0.87 4.21
Table 1: Grid resolution and separation and reattachment positions


To explore the sensitivity to grid resolution, two simulations were performed (Table 1). The separated region (from to ) was resolved using 130 points in the streamwise direction for the coarse mesh, and 280 points in the fine mesh. In both grids, the first grid-cell centre is located below at the bottom wall, whether is scaled with the local or inflow friction velocity. The maximum streamwise cell distance, , occurs close to the computational inlet and is around 35. This value drops to below 20 and 10 for the coarse and fine mesh, respectively, downstream of , with values of order 5 within the separation zone. The maximum spanwise cell dimension, , is 10 for both the coarse and fine simulations (cf. [‌3] for more details).

Downstream of separation, the ratio of maximum cell dimension to the Kolmogorov length (derived from the dissipation rate as part of the budget) is of order 6 for the fine grid and 11 for the coarse grid. The boundary layer on the upper wall was covered by a coarser grid than that on the lower wall, with the conditions between the wall and the wall-nearest location of around 25 being covered by a log-law wall function, and with about 17 grid nodes covering the boundary-layer thickness. This choice is justified by the fact that the upper boundary layer is far away from the region of interest, is only subjected to a weak adverse pressure gradient and is always fully attached.

Turbulent boundary layers were imposed at the domain inlet () on both the upper and lower walls. To this end, data were taken from a precursor simulation, performed with the aid of a recycling method [‌17], for a stretch of a boundary layer developing in the actual channel. At a streamwise position at which , every realization of the velocity field was stored over a period of . This data set was then used to explicitly prescribe the inlet condition. Comparisons between lower-wall boundary-layer data returned by the present precursor simulation and DNS results of Schlatter et al [‌25], for a similar Reynolds number, , are shown on Fig. 3. The present plots pertain to a location immediately downstream of the domain inlet, while Fig. 3 in [‌3] shows distribution derived from the precursor simulation. The slight difference relative to flat-plate boundary-layer data is due to a non-zero streamwise pressure extending from the curved step upstream towards the inlet.


UFR3-31 fig3a.png UFR3-31 fig3b.png
Figure 3: Mean flow and normal Reynolds stress profiles at the inflow. DNS data from [‌25], for .


The time step corresponds to a maximum cell-CFL value of around 0.5. Simulations were run for 300 time units, one time unit being , before statistical sampling commenced. The flow was then computed for further 1030 time units during which data for various statistical properties were gathered. Statistical convergence was increased by spanwise averaging.

Data provided

In the database provided for this test case, all quantities have been non-dimensionalized using the step height and the inflow free-stream velocity . The database contains three types of files, all extracted from the fine grid Large-Eddy simulation:

  1. wallQuantities.dat contains boundary-layer data computed on the bottom wall. The file contains the following columns:
    x, y_wall, C_p, tau_w, delta∗, theta, delta_99
    where x, y_wall is given by Eq. 1),

and



where is taken at the bottom left corner of the domain and delta_99 is the distance to the wall for which . is computed at each streamwise location by , for .


  1. profiles.dat contains the velocity and Reynolds-Stress profiles in the vertical direction:
    x, y, U, V, W, u', v', w', u'v', u'w', v'w'
    for 28 different streamwise locations ( = -7.31 (first cell center from the inflow), -7, -6, -5, -4, -3, -2, -1, 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 6, 7, 8, 9, 10, 11, 12, 13 and 14).

  2. balance.#.dat (with # = k, uu, vv, ww or uv) contains the budgets for , , , and respectively at the same 28 streamwise locations as in profiles.dat. The different terms in these files are given by:


col Abbrev. -equation equations
3 convec.
4 prod.
5 t-diff.
6 p-diff.
7 p-str.
8 v-diff.
9 diss.


Reynolds-Averaged Navier-Stokes computations

The objective of the RANS computations presented in this document is to guide CFD practitioners and modellers for computing the present step curved flow. All RANS computations have been performed by the author at CD-adapco with a commercial CFD software (STAR-CCM+, ver. 7.02), using the default models, and the results should be easily reproduced using other CFD software. Note that better results have been obtained by using specific modification of the models (e.g. low-Re damping modification of the equation, or curvature correction), but the objective here is not to show the best possible results, rather to provide a set of comparison for the most classical formulations.

The four different turbulence models are the standard low-Re model of Launder and Sharma [‌14] but with an alternative damping function [‌15], the realizable model [‌26] used in a two-layer approach (originally introduced in [‌24] as an alternative to low-Reynolds number approach) in which the epsilon equation is not solved near the wall but the turbulent length-scale is specified as a function of the wall-distance, the standard model [‌30], and the SST model [‌19]. Detailed formulation of each model can be found in the STAR-CCM+ documentation and is not repeated here. Note for instance that a low-Re damping of (available as an option), while very much applicable to this flow configuration, has not been selected here as it might not be available as standard in all CFD packages. The purpose of this part is purely to draw attention to some particulars when computing this type of spatially developping flow using standard statistical methods.

The grid is the same structured grid as that used in the fine LES described in the previous section, only in 2D because of the spanwise homogeneity. The resolution is thus quite fine, and several tests with coarser grids have shown that the results are grid-independent. The near-wall flow is fully resolved , with at the inlet. Note that due to the very thin separated region, high-Reynolds number models are not applicable for this test case: computations with high-Reynolds number models return very small or no separation at all.

Spatially developping flows are very sensitive to the definition of the boundary conditions and most notably to the inflow condition. For all computations, the outflow boundary is treated as a pressure outlet, with a negligible impact on the results if the position of the outflow is located at . For the top wall, two options have been tested: a free-slip or a no-slip condition. The only condition which needs specific attention is the inflow. Different specifications of the inflow length-scale have been tested, only the first (and simplest) one has been used in this report, but details on setting-up the second one are given here. The velocity and turbulence-quantity profiles at the inlet are given by:

  1. the velocity vector (only and are necessary) are extracted from the LES database, at ,
  2. the turbulent kinetic energy profile at is also extracted from the same database,
  3. The length-scale quantity is more problematic and cannot be directly extracted from the database (except for a priori testing. The turbulent eddy-viscosity is over-estimated, which in turn leads to a too high production of the turbulent kinetic energy immediately downstream of the inlet. It was thus decided to impose a constant value for both and at the inlet, and to allow the flow to develop downstream. The choice of the inlet value for these parameters will thus be dictated by only one requirement: matching the streamwise evolution of the wall-shear stress (Cf. Fig. 7) from , ie. around upstream of the separation point.

The second method is similar to that used for the large-eddy simulation: it requires a pre-computation of a spatially developing flow on a flat-plate, to determine all the correct quantities at the inflow. A third method would be to extract an eddy-viscosity profile directly from the database and use it at the inflow.





Contributed by: Sylvain Lardeau — CD-adapco, London, UK

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References


© copyright ERCOFTAC 2024