Simulation of flow over a backward-facing step.
The mean velocity profile obtained from a boundary layer simulation by Spalart (1988) is imposed at the inlet at \(Re_{\theta}= 667\). Random velocity fluctuations \(u'\), \(v'\), and \(w'\) are superimposed on this profile according to a variant of the method of Lee et al. (1992). The fluctuations are prescribed such that, at the inlet, the turbulence intensities and Reynolds shear stress of Spalart's data are also duplicated. A convective boundary condition (Pauley et al., 1988) is imposed at the exit. The streamwise domain consists of an entry section of length \(10h\) prior to the step and a \(20h\) post-expansion section, where \(h\) is the step height. The vertical dimensions before and after the expansion are \(H_1 = 5h\) and \(H_2 = 6h\) which give an expansion ratio \(ER\) of 1.20 (see figure 1).
The Reynolds number, based on \(h\) and the mean inlet free stream velocity \(U_0\), is \(Re_h = 5100\).
Contour plots of the instantaneous spanwise vorticity \(\omega_z\) on a typical vertical plane are shown in figure 2. The vorticity is normalised by \(U_0/h\). A free-shear layer spreading from the step and interacting with the lower wall near the mean reattachment location, \(x_R = 6h\), is discernible.
The basic statistical quantities are compared in Le et al (1997) to results from experiments by Jovic and Driver (1994), showing generally good agreement.
The Navier-Stokes equations are discretized using a finite difference method on a staggered mesh. Uniform mesh spacing is applied in the streamwise (\(x\)) and spanwise (\(z\)) directions. In the vertical (\(y\)) direction, non-uniform mesh is employed with mesh refinement at the wall and near the step. The fractional step method from Le and Moin (1991) is used for time advancement. The Navier-Stokes equations are first advanced using a second-order semi-implicit method without the pressure terms. The pressure is then calculated by solving the Poisson equation, and the velocities are then corrected to satisfy the continuity equation.
The spanwise dimension, \(W\), is \(4h\) where periodic boundary conditions are imposed. The simulation uses \(770 \times 194 \times 66\) grid points in the streamwise, wall normal, and spanwise directions, respectively.
The computation used 13 megawords of memory and required approximately 55 CPU seconds per time step on a single processor CRAY Y-MP at a rate of 186 mflops. Statistical quantities are averaged over the spanwise direction and time. About 1100 CPU hours were required to obtain an adequate statistical sample. The total computational time corresponds to approximately 4.5 flow-through times.
Data available includes:
Sample plots of selected quantities are available.
The data can be downloaded as compressed archives from the links below, or as individual files.
Location | Mean velocity and Reynolds stress profiles | \(\overline{u^2}\) Budget | \(\overline{v^2}\) Budget | \(\overline{w^2}\) Budget | \(\overline{uv}\) Budget | \(k\) Budget |
---|---|---|---|---|---|---|
\(x/h=-3\) | x-181.dat | rs11-181.dat | rs22-181.dat | rs33-181.dat | rs12-181.dat | rskk-181.dat |
\(x/h= 4\) | x-360.dat | rs11-360.dat | rs22-360.dat | rs33-360.dat | rs12-360.dat | rskk-360.dat |
\(x/h= 6\) | x-411.dat | rs11-411.dat | rs22-411.dat | rs33-411.dat | rs12-411.dat | rskk-411.dat |
\(x/h=10\) | x-513.dat | rs11-513.dat | rs22-513.dat | rs33-513.dat | rs12-513.dat | rskk-513.dat |
\(x/h=15\) | x-641.dat | rs11-641.dat | rs22-641.dat | rs33-641.dat | rs12-641.dat | rskk-641.dat |
\(x/h=19\) | x-744.dat | rs11-744.dat | rs22-744.dat | rs33-744.dat | rs12-744.dat | rskk-744.dat |
Skin friction and boundary layer data are available in the file stat-inf.dat. The file readme.txt contains further information on the files.
Indexed data:
case031 (dbcase, confined_flow) | |
---|---|
case | 031 |
title | Backward-Facing Step |
author | Le, Moin |
year | 1992 |
type | DNS |
flow_tag | 2d, separated, varying_cross_section |