CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Pre-Processing (https://www.cfd-online.com/Forums/openfoam-pre-processing/)
-   -   Creating baffles in foam-extend (https://www.cfd-online.com/Forums/openfoam-pre-processing/187808-creating-baffles-foam-extend.html)

MarcelK May 15, 2017 15:57

Creating baffles in foam-extend
 
1 Attachment(s)
I am currently working on a three-dimensional channel flow problem with a thin structure, which is finally intended to be integrated into an FSI simulation (that's why I need 0/motionU and constant/dynamicMeshDict).

The thin structure should be located at the face where block 0 and block 2 of my blockMeshDict intersect. Since the structure is very thin, I want to use createBaffles to convert internal faces into a boundary such that I am able to address the two sides of the thin structure separately (such that I am able to request pressure, for example, in the FSI calculation later on). My workflow for this is as follows:

Code:

m4 constant/polyMesh/blockMeshDict.m4 > constant/polyMesh/blockMeshDict
blockMesh
setSet -batch setBatch.batch
setsToZones
createBaffles fsBaffleLeft consoleLeft -overwrite
createBaffles fsBaffleRight consoleRight -overwrite

Apparently, the step to create the baffles does not work since the second invocation of createBaffles overwrites the result of the first invocation. Leaving the -overwrite option apart does not work either. My setBatch.batch file looks as follows:

Code:

faceSet fsBaffleLeft new cellToFace block0 all
faceSet fsBaffleLeft subset cellToFace block2 all
cellSet fsBaffleLeftMasterCells new cellToCell block0

faceSet fsBaffleRight new cellToFace block0 all
faceSet fsBaffleRight subset cellToFace block2 all
cellSet fsBaffleRightMasterCells new cellToCell block2
quit

When using the createBaffles tool in OpenFOAM, it is possible to specify a master and slave patch via a createBafflesDict. Unfortunately, I have not yet found an option in the createBaffles from foam-extend to do this. Any help on this would be much appreciated.

CFD-Henning February 21, 2018 12:05

Creating interface in FE-4.0
 
Hey Marcel,
did you solve the problem? I am sorry to come up as a thread necromancer, but I'm struggeling with the same problem.
I'm a foam-newbie and work also on a FSI-case (Foam-Extend 4.0). Thereby I want to analyze the heave-movement of a floating object in a wave tank. So I want to insert this object as a patch with a wall-bc into an internal part of the pre-meshed domain.

I went as following:
-Meshing with blockMesh
-define cellSet "c0" with the geometrie of the floating Object
-define a faceSet "interface_zone" out of c0 (just all sides)
-create cellZones
-replace pre-defined empty interface-patch with new patch (createBaffles interface_zone interface -overwrite)

Unfortunately this is not working, coming up with the following checkMesh-log

Code:

Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
    all points:          7007
    live points:          7007
    all faces:            17794
    live faces:          17794
    internal faces:      14606
    cells:                5400
    boundary patches:    6
    point zones:          0
    face zones:          1
    cell zones:          2

Overall number of cells of each type:
    hexahedra:    5400
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:    0

Checking topology...
    Boundary definition OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
  *Number of regions: 2
    The mesh has multiple regions which are not connected by any face.
  <<Writing region information to "0/cellToRegion"
Nuumber of cells per region:
    0    5340
    1    60


Checking patch topology for multiply connected surfaces ...
    Patch              Faces    Points  Area [m^2]  Surface topology                 
    inlet              60      77      0.6        ok (non-closed singly connected) 
    outlet              60      77      0.6        ok (non-closed singly connected) 
    bottom              900      1001    9          ok (non-closed singly connected) 
    atmosphere          900      1001    9          ok (non-closed singly connected) 
    frontBack          1080    1274    10.8        ok (non-closed singly connected) 
    interface          188      96      1.88        multiply connected (shared edge) 
  <<Writing 96 conflicting points to set nonManifoldPoints

Checking geometry...
    This is a 3-D mesh
    Overall domain bounding box (0 -0.4 0) (9 0.2 1)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Mesh (non-empty, non-wedge) dimensions 3
    Boundary openness (5.44142e-19 3.55216e-16 6.00188e-17) Threshold = 1e-06 OK.
    Max cell openness = 1.73472e-16 OK.
    Max aspect ratio = 1 OK.
    Minumum face area = 0.01. Maximum face area = 0.01.  Face area magnitudes OK.
    Min volume = 0.001. Max volume = 0.001.  Total volume = 5.4.  Cell volumes OK.
    Mesh non-orthogonality Max: 0 average: 0 Threshold = 70
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 2.51968e-12 OK.

Mesh OK.

End

The interface-patch switched from 0 to 188 faces, whereby the faceSet of the floating Object has got only the half of it - 94 faces. It seems that there are issues with the direction of the faces, as I didn't define a master- and a slave-side.
What am I doing wrong and how is the best-practice-way of getting this solved in foam-extend? Btw. do I have to go with subsetMesh?
I would be very happy if you, or anyone in the community could help me with this .
Thanks in advance,
Henning


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