Turbfoam error

 Register Blogs Members List Search Today's Posts Mark Forums Read

 July 30, 2008, 06:45 Hi foamers I'm studying the f #1 New Member   Daniela Join Date: Mar 2009 Location: Italy Posts: 4 Rep Power: 17 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&, Foam::UList const&, Foam::UList const&) in "/home/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so" #4 void Foam::divide(Foam::GeometricF ield&, Foam::GeometricField const&, Foam::GeometricField const&) in "/home/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libincompressibleRASModels.so" #5 Foam::tmp > Foam::operator/(Foam::tmp > const&, Foam::GeometricField 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

 July 30, 2008, 07:08 Hi Daniela, Somewhere you h #2 Senior Member   Niels Gjoel Jacobsen Join Date: Mar 2009 Location: Copenhagen, Denmark Posts: 1,901 Rep Power: 37 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 __________________ Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.

 July 30, 2008, 07:45 thanks niels... i've run the c #3 New Member   Daniela Join Date: Mar 2009 Location: Italy Posts: 4 Rep Power: 17 thanks niels... i've run the case again right now using your suggestions and... IT WORKS!!!!! thank youuuu