CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[mesh manipulation] StitchMesh sigFpe / bad point

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

Like Tree2Likes
  • 1 Post By wyldckat
  • 1 Post By liquidspoon

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 24, 2015, 18:34
Default StitchMesh sigFpe / bad point
  #1
Member
 
Alex
Join Date: Jun 2011
Posts: 33
Rep Power: 15
liquidspoon is on a distinguished road
Hi all,

I am trying to use mergeMeshes and StitchMesh to attach an inlet port to an annular domain, but am getting stumped on the StitchMesh stage. I would really appreciate assistance on this.

I make the annulus first
Code:
blockMesh
mirrorMesh
and everything looks good (figure):
Step 1.jpg

Then I make the port that connects to the side in another case directory (MeshCase), and that comes out fine:
Code:
blockMesh
(see figure)
Step 2.jpg

Then I run:

Code:
mergeMeshes -overwrite . MeshCase
and both meshes are now combined in the top level mesh. As you can see from the third figure, the grids on both meshes line up perfectly.
Step 3.jpg

Finally, I need to get rid of the dummy patch (inner portion of the inlet port that coincides with the outer boundary of the big mesh). I have tried to use stitchMesh to make this work, but can not figure it out.

Calling:
Code:
stitchMesh -partial Wall Dummy
Throws a floating point exception in the projectPoints() stage.
Code:
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::slidingInterface::projectPoints() const at ??:?
#4  Foam::slidingInterface::changeTopology() const at ??:?
#5  Foam::polyTopoChanger::changeTopology() const at ??:?
#6  Foam::polyTopoChanger::changeMesh(bool, bool, bool, bool) at ??:?
#7  ? at ??:?
#8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9  ? at ??:?
Floating point exception (core dumped)
If I unset FOAM_SIGFPE, and rerun stitchMesh, I get the specific message:

Code:
Bad points:(0.137781 0.980365 0.02) (0.137781 0.980365 0.02) (-0.30265 -1.96664 0.02)

    From function void plane::calcPntAndVec
(
    const point&,
    const point&,
    const point&
)
I'm not really sure what is going on in this stage of the stitchMesh program. I have also tried playing around with toleranceDict, different stich options (perfect, partial, etc.), but no success.

I have attached the case. StitchTest.zip
liquidspoon is offline   Reply With Quote

Old   November 29, 2015, 07:38
Default
  #2
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,981
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Hi Alex,

The problem is that you have a single patch for the main mesh, leading stitchMesh to get confused about which points from the wall to be used for stitching.

The following did the trick:
Code:
blockMesh
mirrorMesh
blockMesh -case MeshCase
autoPatch -overwrite 45
mergeMeshes -overwrite . MeshCase
stitchMesh -overwrite auto1 Dummy
The "-partial" option can't be used in this scenario, because the point correspondence is of type "integral": http://www.openfoam.org/mantisbt/view.php?id=1767#c5338 - see comment #5338

Best regards,
Bruno
liquidspoon likes this.
__________________
wyldckat is offline   Reply With Quote

Old   November 29, 2015, 15:19
Default
  #3
Member
 
Alex
Join Date: Jun 2011
Posts: 33
Rep Power: 15
liquidspoon is on a distinguished road
Hi Bruno,

Thank you very much for the detailed response and assistance. I am eager to try this out when I get back to the office.

You are a real asset to the community. I owe you a beer if we are ever in the same town for a conference :-).
wyldckat likes this.
liquidspoon is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[mesh manipulation] Duplicate point found in cut face - stitchMesh - problem RugbyGandalf OpenFOAM Meshing & Mesh Conversion 5 April 2, 2020 15:45
[mesh manipulation] stitchMesh: multiple meshes GerhardHolzinger OpenFOAM Meshing & Mesh Conversion 3 August 25, 2017 13:43
an odd(at least for me!) reconstructPar error on a field immortality OpenFOAM Running, Solving & CFD 3 June 3, 2013 23:36
[mesh manipulation] stitchMesh Dublicate point ... besto OpenFOAM Meshing & Mesh Conversion 1 October 4, 2011 22:48
Problems of Duns Codes! Martin J Main CFD Forum 8 August 15, 2003 00:19


All times are GMT -4. The time now is 22:01.