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/)
-   -   %230 FoamerrorprintStackFoam%3cIMG SRC%3d%22httpopenfoamcfdonlinecomforumclipartproud gif%22 ALT%3d%22O%22 BORDER%3d0%3estreamamp (http://www.cfd-online.com/Forums/openfoam-solving/59360-230-foamerrorprintstackfoam-3cimg-src-3d-22httpopenfoamcfdonlinecomforumclipartproudgif-22-alt-3d-22o-22-border-3d0-3estreamamp.html)

guimch July 17, 2007 07:18

Hello ! I am a new openFoam u
 
Hello !
I am a new openFoam user and I must admit that it is still much stronger than me ...
I use fluentMeshtoFoam : paraView seems to show me that my mesh is correct.
In my case, I have to model an inlet swirl in a cubic volume. After 5 days of hard working, my inlet condition is ok (according to paraview).
I use the simpleFoam solver (steadystate, turbulent, incompressible).
It works properly as I turn turbulence "off" in turbulenceProperties but as I said before, I need results in the turbulent case ...
My trouble is the follwing error message :


Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model CrossPowerLaw
Selecting turbulence model kEpsilon

Starting time loop

Time = 1

DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.0291, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.0176, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 0.0124, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.00787, No Iterations 37
time step continuity errors : sum local = 1.25e+04, global = 137, cumulative = 137
#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&)
#1 Foam::sigFpe::sigFpeHandler(int)
#2 Uninterpreted: [0x2ee420]
#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 ../sysdeps/i386/elf/start.S:122
/home/lodier/OpenFOAM/OpenFOAM-1.4/bin/runFoamXHB: line 139: 10979 Killed FoamXHostBrowser -ORBNamingAddr $myIOP

I can't sort out this problem. So, it would be most helpfull if you could give me any advices you might have !
Thanks a lot
Guillaume

cedric_duprat July 17, 2007 08:05

hi Guillaume, I can't help
 
hi Guillaume,

I can't help you very well but, go for a look there:
http://www.cfd-online.com/OpenFOAM_D...tml?1184091458
or try to find "continuity error" (which seems to be your pb)
there are small tips to correct your problem.

Good luck,
Cedric

guimch July 17, 2007 09:10

Thank you for your answer Cedr
 
Thank you for your answer Cedric. I had a look overthere, I tried to change the properties in the same way but unfortunatly, it doesn't work better, I get the same error message.
I am really lost. Is there someone to help me to find my way ?
Guillaume

gschaider July 25, 2007 12:52

Hi Guillaume! Are in your I
 
Hi Guillaume!

Are in your Initial/Boundary-conditons k or epsilon 0? That would explain the division by zero you have.

Bernhard

guimch July 26, 2007 04:14

Thanks a lot Bernhard ! That
 
Thanks a lot Bernhard !
That was exactly the case. I have got an other question : the value that I put for k and epsilon are a priori value, aren't they ? And these values will permit to converge quickly if they are well chosen ?
thank you
Guillaume

nishant_hull October 25, 2007 09:26

Hi all.. I am currently fa
 
Hi all..

I am currently facing this error on my openfoaminstallation. I wil appreciate if some body can help me in this regard.

#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0x110420]
#3 void Foam::fvc::surfaceIntegrate<foam::vector<double> >(Foam::Field<foam::vector<double> >&, Foam::GeometricField<foam::vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libfiniteVolume.so"
#4 Foam::tmp<foam::geometricfield<foam::vector<double >, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::surfaceIntegrate<foam::vector<double> >(Foam::GeometricField<foam::vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libfiniteVolume.so"
#5 Foam::tmp<foam::geometricfield<foam::vector<double >, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::surfaceIntegrate<foam::vector<double> >(Foam::tmp<foam::geometricfield<foam::vector<doub le>, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libfiniteVolume.so"
#6 Foam::fv::gaussDivScheme<foam::tensor<double> >::fvcDiv(Foam::GeometricField<foam::tensor<double >, Foam::fvPatchField, Foam::volMesh> const&) in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libfiniteVolume.so"
#7 Foam::tmp<foam::geometricfield<foam::innerproduct< foam::vector<double>, Foam::Tensor<double> >::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div<foam::tensor<double> >(Foam::GeometricField<foam::tensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libincompressibleTurbule nceModels.so"
#8 Foam::tmp<foam::geometricfield<foam::innerproduct< foam::vector<double>, Foam::Tensor<double> >::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div<foam::tensor<double> >(Foam::GeometricField<foam::tensor<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libincompressibleTurbule nceModels.so"
#9 Foam::tmp<foam::geometricfield<foam::innerproduct< foam::vector<double>, Foam::Tensor<double> >::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div<foam::tensor<double> >(Foam::tmp<foam::geometricfield<foam::tensor<doub le>, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libincompressibleTurbule nceModels.so"
#10 Foam::turbulenceModels::kEpsilon::divR(Foam::Geome tricField<foam::vector<double> , Foam::fvPatchField, Foam::volMesh>&) const in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libincompressibleTurbule nceModels.so"
#11 main in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/simpleFoam"
#12 __libc_start_main in "/lib/libc.so.6"
#13 Foam::regIOobject::readIfModified() in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/simpleFoam"
Floating point exception

gabriel October 25, 2007 17:52

To Guillaume, you are right
 
To Guillaume,

you are right, a good first estimation of the internal fields for kappa and epsilon is important. I made the expirience, that using bad values, convergence can becom difficult.

Gabriel


All times are GMT -4. The time now is 04:52.