CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[snappyHexMesh] snappyHexMesh on sharp corners (cube)

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

LinkBack Thread Tools Search this Thread Display Modes
Old   June 4, 2015, 23:47
Default snappyHexMesh on sharp corners (cube)
New Member
Join Date: Jan 2012
Posts: 24
Rep Power: 13
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, 113 views)
File Type: png fig2.png (47.7 KB, 161 views)
Attached Files
File Type: zip (15.0 KB, 15 views)
Regis_ is offline   Reply With Quote


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] Error defining boundary layer around cube snappyHexMesh crizpi21 OpenFOAM Meshing & Mesh Conversion 5 October 16, 2021 10:56
[ICEM] Skewed elements at sharp corners Julian121 ANSYS Meshing & Geometry 3 January 11, 2019 03:55
Does openfoam quite sensitive to the sharp corners of the geometry? jiaodanuma OpenFOAM 0 August 16, 2016 07:53
Unphysical peaks in the pressure field at sharp corners why? jcamilleri STAR-CCM+ 0 March 5, 2016 05:35
[snappyHexMesh] SnappyHexMesh and sharp corners madad2005 OpenFOAM Meshing & Mesh Conversion 2 September 7, 2009 09:26

All times are GMT -4. The time now is 22:00.