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

[ICEM] Export mesh statistics for use in Richardson Extrapolation for grid sensitivity study

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

Like Tree1Likes
  • 1 Post By RodriguezFatz

Reply
 
LinkBack Thread Tools Display Modes
Old   April 16, 2013, 06:33
Question Export mesh statistics for use in Richardson Extrapolation for grid sensitivity study
  #1
Senior Member
 
RicochetJ's Avatar
 
Mr CFD
Join Date: Jun 2012
Location: Britain
Posts: 289
Rep Power: 6
RicochetJ is on a distinguished road
Hi,

Following the post from another user in this thread (Grid sensitivity analysis) I want to do a grid sensitivity study using the Richardson extrapolation method.

The first step is to define a representative cell, mesh or grid size, h. For three-dimensional simulations

h=\left[ \frac{1}{N}\sum^{N}_{i=1}{ \Delta V_{i}} \right]^{ \frac{1}{3}}

where \Delta V_{i} is the volume of the i^{th} cell, N is the total number of cells used for the computations.

My problem is I don't know how to extract the values of \Delta V_{i} from ICEM as ideally I'd like to do the Richardson Extrapolation method in Excel, where I can compare different meshes. In ICEM how do I export mesh statistics such as volume of each cell? Can it be exported as a handy .csv file?

I refer you to the paper I am using: http://journaltool.asme.org/Template...umAccuracy.pdf

Thank you
RicochetJ is offline   Reply With Quote

Old   April 16, 2013, 06:48
Default
  #2
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,097
Rep Power: 16
RodriguezFatz will become famous soon enough
Since you just need the ratio of h(fine)/h(coarse) calculation of an actual value of "h" is not necessary.
How do you want to do the refinement in ICEM?
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   April 16, 2013, 07:01
Default
  #3
Senior Member
 
RicochetJ's Avatar
 
Mr CFD
Join Date: Jun 2012
Location: Britain
Posts: 289
Rep Power: 6
RicochetJ is on a distinguished road
Quote:
Originally Posted by RodriguezFatz View Post
Since you just need the ratio of h(fine)/h(coarse) calculation of an actual value of "h" is not necessary.
How do you want to do the refinement in ICEM?
My variable of interest is wall heat transfer coefficient at the various boundaries. I am using the SST turbulence model. For accurate heat transfer predictions I need a Y+ < 1. So far the maximum Y+ I have on any of my boundaries is 1.6.

I will do the refinement in ICEM by applying more nodes in the thermal boundary layer.
RicochetJ is offline   Reply With Quote

Old   April 16, 2013, 07:03
Default
  #4
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,097
Rep Power: 16
RodriguezFatz will become famous soon enough
Do you use blocking / hexa meshing or unstructured?
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   April 16, 2013, 07:11
Default
  #5
Senior Member
 
RicochetJ's Avatar
 
Mr CFD
Join Date: Jun 2012
Location: Britain
Posts: 289
Rep Power: 6
RicochetJ is on a distinguished road
Quote:
Originally Posted by RodriguezFatz View Post
Do you use blocking / hexa meshing or unstructured?
I use blocking (hexa mesh). Before I export the mesh in .cfx5 format I right click on "Pre-Mesh" and click on "Convert to Unstructured Mesh".
RicochetJ is offline   Reply With Quote

Old   April 16, 2013, 07:13
Default
  #6
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,097
Rep Power: 16
RodriguezFatz will become famous soon enough
Great. Then you have all the numbers...
Do you refine the whole grid by using the mesh->refine utility, or do you change the number of nodes of some relevant edges?

Is this 2d? Can you post some picture?
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   April 16, 2013, 07:24
Default
  #7
Senior Member
 
RicochetJ's Avatar
 
Mr CFD
Join Date: Jun 2012
Location: Britain
Posts: 289
Rep Power: 6
RicochetJ is on a distinguished road
Quote:
Originally Posted by RodriguezFatz View Post
Great. Then you have all the numbers...
Do you refine the whole grid by using the mesh->refine utility, or do you change the number of nodes of some relevant edges?

