[dogan]

Hi everyone,

I am having problems with stitchMesh command in OF 2.1.x. My case a centrifugal pump which has 2 interfaces between the rotating part and the stationary part. I want to run a case with MRFSimpleFoam solver. Mesh was in .msh format, and i converted it to openFoam by using fluent3DMeshToFoam. after that, i used topoSet and i obtained the constant>polyMesh>sets directory with rotor file in it. afterwards, what i know is, i need to stitch the interfaces, which indicates the connection between rotating and the stationary parts, but the problem is, the number of the faces of two patches are not identical, i think because of that the stitchMesh command gives me the following error message:


FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) in "/home/gulacd2k/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/gulacd2k/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so"
#2 Foam::enrichedPatch::calcCutFaces() const in "/home/gulacd2k/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libdynamicMesh.so"
#3 Foam::enrichedPatch::cutFaces() const in "/home/gulacd2k/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libdynamicMesh.so"
#4 Foam::slidingInterface::coupleInterface(Foam:oly TopoChange&) const in "/home/gulacd2k/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libdynamicMesh.so"
#5 Foam:olyTopoChanger::topoChangeRequest() const in "/home/gulacd2k/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libdynamicMesh.so"
#6 Foam:olyTopoChanger::changeMesh(bool, bool, bool, bool) in "/home/gulacd2k/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libdynamicMesh.so"
#7
in "/home/gulacd2k/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/bin/stitchMesh"
#8 __libc_start_main in "/lib64/libc.so.6"
#9
at /home/abuild/rpmbuild/BUILD/glibc-2.15/csu/../sysdeps/x86_64/elf/start.S:116


i couldn't find what should i do,
thanks in advance for your helps
Dogan

[dogan]

Hi,


i am working on a centrifugal pump geometry, and i want to stitch two interfaces with the stitchMesh command. the master interface is GEOM-SIDE-2 has 5481 faces, and the slave interface has 3248 faces.

once i run the command "stitchMesh GEOM-SIDE-2 GEOM-SIDE-1" the following error message comes up:


Create mesh for time = 0

Coupling partially overlapping patches GEOM-SIDE-2 and GEOM-SIDE-1
Resulting internal faces will be in faceZone GEOM-SIDE-2GEOM-SIDE-1CutFaceZone
Any uncovered faces will remain in their patch
Adding pointZone GEOM-SIDE-2GEOM-SIDE-1CutPointZone at index 0
Adding faceZone GEOM-SIDE-2GEOM-SIDE-1MasterZone at index 0
Adding faceZone GEOM-SIDE-2GEOM-SIDE-1SlaveZone at index 1
Adding faceZone GEOM-SIDE-2GEOM-SIDE-1CutFaceZone at index 2
Sliding interface parameters:
pointMergeTol : 0.05
edgeMergeTol : 0.01
nFacesPerSlaveEdge : 5
edgeFaceEscapeLimit : 10
integralAdjTol : 0.05
edgeMasterCatchFraction : 0.4
edgeCoPlanarTol : 0.8
edgeEndCutoffTol : 0.0001
Reading all current volfields
Reading volScalarField p
Reading volScalarField nut
Reading volScalarField k
Reading volScalarField omega
Reading volVectorField U


--> FOAM FATAL ERROR:
Duplicate point found in cut face. Error in the face cutting algorithm for global face 4(476034 466684 466686 476036) local face 4(0 1 2 3)
Slave size: 3248 Master size: 5481 index: 0.
Face: 5(476034 466684 466686 338175 476036)
Cut face:
101
(
476034
466684
466686
644233
338166
466686
466684
559434
559436
........




Those interfaces indicates the connection between rotor and stator of the pump, and my aim is to run a simulation with MRFSimpleFoam after stitching them. I am using OF 2.1.x. I tried whatever i know, but i couldn't find a solution. i hope some of you can help me in this matter.

thanks
Dogan

[wyldckat]

Greetings Dogan,

Without a test case, I can't figure this out myself. But my suggesting is that you try using the "-partial" option with stitchMesh.

There are a few examples given in the following thread and the respective solution for those examples: when can stitchMesh be used?
From it you should be able to derive some information that might help you get closer to the solution.

If you're still not able to figure it out, please create a simple case that can lead to a similar error that you are having and share it with us! This way I or anyone else can help you figure this out!

Best regards,
Bruno