Can I get the geometry for rotor 37 compressor?

[Mohsin-Bashir]

Could you please provide the geometry used to create mesh in tutorial rotor 37 compressor?

Also, where can I get more information about this test case, is there any paper about it?

[friedenhe]

We generated the mesh using TurboGrid. So we got the geometry from TurboGrid's tutorial. We did not save the geometry....

[Mohsin-Bashir]

In Turbogrid tutorial it is mentioned that there is a gap of 2.5% of the total span between blade and shroud. But the mesh given on DAFoam website does not seem to have any gap. Could you please confirm this?

[friedenhe]

Right. We did not add the tip clearance for the R37 tutorial case..

[Mohsin-Bashir]

We generated the mesh using TurboGrid. So we got the geometry from TurboGrid's tutorial. We did not save the geometry....

TurboGrid does not support exporting mesh as .msh. Could you please let me know how you converted the mesh to fluent3d .msh format. I have only two options cgns and ansys cfx mesh files (.gtm .def)

In turboGrid I exported mesh as cgns format then using gmsh saved it to gmsh .msh format. Then I used this command gmshToFoam rotor37.msh to convert it to openfoam format. But the problem is when I run checkMesh my mesh fails 6 mesh checks.

dafoamuser@68e7dbf17200:~/mount/7_dafoam$ checkMesh
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1812                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v1812 OPENFOAM=1812
Arch   : "LSB;label=32;scalar=64"
Exec   : checkMesh
Date   : Jun 19 2024
Time   : 03:05:13
Host   : 68e7dbf17200
PID    : 384
I/O    : uncollated
Case   : /home/dafoamuser/mount/7_dafoam
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Overriding DebugSwitches according to controlDict
    SolverPerformance 0;

Create mesh for time = 0

Time = 0

Mesh stats
    points:           166152
    faces:            475313
    internal faces:   452785
    cells:            154683
    faces per cell:   6
    boundary patches: 7
    point zones:      0
    face zones:       0
    cell zones:       1

Overall number of cells of each type:
    hexahedra:     154683
    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
    hub                 5729     5934     ok (non-closed singly connected)
    shroud              5729     5934     ok (non-closed singly connected)
    Inlet               1053     1120     ok (non-closed singly connected)
    outlet              972      1036     ok (non-closed singly connected)
    blade               4563     4732     ok (non-closed singly connected)
    per1                2241     2352     ok (non-closed singly connected)
    per2                2241     2352     ok (non-closed singly connected)

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

Checking basic cellZone addressing...
    CellZone            Cells        Points       Volume       BoundingBox
    volume              154683       166152       0.0001316654796281227 (0.1763780449527726 -0.01016055771465312 -0.009490815155249004) (0.2522999873927165 0.06228013791976354 0.05798066736143725)

Checking geometry...
    Overall domain bounding box (0.1763780449527726 -0.01016055771465312 -0.009490815155249004) (0.2522999873927165 0.06228013791976354 0.05798066736143725)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
 ***Boundary openness (0.008946865396841775 -0.04427588224865455 0.05411739295007382) possible hole in boundary description.
 ***Open cells found, max cell openness: 1, number of open cells 351
  <<Writing 351 non closed cells to set nonClosedCells
  <<Writing 227 cells with high aspect ratio to set highAspectRatioCells
    Minimum face area = 4.482590898700458e-10. Maximum face area = 1.879119529075913e-05.  Face area magnitudes OK. ***Zero or negative cell volume detected.  Minimum negative volume: -3.084917671882813e-10, Number of negative volume cells: 227
  <<Writing 227 zero volume cells to set zeroVolumeCells
    Mesh non-orthogonality Max: 177.7254724844629 average: 11.96743086930206
   *Number of severely non-orthogonal (> 70 degrees) faces: 277.
 ***Number of non-orthogonality errors: 511.
  <<Writing 788 non-orthogonal faces to set nonOrthoFaces
 ***Error in face pyramids: 1210 faces are incorrectly oriented.
  <<Writing 707 faces with incorrect orientation to set wrongOrientedFaces
 ***Max skewness = 103.7997219375421, 14 highly skew faces detected which may impair the quality of the results
  <<Writing 14 skew faces to set skewFaces
    Coupled point location match (average 0) OK.

Failed 6 mesh checks.

End

[friedenhe]

What we did was saving the mesh in cgns format. Then, we used ICEM-CFD to read the cgns mesh and export it as ansys-fluent format. Then, we use fluent3DMeshToFoam to convert it further to OpennFOAM format.