[sErik]

Hi Foamers,

I'm on simulating a pipe (a quite big mesh, ~4mio cells) flow with pisoFoam and I have some troubles there. The simulation starts fine, but after a certain interations, it just crashes with an floating point error. The Co number seems to be fine. I first tried it with the PCiCG solver but the calculation took an eternity, so I switched to the GAMG solver.

This is the protocol:

...

Calculating averages

Time = 0.00014

Courant Number mean: 0.00608327 max: 1.29077
GAMG: Solving for Ux, Initial residual = 0.16536, Final residual = 1.13479e-06, No Iterations 1
GAMG: Solving for Uy, Initial residual = 0.0842274, Final residual = 1.66113e-06, No Iterations 1
GAMG: Solving for Uz, Initial residual = 0.100704, Final residual = 2.04478e-06, No Iterations 1
GAMG: Solving for p, Initial residual = 0.260381, Final residual = 0.00640364, No Iterations 3
time step continuity errors : sum local = 2.71062e-05, global = -6.4173e-07, cumulative = 2.24495e-05
GAMG: Solving for p, Initial residual = 0.115879, Final residual = 0.00266843, No Iterations 4
time step continuity errors : sum local = 1.75079e-05, global = 1.64369e-06, cumulative = 2.40932e-05
GAMG: Solving for k, Initial residual = 0.0670612, Final residual = 0.00116239, No Iterations 1
bounding k, min: -6.5128 max: 140.59 average: 0.00130018
ExecutionTime = 553.07 s ClockTime = 555 s

Calculating averages

Time = 0.00015

Courant Number mean: 0.00590866 max: 2.04424
GAMG: Solving for Ux, Initial residual = 0.132126, Final residual = 8.87388e-07, No Iterations 1
GAMG: Solving for Uy, Initial residual = 0.065539, Final residual = 1.00676e-06, No Iterations 1
GAMG: Solving for Uz, Initial residual = 0.0772455, Final residual = 1.38296e-06, No Iterations 1
GAMG: Solving for p, Initial residual = 0.544934, Final residual = 0.0172603, No Iterations 2
time step continuity errors : sum local = 9.3184e-05, global = 2.65063e-06, cumulative = 2.67438e-05
GAMG: Solving for p, Initial residual = 0.299918, Final residual = 0.00487727, No Iterations 2
time step continuity errors : sum local = 2.85995e-05, global = 2.48463e-06, cumulative = 2.92284e-05
GAMG: Solving for k, Initial residual = 0.0600043, Final residual = 0.000370476, No Iterations 1
bounding k, min: -0.644196 max: 81.3958 average: 0.00165421
ExecutionTime = 612.18 s ClockTime = 614 s

Calculating averages

Time = 0.00016

Courant Number mean: 0.00596741 max: 2.99101
GAMG: Solving for Ux, Initial residual = 0.10414, Final residual = 0.00127384, No Iterations 1
GAMG: Solving for Uy, Initial residual = 0.0648686, Final residual = 0.00136366, No Iterations 1
GAMG: Solving for Uz, Initial residual = 0.0814198, Final residual = 0.00175266, No Iterations 1
GAMG: Solving for p, Initial residual = 0.888643, Final residual = 0.0321735, No Iterations 3
time step continuity errors : sum local = 0.000540203, global = 8.97908e-06, cumulative = 3.82075e-05
GAMG: Solving for p, Initial residual = 0.370824, Final residual = 0.0159084, No Iterations 30
time step continuity errors : sum local = 0.000694453, global = 8.83416e-05, cumulative = 0.000126549

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib/libc.so.6"
#3 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#5 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::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 in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so"
#8 Foam::incompressible::LESModels:neEqEddy::correc t(Foam::tmp<Foam::GeometricField<Foam::Tensor<doub le>, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleLESModels.so"
#9 Foam::incompressible::LESModel::correct() in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleLESModels.so"
#10 main in "/opt/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/pisoFoam"
#11 __libc_start_main in "/lib/libc.so.6"
#12 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
Gleitkomma-Ausnahme

Any ideas?