[rbbarber]

Hello,

I am trying to model a tidal turbine using the arbitrary mesh interface and pimpleDyMFoam. Right now I have everything meshing nicely (as far as I can tell), and it will run happily (again, as far as I can tell) for ~40-50 time steps, but then it crashes with the error "unable to set source and target faces." This happens right after "AMI: Creating addressing and weights between 122366 source faces and 123756 target faces" in the log, but those face numbers are consistent with every previous step. I am interpreting this as it's having a hard time solving across the rotating AMI boundary... maybe because they aren't meshed identically? But I don't know how to fix this (or even if it's the issue). The .stl files for the AMI and AMI_slave meshes are identical.

Any ideas would be greatly appreciated. Thanks,
R

Update: it seems like the two AMI meshes start out aligned and then start moving apart. Why??

[Elena89]

Hi R,

are you using solidBodyMotionFvMesh? Rotation?
To me it sounds like the rotational axis (if it's a rotation) isn't normal to the AMI faces and that's why they move apart.

regards

Elena

[rbbarber]

Yes, you're right. I had changed that in a previous simulation but forgot to in this one. Thank you,
R

[shantanumandhane]

hi, I am running a case of fluid flow across cylindrical body. Although my rotation axis is is along the AMI-cylinder's axis, it's still showing the same error. Please help.
I am attaching the case folder here.


https://drive.google.com/drive/folde...Em?usp=sharing

[rokasbytautas]

Have you managed to solve this?

[Pratik Wani]

Hello all

I am working on cfd simulation of vertical axis wind turbine. I ran snappyHexMesh and createPatch. Later, I am getting the following error when I run " moveDynamicMesh -checkAMI"

Create time

Create mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: solidBody
Selecting solid-body motion function rotatingMotion
Applying solid body motion to cellZone rotatingZone
Writing VTK files with weights of AMI patches.


PIMPLE: no residual control data found. Calculations will employ 2 corrector loops

Time = 0.001
PIMPLE: iteration 1
PIMPLE: iteration 2
    Point usage OK.
    Upper triangular ordering OK.
    Topological cell zip-up check OK.
    Face vertices OK.
    Face-face connectivity OK.
    Mesh topology OK.
    Boundary openness (-3.41563e-016 -1.46549e-015 -8.0823e-017) OK.
    Max cell openness = 3.41365e-016 OK.
    Max aspect ratio = 8.01026 OK.
    Minimum face area = 2.65006e-009. Maximum face area = 0.000102174.  Face area magnitudes OK.
    Min volume = 5.14543e-013. Max volume = 1.02146e-006.  Total volume = 1.14749.  Cell volumes OK.
    Mesh non-orthogonality Max: 64.993 average: 10.5499
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 3.09662 OK.
    Mesh geometry OK.
Mesh OK.
Calculating AMI weights between owner patch: AMI1 and neighbour patch: AMI2
AMI: Creating addressing and weights between 426887 source faces and 428124 target faces


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

    From function void Foam::faceAreaWeightAMI<SourcePatch, TargetPatch>::setNextFaces(Foam::label&, Foam::label&, Foam::label&, const boolList&, Foam::labelList&, const Foam::DynamicList<int>&, bool) 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> >&>; Foam::label = int; Foam::boolList = Foam::List<bool>; Foam::labelList = Foam::List<int>]
    in file ./AMIInterpolation/AMIInterpolation/AMIMethod/faceAreaWeightAMI/faceAreaWeightAMI.C at line 287.

FOAM aborting

Generating stack trace...


