[nainavinod]

Hello Foamers!

I am trying ot mesh an infinite swept wing geometry such that the sides have a periodic boundary condition. I tried using the "cyclic" condition initially, however, due to the unstructured meshes being different on either side, it didn't work out. I found out that "cyclicAMI" would be a better fit for my case. My current workflow for the final mesh is:
1) blockMesh
2) surfaceFeatureExtract
3) snappyHexMesh
4) topoSet
5) createPatch

When I do a general mesh quality check with the "-allGeometry" flag, I get the following error in the end:

Calculating AMI weights between owner patch: Side1 and neighbour patch: Side2
AMI: Creating addressing and weights between 111996 source faces and 111341 target faces


--> FOAM FATAL ERROR:
Unable to set target face for source face 111995

    From virtual bool Foam::faceAreaWeightAMI::setNextFaces(Foam::label&, Foam::label&, Foam::label&, const Foam::bitSet&, Foam::labelList&, const Foam::DynamicList<int>&, bool) const
    in file AMIInterpolation/AMIInterpolation/faceAreaWeightAMI/faceAreaWeightAMI.C at line 347.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::exitOrAbort(int, bool) at ??:?
#2  Foam::faceAreaWeightAMI::setNextFaces(int&, int&, int&, Foam::bitSet const&, Foam::List<int>&, Foam::DynamicList<int, 16> const&, bool) const at ??:?
#3  Foam::faceAreaWeightAMI::calcAddressing(Foam::List<Foam::DynamicList<int, 16> >&, Foam::List<Foam::DynamicList<double, 16> >&, Foam::List<Foam::DynamicList<Foam::Vector<double>, 16> >&, Foam::List<Foam::DynamicList<int, 16> >&, Foam::List<Foam::DynamicList<double, 16> >&, int, int) at ??:?
#4  Foam::faceAreaWeightAMI::calculate(Foam::PrimitivePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::PrimitivePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::autoPtr<Foam::searchableSurface> const&) at ??:?
#5  Foam::cyclicAMIPolyPatch::resetAMI(Foam::UList<Foam::Vector<double> > const&) const at ??:?
#6  Foam::cyclicAMIPolyPatch::AMI() const at ??:?
#7  ? at ??:?
#8  ? at ??:?
#9  __libc_start_main in /lib64/libc.so.6
#10  ? at ??:?

I tried looking this error up in the forum, but most of the suggestions are for moving/dynamic meshes. While, mine is a stationary case. I couldn't figure out where exactly I am going wrong. My boundary file is written as:

7
(
    Inlet
    {
        type            patch;
        nFaces          2016;
        startFace       35454828;
    }
    Outlet
    {
        type            patch;
        nFaces          2016;
        startFace       35456844;
    }
    Bottom
    {
        type            slip;
        inGroups        1(slip);
        nFaces          3780;
        startFace       35458860;
        faces           ( ( 0 1 5 4 ) );
    }
    Top
    {
        type            slip;
        inGroups        1(slip);
        nFaces          3780;
        startFace       35462640;
        faces           ( ( 3 7 6 2 ) );
    }
    Wing
    {
        type            wall;
        inGroups        1(wall);
        nFaces          145149;
        startFace       35466420;
    }
    Side1
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          111996;
        startFace       35611569;
        matchTolerance  0.0001;
        transform       unknown;
        neighbourPatch  Side2;
        AMIMethod       faceAreaWeightAMI;
        restartUncoveredSourceFace 1;
    }
    Side2
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          111341;
        startFace       35723565;
        matchTolerance  0.0001;
        transform       unknown;
        neighbourPatch  Side1;
        AMIMethod       faceAreaWeightAMI;
        restartUncoveredSourceFace 1;
    }
)

blockMeshDict file:

convertToMeters 1;

vertices
(
    ( -9.00   -7.20  -0.70 )  //0
    (  18.00  -7.20  -0.70 )  //1
    (  18.00   7.20  -0.70 )  //2
    ( -9.00    7.20  -0.70 )  //3

    ( -9.00   -7.20  0.70 )  //4
    (  18.00  -7.20  0.70 )  //5
    (  18.00   7.20  0.70 )  //6
    ( -9.00    7.20  0.70 )  //7
);

blocks
(
    hex (0 1 2 3 4 5 6 7) (270 144 14) simpleGrading (1 1 1)
);

edges
(

);

boundary
(
  
    Inlet
    {
	type    patch;
	faces   ((0 4 7 3));
    }
    Outlet
    {
	type	patch;
	faces	((1 2 6 5));
    }
    Bottom
    {
        type    slip;
        faces   ((0 1 5 4));
    }
    Top
    {
	type	slip;
        faces   ((3 7 6 2));
    }
    Sides
    {
        type    patch;
        //neighbourPatch  Side2;
        faces   ((0 3 2 1)
		 (4 5 6 7));
    }
   
)

topoSetDict file:

actions
(
    {
        name    Side2;
        type    faceSet;
        action  new;
        source  boxToFace;
        sourceInfo
        {
            box (-9 -7.2 0.69991)(18 7.2 0.70001);
        }
    }

    {
        name    Side1;
        type    faceSet;
        action  new;
        source  boxToFace;
        sourceInfo
        {
            box (-9 -7.2 -0.70001)(18 7.2 -0.69991);
        }
    }

);

And finally, the createPatchDict file:

pointSync false;

patches
(
    {
        // Name of new patch
        name Side1;

        // Dictionary to construct new patch from
        patchInfo
        {
            type cyclicAMI;
            neighbourPatch Side2;
     	    //matchTolerance 1e-04;
        }

        // How to construct: either from 'patches' or 'set'
        constructFrom set;

        // If constructFrom = set : name of faceSet
        set  Side1 ;
    }

    {
        // Name of new patch
        name Side2;

        // Dictionary to construct new patch from
        patchInfo
        {
            type cyclicAMI;
            neighbourPatch Side1;
   	    //matchTolerance 1e-04;
        }

        // How to construct: either from 'patches' or 'set'
        constructFrom set;

        // If constructFrom = set : name of faceSet
        set  Side2;
    }
);

Let me know if more details on my case is required. Any help or suggestions would be greatly appreciated.

[nainavinod]

I am still stuck on this issue and it would be great if someone can help me out.

Thanks,
Naina.