|
[Sponsors] |
May 15, 2017, 16:57 |
Creating baffles in foam-extend
|
#1 |
New Member
Join Date: Aug 2013
Posts: 4
Rep Power: 13 |
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 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 |
|
February 21, 2018, 13:05 |
Creating interface in FE-4.0
|
#2 |
New Member
Henning Kuchenbuch
Join Date: Jan 2018
Location: Hannover
Posts: 8
Rep Power: 8 |
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 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 |
|
Tags |
blockmeshdict, createbaffles, foam-extend, internal faces, setset |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Other] mesh airfoil NACA0012 | anand_30 | OpenFOAM Meshing & Mesh Conversion | 13 | March 7, 2022 18:22 |
[CGNS] CGNS converters available | mbeaudoin | OpenFOAM Meshing & Mesh Conversion | 137 | December 14, 2018 05:20 |
[blockMesh] non-orthogonal faces and incorrect orientation? | nennbs | OpenFOAM Meshing & Mesh Conversion | 7 | April 17, 2013 06:42 |
Problem with rhoSimpleFoam | matteo_gautero | OpenFOAM Running, Solving & CFD | 0 | February 28, 2008 07:51 |
[Gmsh] Import gmsh msh to Foam | adorean | OpenFOAM Meshing & Mesh Conversion | 24 | April 27, 2005 09:19 |