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

[snappyHexMesh] snappyHexMesh - cyclic boundaries

Register Blogs Community New Posts Updated Threads Search

Like Tree9Likes
  • 6 Post By Ruli
  • 3 Post By aylalisa

 
 
LinkBack Thread Tools Search this Thread Display Modes
Prev Previous Post   Next Post Next
Old   February 13, 2014, 14:17
Default snappyHexMesh - cyclic boundaries
  #1
Senior Member
 
Join Date: Nov 2012
Location: Bavaria
Posts: 145
Rep Power: 13
aylalisa is on a distinguished road
Hi to All!

I've created a mesh of a simple channel with snappyHexMesh (in parallel).
In principle meshing with SHM worked fine, simulation ran successfully.
Primarily the model had these boundary regions:
inlet, outlet : patch
wallLeft, wallRight, top, bottom: wall

Next I've changed the boundaries of LeftWall and RightWall to cyclic in blockMeshDict and I've adapted the 0/* files accordingly.




log-file 'log.snappyHexMesh_par' discloses this error:

Code:
nProcs : 8
Slaves : 
7
(
"linux-0ze2.14255"
"linux-0ze2.14256"
"linux-0ze2.14258"
"linux-0ze2.14259"
"linux-0ze2.14260"
"linux-0ze2.14261"
"linux-0ze2.14262"
)

Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Read mesh in = 0.14 s

Overall mesh bounding box  : (0 0 0) (3.14 1.57 1)
Relative tolerance         : 1e-06
Absolute matching distance : 3.65027e-06

Reading refinement surfaces.
Read refinement surfaces in = 0.01 s

Reading refinement shells.
Refinement level according to distance to rod
    level 2 for all cells within 0.2 metre.
Read refinement shells in = 0 s

Setting refinement level of surface to be consistent with shells.
Checked shell refinement in = 0 s

Reading features.
Read features in = 0 s


Determining initial surface intersections
-----------------------------------------

Edge intersection testing:
    Number of edges             : 244200
    Number of edges to retest   : 244200
    Number of intersected edges : 542
Calculated surface intersections in = 0.19 s

Initial mesh : cells:80000  faces:244200  points:84240
Cells per refinement level:
    0    80000

Adding patches for surface regions
----------------------------------

Patch    Type    Region
-----    ----    ------
rod:

6    wall    rod_region0

Added patches in = 0.01 s

Selecting decompositionMethod hierarchical

Refinement phase
----------------

[4] Found point (1 0.5 0.5) in cell 203 on processor 4

Surface refinement iteration 0
------------------------------

Marked for refinement due to surface intersection          : 718 cells.
Marked for refinement due to curvature/regions             : 0 cells.
Determined cells to refine in = 0.07 s
Selected for refinement : 718 cells (out of 80000)
[0] Cannot find point in pts1 matching point 0 coord:(0.019625 0 0.02) in pts0 when using tolerance 0.0280204
[0] Searching started from:0 in pts1
.
.
.
[7]     Compared coord: (3.12038 0 0.98) at index 509 with difference to point 1.57051
[7] Cannot find point in pts1 matching point 485 coord:(3.12038 1.57 0.98) in pts0 when using tolerance 0.0280204
[7] Searching started from:509 in pts1
[7]     Compared coord: (3.12038 0 0.98) at index 509 with difference to point 1.57
[7] processorPolyPatch::order : Writing faces to OBJ file procBoundary7to5throughwallLeft_faces.obj
[3] processorPolyPatch::order : Dumping newly found match as lines between corresponding face centres to OBJ file "/home/lisa/OpenFOAM/lisa-2.2.2/run/modelle_2/Anfangsbedingungen_obstacle_rod/1_simpleFoam_cyclic/processor3/procBoundary3to1throughwallLeft_faceCentresConnections.obj"
[3] --> FOAM Serious Error : 
[3]     From function processorPolyPatch::order(const primitivePatch&, labelList&, labelList&) const
[3]     in file meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C at line 1014
[3]     in patch:procBoundary3to1throughwallLeft : Cannot match vectors to faces on both sides of patch
[3]     masterCtrs[0]:(1.58963 0 0.02)
[3]     ctrs[0]:(1.58963 1.57 0.02)
[3]     Check your topology changes or maybe you have multiple separated (from cyclics) processor patches
[3]     Continuing with incorrect face ordering from now on
[7] processorPolyPatch::order : Dumping newly found match as lines between corresponding face centres to OBJ file "/home/lisa/OpenFOAM/lisa-2.2.2/run/modelle_2/Anfangsbedingungen_obstacle_rod/1_simpleFoam_cyclic/processor7/procBoundary7to5throughwallLeft_faceCentresConnections.obj"
[7] --> FOAM Serious Error : 
[7]     From function processorPolyPatch::order(const primitivePatch&, labelList&, labelList&) const
[7]     in file meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C at line 1014
[7]     in patch:procBoundary7to5throughwallLeft : Cannot match vectors to faces on both sides of patch
[7]     masterCtrs[0]:(1.58963 0 0.54)
[7]     ctrs[0]:(1.58963 1.57 0.54)
[7]     Check your topology changes or maybe you have multiple separated (from cyclics) processor patches
[7]     Continuing with incorrect face ordering from now on
Could anybody point out the problem?

Thanks!
Ayla
aylalisa is offline   Reply With Quote

 


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] Problem with using createPatchDict to define cyclic boundaries kaszt OpenFOAM Meshing & Mesh Conversion 0 April 1, 2016 21:18
Possible createPatch/createBaffles bug? simpomann OpenFOAM Bugs 2 July 15, 2014 07:07
[snappyHexMesh] snappyHexMesh and cyclic boundaries Ruli OpenFOAM Meshing & Mesh Conversion 2 December 9, 2013 06:51
Wrong wall distance with cyclic boundaries sebastian OpenFOAM Bugs 4 October 31, 2012 10:24
[snappyHexMesh] snappyHexMesh with cyclic patches? Axel_T OpenFOAM Meshing & Mesh Conversion 1 October 12, 2010 10:03


All times are GMT -4. The time now is 06:46.