[sam.ho]

HI,

I am simulating a radial turbine impeller simulation using openfoam.
I have periodic boundaries as i want a passage to be simulated.
createPatchDict looks like

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      createPatchDict;
}

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

matchTolerance 1E-3;

pointSync true;

patchInfo
(
    {
        name ROTOR_PERIODIC;
        dictionary
        {
            type cyclic;
        }
        constructFrom patches;
        patches ( PERIODIC1 PERIODIC2);
    }

);

When i run createPatch command ..
I will get the following error

sml5kor@BMH301562:~/OpenFOAM/sml5kor-1.6-ext/run/Only_IMPELLER$ createPatch 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM Extend Project: Open source CFD        |
|  \\    /   O peration     | Version:  1.6-ext                               |
|   \\  /    A nd           | Web:      www.extend-project.de                 |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 1.6-ext
Exec   : createPatch
Date   : Apr 14 2014
Time   : 17:17:11
Host   : BMH301562
PID    : 24473
Case   : /home/sml5kor/OpenFOAM/sml5kor-1.6-ext/run/Only_IMPELLER
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

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

Reading createPatchDict.

Create polyMesh for time = 0

Adding new patch ROTOR_PERIODIC as patch 8 from 
{
    type            cyclic;
}


--> FOAM Warning : 
    From function directTopoChange::addMesh(const polyMesh&, const labelList&,const labelList&, const labelList&,const labelList&)
    in file directTopoChange/directTopoChange/directTopoChange.C at line 2281
Segmentation fault (core dumped)

What could be the problem ?

Regards,
Sangamesh Hosur

[sivakumar]

Hello Hosur Boy,
What could be the result from the current Parliament election?

coming to your problem, your problem description is not clear,
I got these kind of problem, in the beginning.

It could be because of your physical memory.

how many GB RAM you are using?
what is the size of your mesh (how many million cells)?

generally 1 million cells needs 1GB ram,

if you are still in problem, please explain your problem clearly.

Have a nice evening,

Sivakumar

[sam.ho]

Hi Siva

Assembly Election results are yet to be announced..

I have 0.9 million cells and my system consists of 2GB RAM.

If i keep pointSync true ;
then the output is

sml5kor@BMH301562:~/OpenFOAM/sml5kor-1.6-ext/run/Turbine_CASE_16.04.2014$ createPatch 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM Extend Project: Open source CFD        |
|  \\    /   O peration     | Version:  1.6-ext                               |
|   \\  /    A nd           | Web:      www.extend-project.de                 |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 1.6-ext
Exec   : createPatch
Date   : Apr 21 2014
Time   : 11:26:16
Host   : BMH301562
PID    : 27877
Case   : /home/sml5kor/OpenFOAM/sml5kor-1.6-ext/run/Turbine_CASE_16.04.2014
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

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

Reading createPatchDict.

Create polyMesh for time = 0

Adding new patch ROTOR_PERIODIC1 as patch 20 from 
{
    type            cyclic;
    transform       rotational;
    rotationAxis    ( 0 0 1 );
    rotationCentre  ( 0 0 0 );
    rotationAngle   32.7272727272727;
}

Adding new patch UPS_PERIODIC as patch 21 from 
{
    type            cyclic;
    transform       rotational;
    rotationAxis    ( 0 0 1 );
    rotationCentre  ( 0 0 0 );
    rotationAngle   32.7272727272727;
}

Adding new patch DWS_PERIODIC as patch 22 from 
{
    type            cyclic;
    transform       rotational;
    rotationAxis    ( 0 0 1 );
    rotationCentre  ( 0 0 0 );
    rotationAngle   32.7272727272727;
}


Moving faces from patch PERIODIC1 to patch 20
Moving faces from patch PERIODIC2 to patch 20
Moving faces from patch PER1_UPS to patch 21
Moving faces from patch PER2_UPS to patch 21
Moving faces from patch PER1_DWS to patch 22
Moving faces from patch PER2_DWS to patch 22

