Hi all, I've been doing com
 

Hi all,

I've been doing computations of a transonic flow around an airfoil with high Re, so I've been using sonicTurbFoam solver for this case. Everything went well while I was using laminar turbulence model or no turbulence at all, but when I've added kEpsilon turbulence solver start to crash after the very first timestep iteration (or even in the very beginig in case of kOmegaSST model).

For the inlet BC I use fixed velocity value (M_infinity = 0.8 - 0.9), fixed small k and epsilon. For the outlet BC I've tried both fixed velocity outlet and zero gradient.

Could anybody help me to find out, where is the possible error?

Error messages for the k-epsilon case look like:
================================================
Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package hThermo<puremixture<consttransport<speciethermo<hc onstthermo<perfectgas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model kEpsilon
Creating field DpDt


Starting time loop

Time = 1e-07

Courant Number mean: 0.000261797 max: 0.115901
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for Ux, Initial residual = 0.998782, Final residual = 0.000437629, No Iterations 1000
GAMG: Solving for Uy, Initial residual = 1, Final residual = 4.432e-09, No Iterations 2
GAMG: Solving for h, Initial residual = 0.538489, Final residual = 0.0813286, No Iterations 1000
GAMG: Solving for p, Initial residual = 1, Final residual = 1.39313e-08, No Iterations 1000
GAMG: Solving for p, Initial residual = 3.71938e-08, Final residual = 1.68472e-12, No Iterations 1
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
GAMG: Solving for p, Initial residual = 1.19098e-10, Final residual = 1.19098e-10, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.7637e-10, global = 8.82502e-12, cumulative = 8.82502e-12
GAMG: Solving for p, Initial residual = 6.59328e-07, Final residual = 3.56788e-11, No Iterations 1
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.27284e-09, Final residual = 1.27284e-09, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 9.11386e-09, global = -1.70428e-10, cumulative = -1.61603e-10
#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/tls/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<foam::fvpatchfield,>(Foam::GeometricF ield<double,>&, Foam::GeometricField<double,> const&, Foam::GeometricField<double,> const&) in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTurbul enceModels.so"
#5 Foam::tmp<foam::geometricfield<double,> > Foam::operator/<foam::fvpatchfield,>(Foam::tmp<foam::geometricfie ld<double,> > const&, Foam::GeometricField<double,> const&) in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTurbul enceModels.so"
#6 Foam::compressible::turbulenceModels::kEpsilon::co rrect() in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTurbul enceModels.so"
#7 main in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/sonicTur bFoam"
#8 __libc_start_main in "/lib64/tls/libc.so.6"
#9 __gxx_personality_v0 in "/wrk/semyonov/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/sonicTur bFoam"
Floating exception
195.862u 6.226s 3:42.68 90.7% 0+0k 0+0io 2515pf+0w
=================================================