Backtrace:
        ZN10StackTraceC1Ev [0x705c1465+0x25]
                 module: C:\PROGRA~1\BLUECF~1\ThirdParty-5.x\platforms\mingw_w64GccDPInt32\lib\libstack_trace.dll
        ZN4Foam5error10printStackERNS_7OstreamE [0x10b1c88+0x218]
                 module: C:\PROGRA~1\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
        ZN4Foam5error5abortEv [0xe65b5d+0x12d]
                 module: C:\PROGRA~1\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
        ZN4Foam17faceAreaWeightAMIINS_14PrimitivePatchINS_4faceENS_7SubListERKNS_5FieldINS_6VectorIdEEEES6_EESA_E14calcAddressingERNS_4ListINS_11DynamicListIiLj0ELj2ELj1EEEEERNSC_INSD_IdLj0ELj2ELj1EEEEESG_SJ_ii [0x63318052+0x1e2]
                 module: C:\PROGRA~1\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libmeshTools.dll
        ZN4Foam17faceAreaWeightAMIINS_14PrimitivePatchINS_4faceENS_7SubListERKNS_5FieldINS_6VectorIdEEEES6_EESA_E9calculateERNS_4ListINSC_IiEEEERNSC_INSC_IdEEEESF_SI_ii [0x63318cb4+0xf4]
                 module: C:\PROGRA~1\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libmeshTools.dll
        ZN4Foam16AMIInterpolationINS_14PrimitivePatchINS_4faceENS_7SubListERKNS_5FieldINS_6VectorIdEEEES6_EESA_E6updateERKSA_SD_ [0x6330e6a9+0x339]
                 module: C:\PROGRA~1\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libmeshTools.dll
        ZN4Foam16AMIInterpolationINS_14PrimitivePatchINS_4faceENS_7SubListERKNS_5FieldINS_6VectorIdEEEES6_EESA_E20constructFromSurfaceERKSA_SD_RKNS_7autoPtrINS_17searchableSurfaceEEE [0x6330d83a+0x5ba]
                 module: C:\PROGRA~1\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libmeshTools.dll
        ZNK4Foam18cyclicAMIPolyPatch8resetAMIERKNS_16AMIInterpolationINS_14PrimitivePatchINS_4faceENS_7SubListERKNS_5FieldINS_6VectorIdEEEES7_EESB_E19interpolationMethodE [0x632cb8ff+0x71f]
                 module: C:\PROGRA~1\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libmeshTools.dll
        ZNK4Foam18cyclicAMIPolyPatch3AMIEv [0x632c6203+0x93]
                 module: C:\PROGRA~1\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libmeshTools.dll
        (No symbol) [0x402935]
                 module: C:\PROGRA~1\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\moveDynamicMesh.exe
        (No symbol) [0x40c0c7]
                 module: C:\PROGRA~1\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\moveDynamicMesh.exe
        (No symbol) [0x4013f7]
                 module: C:\PROGRA~1\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\moveDynamicMesh.exe
        (No symbol) [0x40152b]
                 module: C:\PROGRA~1\BLUECF~1\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\moveDynamicMesh.exe
        BaseThreadInitThunk [0x7ffc88424034+0x14]
                 module: C:\WINDOWS\System32\KERNEL32.DLL
        RtlUserThreadStart [0x7ffc8a293691+0x21]
                 module: C:\WINDOWS\SYSTEM32\ntdll.dll

My dynamicMeshDict is this -

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh   dynamicMotionSolverFvMesh;

motionSolverLibs ("libfvMotionSolvers.so");

motionSolver    solidBody;

cellZone        rotatingZone;

solidBodyMotionFunction  rotatingMotion;

origin      (0 0.361471 0);
axis        (0 1 0);
omega       15.7; // rad/s

My boundary file in Polymesh is this -

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
/*   Windows 32 and 64 bit porting by blueCAPE: http://www.bluecape.com.pt   *\
|  Based on Windows porting (2.0.x v4) by Symscape: http://www.symscape.com   |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      binary;
    class       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

11
(
    walls
    {
        type            wall;
        inGroups        1(wall);
        nFaces          70000;
        startFace       16156900;
    }
    inlet
    {
        type            patch;
        nFaces          4500;
        startFace       16226900;
    }
    outlet
    {
        type            patch;
        nFaces          4500;
        startFace       16231400;
    }
    blade_one_new_meters
    {
        type            wall;
        inGroups        1(wall);
        nFaces          79142;
        startFace       16235900;
    }
    blade_two_new_meters
    {
        type            wall;
        inGroups        1(wall);
        nFaces          78693;
        startFace       16315042;
    }
    blade_three_new_meters
    {
        type            wall;
        inGroups        1(wall);
        nFaces          79990;
        startFace       16393735;
    }
    upper_hub_new_meters
    {
        type            wall;
        inGroups        1(wall);
        nFaces          9309;
        startFace       16473725;
    }
    lower_hub_new_meters
    {
        type            wall;
        inGroups        1(wall);
        nFaces          10003;
        startFace       16483034;
    }
    tower_new_meters
    {
        type            wall;
        inGroups        1(wall);
        nFaces          28261;
        startFace       16493037;
    }
    AMI1
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          426887;
        startFace       16521298;
        matchTolerance  0.0001;
        transform       noOrdering;
        neighbourPatch  AMI2;
    }
    AMI2
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          428124;
        startFace       16948185;
        matchTolerance  0.0001;
        transform       noOrdering;
        neighbourPatch  AMI1;
    }
)

// ************************************************************************* //

My createPatchDict is this -

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      createPatchDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

// 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       noOrdering;
        }
        constructFrom patches;
        patches (AMI);
    }

    {
        //- Slave side patch
        name            AMI2;
        patchInfo
        {
            type            cyclicAMI;
            matchTolerance  0.0001;
            neighbourPatch  AMI1;
            transform       noOrdering;
        }
        constructFrom patches;
        patches (AMI_slave);
    }

/*
    {
        name inlet;
        patchInfo
        {
            type            patch;
        }
        constructFrom set;
        set inletFaces;
    }
    {
        name outlet;
        patchInfo
        {
            type            patch;
        }
        constructFrom set;
        set outletFaces;
    }
*/
);

