[Vinzmann]

Hello everyone,

I am currently trying to save computing power by only simulating one third of my axisymmetric (6-way) case. For this I first tried to add cyclicAMI Boundary conditions in snappyHexMesh.

SnappyHexMesh's output looks something like this when it gets to match the cyclicAMI faces:

AMI: Creating addressing and weights between 19082 source faces and 27777 target faces
--> FOAM Warning :
    From virtual void Foam::faceAreaWeightAMI::calcInterArea(Foam::label, Foam::label, Foam::scalar&, Foam::vector&) const
    in file AMIInterpolation/AMIInterpolation/faceAreaWeightAMI/faceAreaWeightAMI.C at line 424
    Invalid normal for source face 2374 points 4((0.0027848 -0.005 -0.0025) (0.001760275 -0.005 -0.0025) (0.001760275 -0.00625 -0.0025) (0.0027848 -0.00625 -0.0025)) target face 15319 points 4((-0.000573182015 -0.00730039552 -0.0025) (0.000631515189 -0.00696692192 -0.0025) (0.000358193551 -0.00597952799 -0.0025) (-0.000846503653 -0.0063130016 -0.0025))

This goes on for many times and ends with:

--> FOAM FATAL ERROR: (openfoam-2212)
Unable to set target face for source face 19081

    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 361.

This happens for snappyHexMesh in parallel and in serial mode.

If I set my cyclicAMI boundaries as walls in snappyHexMesh, it meshes fine. After reconstructing, setting the patches from wall to cyclicAMI with createPatch (this runs fine and gives no errors), decomposePar throws the exact same error, namely:

AMI: Creating addressing and weights between 19082 source faces and 27777 target faces
--> FOAM FATAL ERROR: (openfoam-2212)
Unable to set target face for source face 19081

    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 361.

When skipping createPatch, decomposing and then changing the processor*/0/* files to cyclicAMI for the two relevant patches, the solver throws the exact same error as decomposePar above.

Can somebody please explain why this is happening and what I am doing wrong? This would be so helpful.

Best regards,
Vinz

[gp24]

Hello, I have encountered the same issue. Did you find a solution?

[Sakun]

Originally Posted by gp24

Hello, I have encountered the same issue. Did you find a solution?

What is your meshing software ? did you define length between 2 faces in SI units correctly ?

[gp24]

Originally Posted by Sakun

What is your meshing software ? did you define length between 2 faces in SI units correctly ?

I am using snappyHexMesh, I defined correctly the length between the faces. I usually start by defining in blockMesh the periodic patches as cyclic and then with createPatch I redefine them as cyclicAMI, I have tried defining it directly as cyclicAMI but the issue persists.