Flow over curved backward-facing step

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. \ref{fig:geometry}. The rounded ramp of height ${H}$ is placed in a high-aspect-ratio duct with upstream height of ${8.52H}$ . In the simulations, the flow is assumed to be spanwise homogeneous, with the spanwise slab being ${3.7H}$ . 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 \cite{zhang2010experimental}, tripped boundary layers were allowed to develop on both walls for a distance of about ${25H}$ . The Reynolds number, based on $H$ and the inlet free-stream velocity ${U_{in}}$ , is 13,700. At ${x/H=-7.4}$ , the computational inlet, the momentum-thickness Reynolds number is ${Re_{\theta }=1192}$ , and the boundary-layer thickness is ${\delta _{99}=0.83H}$ .

The step geometry is based on that used originally by Song and Eaton \cite{song2004reynolds}. 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 \cite{zhang2010experimental}. The step shape is described by the following three relations, with the origin ${x/H=0}$ being the upstream edge of the ramp:

{\begin{aligned}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}&=R_{2}-{\sqrt {R_{2}^{2}-(3-x)^{2}}}\qquad \mathrm {for} \qquad 2.835<x/H\leq 2.937\\\end{aligned}}

with ${R_{1}=4.03}$ , ${\ R_{2}=0.333}$ , ${\ x_{2}=3.449}$ and ${\ y_{2}=1.936}$ and ${\ y_{wall}=1}$ for ${\ x/H<0}$ , ${\ y_{wall}=0}$ for ${\ x/H>2.937}$ .

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 ${x=-400mm}$ 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.

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. \cite{fishpool2009stability} 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 \cite{fishpool2009stability}. 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 \cite{rhie1983numerical} is adopted. Fishpool \& Leschziner \cite{fishpool2009stability} 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 \cite{germanoetal91} and \cite{lilly92} or the mixed-time scale model of \cite{inagakietal05}, 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 (nx,ny,nz)	$x/H$ (separation)	$x/H$ (reattachment)
Fine	(768,160,192)	0.83	4.36
Coarse	(448,128,192)	0.87	4.21

\begin{table} \begin{center} \begin{tabular}{c|ccc} Computation & Grid $(nx,ny,nz)$ & $x/H$ (separation) & $x/H$ (reattachment) \\ \hline \hline Fine & (768,160,192) & 0.83 & 4.36 \\ Coarse & (448,128,192) & 0.87 & 4.21 \\ \end{tabular} \caption{\label{tab:domain} Grid resolution and separation and reattachment positions} \end{center} \end{table}

Contributed by: Sylvain Lardeau — CD-adapco