[giordy]

Dear All.

I got a problem by using laplacian Foam solver. I created the mesh by Gmsh 2.2.4 and convert into OpenFOAM by gmsh2ToFoam command. During the conversion i saw warning messages:

"
File format version: 2

Read nVerts: 3979

Read nElems: 22135
Found a physical region of region number other than 0;
omitting elementary region handling.

Automatically inverting cells with wrong orientation

Removing unused points

Constructing polyMesh
--> FOAM Warning :
From function polyMesh::polyMesh(... construct from shapes...)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 575
Found 4072 undefined faces in mesh; adding to default patch.

Performing repatching:
... repatching done.

Performing renumberMesh matrix bandwidth compression:
Band before renumbering: 17991
Band after renumbering: 288

Performing a simplified checkMesh test:
Number of non-orthogonality errors: 0. Number of severely non-orthogonal faces: 104.
--> FOAM Warning :
From function primitiveMesh::checkFaceSkewness(const bool report, labelHashSet* setPtr) const
in file meshes/primitiveMesh/primitiveMeshCheck.C at line 838
Large face skewness detected. Max skewness = 5614.73 percent.
This may impair the quality of the result.
71 highly skew faces detected.
Failed 1 mesh geometry checks.
Failed some mesh checks.
... checkMesh done.

Patches:
Patch Size Name Base type
0 3888 patch333 patch
1 168 patch334 patch
2 16 defaultFaces patch

CellZones:
Zone Size Name
0 18079 cellZone_338

Writing the final converted mesh and sets:
... done."


Anyway the Geometry seems to be successfully imported. The solver stopped during the run (before the first time step) showing the following message:

"#0 Foam::error::printStack(Foam:stream&)
#1 Foam::sigFpe::sigFpeHandler(int)
#2 Uninterpreted: [0xb7ef9420]
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&)
#4 Foam::operator/(Foam::UList<double> const&, Foam::tmp<foam::field<double> > const&)
#5 Foam::surfaceInterpolation::makeWeights() const
#6 Foam::surfaceInterpolation::weights() const
#7 Foam::surfaceInterpolation::makeDeltaCoeffs() const
#8 Foam::surfaceInterpolation::deltaCoeffs() const
#9 Foam::fvPatch::deltaCoeffs() const
#10 Foam::fixedGradientFvPatchField<double>::evaluate( )
#11 Foam::fixedGradientFvPatchField<double>::fixedGrad ientFvPatchField(Foam::fvPatch const&, Foam::Field<double> const&, Foam::dictionary const&)
#12 Foam::fvPatchField<double>::adddictionaryConstruct orToTable<foam::fixedgradientf vpatchfield<double> >::New(Foam::fvPatch const&, Foam::Field<double> const&, Foam::dictionary const&)
#13 Foam::fvPatchField<double>::New(Foam::fvPatch const&, Foam::Field<double> const&, Foam::dictionary const&)
#14 Foam::GeometricField<double,>::GeometricBoundaryFi eld::GeometricBoundaryField(Fo am::fvBoundaryMesh const&, Foam:: Field<double> const&, Foam::dictionary const&)
#15 Foam::GeometricField<double,>::readField(Foam::Ist ream&)
#16 Foam::GeometricField<double,>::GeometricField(Foam ::IOobject const&, Foam::fvMesh const&)
#17 main
#18 __libc_start_main
#19 __gxx_personality_v0 at /usr/src/packages/BUILD/glibc-2.3/csu/../sysdeps/i386/elf/start.S:122"


Does someone know how to solve this problem?
Thanks in advance,

Stefano

[albcem]

It sounds like your mesh is way off from ideal; your skewness warning is way off and I guarantee OpenFOAM will not run. If you like, post the geo file and I can have a look at it with GMSH.

Jim

[giordy]

Thank you for your kind reply.
In reality i'm not able to upload file in the forum, sorry. If you have private mail i could send you the geo file.

Stefano

[steven]

You have not meshed it right as Jim has pointed out. Its really easy to upload a file here. Have a read through documentation section.

Steve