[Heno]

Hi foamers,
In my program, a control equation only includes a transient term and a source term like that:

ddt(rho,Y) == Sourece

When I run the solver ,the error message shows that:

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam:perator/(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#5 Foam::diagonalSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at /opt/openfoam4/src/finiteVolume/lnInclude/fvMatrixSolve.C:82
#8 Foam::fvMatrix<double>::solve() at /opt/openfoam4/src/finiteVolume/lnInclude/fvMatrixSolve.C:325 (discriminator 6)
#9 ? at ~/OpenFOAM/tutorials/pemfccoldstartModel-4.0.0/pemfccoldstartSolver/applications/../src/myFoam/myFoam12.H:12 (discriminator 42)
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11 ? at ??:?

Please tell me how can I do to solve the problem, thanks a lot.

[Tobermory]

Firstly, you have given us almost no information on what your coding is, so it's not very easy to help you out. Have a read of Guide: How to ask a question on the forums and then have a think about how you might phrase your post to give readers the best chance of helping you.

In the mean time, one thing that we can do is help you understand the error message - it is complaining about a floating point error (#1) as a result of a divide operation (#3, #4) in the diagonal solver (#5). Perhaps check that rho is properly defined and non-zero?

[Heno]

Originally Posted by Tobermory

Firstly, you have given us almost no information on what your coding is, so it's not very easy to help you out. Have a read of Guide: How to ask a question on the forums and then have a think about how you might phrase your post to give readers the best chance of helping you.

In the mean time, one thing that we can do is help you understand the error message - it is complaining about a floating point error (#1) as a result of a divide operation (#3, #4) in the diagonal solver (#5). Perhaps check that rho is properly defined and non-zero?

Thanks for your help Tobermory! I will read the how-ask-question-forums carefully.

I cannot provide my code for some reasons, and I sure the rho is non-zero. But I added a laplacian term in my control equation and my code ran(The value of the laplacian term is relatively small). I tried to delete the laplacian term in another equation and kept the transit term and source term only. It also didn't work and apear the same error message. So I gueess the problem is in the matrix solver. I am not sure if the matrix solver in OpenFOAM can solve the equations with only transient and source terms.