[keishawillstone]

Hi there!
I'm dealing with a porous media modeling,after solving the pressure I have the velocity which is proportional to thi latter according to this:
U = -C*fvc::grad(p);
But when it comes to the equation that account for the generation(within the incinerator) of each species like O2,CO2..it seems like it crashes,here is the part of it that crashes,
solve
(
fvm::div(phirhoesp,yCO) == SCO
);
about the solver I'm using:
in the fvSchemes file
div(phirhoesp,yCO) Gauss upwind;
(I also tried "div(phirhoesp,yCO) Gauss limitedLinear 1")
in the fvSolution file
yCO PBiCG
{
preconditioner DILU;
tolerance 1e-05;
relTol 0;
};
I even tried this:
yCO smoothSolver
{
smoother GaussSeidel;
tolerance 1e-05;
relTol 0;
};
Then it sends this error message!!!!
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xb8052400]
#3 Foam:ILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam:ILUPreconditioner:ILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::Istream&) in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#5 Foam::lduMatrix:reconditioner::addasymMatrixCons tructorToTable<Foam:ILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::Istream&) in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#6 Foam::lduMatrix:reconditioner::New(Foam::lduMatr ix::solver const&, Foam::Istream&) in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#8 Foam::fvMatrix<double>::solve(Foam::Istream&) in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libfiniteVolume.so"
#9 Foam::fvMatrix<double>::solve() at OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fvMatrixSolve.C:162
#10 Foam::lduMatrix::solverPerformance Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/fvMatrix.C:1301
#11 main at OpenFOAM/energy-1.4.1/applications/solvers/ironwill/rhohEqn.H:35
#12 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#13 _start in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/ironwill"
Floating point exception
Can somebody tell me how to handle this?please,I've already tried so many things but it didn't work,I'm quite stuck!!!!
Best Regards
Thanks in advance