[openfoam_user]

Dear users,

I have read many threads about 'bounding problems for k and epsilon', but without finding how to solve it.

I am trying to compute a free surface calculation with a boat.

The mesh was generated using ICEMCFD Hexa.

The cae crashes after 2-4 iterations with the following message :

MULES: Solving for gamma
Liquid phase volume fraction = 0.499998 Min(gamma) = -2.96063e-12 Max(gamma) = 1.01884
MULES: Solving for gamma
Liquid phase volume fraction = 0.499998 Min(gamma) = -2.96285e-12 Max(gamma) = 1.01883
MULES: Solving for gamma
Liquid phase volume fraction = 0.499998 Min(gamma) = -2.96507e-12 Max(gamma) = 1.01883
MULES: Solving for gamma
Liquid phase volume fraction = 0.499998 Min(gamma) = -2.96729e-12 Max(gamma) = 1.01882
DILUPBiCG: Solving for Ux, Initial residual = 0.611633, Final residual = 8.16394e-07, No Iterations 625
DILUPBiCG: Solving for Uy, Initial residual = 0.673241, Final residual = 8.43185e-07, No Iterations 631
DILUPBiCG: Solving for Uz, Initial residual = 0.927937, Final residual = 8.75002e-07, No Iterations 564
DICPCG: Solving for pd, Initial residual = 1, Final residual = 0.038253, No Iterations 4
DICPCG: Solving for pd, Initial residual = 6.54561e-09, Final residual = 6.54561e-09, No Iterations 0
DICPCG: Solving for pd, Initial residual = 6.54561e-09, Final residual = 6.54561e-09, No Iterations 0
DICPCG: Solving for pd, Initial residual = 1.36844e-07, Final residual = 1.2005e-08, No Iterations 1
DICPCG: Solving for pd, Initial residual = 3.54892e-08, Final residual = 3.54892e-08, No Iterations 0
DICPCG: Solving for pd, Initial residual = 3.54892e-08, Final residual = 3.54892e-08, No Iterations 0
DICPCG: Solving for pd, Initial residual = 1.63417e-07, Final residual = 1.93669e-08, No Iterations 1
DICPCG: Solving for pd, Initial residual = 4.80168e-08, Final residual = 4.80168e-08, No Iterations 0
DICPCG: Solving for pd, Initial residual = 4.80168e-08, Final residual = 4.80168e-08, No Iterations 0
time step continuity errors : sum local = 0.000631323, global = 8.45631e-07, cumulative = 8.45629e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.870335, Final residual = 32077.1, No Iterations 2001
bounding epsilon, min: -6.62374e+16 max: 1.14934e+17 average: -7.2268e+11
DILUPBiCG: Solving for k, Initial residual = 0.995697, Final residual = 8.46144e-05, No Iterations 2001
bounding k, min: -2.04339e+12 max: 8.30451e+14 average: 2.41778e+11
ExecutionTime = 282.88 s ClockTime = 284 s

Courant Number mean: 0.252986 max: 452593
deltaT = 3.58683e-11
Time = 0.00503247

MULES: Solving for gamma
[1] [3] ##00 Foam::error::printStack(Foam:stream&)Foam::error::printStack(Foam:stream&) in "/shared/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #1 Foam::sigFpe::sigFpeHandler(int) in "/shared/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #1 Foam::sigFpe::sigFpeHandler(int)[2] #0 Foam::error::printStack(Foam:stream&)[0] in "/shared/OpenF#O0 AFoam::error::printStack(Foam:stream&)M/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #2 in [5] "#/0share dFoam::error::printStack(Foam:stream&)/OpenFOAM/Ope[4] n#F0O AFoam::error::printStack(Foam:stream&)M-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #2 ?? in "/shared/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[2] #1 Foam::sigFpe::sigFpeHandler(int) in "/shared/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[0] #1 Foam::sigFpe::sigFpeHandler(int)?? in in "/shared/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[5] #1 Foam::sigFpe::sigFpeHandler(int)"/shared/OpenFOAM in /O"p/elniFbO6A4M/-l1i.b5c/.lsiob./6l"i
n[1] u#x36 4Gvoid Foam::MULES::limiter<foam::onefield,>(Foam::Field< double>&, Foam::oneField const&, Foam::GeometricField<double,> const&, Foam::GeometricField<double,> const&, Foam::GeometricField<double,> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int)ccDPOpt/libOpenFOAM.so"
[4] #1 Foam::sigFpe::sigFpeHandler(int) in "/shared/OpenFOAM/OpenFOAM-1.5 in /"l/isbh/alriendu/xO6p4eGncFcODAPMO/pOtp/elniFbOOApMe-n1F.O5A/Ml.isbo/"l
i[2] n#u2x 64GccDPOpt/libOpenFOAM.so"
[0] #2 in "/lib64/libc.so.6"

The error message is much longer. This is only the beginning.

I have tried to play with many parameters (deltaT, maxDeltaT, k and setFields, ...) but without success !!!!!!!!!!!!

Stephane.