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

transient case diverged all of a sudden after already convergence

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 11, 2012, 10:39
Default transient case diverged all of a sudden after already convergence
  #1
Member
 
Cong
Join Date: Nov 2012
Posts: 66
Rep Power: 3
wc34071209 is on a distinguished road
Hi everyone,

I am running an internal flow in a duct. Firstly I ran the steady simulation (3500 iterations ) using simpleFoam and everything is OK. Then I started running transient simulation using transient SIMPLE method. In the first 0.9958 s, everything is also OK. However, when the 0.99585 time step starts, the pressure calculation gets suddenly divergence. I tried to modify convergence tolerance control and other linear system solvers and ran from a previous time step, but the divergence exactly recurred at 0.99585 s.

Although the pressure field diverged, the velocity field not. It is really strange.

The log is as follows.

I hope someone could give me a hint. Thank you everyone in advance.

Courant Number mean: 0.0002247080246 max: 7.268116999
Time = 3500.99585
Max Inner Iteration: 40
ConvergenceCriterion: 9e-08
-- Inner loop: 0 ---------------------------------
smoothSolver: Solving for Ux, Initial residual = 3.722149168e-06, Final residual = 4.894892661e-10, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 1.756301862e-05, Final residual = 1.469450516e-09, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 1.01839746e-05, Final residual = 1.132826601e-09, No Iterations 1
GAMG: Solving for p, Initial residual = 4.942934403e-06, Final residual = 4.204448761e-08, No Iterations 7
time step continuity errors : sum local = 1.081550112e-13, global = 9.481328711e-15, cumulative = -2.379358346e-11
smoothSolver: Solving for nuTilda, Initial residual = 9.707065056e-06, Final residual = 2.616585371e-09, No Iterations 1
-> maxResidual(U,p) 4.942934403e-06
-- Inner loop: 1 ---------------------------------
smoothSolver: Solving for Ux, Initial residual = 1.117159902e-06, Final residual = 1.819146731e-10, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 5.270313626e-06, Final residual = 4.919923245e-10, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 3.05667921e-06, Final residual = 4.055659424e-10, No Iterations 1
GAMG: Solving for p, Initial residual = 4.934697876e-05, Final residual = 1.730930522e-07, No Iterations 4
time step continuity errors : sum local = 4.452616788e-13, global = -1.255549445e-13, cumulative = -2.39191384e-11
smoothSolver: Solving for nuTilda, Initial residual = 2.915621506e-06, Final residual = 9.895764204e-10, No Iterations 1
-> maxResidual(U,p) 4.934697876e-05
-- Inner loop: 2 ---------------------------------
smoothSolver: Solving for Ux, Initial residual = 3.502918942e-07, Final residual = 7.495885771e-11, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 1.675695577e-06, Final residual = 1.964754382e-10, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 9.391553624e-07, Final residual = 1.568032312e-10, No Iterations 1
GAMG: Solving for p, Initial residual = 0.3252932271, Final residual = 0.000827214952, No Iterations 3
time step continuity errors : sum local = 2.924368657e-09, global = -1.922521341e-10, cumulative = -2.161712725e-10
smoothSolver: Solving for nuTilda, Initial residual = 1.225521515e-06, Final residual = 4.562680108e-10, No Iterations 1
-> maxResidual(U,p) 0.3252932271
wc34071209 is offline   Reply With Quote

Old   November 11, 2012, 17:31
Default
  #2
Senior Member
 
sail's Avatar
 
Vieri Abolaffio
Join Date: Jul 2010
Location: Always on the move.
Posts: 308
Rep Power: 7
sail is on a distinguished road
your courant number is too large, try using a timestep equal to 1/10 of what you are actually using
__________________
http://www.leadingedge.it/
Naval architecture and CFD consultancy
sail is offline   Reply With Quote

Old   November 12, 2012, 06:55
Default
  #3
Member
 
Cong
Join Date: Nov 2012
Posts: 66
Rep Power: 3
wc34071209 is on a distinguished road
Hi,

Thank you for your reply.

I know that CFL < 1 is a common recommendation for transient calculation, but I am also told that for URANS, CFL between 1 and 10 can be also acceptable. Is it true ?

And even if CFL is too big, but why for the first 0.9958 s does it work 'well' ?
wc34071209 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
Case comparison of steady state and time-averaged transient solutions k.vafiadis CFX 2 October 20, 2012 05:37
convergence of pitzdailyExptinlet case openfoambeginner OpenFOAM Running, Solving & CFD 0 July 23, 2012 04:33
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 12:24
Fluidisation transient case Vikas Kumar Main CFD Forum 0 April 24, 2003 03:29
ICCG convergence case hong Main CFD Forum 7 December 30, 2000 05:44


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