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

why divergence occures in these cases

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By lfgmarc

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 28, 2012, 10:51
Default why divergence occures in these cases
  #1
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26
immortality is on a distinguished road
the runs were performing for 2 days but eventually stoped with these errors.
I use constant thermophysical properties as the case of shock wave.what may be the cause?
Code:
Time = 0.002893253
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 1.12046e-05, Final residual = 7.31173e-17, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 4.25274e-05, Final residual = 3.62861e-17, No Iterations 3
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 4.75901e-06, Final residual = 3.99706e-15, No Iterations 3
smoothSolver: Solving for omega, Initial residual = 0.000196995, Final residual = 2.13054e-13, No Iterations 4
smoothSolver: Solving for k, Initial residual = 0.000793458, Final residual = 6.24992e-12, No Iterations 4
bounding k, min: -34.3897 max: 99716.7 average: 95.1909
ExecutionTime = 9947.03 s ClockTime = 9976 s
Mean and max Courant Numbers = 0.00838603 0.0499914
deltaT = 3.77558e-09
Time = 0.002893257
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 1.11839e-05, Final residual = 7.19837e-17, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 4.24677e-05, Final residual = 3.78725e-17, No Iterations 3
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 4.75475e-06, Final residual = 3.99285e-15, No Iterations 3
smoothSolver: Solving for omega, Initial residual = 0.000195555, Final residual = 2.08341e-13, No Iterations 4
smoothSolver: Solving for k, Initial residual = 0.000784559, Final residual = 6.17617e-12, No Iterations 4
bounding k, min: -298.288 max: 99714.7 average: 95.0889
ExecutionTime = 9947.31 s ClockTime = 9976 s
Mean and max Courant Numbers = 0.00838603 0.0499914
deltaT = 3.77558e-09
Time = 0.002893261
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 1.11625e-05, Final residual = 7.20281e-17, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 4.24064e-05, Final residual = 3.99302e-17, No Iterations 3
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 4.75018e-06, Final residual = 3.99205e-15, No Iterations 3
smoothSolver: Solving for omega, Initial residual = 0.000194523, Final residual = 2.07203e-13, No Iterations 4
smoothSolver: Solving for k, Initial residual = 0.000776509, Final residual = 6.20535e-12, No Iterations 4
bounding k, min: -72.8202 max: 99710.3 average: 95.0768
ExecutionTime = 9947.6 s ClockTime = 9976 s
Mean and max Courant Numbers = 0.00838603 0.0499915
deltaT = 3.77558e-09
Time = 0.002893264
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 1.11423e-05, Final residual = 7.28519e-17, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 4.23456e-05, Final residual = 3.49586e-17, No Iterations 3
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 4.74569e-06, Final residual = 3.99724e-15, No Iterations 3
 
--> FOAM FATAL ERROR: 
Maximum number of iterations exceeded
From function specieThermo<Thermo>::T(scalar f, scalar T0, scalar (specieThermo<Thermo>::*F)(const scalar) const, scalar (specieThermo<Thermo>::*dFdT)(const scalar) const) const
in file /home/opencfd/OpenFOAM/OpenFOAM-2.1.0/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 69.
FOAM aborting
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::T(double, double, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#3 Foam::ePsiThermo<Foam::pureMixture<Foam::constTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::calculate() in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#4 Foam::ePsiThermo<Foam::pureMixture<Foam::constTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::correct() in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#5 
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7 
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
thesis@thesis-X58A-UD7:~/Desktop/method_4_2_2(revised)$ 
 
----------------------------------------------------------------------
and at this:
----------------------------------------------------------------------
Code:
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 0.000123841, Final residual = 5.66411e-17, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.000169952, Final residual = 3.73178e-17, No Iterations 3
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 7.07211e-05, Final residual = 5.33129e-15, No Iterations 3
smoothSolver: Solving for omega, Initial residual = 0.000158878, Final residual = 7.99203e-14, No Iterations 4
smoothSolver: Solving for k, Initial residual = 0.000316136, Final residual = 6.48863e-13, No Iterations 4
ExecutionTime = 7423.7 s ClockTime = 7443 s
Mean and max Courant Numbers = 0.00916482 0.049979
deltaT = 4.26603e-09
Time = 0.000514313
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 0.000123838, Final residual = 5.72206e-17, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.000170605, Final residual = 4.01179e-17, No Iterations 3
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 7.09006e-05, Final residual = 5.40426e-15, No Iterations 3
smoothSolver: Solving for omega, Initial residual = 0.000159318, Final residual = 8.0518e-14, No Iterations 4
smoothSolver: Solving for k, Initial residual = 0.000316903, Final residual = 6.50375e-13, No Iterations 4
ExecutionTime = 7423.79 s ClockTime = 7444 s
Mean and max Courant Numbers = 0.00916483 0.0499776
deltaT = 4.26603e-09
Time = 0.000514318
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 0.000123835, Final residual = 6.2299e-17, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.000171294, Final residual = 3.69854e-17, No Iterations 3
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 7.10944e-05, Final residual = 5.32607e-15, No Iterations 3
 
--> FOAM FATAL ERROR: 
Maximum number of iterations exceeded
From function specieThermo<Thermo>::T(scalar f, scalar T0, scalar (specieThermo<Thermo>::*F)(const scalar) const, scalar (specieThermo<Thermo>::*dFdT)(const scalar) const) const
in file /home/opencfd/OpenFOAM/OpenFOAM-2.1.0/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 69.
FOAM aborting
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::T(double, double, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#3 Foam::ePsiThermo<Foam::pureMixture<Foam::constTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::calculate() in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#4 Foam::ePsiThermo<Foam::pureMixture<Foam::constTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::correct() in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#5 
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7 
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
thesis@thesis-X58A-UD7:~/Desktop/method_4_2_2(revised)-.042$ paraFoam 
created temporary 'method_4_2_2(revised)-.042.OpenFOAM'
/opt/openfoam210/bin/paraFoam: 1: Syntax error: "(" unexpected
thesis@thesis-X58A-UD7:~/Desktop/method_4_2_2(revised)-.042$ 
 
 
immortality is offline   Reply With Quote

Old   December 30, 2012, 02:27
Default
  #2
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26
immortality is on a distinguished road
does errors come from thermophyisical model?
immortality is offline   Reply With Quote

Old   January 25, 2013, 10:21
Default
  #3
Member
 
lfgmarc's Avatar
 
Luis Felipe Gutierrez Marcantoni
Join Date: Oct 2010
Location: Cordoba-Argentina
Posts: 47
Rep Power: 15
lfgmarc is on a distinguished road
Send a message via MSN to lfgmarc
Hi, this error comes from the iterative method used to convert the energetic variable to Temperature (see OpenFOAM-1.6-ext/src/thermophysicalModels/specie/thermo/specieThermo/specieThermoI.H).


Try to use a small time step .
Regards
immortality likes this.
__________________
Felipe G
lfgmarc is offline   Reply With Quote

Reply


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
Fluidized Bed: Error: Divergence detected in AMG solver: pressure correction Error Ob Mole89 FLUENT 5 April 12, 2014 09:32
Divergence problem for species transport model MY FLUENT 3 January 11, 2014 04:46
Quarter Burner mesh with periosic condition SamCanuck FLUENT 2 August 31, 2011 11:34
How to speed up the simulation of a large amount of similar steady state cases tss Main CFD Forum 4 November 5, 2010 14:36
Directory constructions of complex cases 7islands OpenFOAM 2 August 26, 2008 05:25


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