CFD Online Discussion Forums - [blockMesh] stitching on arcs

Dear All,

For a multi-phase flow simulation I have a spider-web mesh, as can be seen in the picture. I create the blockmesh file with a matlab script. The mesh exists of a couple of wedge-type blocks in circumferential direction. I have noticed for a test-case simulation that other meshes (like o-type) don't give correct results (not even close!).

At a certain location in axial direction, say at the arrow, I want fluid dropping into the domain. Therefore I thought of making an inlet at a certain location in axial direction. I also want to make a block for this with new nodes and an arc on the downside going through the mesh' radius.

The total outside boundary of the cylinder is called 'freeStream' I have tried to do stitchmesh freestream inletDownSide -partial, but all I got was errors like

Code:

#0  Foam::error::printStack(Foam::Ostream&) in "/uhome/koldewei/OpenFOAM/OpenFOA

M-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"

#1  Foam::error::abort() in "/uhome/koldewei/OpenFOAM/OpenFOAM-1.6/lib/linux64Gc

cDPOpt/libOpenFOAM.so"

#2  Foam::enrichedPatch::calcCutFaces() const in "/uhome/koldewei/OpenFOAM/OpenF

OAM-1.6/lib/linux64GccDPOpt/libdynamicMesh.so"

#3  Foam::enrichedPatch::cutFaces() const in "/uhome/koldewei/OpenFOAM/OpenFOAM-

1.6/lib/linux64GccDPOpt/libdynamicMesh.so"

#4  Foam::slidingInterface::coupleInterface(Foam::polyTopoChange&) const in "/uh

ome/koldewei/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libdynamicMesh.so"

#5  Foam::polyTopoChanger::topoChangeRequest() const in "/uhome/koldewei/OpenFOA

M/OpenFOAM-1.6/lib/linux64GccDPOpt/libdynamicMesh.so"

#6  Foam::polyTopoChanger::changeMesh(bool, bool, bool, bool) in "/uhome/koldewe

i/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libdynamicMesh.so"

#7  main in "/uhome/koldewei/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDP

Opt/stitchMesh"

#8  __libc_start_main in "/lib64/libc.so.6"

#9  _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.

S:116

 

    From function void enrichedPatch::calcCutFaces() const

    in file slidingInterface/enrichedPatch/enrichedPatchCutFaces.C at line 567.

FOAM aborting

Does anyone have any clues on how to perform this?

Of course it would be possible to manually select the nodes on which the inlet should be added, but I'd like to make it more an automatic procedure: change some variables in my Matlab script -> blockMesh > stitchmesh > run simulation