|
[Sponsors] |
December 28, 2012, 11:51 |
why divergence occures in these cases
|
#1 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27 |
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$ |
|
December 30, 2012, 03:27 |
|
#2 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27 |
does errors come from thermophyisical model?
|
|
January 25, 2013, 11:21 |
|
#3 |
Member
|
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
__________________
Felipe G |
|
|
|
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 10:32 |
Divergence problem for species transport model | MY | FLUENT | 3 | January 11, 2014 05:46 |
Quarter Burner mesh with periosic condition | SamCanuck | FLUENT | 2 | August 31, 2011 12:34 |
How to speed up the simulation of a large amount of similar steady state cases | tss | Main CFD Forum | 4 | November 5, 2010 15:36 |
Directory constructions of complex cases | 7islands | OpenFOAM | 2 | August 26, 2008 06:25 |