CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Who can help me with the OpenFoam errors? (http://www.cfd-online.com/Forums/openfoam/91807-who-can-help-me-openfoam-errors.html)

waynezoon August 23, 2011 08:18

Who can help me with the OpenFoam errors?
 
I'm create a new and simple OpenFoam case. But when i running the interFoam,i get the error below.Who can tell me what cause the errors,and how can i fix them! Thanks!! a lot!!
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.7.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.7.1-03e7e056c215
Exec : interFoam
Date : Aug 23 2005
Time : 04:55:44
Host : ubuntu
PID : 8399
Case : /home/cfd/OpenFOAM/OpenFOAM-1.7.1/tutorials1/multiphase/interFoam/les/tankcase
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Reading field p_rgh
Reading field alpha1
Reading field U
Reading/calculating face flux field phi
Reading transportProperties
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
#0 Foam::error::printStack(Foam::Ostream&) at ~/OpenFOAM/OpenFOAM-1.7.1/src/OSspecific/POSIX/printStack.C:202
#1 Foam::sigFpe::sigFpeHandler(int) at ~/OpenFOAM/OpenFOAM-1.7.1/src/OSspecific/POSIX/signals/sigFpe.C:127
#2 Uninterpreted:
#3 Foam::plusEqMagSqrOp2<double, double>::operator()(double&, double const&) const at ~/OpenFOAM/OpenFOAM-1.7.1/src/OpenFOAM/lnInclude/ops.H:76
#4 void VectorSpaceOps<3, 1>::SeqOp<double, Foam::VectorSpace<Foam::Vector<double>, double, 3>, Foam::plusEqMagSqrOp2<double, double> >(double&, Foam::VectorSpace<Foam::Vector<double>, double, 3> const&, Foam::plusEqMagSqrOp2<double, double>) at ~/OpenFOAM/OpenFOAM-1.7.1/src/OpenFOAM/lnInclude/VectorSpaceM.H:20
#5 double Foam::magSqr<Foam::Vector<double>, double, 3>(Foam::VectorSpace<Foam::Vector<double>, double, 3> const&) at ~/OpenFOAM/OpenFOAM-1.7.1/src/OpenFOAM/lnInclude/VectorSpaceI.H:311
#6 double Foam::mag<Foam::Vector<double>, double, 3>(Foam::VectorSpace<Foam::Vector<double>, double, 3> const&) at ~/OpenFOAM/OpenFOAM-1.7.1/src/OpenFOAM/lnInclude/VectorSpaceI.H:321
#7 void Foam::mag<Foam::Vector<double> >(Foam::Field<double>&, Foam::UList<Foam::Vector<double> > const&) at ~/OpenFOAM/OpenFOAM-1.7.1/src/OpenFOAM/lnInclude/FieldFunctions.C:174
#8 void Foam::mag<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ~/OpenFOAM/OpenFOAM-1.7.1/src/OpenFOAM/lnInclude/GeometricFieldFunctions.C:319
#9 Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::mag<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<Foam::Vect or<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ~/OpenFOAM/OpenFOAM-1.7.1/src/OpenFOAM/lnInclude/GeometricFieldFunctions.C:345
#10 Foam::interfaceProperties::calculateK() at ~/OpenFOAM/OpenFOAM-1.7.1/src/transportModels/interfaceProperties/interfaceProperties.C:121
#11 Foam::interfaceProperties::interfaceProperties(Foa m::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::IOdictionary const&) at ~/OpenFOAM/OpenFOAM-1.7.1/src/transportModels/interfaceProperties/interfaceProperties.C:197
#12
at ~/OpenFOAM/OpenFOAM-1.7.1/applications/solvers/multiphase/interFoam/createFields.H:94
#13 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#14
in "/home/cfd/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linuxGccDPDebug/interFoam"

wyldckat August 23, 2011 15:49

Greetings Wayne and welcome to the forum!

Well, you're asking a very generic question... therefore, I'm going to use a pre-packaged answer:
  1. Don't panic.
  2. Start reading here: Foam::error::PrintStack - post #2
You'll understand why it's a very generic question once you start reading that post ;)

Best regards and good luck!
Bruno

waynezoon August 23, 2011 21:56

Quote:

Originally Posted by wyldckat (Post 321252)
Greetings Wayne and welcome to the forum!

Well, you're asking a very generic question... therefore, I'm going to use a pre-packaged answer:
  1. Don't panic.
  2. Start reading here: Foam::error::PrintStack - post #2
You'll understand why it's a very generic question once you start reading that post ;)

Best regards and good luck!
Bruno

