CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Native Meshers: blockMesh (
-   -   stitching on arcs (

RBJ March 15, 2010 10:48

stitching on arcs
1 Attachment(s)
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


#0  Foam::error::printStack(Foam::Ostream&) in "/uhome/koldewei/OpenFOAM/OpenFOA
#1  Foam::error::abort() in "/uhome/koldewei/OpenFOAM/OpenFOAM-1.6/lib/linux64Gc
#2  Foam::enrichedPatch::calcCutFaces() const in "/uhome/koldewei/OpenFOAM/OpenF
#3  Foam::enrichedPatch::cutFaces() const in "/uhome/koldewei/OpenFOAM/OpenFOAM-
#4  Foam::slidingInterface::coupleInterface(Foam::polyTopoChange&) const in "/uh
#5  Foam::polyTopoChanger::topoChangeRequest() const in "/uhome/koldewei/OpenFOA
#6  Foam::polyTopoChanger::changeMesh(bool, bool, bool, bool) in "/uhome/koldewe
#7  main in "/uhome/koldewei/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDP
#8  __libc_start_main in "/lib64/"
#9  _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.
    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

All times are GMT -4. The time now is 21:47.