[duguyoudi]

Dear all,


I generated a mesh using blockMesh and created cyclicAMI boundaries. The process of creating cyclicAMI seems fine. The results is

Create polyMesh for time = 0

Reading createPatchDict

Adding new patch left as patch 6 from 
{
    type            cyclicAMI;
    neighbourPatch  right;
    transform       translational;
    separationVector ( 0 1 0 );
    matchTolerance  1e-06;
}

Adding new patch right as patch 7 from 
{
    type            cyclicAMI;
    neighbourPatch  left;
    transform       translational;
    separationVector ( 0 -1 0 );
    matchTolerance  1e-06;
}


Moving faces from patch left1 to patch 6
Moving faces from patch right1 to patch 7

Doing topology modification to order faces.

Not synchronising points.

Removing patches with no faces in them.

Removing zero-sized patch left1 type patch at position 1
Removing zero-sized patch right1 type patch at position 2
Removing patches.
Writing repatched mesh to 0

End

But when I tried to check the mesh using paraFoam, I got an error message:

AMI: Creating addressing and weights between 3600 source faces and 3600 target faces
--> FOAM Warning : 
    From function void Foam::AMIMethod<SourcePatch, TargetPatch>::checkPatches() const [with SourcePatch = Foam::PrimitivePatch<Foam::face, Foam::SubList, const Foam::Field<Foam::Vector<double> >&>; TargetPatch = Foam::PrimitivePatch<Foam::face, Foam::SubList, const Foam::Field<Foam::Vector<double> >&>]
    in file lnInclude/AMIMethod.C at line 57
    Source and target patch bounding boxes are not similar
    source box span     : (1.6 3.33067e-16 0.9)
    target box span     : (1.6 3.33067e-16 0.9)
    source box          : (-0.4 -1 0) (1.2 -1 0.9)
    target box          : (-0.4 -1.11022e-16 0) (1.2 2.22045e-16 0.9)
    inflated target box : (-0.491788 -0.0917878 -0.0917878) (1.29179 0.0917878 0.991788)


--> FOAM FATAL ERROR: 
Unable to find initial target face

Here is blockMeshDict

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.1;

vertices
(
    (-4 -10 0)  //vertex 0
    (12 -10 0)  //vertex 1
    (12  10 0)  //vertex 2
    (-4  10 0)  //vertex 3

    (-4 -10 9)  //vertex 5
    (12 -10 9)  //vertex 6
    (12  10 9)  //vertex 7
    (-4  10 9)  //vertex 8
);

blocks
(
    hex (0 1 2 3 4 5 6 7) (80 100 45) 
    simpleGrading (1 1 1)
);
edges
(
);

boundary
(
    walls
    {
        type wall;
        faces
        (
            (0 1 2 3)
        );
    }
    left1
    {
        //type cyclic;
        //neighbourPatch right;
        type patch;
        faces
        (
            (0 1 5 4)
        );
    }
    right1
    {
        //type cyclic;
        //neighbourPatch left;
        type patch;
        faces
        (
            (3 2 6 7)
        );
    }
    inlet
    {
        type patch;
        faces
        (
            (0 3 7 4)
        );
    }
    outlet
    {
        type patch;
        faces
        (
            (1 2 6 5)
        );
    }
    top
    {
        type patch;
        faces
        (
            (4 5 6 7)
        );
    }
);

and createPatchDict

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

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

pointSync false;

patches
(
   {
      name left;
      patchInfo
      {
         type cyclicAMI;
         neighbourPatch right;

         transform translational;
         separationVector (0 1 0);
         matchTolerance 1e-6;
      }
      constructFrom patches;
      patches (left1);
   }
   {
      name right;
      patchInfo
      {
         type cyclicAMI;
         neighbourPatch left;

         transform translational;
         separationVector (0 -1 0);
         matchTolerance 1e-6;
      }
      constructFrom patches;
      patches (right1);
   }
);

The strange thing is that when I changed the parameter in blockMeshDict like this:

vertices
(
    (-4 -5 0)  //vertex 0
    (12 -5 0)  //vertex 1
    (12  5 0)  //vertex 2
    (-4  5 0)  //vertex 3

    (-4 -5 9)  //vertex 5
    (12 -5 9)  //vertex 6
    (12  5 9)  //vertex 7
    (-4  5 9)  //vertex 8
);

blocks
(
    hex (0 1 2 3 4 5 6 7) (80 50 45) 
    simpleGrading (1 1 1)
);

This time every thing is fine.


Does anyone have any ideas about this problem?


Best regards,


Leo

[duguyoudi]

For anyone who is interested.


The " separationVector" should be (0 2 0) and (0 -2 0) for left and right respectively. The vector is not normalized, i.e. ||v||=1. The length of the vector should be the distance between two planes.