[blockMesh] Issue with mergePatchPairs

JohnnyMercu · May 11, 2021, 09:08

Hi everyone,
As a beginner with OpenFOAM, I wanted to realize a mesh composed of 9 different zones with different mesh refinements. To avoid the necessity of continuous mesh between the zones, I've read I could use mergePatchPairs. The thing is I am not sure I am coding it right and I don't really know which zone should I define as Master or as Slave. Here is the blockMeshDict I use and the error message :

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  7                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#include "../../include/initialConditions"

scale   1;

vertices
(

    //Block 0
    ( -60   6.7      -60)
    (-7.2   6.7      -60)
    (-7.2   6.8      -60)
    ( -60   6.8      -60)
    ( -60   6.7     -7.2)
    (-7.2   6.7     -7.2)
    (-7.2   6.8     -7.2)
    ( -60   6.8     -7.2)

    //Block 1
    (-7.2   6.7      -60)
    ( 7.2   6.7      -60)
    ( 7.2   6.8      -60)
    (-7.2   6.8      -60)
    (-7.2   6.7     -7.2)
    ( 7.2   6.7     -7.2)
    ( 7.2   6.8     -7.2)
    (-7.2   6.8     -7.2)

    //Block 2

    ( 7.2   6.7      -60)
    ( 120   6.7      -60)
    ( 120   6.8      -60)
    ( 7.2   6.8      -60)
    ( 7.2   6.7     -7.2)
    ( 120   6.7     -7.2)
    ( 120   6.8     -7.2)
    ( 7.2   6.8     -7.2)

    //Block 3

    ( 7.2   6.7     -7.2)
    ( 120   6.7     -7.2)
    ( 120   6.8     -7.2)
    ( 7.2   6.8     -7.2)
    ( 7.2   6.7      7.2)
    ( 120   6.7      7.2)
    ( 120   6.8      7.2)
    ( 7.2   6.8      7.2)

    //Block 4

    (-7.2   6.7     -7.2)
    ( 7.2   6.7     -7.2)
    ( 7.2   6.8     -7.2)
    (-7.2   6.8     -7.2)
    (-7.2   6.7      7.2)
    ( 7.2   6.7      7.2)
    ( 7.2   6.8      7.2)
    (-7.2   6.8      7.2)

    //Block 5

    ( -60   6.7     -7.2)
    (-7.2   6.7     -7.2)
    (-7.2   6.8     -7.2)
    ( -60   6.8     -7.2)
    ( -60   6.7      7.2)
    (-7.2   6.7      7.2)
    (-7.2   6.8      7.2)
    ( -60   6.8      7.2)

    //Block 6

    ( -60   6.7      7.2)
    (-7.2   6.7      7.2)
    (-7.2   6.8      7.2)
    ( -60   6.8      7.2)
    ( -60   6.7       60)
    (-7.2   6.7       60)
    (-7.2   6.8       60)
    ( -60   6.8       60)

    //Block 7

    (-7.2   6.7      7.2)
    ( 7.2   6.7      7.2)
    ( 7.2   6.8      7.2)
    (-7.2   6.8      7.2)
    (-7.2   6.7       60)
    ( 7.2   6.7       60)
    ( 7.2   6.8       60)
    (-7.2   6.8       60)

    //Block 8

    ( 7.2   6.7      7.2)
    ( 120   6.7      7.2)
    ( 120   6.8      7.2)
    ( 7.2   6.8      7.2)
    ( 7.2   6.7       60)
    ( 120   6.7       60)
    ( 120   6.8       60)
    ( 7.2   6.8       60)



/*
    ($min_box  1.7  $min_box)
    ($max_box  1.7  $min_box)
    ($max_box  1.8  $min_box)
    ($min_box  1.8  $min_box)
    ($min_box  1.7  $max_box)
    ($max_box  1.7  $max_box)
    ($max_box  1.8  $max_box)
    ($min_box  1.8  $max_box)
    */
);

blocks
(
    hex (0 1 2 3 4 5 6 7)           (10 1 10)     simpleGrading (1 1 1)
    hex (8 9 10 11 12 13 14 15)     (29 1 10)     simpleGrading (1 1 1)
    hex (16 17 18 19 20 21 22 23)   (10 1 10)     simpleGrading (1 1 1)
    hex (24 25 26 27 28 29 30 31)   (10 1 29)     simpleGrading (1 1 1)
    hex (32 33 34 35 36 37 38 39)   (29 1 29)     simpleGrading (1 1 1)
    hex (40 41 42 43 44 45 46 47)   (10 1 29)     simpleGrading (1 1 1)
    hex (48 49 50 51 52 53 54 55)   (10 1 10)     simpleGrading (1 1 1)
    hex (56 57 58 59 60 61 62 63)   (29 1 10)     simpleGrading (1 1 1)
    hex (64 65 66 67 68 69 70 71)   (10 1 10)     simpleGrading (1 1 1)
);

edges
(
);