Is this 2d? Can you post some picture?
It is 2D axisymmetric. I change the number of nodes of relevant edges. The geometry is shown below:



Do you want a picture of the mesh?
RicochetJ is offline   Reply With Quote

Old   April 16, 2013, 07:31
Default
  #8
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,097
Rep Power: 16
RodriguezFatz will become famous soon enough
How large is your coarsest mesh?
In 2d (and such a simple geometry) it is often affordable to just refine the whole mesh, which means you will end up with a 4 times larger mesh (2 x in "x" and 2 x in "y" direction). Would that be ok for you?
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   April 16, 2013, 07:48
Default
  #9
Senior Member
 
RicochetJ's Avatar
 
Mr CFD
Join Date: Jun 2012
Location: Britain
Posts: 289
Rep Power: 6
RicochetJ is on a distinguished road
Quote:
Originally Posted by RodriguezFatz View Post
How large is your coarsest mesh?
In 2d (and such a simple geometry) it is often affordable to just refine the whole mesh, which means you will end up with a 4 times larger mesh (2 x in "x" and 2 x in "y" direction). Would that be ok for you?
The coarsest mesh has 6500 elements and 3072 nodes. I do have the computational resources to double, triple and quadruple the number of elements and nodes.

However what I am interested in is the Richardson extrapolation method I mentioned in my original post.
RicochetJ is offline   Reply With Quote

Old   April 16, 2013, 08:05
Default
  #10
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,097
Rep Power: 16
RodriguezFatz will become famous soon enough
It's always the same method:
In your case, save your coarse mesh. Then refine the (unstructured) mesh, by clicking on "Edit Mesh -> Adjust Mesh Density". Here you can use "Refine All Mesh", method "pure refinement", steps "1". ->apply.
You will get a mesh, that is 2 times finer in each dimension.
Save that mesh. Do the same refinement again and you will end up with a even finer mesh.

Now run your simulation on all three meshes. Use the coarse results as interpolation for initializing the other runs.
For each simulation you need one single "important" value, such as you heat transfer coefficient or whatever. Using the nomenclature of your pdf, these values are the phi 3, 2 an 1. Since you doubled the number of points in each dimension the value of r21 and r32 is "2".
Here we go!
RicochetJ likes this.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   April 16, 2013, 08:45
Default
  #11
Senior Member
 
RicochetJ's Avatar
 
Mr CFD
Join Date: Jun 2012
Location: Britain
Posts: 289
Rep Power: 6
RicochetJ is on a distinguished road
Thank you! That actually does help.
RicochetJ is offline   Reply With Quote

Old   April 16, 2013, 09:34
Default
  #12
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,097
Rep Power: 16
RodriguezFatz will become famous soon enough
If you did that, you could post the values here for others to see how it works...
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   January 18, 2014, 19:57
Default
  #13
Member
 
don lenardo
Join Date: Nov 2010
Posts: 75
Rep Power: 6
lentschi is on a distinguished road
Hello,

I have one additional question regarding the Richardson Interpolation: What is about r21 and r32 if only the total number of cells is doubled (by splitting one direction)?

Example:

Coarse: 3 Million cells
Reference: 6 Million cells
Fine: 12 Million cells

I then r21=r32=1.25??

Thanks in advance.

//Markus
lentschi is offline   Reply With Quote

Reply

Tags
cfx, grid sensitivity, icem, richardson extrapolation

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
2D mesh export of an original 3D model Yogendra Shastri Main CFD Forum 0 February 22, 2008 14:14
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 12:55
Icemcfd 11: Loss of mesh from surface mesh option? Joe CFX 2 March 26, 2007 18:10
How to control Minximum mesh space? hung FLUENT 7 April 18, 2005 09:38
Gambit problems Althea FLUENT 21 February 6, 2001 08:05


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