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

Cyclic BC for a twisted hexagonal cylinder

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 29, 2022, 20:39
Default Cyclic BC for a twisted hexagonal cylinder
  #1
New Member
 
Yicong
Join Date: Nov 2022
Posts: 2
Rep Power: 0
y.liu is on a distinguished road
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:

Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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:

Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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:

Code:
--> 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!
y.liu is offline   Reply With Quote

Old   December 26, 2022, 12:10
Default
  #2
Member
 
Join Date: Jun 2020
Posts: 49
Rep Power: 5
Swirl is on a distinguished road
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.
Swirl is offline   Reply With Quote

Reply

Tags
cyclic bc, twisted geometry

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
Flow pass 2D square cylinder using LES-channel 395 question Bashar OpenFOAM 1 January 8, 2021 01:53
How to implement time-varying boundary conditions in an enclosed cylinder Alejandro-FA OpenFOAM Pre-Processing 0 May 14, 2020 06:29
cyclic BC, cylinder array el_nando_97 OpenFOAM Running, Solving & CFD 1 February 25, 2020 02:34
flow over a cylinder urgent! kevin FLUENT 8 August 11, 2015 13:00
[blockMesh] Specifying boundary faces failes in blockMesh blaise OpenFOAM Meshing & Mesh Conversion 0 May 10, 2010 03:56


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