Doing topology modification to order faces.

cyclicPolyPatch::order : Writing half0 faces to OBJ file "UPS_PERIODIC_half0_faces.obj"
cyclicPolyPatch::order : Writing half1 faces to OBJ file "UPS_PERIODIC_half1_faces.obj"
cyclicPolyPatch::order : Dumping currently found cyclic match as lines between corresponding face centres to file "/home/sml5kor/OpenFOAM/sml5kor-1.6-ext/run/Turbine_CASE_16.04.2014/UPS_PERIODIC_faceCentres.obj"
--> FOAM Serious Error : 
    From function cyclicPolyPatch::order(const primitivePatch&, labelList&, labelList&) const
    in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 1523
    Patch UPS_PERIODIC : Cannot match vectors to faces on both sides of patch UPS_PERIODIC
.  Perhaps your faces do not match? The obj files written contain the current match.
    Continuing with incorrect face ordering from now on!
cyclicPolyPatch::order : Writing half0 faces to OBJ file "DWS_PERIODIC_half0_faces.obj"
cyclicPolyPatch::order : Writing half1 faces to OBJ file "DWS_PERIODIC_half1_faces.obj"
cyclicPolyPatch::order : Dumping currently found cyclic match as lines between corresponding face centres to file "/home/sml5kor/OpenFOAM/sml5kor-1.6-ext/run/Turbine_CASE_16.04.2014/DWS_PERIODIC_faceCentres.obj"
--> FOAM Serious Error : 
    From function cyclicPolyPatch::order(const primitivePatch&, labelList&, labelList&) const
    in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 1523
    Patch DWS_PERIODIC : Cannot match vectors to faces on both sides of patch DWS_PERIODIC
.  Perhaps your faces do not match? The obj files written contain the current match.
    Continuing with incorrect face ordering from now on!

    From function void polyMesh::initMesh()
    in file meshes/polyMesh/polyMeshInitMesh.C at line 82
    Truncating neighbour list at 1813845 for backward compatibility
Dumping ROTOR_PERIODIC1 half0 faces to "coupled_ROTOR_PERIODIC1_half0.obj"
Dumping ROTOR_PERIODIC1 half1 faces to "coupled_ROTOR_PERIODIC1_half1.obj"
Dumping cyclic match as lines between face centres to "coupled_ROTOR_PERIODIC1_match.obj"
Dumping UPS_PERIODIC half0 faces to "coupled_UPS_PERIODIC_half0.obj"
Dumping UPS_PERIODIC half1 faces to "coupled_UPS_PERIODIC_half1.obj"
Dumping cyclic match as lines between face centres to "coupled_UPS_PERIODIC_match.obj"
Dumping DWS_PERIODIC half0 faces to "coupled_DWS_PERIODIC_half0.obj"
Dumping DWS_PERIODIC half1 faces to "coupled_DWS_PERIODIC_half1.obj"
Dumping cyclic match as lines between face centres to "coupled_DWS_PERIODIC_match.obj"
Synchronising points.

On coupled patch ROTOR_PERIODIC1 forcing uniform rotation of (0.841253532831181 -0.540640817455598 0 0.540640817455598 0.841253532831181 0 0 0 1)
On coupled patch ROTOR_PERIODIC1 forcing uniform rotation of 1((0.841253532831181 -0.540640817455598 0 0.540640817455598 0.841253532831181 0 0 0 1))
On coupled patch UPS_PERIODIC forcing uniform rotation of (0.841253532831181 -0.540640817455598 0 0.540640817455598 0.841253532831181 0 0 0 1)
On coupled patch UPS_PERIODIC forcing uniform rotation of 1((0.841253532831181 -0.540640817455598 0 0.540640817455598 0.841253532831181 0 0 0 1))
On coupled patch DWS_PERIODIC forcing uniform rotation of (0.841253532831181 -0.540640817455598 0 0.540640817455598 0.841253532831181 0 0 0 1)
On coupled patch DWS_PERIODIC forcing uniform rotation of 1((0.841253532831181 -0.540640817455598 0 0.540640817455598 0.841253532831181 0 0 0 1))
Synchronising points.
--> FOAM Warning : 
    From function syncPoints(const polyMesh&, pointField&, const CombineOp&, const point&)
    in file createPatch.C at line 484
    There are decomposed cyclics in this mesh with transformations.
