CFD Assignment 3 Lid Driven cavity flow -Due date: Monday, 7th January, 2020
1. The present task is to provide numerical solution of a steady flow within a square cavity, where the
top lid is moving to the right at a constant velocity. The following figure shows a sample solution for
Re=100 with central difference scheme. The Reynolds number is defined as Re = ρUlid L/µ, where Ulid
is the top lid velocity (1m/s) and L is the cavity height.
V=0, U=1 1 m/s
0.8
U=V=0
U=V=0
0.6
Y
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2
U=V=0 X
The flow field can be solved by the continuity and Navier-Stokes equations (ρ = constant):
∂ρU ∂ρV
+ =0
∂x ∂y
∂ρU ∂ρU U ∂ρV U ∂P ∂ 2U ∂ 2U
+ + =− + µ[ 2 + ]
∂t ∂x ∂y ∂x ∂x ∂y 2
∂ρV ∂ρU V ∂ρV V ∂P ∂ 2V ∂ 2V
+ + =− + µ[ 2 + ]
∂t ∂x ∂y ∂y ∂x ∂y 2
(a) Please compute the cavity flows with Reynolds number being 100, 1000 and 5000.
(b) Please use mesh size at least 81x81 and 161x161.
(c) Please also compare your results with the benchmark solutions from Gihia et al. (1982) for u(y) at
x=0.5 and v(x) at y=0.5.
Ghia U, Ghia KN, Shin CT. High-Resolutions for incompressible flow using the Navier-Stokes
equations and a multigrid method. J Comput Phys 1982;48:387.
(d) Please compare the accuracy of the central difference, QUICK and MUSCL schemes.
2. Programming languages
(a) Please use programming language at your own choice.
(b) Please list your program in the report
(c) You can use artificial compressibility method, fractional step or SIMPLE method.
3. Please compile your results into a report. No hand written material should be included in the report.
The report should be submitted through the turnitin website.
The report should contain
(a) Problem descriptions.
(b) Introduction of methodology adopted
(c) Results and discussions
(d) Conclusion
(e) List of programs