[JBeilke]

Originally Posted by wyldckat

Quick question @JBeilke: Can you please provide more details, namely which tutorial case I can use to test this, so that I can reproduce the problem?

The reason that it worked in 3.0 was, that the library was not compiled there. Otherwise I have to delete the calls to any additional libraries in controlDict, like it is described in the related mantis error entry..

I used the ERCOFTAC tutorial from twcoPhaseMixingCentralFoam. It uses the libCompressiblrTools Library from the same server.

Jörn

[wyldckat]

Hi Jörn,

Originally Posted by JBeilke

The reason that it worked in 3.0 was, that the library was not compiled there. Otherwise I have to delete the calls to any additional libraries in controlDict, like it is described in the related mantis error entry..

Sorry, but now I'm confused... do you mean that the problem is now solved, after you've rebuilt the custom library?

Best regards,
Bruno

[JBeilke]

Originally Posted by wyldckat

Hi Jörn,


Sorry, but now I'm confused... do you mean that the problem is now solved, after you've rebuilt the custom library?

Best regards,
Bruno

No. The reconstruction only works, when the library is missing or the call is removed from controlDict. It is the same problem like it was reported with swak.

[rotschipf]

is this problem solved yet. I can not use cyclicAMI without getting the error message:
Unable to find initial target face
It occurs even with a simple meshed cube with two sides of cyclicAMI.
Please help!

[SirIsaac90]

Hi!

I had the same problem with reconstructPar in OpenFOAM 2.4.0.

The problem was with a *.lib in controlDict, which I had to use to define a porousBafflePressure BC. After deleting it and redefining the BC it worked.

Maybe this helps other Foamers.

Regards
JW

[majid pourdian]

Hi every body
I simulated a circular channel in blockMesh and want to solve it with sonicFoam but it give me the following error:
unable to find initial target face
I used two pair of cyclicAMI Bc's and I attached the picture of my geometry.

[jbcourt]

Hi everyone,
I seem to be experiencing the same issue on OF3.0.1 (full build 3.0.1-d8a290b55d28) - I'm developing a numerical wave tank using the olaFoam toolbox and a transverse AMI cylinder to deal with wave energy device motion. I hadn't been experiencing this issue until I changed my createPatchDict to explicitly state the transform for setup of the AMI patches (see the updated createPatchDict below):

// Do a synchronisation of coupled points after creation of any patches.
// Note: this does not work with points that are on multiple coupled patches
//       with transformations (i.e. cyclics).
pointSync false;

