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

[ICEM] Exporting Grid Volumes/Areas from ICEM to Calculate Discretization Error

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 18, 2010, 14:12
Default Exporting Grid Volumes/Areas from ICEM to Calculate Discretization Error
  #1
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 9
Josh is on a distinguished road
Hi -

I'm trying to calculate the discretization error of my CFD simulations using Richardson's method as outlined in the Journal of Fluids Engineering:
http://journaltool.asme.org/Template...umAccuracy.pdf

I created my "2D" (one element deep) structured mesh using ICEM 12. If you observe Equation 2 (fourth page of the above link), you'll see that a representative cell size is required. Basically, I need the following parameters to solve the equation:

N = total number of cells
Ai = area of each cell (the i'th cell)
or...
Vi = volume of each cell

I know how to get N, but I don't know how to export a file from ICEM with the area (or volume) of each cell.

Is it possible to export a file from ICEM with the volume of each cell? If so, how?

Thanks for any help!
Josh is offline   Reply With Quote

Old   May 19, 2010, 12:05
Default Hmm...
  #2
Senior Member
 
PSYMN's Avatar
 
Simon Pereira
Join Date: Mar 2009
Location: Ann Arbor, MI
Posts: 2,662
Blog Entries: 1
Rep Power: 36
PSYMN has a spectacular aura aboutPSYMN has a spectacular aura about
There are diagnostics that create a historgram for a property such as volume, but I don't know how to do exactly what you want. I will look into it and get back to you...
PSYMN is offline   Reply With Quote

Old   May 19, 2010, 19:48
Default
  #3
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 9
Josh is on a distinguished road
Thanks, as always, Simon.
Josh is offline   Reply With Quote

Old   May 20, 2010, 11:58
Default Scripting...
  #4
Senior Member
 
PSYMN's Avatar
 
Simon Pereira
Join Date: Mar 2009
Location: Ann Arbor, MI
Posts: 2,662
Blog Entries: 1
Rep Power: 36
PSYMN has a spectacular aura aboutPSYMN has a spectacular aura about
Not sure how comfortable you are with scripting, but I think it will be the only way to produce a file with this info for each element.

You will need to look into the programmers guide and tcl scripting in general, but I asked one of my friendly neighborhood developers (actually he is on a different continent) for some tips and he sent me this...

At least it looks easy for the surface elements.

Enjoy...

There exists 2 commands to get volume/area information.

set area [ic_uns_area $subset_name]
set vol [ic_uns_volume $subset_name]

What the user needs to do is to create a subset ($subset_name) for each element in the mesh and to write this information into a file.

The following existing ic_ command writes the area information (only for shell elements) into a file:

ic_uns_write_elem_list all ./my_file_name

Output format of the file is

"element number" "area" "element node numbers"

Hope that helps a little.
Jari
PSYMN is offline   Reply With Quote

Old   May 20, 2010, 13:40
Default
  #5
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 9
Josh is on a distinguished road
Thanks, Simon. Chances are I'll never get around to doing this, but it's good to have this as a reference should I have the time to try it.
Josh 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
UDF: DEFINE_CG_MOTION for vertical jump motion of an electrode! alban Fluent UDF and Scheme Programming 2 June 8, 2010 18:54
POSDAT problem piotka STAR-CD 4 June 12, 2009 08:43
Errors running allwmake in OpenFOAM141dev with WM_COMPILE_OPTION%3ddebug unoder OpenFOAM Installation 11 January 30, 2008 21:30
Compiling OpenFOAM13 on AMD64 with OpenSUSE 101 silent_missile OpenFOAM Installation 5 August 10, 2007 07:31
Combustion Convergence problems Art Stretton Phoenics 5 April 2, 2002 05:59


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