[CST]

Dear Forum,

I'm testing the coalChemistryFoam-Solver on OpenFOAM-1.7. There i got a small case, i.e. a tube with only one Particle of coal, i want to track devolating and combusting. Therefore I'd like to simulate coldflow first. But everytime is start the calculation, i get this errormessage:

==========================
Solving chemistry
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.9999999999, Final residual = 2.357901633e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 2.475990053e-07, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 8.867723085e-06, No Iterations 1
DILUPBiCG: Solving for O2, Initial residual = 1, Final residual = 9.395203123e-08, No Iterations 2
DILUPBiCG: Solving for CH4, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for hs, Initial residual = 1, Final residual = 3.923674927e-09, No Iterations 2

--> FOAM FATAL ERROR:
attempt to use janafThermo<equationOfState> out of temperature range 0 -> 5000; T = -38253.3
From function janafThermo<equationOfState>::checkT(const scalar T) const
in file /home/nini/OpenFOAM/OpenFOAM-1.7.0/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 63.
FOAM aborting
#0 Foam::error: :PrintStack(Foam::Ostream&) in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam::specieThermo<Foam::janafThermo<Foam: :PerfectGas> >::Hs(double) const in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#3 Foam::specieThermo<Foam::janafThermo<Foam: :PerfectGas> >::THs(double, double) const in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#4 Foam::hsPsiMixtureThermo<Foam::reactingMixture<Foa m::sutherlandTransport<Foam::specieThermo<Foam::ja nafThermo<Foam: :PerfectGas> > > > >::calculate() in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#5 Foam::hsPsiMixtureThermo<Foam::reactingMixture<Foa m::sutherlandTransport<Foam::specieThermo<Foam::ja nafThermo<Foam: :PerfectGas> > > > >::correct() in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/libreactionThermophysicalModels.so"
#6 main in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/applications/bin/linux64GccDPOpt/coalChemistryFoam"
#7 __libc_start_main in "/lib64/libc.so.6"
#8 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/openfoam/OpenFOAM/OpenFOAM-1.7.0/applications/bin/linux64GccDPOpt/coalChemistryFoam"
Abgebrochen
===============================


So I'm "out of range". My Temperature seems to be not within the borders of 0 -> 5000K. But my T-file is set-up with an internal field of 1300K and zeroGradient at walls/patches.


So where could be my error? Setting the calculation/fv-Schemes to Gauss linear uncorrected or limited doesn't seem to gain some success.


The code seems to read the "scalar T" (I guess from my T-file?) and compares with the borders - but where can be the fault?


With best regards,
STefan