[worldsheep]

Hello everyone,

I am simulating a rotor-stator system using different methods starting from the examples provided in openfoam (ex. mixerVesselAMI2D).

Every time I use an AMI interface I get the following warings:

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 414    
Invalid normal for source face 0 points 3((0.00193882 0.0019985 -0.00016) (0.00188455 0.00193079 -0.00016) (0.00185038 0.00199937 -0.00016)) target face
41675 points 3((-0.00230661 0.00266591 0.00065) (-0.00225043 0.00274269 0.00065) (-0.00222007 0.00264246 0.00065))

This is independent of the solver. I get the same error for:

• Stationary MRF simulations using AMI solved with simpleFoam

• Transient sumulations using AMI solved with pimpleFoam

My mesh is generated using Salome. I mesh the rotor region and the stator region separately with netgen. To make it conformal at time zero, I import the 1D-2D elements from the rotor to mesh the stator and then make a compound with the options below:

• Processing identical groups -> Unite

• Create groups from input objects -> unchecked

• Merge coincident nodes and elements -> unckeked

The mesh is then exported using: https://github.com/nicolasedh/salomeToOpenFOAM

In the boundary file:

    AMI_in
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          41676;
        startFace       4200354;
        matchTolerance  0.0001;
        transform       noOrdering;
        neighbourPatch  AMI_out;
        AMIMethod       faceAreaWeightAMI;
        restartUncoveredSourceFace 1;
    }
    AMI_out
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          41676;
        startFace       4459504;
        matchTolerance  0.0001;
        transform       noOrdering;
        neighbourPatch  AMI_in;
        AMIMethod       faceAreaWeightAMI;
        restartUncoveredSourceFace 1;
    }

In my velocity boundary conditions:

    AMI_in
    {
        type            cyclicAMI;
        value           uniform (0 0 0);
    }
    AMI_out
    {
        type            cyclicAMI;
        value           uniform (0 0 0);
    }

In my pressure boundary condition:

    AMI_in
    {
        type            cyclicAMI;
        value           uniform 0;
    }
    AMI_out
    {
        type            cyclicAMI;
        value           uniform 0;
    }

Has anyone an idea why is this warining generated?