// Patches to create.
patches
(
    {
        //- Master side patch
        name            AMI1;
        patchInfo
        {
            type            cyclicAMI;
            matchTolerance  0.0001;
            neighbourPatch  AMI2;
            transform rotational;
            rotationAxis (0 1 0);
            rotationCentre (0 0 $hingeDepth);
        }
        constructFrom patches;
        patches (AMICylinder);
    }

    {
        //- Slave side patch
        name            AMI2;
        patchInfo
        {
            type            cyclicAMI;
            matchTolerance  0.0001;
            neighbourPatch  AMI1;
            transform rotational;
            rotationAxis (0 1 0);
            rotationCentre (0 0 $hingeDepth);
        }
        constructFrom patches;
        patches (AMICylinder_slave);
    }

As an aside, I'm only about 50% confident that I've defined my rotationAxis and rotationCentre correctly for rotation about a hinge parallel to the y axis, on a line from (0 0 $hingeDepth) to (0 2.5 $hingeDepth).

Anyway, previously I'd set transform to noOrdering as it is in the propeller tutorial and had been experiencing issues whenever I tried to run the simulation to solve for device motion due to wave-structure interaction (another issue which I've started a separate thread for http://www.cfd-online.com/Forums/ope...mi-motion.html), so was trying different things in an attempt to get the thing to run.

Back to the matter in hand, I'm going to try editing my Allrun script to comment out the olaFoam libraries in controlDict just while setFields - which is where the problem starts - is running and then bring it back for decomposePar and the solver. I'll update with how I get on (I'm a bit of a novice in bash, so it might take me a while!).

The output from log.setFields is below, just for reference.

Create time

Create mesh for time = 0

Reading setFieldsDict

Setting field default values
    Setting internal values of volScalarField alpha.water
AMI: Creating addressing and weights between 33904 source faces and 34112 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     : (5.20002 2.499 2.92903)
    target box span     : (4.58053 2.499 4.72778)
    source box          : (-2.60001 0 -1.8) (2.60001 2.499 1.12903)
    target box          : (-2.59996 0 -3.59879) (1.98057 2.499 1.12899)
    inflated target box : (-2.95202 -0.352059 -3.95085) (2.33263 2.85106 1.48105)


--> FOAM FATAL ERROR: 
Unable to set source and target faces

    From function void Foam::faceAreaWeightAMI<SourcePatch, TargetPatch>::setNextFaces(label&, label&, label&, const boolList&, labelList&, const DynamicList<label>&, bool) const
    in file lnInclude/faceAreaWeightAMI.C at line 300.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/lib/libOpenFOAM.so"
#2  Foam::faceAreaWeightAMI<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> > >::calcAddressing(Foam::List<Foam::DynamicList<long, 0u, 2u, 1u> >&, Foam::List<Foam::DynamicList<double, 0u, 2u, 1u> >&, Foam::List<Foam::DynamicList<long, 0u, 2u, 1u> >&, Foam::List<Foam::DynamicList<double, 0u, 2u, 1u> >&, long, long) in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/lib/libmeshTools.so"
#3  Foam::faceAreaWeightAMI<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> > >::calculate(Foam::List<Foam::List<long> >&, Foam::List<Foam::List<double> >&, Foam::List<Foam::List<long> >&, Foam::List<Foam::List<double> >&, long, long) in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/lib/libmeshTools.so"
#4  Foam::AMIInterpolation<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> > >::update(Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > const&, Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > const&) in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/lib/libmeshTools.so"
#5  Foam::AMIInterpolation<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> > >::constructFromSurface(Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > const&, Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > const&, Foam::autoPtr<Foam::searchableSurface> const&) in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/lib/libmeshTools.so"
#6  Foam::cyclicAMIPolyPatch::resetAMI(Foam::AMIInterpolation<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> > >::interpolationMethod const&) const in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/lib/libmeshTools.so"
#7  Foam::cyclicAMIPolyPatch::AMI() const in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/lib/libmeshTools.so"
#8  Foam::cyclicAMIPolyPatch::applyLowWeightCorrection() const in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/lib/libmeshTools.so"
#9  Foam::cyclicAMIFvPatchField<double>::patchNeighbourField() const in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/lib/libfiniteVolume.so"
#10  Foam::coupledFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/lib/libfiniteVolume.so"
#11  Foam::cyclicAMIFvPatchField<double>::cyclicAMIFvPatchField(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/lib/libfiniteVolume.so"
#12  Foam::fvPatchField<double>::adddictionaryConstructorToTable<Foam::cyclicAMIFvPatchField<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/lib/libfiniteVolume.so"
#13  Foam::fvPatchField<double>::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/bin/setFields"
#14  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::readField(Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/bin/setFields"
#15  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/bin/setFields"
#16  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields() in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/bin/setFields"
#17  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&) in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/bin/setFields"
#18  ? in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/bin/setFields"
#19  ? in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/bin/setFields"
#20  ? in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/bin/setFields"
#21  ? in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/bin/setFields"
#22  __libc_start_main in "/lib64/libc.so.6"
#23  ? in "/apps/openfoam/gcc/3.0.1/OpenFOAM-3.0.1/platforms/linux64GccDPInt64Opt/bin/setFields"

[zordiack]

--snip--

Create time

Create mesh for time = 0

Reading setFieldsDict

Setting field default values
    Setting internal values of volScalarField alpha.water
AMI: Creating addressing and weights between 33904 source faces and 34112 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     : (5.20002 2.499 2.92903)
    target box span     : (4.58053 2.499 4.72778)
    source box          : (-2.60001 0 -1.8) (2.60001 2.499 1.12903)
    target box          : (-2.59996 0 -3.59879) (1.98057 2.499 1.12899)
    inflated target box : (-2.95202 -0.352059 -3.95085) (2.33263 2.85106 1.48105)


--> FOAM FATAL ERROR: 
Unable to set source and target faces

--snip--

You have a geometrical problem with your AMI patches. As you can see, AMI is telling you that the target and source patches do not match (size is different in X and Z co-ordinates). Check your geometrical patch assingments, as for AMI the patches have to match geometrically exactly (or at least very very closely).

Maybe you should try specifying also rotationAngle? Check this: http://openfoam.org/release/2-3-0/non-conforming-ami/. Documentation or examples are pretty hard to find for how to use rotationAngle, I strugled with a similar issue myself some time ago.

[jbcourt]

Thanks Pekka, that link does mention the rotationCentre and rotationAxis keywords that I'm unsure about - I generate the patches that are neighbours for my AMI cylinder (i.e. AMI1 and AMI2) from a user-defined cylinder in snappyHexMesh - AMICylinder - from which another patch, AMICylinder_slave is derived. This pair are then used as the two AMI patches, so they should conform perfectly to one-another.

I've successfully (I think, as the solver log file seemed sensible - but I'm still transferring over from the supercomputer to check the results) run this with prescribed motion for the AMI cylinder: the only difference being that I'd set the transform to noOrdering. I'd been running my case with calculated motion (sixDoFRigidBodyMotion solver) and with the transform set to noOrdering too, and only getting errors once the solver had run to about 0.02s - so the meshing was fine with that.

This setFields error has only arisen since I changed the transform setting in createPatchDict - but there must have been something wrong with the way I'd implemented my AMI before for the solver not to be able to run the simulation to completion.

I think therefore that the issue here could lie in how I've set the rotationCentre and rotationAxis parameters. Either that or it's the bug mentioned earlier in the thread that leads to this same error.
Thanks again,
James

[jbcourt]

To be honest I'm a little unclear as to what rotationAngle does. I've had a quick grep through the tutorials that come with OF3.0.1 and I can't see it used anywhere - I definitely agree it's hard to find documentation for it!

[majid pourdian]

Hi Pekka
my patch was rotated so it's bounding box is different in all direction. the dx and dy of my and area of my patches is equal but I have error ani I cant resolve it.

[jbcourt]

Originally Posted by jbcourt

Back to the matter in hand, I'm going to try editing my Allrun script to comment out the olaFoam libraries in controlDict just while setFields - which is where the problem starts - is running and then bring it back for decomposePar and the solver. I'll update with how I get on (I'm a bit of a novice in bash, so it might take me a while!).

I've just implemented this by creating two versions of controlDict - one with the external libraries and one without - and switching to the version without external libraries for setFields: the error still cropped up so I think mine is an issue with how I've entered rotationCentre and rotationAxis. Guess I'll have to start trying different possibilities for those entries to try to get it to work!

[Sugajen]

Originally Posted by dhruv

Thanks for the reply.

1. The case is stationary. So, I dont have a dynamicMeshDict.
2. I tried removing the faces, and making the boundary file like you have posted above, but still I get the same error.
3. I will now try to refine the mesh better to see if it works

Thanks,

Dhruv.

Hi all,

I am facing a similar error "unable to find initial target face" using cyclicAMI on snappyHexMesh. I have refined my mesh multiple times but that does not help.

@Dhruv were you able sort out your issue ?

[hooman.4028]

Originally Posted by wyldckat

Hi Jörn,


Sorry, but now I'm confused... do you mean that the problem is now solved, after you've rebuilt the custom library?

Best regards,
Bruno

Hello Bruno,

I read pretty much everything in this forum on this issue but have not found any solid solution. I have a simple geometry and would like to assign cyclicAMI after meshing with SHM. First I got the error "Unable to set source and target faces
" and managed to solve it by refining my surface mesh on these two cyclic BCs. Now I am getting "unable to find initial target face" and have no clue how to resolve it. Do you guys have any hint on that? There are just two sides of a box and have the exact same geometry .

[akin1078]

I have the same problem. The all day I was trying to solve where this problem comes from. And unfortunately I could not find. So my case is;
1-Flow in a channel
2-Using PimpleDyMFoam
3-There is a rotating disc inside of channel and I am using cyclicAMI
4-The inlet and outlet of channel are cyclic
The decomposePar is running and the simulation also. But when I reconstruct the processors, I am geting below error;


[QUOTE]
Reconstructing FV fields

Reconstructing volScalarFields

p
AMI: Creating addressing and weights between 2224 source faces and 1364 target faces


--> FOAM FATAL ERROR:
Unable to set source and target faces
[QUOTE]


Below solutions was already tried and could not fix:
1-change decomposeParDict from hierarchical to simple. So it is not related with the method of decomposition.
2-use preservePatch in order to put the cyclic boundaries in to one processors. But it gave me the same problem.
3-I am using HPC and I need at least 160 processors. So I deacreased the number of processors from (1 1 160) to
(1 1 4) -> okay
(1 1 40) -> okay
(1 2 40) -> reconstractPar is not running
(1 2 40) -> reconstractPar is not running
(2 1 40) -> okay
(1 1 80) -> reconstractPar is not running
(2 1 60) -> reconstractPar is not running
So it changes according to number of processors obviously but how???
4-I changed the cyclic boundary condition to normal boundary condition with keeping number of processors same as (1 1 160) and reconstructPar works without any error. So it means also that it depends on the cyclic boundary condition.
5- change the OpenFOAM version. I used 2.3 , 4 and 5 and all of them are the same. So it is not version dependant.
6-I used another computer and same problem.
So how it is possible that reconstractPar depends on both cyclic B.C. and number of processors. Even I used the preservePatch, it means that it should not be related with cyclic B.C.
Any ideas???

[Bashar]

[QUOTE=akin1078;662166]I have the same problem. The all day I was trying to solve where this problem comes from. And unfortunately I could not find. So my case is;
1-Flow in a channel
2-Using PimpleDyMFoam
3-There is a rotating disc inside of channel and I am using cyclicAMI
4-The inlet and outlet of channel are cyclic
The decomposePar is running and the simulation also. But when I reconstruct the processors, I am geting below error;


[QUOTE]
Reconstructing FV fields

Reconstructing volScalarFields

p
AMI: Creating addressing and weights between 2224 source faces and 1364 target faces


--> FOAM FATAL ERROR:
Unable to set source and target faces

Below solutions was already tried and could not fix:
1-change decomposeParDict from hierarchical to simple. So it is not related with the method of decomposition.
2-use preservePatch in order to put the cyclic boundaries in to one processors. But it gave me the same problem.
3-I am using HPC and I need at least 160 processors. So I deacreased the number of processors from (1 1 160) to
(1 1 4) -> okay
(1 1 40) -> okay
(1 2 40) -> reconstractPar is not running
(1 2 40) -> reconstractPar is not running
(2 1 40) -> okay
(1 1 80) -> reconstractPar is not running
(2 1 60) -> reconstractPar is not running
So it changes according to number of processors obviously but how???
4-I changed the cyclic boundary condition to normal boundary condition with keeping number of processors same as (1 1 160) and reconstructPar works without any error. So it means also that it depends on the cyclic boundary condition.
5- change the OpenFOAM version. I used 2.3 , 4 and 5 and all of them are the same. So it is not version dependant.
6-I used another computer and same problem.
So how it is possible that reconstractPar depends on both cyclic B.C. and number of processors. Even I used the preservePatch, it means that it should not be related with cyclic B.C.
Any ideas???

Hi,

Any chance that you solve this issues? I am facing exactly the same issue when working with AMI!

Bashar

[anhkenyt]

Hi all,
l created a turbine mesh in turborGrid and imported it to OF 6.0, checkMesh is Ok, but when l run PimpleFoam I am facing a similar error :

AMI: Creating addressing and weights between 6625 source faces and 6625 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     : (0.041425 0.0215107 0.096)
    target box span     : (0.0419304 0.0195127 0.096)
    source box          : (0.0874741 0.04003 -0.016) (0.128899 0.0615407 0.08)
    target box          : (0.0897394 0.035125 -0.016) (0.13167 0.0546376 0.08)
    inflated target box : (0.0844115 0.029797 -0.0213279) (0.136998 0.0599656 0.085328)


--> FOAM FATAL ERROR: 
Unable to set source and target faces

My case is pretty similar to the tutorial Impeller in OF 6.0 and I also follow the setup in that tutorial for boundary conditions.

Please have a look at my geometry in the attached images. My boundaries is like below:

12
(
    outlet
    {
        type            patch;
        nFaces          6360;
        startFace       4074577;
    }
    movingWalls
    {
        type            wall;
        inGroups        1(wall);
        nFaces          40680;
        startFace       4080937;
    }
    rotorBlade
    {
        type            wall;
        inGroups        1(wall);
        nFaces          33231;
        startFace       4121617;
    }
    statorBlade
    {
        type            wall;
        inGroups        1(wall);
        nFaces          6976;
        startFace       4154848;
    }
    inlet
    {
        type            patch;
        nFaces          1280;
        startFace       4161824;
    }
    stationnaryWalls
    {
        type            wall;
        inGroups        1(wall);
        nFaces          18920;
        startFace       4163104;
    }
    cyclicRepeatAMIInterface1
    {
        type            cyclicRepeatAMI;
        inGroups        1(RepeatAMI1);
        nFaces          1792;
        startFace       4182024;
        inGroups        1 ( RepeatAMI1 );
        name            cyclicRepeatAMIInterface1;
        neighbourPatch  cyclicRepeatAMIInterface2;
        transformPatch  cyclic_Out1;
    }
    cyclicRepeatAMIInterface2
    {
        type            cyclicRepeatAMI;
        inGroups        1(RepeatAMI2);
        nFaces          5565;
        startFace       4183816;
        inGroups        1 ( RepeatAMI2 );
        name            cyclicRepeatAMIInterface2;
        neighbourPatch  cyclicRepeatAMIInterface1;
        transformPatch  cyclic_Out1;
    }
    cyclic_In1
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          3616;
        startFace       4189381;
        matchTolerance  0.0001;
        transform       rotational;
        neighbourPatch  cyclic_In2;
        rotationAxis    (0 0 1);
        rotationCentre  (0 0 0);
    }
    cyclic_In2
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          3616;
        startFace       4192997;
        matchTolerance  0.0001;
        transform       rotational;
        neighbourPatch  cyclic_In1;
        rotationAxis    (0 0 1);
        rotationCentre  (0 0 0);
    }
    cyclic_Out1
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          6625;
        startFace       4196613;
        matchTolerance  0.0001;
        transform       rotational;
        neighbourPatch  cyclic_Out2;
        rotationAxis    (0 0 1);
        rotationCentre  (0 0 0);
    }
    cyclic_Out2
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          6625;
        startFace       4203238;
        matchTolerance  0.0001;
        transform       rotational;
        neighbourPatch  cyclic_Out1;
        rotationAxis    (0 0 1);
        rotationCentre  (0 0 0);
    }
)

