# why divergence occures in these cases

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

 December 28, 2012, 10:51 why divergence occures in these cases #1 Senior Member     Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,208 Rep Power: 26 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::T(scalar f, scalar T0, scalar (specieThermo::*F)(const scalar) const, scalar (specieThermo::*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 >::T(double, double, double (Foam::specieThermo >::*)(double) const, double (Foam::specieThermo >::*)(double) const, double (Foam::specieThermo >::*)(double) const) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so" #3 Foam::ePsiThermo > > > >::calculate() in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so" #4 Foam::ePsiThermo > > > >::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::T(scalar f, scalar T0, scalar (specieThermo::*F)(const scalar) const, scalar (specieThermo::*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 >::T(double, double, double (Foam::specieThermo >::*)(double) const, double (Foam::specieThermo >::*)(double) const, double (Foam::specieThermo >::*)(double) const) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so" #3 Foam::ePsiThermo > > > >::calculate() in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so" #4 Foam::ePsiThermo > > > >::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\$ 　 　```

 December 30, 2012, 02:27 #2 Senior Member     Ehsan Join Date: Oct 2012 Location: Iran Posts: 2,208 Rep Power: 26 does errors come from thermophyisical model?

 January 25, 2013, 10:21 #3 Member     Luis Felipe Gutierrez Marcantoni Join Date: Oct 2010 Location: Cordoba-Argentina Posts: 47 Rep Power: 15 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