Limitations on time the step for lid driven cavity flow
 

I'm following these slides

http://www3.nd.edu/~gtryggva/CFD-Cou...-Lecture-5.pdf

They formulate N-S equations usin stream, vorticity functions, I'm want to know what kind of analysis they do to conclude the time step restriction showed in the page 5. Its very important for me to know how it is done, thanks!

https://lh3.googleusercontent.com/mG...2=w367-h157-no

It stems from the numerical stability analysis for the explicit Forward Time Central Space (FTCS) applied on a 2D convection-diffusion equation. Note that the first constraint is exact only for the pure diffusive case (no convection). The method is unconditionally unstable for the pure convective case. In the mixed case the constraint involves a more complex functional law in the (cfl, Reh) plane. You can find details in any good CFD textbook.

I'm curious about your answer, because as I know it is a Von Newman stability analysis, however this problem is not linear so the Von Newman analysis can not be applied here. I'm missing something?

Yes. The Von Neumann analysis is applied for the linear case (as well as the matrix analysis you can perform by using the Gershgorin locus method). The constraint in those slides are valid for the linear equation. Usually, they are roughly used also for the non linear equation.

P.S.:the second constraint must have the ratio to the mesh size h but, as I wrote, the case for the pure convective case is unstable. The key for a wiggles-free solution is to work with the Reh<2 condition.

sorry to be annoying, but do you know any paper or article where this specific analysis is done?

Fletcher: Computational Techniques for Fluid Dynamics, Vol. 2, pp. 335 cites a more general condition for the 2D equations. He actually cites this:

http://link.springer.com/book/10.100...-3-642-85952-6

pp. 183 as original source. If you don't trust them, I swear, I did it by myself (straightforward Neumann on the equivalent scalar equation) and it is correct ;).

It is ():

Thanks a lot! :) :) :)

you can read any good textbook of both CFD and numerical anlysis...
as example, http://dl.iranidata.com/book/daneshg...idata.com).pdf

Sbafiny, I've checked the books you mentioned but again I don't see they do any analysis to arrive to the time step condition.

You mentioned something about "straightforward Neumann on the equivalent scalar equation", may you point me to some source about these formulation and the analysis done over it?

That's what i meant by all those "cites" and "swear". If you can handle the italian language between the 14 equations on 2 pages (i guess so), here it is attached my old analysis. As you can see, it is pretty straightforward Neumann analysis.

I strongly suggest to plot the amplification factor for the value = 1 in the (Reh, cfl) plane.

This issue is not as simple as it would appear...


https://estudogeral.sib.uc.pt/bitstr...20analysis.pdf