immortality 
December 28, 2012 11:51 
why divergence occures in these cases
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.12046e05, Final residual = 7.31173e17, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 4.25274e05, Final residual = 3.62861e17, No Iterations 3
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 4.75901e06, Final residual = 3.99706e15, No Iterations 3
smoothSolver: Solving for omega, Initial residual = 0.000196995, Final residual = 2.13054e13, No Iterations 4
smoothSolver: Solving for k, Initial residual = 0.000793458, Final residual = 6.24992e12, 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.77558e09
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.11839e05, Final residual = 7.19837e17, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 4.24677e05, Final residual = 3.78725e17, No Iterations 3
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 4.75475e06, Final residual = 3.99285e15, No Iterations 3
smoothSolver: Solving for omega, Initial residual = 0.000195555, Final residual = 2.08341e13, No Iterations 4
smoothSolver: Solving for k, Initial residual = 0.000784559, Final residual = 6.17617e12, 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.77558e09
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.11625e05, Final residual = 7.20281e17, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 4.24064e05, Final residual = 3.99302e17, No Iterations 3
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 4.75018e06, Final residual = 3.99205e15, No Iterations 3
smoothSolver: Solving for omega, Initial residual = 0.000194523, Final residual = 2.07203e13, No Iterations 4
smoothSolver: Solving for k, Initial residual = 0.000776509, Final residual = 6.20535e12, 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.77558e09
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.11423e05, Final residual = 7.28519e17, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 4.23456e05, Final residual = 3.49586e17, No Iterations 3
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 4.74569e06, Final residual = 3.99724e15, 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/OpenFOAM2.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_64linuxgnu/libc.so.6"
#7
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
thesis@thesisX58AUD7:~/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.66411e17, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.000169952, Final residual = 3.73178e17, No Iterations 3
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 7.07211e05, Final residual = 5.33129e15, No Iterations 3
smoothSolver: Solving for omega, Initial residual = 0.000158878, Final residual = 7.99203e14, No Iterations 4
smoothSolver: Solving for k, Initial residual = 0.000316136, Final residual = 6.48863e13, No Iterations 4
ExecutionTime = 7423.7 s ClockTime = 7443 s
Mean and max Courant Numbers = 0.00916482 0.049979
deltaT = 4.26603e09
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.72206e17, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.000170605, Final residual = 4.01179e17, No Iterations 3
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 7.09006e05, Final residual = 5.40426e15, No Iterations 3
smoothSolver: Solving for omega, Initial residual = 0.000159318, Final residual = 8.0518e14, No Iterations 4
smoothSolver: Solving for k, Initial residual = 0.000316903, Final residual = 6.50375e13, No Iterations 4
ExecutionTime = 7423.79 s ClockTime = 7444 s
Mean and max Courant Numbers = 0.00916483 0.0499776
deltaT = 4.26603e09
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.2299e17, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.000171294, Final residual = 3.69854e17, No Iterations 3
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 7.10944e05, Final residual = 5.32607e15, 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/OpenFOAM2.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_64linuxgnu/libc.so.6"
#7
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
thesis@thesisX58AUD7:~/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@thesisX58AUD7:~/Desktop/method_4_2_2(revised).042$