This is not supported. The result will be incorrect
Points changed by average:0.180625681331992 max:26.2926718631749



--> FOAM FATAL ERROR: 
face 1301 area does not match neighbour 2603 by 6.39394% -- possible face ordering problem.
patch:UPS_PERIODIC my area:0.000429092 neighbour area:0.000402506 matching tolerance:0.0001
Mesh face:1878732 vertices:4((28.9613 8.5167 0.0772216) (28.9613 8.5167 0.0685689) (28.9137 8.50274 0.0685642) (28.9137 8.50275 0.0772158))
Neighbour face:1880034 vertices:4((28.9639 -8.50583 0.0685689) (28.9639 -8.50583 0.0772216) (28.9683 -8.49296 0.0772216) (28.9163 -8.49183 0.0685642))
Other errors also exist, only the largest is reported. Please rerun with cyclic debug flag set for more information.

    From function cyclicPolyPatch::calcTransforms()
    in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 294.

FOAM aborting

Aborted (core dumped)

If i Keep pointSync false;
then the output is

sml5kor@BMH301562:~/OpenFOAM/sml5kor-1.6-ext/run/Turbine_CASE_16.04.2014$ createPatch 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM Extend Project: Open source CFD        |
|  \\    /   O peration     | Version:  1.6-ext                               |
|   \\  /    A nd           | Web:      www.extend-project.de                 |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 1.6-ext
Exec   : createPatch
Date   : Apr 21 2014
Time   : 11:28:48
Host   : BMH301562
PID    : 27903
Case   : /home/sml5kor/OpenFOAM/sml5kor-1.6-ext/run/Turbine_CASE_16.04.2014
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

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

Reading createPatchDict.

Create polyMesh for time = 0

Adding new patch ROTOR_PERIODIC1 as patch 20 from 
{
    type            cyclic;
    transform       rotational;
    rotationAxis    ( 0 0 1 );
    rotationCentre  ( 0 0 0 );
    rotationAngle   32.7272727272727;
}

Adding new patch UPS_PERIODIC as patch 21 from 
{
    type            cyclic;
    transform       rotational;
    rotationAxis    ( 0 0 1 );
    rotationCentre  ( 0 0 0 );
    rotationAngle   32.7272727272727;
}

Adding new patch DWS_PERIODIC as patch 22 from 
{
    type            cyclic;
    transform       rotational;
    rotationAxis    ( 0 0 1 );
    rotationCentre  ( 0 0 0 );
    rotationAngle   32.7272727272727;
}


Moving faces from patch PERIODIC1 to patch 20
Moving faces from patch PERIODIC2 to patch 20
Moving faces from patch PER1_UPS to patch 21
Moving faces from patch PER2_UPS to patch 21
Moving faces from patch PER1_DWS to patch 22
Moving faces from patch PER2_DWS to patch 22

Doing topology modification to order faces.

cyclicPolyPatch::order : Writing half0 faces to OBJ file "UPS_PERIODIC_half0_faces.obj"
cyclicPolyPatch::order : Writing half1 faces to OBJ file "UPS_PERIODIC_half1_faces.obj"
cyclicPolyPatch::order : Dumping currently found cyclic match as lines between corresponding face centres to file "/home/sml5kor/OpenFOAM/sml5kor-1.6-ext/run/Turbine_CASE_16.04.2014/UPS_PERIODIC_faceCentres.obj"
--> FOAM Serious Error : 
    From function cyclicPolyPatch::order(const primitivePatch&, labelList&, labelList&) const
    in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 1523
    Patch UPS_PERIODIC : Cannot match vectors to faces on both sides of patch UPS_PERIODIC
