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

Hi

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

Thanks
Anant

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?

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

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.