CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   why divergence occures in these cases (https://www.cfd-online.com/Forums/openfoam-solving/111135-why-divergence-occures-these-cases.html)

 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.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 December 30, 2012 03:27

does errors come from thermophyisical model?

 lfgmarc January 25, 2013 11:21

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

 All times are GMT -4. The time now is 10:15.