CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Native Meshers: blockMesh

mergePatchPairs vs stitchMesh

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 5, 2010, 09:43
Default mergePatchPairs vs stitchMesh
  #1
New Member
 
Paul Garlick
Join Date: Mar 2009
Location: Bournemouth, UK
Posts: 25
Rep Power: 8
pg22 is on a distinguished road
I am using the mergePatchPairs feature of blockMesh with OF16.

I have noticed a limitation in its capability of merging two cylindrical patches. This is for the case of merging a fine grid patch to a coarse grid patch. When the cylindrical patches are short (approx. 2*diameter) the merging works; when the patches are longer the merging fails, with the following error:

Zero length edge detected. Probable projection error: slave patch probably does not project onto master. Please switch on enriched patch debug for more info
.
.
.
From function void enrichedPatch::calcCutFaces() const
in file slidingInterface/enrichedPatch/enrichedPatchCutFaces.C at line 237.

FOAM aborting

I have found that it is possible to use stitchMesh to connect the two regions of the mesh (using the same long cylindrical patches). However, the zone information from blockMesh must be deleted first, so this is not ideal.

From a quick look at the code I can see that mergePatchPairs and stitchMesh use different methods for connecting the patches. However, I cannot quite see where the error comes from.

Is this a known limitation of mergePatchPairs? Could there be a fix?
pg22 is offline   Reply With Quote

Old   March 10, 2010, 06:50
Default
  #2
New Member
 
Paul Garlick
Join Date: Mar 2009
Location: Bournemouth, UK
Posts: 25
Rep Power: 8
pg22 is on a distinguished road
Some more information on this case:

Following the message to 'switch on enriched patch debug' from above I ran blockMesh again (with the long cylindrical patches) and received the following error message:
Duplicate point found in cut face. Error in the face cutting algorithm for global face 4(221561 221576 221771 221756) local face 4(1606 1607 1620 1619)
Slave size: 5952 Master size: 11904 index: 1483.
Face: 6(171216 171219 171258 221771 171294 171255)
Cut face: 6(171258 221771 171294 171297 171300 221771) local: 6(3219 13572 3231 3232 3233 13572)
Points: 6((0.00225148 -0.000603282 0.0497984) (0.00225148 -0.000603282 0.05) (0.00231096 -0.000304244 0.05) (0.00225148 -0.000603282 0.05) (0.00215347 -0.000891997 0.05) (0.00225148 -0.000603282 0.05))
This message has been recently linked to a bug Change in triangle::intersection breaks stitchMesh in 1.5 and 1.6
However, this bug was reported for stitchMesh rather than mergePatchPairs, so it may not be directly related.

Can anyone explain the conditions that generate these errors?
pg22 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
mergePatchPairs master with few slaves piprus OpenFOAM Native Meshers: blockMesh 3 November 8, 2012 10:32
Change in triangle::intersection breaks stitchMesh in 1.5 and 1.6 gschaider OpenFOAM Bugs 7 February 12, 2010 10:21
BlockMeshDict with multiple hex blocks mergePatchPairs prousa OpenFOAM Native Meshers: blockMesh 1 May 14, 2009 23:42
StitchMesh on OF15 error any fix flo OpenFOAM Bugs 0 August 20, 2008 10:27
BlockMesh documentation User Guide clearance of mergePatchPairs proposals benru OpenFOAM Native Meshers: blockMesh 0 June 9, 2008 02:22


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