CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Native Meshers: snappyHexMesh and Others

snappyHexMesh on sharp corners (cube)

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

LinkBack Thread Tools Display Modes
Old   June 4, 2015, 23:47
Default snappyHexMesh on sharp corners (cube)
New Member
Join Date: Jan 2012
Posts: 14
Rep Power: 7
Regis_ is on a distinguished road
I'm trying to mesh a *really simple* geometry with snappyHexMesh but I'm getting weird results upon viewing it in Paraview.
Before getting into details, I should say that I'm using OpenFOAM 2.0.1. I know this is a relatively old release, but right now I need to use some libraries that currently doesn't work with very recent OpenFOAM releases and I won't have time to port them in the next couple of months. Since what I'm trying to achieve is very simple, 2.0.1 might be enough for now.

This is the geometry I'm looking for: a 3000x3000x1000m domain with a cube on the floor, centered. Cube is 150m in length. The bottom part of the cube should lie on the domain "floor". The cube is just a representation of the geometry I'm going to have there (yes, the geometry is still relatively simple and will contain sharp corners).

See the attached images of what I'm obtaining. The translucent cube is where my mesh should be; it represents the 150x150x150 cube. It is like the mesh is snapping to something smaller than my geometry. And since the cube geometry I'm getting is smaller that what I want, the cube end up not touching the floor and I get cells in between the cube and floor.

These are the steps I'm taking:

1. run blockMesh on the outside domain

2. use decomposePar (scotch method)

3. use surfaceFeatureExtract by using the following command:
surfaceFeatureExtract -includedAngle 150 -writeObj constant/triSurface/cube_ig.stl cubeSFE
(in OpenFOAM 2.0.1 we don't need the dictionary)

4. run snappyHexMesh (now using ptscotch in decomposeParDict)
mpirun -np 32 snappyHexMesh -noFunctionObjects -overwrite -parallel

5. Reconstruct the mesh, remove processor* dirs, decompose again (scotch method)

Some observations:
- The STL file looks ok. I have used surfaceCheck to verify. It is so simple that I had even verified it manually.

- I'm running checkMesh as I go. Everything looks fine. The cube boundary is being added as expected.

- I'm aware that the 5th step can be replaced by using changeDictionary, but I still haven't got to that point yet (any tips are appreciated, though).

- Since my geometry can be meshed with just regular hexahedrons, I tried turning off the snap feature on snappyHexMeshDict. But that didn't really help. It doesn't matter if is on or off, still doesn't work with my current settings

- I did notice I'm getting a warning during sHM, but don't really know what it means:
- --> FOAM Warning : Displacement (0 0 0.694480403) at mesh point 265 coord (1425.69448 1460 149.30552) points through the surrounding patch faces
Smoothing displacement ...
I have attached the dictionaries and a simple script that goes through steps 1--5 outlined above. The output of the script is the defective mesh.

Any ideas on what I'm missing?

Attached Images
File Type: jpg fig1.jpg (50.8 KB, 64 views)
File Type: png fig2.png (47.7 KB, 80 views)
Attached Files
File Type: zip (15.0 KB, 6 views)
Regis_ is offline   Reply With Quote


Thread Tools
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 On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
SnappyHexMesh: Layer on corners Pj. OpenFOAM Native Meshers: snappyHexMesh and Others 9 July 17, 2014 04:26
Meshing sharp corners - mesh for a MVG fportela Mesh Generation & Pre-Processing 0 March 17, 2013 12:33
Sharp corners Ford Prefect Main CFD Forum 5 January 31, 2012 13:30
SnappyHexMesh and sharp corners madad2005 OpenFOAM 2 September 7, 2009 09:26
Flow near sharp corners Harish Main CFD Forum 4 February 21, 2007 22:55

All times are GMT -4. The time now is 06:46.