[Hanseny]

Hi,

I am running a 2-D laminar airfoil S809 in Simplefoam with S-A turbulent model and debug mode, It just crashed when 68 steps. I copy couples of steps as below:

Time = 65

smoothSolver: Solving for Ux, Initial residual = 0.15809, Final residual = 0.00202525, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.238226, Final residual = 0.00310446, No Iterations 2
GAMG: Solving for p, Initial residual = 0.951543, Final residual = 0.0731533, No Iterations 1
time step continuity errors : sum local = 1.07535e+26, global = -8.07586e+10, cumulative = -7.80351e+10
smoothSolver: Solving for nuTilda, Initial residual = 0.221721, Final residual = 0.00212716, No Iterations 6
bounding nuTilda, min: -0.0282978 max: 6.34858e+19 average: 2.00363e+15
ExecutionTime = 29.77 s ClockTime = 30 s

Time = 66

smoothSolver: Solving for Ux, Initial residual = 0.247476, Final residual = 0.0112385, No Iterations 6
smoothSolver: Solving for Uy, Initial residual = 0.212603, Final residual = 0.00806766, No Iterations 6
GAMG: Solving for p, Initial residual = 0.995719, Final residual = 0.0864595, No Iterations 1
time step continuity errors : sum local = 3.60256e+32, global = 2.21104e+16, cumulative = 2.21103e+16
smoothSolver: Solving for nuTilda, Initial residual = 0.000146617, Final residual = 2.40179e-06, No Iterations 10
bounding nuTilda, min: -2355.47 max: 7.5393e+22 average: 1.78911e+18
ExecutionTime = 29.95 s ClockTime = 30 s


Time = 67

smoothSolver: Solving for Ux, Initial residual = 0.193772, Final residual = 0.00562729, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.161721, Final residual = 0.00563807, No Iterations 2
GAMG: Solving for p, Initial residual = 0.99604, Final residual = 2.502e+16, No Iterations 1000
time step continuity errors : sum local = 1.81752e+53, global = 5.00108e+37, cumulative = 5.00108e+37
smoothSolver: Solving for nuTilda, Initial residual = 0.472311, Final residual = 1.12189e+12, No Iterations 1000
ExecutionTime = 40.09 s ClockTime = 40 s

Time = 68

smoothSolver: Solving for Ux, Initial residual = 0.202715, Final residual = 0.00219088, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.0602311, Final residual = 0.00212992, No Iterations 2
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const at ??:?
#6 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::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#7 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#10
at ??:?
#11
at ??:?
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13
at ??:?
Floating point exception (core dumped)

It looks like always got "time step continuity errors" and nuTilda blowed up. I have no idea why.

thank you very much for your help!