CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Main CFD Forum (https://www.cfd-online.com/Forums/main/)
-   -   2D Lid Driven Cavity Problem in C (https://www.cfd-online.com/Forums/main/133621-2d-lid-driven-cavity-problem-c.html)

satty_00 April 17, 2014 12:13

2D Lid Driven Cavity Problem in C
 
3 Attachment(s)
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.

Attachment 30235

Attachment 30236

Attachment 30237

tas38 April 19, 2014 10:03

more information would be helpful
 
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 April 19, 2014 10:49

Details
 
1 Attachment(s)
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.

Attachment 30262

And thank you very much for the quick response mate !


All times are GMT -4. The time now is 00:37.