CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Mesh generation

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   October 29, 2008, 06:31
Default Dear All. I got a problem b
  #1
New Member
 
stefano zinna
Join Date: Mar 2009
Posts: 14
Rep Power: 7
giordy is on a distinguished road
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
giordy is offline   Reply With Quote

Old   October 29, 2008, 21:47
Default It sounds like your mesh is wa
  #2
Member
 
Cem Albukrek
Join Date: Mar 2009
Posts: 50
Rep Power: 7
albcem is on a distinguished road
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
albcem is offline   Reply With Quote

Old   October 30, 2008, 15:46
Default Thank you for your kind reply.
  #3
New Member
 
stefano zinna
Join Date: Mar 2009
Posts: 14
Rep Power: 7
giordy is on a distinguished road
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
giordy is offline   Reply With Quote

Old   October 30, 2008, 17:13
Default You have not meshed it right a
  #4
New Member
 
Steven Parole
Join Date: Mar 2009
Posts: 3
Rep Power: 7
steven is on a distinguished road
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
steven is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
mesh generation srinivas CD-adapco 1 September 13, 2007 01:46
mesh generation srikanth Main CFD Forum 1 January 26, 2006 06:02
RE. Mesh Generation Help Patel Amit R. FLUENT 0 July 12, 2003 03:03
mesh generation Matthias Main CFD Forum 0 February 4, 2002 10:22
MESH GENERATION SANJIB DAS SHARMA Main CFD Forum 7 May 13, 2000 08:16


All times are GMT -4. The time now is 14:38.