CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

"Unable to set source and target faces" in pimpleDyMFoam

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 3 Post By Elena89

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 27, 2016, 15:26
Default "Unable to set source and target faces" in pimpleDyMFoam
  #1
New Member
 
R
Join Date: Jan 2016
Posts: 2
Rep Power: 0
rbbarber is on a distinguished road
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??

Last edited by rbbarber; January 27, 2016 at 16:23. Reason: update
rbbarber is offline   Reply With Quote

Old   January 28, 2016, 03:12
Default
  #2
New Member
 
Join Date: Jul 2013
Posts: 3
Rep Power: 13
Elena89 is on a distinguished road
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
Swift, ellenripley and hogsonik like this.
Elena89 is offline   Reply With Quote

Old   January 28, 2016, 20:16
Default
  #3
New Member
 
R
Join Date: Jan 2016
Posts: 2
Rep Power: 0
rbbarber is on a distinguished road
Yes, you're right. I had changed that in a previous simulation but forgot to in this one. Thank you,
R
rbbarber is offline   Reply With Quote

Old   February 28, 2019, 10:33
Default Unable to set source and target faces
  #4
New Member
 
shantanu mandhane
Join Date: Jun 2018
Posts: 10
Rep Power: 8
shantanumandhane is on a distinguished road
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
__________________
Regards!
Shantanu
shantanumandhane is offline   Reply With Quote

Old   February 4, 2021, 05:40
Default
  #5
New Member
 
Rokas Bytautas
Join Date: Jan 2021
Posts: 3
Rep Power: 5
rokasbytautas is on a distinguished road
Have you managed to solve this?
rokasbytautas is offline   Reply With Quote

Old   March 2, 2021, 02:30
Default Unable to set source and target faces
  #6
New Member
 
Pratik
Join Date: Apr 2020
Location: Pune
Posts: 3
Rep Power: 6
Pratik Wani is on a distinguished road
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"

Code:
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 -

Code:
/*--------------------------------*- 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 -

Code:
/*--------------------------------*- 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 -

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

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

Last edited by Pratik Wani; March 3, 2021 at 03:41.
Pratik Wani is offline   Reply With Quote

Old   May 26, 2023, 10:51
Default Same AMI problem
  #7
New Member
 
Helloy
Join Date: Jul 2012
Posts: 4
Rep Power: 14
eloy_785 is on a distinguished road
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
eloy_785 is offline   Reply With Quote

Old   July 29, 2023, 07:29
Default
  #8
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,267
Blog Entries: 1
Rep Power: 25
nimasam is on a distinguished road
check the grid motion with the command

moveMesh


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





Quote:
Originally Posted by eloy_785 View Post
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
__________________
My Personal Website (http://nimasamkhaniani.ir/)
Telegram channel (https://t.me/cfd_foam)
nimasam is offline   Reply With Quote

Old   August 11, 2023, 11:18
Default
  #9
New Member
 
Darya Ilina
Join Date: Aug 2023
Posts: 1
Rep Power: 0
toonco17 is on a distinguished road
Quote:
Originally Posted by Elena89 View Post
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
toonco17 is offline   Reply With Quote

Reply

Tags
arbitrary mesh interface, pimpledymfoam


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Building OpenFOAM1.7.0 from source ata OpenFOAM Installation 46 March 6, 2022 14:21
AMI speed performance danny123 OpenFOAM 21 October 24, 2020 05:13
OpenFOAM without MPI kokizzu OpenFOAM Installation 4 May 26, 2014 10:17
centOS 5.6 : paraFoam not working yossi OpenFOAM Installation 2 October 9, 2013 02:41
mixerVesselAMI2D's mass is not balancing sharonyue OpenFOAM Running, Solving & CFD 6 June 10, 2013 10:34


All times are GMT -4. The time now is 09:56.