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

Question on courant number

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 18, 2012, 12:14
Default Question on courant number
  #1
Senior Member
 
Daniele Vicario
Join Date: Mar 2009
Location: Novara, Italy
Posts: 142
Rep Power: 8
danvica is on a distinguished road
I'm running a pisofoam with k-OmegaSST turbolence model case with an average courant number of about 0.04 but a max value of 3.2.

The residual for all the variables are ok, so far the solver had calculated thousand steps without any problem, but...

How can I verify whether the results are trustable ?

Is courant number just an indication of the convergence or something more ?

Thanks for any comment.
Daniele
danvica is offline   Reply With Quote

Old   March 19, 2012, 03:48
Default
  #2
Senior Member
 
Suresh kumar Kannan
Join Date: Mar 2009
Location: Luxembourg, Luxembourg, Luxembourg
Posts: 129
Rep Power: 8
kumar is on a distinguished road
HI,
From my experience to check if your solution is right or not you should also check the residuals, especially the Final residuals.

I think courant number is mostly used to check the stability and to keep your solver with in the limits of the deltaT, to achieve a reliable and stable solution.

regards
K.Suresh kumar
kumar is offline   Reply With Quote

Old   March 19, 2012, 03:55
Default
  #3
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 919
Rep Power: 17
akidess will become famous soon enough
The Courant number doesn't even indicate convergence. All it says is that for the PISO algorithm to work, you'll probably get into troubles if your maximum Courant number is larger than 1.

Does it make sense for your case that the maximum and the average velocity are so two orders of magnitude apart? Do you have any (experimental) data you can use as a guide?

- Anton
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
*Help define the OpenFOAM stackexchange Q&A site: http://area51.stackexchange.com/prop...oam-technology
akidess is offline   Reply With Quote

Old   March 19, 2012, 07:19
Default
  #4
Senior Member
 
Daniele Vicario
Join Date: Mar 2009
Location: Novara, Italy
Posts: 142
Rep Power: 8
danvica is on a distinguished road
Thanks for the feedback.

I'm simulation the flow of water into a valve. As you (Anton) said, two order of magnitude is a lot.
So far, from Paraview, all I can say is that the flow velocity in the valve is nowhere more that three times the one at the inlet.

Basing on the definition of Courant number this means there's some "problem" with the mesh. But is it a problem (if the solver converges) ?

Daniele
danvica is offline   Reply With Quote

Old   March 19, 2012, 11:11
Default
  #5
Senior Member
 
kmooney's Avatar
 
Kyle Mooney
Join Date: Jul 2009
Location: Amherst, MA USA - San Diego, CA USA
Posts: 268
Rep Power: 9
kmooney is on a distinguished road
If your Co is getting that high I'm guessing that you don't have adjustable time stepping enabled and a maxCo defined. Look into the controlDict options available and try it out.
kmooney is offline   Reply With Quote

Old   March 19, 2012, 11:32
Default
  #6
Senior Member
 
Daniele Vicario
Join Date: Mar 2009
Location: Novara, Italy
Posts: 142
Rep Power: 8
danvica is on a distinguished road
No, your're right.
But I don't want to limit the timestep just because some cells are too small. I'm not worried about them... or should I ? That's the question.

In the meanwhile the solver is still parallel running, residuals are fine (convergence is reached within tollerance in 1-2 iteration), and mean/max Courant num are oscillating at about 0.04/3.0...

Is there any way to display Courant number in Parafoam ? As it was a physical field.

Daniele
danvica is offline   Reply With Quote

Old   March 19, 2012, 11:38
Default
  #7
Senior Member
 
kmooney's Avatar
 
Kyle Mooney
Join Date: Jul 2009
Location: Amherst, MA USA - San Diego, CA USA
Posts: 268
Rep Power: 9
kmooney is on a distinguished road
You could instantiate a new volScalarField, populate it with the local Courant No and have it print out with the rest of your results. Take a look at CourantNo.H for some hints on how to calculate the Co field.
kmooney 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
Courant Number Problems wschosta OpenFOAM Running, Solving & CFD 4 July 15, 2011 15:57
external flow with snappyHexMesh chelvistero OpenFOAM 11 January 15, 2010 20:43
Problems with Courant number (LaunderGibsonTurbulence Model) sven OpenFOAM 3 August 10, 2009 03:12
Courant number, patches, etc oort OpenFOAM 1 July 24, 2009 18:05
Fluent 6.4 courant number Aris Nikolopoulos FLUENT 0 May 6, 2008 08:52


All times are GMT -4. The time now is 01:38.