floating point exception running pimpleFoam

elisa_france · July 3, 2023, 12:20

Hello everyone,
I'm working on the CFD simulation of a turbine. I use a rotational mesh obtained with snappyHexMesh. The checkMesh is ok and the orhtogonality of the cells does not exceed 60.

I use the k-epsilon model with standard parameters and the pimpleFoam solver.
PIMPLE
{
// correctPhi yes;
nOuterCorrectors 2;
nCorrectors 1;
nNonOrthogonalCorrectors 2;

}

My ICs and BCs are the following :

epsilon :

Code:

internalField   uniform 0.000773;
boundaryField
{
    #includeEtc "caseDicts/setConstraintTypes"

    inlet
    {
        type            fixedValue;
        value           $internalField;
    }

    outlet
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }


    OuterCylinder
    {
        type            epsilonWallFunction;
        value           $internalField;
    }
    rotor
    {
        type            epsilonWallFunction;
        value           $internalField;
    }
}

k:

Code:

InternalField   uniform  3.75e-3; 
boundaryField
{
    #includeEtc "caseDicts/setConstraintTypes"

    inlet
    {
        type            fixedValue;
        value           $internalField;
    }
    outlet
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }
    rotor
    {
        type            kqRWallFunction;
        value           $internalField;
    }
    OuterCylinder
    {
        type            kqRWallFunction;
        value           $internalField;
    }
}

nut:

Code:

internalField   uniform 0;

boundaryField
{
    #includeEtc "caseDicts/setConstraintTypes"

    inlet
    {
        type             calculated; 
        value           $internalField;
    }
    outlet
    {
        type             calculated; 
        value           $internalField;
    }
    rotor
    {
        type            nutkWallFunction;
        value           $internalField;
    }
    OuterCylinder
    {
        type            nutkWallFunction;
        value           $internalField;
    }
}

p:

Code:

internalField   uniform 0;

boundaryField
{
    #includeEtc "caseDicts/setConstraintTypes"

    inlet
    {
        type            zeroGradient;
    }

    outlet
    {
        type            fixedValue;
        value           uniform 0;
    }

   wall
    {
        type            zeroGradient;
    }
}

U:

Code:

internalField   uniform (0 0 0);

boundaryField
{
    #includeEtc "caseDicts/setConstraintTypes"

    inlet
    {
        type            fixedValue;
        value           uniform (5 0 0);
    }

    outlet
    {
        type            inletOutlet;
        inletValue      uniform (0 0 0); 
        value           $internalField;
    }

    OuterCylinder
    {
        type            noSlip;
    }

    rotor
    {
        type            movingWallVelocity;
        value           uniform (0 0 0);
    }
}

I tried to change multiple times the schemes of discretization but I obtain the same error when I run my simulation in parallel (with the ultimate iteration):

Code:

PIMPLE: Iteration 1
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.00518303, No Iterations 2
GAMG:  Solving for pcorr, Initial residual = 0.110429, Final residual = 0.00686751, No Iterations 1
GAMG:  Solving for pcorr, Initial residual = 0.0412413, Final residual = 0.00282683, No Iterations 1
time step continuity errors : sum local = 6.68413e-12, global = -3.13249e-12, cumulative = -3.81649e-08
smoothSolver:  Solving for Ux, Initial residual = 0.000167296, Final residual = 1.61408e-09, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.00142363, Final residual = 5.53925e-08, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.00142653, Final residual = 3.14146e-07, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.529876, Final residual = 0.00152502, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.0134142, Final residual = 9.01409e-05, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00346595, Final residual = 3.82134e-07, No Iterations 10
time step continuity errors : sum local = 6.23079e-13, global = -2.3497e-13, cumulative = -3.81651e-08
PIMPLE: Iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.000114029, Final residual = 3.50674e-09, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.000993846, Final residual = 1.42656e-07, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.000997294, Final residual = 7.1382e-07, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.769851, Final residual = 0.00236002, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.111126, Final residual = 0.000683643, No Iterations 4
GAMG:  Solving for p, Initial residual = 0.0215559, Final residual = 5.31638e-07, No Iterations 14
time step continuity errors : sum local = 1.15072e-13, global = 3.08406e-14, cumulative = -3.81651e-08
smoothSolver:  Solving for epsilon, Initial residual = 9.75282e-07, Final residual = 9.75282e-07, No Iterations 0
smoothSolver:  Solving for k, Initial residual = 4.45736e-06, Final residual = 8.17555e-07, No Iterations 2
bounding k, min: -88.1444 max: 255.752 average: 0.30806
ExecutionTime = 161.49 s  ClockTime = 161 s

