CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[snappyHexMesh] Mesh embraces only one of four specified domain sub-volumes

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 30, 2020, 18:28
Default Mesh embraces only one of four specified domain sub-volumes
  #1
New Member
 
J.D. Wilson
Join Date: Nov 2020
Location: Edmonton, Canada
Posts: 15
Rep Power: 2
JayDeeUU is on a distinguished road
I'm new to OpenFoam, and attempting to simulate buoyancy-driven internal circulation within a vessel whose geometry is that of a round (but flat-bottomed) flask with a neck. I want to impose a steady but non-uniform temperature distribution on the flow boundary, and have inferred (rightly or wrongly) that to accomplish that I have to build the computational mesh in such a manner that distinct boundary patches or regions are recognised.

I've begun by dividing the vessel into four parts, with the divisions occurring by slicing perpendicular to the vertical axis. Each part is azimuthally-symmetric: a lower partial sphere, an upper partial sphere, a lower cylinder and an upper cylinder. I also created two discs, one configured for the bounding bottom plane, the other for the top.

I created these as .stl files in FreeCAD, converted them to ascii file format, and renamed the parts in the first line (e.g. "solid cyl1") and last line ("endsolid cyl1") of each .stl file. I then created the blockMesh, ran "surfaceFeatures" and finally "snappyHexMesh" (within which I set "scale 0.001;" to rescale FreeCAD's [mm] to OpenFoam's expected [m] dimensions).

The problem is that whatever specification I make for "locationInMesh" within snappyHexMeshDict, the outcome is that I see only one of the four parts (a different part according to where I place that "locationInMesh"). I've been following Calum Douglas's www.youtube.com/watch?v=ObsFQUiVi1U to some extent, so have experimented with both a single .stl containing all parts, or, listing the parts individually. I get the same outcome either way. The attached image (from MeshLab) shows some of the "parts" as represented by the .stl files (other parts suppressed for display).

I'm running OpenFoam Version 8 under Ubuntu 20.04. As an aside, I have made low resolution simulations of this flow with my own Fortran code (collocated grid, SIMPLE method) - which however is grossly compromised by its use of a uniform mesh (of cubes) resulting in a stepwise representation of the flask's wall, and inadequate resolution of the near-wall zone.
Attached Images
File Type: png illustration.png (102.4 KB, 13 views)

Last edited by JayDeeUU; December 2, 2020 at 17:23.
JayDeeUU is offline   Reply With Quote

Old   December 9, 2020, 15:46
Default
  #2
New Member
 
J.D. Wilson
Join Date: Nov 2020
Location: Edmonton, Canada
Posts: 15
Rep Power: 2
JayDeeUU is on a distinguished road
The cause of this problem was that each of the individual .stl files was a closed shell. This whatever value was used for "locationInMesh" (set in snappyHexMeshDict), only one of my (wanted) four sub-volumes was meshed.
JayDeeUU is offline   Reply With Quote

Old   December 9, 2020, 15:50
Default
  #3
New Member
 
J.D. Wilson
Join Date: Nov 2020
Location: Edmonton, Canada
Posts: 15
Rep Power: 2
JayDeeUU is on a distinguished road
The cause of this problem was that each of my four .stl files described a closed shell. Therefore the mesh covered only the subvolume for which the "locationInMesh" setting was an interior point.
JayDeeUU is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[snappyHexMesh] SnappyHexMesh/splitMeshRegion : region1 in zone "-1" GuiMagyar OpenFOAM Meshing & Mesh Conversion 1 September 10, 2019 10:59
Gambit problems Althea FLUENT 22 January 4, 2017 04:19
sliding mesh problem in CFX Saima CFX 45 September 22, 2015 11:53
[snappyHexMesh] how to mesh a domain splitted in two? aCFDuser OpenFOAM Meshing & Mesh Conversion 6 August 19, 2015 09:27
injection problem Mark New FLUENT 0 August 4, 2013 02:30


All times are GMT -4. The time now is 18:07.