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/)
-   -   janafThermo Error on coalChemistryFoam (http://www.cfd-online.com/Forums/openfoam-solving/85520-janafthermo-error-coalchemistryfoam.html)

CST February 28, 2011 09:48

janafThermo Error on coalChemistryFoam
 
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

yingkun June 13, 2011 23:37

Hi,Stefan
I have the same problem with you,have you solved it?

CST June 14, 2011 03:00

Hey you,

sorry, i do not have an answer to that problem above - but for me worked another case:

i did just turn off the radiation (my theme of work was able to be changed to no radiation) - then the error didn't occur another time.

Good luck,

CST

gregor June 14, 2011 03:53

Hi,

you should provide some more informations about your case because the reason could be anywhere. What are your boundary conditions ? Which models are you using ...

gregor

yingkun June 14, 2011 05:02

hi
the model is just as the tutorial case simplifiedSiwek, I just change the mesh and boundary conditions and relational files in constant, my case has two inlets(inlet1 jets coal and primary air at 300K,inlet2 jets preheated air at 1500K)

the BC like follows:
http://cg5a39.mail.163.com/app/wp/fi...ction=showfile
(提取码:9sf3ezow ~~sorry,I don't know another way to upload my BC file)

lfgmarc June 14, 2011 21:46

I have the same problem with XiFOAM
 
Hi, foamers, I have the same problem with XiFoam solver, I used the pitzdialy case how starting point, I changed the geometry according to :

https://picasaweb.google.com/1009028...3/XiFoam_case#

The conditions for my case are:

http://ifile.it/j1mgu7h/conditions.tar.gz

If any one can help me, I was very thankful

best regards

Felipe


All times are GMT -4. The time now is 17:50.