[Sorabh]

I am using pimpleDyMFoam to simulate a stirrer. The Courant number is within the specified limits and the residuals seem to be okay too. However the following error is obtained.

Courant Number mean: 0.104544 max: 0.825854
deltaT = 5.1058e-05
Time = 0.000247919

AMI: Creating addressing and weights between 22980 source faces and 23116 target faces
AMI: Patch source sum(weights) min/max/average = 0.618253, 1.38838, 1.00215
AMI: Patch target sum(weights) min/max/average = 0, 1.12629, 0.994905
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 0.0188878, Final residual = 0.000519529, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.0161195, Final residual = 0.0004383, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.0185892, Final residual = 0.000509428, No Iterations 1
#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 "/home/mslab02/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/bin/pimpleDyMFoam"
#5  ? in "/home/mslab02/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/bin/pimpleDyMFoam"
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  ? in "/home/mslab02/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/bin/pimpleDyMFoam"
Floating point exception (core dumped)

I used the checkMesh command to see if there were any errors, and these are the 2 errors I found:

 ***Cells with small determinant (< 0.001) found, number of cells: 56
  <<Writing 56 under-determined cells to set underdeterminedCells
 ***Concave cells (using face planes) found, number of cells: 2773
  <<Writing 2773 concave cells to set concaveCells
    Face interpolation weight : minimum: 0.14992 average: 0.478769
    Face interpolation weight check OK.
    Face volume ratio : minimum: 0.0380045 average: 0.908863
    Face volume ratio check OK.
Calculating AMI weights between owner patch: AMI1 and neighbour patch: AMI2
AMI: Creating addressing and weights between 22980 source faces and 23116 target faces
AMI: Patch source sum(weights) min/max/average = 0.602327, 1.38482, 1.00254
AMI: Patch target sum(weights) min/max/average = 0.216937, 1.06614, 0.995235

Failed 2 mesh checks.

End

Could someone please help me out, as to how should I deal with this error

[Sorabh]

No one? I could gladly use some help..

[JM27]

Hi Sorabh,

I am also encountering a similar floating point exception error and I suspect this is due to a failed mesh check. Apparently there are some concave cells in my mesh. Strangely, I did not encounter this error before using a similar mesh.

Just wondering whether you had managed to solve this issue and, if so, how you did this?

I am also using OpenFOAM 5.x.

Any help will be appreciated. Thank you.

[louisgag]

Hello Saurabh,

AMI: Patch target sum(weights) min/max/average = 0, 1.12629, 0.994905

the 0 means that a cell from your AMI interface cannot find any neighbor. I suggest you use paraview to see what goes wrong with your rotation and if need be run something like

mergeOrSplitBaffles -split -overwrite

before your simulation.