Lab 7
Lab 7
Lab assignment 7:
Total points: 20
Date of release: March 19th; Lab date: March 21st
Due date for report and code (please submit the matlab codes and short report with results on
moodle only): April 2nd; Evaluation date: April 4th
In this lab, we will continue development of SIMPLE scheme on a staggered grid. The 2D flow
domain is given by a rectangular box of size where (See Figure 1 below).
The number of grid cells are and in both directions, and uniform grid is being
used. The flow field satisfies the discretized steady Stokes equation (i.e. no unsteady or convection
terms are present):
and (1)
in which Reynolds number is . The boundary conditions are:
, (inlet condition), and (at outlet) for all
. Also (no slip at top and bottom walls) for all . Here:
, and
The variables are located on a staggered grid, where
,
are the location of the nodes at cell centers. The cell faces are located at:
,
For a cell with given index , we denote discretized guess and predicted fields as:
, ,
and the corresponding discretized corrected fields as .
Q1 (2 points): Clearly state the discretized prediction equations in terms of for
interior nodes.
Q2 (2 points): Clearly state the correction equations for in
terms of pressure correction for interior nodes.
Q3 (2 points) Clearly derive and state the discretized Poisson equation for
Q4 (4 points) Clearly state the transformations required for the coefficients
of the predictor equations (for ) at the 4 boundaries. (Note:
We have already derived this for pressure Poisson equation in Lab 6. At , we can think of
“extrapolating” the velocities to ghost cell values, e.g. to transform
the coefficients.)
Q5 (10=1+4+2+2+1 points) Solve for discretized fields which satisfy Eqn 1, using SIMPLE
scheme:
(a) You will have to use very small under-relaxation factor (i.e. around 0.001) in the outer iteration.
Clearly state the values of the under-relaxation factor you are using, along with the number of
iterations required for convergence.
(b) Plot isocontours of (3 graphs)
(c) Make a vector plot of the velocity field. You have to interpolate velocity components at the
pressure nodes for this (1 graph)
(d) Plot isocontours of (1 graph)
(e) Report the value of residuals in Eqn 1 after convergence has been achieved
Figure 1: Flow domain and boundary conditions