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

OpenFOAM v.2.x Courant Number 4 time higher then in previous versions

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

Like Tree2Likes
  • 1 Post By makaveli_lcf
  • 1 Post By makaveli_lcf

Reply
 
LinkBack Thread Tools Display Modes
Old   February 16, 2012, 05:57
Default OpenFOAM v.2.x Courant Number 4 time higher then in previous versions
  #1
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
Hi Foamers!

I got and issue using OpenFOAM v2.x:

Courant Number jumps 4 times if compare to OpenFOAM v1.6 and 1.7.x.
I tried cases for different geometries ofcouse with the same setting between OF versions and result is always the same!
Actually this issue results in a calculation time increasing 4 times. Unappropriated!

Please, comment on that issue, what is your experience with v2.x usage?

Cheer,
Alex
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   February 16, 2012, 08:54
Default
  #2
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
Reported at Mantis:

http://www.openfoam.com/mantisbt/view.php?id=424
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   February 16, 2012, 10:05
Default
  #3
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
Issue closed (http://www.openfoam.com/mantisbt/view.php?id=424#c1039):

Quote:
The method for calculating the Courant number has been changed since version 1.7.x to be better represent the stability requirement of collocated meshes. In general the change hase resulted in slightly lower values.

A change in the method to calculate the Courant number does not in itself change the stability of the algorithm or lead to divergence as it is not used in the algorithm; only to set the time-step. If you have a case which behaved differently for the same fixed time-step between different versions of OpenFOAM the issue does not relate to a change in the way in which the Courant number is calculated.
So then it means that OF 2.x is less stable with the same fixed time step...
Will do more checks and report regarding actual calculation time.
chegdan likes this.
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   February 17, 2012, 12:02
Default
  #4
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Quote:
Originally Posted by makaveli_lcf View Post
Issue closed (http://www.openfoam.com/mantisbt/view.php?id=424#c1039):



So then it means that OF 2.x is less stable with the same fixed time step...
Will do more checks and report regarding actual calculation time.
What solver? I was aware of the change in the definition of Co, which was made to remove inconsistency at boundaries, but I haven't experienced major issues with 2.1.x.
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   February 17, 2012, 16:20
Default
  #5
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
Hi Alberto!

pimpleFoam. According to my test to keep the same Co number (with which solution does not diverge) actual calculation time increased 2 times (I remind with the same constant dt I've got Co_new_OF = 4*Co_old_OF).
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   February 17, 2012, 16:22
Default
  #6
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
Slow down, to my opinion, is not 4 times because now turbulence is solved at last PIMPLE iteration. But of course I tried with turbulence updated each iteration: the same result, solution diverges with the dt I used before with old OF versions..
anmartin likes this.
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   February 17, 2012, 22:14
Default
  #7
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Quote:
Originally Posted by makaveli_lcf View Post
Slow down, to my opinion, is not 4 times because now turbulence is solved at last PIMPLE iteration. But of course I tried with turbulence updated each iteration: the same result, solution diverges with the dt I used before with old OF versions..
OK, let's try to do some troubleshooting, if you like

Do you use any particular boundary condition? For example mapped, or time-varying, or custom BC? Some of these were updated in 2.x.
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   February 18, 2012, 07:55
Default
  #8
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
Nope, standard k-epsilon (also tried with realizable which I mostly use) with standard wall functions for k, epsilon, nut, fixed values for U at inlet and zero gradient at outlet, pressure 0 at outlet and zero gradient at the inlet. I think it would be better if I upload here some test case. Will do it on Monday.
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   February 18, 2012, 22:28
Default
  #9
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Yes, that's a good idea. The BC's are surely not the reason.
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Reply

Tags
courant number, openfoam

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
Very high courant number, time step continuity errors and bounding epsilon erncyc OpenFOAM Running, Solving & CFD 1 March 24, 2011 12:00
Error while running rhoPisoFoam.. nileshjrane OpenFOAM Running, Solving & CFD 8 August 26, 2010 12:50
calculation diverge after continue to run zhajingjing OpenFOAM 0 April 28, 2010 04:35
Differences between serial and parallel runs carsten OpenFOAM Bugs 11 September 12, 2008 11:16
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07


All times are GMT -4. The time now is 08:15.