CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Main CFD Forum (https://www.cfd-online.com/Forums/main/)
-   -   Lid-Driven Cavity Fortran Code works for coarse mesh but not fine mesh (https://www.cfd-online.com/Forums/main/255789-lid-driven-cavity-fortran-code-works-coarse-mesh-but-not-fine-mesh.html)

praveenu April 29, 2024 13:59
Lid-Driven Cavity Fortran Code works for coarse mesh but not fine mesh
 
1 Attachment(s)
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 April 29, 2024 15:29
Quote:
Originally Posted by praveenu (Post 868490)
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 April 29, 2024 16:03
1 Attachment(s)
Quote:
Originally Posted by FMDenaro (Post 868498)
Re number? Solution on the 11x11 grid?

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

FMDenaro April 29, 2024 16:52
Quote:
Originally Posted by praveenu (Post 868500)
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.


All times are GMT -4. The time now is 17:04.