forces forces write:
    sum of forces:
        pressure : (4.1351 0.000648792 0.00196503)
        viscous  : (0.0152723 0.000181542 -0.000325192)
        porous   : (0 0 0)
    sum of moments:
        pressure : (-0.144806 0.000382245 0.000819488)
        viscous  : (0.00404579 2.59582e-05 2.95074e-05)
        porous   : (0 0 0)

Courant Number mean: 0.000618101 max: 1.0822
deltaT = 1.53306e-05
Time = 0.00178537

PIMPLE: Iteration 1
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.00487637, No Iterations 2
GAMG:  Solving for pcorr, Initial residual = 0.109329, Final residual = 0.0069315, No Iterations 1
GAMG:  Solving for pcorr, Initial residual = 0.0414208, Final residual = 0.00311275, No Iterations 1
time step continuity errors : sum local = 6.03569e-12, global = -3.01038e-12, cumulative = -3.81681e-08
smoothSolver:  Solving for Ux, Initial residual = 0.000791879, Final residual = 7.50408e-05, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.0333141, Final residual = 0.00294954, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.244911, Final residual = 0.018558, No Iterations 7
[1] #0  Foam::error::printStack(Foam::Ostream&)[2] #0  Foam::error::printStack(Foam::Ostream&)[3] #0  Foam::error::printStack(Foam::Ostream&)[5] #0  Foam::error::printStack(Foam::Ostream&)[6] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[0] #1  Foam::sigFpe::sigHandler(int) at ??:?
 at ??:?
 at ??:?
[7] #1  Foam::sigFpe::sigHandler(int) at ??:?
[4] #1  [5] #1  Foam::sigFpe::sigHandler(int)Foam::sigFpe::sigHandler(int)[1] #1  Foam::sigFpe::sigHandler(int) at ??:?
[6] #1   at ??:?
Foam::sigFpe::sigHandler(int)[2] #1  Foam::sigFpe::sigHandler(int) at ??:?
[0] #2  ? at ??:?
[3] #2  ? at ??:?
[5] #2  ? at ??:?
[2] #2  ? in "/lib/x86_64-linux-gnu/ at ??:?
[6] #2  libc.so.6"
[0] #3  ?? in "/lib/x86_64-linux-gnu/libc.so.6"
 in [5] #3  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3  ?"/lib/x86_64-linux-gnu/libc.so.6"