.  Perhaps your faces do not match? The obj files written contain the current match.
    Continuing with incorrect face ordering from now on!
cyclicPolyPatch::order : Writing half0 faces to OBJ file "DWS_PERIODIC_half0_faces.obj"
cyclicPolyPatch::order : Writing half1 faces to OBJ file "DWS_PERIODIC_half1_faces.obj"
cyclicPolyPatch::order : Dumping currently found cyclic match as lines between corresponding face centres to file "/home/sml5kor/OpenFOAM/sml5kor-1.6-ext/run/Turbine_CASE_16.04.2014/DWS_PERIODIC_faceCentres.obj"
--> FOAM Serious Error : 
    From function cyclicPolyPatch::order(const primitivePatch&, labelList&, labelList&) const
    in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 1523
    Patch DWS_PERIODIC : Cannot match vectors to faces on both sides of patch DWS_PERIODIC
.  Perhaps your faces do not match? The obj files written contain the current match.
    Continuing with incorrect face ordering from now on!

    From function void polyMesh::initMesh()
    in file meshes/polyMesh/polyMeshInitMesh.C at line 82
    Truncating neighbour list at 1813845 for backward compatibility
Dumping ROTOR_PERIODIC1 half0 faces to "coupled_ROTOR_PERIODIC1_half0.obj"
Dumping ROTOR_PERIODIC1 half1 faces to "coupled_ROTOR_PERIODIC1_half1.obj"
Dumping cyclic match as lines between face centres to "coupled_ROTOR_PERIODIC1_match.obj"
Dumping UPS_PERIODIC half0 faces to "coupled_UPS_PERIODIC_half0.obj"
Dumping UPS_PERIODIC half1 faces to "coupled_UPS_PERIODIC_half1.obj"
Dumping cyclic match as lines between face centres to "coupled_UPS_PERIODIC_match.obj"
Dumping DWS_PERIODIC half0 faces to "coupled_DWS_PERIODIC_half0.obj"
Dumping DWS_PERIODIC half1 faces to "coupled_DWS_PERIODIC_half1.obj"
Dumping cyclic match as lines between face centres to "coupled_DWS_PERIODIC_match.obj"
Not synchronising points.

Removing patches with no faces in them.

Removing empty patch PER1_UPS at position 3
Removing empty patch PER2_UPS at position 4
Removing empty patch PERIODIC1 at position 6
Removing empty patch PERIODIC2 at position 7
Removing empty patch PER1_DWS at position 18
Removing empty patch PER2_DWS at position 19
Removing patches.
Dumping ROTOR_PERIODIC1 half0 faces to "final_ROTOR_PERIODIC1_half0.obj"
Dumping ROTOR_PERIODIC1 half1 faces to "final_ROTOR_PERIODIC1_half1.obj"
Dumping cyclic match as lines between face centres to "final_ROTOR_PERIODIC1_match.obj"
Dumping UPS_PERIODIC half0 faces to "final_UPS_PERIODIC_half0.obj"
Dumping UPS_PERIODIC half1 faces to "final_UPS_PERIODIC_half1.obj"
Dumping cyclic match as lines between face centres to "final_UPS_PERIODIC_match.obj"
Dumping DWS_PERIODIC half0 faces to "final_DWS_PERIODIC_half0.obj"
Dumping DWS_PERIODIC half1 faces to "final_DWS_PERIODIC_half1.obj"
Dumping cyclic match as lines between face centres to "final_DWS_PERIODIC_match.obj"
Writing repatched mesh to 1

End

Is that necessary to keep pointSync true ?
How i need to take a call on this ?