[Metti]

Dear all,

I'm trying to run an internal flow through a turbine cascade geometry.

I use "createPach" to change 4 of my BCs to "cyclicAMI" for the periodic BC.

besides, when i try to run using any solver, it crashes out from the begining like below.

Running icoFoam:

Reading transportProperties

Reading field p

AMI: Creating addressing and weights between 423 source faces and 476 target faces
AMI: Patch source weights min/max/average = 0.767896, 1.02016, 0.998634
AMI: Patch target weights min/max/average = 0, 1.03577, 0.978241
AMI: Creating addressing and weights between 808 source faces and 807 target faces
AMI: Patch source weights min/max/average = 0.954999, 1.03795, 1.00005
AMI: Patch target weights min/max/average = 0.0262439, 1.00738, 0.990712
Reading field U

Reading/calculating face flux field phi


Starting time loop

Time = 0.025

Courant Number mean: 1.39324e-05 max: 0.0057732
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 Uninterpreted:
#3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
#4
at ??:?
#5
at ??:?
#6 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#7
at ??:?
Floating point exception (core dumped)

boundary file:

10
(
Front
{
type empty;
nFaces 51859;
startFace 787374;
}
Back
{
type empty;
nFaces 51861;
startFace 839233;
}
Inlet
{
type wall;
nFaces 504;
startFace 891094;
}
Outlet
{
type wall;
nFaces 496;
startFace 891598;
}
P-Wall
{
type wall;
nFaces 2576;
startFace 892094;
}
S-Wall
{
type wall;
nFaces 2968;
startFace 894670;
}
ID-Cyclic1
{
type cyclicAMI;
inGroups 1(cyclicAMI);
nFaces 423;
startFace 897638;
matchTolerance 0.001;
transform translational;
neighbourPatch IU-Cyclic1;
separationVector (0 158 0);
}
IU-Cyclic1
{
type cyclicAMI;
inGroups 1(cyclicAMI);
nFaces 476;
startFace 898061;
matchTolerance 0.001;
transform translational;
neighbourPatch ID-Cyclic1;
separationVector (0 -158 0);
}
OD-Cyclic1
{
type cyclicAMI;
inGroups 1(cyclicAMI);
nFaces 808;
startFace 898537;
matchTolerance 0.001;
transform translational;
neighbourPatch OU-Cyclic1;
separationVector (0 158 0);
}
OU-Cyclic1
{
type cyclicAMI;
inGroups 1(cyclicAMI);
nFaces 807;
startFace 899345;
matchTolerance 0.001;
transform translational;
neighbourPatch OD-Cyclic1;
separationVector (0 -158 0);
}
)

checkMesh:

Mesh stats
points: 339811
faces: 900152
internal faces: 787374
cells: 280591
faces per cell: 6.01418
boundary patches: 10
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 273575
prisms: 3985
wedges: 0
pyramids: 0
tet wedges: 2
tetrahedra: 0
polyhedra: 3029
Breakdown of polyhedra by number of faces:
faces number of cells
4 111
5 10
6 768
7 556
8 9
9 819
10 10
12 537
15 200
18 9

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 51859 52629 ok (non-closed singly connected)
Back 51861 52633 ok (non-closed singly connected)
Inlet_solid 504 635 ok (non-closed singly connected)
Outlet_solid 496 625 ok (non-closed singly connected)
P-Wall_solid 2576 2929 ok (non-closed singly connected)
S-Wall_solid 2968 3383 ok (non-closed singly connected)
ID-Cyclic1 423 529 ok (non-closed singly connected)
IU-Cyclic1 476 584 ok (non-closed singly connected)
OD-Cyclic1 808 1019 ok (non-closed singly connected)
OU-Cyclic1 807 1030 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-105 -308.338 -0.0270165) (275.001 173.067 5.03104)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (3.17921e-17 5.13671e-18 4.65181e-15) OK.
Max cell openness = 3.0873e-16 OK.
Max aspect ratio = 9.84905 OK.
Minimum face area = 0.0156956. Maximum face area = 7.20545. Face area magnitudes OK.
Min volume = 0.014621. Max volume = 3.59492. Total volume = 282960. Cell volumes OK.
Mesh non-orthogonality Max: 84.028 average: 4.41006
*Number of severely non-orthogonal (> 70 degrees) faces: 1.
Non-orthogonality check OK.
<<Writing 1 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
Max skewness = 2.46165 OK.
Coupled point location match (average 0) OK.

Mesh OK.

End

whether without changing the BCs to "cyclicAMI" it works fine.
I'm really messed up of where does the error originate from, but it seems that it does not belong to my mesh. since the same mesh works fine without "cyclicAMI".(correct me if I'm out of picture)

my steps are:
1. blockMesh
2. surfaceFeatureExtract
3. snappyHexMesh
4. extrudeMesh
5. createPatch
6. run

any suggestions..?

best regards
Metti,