![]() |
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 |
Hi,Stefan
I have the same problem with you,have you solved it? |
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 |
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 |
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) |
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 02:19. |