[satty_00]

Hi everybody,

I've been trying to solve the 2D Lid Driven cavity problem in C using SIMPLEC algorithm. I'm able to attain a converged solution but not Benchmark solutions for the given Reynolds number.
Can anyone help me in correcting the Bug, or identify the problem in my Code.

V-Velocity.jpg

U-Velocity_RNO100.jpg

C Code.txt

[tas38]

It would be helpful if you could elaborate on the code (finite difference vs. finite volume, staggered vs. cell-centered variable arrangement, method of convective flux discretization, method of gradient computation etc.). Additionally, if you could narrow down possible locations in the code that are more suspect than others, that would be helpful.

Also, it would be helpful if you could show the velocity and pressure fields following one single outer iteration and without the first velocity correction applied. Please also include the color legends on your contour plots.

[satty_00]

I'm Using FVM for calculation. I've attached an Paint image of how i've discretized the Central and Staggered Meshes. I've used SOR (TDMA + Gauss Seidal Method) to solve for U- Velocity, V- Velocity and Pressure corrections.

I've used the power law scheme to calculate the coefficients. I've previously solved Convection Diffusion equations using the same scheme and have obtained the benchmark solutions.

And I've used Node based Computation.

So I suspect the errors to be in the pressure correction part or Velocity correction part or the way I iterate.

Grid.png

And thank you very much for the quick response mate !