CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

blockMesh problem when grading applied with cyclic patch

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 3, 2021, 05:54
Default blockMesh problem when grading applied with cyclic patch
  #1
New Member
 
Henrik
Join Date: Jan 2016
Posts: 5
Rep Power: 10
broccolibadger is on a distinguished road
blockMesh with grading works like a charm without grading as per attached images A & B.

With grading specified and changing cyclic to wall, the mesh is generated beautifully. Problem starts when cyclic and grading is combined, see image C.

blockMeshDict is attached.

Gradings:
Code:
/*
gz0 1;
gz1 1;
gr 1;
*/

gz0 .1;
gz1
(
    (.3 3 1)
    (.7 2 5)
);
gr 10;
Patch specifications:
Code:
cyclicA
    {
        type            cyclic;
        neighbourPatch  cyclicB;

        transform       rotational;
        rotationAxis    (0 0 1);
        rotationCentre  (0 0 0);
//        rotationAngle   270;
        matchTolerance  1e-5;

        faces
            (
                (7 32 29 0)
                (32 10 3 29)
                (10 12 5 3)
                (10 17 19 12)
                (17 24 26 19)
                (35 38 24 17)
                (14 21 38 35)
            );
    }

    cyclicB
    {
        type            cyclic;
        neighbourPatch  cyclicA;

        transform       rotational;
        rotationAxis    (0 0 1);
        rotationCentre  (0 0 0);
//        rotationAngle   90;
        matchTolerance  1e-5;

        faces
            (
                (0 28 31 7)
                (28 1 8 31)
                (1 2 9 8)
                (9 16 15 8)
                (15 16 23 22)
                (22 37 15 34)
                (34 37 21 14)
            );
    }
Error message:
Code:
cyclicPolyPatch::calcTransforms : Writing cyclicA faces to OBJ file "cyclicA_faces.obj"
cyclicPolyPatch::calcTransforms : Writing cyclicB faces to OBJ file "cyclicB_faces.obj"
cyclicPolyPatch::calcTransforms : Writing coupled face centres as lines to "cyclicA_to_cyclicB.obj"


--> FOAM FATAL ERROR: 
face 1 area does not match neighbour by 25.44563% -- possible face ordering problem.
patch:cyclicA my area:4.1126e-05 neighbour area:3.184237e-05 matching tolerance:1e-05
Mesh face:6011 fc:(0 0.002292 0.01345746)
Neighbour fc:(0.000764 0 0.03733454)
If you are certain your matching is correct you can increase the 'matchTolerance' setting in the patch dictionary in the boundary file.
Rerun with cyclic debug flag set for more information.

    From void Foam::cyclicPolyPatch::calcTransforms(const primitivePatch&, const pointField&, const vectorField&, const pointField&, const vectorField&)
    in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 215.

FOAM exiting
I'm using .com-version v2006

To me this looks like a bug in how cyclic is implemented. I can't seem to find any other reason for this not to work.
Attached Images
File Type: png A.png (180.6 KB, 4 views)
File Type: png B.png (113.0 KB, 5 views)
File Type: png C.png (104.6 KB, 5 views)
Attached Files
File Type: txt blockMeshDict.txt (5.3 KB, 2 views)
broccolibadger is offline   Reply With Quote

Reply


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
openfoam ext4.0 and problem in cyclic condition Hgholami OpenFOAM Pre-Processing 2 February 19, 2019 11:53
[Commercial meshers] Mesh conversion problem (fluent3DMeshToFoam) Aadhavan OpenFOAM Meshing & Mesh Conversion 2 March 8, 2018 02:47
[OpenFOAM.org] Compile OF 2.3 on Mac OS X .... the patch gschaider OpenFOAM Installation 225 August 25, 2015 20:43
createPatch Segmentation Fault (CORE DUMPED) sam.ho OpenFOAM Pre-Processing 2 April 21, 2014 03:01
Problem in defining patch deformation paul b OpenFOAM Programming & Development 3 April 27, 2010 01:31


All times are GMT -4. The time now is 16:42.