[gisc]

Hello there,
I'm using topoSet to define a cellSet in which I'd like to apply trip later on.


I've created a 2D mesh in ICEM, so I run

fluentMeshToFoam

and subsequently the command

topoSet

.

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  21.06                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      faceSetDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

actions
(
    {
        name    tripSet;
        action  new;
        type    cellSet;
        source  rotatedBoxToCell;
        sourceInfo
        {
                    origin      (0.09977 0.05577 -0.5);
                    i            (4.93029e-3 1.09302e-3 0);
                    j            (-1.29864e-4 5.85778e-4 0);
                    k           (0 0 0.5);
        }
    }
);

I have defined the tripping box from Paraview, and I'm expecting to find 50 cells after the topoSet command (5 cells high, 10 cells length, second image for reference, red box).


The problem is that topoSet just select some random cells and not all the cells inside the rectangle. The same code works for the LES 3d case (which is the second image), but when I apply it to the 2D RANS, it gives me this random output (first image).


I don't understand what I'm missing. At the beginning I thought that it was a problem due to the definiton of RANS's patch

frontAndBackPlanes  empty

, therefore I splitted it in

left

and

right

, and defined them as

cyclic

. However, still I got the same random output.


If anybody might have a clue or a hint to give it to me, it'd be very appreciated. Below you can find the

checkMesh

output for reference.

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2106                                  |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Time = 0

Mesh stats
    points:           942026
    internal points:  0
    faces:            1874296
    internal faces:   932270
    cells:            467761
    faces per cell:   6
    boundary patches: 6
    point zones:      0
    face zones:       0
    cell zones:       1

Overall number of cells of each type:
    hexahedra:     467761
    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
    OUTLET              266      534      ok (non-closed singly connected)
    INLET               3517     7036     ok (non-closed singly connected)
    SS                  1361     2724     ok (non-closed singly connected)
    PS                  1360     2722     ok (non-closed singly connected)
    front               467761   471013   ok (non-closed singly connected)
    back                467761   471013   ok (non-closed singly connected)

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

Checking basic cellZone addressing...
    CellZone            Cells        Points       Volume       BoundingBox
    tripSet             34           128          6.40191e-07  (0.0997363 0.0558234 -0.180277) (0.104675 0.0573815 0.180277)

Checking geometry...
    Overall domain bounding box (-4.99999 -5 -0.180277) (10 5 0.180277)
    Mesh has 2 geometric (non-empty/wedge) directions (1 1 0)
    Mesh has 2 solution (non-empty) directions (1 1 0)
    All edges aligned with or perpendicular to non-empty directions.
    Boundary openness (-1.24833e-17 2.34241e-17 -8.11459e-15) OK.
    Max cell openness = 4.27054e-15 OK.
    Max aspect ratio = 83.5961 OK.
    Minimum face area = 1.49745e-08. Maximum face area = 0.0510004.  Face area magnitudes OK.
    Min volume = 5.39913e-09. Max volume = 0.00176845.  Total volume = 50.1848.  Cell volumes OK.
    Mesh non-orthogonality Max: 62.4138 average: 22.3492
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 2.11943 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End