[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;
    }
*/
);

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