The answers above say that this is a bug but it has been resolved in new version of openfoam. Is there any chance that the bug re-appear in OF 6?

Could any body help me on this problems?
1.png

2.png

3.png

[homayoun9604]

Originally Posted by Kaskade

Hello.

I've set up a case using 3 AMIs and MRFSimpleFOAM. It runs fine on one processor, it even runs fine on multiple processors. But I am running into trouble when I try to piece it back together after a parrallel run.

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

Create mesh for time = 100

Time = 100

Reconstructing FV fields

    Reconstructing volScalarFields

        p
AMI: Creating addressing and weights between 1628 source faces and 1071 target faces
--> FOAM Warning : 
    From function AMIInterpolation<SourcePatch, TargetPatch>::checkPatches(const primitivePatch&, const primitivePatch&)
    in file lnInclude/AMIInterpolation.C at line 146
    Source and target patch bounding boxes are not similar
    source box span     : (8.5e-10 0.55 0.55)
    target box span     : (1.07e-09 0.0576107 0.238538)
    source box          : (0.0423665 -0.275 -0.275) (0.0423665 0.275 0.275)
    target box          : (0.042366 0.217361 -0.136582) (0.042366 0.274972 0.101956)
    inflated target box : (0.0300962 0.205091 -0.148852) (0.0546358 0.287242 0.114226)


