immortality |
January 6, 2013 10: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;
}
}
|