[4] #3  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[7] #3  ? in "/lib/x86_64-linux-gnu/libmpi.so.40"
[0] #4  ompi_coll_base_allreduce_intra_recursivedoubling in "/lib/x86_64-linux-gnu/libc.so.6"
[6] #3  ? in "/lib/x86_64-linux-gnu/libmpi.so.40"
[4] #4  ompi_coll_base_allreduce_intra_recursivedoubling in "/lib/x86_64-linux-gnu/libmpi.so.40"
[1] #4  ompi_coll_base_allreduce_intra_recursivedoubling in "/lib/x86_64-linux-gnu/libmpi.so.40"
[0] #5  PMPI_Allreduce[2] #4  ompi_coll_base_allreduce_intra_recursivedoubling in "/lib/x86_64-linux-gnu/libmpi.so.40"
[7] #4  ompi_coll_base_allreduce_intra_recursivedoubling in "/lib/x86_64-linux-gnu/libmpi.so.40"
[6] #4  ompi_coll_base_allreduce_intra_recursivedoubling in "/lib/x86_64-linux-gnu/libmpi.so.40"
[1] #5  PMPI_Allreduce in "/lib/x86_64-linux-gnu/libmpi.so.40"
[4] #5  PMPI_Allreduce in "/lib/x86_64-linux-gnu/libmpi.so.40"
[3] #5  PMPI_Allreduce in "/lib/x86_64-linux-gnu/libmpi.so.40"
[2] #5  PMPI_Allreduce in "/lib/x86_64-linux-gnu/libmpi.so.40"
[5] #5  PMPI_Allreduce in "/lib/x86_64-linux-gnu/libmpi.so.40"
[6] #5  PMPI_Allreduce in "/lib/x86_64-linux-gnu/libmpi.so.40"
[7] #5  PMPI_Allreduce in "/lib/x86_64-linux-gnu/libmpi.so.40"
[0] #6  void Foam::allReduce<double, Foam::sumOp<double> >(double&, int, ompi_datatype_t*, ompi_op_t*, Foam::sumOp<double> const&, int, int) in "/lib/x86_64-linux-gnu/libmpi.so.40"
[3] #6   in "/lib/x86_64-linux-gnu/libmpi.so.40"
[4] #6  void Foam::allReduce<double, Foam::sumOp<double> >(double&, int, ompi_datatype_t*, ompi_op_t*, Foam::sumOp<double> const&, int, int)void Foam::allReduce<double, Foam::sumOp<double> >(double&, int, ompi_datatype_t*, ompi_op_t*, Foam::sumOp<double> const&, int, int) in "/lib/x86_64-linux-gnu/libmpi.so.40"
[1] #6  void Foam::allReduce<double, Foam::sumOp<double> >(double&, int, ompi_datatype_t*, ompi_op_t*, Foam::sumOp<double> const&, int, int) in "/lib/x86_64-linux-gnu/libmpi.so.40"
[7] #6  void Foam::allReduce<double, Foam::sumOp<double> >(double&, int, ompi_datatype_t*, ompi_op_t*, Foam::sumOp<double> const&, int, int) in "/lib/x86_64-linux-gnu/libmpi.so.40"
[2] #6  void Foam::allReduce<double, Foam::sumOp<double> >(double&, int, ompi_datatype_t*, ompi_op_t*, Foam::sumOp<double> const&, int, int) in "/lib/x86_64-linux-gnu/libmpi.so.40"
[5] #6  void Foam::allReduce<double, Foam::sumOp<double> >(double&, int, ompi_datatype_t*, ompi_op_t*, Foam::sumOp<double> const&, int, int) in "/lib/x86_64-linux-gnu/libmpi.so.40"
[6] #6  void Foam::allReduce<double, Foam::sumOp<double> >(double&, int, ompi_datatype_t*, ompi_op_t*, Foam::sumOp<double> const&, int, int) at ??:?
[0] #7  Foam::reduce(double&, Foam::sumOp<double> const&, int, int) at ??:?
[4] #7  Foam::reduce(double&, Foam::sumOp<double> const&, int, int) at ??:?
[3] #7  Foam::reduce(double&, Foam::sumOp<double> const&, int, int) at ??:?
[1] #7  Foam::reduce(double&, Foam::sumOp<double> const&, int, int) at ??:?
[5] #7  Foam::reduce(double&, Foam::sumOp<double> const&, int, int) at ??:?
[6] #7  Foam::reduce(double&, Foam::sumOp<double> const&, int, int) at ??:??
[0] #8  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[7] #7  Foam::reduce(double&, Foam::sumOp<double> const&, int, int) at ??:?
[4] #8   at ??:?
Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const[3] #8  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[1] #8  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
 at ??:?
[6] #8  [2] #8  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) constFoam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[7] #8  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[0] #9  Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const at ??:?
 at ??:?
[3] #9[4] #9    Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) constFoam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const at ??:?
[2] #9  Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const at ??:?
 at ??:?
[6] #9  Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const[1] #9  Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const at ??:?
[5] #9  Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const at ??:?
[7] #9  Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const at ??:?
 at ??:?
