2D Backward-Facing Step: Validation (SA)

validation

Velocity field over the backward-facing step
Velocity-magnitude field through the step region at ReH=36,000, M=0.128; the separated shear layer springs from the step lip and the recirculation bubble fills the corner before reattachment on the lower floor.

Overview

The backward-facing step is a canonical test of how a turbulence model handles separation and reattachment. As the incoming boundary layer passes the step lip it detaches into a free shear layer, a recirculation bubble forms in the corner, and the flow reattaches on the floor some distance downstream. The length of that bubble, together with the recovery of wall pressure and skin friction behind it, probes how a RANS closure predicts both the spreading of the separated shear layer and the adverse-pressure recovery that brings the flow back to the wall.

This case reproduces the Driver & Seegmiller (1985) experiment, a long-standing benchmark on the NASA Turbulence Modeling Resource. Luminary's Spalart–Allmaras solution is compared to the measured lower-wall pressure and skin friction and to the NASA CFL3D reference code computed on the same grid.

Problem Setup

The step height is H=1, the free-stream Mach number is M=0.128, and the Reynolds number based on step height is ReH=36,000. The fluid is ideal-gas air with Sutherland viscosity. The inlet is placed 130H upstream of the step so that a uniform inflow develops into the correct turbulent boundary layer by the time it reaches the lip, no inflow profile is prescribed. The outlet uses a back-pressure outflow set slightly above free-stream to match the mild streamwise pressure gradient of the tunnel. The step face, the floor upstream and downstream of the step, and the opposite tunnel wall are all no-slip adiabatic walls; the spanwise faces are symmetry planes for this quasi-two-dimensional computation.

The mesh is the NASA TMR's own published structured grid for this geometry, a point-matched multiblock family stitched here into a single conformal grid. Using that grid means Luminary and CFL3D are solving on identical points, so any difference is a genuine model-and-solver difference rather than a grid effect. The discretization is second order and the solution is advanced until the mean flow reaches a statistically steady state.

Quantities of Interest

The validation quantities are the lower-wall pressure coefficient Cp(x/H) and skin-friction coefficient Cf(x/H), both reported in step heights along the floor. The pressure coefficient is re-datumed to zero far downstream (near x/H=40), the same convention used for the experiment and CFL3D, so the curves are directly comparable. The skin friction is especially informative: the point where Cf changes sign locates flow reattachment, which the experiment places at xr/H6.26. Together these distributions test whether the turbulence model captures the size of the recirculation zone and the strength of the downstream pressure recovery.

Sources

NASA Langley Turbulence Modeling Resource, "2D Backward Facing Step" (SA model).: https://tmbwg.github.io/turbmodels/backstep_val_sa.html
Reference lower-wall Cp and Cf from the NASA CFL3D code at ReH=36,000, M=0.128.
D. M. Driver and H. L. Seegmiller, "Features of a Reattaching Turbulent Shear Layer in Divergent Channel Flow," AIAA Journal, 23(2):163–171, 1985.: https://tmbwg.github.io/turbmodels/backstep_val.html
Measured lower-wall Cp, Cf and reattachment length xr/H6.26.

Results

Lower-wall pressure coefficient Cp(x/H)

Lower-wall pressure coefficient \( C_p(x/H) \)
Lower-wall Cp versus x/H. Luminary's pressure coefficient is shifted to zero near x/H=40 to match the reference datum.

Reference: NASA Langley Turbulence Modeling Resource, "2D Backward Facing Step" (SA model). (CFL3D (SA), lower-wall Cp): https://tmbwg.github.io/turbmodels/backstep_val_sa.html

Lower-wall skin-friction coefficient Cf(x/H)

Lower-wall skin-friction coefficient \( C_f(x/H) \)
Lower-wall Cf versus x/H. The sign change marks flow reattachment; the experiment reports xr/H6.26.

Reference: NASA Langley Turbulence Modeling Resource, "2D Backward Facing Step" (SA model). (CFL3D (SA), lower-wall Cf): https://tmbwg.github.io/turbmodels/backstep_val_sa.html

Solver configuration

features exercised
VerificationValidation2DRANSSpalart–AllmarasIdeal gasSteadyEnergy equation