CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Pre-Processing

Problem : Cell expansion ratio (blockMeshDict)

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 16, 2012, 06:16
Default Problem : Cell expansion ratio (blockMeshDict)
  #1
New Member
 
Seon Hye
Join Date: Feb 2012
Location: South Korea
Posts: 18
Rep Power: 5
ysh1227 is on a distinguished road
Hi, guys.

I try to make my mesh for simulating 2D flow over the airfoil using blockMesh, so I compiled my blockMeshDict file. But when I check in paraView, I got a problem that the mesh is twisted.

First, I inputted small expansion ratio about 20 in normal direction to the airfoil surface, and there is no problem.

And then, I increased the expansion ratio (about 100). blockMesh was good working and no error existed. However, in paraView, The mesh had big trouble.

I attach a part of my blockMeshDict (related to definition of block) and picture about bad mesh.

How can I resolve it?


PHP Code:
blocks
(
    
hex (11 10 0 1 27 26 16 17) (             60            256              1simpleGrading (1 1000 1)
    
hex (12 11 1 2 28 27 17 18) (             32            256              1simpleGrading (1 1000 1)
    
hex (13 12 2 3 29 28 18 19) (             36            256              1simpleGrading (1 1000 1)
    
hex (5 13 3 4 21 29 19 20) (            256            256              1simpleGrading (0.001 1000 1)
    
hex (10 15 9 0 26 31 25 16) (             60            256              1simpleGrading (1 1000 1)
    
hex (15 14 8 9 31 30 24 25) (             32            256              1simpleGrading (1 1000 1)
    
hex (14 13 7 8 30 29 23 24) (             36            256              1simpleGrading (1 1000 1)
    
hex (13 5 6 7 29 21 22 23) (            256            256              1simpleGrading (1000 1000 1)
); 
Attached Images
File Type: jpg mesh1.jpg (98.0 KB, 28 views)
File Type: jpg mesh2.jpg (86.2 KB, 29 views)
ysh1227 is offline   Reply With Quote

Old   July 16, 2012, 08:53
Default
  #2
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
What is the output of checkMesh?
Bernhard is offline   Reply With Quote

Old   July 16, 2012, 09:28
Default
  #3
New Member
 
Seon Hye
Join Date: Feb 2012
Location: South Korea
Posts: 18
Rep Power: 5
ysh1227 is on a distinguished road
The output of checkMesh is

PHP Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  2.1.0-0bc225064152
Exec   
checkMesh
Date   
Jul 16 2012
Time   
22:26:18
Host   
"wtcfd"
PID    12055
Case   : /home/wtcfd/MSthesis/MATLABcodetest
nProcs 
1
sigFpe 
Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking Monitoring run-time modified files using timeStampMaster
allowSystemOperations 
Disallowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create polyMesh 
for time 0
Time 
0
Mesh stats
    points
:           394752
    internal points
:  0
    faces
:            787200
    internal faces
:   392448
    cells
:            196608
    boundary patches
4
    point zones
:      0
    face zones
:       0
    cell zones
:       0
Overall number of cells of each type
:
    
hexahedra:     196608
    prisms
:        0
    wedges
:        0
    pyramids
:      0
    tet wedges
:    0
    tetrahedra
:    0
    polyhedra
:     0
Checking topology
...
    
Boundary definition OK.
    
Cell to face addressing OK.
    
Point usage OK.
    
Upper triangular ordering OK.
    
Face vertices OK.
    
Number of regions(OK).
Checking patch topology for multiply connected surfaces ...
    
Patch               Faces    Points   Surface topology                  
    inlet               768      1538     ok 
(non-closed singly connected)  
    
outlet              512      1026     ok (non-closed singly connected)  
    
airfoil             256      512      ok (non-closed singly connected)  
    
frontAndBack        393216   394752   ok (non-closed singly connected)  
Checking geometry...
    
Overall domain bounding box (-30 -30 0) (30 30 0.2)
    
Mesh (non-empty, non-wedgedirections (1 1 0)
    
Mesh (non-empty) directions (1 1 0)
    
All edges aligned with or perpendicular to non-empty directions.
    
Boundary openness (-2.14328e-19 -9.24849e-20 2.84901e-16OK.
 ***
High aspect ratio cells foundMax aspect ratio1000number of cells 4
  
<<Writing 4 cells with high aspect ratio to set highAspectRatioCells
    Minumum face area 
9.46298e-08Maximum face area 0.644086.  Face area magnitudes OK.
    
Min volume 2e-300Max volume 0.128817.  Total volume 643.668.  Cell volumes OK.
    
Mesh non-orthogonality Max179.984 average20.1941
   
*Number of severely non-orthogonal faces914.
 
***Number of non-orthogonality errors2500.
  
<<Writing 3414 non-orthogonal faces to set nonOrthoFaces
 
***Error in face pyramids8016 faces are incorrectly oriented.
  <<
Writing 5730 faces with incorrect orientation to set wrongOrientedFaces
 
***Max skewness 23.752771 highly skew faces detected which may impair the quality of the results
  
<<Writing 71 skew faces to set skewFaces
    Coupled point location match 
(average 0OK.
Failed 4 mesh checks
ysh1227 is offline   Reply With Quote

Old   July 16, 2012, 10:17
Default
  #4
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
Only thing I get think of is the writePrecision in controlDict, maybe you want to set that to a larger value.
Bernhard is offline   Reply With Quote

Old   July 16, 2012, 10:28
Default
  #5
New Member
 
Seon Hye
Join Date: Feb 2012
Location: South Korea
Posts: 18
Rep Power: 5
ysh1227 is on a distinguished road
Thank you for your advice
Unfortunately, I can not correct my problem with your suggestion.
I increased writePrecision till 20, but results of the mesh that I can see in paraView are all the same.

Do you need more information to resolve this trouble?
ysh1227 is offline   Reply With Quote

Old   July 16, 2012, 10:33
Default
  #6
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
Than you should share everything one needs to reproduce the problem, so that we can try ourselves.
Bernhard is offline   Reply With Quote

Old   July 16, 2012, 10:47
Default
  #7
New Member
 
Seon Hye
Join Date: Feb 2012
Location: South Korea
Posts: 18
Rep Power: 5
ysh1227 is on a distinguished road
I can not understand this problem.
SAME blockMeshDict file, but Different expension ratio.
For the small expansion ratio, there is no matter.
For the bigger expansion ratio, there is bad results.

First, I'll try more things by myself.
If I do not find a solution, I'll upload blockMeshDict.

Once again, I heartily thank you for your advice
ysh1227 is offline   Reply With Quote

Old   July 17, 2012, 23:31
Default
  #8
New Member
 
Seon Hye
Join Date: Feb 2012
Location: South Korea
Posts: 18
Rep Power: 5
ysh1227 is on a distinguished road
Bernhard, I decided to make a blockMeshDict file again with different structure for blocks.

Last edited by ysh1227; July 18, 2012 at 03:35.
ysh1227 is offline   Reply With Quote

Old   July 22, 2012, 08:39
Default
  #9
New Member
 
Seon Hye
Join Date: Feb 2012
Location: South Korea
Posts: 18
Rep Power: 5
ysh1227 is on a distinguished road
.Bernhard, I decided to make a blockMeshDict file again with different structure for blocks.

Last edited by ysh1227; July 22, 2012 at 10:01.
ysh1227 is offline   Reply With Quote

Old   September 10, 2012, 03:56
Default
  #10
Senior Member
 
sivakumar selvaraju
Join Date: Mar 2009
Location: Cape Town - South Africa
Posts: 186
Rep Power: 8
sivakumar is on a distinguished road
Send a message via Skype™ to sivakumar
Can you please attach your whole blockMeshDict, so that we can have a look of that.
It can be problem with your blockMeshDict.

Siva
sivakumar is offline   Reply With Quote

Old   September 10, 2012, 04:12
Default
  #11
New Member
 
Seon Hye
Join Date: Feb 2012
Location: South Korea
Posts: 18
Rep Power: 5
ysh1227 is on a distinguished road
I resolve my problem.
Output of my blockMeshDict shows some errors about the order of definition of block. But the mesh has no errors in paraView.
So I thought that is not critical and continued my calculation.
When I use simpleFoam, the calculation is perfect.
Thank you for your concern.
ysh1227 is offline   Reply With Quote

Reply

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
equivalence ratio calculation on star-cd cell kri321shna STAR-CD 4 April 11, 2012 05:34
Problem on Star-CD deactive cell command willy511 STAR-CD 2 March 27, 2012 21:19
Cells with t below lower limit Purushothama CD-adapco 2 May 31, 2010 21:58
Problem on high density ratio in Level Set method Kai Yan Main CFD Forum 10 December 25, 2007 07:12
v4 trimmed cell problem work CD-adapco 3 December 18, 2007 18:07


All times are GMT -4. The time now is 15:18.