Hi everybody, I'm strugglin
 

Hi everybody,

I'm struggling with a rae profile analysed for Ma=0.25 (aoa=0 degree). Although the boundary conditions seem to be correct (indeed potentialfoam did not have any problem with the same BC) this is the message I get:

------------------------------------
Time = 0.066

DILUPBiCG: Solving for Ux, Initial residual = 0.439752, Final residual = 4.10432e-06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.453736, Final residual = 5.53212e-07, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.987358, Final residual = 0.000979867, No Iterations 192
time step continuity errors : sum local = 2.13623e+86, global = -6.30715e+72, cumulative = -5.9625e+72
DILUPBiCG: Solving for epsilon, Initial residual = 0.0513206, Final residual = 8.37324e-05, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 2.87352e-09, Final residual = 2.87352e-09, No Iterations 0
ExecutionTime = 39.91 s ClockTime = 41 s

Time = 0.067

DILUPBiCG: Solving for Ux, Initial residual = 0.373545, Final residual = 9.78289e-05, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.435162, Final residual = 1.79789e-05, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.98127, Final residual = 0.000992254, No Iterations 197
time step continuity errors : sum local = 5.44038e+87, global = 1.91227e+74, cumulative = 1.85265e+74
DILUPBiCG: Solving for epsilon, Initial residual = 0.303376, Final residual = 3.81456e-05, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 3.13267e-08, Final residual = 3.13267e-08, No Iterations 0
ExecutionTime = 40.49 s ClockTime = 42 s

Time = 0.068

DILUPBiCG: Solving for Ux, Initial residual = 0.434953, Final residual = 3.85775e-06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.45262, Final residual = 7.94863e-07, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.987109, Final residual = 0.000802516, No Iterations 208
time step continuity errors : sum local = 1.21666e+89, global = -2.62405e+75, cumulative = -2.43879e+75
#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/home/nettis/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/nettis/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/nettis/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 Foam::fvMatrix<double>::solve(Foam::Istream&) in "/home/nettis/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::lduMatrix::solverPerformance Foam::solve<double>(Foam::tmp<foam::fvmatrix<doubl e> > const&) in "/home/nettis/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libincompressibleTurbu lenceModels.so"
#6 Foam::turbulenceModels::LaunderSharmaKE::correct() in "/home/nettis/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libincompressibleTurbu lenceModels.so"
#7 main in "/home/nettis/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/simpleFoa m"
#8 __libc_start_main in "/lib64/libc.so.6"
#9 Foam::regIOobject::readIfModified() in "/home/nettis/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/simpleFoa m"
Errore di virgola mobile

-----------------------------------

I'm using a Low Re model (LaunderSharmaKE) with very low starting value of k,e and nutilda (1e-10) in the whole domain. Anyway I don't think these value could influence the problem since I get approximately the same message with turbulence off.
Finally I do not think it is due to the solver since it is a classic case of incompressible steady-state flow.
Therefore I hope somebody can give me an advice.
thanks in advance.

dino