Thank you very much,wyldckat.
But when i use 'checkMesh',i also have some errors below. Do i have some problems in my Mesh?
Quote:

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface topology
frontwall 2871 3000 ok (non-closed singly connected)
lowwall 2871 3000 ok (non-closed singly connected)
mass-flow-inlet 841 900 ok (non-closed singly connected)
atmosphere 2871 3000 ok (non-closed singly connected)
back 2871 3000 ok (non-closed singly connected)
outwall 841 900 ok (non-closed singly connected)
Checking geometry...
Overall domain bounding box (0 0 0) (10 2 2)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (7.3134365e-17 -3.4950735e-17 8.8586217e-16) OK.
***High aspect ratio cells found, Max aspect ratio: 4.5156526e+197, number of cells 83259
<<Writing 83259 cells with high aspect ratio to set highAspectRatioCells
Minumum face area = 0.0045853684. Maximum face area = 0.0086395093. Face area magnitudes OK.
Min volume = 2e-300. Max volume = 2e-300. Total volume = 1.66518e-295. Cell volumes OK.
#0 Foam::error::printStack(Foam::Ostream&) at ~/OpenFOAM/OpenFOAM-1.7.1/src/OSspecific/POSIX/printStack.C:202
#1 Foam::sigFpe::sigFpeHandler(int) at ~/OpenFOAM/OpenFOAM-1.7.1/src/OSspecific/POSIX/signals/sigFpe.C:127
#2 Uninterpreted:
#3 in "/lib/tls/i686/cmov/libm.so.6"
#4 acos in "/lib/tls/i686/cmov/libm.so.6"
#5 Foam::primitiveMesh::checkFaceOrthogonality(bool, Foam::HashSet<int, Foam::Hash<int> >*) const at ~/OpenFOAM/OpenFOAM-1.7.1/src/OpenFOAM/meshes/primitiveMesh/primitiveMeshCheck/primitiveMeshCheck.C:503
#6
at ~/OpenFOAM/OpenFOAM-1.7.1/applications/utilities/mesh/manipulation/checkMesh/checkGeometry.C:395
#7
at ~/OpenFOAM/OpenFOAM-1.7.1/applications/utilities/mesh/manipulation/checkMesh/checkMesh.C:94
#8 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#9
in "/home/cfd/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linuxGccDPDebug/checkMesh"
Floating point exception

akidess August 24, 2011 02:48

Quote:

Originally Posted by waynezoon (Post 321273)
Do i have some problems in my Mesh?

Yes, you clearly do.

waynezoon August 24, 2011 03:09

Could you analysis it more clearly? What's the problem with the Mesh?

waynezoon August 24, 2011 03:21

Quote:

Originally Posted by akidess (Post 321297)
Yes, you clearly do.

Could you analysis it more clearly? What's the problem with the Mesh?
I have a very simple Mesh. Just a tank.So i cann't know the specific problem with my Mesh.

wyldckat August 24, 2011 04:33

Greetings to all!

@waynezoon:
Quote:

Originally Posted by waynezoon (Post 321308)
Could you analysis it more clearly? What's the problem with the Mesh?
I have a very simple Mesh. Just a tank.So i cann't know the specific problem with my Mesh.

Allow me to point out the issues...

This seems to be a clear indication of what's going wrong:
Quote:

Originally Posted by waynezoon (Post 321273)
High aspect ratio cells found, Max aspect ratio: 4.5156526e+197, number of cells 83259

Imagine that your mesh was a cube of 1x1x1 m^3. Then imagine that you stretch only one of the corners to 4.5156526e+197 m! One light year is 1e+16 m... so, that point would be 4.5e+181 light years away from the original cube...

Ironically, the volume indicated is:
Quote:

Originally Posted by waynezoon (Post 321273)
Min volume = 2e-300. Max volume = 2e-300. Total volume = 1.66518e-295

Which is... strange. Is your tank the size of molecule, atom or electron?
It almost looks like you don't even have a mesh :( Or that you are trying to simulate the whole universe...

If you give some details about the mesher you are using, then it might be easier to point out what is really going on :(

Best regards,
Bruno

waynezoon August 24, 2011 07:52

Dear Bruno:
Thanks for your precious time to analysis my poor problems.
My tank is 10*2*2 m^3. And i use gridgen to make a fluent mesh. Then using the tool 'fluentMeshToFoam' to transform the mesh .
I have known the problem!
My hard disk is full.So when i copy the Mesh file,maybe some error happened. And my Mesh is broken.
Thank you very much man . You are the frist guy who reply me in OpenFoam forume.
But when i running the case ,there some error happens.
[error]-> FOAM FATAL IO ERROR:
keyword backwall is undefined in dictionary "/home/cfd/OpenFOAM/OpenFOAM-1.7.1/tutorials1/multiphase/interFoam/ras/tankcase/0/p_rgh::boundaryField"
file: /home/cfd/OpenFOAM/OpenFOAM-1.7.1/tutorials1/multiphase/interFoam/ras/tankcase/0/p_rgh::boundaryField from line 25 to line 57.
From function dictionary::subDict(const word& keyword) const
in file db/dictionary/dictionary.C at line 456.
FOAM exiting
[/error]
What the matter with the p_rgh?

joel.lehikoinen August 24, 2011 10:39

You have a patch called backwall on which you have not defined a boundary condition for p_rgh. So you must add (or correct, if there is a typo) an entry for backwall in the 0/p_rgh file.

waynezoon August 24, 2011 21:08

Dear joel.lehikoinen
Greet to meet you and thank you for your time to helping me with the problem.
With your analysis ,i have fixed the error and my case is running now .
Thanks guys above. You are so kind!

umer.chaudrey September 5, 2011 05:00

Floating Point Exception
 
Guys,

I am experiencing a similar problem. I posted my error details on another thread in this forum, but my problem has not resolved. I was told to try changing all "0" value boundary conditions to a value such as "1e10-16". But it didn't work.

My mesh is alright, and my boundary conditions worked well for incompressible case. Can I get some help by posting my error here?

Regars,

Umer

wyldckat September 5, 2011 14:16

Hi Umer,

You could have posted the link to said thread ;)

I've answered you here: http://www.cfd-online.com/Forums/ope...ion-1-6-a.html

Best regards,
Bruno


All times are GMT -4. The time now is 10:18.