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

sudden overflow of the courant number

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 3, 2013, 08:54
Default sudden overflow of the courant number
  #1
Senior Member
 
Join Date: Jan 2013
Posts: 134
Rep Power: 13
kkpal is on a distinguished road
dear everyone
In my simulation with pimpleFoam, I firstly set the MaxCo at 0.5, it went smoothly. To slightly save the computation time I increased the MaxCo to 0.85, then bad things happened.
The first few hundreds iterations after the MaxCo change were fine, with the courant number varying a little below and above 0.85. But at a precise moment, a sudden increase in the courant number(from 0.85 to 3) occured, here is what it looks like.
Code:
Courant Number mean: 0.0554198 max: 0.834122
deltaT = 0.0055491
Time = 7.717

PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 0.000262906, Final residual = 3.06265e-07, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00444772, Final residual = 5.03058e-06, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.000323524, Final residual = 3.89366e-07, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0589746, Final residual = 0.000576175, No Iterations 28
time step continuity errors : sum local = 1.20116e-09, global = 4.29182e-10, cumulative = 3.79531e-10
PIMPLE: iteration 2
DILUPBiCG:  Solving for Ux, Initial residual = 6.00503e-07, Final residual = 6.00503e-07, No Iterations 0
DILUPBiCG:  Solving for Uy, Initial residual = 1.21582e-05, Final residual = 1.39211e-07, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 5.54895e-07, Final residual = 5.54895e-07, No Iterations 0
GAMG:  Solving for p, Initial residual = 0.189246, Final residual = 0.00174058, No Iterations 30
time step continuity errors : sum local = 4.25715e-09, global = 1.74154e-09, cumulative = 2.12107e-09
PIMPLE: iteration 3
DILUPBiCG:  Solving for Ux, Initial residual = 1.23274e-06, Final residual = 1.80222e-08, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 4.09755e-05, Final residual = 1.73558e-08, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 1.13289e-07, Final residual = 1.79303e-09, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.589279, Final residual = 9.8303e-10, No Iterations 412
time step continuity errors : sum local = 4.74361e-15, global = 1.39072e-15, cumulative = 2.12108e-09
ExecutionTime = 162003 s  ClockTime = 165697 s

forceCoeffs output:
    Cm    = -0.0148137
    Cd    = 1.20736
    Cl    = 0.31442
    Cl(f) = 0.142396
    Cl(r) = 0.172023

Calculating averages

Courant Number mean: 0.055421 max: 3.13703
deltaT = 0.00149738
Time = 7.71849

PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 7.40663e-05, Final residual = 1.49712e-07, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00127949, Final residual = 2.6649e-06, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 8.75115e-05, Final residual = 5.88664e-08, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.92039, Final residual = 0.00902539, No Iterations 25
time step continuity errors : sum local = 2.01598e-08, global = 5.18475e-09, cumulative = 7.30583e-09
PIMPLE: iteration 2
DILUPBiCG:  Solving for Ux, Initial residual = 8.33085e-06, Final residual = 8.33085e-06, No Iterations 0
DILUPBiCG:  Solving for Uy, Initial residual = 0.000220168, Final residual = 2.06308e-06, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 5.25244e-07, Final residual = 5.25244e-07, No Iterations 0
GAMG:  Solving for p, Initial residual = 0.950336, Final residual = 0.00888994, No Iterations 26
time step continuity errors : sum local = 4.3503e-08, global = 1.43381e-08, cumulative = 2.1644e-08
PIMPLE: iteration 3
DILUPBiCG:  Solving for Ux, Initial residual = 1.29058e-05, Final residual = 1.3311e-08, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.000512759, Final residual = 1.47315e-08, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 7.73346e-07, Final residual = 8.57267e-09, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.980802, Final residual = 9.85576e-10, No Iterations 546
time step continuity errors : sum local = 1.26467e-14, global = 3.30246e-15, cumulative = 2.1644e-08
ExecutionTime = 162067 s  ClockTime = 165761 s

forceCoeffs output:
    Cm    = -1.12254
    Cd    = 6.38589
    Cl    = 3.02313
    Cl(f) = 0.389023
    Cl(r) = 2.6341
This sudden change seemed to take place from nowhere
anyone knows why this is happening?
kkpal is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Courant number fireman FLUENT 7 September 11, 2021 12:33
Sudden jump in Courant number NJG OpenFOAM Running, Solving & CFD 7 May 15, 2014 14:52
Stable boundaries marcoymarc CFX 33 March 13, 2013 07:39
LES near wall model & courant number kasim CFX 5 March 16, 2008 19:23
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 03:58


All times are GMT -4. The time now is 12:05.