[Kappa]

Hey, I am using interfoam and I got a Floating Point Exception but I am not able to find its source. I am using the k-omega-SST model for a multiphase flow.
The error message reads as follows:

[3] #0  Foam::error::printStack(Foam::Ostream&) in "/appl/openfoam/2.3.1/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #1  Foam::sigFpe::sigHandler(int) in "/appl/openfoam/2.3.1/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #2  
[3]  at sigaction.c:0
[3] #3  double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) in "/appl/openfoam/2.3.1/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #4  Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/appl/openfoam/2.3.1/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #5  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/appl/openfoam/2.3.1/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[3] #6  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "myFolder/myApplication"
[3] #7  Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) in "/appl/openfoam/2.3.1/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[3] #8  Foam::incompressible::RASModels::kOmegaSST::correct() in "/appl/openfoam/2.3.1/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[3] #9  
[3]  in "myFolder/myApplication"
[3] #10  __libc_start_main in "/lib64/libc.so.6"
[3] #11  
[3]  in "myFolder/myApplication"
[myPC] *** Process received signal ***
[myPC] Signal: Floating point exception (8)
[myPC] Signal code:  (-6)

Thanks everybody for help. I really do not know how to go on. The CFL number should be smaller than zero and the mesh quality should be good also. I cannot find where I may divide by zero :/

Btw, at the end it says:

--------------------------------------------------------------------------
mpirun noticed that process rank 3 with PID 14925 on node pbplnx10 exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

[Antimony]

Hi,

In your k and/or omega files, is there anything with a value 0? Division by zero will result in floating point error....

Cheers,
Antimony

[Kappa]

Hey,

thanks for your answer.
Neither in the initial data, nor in the timesteps I'm writing out...Omega is between 200-300 in average, k is between 2e-10 - 1e-3 in average.

So I don't think that they get zero, the minimal value is always larger than 0.

I just compared my kOmegaSST.H and ~.C file and I did not change anything in these. They are exactly the same as in version 2.3.1.

[Kappa]

Okay I was wrong.
the values for k and omega start diverging somewhen.
It says that 1001 iterations were used to solve for k but k ranges from 1e-64 to 1e+64.

So this is where the error comes from, but why did it need so many iterations without any effect?

[vasava]

Can you post your case setup?