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

[mesh manipulation] Does extrudeMesh affect mesh quality in non extruded cells?

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 13, 2015, 12:04
Default Does extrudeMesh affect mesh quality in non extruded cells?
  #1
Senior Member
 
Blanco
Join Date: Mar 2009
Location: Torino, Italy
Posts: 193
Rep Power: 17
Blanco is on a distinguished road
Hello everybody,

I have a strange problem with a mesh I created through blockMesh.
The mesh is for a 2D case, you can find the blockMeshDict here attached.

After creating the mesh if I check its quality, using checkMesh utility, I get:

Code:
Checking geometry...
    Overall domain bounding box (-0.2 -0.9 0) (0.8 1.1 0.1)
    Mesh (non-empty, non-wedge) directions (1 1 0)
    Mesh (non-empty) directions (1 1 0)
    All edges aligned with or perpendicular to non-empty directions.
    Boundary openness (-4.47678e-17 -3.69151e-18 -1.16236e-15) OK.
    Max cell openness = 1.76034e-14 OK.
    Max aspect ratio = 169.31 OK.
    Minimum face area = 3.27985e-09. Maximum face area = 0.00622355.  Face area magnitudes OK.
    Min volume = 3.27985e-10. Max volume = 0.000114482.  Total volume = 0.193475.  Cell volumes OK.
    Mesh non-orthogonality Max: 62.6561 average: 27.9759
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 3.23037 OK.
    Coupled point location match (average 0) OK.
...and everything seems ok.
After this, I want to extrude the mesh starting from the top and bottom boundaries, therefore I issue an extrudeMesh command using the extrudeMeshDict here attached. Everything seems to go well but if then I issue again the checkMesh command I get the following problems:

Code:
Checking geometry...
    Overall domain bounding box (-0.2 -10.9 0) (0.8 11.1 0.1)
    Mesh (non-empty, non-wedge) directions (1 1 0)
    Mesh (non-empty) directions (1 1 0)
    All edges aligned with or perpendicular to non-empty directions.
    Boundary openness (2.07326944578e-16 5.89928008973e-19 -4.02186033899e-15) OK.
 ***High aspect ratio cells found, Max aspect ratio: 13192.513369, number of cells 5577
  <<Writing 5577 cells with high aspect ratio to set highAspectRatioCells
    Minimum face area = 3.27984629532e-09. Maximum face area = 0.02467.  Face area magnitudes OK.
    Min volume = 3.27984629533e-10. Max volume = 0.00087830134.  Total volume = 2.19347456322.  Cell volumes OK.
    Mesh non-orthogonality Max: 70.4886475638 average: 22.9159292824
   *Number of severely non-orthogonal (> 70 degrees) faces: 2.
    Non-orthogonality check OK.
  <<Writing 2 non-orthogonal faces to set nonOrthoFaces
    Face pyramids OK.
    Max skewness = 3.23049530527 OK.
    Coupled point location match (average 0) OK.

Failed 1 mesh checks.

End
If I check where are the 2 severe non-orthogonal cells...I find that they are placed inside the non-extruded mesh! No problematic cell is placed in the newly extruded zone (apart high aspect ratio cells which were expected).

Does anyone have an idea why this is happening? Or does anyone have already seen this behavior?

Regards,

Blanco
Attached Files
File Type: txt blockMeshDict.txt (16.1 KB, 6 views)
File Type: txt extrudeMeshDict.txt (3.0 KB, 26 views)

Last edited by Blanco; January 13, 2015 at 16:49.
Blanco is offline   Reply With Quote

Old   January 23, 2015, 11:54
Default Update -- SOLVED
  #2
Senior Member
 
Blanco
Join Date: Mar 2009
Location: Torino, Italy
Posts: 193
Rep Power: 17
Blanco is on a distinguished road
I update the status for anyone who could have the same problem:

It works if you write in binary or increase your writePrecision.
In particular, the problem still appear even if you use ascii with 12 as write precision...therefore binary is strongly suggested.
Blanco is offline   Reply With Quote

Old   January 23, 2015, 15:44
Default
  #3
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,974
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Greetings Blanco,

Many thanks for providing the solution to your issue!

For future reference, this was also reported as a bug here: http://www.openfoam.org/mantisbt/view.php?id=1497

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   February 25, 2023, 09:19
Default
  #4
New Member
 
Taihang Zhu
Join Date: May 2018
Location: London
Posts: 5
Rep Power: 7
zhutaihang is on a distinguished road
For someone having the same problem, if the suggestions from the above doesn't work, please try changing the extrudeModel into linearNormal.

I have been surfering from the same problem until I changed linearDirection into linearNormal.

Good luck!
zhutaihang 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
cellZone not taking all the cells inside rahulksoni OpenFOAM Running, Solving & CFD 6 January 25, 2019 01:11
cellZone not taking all the cells inside rahulksoni OpenFOAM 0 January 16, 2019 02:16
[snappyHexMesh] SHM Layer Addition Phase dickcruz OpenFOAM Meshing & Mesh Conversion 4 November 1, 2018 08:05
[snappyHexMesh] sHM layer process keeps getting killed MBttR OpenFOAM Meshing & Mesh Conversion 4 August 15, 2016 04:21
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 12:55


All times are GMT -4. The time now is 21:48.