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

[blockMesh] blockMesh not producing straight lines when using edgeGrading

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 1, 2019, 09:21
Default blockMesh not producing straight lines when using edgeGrading
  #1
New Member
 
Join Date: Feb 2012
Posts: 14
Rep Power: 14
ivan_CFD is on a distinguished road
Hi everyone,


I would like to request help with an issue when using edgeGrading in blockMesh. I attach some pictures to illustrate it.


I have a mesh consisting of two blocks, one of them is rectangular and the other one is trapezoidal (see blocks.jpg). There is only one cell in the z direction. The rectangular block has a uniform spacing in the vertical direction (y), whereas the trapezoidal block has a grading along the vertical direction.


To avoid jumps in cell size, I want the vertical cell size at each side of the interface between the two blocks to be the same. In other words, I want the first cell size of the trapezoidal block to be equal to the last cell size of the rectangular block (in the y direction). To achieve this, I use edgeGrading on the trapezoidal block, computing the respective expansion ratio that is needed at the left (inlet) and right (outlet) vertical edges of the block.



The problem is the following: at the left and right edges of the trapezoidal block, where the grading is imposed in the blockMeshDict (also in attachment), the grading is working exactly as expected and the cell size is matching (see mesh_inlet.jpg and mesh_outlet.jpg). On the other hand, in the vertical mesh divisions inside the block this is not being respected. This implies that for instance, the first horizontal grid line of the trapezoidal block is not straight but has a small curvature, causing a jump in cell size at the interface of the two blocks (see mesh_middle.jpg). This jump is maximum at the middle of the block and decreases towards the left and right vertical edges, where the grading is imposed. The longer is the domain, the more this problem accentuates.



I would expect blockMesh to join the distribution of points at the left and right vertical edges with straight lines, but for some reason this is not the case. It is not an optical artifact as can be confirmed by looking at the y coordinate of the points.



I have not been able to find another thread related to this problem. Only the following one https://www.cfd-online.com/Forums/op...s-geomety.html in which a similar problem is reported.



I have tried blockMesh with OF versions 3 to 6 with the same result.


Many thanks in advance.


Ivan
Attached Images
File Type: jpg blocks.jpg (22.6 KB, 28 views)
File Type: jpg mesh_inlet.jpg (88.2 KB, 29 views)
File Type: jpg mesh_middle.jpg (64.8 KB, 27 views)
File Type: jpg mesh_outlet.jpg (63.0 KB, 22 views)
Attached Files
File Type: txt blockMeshDict.txt (1.3 KB, 7 views)
ivan_CFD is offline   Reply With Quote

Reply

Tags
bent, blockmesh, edgegrading, lines, straight

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
[Tutorials] Tutorial of how to plot residuals ! wolle1982 OpenFOAM Community Contributions 171 February 20, 2024 02:55
[blockMesh] blockmesh arc, no curves, no arc, only straight lines heavy_user OpenFOAM Meshing & Mesh Conversion 9 January 26, 2017 20:48
GMSH producing lines through mesh APCoding Mesh Generation & Pre-Processing 1 January 23, 2017 10:46
[ANSYS Meshing] How to align nodes in straight (vertical) lines Judith ANSYS Meshing & Geometry 5 July 25, 2014 04:55
Plotting Residuals ata OpenFOAM 12 July 8, 2011 02:50


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