Skip to content

Sloshing2D tutorial crashing in openfoam2012

Summary

The sloshing2D tutorial case is crashing with a floating point exception.

Steps to reproduce

Run the sloshing2D tutorial case.

Example case

See the sloshing2D case included in the openfoam2012 release, located at tutorials/incompressible/pimpleFoam/laminar/sloshing2D

What is the current bug behaviour?

Run crashes with floating point exception. When run in the 2006 version the tutorial runs with no issue.

What is the expected correct behavior?

Tutorial is expected to run with no issue.

Relevant logs and/or images

tail of output of pimpleFoam

Courant Number mean: 0.0219597 max: 1.05238
Time = 0.48

PIMPLE: iteration 1
Maximal capillary Courant number: 4.01058
Free surface curvature: min = -0.103997, max = 0.0975977, average = -3.99992e-05
Free surface continuity error : sum local = 0.0197539, global = 1.75478e-11
DICPCG:  Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG:  Solving for cellMotionUy, Initial residual = 0.457312, Final residual = 5.64373e-09, No Iterations 48
smoothSolver:  Solving for Ux, Initial residual = 0.0427553, Final residual = 1.82027e-10, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 0.0483463, Final residual = 8.64354e-11, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.208065, Final residual = 8.26681e-09, No Iterations 21
time step continuity errors : sum local = 8.53615e-12, global = -3.40791e-12, cumulative = -1.51755e-10
GAMG:  Solving for p, Initial residual = 0.0254321, Final residual = 7.7193e-09, No Iterations 16
time step continuity errors : sum local = 9.91592e-12, global = 4.01099e-12, cumulative = -1.47744e-10
GAMG:  Solving for p, Initial residual = 0.00370436, Final residual = 4.97571e-09, No Iterations 16
time step continuity errors : sum local = 6.30442e-12, global = -2.5717e-12, cumulative = -1.50316e-10
PIMPLE: iteration 2
Free surface continuity error : sum local = 0.0314803, global = 7.36567e-10
DICPCG:  Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG:  Solving for cellMotionUy, Initial residual = 0.74463, Final residual = 9.73061e-09, No Iterations 49
smoothSolver:  Solving for Ux, Initial residual = 0.017863, Final residual = 1.36831e-10, No Iterations 4
smoothSolver:  Solving for Uy, Initial residual = 0.0240514, Final residual = 2.23125e-10, No Iterations 4
GAMG:  Solving for p, Initial residual = 0.129237, Final residual = 9.01299e-09, No Iterations 21
time step continuity errors : sum local = 9.515e-12, global = 3.51329e-12, cumulative = -1.46802e-10
GAMG:  Solving for p, Initial residual = 0.0444841, Final residual = 9.74374e-09, No Iterations 18
time step continuity errors : sum local = 7.47935e-12, global = -2.77065e-12, cumulative = -1.49573e-10
GAMG:  Solving for p, Initial residual = 0.00602138, Final residual = 7.04264e-09, No Iterations 16
time step continuity errors : sum local = 5.44304e-12, global = 2.03037e-12, cumulative = -1.47543e-10
PIMPLE: iteration 3
Free surface continuity error : sum local = 0.0307808, global = 7.89252e-10
DICPCG:  Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG:  Solving for cellMotionUy, Initial residual = 0.827931, Final residual = 6.51922e-09, No Iterations 49
smoothSolver:  Solving for Ux, Initial residual = 0.00918507, Final residual = 1.91392e-10, No Iterations 4
smoothSolver:  Solving for Uy, Initial residual = 0.011106, Final residual = 1.70839e-10, No Iterations 4
GAMG:  Solving for p, Initial residual = 0.188214, Final residual = 5.6468e-09, No Iterations 21
time step continuity errors : sum local = 5.19288e-12, global = 3.17129e-12, cumulative = -1.44371e-10
GAMG:  Solving for p, Initial residual = 0.0291086, Final residual = 8.36559e-09, No Iterations 16
time step continuity errors : sum local = 8.81108e-12, global = -5.34691e-12, cumulative = -1.49718e-10
GAMG:  Solving for p, Initial residual = 0.00488945, Final residual = 7.43674e-09, No Iterations 15
time step continuity errors : sum local = 7.90801e-12, global = 4.81474e-12, cumulative = -1.44903e-10
PIMPLE: iteration 4
Free surface continuity error : sum local = 0.00477103, global = 4.35966e-11
DICPCG:  Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG:  Solving for cellMotionUy, Initial residual = 0.0896289, Final residual = 6.78209e-09, No Iterations 47
smoothSolver:  Solving for Ux, Initial residual = 0.00788484, Final residual = 2.09289e-12, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 0.00482004, Final residual = 1.00621e-12, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.149903, Final residual = 7.07337e-09, No Iterations 18
time step continuity errors : sum local = 1.05058e-11, global = 6.34139e-12, cumulative = -1.38562e-10
GAMG:  Solving for p, Initial residual = 0.027082, Final residual = 8.70494e-09, No Iterations 16
time step continuity errors : sum local = 1.45241e-11, global = 8.70469e-12, cumulative = -1.29857e-10
GAMG:  Solving for p, Initial residual = 0.00467499, Final residual = 8.92156e-09, No Iterations 13
time step continuity errors : sum local = 1.51271e-11, global = -9.2286e-12, cumulative = -1.39086e-10
PIMPLE: iteration 5
Free surface continuity error : sum local = 0.0287436, global = -8.29794e-10
DICPCG:  Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG:  Solving for cellMotionUy, Initial residual = 0.551951, Final residual = 7.17263e-09, No Iterations 51
smoothSolver:  Solving for Ux, Initial residual = 0.0113926, Final residual = 4.32554e-10, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.0139625, Final residual = 1.39168e-10, No Iterations 7
GAMG:  Solving for p, Initial residual = 0.0598591, Final residual = 4.9069e-09, No Iterations 20
time step continuity errors : sum local = 7.71824e-12, global = -4.4745e-12, cumulative = -1.4356e-10
GAMG:  Solving for p, Initial residual = 0.0327082, Final residual = 5.9765e-09, No Iterations 16
time step continuity errors : sum local = 9.03766e-12, global = 5.26532e-12, cumulative = -1.38295e-10
GAMG:  Solving for p, Initial residual = 0.0268987, Final residual = 8.39761e-09, No Iterations 15
time step continuity errors : sum local = 1.28788e-11, global = -7.44382e-12, cumulative = -1.45739e-10
PIMPLE: iteration 6
Free surface continuity error : sum local = 0.0346845, global = -7.9425e-10
DICPCG:  Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG:  Solving for cellMotionUy, Initial residual = 0.809593, Final residual = 7.7121e-09, No Iterations 52
smoothSolver:  Solving for Ux, Initial residual = 0.0133566, Final residual = 8.67802e-10, No Iterations 17
smoothSolver:  Solving for Uy, Initial residual = 0.0189276, Final residual = 5.5088e-10, No Iterations 17
GAMG:  Solving for p, Initial residual = 0.315144, Final residual = 5.3234e-09, No Iterations 6
time step continuity errors : sum local = 15973, global = -27168.6, cumulative = -27168.6
GAMG:  Solving for p, Initial residual = 0.00485725, Final residual = 6.75073e-09, No Iterations 17
time step continuity errors : sum local = 42937.3, global = 26776.7, cumulative = -391.905
GAMG:  Solving for p, Initial residual = 2.02096e-15, Final residual = 2.02096e-15, No Iterations 0
time step continuity errors : sum local = 4.27168e+10, global = -1.42809e+10, cumulative = -1.42809e+10
PIMPLE: iteration 7
Free surface continuity error : sum local = 8.3222e+10, global = -7.14045e+10
DICPCG:  Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG:  Solving for cellMotionUy, Initial residual = 1, Final residual = 7.14591e-09, No Iterations 52
smoothSolver:  Solving for Ux, Initial residual = 0.70222, Final residual = 2.92166e-10, No Iterations 15
smoothSolver:  Solving for Uy, Initial residual = 0.604349, Final residual = 5.74197e-10, No Iterations 13
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib/x86_64-linux-gnu/libpthread.so.0
#3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4  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 ??:?
#5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?
#8  Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
#9  ? in ~/OpenFOAM/OpenFOAM-v2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam
#10  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#11  ? in ~/OpenFOAM/OpenFOAM-v2012/platforms/linux64GccDPInt32Opt/bin/pimpleFoam
Floating point exception (core dumped)

Environment information

  • OpenFOAM version :v2012
  • Operating system :Ubuntu
  • Hardware info :Intel(R) Core(TM) i5-8265U CPU @ 1.60GHz
  • Compiler :gcc

Possible fixes