[4] #10  [0] #10  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const[3] #10  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&,
Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) constFoam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
[2] #[1] #10  10  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, unsigned char) constFoam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const[6] #10  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
[5] #10  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
[7] #10  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
 at ??:?
 at ??:?
[4] #11  [3] #11  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const[0] #11  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) constFoam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
 at ??:?
[6] #11  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const[2] #11  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[1] #11  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[5] #11  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[7] #11  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const[0] #12  [4] #12  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&)Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&)[3] #12  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
 at ??:?
[6] #12  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&)[2] #12  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[1] #12  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[5] #12  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[7] #12  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[4] #13   at ??:?
[0] #13  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
[3] #13  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
[2] #13   at ??:?
[6] #13  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
[1] #13  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[4] #14  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[0] #14   at ??:?
[7] #13  Foam::fvMatrix<double>::solve(Foam::dictionary const&)Foam::fvMatrix<double>::solve() at ??:?
[5] #13   in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[2] #14   in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[3] #14  Foam::fvMatrix<double>::solve() in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[6] #14  Foam::fvMatrix<double>::solve(Foam::dictionary const&)Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[1] # in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[4] #15  14  Foam::fvMatrix<double>::solve()Foam::fvMatrix<double>::solve() in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[0] #Foam::fvMatrix<double>::solve()15  Foam::fvMatrix<double>::solve()? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[7] #14   in  in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[2] #15  "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[3] #15   in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[5] #14  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[6] #15  Foam::fvMatrix<double>::solve() in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[1] #15   in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[4] #16  __libc_start_mainFoam::fvMatrix<double>::solve()?? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[0] #16  ?__libc_start_main? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[7] #15   in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[5] #15   in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[2] #16  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[4] #17   in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #17  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
 in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[6] #16  [3] #16  __libc_start_main__libc_start_main? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[1] #16  __libc_start_main? in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #17  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[5] #16  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #17   in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"

I don't know what to do any more, I've changed the initial conditions, the boundary conditions, the mesh, the discretisation schemes and nothing helps. Can anyone give me some advice?

Alczem · July 4, 2023, 07:11

Hey

I am running a somewhat similar simulation. Did you make sure your mesh is properly split up at the interface between the moving zone and the static zone? You can check with moveDynamicMesh for a few timesteps and see if the mesh is deformed with Paraview. If the cells are warped, it means that the cells are not split up. You can use mergeOrSplitBaffles -split to split the mesh at the faceZone used to create the sliding mesh (I assume that you used a cellzone and facezone definition in your snappyHexMeshDict?).

Can you provide the meshing setup and your dynamicMeshDict to be sure?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Floating point exception (PLEASE HELP :) )	mn14cat	OpenFOAM CC Toolkits for Fluid-Structure Interaction	9	December 8, 2021 12:24
InterFoam based solver running into floating point error on restarting simulation	Venky_94	OpenFOAM Running, Solving & CFD	9	November 23, 2021 17:53
ERORR floating point exception fluent udf	cudau.95	Fluent UDF and Scheme Programming	0	August 4, 2021 23:11
[blockMesh] Floating Point Exception while generating wedge based mesh	jns-v	OpenFOAM Meshing & Mesh Conversion	9	July 8, 2021 06:36
Floating point exception for pimpleFoam	gharek	OpenFOAM Running, Solving & CFD	8	May 10, 2021 09:55

July 4, 2023, 07:11		#2
Alczem Senior Member Join Date: Dec 2021 Posts: 251 Rep Power: 5	Hey I am running a somewhat similar simulation. Did you make sure your mesh is properly split up at the interface between the moving zone and the static zone? You can check with moveDynamicMesh for a few timesteps and see if the mesh is deformed with Paraview. If the cells are warped, it means that the cells are not split up. You can use mergeOrSplitBaffles -split to split the mesh at the faceZone used to create the sliding mesh (I assume that you used a cellzone and facezone definition in your snappyHexMeshDict?). Can you provide the meshing setup and your dynamicMeshDict to be sure?