CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   why divergence occures in these cases (http://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 13:55.