CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > ANSYS Meshing & Geometry

[ICEM] courant number <-> mesh quality

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 20, 2012, 05:20
Default courant number <-> mesh quality
  #1
Member
 
Join Date: Jul 2012
Posts: 31
Rep Power: 5
kpax is on a distinguished road
hey folks,

i am currently encountering the problem that even for very small time-steps (0.0001), the courant number of my problem explodes at some point of the simulation (in icoFoam). Since i have to simulate ~50s in total, I really dont wanna go to even smaller time steps.
I think the problem is mesh quality, because the simulation stayed stable a bit longer when I used more non-orthogonal corrector steps. plus i could see in the results that the velocity exploded in regions of low mesh quality...

so my question is this: what kind of mesh errors have the biggest influence on the courant number, according to your experience?
so far, i have mainly focused on improving equiangle skewness, but i'm not sure if that is really the best choice.


btw: both "smooth elements globally" and "smooth hexahedral mesh - orthogonal" did not work for me at all (made it all worse) - is there someone else who experienced that?


best,
kpax
kpax is offline   Reply With Quote

Old   November 20, 2012, 05:43
Default
  #2
Senior Member
 
AB
Join Date: Sep 2009
Location: France
Posts: 323
Rep Power: 12
BrolY will become famous soon enough
First of all, you have to know that if you increase the determinant too much, it could decrease the angle, and vice versa .. But, if you increase the aspect ratio, it could increase the determinant. Most of the criteria are linked.

From my experience, angle is worth for hexa. A very bad angle could crash your simulation. Could you post some pictures of your problem. For example, a picture of the region where your simulation crashes with the picture of your mesh at the same location.

Also, it depends on the conditions of your simulation : very high velocity in area where your mesh is not refinned is not good at all !
BrolY is offline   Reply With Quote

Old   November 20, 2012, 09:03
Default
  #3
Member
 
Join Date: Jul 2012
Posts: 31
Rep Power: 5
kpax is on a distinguished road
ok, attached you find screenshots of the areas of high velocity (at the outlet), the mesh at the outlet, and areas of low skew values (0.4-0.5, i thought that would not be too bad?)..
obviously, the high u values appear at the corners of the o-grid. i tried moving the vertices around, but it didn't improve results. weird, isnt it?



Quote:
Originally Posted by BrolY View Post
Also, it depends on the conditions of your simulation : very high velocity in area where your mesh is not refinned is not good at all !
i thought a fine mesh and high velocity would be the worst that can happen in terms of the courant number? because c = u*dt/dx then becomes large...
Attached Images
File Type: jpg outlet.jpg (51.1 KB, 46 views)
File Type: png skew.png (67.8 KB, 38 views)
File Type: png u.png (43.9 KB, 42 views)
kpax is offline   Reply With Quote

Old   November 20, 2012, 10:46
Default
  #4
Senior Member
 
AB
Join Date: Sep 2009
Location: France
Posts: 323
Rep Power: 12
BrolY will become famous soon enough
About the high velocity and fine mesh, the problem doesn't come from the Courant, but comes from the solver. You need more elements to capture more information.

About the velocity profile at the outlet, it looks weird.
What solver do you use ? Are you sure everything is correctly configured with your solver ? Because, from the picture of your mesh, it looks like there is no trouble with the mesh !
Unless other users see something I do not, maybe you should ask your question in the solver forum.
BrolY is offline   Reply With Quote

Old   November 21, 2012, 02:07
Default
  #5
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,914
Blog Entries: 6
Rep Power: 38
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
did you try the elliptical smoother? It comes under blocking and use the option orthogonality.
Far is offline   Reply With Quote

Old   November 21, 2012, 10:02
Default
  #6
Member
 
Join Date: Jul 2012
Posts: 31
Rep Power: 5
kpax is on a distinguished road
Quote:
Originally Posted by BrolY View Post
About the velocity profile at the outlet, it looks weird.
What solver do you use ? Are you sure everything is correctly configured with your solver ? Because, from the picture of your mesh, it looks like there is no trouble with the mesh !
Unless other users see something I do not, maybe you should ask your question in the solver forum.
broly, you were absolutely right. some of my BC were not adjusted correctly... once if fixed that the velocity profile looked alright.
thx a lot for pointing me in the right direction!!


Quote:
Originally Posted by Far View Post
did you try the elliptical smoother? It comes under blocking and use the option orthogonality.
yep, made things worse. but maybe that is because my mesh is in fact not too bad to begin with? or i might not have chosen the best configuration, there are a lot of options...
kpax is offline   Reply With Quote

Old   November 26, 2012, 09:22
Default
  #7
Member
 
Join Date: Jul 2012
Posts: 31
Rep Power: 5
kpax is on a distinguished road
hey everybody,

turns out i still have some problems with the courant number. after around 50 time steps (each with dt = 0.0001), the simulation crashes due to very high U and p values. i attached a picture of the area where p values get very high - obviously not a reasonable result.
strangely though, this area does not seem to have worse-than-average mesh elements, as you can see in the other picture, which shows elements with determinant between 0.45 and 0.65 (all other elements have higher determinant).

in some publications, i have seen reasonable results for my kind of simulation with roughly the same number of mesh elements, but dt = 0.01. i wonder what's the difference?
Attached Images
File Type: jpg pressure.jpg (17.0 KB, 17 views)
File Type: jpg meshErrors.jpg (14.0 KB, 21 views)
kpax is offline   Reply With Quote

Old   December 20, 2012, 12:06
Default
  #8
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
Hi kpax,
I am using OpenFOAM too, and I am facing some similary problem with you.
same geometry,structured and unstructured mesh,different behaviour.
in unstructured mesh the quality is not good. in structured mesh the quality is above 0.8. I think thats the mesh quality which incur the result explode on twophaseeulerfoam. If you handle your problem, I can get it a hint.~
sharonyue 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
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20
How to Do External Mesh for Airfoil sHM msuaeronautics OpenFOAM Native Meshers: snappyHexMesh and Others 1 September 23, 2012 04:00
Add Mesh Layers doesnt work on the whole surface Kryo OpenFOAM Native Meshers: snappyHexMesh and Others 8 September 13, 2012 09:28
snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Native Meshers: snappyHexMesh and Others 2 March 27, 2011 21:11
BlockMeshmergePatchPairs hjasak OpenFOAM Native Meshers: blockMesh 11 August 15, 2008 07:36


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