[danie]

Hi foamers
I'm studying the flow around a circular cylinder at Re=3900 using turbFoam.
Can somebody help me?
This is the error message i got (The mesh is ok)
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
alphaEps 0.76923;
}


Starting time loop

Time = 1e-05

Courant Number mean: 3.24456e-05 max: 0.000221393
DILUPBiCG: Solving for Ux, Initial residual = 0.998615, Final residual = 5.2295e-07, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 2.50604e-08, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.999968, Final residual = 8.51587e-09, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 9.49171e-07, No Iterations 437
time step continuity errors : sum local = 2.2088e-14, global = -5.07405e-17, cumulative = -5.07405e-17
DICPCG: Solving for p, Initial residual = 0.0544752, Final residual = 9.84023e-07, No Iterations 366
time step continuity errors : sum local = 3.06552e-13, global = 6.07669e-17, cumulative = 1.00264e-17
#0 Foam::error::printStack(Foam:stream&) in "/home/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<foam::fvpatchfield,>(Foam::GeometricF ield<double,>&, Foam::GeometricField<double,> const&, Foam::GeometricField<double,> const&) in "/home/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#5 Foam::tmp<foam::geometricfield<double,> > Foam::operator/<foam::fvpatchfield,>(Foam::tmp<foam::geometricfie ld<double,> > const&, Foam::GeometricField<double,> const&) in "/home/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#6 Foam::incompressible::RASModels::kEpsilon::correct () in "/home/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#7 main in "/home/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/turbFoam"
#8 __libc_start_main in "/lib64/libc.so.6"
#9 Foam::regIOobject::readIfModified() in "/home/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/turbFoam"
Errore di virgola mobile

[ngj]

Hi Daniela,

Somewhere you have set k or epsilon equal to 0. As you are dividing by either of them you get problems. If I have a boundary where I need zero quantity I usually use 1e-11 instead.

Best regards,

Niels

P.S. I like the non-english error message in the very bottom

[danie]

thanks niels... i've run the case again right now using your suggestions and... IT WORKS!!!!! thank youuuu