boundary
(   
    topAndBottom
    {
        type patch;
        faces
        (
            (52 53 54 55)
            (60 61 62 63)
            (68 69 70 71)
            (0 1 2 3)
            (8 9 10 11)
            (16 17 18 19)
        );
    }

    inlet
    {
        type patch;
        faces
        (            
            (0 4 7 3)
            (40 44 47 43)
            (48 52 55 51)
        );
    }

    outlet
    {
        type patch;
        faces
        (
            (17 21 22 18)
            (25 29 30 26)
            (65 69 70 66)

        );
    }

    front
    {
        type empty;
        faces
        (
            (3 2 6 7)
            (11 10 14 15)
            (19 18 22 23)
            (27 26 30 31)
            (35 34 38 39)
            (43 42 46 47)
            (51 50 54 55)
            (59 58 62 63)
            (67 66 70 71)
        );
    }

    back
    {
        type empty;
        faces
        (
            (0 1 5 4)
            (8 9 13 12)
            (16 17 21 20)
            (24 25 29 28)
            (32 33 37 36)
            (40 41 45 44)
            (48 49 53 52)
            (56 57 61 60)
            (64 65 69 68)
        );
    }

    face05
    {
        type patch;
        faces
        (
            (4 5 6 7)
        );
    }

    face50
    {
        type patch;
        faces
        (
            (40 41 42 43)
        );
    }

    face01
    {
        type patch;
        faces
        (
            (1 5 6 2)
        );
    }

    face10
    {
        type patch;
        faces
        (
            (8 12 15 11)
        );
    }

    face14
    {
        type patch;
        faces
        (
            (12 13 14 15)
        );
    }

    face41
    {
        type patch;
        faces
        (
            (32 33 34 35)
        );
    }

    face12
    {
        type patch;
        faces
        (
            (9 13 14 10)
        );
    }

    face21
    {
        type patch;
        faces
        (
            (16 20 23 19)
        );
    }

    face23
    {
        type patch;
        faces
        (
            (20 21 22 23)
        );
    }

    face32
    {
        type patch;
        faces
        (
            (24 25 26 27)
        );
    }

    face34
    {
        type patch;
        faces
        (
            (24 28 31 27)
        );
    }

    face43
    {
        type patch;
        faces
        (
            (33 37 38 34)
        );
    }

    face45
    {
        type patch;
        faces
        (
            (32 36 39 35)
        );
    }

    face54
    {
        type patch;
        faces
        (
            (41 45 46 42)
        );
    }

    face56
    {
        type patch;
        faces
        (
            (44 45 46 47)
        );
    }

    face65
    {
        type patch;
        faces
        (
            (48 49 50 51)
        );
    }

    face67
    {
        type patch;
        faces
        (
            (49 53 54 50)
        );
    }

    face76
    {
        type patch;
        faces
        (
            (56 60 63 59)
        );
    }

    face78
    {
        type patch;
        faces
        (
            (57 61 62 58)
        );
    }

    face87
    {
        type patch;
        faces
        (
            (64 68 71 67)
        );
    }

    face47
    {
        type patch;
        faces
        (
            (36 37 38 39)
        );
    }

    face74
    {
        type patch;
        faces
        (
            (56 57 58 59)
        );
    }

    face38
    {
        type patch;
        faces
        (
            (28 29 30 31)
        );
    }

    face83
    {
        type patch;
        faces
        (
            (64 65 66 67)
        );
    }
   
);

mergePatchPairs
(
    (face01 face10)
    (face12 face21)
    (face23 face32)
    (face41 face14)
    (face05 face50)
    (face45 face54)
    (face43 face34)
    (face56 face65)
    (face47 face74)
    (face38 face83)
    (face67 face76)
    (face78 face87)
);

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

Code:

--> FOAM FATAL ERROR: (openfoam-2012)
Face 7504 reduced to less than 3 points.  Topological/cutting error A.
Old face: 2(1101 1123) new face: 2(1101 1123)

    From void Foam::slidingInterface::coupleInterface(Foam::polyTopoChange&) const
    in file slidingInterface/coupleSlidingInterface.C at line 1450.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::exitOrAbort(int, bool) at ??:?
#2  Foam::slidingInterface::coupleInterface(Foam::polyTopoChange&) const at ??:?
#3  Foam::slidingInterface::setRefinement(Foam::polyTopoChange&) const at ??:?
#4  Foam::polyTopoChanger::topoChangeRequest() const at ??:?
#5  Foam::polyTopoChanger::changeMesh(bool, bool, bool, bool) at ??:?
#6  Foam::attachPolyTopoChanger::attach(bool) at ??:?
#7  ? at ??:?
#8  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#9  ? at ??:?
./run_full: line 102: 17686 Aborted                 blockMesh > log.blockMesh


--> FOAM FATAL ERROR: (openfoam-2012)
Cannot find file "points" in directory "polyMesh" in times "0" down to constant

    From virtual Foam::IOobject Foam::fileOperation::findInstance(const Foam::IOobject&, Foam::scalar, const Foam::word&) const
    in file global/fileOperations/fileOperation/fileOperation.C at line 1003.

FOAM exiting



