CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Severe nonorthogonality and severe skewness problem

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

Like Tree1Likes
  • 1 Post By andrewburns

Reply
 
LinkBack Thread Tools Display Modes
Old   January 22, 2008, 18:28
Default Dear all, I did high reynol
  #1
Member
 
Quinn Tian
Join Date: Mar 2009
Posts: 62
Rep Power: 8
qtian is on a distinguished road
Dear all,

I did high reynolds number RANS simulation case one month ago and the results seems ok. Now I am trying to do a low reynolds number RANS simualtion with LaunderSharmaKE turbulence model for comparison. I just refined the coarse mesh which I used in the high reynolds number simulation only along the vertical direction in order to resolve yplus=1. For some reason, my simulation for the low reynolds number case has convergency problem and with "unbound k and epsilon" and also has large continuity error. I also tried to turn off the turbulence first and no help also.

I generated mesh in gridgen with fluint format. Here is the results from fluentMeshToFoam and checkMesh. I don't know how to improve these severe non-orthogonality and skewed meshes? Does anyone know how to visualze these severe non-orthogonality and skewed meshes? I don't even know where they are. Thanks for your help.


Default patch type set to empty

Checking mesh...Number of non-orthogonality errors: 0. Number of severely non-orthogonal faces: 149.
--> FOAM Warning :
From function primitiveMesh::checkFaceSkewness(const bool report, labelHashSet* setPtr) const
in file meshes/primitiveMesh/primitiveMeshCheck.C at line 838
Large face skewness detected. Max skewness = 531.581 percent.
This may impair the quality of the result.
105 highly skew faces detected.
Failed 1 mesh geometry checks.
Failed some mesh checks.
done.
Writing mesh... to "constant/polyMesh" done.


checkMesh

Domain bounding box: min = (-0.469 -0.000233522 0) max = (0.78 0.254 0.468) meters.

Checking geometry...
Boundary openness in x-direction = -2.03288e-20
Boundary openness in y-direction = -4.14707e-18
Boundary openness in z-direction = 3.01544e-19
Boundary closed (OK).
Max cell openness = 2.26317e-19 Max aspect ratio = 45.748. All cells OK.

Minumum face area = 4.0582e-08. Maximum face area = 0.00063478. Face area magnitudes OK.

Min volume = 1.66625e-10. Max volume = 2.40859e-06. Total volume = 0.147875. Cell volumes OK.

Severe non-orthogonality for face 1 between cells 0 and 177: Angle = 76.8354 deg.
Severe non-orthogonality for face 526 between cells 175 and 352: Angle = 74.5599 deg.
Severe non-orthogonality for face 3242694 between cells 1086602 and 1086779: Angle = 77.6354 deg.
Severe non-orthogonality for face 3294987 between cells 1104125 and 1104302: Angle = 75.4307 deg.
Severe non-orthogonality for face 3347280 between cells 1121648 and 1121825: Angle = 72.0853 deg.
Number of non-orthogonality errors: 0. Number of severely non-orthogonal faces: 149.
Mesh non-orthogonality Max: 86.6762 average: 8.18947
Non-orthogonality check OK.

Writing 149 non-orthogonal faces to set nonOrthoFaces

Face pyramids OK.

Severe skewness for face 1058 skewness = 5.31581
Severe skewness for face 53351 skewness = 5.30193
Severe skewness for face 105644 skewness = 5.28783
Severe skewness for face 314816 skewness = 5.22916
Severe skewness for boundary face 3550338 skewness = 2.75832
Severe skewness for boundary face 3550339 skewness = 2.56192
Severe skewness for boundary face 3550340 skewness = 2.33748
Severe skewness for boundary face 3550341 skewness = 2.09081
--> FOAM Warning :
From function primitiveMesh::checkFaceSkewness(const bool report, labelHashSet* setPtr) const
in file meshes/primitiveMesh/primitiveMeshCheck.C at line 838
Large face skewness detected. Max skewness = 531.581 percent.
This may impair the quality of the result.
105 highly skew faces detected.
Writing 105 skew faces to set skewFaces

Minumum edge length = 1.54184e-05. Maximum edge length = 0.0312349.

All angles in faces are convex or less than 10 degrees concave.

Face flatness (1 = flat, 0 = butterfly) : average = 0.999998 min = 0.999867
All faces are flat in that the ratio between projected and actual area is > 0.8

Geometry check done.

Number of cells by type:
hexahedra: 1191564
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0
Number of regions: 1 (OK).
Failed 1 mesh checks.
qtian is offline   Reply With Quote

Old   January 22, 2008, 18:36
Default You need to adjust your discre
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
You need to adjust your discretisation on k and epsilon. Do upwind on convection and Gauss linear limited 0.7 on the laplacian. Apart from that you may need to tighten the pressure tolerance.

The above will remove the unboundedness problem and the code should behave OK.

If you want to visualise the bad cells, play around with
subsetMesh, which will allow you to see the bad cells.

Enjoy,

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   January 22, 2008, 19:47
Default I almost always have similar p
  #3
Member
 
Andrew Burns
Join Date: Mar 2009
Posts: 36
Rep Power: 8
andrewburns is on a distinguished road
I almost always have similar problems with k and epsilon, here's what I usually do.

Run a laminar solution for a few hundred iterations, or however long it takes for the residuals to drop to about 1e-3.

Turn turbulence back on, usually it still goes crazy at this point so then...

As Hrvoje Jasak said change the laplacian for k and epsilon to linear limited 0.7 and I also change the relaxation factors on k and e to either 0.5 or 0.3.

This usually works for me.
atg likes this.
andrewburns 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
Severe skewness Where dmoroian OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 7 August 21, 2014 06:46
I have a severe trouble with the compressible flow jinwon Main CFD Forum 1 June 22, 2007 04:57
Blockmesh problem about skewness aderliner OpenFOAM Native Meshers: blockMesh 0 December 1, 2006 22:04
severe error in cfx-post process chang CFX 0 May 18, 2006 11:43
SEVERE Operation 'calculate' failed with unrecogni Mohamed Musthafa CFX 0 September 8, 2005 10:38


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