[y.liu]

Hi guys,

I have created an O-grid for a twisted hexagonal cylinder in Pointwise v18.5 and exported the mesh as a polyMesh folder into openfoam-v2112 (the .com version).

cylinder_spanwise.PNG

cylinder.PNG

O-grid.PNG

When I try and make the z-axis spanwise endwalls cyclic, I get the following error:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2112                                  |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : _14aeaf8dab-20211220 OPENFOAM=2112 version=v2112
Arch   : "LSB;label=32;scalar=64"
Exec   : checkMesh
Date   : Nov 30 2022
Time   : 00:17:56
Host   : 
PID    : 
I/O    : uncollated
Case   : 
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


--> FOAM FATAL ERROR: (openfoam-2112)
face 647 area does not match neighbour by 0.013587605% -- possible face ordering problem.
patch:back my area:0.00030710446 neighbour area:0.00030714619 matching tolerance:0.0001
Mesh face:13372775 fc:(-0.061516008 0.38151984 -0.12566371)
Neighbour fc:(-0.061516383 -0.3815217 -1.5549985e-11)
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

The endwalls are geometrically the same, so I increased the matchTolerance to 0.01, but then I get the following error:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2112                                  |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : _14aeaf8dab-20211220 OPENFOAM=2112 version=v2112
Arch   : "LSB;label=32;scalar=64"
Exec   : checkMesh
Date   : Nov 30 2022
Time   : 01:19:14
Host   : 
PID    : 
I/O    : uncollated
Case   : 
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

Time = 0

Mesh stats
    points:           4543992
    faces:            13501728
    internal faces:   13372128
    cells:            4478976
    faces per cell:   6
    boundary patches: 5
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     4478976
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     0

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
    Patch               Faces    Points   Surface topology
    front               23328    23544    ok (non-closed singly connected)
    cylinder            41472    41688    ok (non-closed singly connected)
    back                23328    23544    ok (non-closed singly connected)
    inlet               20736    21037    ok (non-closed singly connected)
    outlet              20736    21037    ok (non-closed singly connected)

Checking faceZone topology for multiply connected surfaces...
    No faceZones found.

Checking basic cellZone addressing...
    No cellZones found.

Checking geometry...
    Overall domain bounding box (-0.4 -0.4 -0.12566371) (0.4 0.4 2.5839636e-10)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (8.5138466e-17 7.1580333e-18 -3.4146196e-15) OK.
    Max cell openness = 6.4200843e-16 OK.
    Max aspect ratio = 41.659087 OK.
    Minimum face area = 2.466335e-09. Maximum face area = 0.00031118574.  Face area magnitudes OK.
    Min volume = 1.618582e-12. Max volume = 2.0366275e-07.  Total volume = 0.063112991.  Cell volumes OK.
    Mesh non-orthogonality Max: 53.984103 average: 10.653729
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 1.5850591 OK.
  **Error in coupled point location: 14398 faces have their 0th or consecutive vertex not opposite their coupled equivalent. Average mismatch 0.00099244962.
  <<Writing 14398 faces with incorrectly matched 0th (or consecutive) vertex to set coupledFaces

Failed 1 mesh checks.

End

I have also tried cyclicAMI instead of cyclic, and it passes the checkMesh, but I get an error when running the case:

--> FOAM FATAL ERROR: (openfoam-2112)
Unable to set target face for source face 15178

    From virtual bool Foam::faceAreaWeightAMI::setNextFaces(Foam::label&, Foam::label&, Foam::label&, const Foam::bitSet&, Foam::labelList&, const Foam::DynamicList<int>&, bool) const
    in file AMIInterpolation/AMIInterpolation/faceAreaWeightAMI/faceAreaWeightAMI.C at line 363.

FOAM aborting

I have provided a link to the entire case file. Does anyone know what the issue is?

case file

The mesh is very fine, so I am running the case in parallel, but I haven't included this in the case file link, since I don't think this is affecting the cyclic BC issue.

Thanks!

[Swirl]

General question, did you use a createPatchDict and the createPatch command, to set up the polyMesh/boundary file correctly? because i couldn´t find the dict in your case files.