[frobaux]

I everyone,

I already read pretty much all tread about meshToMesh (If I missed some, let me know =) ) but there is something I don't understand. I want to map a field which is defined on mesh externalPhiRegions[0] on a mesh completely included, fluidRegions[i].
I use the meshToMesh toolbox.

meshToMesh mapper(
                externalPhiRegions[0], 
                fluidRegions[i], 
                meshToMesh::imCellVolumeWeight,
                patchMap,        // patchMap,
                cuttingPatches,  // cuttingPatches
                true 
            );

If I got it right, cuttingPatches are all the patches which needs to be interpolated and are completely inside the external region. So they should be defined in fluidRegions[i], right?
The problem is I got a Fatal error as the cuttingPatch is not inside the externalPhiRegions[0] mesh (which name is of_GlobMesh).

--> FOAM FATAL ERROR: 
Unable to find patch 'inlet' in mesh 'of_GlobMesh'.  Available patches include:
6
(
front
back
LeftWall
BottomWall
RightWall
athmosphere
 )
 From function void Foam::meshToMesh::constructFromCuttingPatches(const Foam::word&, const Foam::word&, const Foam::HashTable<Foam::word>&, const wordList&, bool)
    in file meshToMesh/meshToMesh.C at line 789.

Why does the code need the cutting patches to be in the source mesh??

[frobaux]

I tried to do it manually, using the mapFieldsPar utility

(which uses the meshToMesh library, where mapFields uses the meshToMesh0)
I obtain the same output (this time of_boxMesh3 is the source mesh)

mapFieldsPar -sourceRegion of_boxMesh3 -targetRegion Hsl35 -fields '(Phi_interp)'  ./

output:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Source: "." ""
Source region: of_boxMesh3
Target: "/home/robaux/These/HPC/IBMfs/CylinderVenugopal/TestConversionOF/nx60CFL200" "mergedboxMesh3"
Target region: Hsl35

Create databases as time
I/O    : uncollated
Case   : /home/robaux/These/HPC/IBMfs/CylinderVenugopal/TestConversionOF/nx60CFL200/mergedboxMesh3
nProcs : 1
Overriding DebugSwitches according to controlDict
    meshToMesh      1;

Overriding DebugSwitches according to controlDict
    meshToMesh      1;


Source time: 6.7
Target time: 6.7

Create meshes

Source mesh size: 1800    Target mesh size: 1341


Creating and mapping fields for time 6.7



--> FOAM FATAL ERROR: 
Unable to find patch 'inlet' in mesh 'of_boxMesh3'.  Available patches include:
7
(
front
back
extInterpLeft
extInterpBottom
extInterpRight
extInterpTop
wallbox
)


    From function void Foam::meshToMesh::constructFromCuttingPatches(const Foam::word&, const Foam::word&, const Foam::HashTable<Foam::word>&, const wordList&, bool)
    in file meshToMesh/meshToMesh.C at line 789.

FOAM exiting

I don't understand why the cuttingPatches should be in the sourceMesh. Is there a bug here?

[frobaux]

Well, I think I may have find the answer: A problem in the code of my OF version 1712, as well as in 1812:
lines are in src/sampling/meshToMesh/meshToMesh.C

 cuttingPatches_.setSize(cuttingPatches.size());

     forAll(cuttingPatches_, i)

     {

         const word& patchName = cuttingPatches[i];

         label cuttingPatchi = srcBm.findPatchID(patchName);



         if (cuttingPatchi == -1)

         {

             FatalErrorInFunction

                 << "Unable to find patch '" << patchName

                 << "' in mesh '" << srcRegion_.name() << "'. "

                 << " Available patches include:" << srcBm.names()

                 << exit(FatalError);

         }



         cuttingPatches_[i] = cuttingPatchi;

     }

I tried to replace the "src" by "tgt" and it seems to work like a charm.
In OpenFoam v6, corresponding line takes the mesh target, as well as in foam-extend 3.2.


Should I do something about it? Post that somewhere or inform someone?

Thanks a lot