CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Mistake when i calculate (http://www.cfd-online.com/Forums/openfoam-solving/58590-mistake-when-i-calculate.html)

ivanyao August 31, 2008 09:04

hi everyone, when i calculate
 
hi everyone,
when i calculate a case,there are some errors:
Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model kEpsilon
#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&)
#1 Foam::sigFpe::sigFpeHandler(int)
#2 Uninterpreted: [0xb7f7f420]
#3 Foam::turbulenceModels::kEpsilon::kEpsilon(Foam::G eometricField<foam::vector<dou ble>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::transportModel&)
#4 Foam::turbulenceModel::adddictionaryConstructorToT able<foam::turbulencemodels::k epsilon>::New(Foam::GeometricField<foam::vector<do uble>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::transportModel&)
#5 Foam::turbulenceModel::New(Foam::GeometricField<fo am::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::transportModel&)
#6 main
#7 __libc_start_main
#8 __gxx_personality_v0 at /usr/src/packages/BUILD/glibc-2.3/csu/../sysdeps/i386/elf/start.S:122
what is wrong?thank!

smehdi609 August 31, 2008 13:57

Hi, You might have initialize
 
Hi,
You might have initialized both k and epsilon with zero somewhere.
Cheers,
Mahdi

ivanyao August 31, 2008 23:06

hi, i have fix the problem,bu
 
hi,
i have fix the problem,but there is another one.when i when i calculate a case,and show:
Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model kEpsilon

Starting time loop

Time = 0.005

Courant Number mean: 0 max: 0.0777778
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 1.21994e-12, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 1.11438e-12, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 1.03725e-12, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 9.37245e-07, No Iterations 181
time step continuity errors : sum local = 8.8762e-13, global = -5.54212e-14, cumulative = -5.54212e-14
DICPCG: Solving for p, Initial residual = 4.97522e-08, Final residual = 4.97522e-08, No Iterations 0
time step continuity errors : sum local = 1.3379e-11, global = -4.27122e-12, cumulative = -4.32665e-12
#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&)
#1 Foam::sigFpe::sigFpeHandler(int)
#2 Uninterpreted: [0xb7faa420]
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&)
#4 void Foam::divide<foam::fvpatchfield,>(Foam::GeometricF ield<double,>&, Foam::GeometricField<double,> const&, Foam::GeometricField<double,> const&)
#5 Foam::tmp<foam::geometricfield<double,> > Foam::operator/<foam::fvpatchfield,>(Foam::tmp<foam::geometricfie ld<double,> > const&, Foam::GeometricField<double,> const&)
#6 Foam::turbulenceModels::kEpsilon::correct()
#7 main
#8 __libc_start_main
#9 __gxx_personality_v0 at /usr/src/packages/BUILD/glibc-2.3/csu/../sysdeps/i386/elf/start.S:122
would someone know that?

ngj September 1, 2008 02:17

Hi You are still dividing b
 
Hi

You are still dividing by zero in your turbulence model. Check that no k or epsilon is initialized to 0.

- Niels

fshak92 January 8, 2013 20:16

Hi

I'm facing to these errors after modifying my mesh(which are generated in another software).
Im working with K omega and I've checked that for zero K or omega but they are ok.Moreover, they work for other meshes.
I would appreciate if you have any suggestion.
----------------------------------------------------------------
Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::surfaceInterpolation::makeWeights() const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 Foam::surfaceInterpolation::weights() const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/simpleFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/simpleFoam"
Floating point exception (core dumped)
-------------------------------------------------------------------


All times are GMT -4. The time now is 13:56.