|
[Sponsors] | |||||
|
|
|
#1 |
|
New Member
Join Date: Aug 2013
Posts: 4
Rep Power: 14 ![]() |
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 |
|
|
|
|
|
|
|
|
#2 |
|
New Member
Henning Kuchenbuch
Join Date: Jan 2018
Location: Hannover
Posts: 8
Rep Power: 10 ![]() |
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 |
| Thread Tools | Search this Thread |
| Display Modes | |
|
|
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 |