VVF15 - Incompressible laminar lid-driven cavity

Test case
SVTEST193

Description

This test case examines the steady, incompressible, laminar lid-driven cavity flow. Despite its simple geometry, the lid-driven cavity flow retains a rich fluid flow structure manifested by multiple counter rotating recirculation regions near the corners of the cavity depending on the Reynolds number.

Geometry

The geometry is an extruded square: 100 mm (+ X-dir) x 100 mm (+ Y-dir) x 2 mm (+Z-dir).

Simulation Model

This model uses the Advanced Thermal-Flow solution type.

The swept mesh is made of hexahedral elements. The following mesh controls are used without any bias:

  • 200 elements along the + Y-dir
  • 200 elements along the + X-dir
  • 2 elements along the + Z-dir

Total number of elements: 80,000

Total number of nodes: 121,203

The fluid is modeled using incompressible air with the following properties:

  • Mass density: ρ = 1.21 kg/m3
  • Thermal conductivity: k = 1 W/mm·C
  • Dynamics viscosity: µ = 1.85 x 10−5 Pa·s
  • Specific heat at constant pressure: Cp = 1 J/kg·K

The gas constant is kept blank in order to simulate a truly incompressible flow.

The following boundary conditions are applied:

  • Symmetry Plane on the +Z and −Z walls
  • Flow Surface: Boundary Flow Surface on the +X, −X and −Y walls using No Slip Wall conditions
  • Flow Surface: Boundary Flow Surface on the + Y wall with specified velocity Ulid = 15.29 mm/s

The following solution option is set:

  • Turbulence Model: None (Laminar Flow)

The following solver parameters are selected:

  • 3D Flow Solver: Maximum Residuals = 1e-6
  • Advection Schemes Momentum : Second-order (CDS)
  • Limiter : Automatic; Momentum Limiter Stabilization : General Convective Boundedness Condition

Theory

The laminar incompressible flow in a square cavity whose top wall moves with a uniform velocity in its own plane has served extensively as a model problem for testing and evaluating numerical techniques. The multiple recirculation regions in the vicinity of the corners of the cavity provides for interesting flow features to test and benchmark numerical methods.

Results

The following section presents the results of a dimensionalized u-velocity and v-velocity profiles as a function of respectively the y and x centerline of the fluid domain.

Steady-state centerline u-velocity at Re = 100:

Steady-state centerline v-velocity at Re = 100:

The figures above compare the theoretical data with computed data. A grid size of 200x200x2 was selected yielding a maximum percent difference of 0.69% in streamwise direction and 0.37% in the crosswise direction with the steady state results of Ghia et al.[18].