|
[Sponsors] |
blockMesh problem when grading applied with cyclic patch |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
January 3, 2021, 05:54 |
blockMesh problem when grading applied with cyclic patch
|
#1 |
New Member
Henrik
Join Date: Jan 2016
Posts: 5
Rep Power: 10 |
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; 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) ); } 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 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. |
|
|
|
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 |