[Alasir]

Hi

Im having a channel flow with a repeating geometry inside. To simplify this I'm using a LES model where I try to use the cyclic BC.

Im using FreeCAD and SnappyHexMesh to generate the mesh, thus causing the iunlet and outlet to have a slight difference in faces. I try to compensate by using cyclicAMI.

However, I keep getting the error "Unable to find initial target face". Does anyone know what this means?

Boundary:
    inlet
    {
        type            cyclicAMI;
        nFaces          3128;
        startFace       400334;
        matchTolerance  0.01;
        transform       noOrdering;
        neighbourPatch  outlet;
    }
    outlet
    {
        type            cyclicAMI;
        nFaces          3172;
        startFace       403462;
        matchTolerance  0.01;
        transform       noOrdering;
        neighbourPatch  inlet;

U:
    inlet
    {
        type            fixedValue;
        value           uniform (0.01 0 0);

    }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform (0 0 0);
        value           uniform (0.01 0 0);
    }

P:
    inlet
    {
        type            cyclicAMI;
    }
    outlet
    {
        type            cyclicAMI;
    }

[Alasir]

Resolved.

Turns out that my error was using noOrdering when I should use translational.

Boundary:
    inlet
    {
        type            cyclicAMI;
        nFaces          3128;
        startFace       400334;
        matchTolerance  0.001;
        transform       translational;
        neighbourPatch  outlet;
        separationVector (0.000480 0 0);
    }
    outlet
    {
        type            cyclicAMI;
        nFaces          3172;
        startFace       403462;
        matchTolerance  0.001;
        transform       translational;
        neighbourPatch  inlet;
        separationVector (-0.000480 0 0);


P:
    inlet
    {
        type            cyclicAMI;

    }
    outlet
    {
        type            cyclicAMI;

    }

[minzhang]

Hello Alasir,

I am having similar issues with the implementation of cyclicAMI boundary condition.
Your kind help would be greatly appreciated!

My questions are as follows:
(1) My patch2 and patch3 are quite different, does it matter?
(2) What does "translational" mean here?
(3) How to set the "separationVector" value? My patch2 position in y-direction is y=0 and patch3 position is y=2.19.

Thanks in advance!!

My CFD/constant/polyMesh/boundary file is attached below:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/cellZone_0/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

5
(
patch2
{
type cyclicAMI;//symmetry;
physicalType cyclicAMI;//symmetry;
neighbourPatch patch3;
nFaces 406;
startFace 91945;
transform translational;
separationVector (0 -2.19 0);
}
patch3
{
type cyclicAMI;//symmetry;
physicalType cyclicAMI;//symmetry;
neighbourPatch patch2;
nFaces 351;
startFace 92351;
transform translational;
separationVector (0 2.19 0);
}
inlet
{
type patch;
physicalType patch;
nFaces 316;
startFace 92702;
}
outlet
{
type patch;
physicalType patch;
nFaces 160;
startFace 93018;
}
walls
{
type wall;
nFaces 42871;
startFace 93178;
}
)

My errors are here:
Number of processor faces = 3643
Max number of cells = 1156 (46.0249% above average 791.646)
Max number of processor patches = 5 (41.1765% above average 3.54167)
Max number of faces between processors = 229 (50.8647% above average 151.792)

Time = 0
AMI: Creating addressing and weights between 406 source faces and 351 target faces
--> FOAM Warning :
From function AMIMethod<SourcePatch, TargetPatch>::checkPatches()
in file lnInclude/AMIMethod.C at line 57
Source and target patch bounding boxes are not similar
source box span : (3.718 0 0.1752)
target box span : (3.9 0 0.219)
source box : (0 0 0.0219) (3.718 0 0.1971)
target box : (0 4.38 0.0219) (3.9 4.38 0.2409)
inflated target box : (-0.195307 4.18469 -0.173407) (4.09531 4.57531 0.436207)


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

From function void Foam::AMIMethod<SourcePatch, TargetPatch>::initialise(labelListList&, scalarListList&, labelListList&, scalar\
ListList&, label&, label&)
in file lnInclude/AMIMethod.C at line 149.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::AMIMethod<Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> \
>, Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > >::initialise(F\
oam::List<Foam::List<int> >&, Foam::List<Foam::List<double> >&, Foam::List<Foam::List<int> >&, Foam::List<Foam::List<double> >&, int\
&, int&) at ??:?

[Vishsel]

Originally Posted by Alasir

Resolved.

Turns out that my error was using noOrdering when I should use translational.

Boundary:
    inlet
    {
        type            cyclicAMI;
        nFaces          3128;
        startFace       400334;
        matchTolerance  0.001;
        transform       translational;
        neighbourPatch  outlet;
        separationVector (0.000480 0 0);
    }
    outlet
    {
        type            cyclicAMI;
        nFaces          3172;
        startFace       403462;
        matchTolerance  0.001;
        transform       translational;
        neighbourPatch  inlet;
        separationVector (-0.000480 0 0);


P:
    inlet
    {
        type            cyclicAMI;

    }
    outlet
    {
        type            cyclicAMI;

    }

How to find the value #(0 0.000480 0) for separation vector??

Anyone ?

Thankyou in advance
Vishsel

[ameyer]

For the sake of completeness and because I had the same question:


As far as I understood (and it seems to work in my case) the


separationVector (x y z)


is just the offset-vector (or translation vector) between your 2 patches - as seen from the patch where it is defined (as in the example above you have to specify it for both patches, so the second
separationVector has the opposite dir. as the first one)

[openfoam_aero]

Originally Posted by minzhang

Hello Alasir,

I am having similar issues with the implementation of cyclicAMI boundary condition.
Your kind help would be greatly appreciated!

My questions are as follows:
(1) My patch2 and patch3 are quite different, does it matter?
(2) What does "translational" mean here?
(3) How to set the "separationVector" value? My patch2 position in y-direction is y=0 and patch3 position is y=2.19.

Thanks in advance!!

My CFD/constant/polyMesh/boundary file is attached below:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/cellZone_0/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

5
(
patch2
{
type cyclicAMI;//symmetry;
physicalType cyclicAMI;//symmetry;
neighbourPatch patch3;
nFaces 406;
startFace 91945;
transform translational;
separationVector (0 -2.19 0);
}
patch3
{
type cyclicAMI;//symmetry;
physicalType cyclicAMI;//symmetry;
neighbourPatch patch2;
nFaces 351;
startFace 92351;
transform translational;
separationVector (0 2.19 0);
}
inlet
{
type patch;
physicalType patch;
nFaces 316;
startFace 92702;
}
outlet
{
type patch;
physicalType patch;
nFaces 160;
startFace 93018;
}
walls
{
type wall;
nFaces 42871;
startFace 93178;
}
)

My errors are here:
Number of processor faces = 3643
Max number of cells = 1156 (46.0249% above average 791.646)
Max number of processor patches = 5 (41.1765% above average 3.54167)
Max number of faces between processors = 229 (50.8647% above average 151.792)

Time = 0
AMI: Creating addressing and weights between 406 source faces and 351 target faces
--> FOAM Warning :
From function AMIMethod<SourcePatch, TargetPatch>::checkPatches()
in file lnInclude/AMIMethod.C at line 57
Source and target patch bounding boxes are not similar
source box span : (3.718 0 0.1752)
target box span : (3.9 0 0.219)
source box : (0 0 0.0219) (3.718 0 0.1971)
target box : (0 4.38 0.0219) (3.9 4.38 0.2409)
inflated target box : (-0.195307 4.18469 -0.173407) (4.09531 4.57531 0.436207)


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

From function void Foam::AMIMethod<SourcePatch, TargetPatch>::initialise(labelListList&, scalarListList&, labelListList&, scalar\
ListList&, label&, label&)
in file lnInclude/AMIMethod.C at line 149.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::AMIMethod<Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> \
>, Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > >::initialise(F\
oam::List<Foam::List<int> >&, Foam::List<Foam::List<double> >&, Foam::List<Foam::List<int> >&, Foam::List<Foam::List<double> >&, int\
&, int&) at ??:?

Hello Min,


Did you find a solution to this? I am facing a similar problem.

[Gang Wang]

Dear Mr. Zhang,

Have you solved this issue? I've met the same error with urs for several days but I cannot find the solution.


Best regards,
Gang

Originally Posted by minzhang

Hello Alasir,

I am having similar issues with the implementation of cyclicAMI boundary condition.
Your kind help would be greatly appreciated!

My questions are as follows:
(1) My patch2 and patch3 are quite different, does it matter?
(2) What does "translational" mean here?
(3) How to set the "separationVector" value? My patch2 position in y-direction is y=0 and patch3 position is y=2.19.

Thanks in advance!!

My CFD/constant/polyMesh/boundary file is attached below:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/cellZone_0/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

5
(
patch2
{
type cyclicAMI;//symmetry;
physicalType cyclicAMI;//symmetry;
neighbourPatch patch3;
nFaces 406;
startFace 91945;
transform translational;
separationVector (0 -2.19 0);
}
patch3
{
type cyclicAMI;//symmetry;
physicalType cyclicAMI;//symmetry;
neighbourPatch patch2;
nFaces 351;
startFace 92351;
transform translational;
separationVector (0 2.19 0);
}
inlet
{
type patch;
physicalType patch;
nFaces 316;
startFace 92702;
}
outlet
{
type patch;
physicalType patch;
nFaces 160;
startFace 93018;
}
walls
{
type wall;
nFaces 42871;
startFace 93178;
}
)

My errors are here:
Number of processor faces = 3643
Max number of cells = 1156 (46.0249% above average 791.646)
Max number of processor patches = 5 (41.1765% above average 3.54167)
Max number of faces between processors = 229 (50.8647% above average 151.792)

Time = 0
AMI: Creating addressing and weights between 406 source faces and 351 target faces
--> FOAM Warning :
From function AMIMethod<SourcePatch, TargetPatch>::checkPatches()
in file lnInclude/AMIMethod.C at line 57
Source and target patch bounding boxes are not similar
source box span : (3.718 0 0.1752)
target box span : (3.9 0 0.219)
source box : (0 0 0.0219) (3.718 0 0.1971)
target box : (0 4.38 0.0219) (3.9 4.38 0.2409)
inflated target box : (-0.195307 4.18469 -0.173407) (4.09531 4.57531 0.436207)


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

From function void Foam::AMIMethod<SourcePatch, TargetPatch>::initialise(labelListList&, scalarListList&, labelListList&, scalar\
ListList&, label&, label&)
in file lnInclude/AMIMethod.C at line 149.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::AMIMethod<Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> \
>, Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > >::initialise(F\
oam::List<Foam::List<int> >&, Foam::List<Foam::List<double> >&, Foam::List<Foam::List<int> >&, Foam::List<Foam::List<double> >&, int\
&, int&) at ??:?