[mulchah]

Hi!

I am having some trouble completing my Matlab code for Couette flow (parallel flow between two plates, the top one moving at a specified velocity) using the pressure correction method, as described in section 9 of John Anderson's book "CFD - Basics with Applications".

The trouble is that while my code manages to complete the suggested 300 iterations to solve for the velocity and pressure fields, the mass source term I get at the end of the iterations is d = -0.350 x 10^-5 slug/ft^3.s, different from Anderson's value of d = -0.172 x 10^-5 slug/ft^3.s. Related to this, the horizontal component of velocity, u, looks slightly different after 4 iterations when compared with Anderson's values (there is a graph on page 444 of his book). The vertical component of velocity, v, looks similar after 4 iterations though. This makes me think the horizontal component is incorrect. However, I can't find anything wrong with it in my code!

Having said all that, none of this makes a huge difference to the flowfield as the horizontal component of velocity, u, still looks very linear when it has converged but nevertheless, I would like to know what (small) thing I am doing wrong.

Can you please let me know the error of my ways? My code is attached for your perusal. Thanks in advance!