Error after running pisoFoam
 

Hi everybody,
I'm a new comer and new user of OpenFoam within just one month. Currently I'm working on simple simulation using a solid block sets on the middle of the domain. After running the simulation I get below errors and I couldn't understand the meaning. Could anybody help me on this matter?Your cooperation is highly appreciated. Thanks in advance.


#0 Foam::error::printStack(Foam::Ostream&) in "/opt/usrx/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/usrx/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 __restore_rt at sigaction.c:0
#3 void Foam::fvc::surfaceIntegrate<double>(Foam::Field<do uble>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/usrx/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/pisoFoam"
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::surfaceIntegrate<double>(Foam::Geometri cField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/usrx/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/pisoFoam"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div<double>(Foam::GeometricField<double , Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/usrx/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/pisoFoam"
#6 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/usrx/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7 Foam::fv::laplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/usrx/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/opt/usrx/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/pisoFoam"
#9 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/usrx/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/pisoFoam"
#10 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/usrx/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#11 Foam::incompressible::RASModels::kEpsilon::correct () in "/opt/usrx/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#12 main in "/opt/usrx/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/pisoFoam"
#13 __libc_start_main in "/lib64/libc.so.6"
#14 Foam::UOPstream::write(char) in "/opt/usrx/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/pisoFoam"
Floating point exception

Hello,

do you have a division by zero or for any of your variables very high values?

Another thing could be wrong boundary conditions, did you check them?

These two things came in my mind first for a "floating point exception".


regards
treima

This kind of error can be several things and is indicative of a divide by zero. if you initialize your k or epsilon field with 0 or have a 0 boundary condition...then you will get this. To avoid this, try to change your zero entries to 1e-20 and provide a good initial guess for k and epsilon fields.

Dear Treima,
Thank for your prompt reply.
1)I dont have any division zero and variables with high values
2)For BC, I set U and p as below;
U;
inlet
{ type fixedValue;
value (12 0 0);
}
outlet
{ type zeroGradient;}
for surrounding walls I set as fixedValue, uniform (0 0 0)

p;
inlet
{type zeroGradient;}
outlet
{type fixedValue;
value uniform 0;
}
for surrounding walls I set as zeroGradient;
I have calculated the values of k and epsilon.
Need your advice on this matter.

While I don´t know your case exactly I can imagine that your boundary condition for p at outlet, p = 0, might cause some problems.

Which solver do you use? Did you have a look in the tutorials for some example cases and set up your case in a similar form?

The boundary conditions for U and p seem ok.

The question is not whether there is one or not, the question is really where it is and how it can be found. Foam::sigFpe indicates a floating point exception which almost always indicates a division by zero.

When does the error occur exactly? Is it before the start of the iterations, after a few iterations, or at the end of the full simulation? Or at an other point in time?

Regards,

L

Dear Lieven,
Thanks for your explanation.
Actually this occurred after a few iterations.