[immortality]

Hi.the error persists.I included the thermophysical file.why this error ccures?how can I fix it?
any idea is valuable to me.
thanks.

diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 1.89761e-07, Final residual = 4.15927e-17, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 1.36181e-07, Final residual = 3.78343e-17, No Iterations 3
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for e, Initial residual = 5.14607e-08, Final residual = 4.03473e-15, No Iterations 3
ExecutionTime = 4239.85 s  ClockTime = 4250 s

Mean and max Courant Numbers = 0.00877655 0.0499934
deltaT = 3.73726e-09
Time = 0.000429391

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.89779e-07, Final residual = 4.36025e-17, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 1.36267e-07, Final residual = 4.10536e-17, No Iterations 3
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for e, Initial residual = 5.15161e-08, Final residual = 3.99024e-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$

thermoType      ePsiThermo<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>;

mixture
{
    specie
    {
        nMoles          1;
        molWeight       28.96;
    }
    thermodynamics
    {
        Cp              1004.5;
        Hf              2.544e+06;
    }
    transport
    {
        mu              0;
        Pr              1;
    }
}

The T BC is as below:

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 288;

boundaryField
{
   right
    {
        type zeroGradient;
    }

    left
    {
       type groovyBC;

        variables (

      //"pi=3.1415926535;"
                  "ymax=max(pos().y);"
                  "ymin=min(pos().y);"
                  "r=0.06;"
                  "rpm=7200;"
                  "omega=rpm*pi/30;"
                  "v_r=r*omega;"
                  "w_cell=.003;"
                  "w_w1=.003;"
                  "w_p1=.06*11*pi/180;"
                  "w_w2=.06*5*pi/180;"
                  "w_p2=.06*74*pi/180;"
                  "w_w3=.003;"
                  "p0_1=250000;"
                  "T0_1=288;"
                  "p0_2=100000;"
                  "T0_2=288;"
                  "gamma=1.4;"
                  "R=287.14;"
                  
);

     
         
       fractionExpression "((w_w1-pos().y)/v_r<time() && time()<(w_w1+w_p1-pos().y)/v_r)?1:0";
       valueExpression "T0_1-(gamma-1)/(2*gamma*R)*magSqr(internalField(U))";
       value uniform 288;
       gradientExpression "0";
       
    }

    walls
    {
        type zeroGradient;
        
    }

    empty
    {
        type empty;
        
    }
}