--> FOAM FATAL ERROR: 
Unable to set source and target faces

    From function void Foam::cyclicAMIPolyPatch::setNextFaces(label&, label&, const primitivePatch&, const primitivePatch&, const boolList&, labelList&, const DynamicList<label>&) const
    in file lnInclude/AMIInterpolation.C at line 878.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   at cyclicAMIPolyPatch.C:0
#3  Foam::AMIInterpolation<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> > >::update(Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > const&, Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > const&) in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libmeshTools.so"
#4  Foam::AMIInterpolation<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> > >::AMIInterpolation(Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > const&, Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > const&, Foam::autoPtr<Foam::searchableSurface> const&, Foam::faceAreaIntersect::triangulationMode const&, bool) in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libmeshTools.so"
#5  Foam::cyclicAMIPolyPatch::resetAMI() const in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libmeshTools.so"
#6  Foam::cyclicAMIPolyPatch::AMI() const in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libmeshTools.so"
#7  Foam::tmp<Foam::Field<double> > Foam::cyclicAMIPolyPatch::interpolate<double>(Foam::tmp<Foam::Field<double> > const&) const in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8  Foam::cyclicAMIFvPatch::makeWeights(Foam::Field<double>&) const in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#9  Foam::surfaceInterpolation::makeWeights() const in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#10  Foam::surfaceInterpolation::weights() const in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#11  Foam::fvPatch::weights() const in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#12  Foam::coupledFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#13  Foam::cyclicFvPatchField<double>::cyclicFvPatchField(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#14  Foam::fvPatchField<double>::adddictionaryConstructorToTable<Foam::cyclicFvPatchField<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#15  Foam::fvPatchField<double>::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/reconstructPar"
#16  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::GeometricBoundaryField(Foam::fvBoundaryMesh const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/reconstructPar"
#17  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readField(Foam::dictionary const&) in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/reconstructPar"
#18  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readField(Foam::Istream&) in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/reconstructPar"
#19  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&) in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/reconstructPar"
#20  
 in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/reconstructPar"
#21  
 in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/reconstructPar"
#22  
 in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/reconstructPar"
#23  __libc_start_main in "/lib64/libc.so.6"
#24  
 in "/software/openfoam-2.1.0/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/reconstructPar"

I used scotch to decompose it, which worked on a previous case, set up analogous to the propeller-tutorial. This mesh however has been generated using ICEM and CFX. My first guess was that OF210 can't handle multiple AMIs, but fusing them using createPatch led to the pretty much the same error. Can anyone point me in the right direction?

Hello,
How can I fix this problem?

[jvinuales]

edit: reply on wrong post