[cyln]

Hello, I am having error due to division by zero. I checked all my BCs. Everything looks fine. Any idea about what the problem is?

PIMPLE: iteration 1
--> FOAM Warning : 
    From function void Foam::lduMatrix::operator-=(const Foam::lduMatrix&)
    in file matrices/lduMatrix/lduMatrix/lduMatrixOperations.C at line 286
    Unknown matrix type combination
    this : diagonal:0 symmetric:0 asymmetric:1
    A    : diagonal:0 symmetric:0 asymmetric:0
   Normalisation factor = 87582.7
smoothSolver:  Iteration 0 residual = 1
smoothSolver:  Iteration 4 residual = 8.65356e-09
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 8.65356e-09, No Iterations 4
   Normalisation factor = 0.0721236
smoothSolver:  Iteration 0 residual = 0.743557
smoothSolver:  Iteration 4 residual = 5.10218e-05
smoothSolver:  Solving for Uy, Initial residual = 0.743557, Final residual = 5.10218e-05, No Iterations 4
   Normalisation factor = 0.0705167
smoothSolver:  Iteration 0 residual = 0.802723
smoothSolver:  Iteration 4 residual = 4.89182e-05
smoothSolver:  Solving for Uz, Initial residual = 0.802723, Final residual = 4.89182e-05, No Iterations 4

#0  Foam::error::printStack(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>&, double const&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/openfoam4/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
#5  ? in "/opt/openfoam4/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  ? in "/opt/openfoam4/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
Floating point exception (core dumped)

[kindle]

As indicated, there should be a "divide by zero" for fields algebra.

You can try the full debug build to get the exact function "Foam::divide" taking arg (Field<T>, T, UList<T>). Or try turn on debugSwitch to help locate the function.

[cyln]

Thank you kindle,

When I type "unset FOAM_SIGFPE", I get rid of this error.

However, I cannot manage to converge my simulation. My momentum equations are not converging. May this have an effect on converge?

[kindle]

Well, you don't get rid of the problem by unsetting/ignoring the Error. The Info you get above means that there is a division by zero happening and you have to deal with that.