[maysmech]

Dear Foamers,
i used icoFoam solver for cavity. when the reynolds number is below 5000 the results are good. but for reynolds 5000 and 10000 results are different from Ghia benchmark. i promoted orders of discretisations in fvSchemes but results were same.
it is wondering because for reynolds below Re=10000 the cavity case is laminar and icoFoam should leads to true results.
WHAT IS THE PROBLEM?

[stevenvanharen]

Hi maysam,

I am assuming you are referring to the lid-driven cavity flow.

i think both Reynolds number will not result in steady laminar flow. The OpenFoam user guide switches to pisoFoam for Re = 10^4. Botella and Peyret found unsteady behaviour at least for Re=9000.

What is your source for assuming laminar flow for Re < 10^4?

[maysmech]

I examined cavity lid driven with pisoFoam RAS but results are same as icoFoam. for Re= 10000 and 5000 results are far from Ghia benchmark

[stevenvanharen]

Ok, I will look at the Ghia paper this afternoon.

In the mean time, in what way do your results differ from the Ghia paper? In what way are velocity plots and pressure plots different?

[maysmech]

Results for Re=10000.
velocity diagram over horizontal and vertical centreline.
U=1, D=1, nu=10^-4, time 100 sec


[IMG]file:///home/maysam/Desktop/U.gif[/IMG]

[stevenvanharen]

Ok, I am not a 100 percent sure about what it is but I have some thoughts:

- Ghai et al use upwinding, are you using upwinding too in icoFoam? Could it be that your numerical diffusion is less strong than theirs? The fact that your velocity profile is less smooth suggests this.
- the fact that you use upwinding could explain the laminar solution for Re=10000.
- numerical diffusion is less dominant at low Reynolds numbers, therefore your error increases for large Reynolds number
- I suggest you use the paper by Botella & Peyret (1998), their method is much more accurate (spectral) and they also compare to Ghai et. al., unfortunately only up to Re = 1000.

Hope this helps!

[maysmech]

Originally Posted by stevenvanharen

Ok, I am not a 100 percent sure about what it is but I have some thoughts:

- Ghai et al use upwinding, are you using upwinding too in icoFoam? Could it be that your numerical diffusion is less strong than theirs? The fact that your velocity profile is less smooth suggests this.
- the fact that you use upwinding could explain the laminar solution for Re=10000.
- numerical diffusion is less dominant at low Reynolds numbers, therefore your error increases for large Reynolds number
- I suggest you use the paper by Botella & Peyret (1998), their method is much more accurate (spectral) and they also compare to Ghai et. al., unfortunately only up to Re = 1000.

Hope this helps!

Dear Steven,
Thanks for your suggestions.
i examined icoFoam for Re=100, 400 and 1000 and there was no difference between them and Ghia et al.
but there was problem in 5000 and 10000 and it is solved by considering two things:
- I should use higher order discretisation method i used vanleer.
- I shouldn't use initial value u=0 for all domain. i used steady state velocity and pressure of Re=1000 for Re= 5000 and also used steady satate values of Re=5000 for 0 folder of 10000. i think it is an important work for high reynolds problems.
Best wishes,
Maysam

[stevenvanharen]

ok, first of all congratulations on solving the problem!

But your solution should be independent of the initial velocity field!

In the previous simulation which yielded the erroneous results how did you stop the simulation? A simulation with a higher Reynolds number should take longer to converge to a steady state solution (in physical time).

[maysmech]

for Re=1000 after 20 sec the problem became steady and i let 50 sec for 5000 and 100sec for 10000. i thought the time was sufficient to become steady because i didn't see any changes in several last time steps by eye.
maybe it was not enough i will try it with more time because a sudden shock in high reynolds leads to late steady and putting steady state values of lower Re to 0 values of higher Re makes less run time.

[Ali Raza Abid]

Dear Foamers

I am unable to get comparable resutls for lid driven cavity case at high Re 10,000
for openFoam vs Ghia et al.

I tried using pisoFoam and using mesh 257x257 with Courant no. kept below 1.

I even tried using vanLeer scheme which some friend mentioned here but the results are not improved.

looing for suggestions please.