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/)
-   -   Problem in solving (http://www.cfd-online.com/Forums/openfoam-solving/110068-problem-solving.html)

immortality December 3, 2012 09:19

Problem in solving
 
this error tells the maximum number of iterations the solver is reached probably in e equation solving.why this is occuring?
Im engaged in this problem but it don't want to be resolved!



fluxScheme: Kurganov

Starting time loop

Mean and max Courant Numbers = 0.10019 0.211345
Time = 1e-09

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
swak4Foam: Allocating new repository for sampledGlobalVariables
smoothSolver: Solving for Ux, Initial residual = 0.878819, Final residual = 9.93747e-10, No Iterations 438
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 9.87651e-10, No Iterations 693
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 0.94115, Final residual = 5.58403e-06, No Iterations 1002


--> 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.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 69.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::specieThermo<Foam::hConstThermo<Foam::perfec tGas> >::T(double, double, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfe ctGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfe ctGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfe ctGas> >::*)(double) const) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#3 Foam::ePsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::calculate() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#4 Foam::ePsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::correct() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#5
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
Aborted

diwakaranant January 29, 2013 02:10

Hi

I am also having the same problem.
Did you find the cause of its occurrence ?

Thanks
Anant

immortality January 29, 2013 10:51

hi dear Anant.I did some changes in fvSchemes and fvSolutions.using CrankNicholson .5;and div(phi,U) Gauss upwind; is better than Euler and linear respectively.
a very important thing is boundary conditions and initial conditions.I changed initial condition and increased it because it had made a supersonic flow at inlet but it must be subsonic in my problem.
and also try another solver and change maxIter for e in specieThermo.c (I think it was there) can lead to more stability.
whats the problem you want to model?comp or incomp?

diwakaranant January 30, 2013 01:45

Dear Ehsan

Thanks a lot for your help. I tweaked the initial conditions and it's
working now.

I am solving compressible flow between rotating cylinders using rhoCentralFoam.

Regards
Anant

immortality January 30, 2013 05:57

Dear Anant
Im working with rhoCentralFoam and rhoPimpleFoam.
Im interested in the problem you are working on.could you send me your case please?
Thanks.


All times are GMT -4. The time now is 19:14.