--> FOAM FATAL ERROR: (openfoam-2012)
Cannot find file "points" in directory "polyMesh" in times "0" down to constant

    From virtual Foam::IOobject Foam::fileOperation::findInstance(const Foam::IOobject&, Foam::scalar, const Foam::word&) const
    in file global/fileOperations/fileOperation/fileOperation.C at line 1003.

FOAM exiting



--> FOAM FATAL ERROR: (openfoam-2012)
Cannot find file "points" in directory "polyMesh" in times "0" down to constant

    From virtual Foam::IOobject Foam::fileOperation::findInstance(const Foam::IOobject&, Foam::scalar, const Foam::word&) const
    in file global/fileOperations/fileOperation/fileOperation.C at line 1003.

FOAM exiting



--> FOAM FATAL ERROR: (openfoam-2012)
Cannot find file "points" in directory "polyMesh" in times "0" down to constant

    From virtual Foam::IOobject Foam::fileOperation::findInstance(const Foam::IOobject&, Foam::scalar, const Foam::word&) const
    in file global/fileOperations/fileOperation/fileOperation.C at line 1003.

FOAM exiting



--> FOAM FATAL ERROR: (openfoam-2012)
Cannot find file "points" in directory "polyMesh" in times "0" down to constant

    From virtual Foam::IOobject Foam::fileOperation::findInstance(const Foam::IOobject&, Foam::scalar, const Foam::word&) const
    in file global/fileOperations/fileOperation/fileOperation.C at line 1003.

FOAM exiting



--> FOAM FATAL ERROR: (openfoam-2012)
Cannot find file "points" in directory "polyMesh" in times "0" down to constant

    From virtual Foam::IOobject Foam::fileOperation::findInstance(const Foam::IOobject&, Foam::scalar, const Foam::word&) const
    in file global/fileOperations/fileOperation/fileOperation.C at line 1003.

FOAM exiting



--> FOAM FATAL ERROR: (openfoam-2012)
Cannot find file "points" in directory "polyMesh" in times "0" down to constant

    From virtual Foam::IOobject Foam::fileOperation::findInstance(const Foam::IOobject&, Foam::scalar, const Foam::word&) const
    in file global/fileOperations/fileOperation/fileOperation.C at line 1003.

FOAM exiting



--> FOAM FATAL ERROR: (openfoam-2012)
Cannot find file "points" in directory "polyMesh" in times "0" down to constant

    From virtual Foam::IOobject Foam::fileOperation::findInstance(const Foam::IOobject&, Foam::scalar, const Foam::word&) const
    in file global/fileOperations/fileOperation/fileOperation.C at line 1003.

FOAM exiting



--> FOAM FATAL ERROR: (openfoam-2012)
Cannot find file "points" in directory "polyMesh" in times "0" down to constant

    From virtual Foam::IOobject Foam::fileOperation::findInstance(const Foam::IOobject&, Foam::scalar, const Foam::word&) const
    in file global/fileOperations/fileOperation/fileOperation.C at line 1003.

FOAM exiting



--> FOAM FATAL ERROR: (openfoam-2012)
Cannot find file "points" in directory "polyMesh" in times "0" down to constant

    From virtual Foam::IOobject Foam::fileOperation::findInstance(const Foam::IOobject&, Foam::scalar, const Foam::word&) const
    in file global/fileOperations/fileOperation/fileOperation.C at line 1003.

FOAM exiting



--> FOAM FATAL ERROR: (openfoam-2012)
Cannot find file "points" in directory "polyMesh" in times "0" down to constant

    From virtual Foam::IOobject Foam::fileOperation::findInstance(const Foam::IOobject&, Foam::scalar, const Foam::word&) const
    in file global/fileOperations/fileOperation/fileOperation.C at line 1003.

FOAM exiting

I tried to find a solution, and I am sorry if the subject has already been solved before.

If anyone could help me I would be very grateful !

anaspauzi · May 30, 2021, 01:24

Same, I'm having the same problem and found the solution.

Quote:

--> FOAM FATAL ERROR: (openfoam-2012)
Face XXX reduced to less than 3 points. Topological/cutting error A.
Old face: 2(XXX XXX) new face: 2(XXX XXX)

From void Foam::slidingInterface::coupleInterface(Foam:

oly TopoChange&) const
in file slidingInterface/coupleSlidingInterface.C at line XXX.

FOAM aborting

I trace back every vertices and make sure the pair of vertices in the internal faces are exactly oppose to each other, only then, the mergePatchPairs works fine

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Issue on zeroGradient BC for system rotation	magnushaese	OpenFOAM Running, Solving & CFD	2	January 28, 2022 14:43
rhoCentralFoam inviscid airfoil issue (Foam::error::printStack(Foam::Ostream&) at ??)	kmkb21	OpenFOAM Running, Solving & CFD	1	March 1, 2018 02:07
Prism Layer - Core Mesh Boundary Issue	InfernoxCJC	STAR-CCM+	7	November 21, 2017 16:46
foamToTecplot360 issue with multiregion solvers	manuc	OpenFOAM Post-Processing	2	November 21, 2016 14:51
Pressure boundary condition issue	Vijay	FLUENT	0	April 6, 2012 14:35