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

Easy way to take surface average at different planes?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 29, 2011, 15:54
Question Easy way to take surface average at different planes?
  #1
Member
 
Charlie
Join Date: Dec 2010
Location: 415 Kinross Dr. Newark, DE 19711
Posts: 72
Rep Power: 5
cheng1988sjtu is on a distinguished road
Hi Foamers:

Now I'm working on channelFoam to reproduce the log-law and turbulence intensity. I just finished my first try using dynSmagorinsky. However, I've got some problems about the post-Processing.

1.
in the log file, I found that
************************************************** ************************************************** *
Time = 10

............
Uncorrected Ubar = 0.1335 pressure gradient = 0.000133009
ExecutionTime = 172139 s ClockTime = 345743 s

Calculating averages
************************************************** ************************************************** **

Is the 'pressure gradient' normalized by 'rho' already? I've seem some threads about it and got confused. Basically I need to use it to calculate the wall shear velocity, so it's very important to understand that.

2. Although in the 'controlDict' file, we have fieldAverage, but It seems that it is just a time average in the UMean, Pmean files, so I have 2 questions on that: a.> iis it averaged from t=0? even if it's not steady state? b>. Is there any existing utility that I can use to take a surface average of U and P? BTW, I may need to take surface average at several planes, so maybe sampling is not such a good choice for me. Any suggestions?

Thank you, foamers!

Zhen
cheng1988sjtu is offline   Reply With Quote

Old   January 31, 2011, 14:52
Default
  #2
Senior Member
 
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 5
stevenvanharen is on a distinguished road
Quote:
Originally Posted by cheng1988sjtu View Post
Hi Foamers:

Now I'm working on channelFoam to reproduce the log-law and turbulence intensity. I just finished my first try using dynSmagorinsky. However, I've got some problems about the post-Processing.

1.
in the log file, I found that
************************************************** ************************************************** *
Time = 10

............
Uncorrected Ubar = 0.1335 pressure gradient = 0.000133009
ExecutionTime = 172139 s ClockTime = 345743 s

Calculating averages
************************************************** ************************************************** **

Is the 'pressure gradient' normalized by 'rho' already? I've seem some threads about it and got confused. Basically I need to use it to calculate the wall shear velocity, so it's very important to understand that.
Yes, I think all the incompressible solvers assume rho = 1. That is why the density does not appear in the momentum equation and you have to specify kinematic viscosity instead of dynamic viscosity.

Quote:
Originally Posted by cheng1988sjtu View Post

2. Although in the 'controlDict' file, we have fieldAverage, but It seems that it is just a time average in the UMean, Pmean files, so I have 2 questions on that: a.> iis it averaged from t=0? even if it's not steady state? b>. Is there any existing utility that I can use to take a surface average of U and P? BTW, I may need to take surface average at several planes, so maybe sampling is not such a good choice for me. Any suggestions?

Thank you, foamers!

Zhen
a) Yes, do not expect the code to identify the steady state for you
b) If you play around with the 'sample' utility you can do quite a lot. See your source code for an example of the sampleDict file.
stevenvanharen is offline   Reply With Quote

Old   January 31, 2011, 23:35
Default thank you!
  #3
Member
 
Charlie
Join Date: Dec 2010
Location: 415 Kinross Dr. Newark, DE 19711
Posts: 72
Rep Power: 5
cheng1988sjtu is on a distinguished road
Hi Steven,

Thank you very much, it's really helpful.

I've looked into the source file for the field average, and yes, the average is taken from t=0 if you didn't specify the starting point of averaging.

Also the postChannel utility turns out to be a good choice for analyzing the channel flow.

Thank you!

Zhen
cheng1988sjtu is offline   Reply With Quote

Old   October 19, 2011, 17:46
Unhappy postChannel
  #4
New Member
 
Join Date: Mar 2011
Posts: 12
Rep Power: 5
tonkiplis is on a distinguished road
Hi Charlie,

I'm trying to run the channel395 tutorial. I encounter troubles at the post processing step. I added k in the fieldAverage function present in the controlDict.

My questions is: I obtain for k after postChannel utility:

k

kMean

k.xy
-----------------------

I assume k is the SGS turbulent kinetic energy (based on previous posts concerning post processing of LES) kMean is therefore <k> "time averaged value". k.xy would therefore be the <k> value averaged over the homogenous directions?

-----------------------

What is exaclty UPrime2Mean and u.xy, v.xy, w.xy??

I assume UPrime2Mean=<u_i' u_j'> and u.xy to be u_rms--> Sqrt(<u'u'>)?

Am I right?
-----------------------

Any inside would be apreciable..

How can I compare my k to k obtained by DNS?? I don't know how to account for the SGS turbulent kinetic energy.


Question solved.

Last edited by tonkiplis; October 25, 2011 at 02:08.
tonkiplis is offline   Reply With Quote

Reply

Tags
pressure gradient, surface average

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
[ICEM] How to generate sunstructured "all-tri patch-dependant" surface mesh in ICEM? jash ANSYS Meshing & Geometry 19 July 23, 2013 18:48
snappyHexMesh aborting Tobi OpenFOAM Native Meshers: snappyHexMesh and Others 0 November 10, 2010 03:23
surface planes in gambit falak FLUENT 9 October 14, 2005 23:59
Equilibrium of a free surface under surface tensio Ryan Main CFD Forum 1 August 7, 2001 16:14
free convection heat transfer from a heated horizontal surface through a liquid to a thin cooled fin Kaushik FLUENT 1 May 8, 2000 06:47


All times are GMT -4. The time now is 01:00.