CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Pre-Processing (
-   -   Create an internal face (

aylalisa June 19, 2013 09:42

Create an internal face
Hi to All,

I simulate a flow in a channel.

-fluid: water
-transport model: Newtonian
-very small velocity at inlet face that results in a laminar flow in the whole channel

The channel has a length of 3,1m (x-direction, width = y-direction = 0,250m, height = z-direction = 0,02m).

Attempt 1

I would like to define an internal face in the blockMeshDict to be able to map the final velocity field of that internal face on the inlet face of another channel model as boundary condition with help of mapFields utility.

The channel consists of three blocks. The desired internal face connects block 1 (delta x = 0m - 1m = 1m) with block 2 (delta x = 1m - 2,5m = 1,5m).

If I add two additional patches (outlet_block1, inlet_block2) in the boundary () subdictionary and list them up in subdict mergePatchPairs () I get an error message after starting blockMesh.


Trying to specify a boundary face 4(0 3 7 4) on the face on cell 0 which is either an internal face or already belongs to some other patch. This is face 0 of patch 5 named Developing_outlet.

From function polyMesh::setTopology
const cellShapeList& cellsAsShapes,
const faceListList& boundaryFaces,
const wordList& boundaryPatchNames,
labelList& patchSizes,
labelList& patchStarts,
label& defaultPatchStart,
label& nFaces,
cellList& cells
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 360.

FOAM aborting
Attempt 2

I've created the velocity field of that internal face with sample utility --> folder surfaces is created and contains for each time step the required data. But then I am not able to map that data on the inlet boundary of the other model with help of mapFields.

MapFields maps the complete volume velocity field from model 1 to model 2. But I would like to map only the velocity field of the internal surface to the inlet surface (inlet boundary condition) of model 2.

Can anybody help???


mturcios777 June 19, 2013 12:02

I don't believe you can do such a thing with blockMesh. Why not define a faceZone (using topoSet or similar tools), then sample off the faceZone and map. Though I'm not sure how I would map from a faceZone to a boundary; I need to think harder about this...

All times are GMT -4. The time now is 16:48.