CFD Online Discussion Forums - data.solverPerformance.tg at line 0.

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.

Code:

#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).

Code:

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;

Code:

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