CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   how to handle this crash (

keishawillstone July 9, 2009 08:49

how to handle this crash
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, seems like it crashes,here is the part of it that crashes,
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
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::printStack(Foam::Ostream&) in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/"
#2 Uninterpreted: [0xb8052400]
#3 Foam::DILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/"
#4 Foam::DILUPreconditioner::DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::Istream&) in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/"
#5 Foam::lduMatrix::preconditioner::addasymMatrixCons tructorToTable<Foam::DILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::Istream&) in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/"
#6 Foam::lduMatrix::preconditioner::New(Foam::lduMatr ix::solver const&, Foam::Istream&) in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/"
#7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/"
#8 Foam::fvMatrix<double>::solve(Foam::Istream&) in "/home/energy/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/"
#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/"
#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

All times are GMT -4. The time now is 04:10.