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/)
-   -   the persistent printStack error (http://www.cfd-online.com/Forums/openfoam-solving/111407-persistent-printstack-error.html)

immortality January 6, 2013 11:59

the persistent printStack error
 
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.
Code:

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$

Code:

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:
Code:

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;
       
    }
}



All times are GMT -4. The time now is 08:44.