extracting kinetic energy

argse2 · November 12, 2009, 03:32

hi foamers,

I am new to openfoam. i created a mesh and modified pisoFoam to my case.
now I need to know the kinetic energy m/2 * v of each control volume.
Is there a post processing tool, a function in paraFoam which can do that task?
Otherwise I like to know how to read the solution domain mesh and iterator through all cells.

best regards and thx

kumar · November 13, 2009, 04:44

Hello,
I think there is no direct utility to extract kinetic energy, but it is possible by simple coding and compiling.
It is also a good excercise for understanding the way to make a new application and compile it.

I looked in to this post and did exactly as it is specified in this post and my utility to extract KE was ready. So i want you to check this link

http://www.cfd-online.com/Forums/ope...ic-energy.html

hope it helps.
bye
K.Suresh kumar

argse2 · November 13, 2009, 10:18

Hi Kumar,

thx alot. coding seems to be quite easy. but according to linked post, shouldn't it be
0.5 * m * U & U or am i completely wrong?

best regards

muellea · November 13, 2009, 13:24

Well, it depend on what you wan't to know. If you wanna know the kinetic energy of your entire computational domain, then you are right (you have to take the total mass then).
But in most cases you are interested in a kinetic energy field, so it is convenient to represented:

either mass-specific: 0.5*U*U
or volume-specifc: 0.5*rho*U*U

The other representation with the mass is not used in continuums mechanics, since the mass of a particle (cell) is grid-dependent.

do you agree?

argse2 · November 13, 2009, 14:25

hi muellea,

i agree with you. thx.

maybe my previous post have been to vague.

i need to compare the solution of our inhouse fv solver with those of openfoam. our solver computes the volume of each control volume, multiplies it with the fluid's density and the mass specific kinetic energy as i (almost) wrote in my initial question.

so 0.5 * rho * mesh.V() & U & U should do the job if openfoam computes on the same grid than ours?

muellea · November 13, 2009, 17:23

Hi Argse2,

sorry i should have read your first post properly. I think it should work like that, except that i would take 0.5*rho*mesh.V()*U&U, since the & stands for the inner product and Mesh.V() and rho are as far as i know scalar-fields, whereas U is a vectorfield.

Hopefully you get it run. And let's see whether Openfoam beats your inhouse solver

.

November 12, 2009, 03:32	extracting kinetic energy	#1
argse2 New Member Argse Join Date: Nov 2009 Posts: 3 Rep Power: 17	hi foamers, I am new to openfoam. i created a mesh and modified pisoFoam to my case. now I need to know the kinetic energy m/2 * v of each control volume. Is there a post processing tool, a function in paraFoam which can do that task? Otherwise I like to know how to read the solution domain mesh and iterator through all cells. best regards and thx

November 13, 2009, 13:24		#4
muellea New Member anonymous Join Date: Nov 2009 Location: Switzerland Posts: 8 Rep Power: 17	Well, it depend on what you wan't to know. If you wanna know the kinetic energy of your entire computational domain, then you are right (you have to take the total mass then). But in most cases you are interested in a kinetic energy field, so it is convenient to represented: either mass-specific: 0.5UU or volume-specifc: 0.5rhoU*U The other representation with the mass is not used in continuums mechanics, since the mass of a particle (cell) is grid-dependent. do you agree?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
LES and total turbulent kinetic energy	Boerge	FLUENT	1	September 8, 2012 12:41
How to view subgrid kinetic energy	Craig	FLUENT	3	September 8, 2012 12:34
ATTENTION! Reliability problems in CFX 5.7	Joseph	CFX	14	April 20, 2010 16:45
Mean flow kinetic energy	Jenny	CFX	3	November 19, 2008 21:12
what's the macro of SGS turbulent kinetic energy?	lcw	FLUENT	1	June 13, 2006 04:08

November 13, 2009, 04:44		#2
kumar Senior Member Suresh kumar Kannan Join Date: Mar 2009 Location: Luxembourg, Luxembourg, Luxembourg Posts: 129 Rep Power: 17	Hello, I think there is no direct utility to extract kinetic energy, but it is possible by simple coding and compiling. It is also a good excercise for understanding the way to make a new application and compile it. I looked in to this post and did exactly as it is specified in this post and my utility to extract KE was ready. So i want you to check this link http://www.cfd-online.com/Forums/ope...ic-energy.html hope it helps. bye K.Suresh kumar

November 13, 2009, 10:18		#3
argse2 New Member Argse Join Date: Nov 2009 Posts: 3 Rep Power: 17	Hi Kumar, thx alot. coding seems to be quite easy. but according to linked post, shouldn't it be 0.5 * m * U & U or am i completely wrong? best regards

November 13, 2009, 14:25		#5
argse2 New Member Argse Join Date: Nov 2009 Posts: 3 Rep Power: 17	hi muellea, i agree with you. thx. maybe my previous post have been to vague. i need to compare the solution of our inhouse fv solver with those of openfoam. our solver computes the volume of each control volume, multiplies it with the fluid's density and the mass specific kinetic energy as i (almost) wrote in my initial question. so 0.5 * rho * mesh.V() & U & U should do the job if openfoam computes on the same grid than ours?

November 13, 2009, 17:23		#6
muellea New Member anonymous Join Date: Nov 2009 Location: Switzerland Posts: 8 Rep Power: 17	Hi Argse2, sorry i should have read your first post properly. I think it should work like that, except that i would take 0.5rhomesh.V()*U&U, since the & stands for the inner product and Mesh.V() and rho are as far as i know scalar-fields, whereas U is a vectorfield. Hopefully you get it run. And let's see whether Openfoam beats your inhouse solver .