[waynezoon]

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:rintStack(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:lusEqMagSqrOp2<double, double>:perator()(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:lusEqMagSqrOp2<double, double> >(double&, Foam::VectorSpace<Foam::Vector<double>, double, 3> const&, Foam:lusEqMagSqrOp2<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]

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]

Originally Posted by wyldckat

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?

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:rintStack(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:rimitiveMesh::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]

Originally Posted by waynezoon

Do i have some problems in my Mesh?

Yes, you clearly do.

[waynezoon]

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

[waynezoon]

Originally Posted by akidess

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]

Greetings to all!

@waynezoon:

Originally Posted by waynezoon

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:

Originally Posted by waynezoon

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:

Originally Posted by waynezoon

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]

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]

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]

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]

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]

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