CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Why icoFoam solver results are not true for cavity (10000>Re>5000) ?

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   August 9, 2010, 11:57
Question Why icoFoam solver results are not true for cavity (10000>Re>5000) ?
  #1
Senior Member
 
maysmech's Avatar
 
Join Date: Jan 2010
Posts: 347
Blog Entries: 2
Rep Power: 8
maysmech is on a distinguished road
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?
maysmech is offline   Reply With Quote

Old   August 10, 2010, 10:44
Default
  #2
Senior Member
 
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 6
stevenvanharen is on a distinguished road
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?
stevenvanharen is offline   Reply With Quote

Old   August 10, 2010, 15:25
Exclamation ???
  #3
Senior Member
 
maysmech's Avatar
 
Join Date: Jan 2010
Posts: 347
Blog Entries: 2
Rep Power: 8
maysmech is on a distinguished road
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
maysmech is offline   Reply With Quote

Old   August 11, 2010, 04:48
Default
  #4
Senior Member
 
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 6
stevenvanharen is on a distinguished road
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?
stevenvanharen is offline   Reply With Quote

Old   August 11, 2010, 05:02
Default
  #5
Senior Member
 
maysmech's Avatar
 
Join Date: Jan 2010
Posts: 347
Blog Entries: 2
Rep Power: 8
maysmech is on a distinguished road
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]
maysmech is offline   Reply With Quote

Old   August 11, 2010, 08:57
Default
  #6
Senior Member
 
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 6
stevenvanharen is on a distinguished road
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!
stevenvanharen is offline   Reply With Quote

Old   August 12, 2010, 04:10
Default
  #7
Senior Member
 
maysmech's Avatar
 
Join Date: Jan 2010
Posts: 347
Blog Entries: 2
Rep Power: 8
maysmech is on a distinguished road
Quote:
Originally Posted by stevenvanharen View Post
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


maysmech is offline   Reply With Quote

Old   August 12, 2010, 16:12
Default
  #8
Senior Member
 
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 6
stevenvanharen is on a distinguished road
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).
stevenvanharen is offline   Reply With Quote

Old   August 12, 2010, 16:42
Default
  #9
Senior Member
 
maysmech's Avatar
 
Join Date: Jan 2010
Posts: 347
Blog Entries: 2
Rep Power: 8
maysmech is on a distinguished road
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.
maysmech is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Export CFD Solver results to text file format Maik Main CFD Forum 2 May 23, 2008 01:34
Could you comare StarCD with CFX 5?Help, please... Suteh CFX 54 November 7, 2001 21:12
CFX 5.5 Roued CFX 1 October 2, 2001 16:49
Setting a B.C using UserFortran in 4.3 tokai CFX 10 July 17, 2001 16:25
results of my MAC solver Maciej Matyka Main CFD Forum 4 January 5, 2001 18:38


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