// ************************************************************************* //

[eloy_785]

Hello,

I'm having the same AMI ERROR:
[101] --> FOAM FATAL ERROR:
[101] Unable to set source and target faces
[101]
[101] From function void Foam::faceAreaWeightAMI<SourcePatch, TargetPatch>::setNextFaces(Foam::label&, Foam::label&, Foam::label&, const boolList&, Foam::labelList&, const Foam:ynamicList<int>&, bool) 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> >&>; Foam::label = int; Foam::boolList = Foam::List<bool>; Foam::labelList = Foam::List<int>]
[101] in file lnInclude/faceAreaWeightAMI.C at line 287.
[101]
FOAM parallel run aborting


This is my createPatchDict
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location system;
object createPatchDict;
}
pointSync false;
patches
(
{
name cm_fan_region_00_940w_brose_interface_ami;
patchInfo
{
type cyclicAMI;
matchTolerance 1.0E-4;
neighbourPatch cm_fan_region_00_940w_brose_interface_slave_ami;
transform noOrdering;
}
constructFrom patches;
patches (cm_fan_region_00_940w_brose_interface);
}
{
name cm_fan_region_00_940w_brose_interface_slave_ami;
patchInfo
{
type cyclicAMI;
matchTolerance 1.0E-4;
neighbourPatch cm_fan_region_00_940w_brose_interface_ami;
transform noOrdering;
}
constructFrom patches;
patches (cm_fan_region_00_940w_brose_interface_slave);
}
);




Does anyone has a solution or suggestion?


Regards

[nimasam]

check the grid motion with the command

moveMesh


then visualize it with paraview, then you can understand the reason behind the error.

Originally Posted by eloy_785

Hello,

I'm having the same AMI ERROR:
[101] --> FOAM FATAL ERROR:
[101] Unable to set source and target faces
[101]
[101] From function void Foam::faceAreaWeightAMI<SourcePatch, TargetPatch>::setNextFaces(Foam::label&, Foam::label&, Foam::label&, const boolList&, Foam::labelList&, const Foam:ynamicList<int>&, bool) 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> >&>; Foam::label = int; Foam::boolList = Foam::List<bool>; Foam::labelList = Foam::List<int>]
[101] in file lnInclude/faceAreaWeightAMI.C at line 287.
[101]
FOAM parallel run aborting


This is my createPatchDict
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location system;
object createPatchDict;
}
pointSync false;
patches
(
{
name cm_fan_region_00_940w_brose_interface_ami;
patchInfo
{
type cyclicAMI;
matchTolerance 1.0E-4;
neighbourPatch cm_fan_region_00_940w_brose_interface_slave_ami;
transform noOrdering;
}
constructFrom patches;
patches (cm_fan_region_00_940w_brose_interface);
}
{
name cm_fan_region_00_940w_brose_interface_slave_ami;
patchInfo
{
type cyclicAMI;
matchTolerance 1.0E-4;
neighbourPatch cm_fan_region_00_940w_brose_interface_ami;
transform noOrdering;
}
constructFrom patches;
patches (cm_fan_region_00_940w_brose_interface_slave);
}
);




Does anyone has a solution or suggestion?


Regards

[toonco17]

Originally Posted by Elena89

Hi R,

are you using solidBodyMotionFvMesh? Rotation?
To me it sounds like the rotational axis (if it's a rotation) isn't normal to the AMI faces and that's why they move apart.

regards

Elena

Hello, Elena,
would you be so king to tell me what you meant here in detail, please? If you still remember that case, ofc. Having the same problem as mr.R. I'm new to CFD and don't really understand how to how fix that problem with those normals you spoke about. Could you, please, tell me what files to fix and how exactly may I do it? Thanks in advance!
Daria