[praveenu]

Hi all,
I've written a FORTRAN code to solve the classical lid driven cavity problem using SMAC algorithm. The time marching is done using first order Euler scheme and spatial discretized using second order central difference.



The code seems to work fine for a gird of 11x11 but blows up for 51x51. I've tried reducing the CFL number but it doesn't seem to help.


Any advise?


thanks in advance!


PS: the code is attached

[FMDenaro]

Originally Posted by praveenu

Hi all,
I've written a FORTRAN code to solve the classical lid driven cavity problem using SMAC algorithm. The time marching is done using first order Euler scheme and spatial discretized using second order central difference.



The code seems to work fine for a gird of 11x11 but blows up for 51x51. I've tried reducing the CFL number but it doesn't seem to help.


Any advise?


thanks in advance!


PS: the code is attached

Re number? Solution on the 11x11 grid?

[praveenu]

Originally Posted by FMDenaro

Re number? Solution on the 11x11 grid?

Hi, the Re = 100, I've attached the tecplot file of the last time step.

[FMDenaro]

Originally Posted by praveenu

Hi, the Re = 100, I've attached the tecplot file of the last time step.

This picture means nothing. You need to check at the final steady state. Use the quantitative comparison, velocity profiles, vorticity and so on.