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

[snappyHexMesh] snappyHexMesh on slice of cylinder

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 13, 2018, 11:46
Default snappyHexMesh on slice of cylinder
  #1
New Member
 
Loek
Join Date: Oct 2018
Posts: 14
Rep Power: 4
Loekatoni is on a distinguished road
Dear users of cfd-online.com


Here I m again, with a different question. So I want to make a bubble inside a channel. I succeeded in making the channel in blockMesh (bM) and I have manually made a bubble.stl file. Now it is axisymmetric, so my slices goes from -2.5 until 2.5 degrees (for my stl file, the angle is exactly the same) to make the system not have many cells.
_______________________________________________

The line I do are as follows:


surfaceFeatureExtract
blockMesh
decomposePar
mpirun -np 4 snappyHexMesh -parallel
reconstructParMesh -constant
reconstructPar -latestTime
_______________________________________________


But it stops during the parallel run of snappyHexMesh (sHM). the error I get is as follows:


Marked for refinement due to explicit features : 106 cells.
Determined cells to refine in = 0.02 s
Selected for feature refinement : 106 cells (out of 1920)
[0] hexRef8 : Dumping cell as obj to "/home/s164072/OpenFOAM/s164072-v1806/run/tutorials/Loek/graduation/scaleR0005/processor0/cell_114.obj"
[0]
[0]
[0] --> FOAM FATAL ERROR:
[0] cell 114 of level 0 does not seem to have 8 points of equal or lower level
cellPoints:6(268 269 284 237 238 253)
pointLevels:6{0}
[0]
[0]
[0] From function Foam::labelListList Foam::hexRef8::setRefinement(const labelList&, Foam:olyTopoChange&)
[0] in file polyTopoChange/polyTopoChange/hexRef8/hexRef8.C at line 3679.
[0]
FOAM parallel run aborting
[0]
[0] #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[0] #1 Foam::error::abort() at ??:?
[0] #2 Foam::hexRef8::setRefinement(Foam::List<int> const&, Foam:olyTopoChange&) at ??:?
[0] #3 Foam::meshRefinement::refine(Foam::List<int> const&) at ??:?
[0] #4 Foam::meshRefinement::refineAndBalance(Foam::strin g const&, Foam::decompositionMethod&, Foam::fvMeshDistribute&, Foam::List<int> const&, double) at ??:?
[0] #5 Foam::snappyRefineDriver::featureEdgeRefine(Foam:: refinementParameters const&, int, int) at ??:?
[0] #6 Foam::snappyRefineDriver::doRefine(Foam::dictionar y const&, Foam::refinementParameters const&, Foam::snapParameters const&, bool, bool, Foam::dictionary const&) at ??:?
[0] #7 ? at ??:?
[0] #8 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
[0] #9 ? at ??:?
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.


__________________________________________________ __________


To me it seems that the error is about my wedge mesh. I have points overlapping in my bM, so to get a triangular mesh, or slice of the channel. These cells have an input of 8 vertices, but 2 overlap, which could result in openFOAM seeing it as a cell with 6 vertices.


My question: what can I do against this problem?


Sincerely,
Loek
Loekatoni is offline   Reply With Quote

Old   November 14, 2018, 04:20
Default
  #2
Senior Member
 
Zander Meiring
Join Date: Jul 2018
Posts: 125
Rep Power: 4
yambanshee is on a distinguished road
I would avoid the problem by reconstructing your blockmesh with a diamond shape at the wedge (see attached image)
Attached Images
File Type: jpg index.jpg (25.3 KB, 7 views)
yambanshee is offline   Reply With Quote

Old   November 14, 2018, 06:53
Default
  #3
New Member
 
Loek
Join Date: Oct 2018
Posts: 14
Rep Power: 4
Loekatoni is on a distinguished road
Dear yambanshee,


Thanks for the reply.


I applied your strategy, but slightly different. Your idea gives two gridcells in the direction where I want just one. So i did it as shown in the picture. Now lets see if sHM will work!


EDIT: so this seems to work fine for this kind of Mesh
Attached Images
File Type: jpg IMG_4818.jpg (86.6 KB, 9 views)

Last edited by Loekatoni; November 16, 2018 at 05:29. Reason: Fixed a typo and conclusion
Loekatoni is offline   Reply With Quote

Reply

Tags
snappyhexmesh, stl file, vertices, wedge

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
Drag force coefficient too low for a flow past cylinder at Re= 1e05 Scabbard STAR-CCM+ 2 June 5, 2020 15:44
Flow past a cylinder at Re 1e05 using LES, drag force coefficient is to low Scabbard Main CFD Forum 21 June 19, 2018 14:58
Lift Cylinder in Liquid using ANSYS MikeShuo Main CFD Forum 0 October 24, 2016 02:46
flow over a cylinder urgent! kevin FLUENT 8 August 11, 2015 14:00
[blockMesh] Specifying boundary faces failes in blockMesh blaise OpenFOAM Meshing & Mesh Conversion 0 May 10, 2010 04:56


All times are GMT -4. The time now is 10:53.