CFD Online URL
[Sponsors]
Home > Forums > Main CFD Forum

Energy spectra

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

Like Tree2Likes
  • 2 Post By sbaffini

Reply
 
LinkBack Thread Tools Display Modes
Old   February 29, 2012, 20:27
Default Energy spectra
  #1
Member
 
Yong Wang
Join Date: Apr 2009
Posts: 32
Rep Power: 7
ywang is on a distinguished road
Hi, everyone

I have simulated turbulent channel flow with Re_tau=180. The mean velocity and stresses profiles agree with those in Kim, Moin and Moser's paper (JFM 1987, vol 177, pp 133-166). Now, I have problem to plot the one-dimensional spectra.

Taking the spectrum E_uu (k_x) for example. I calculated it with the following steps:
1. Output the 2D velocity field at y+ = 5, and then obtained uu(x,z).
2. With a given z, such as z=0, I performed FFT for uu(x,0) and obtained uu(k_x) at z=0.
3. Averaged uu(k_x) over all z.

The obtained spectrum was much different from the one in reference. However, I think the sharp of the profile is right.

I know there should be something wrong in my steps. And, the E_uu and k_x still need to be normalized by some reference variables, but I don't know how to choose them.

Any help from you would be appreciated. Thanks.


Yong
ywang is offline   Reply With Quote

Old   March 1, 2012, 07:43
Default
  #2
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 466
Blog Entries: 13
Rep Power: 14
sbaffini will become famous soon enough
Hi Yong,

the procedure you described is correct, for a given distance from the wall (say y*) you take your Nz (taken in the transverse coordinate direction) velocity vectors (say x velocity component along the x coordinate), FFT and average them (each coefficient is the average over the Nz ones).

It seems to me that the problem is in the proper "normalization".

I suggest you to check your spectral routine toward some test sin functions
sbaffini is offline   Reply With Quote

Old   March 1, 2012, 15:04
Default
  #3
Member
 
Yong Wang
Join Date: Apr 2009
Posts: 32
Rep Power: 7
ywang is on a distinguished road
Quote:
Originally Posted by sbaffini View Post
Hi Yong,

the procedure you described is correct, for a given distance from the wall (say y*) you take your Nz (taken in the transverse coordinate direction) velocity vectors (say x velocity component along the x coordinate), FFT and average them (each coefficient is the average over the Nz ones).

It seems to me that the problem is in the proper "normalization".

I suggest you to check your spectral routine toward some test sin functions
Hi Sbaffini,

Thanks for your reply. Yes, the key point is to find out the way to normalize k_x and E_uu.

I think there should be some ones who have simulated the same flow in this forum, as Kim's paper is seen as a benchmark for turbulent flow. I searched and read corresponding threads, but didn't find out the way.
ywang is offline   Reply With Quote

Old   March 1, 2012, 15:26
Default
  #4
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 466
Blog Entries: 13
Rep Power: 14
sbaffini will become famous soon enough
Ok, this morning i was too lazy to check my routine. Now i did it, you can find it attached. It works for MATLAB and, given the velocity field along a direction, it provides you with the correct abscissa and the coefficients to be averaged.

Just to be sure, this is the database we are talking about:

http://turbulence.ices.utexas.edu/da...n180.xspec.178
Attached Files
File Type: txt spectra.txt (416 Bytes, 95 views)
cfdonline2mohsen and mgg like this.
sbaffini is offline   Reply With Quote

Old   March 1, 2012, 15:39
Default
  #5
Member
 
Yong Wang
Join Date: Apr 2009
Posts: 32
Rep Power: 7
ywang is on a distinguished road
Quote:
Originally Posted by sbaffini View Post
Ok, this morning i was too lazy to check my routine. Now i did it, you can find it attached. It works for MATLAB and, given the velocity field along a direction, it provides you with the correct abscissa and the coefficients to be averaged.

Just to be sure, this is the database we are talking about:

http://turbulence.ices.utexas.edu/da...n180.xspec.178
Hi,

Thanks very much. Yes, I am using the data from the same website as the reference. I am comparing the spectra at y+=5.

I am using Qtiplot to perform FFT and don't have Matlab now. I will have a try. What are the characteristic parameters used for normalization?
ywang is offline   Reply With Quote

Old   March 1, 2012, 17:37
Default
  #6
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 466
Blog Entries: 13
Rep Power: 14
sbaffini will become famous soon enough
According to my memory there are no parameters used for normalization. However, the code you are using to perform the FFT may use some kind of windowing (like the Pwelch mentioned in the file). While this can be useful or even required for several reasons, this does not apply in this case because would invalidate the comparison (and you already average your coefficients in time and transverse directions, so it is not needed).
sbaffini is offline   Reply With Quote

Old   March 8, 2012, 21:24
Default
  #7
Member
 
Yong Wang
Join Date: Apr 2009
Posts: 32
Rep Power: 7
ywang is on a distinguished road
Quote:
Originally Posted by sbaffini View Post
According to my memory there are no parameters used for normalization. However, the code you are using to perform the FFT may use some kind of windowing (like the Pwelch mentioned in the file). While this can be useful or even required for several reasons, this does not apply in this case because would invalidate the comparison (and you already average your coefficients in time and transverse directions, so it is not needed).
Hi Sbaffini,

Thanks very much. Finally, I didn't compare my sepctra with those from Kim's paper. I just calculated the Komlogorov length scale eta, and plotted E_uu via k*eta.
ywang 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
energy spectra from LES simulations kumar OpenFOAM Post-Processing 0 October 29, 2009 11:42
Obtaining Energy Spectra from Physical variables CFDtoy Main CFD Forum 0 September 19, 2006 02:44
Energy Spectra! Murthy FLUENT 0 December 26, 2005 08:42
Energy Spectra! Murthy Main CFD Forum 0 December 26, 2005 08:41
turbulent energy spectra ulrich bieder Main CFD Forum 8 June 25, 1999 11:27


All times are GMT -4. The time now is 03:10.