[Thamali]

Dear Foamers,

I am solving a case using an Euler-Euler solver developed for combustion of wood chips.
When I run the case using the solver,it stops at calculating "rCH4" and FPE error appears.

#0  Foam::error::printStack(Foam::Ostream&) at ~/OpenFOAM/OpenFOAM-2.2.2/src/OSspecific/POSIX/printStack.C:221
#1  Foam::sigFpe::sigHandler(int) at ~/OpenFOAM/OpenFOAM-2.2.2/src/OSspecific/POSIX/signals/sigFpe.C:117
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3   in "/lib/x86_64-linux-gnu/libm.so.6"
#4  Foam::pow(double, double) at ~/OpenFOAM/OpenFOAM-2.2.2/src/OpenFOAM/lnInclude/doubleFloat.H:78
#5  Foam::pow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ~/OpenFOAM/OpenFOAM-2.2.2/src/OpenFOAM/fields/Fields/scalarField/scalarField.C:118 (discriminator 2)
#6  void Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<double> const&) at ~/OpenFOAM/OpenFOAM-2.2.2/src/OpenFOAM/lnInclude/GeometricScalarField.C:274
#7  
 at ~/OpenFOAM/OpenFOAM-2.2.2/src/OpenFOAM/lnInclude/GeometricScalarField.C:328
#8  
 at ~/OpenFOAM/OpenFOAM-2.2.2/src/OpenFOAM/lnInclude/GeometricScalarField.C:350 (discriminator 1)
#9  
 at ~/OpenFOAM/thamali-2.2.2/applications/solvers/latest/testonlymoistandcharnpyrovolat/UEqn1.H:95 (discriminator 1)
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11  
 at ??:?
Floating point exception (core dumped)

But,when I check the parameters YCH4,YO2 or tg of previous time step necessary for calculation,all are positive values.
When I,unset FOAM_SIGFPE,
The "rCH4" become negative at 2 points(1 point in a boundary and 1point in internalField).

wrong token type - expected Scalar, found on line 0 the word 'nan'

file: /home/thamali/OpenFOAM/thamali-2.2.2/run/tutorials/combustion/woodchipcombustion/new/latest/withwallradiation/test/drypyrochacovolat/system/data.solverPerformance.tg at line 0.

    From function operator>>(Istream&, Scalar&)
    in file lnInclude/Scalar.C at line 91.

FOAM exiting

I checked case log and for tg;

smoothSolver:  Solving for tg, Initial residual = nan, Final residual = nan, No Iterations 1000

I tried by changing "fvSchemes " and they did not work.
I am attaching the files calculating "rCH4" and file calculating Gasphase mas fractions(Y) and tg.

I highly appreciate your help.
If you want to look at other files of the solver,I am happy put upload.

Thanks in advance.
Regards,
Thamali