|
[Sponsors] |
%230 FoamerrorprintStackFoam%3cIMG SRC%3d%22httpopenfoamcfdonlinecomforumclipartproud gif%22 ALT%3d%22O%22 BORDER%3d0%3estreamamp |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 17, 2007, 07:18 |
Hello !
I am a new openFoam u
|
#1 |
New Member
Guillaume Lodier
Join Date: Mar 2009
Location: Rouen, France
Posts: 4
Rep Power: 17 |
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:stream&) #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 |
|
July 17, 2007, 08:05 |
hi Guillaume,
I can't help
|
#2 |
Senior Member
Cedric DUPRAT
Join Date: Mar 2009
Location: Nantes, France
Posts: 195
Rep Power: 17 |
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 |
|
July 17, 2007, 09:10 |
Thank you for your answer Cedr
|
#3 |
New Member
Guillaume Lodier
Join Date: Mar 2009
Location: Rouen, France
Posts: 4
Rep Power: 17 |
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 |
|
July 25, 2007, 12:52 |
Hi Guillaume!
Are in your I
|
#4 |
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 |
Hi Guillaume!
Are in your Initial/Boundary-conditons k or epsilon 0? That would explain the division by zero you have. Bernhard
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request |
|
July 26, 2007, 04:14 |
Thanks a lot Bernhard !
That
|
#5 |
New Member
Guillaume Lodier
Join Date: Mar 2009
Location: Rouen, France
Posts: 4
Rep Power: 17 |
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 |
|
October 25, 2007, 09:26 |
Hi all..
I am currently fa
|
#6 |
Senior Member
Nishant
Join Date: Mar 2009
Location: Glasgow, UK
Posts: 166
Rep Power: 17 |
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:stream&) 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
__________________
Thanks and regards, Nishant |
|
October 25, 2007, 17:52 |
To Guillaume,
you are right
|
#7 |
New Member
Gabriel Barroso
Join Date: Mar 2009
Posts: 23
Rep Power: 17 |
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 |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|