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

Energy Spectrum

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 11, 2012, 17:11
Default Energy Spectrum
  #1
Member
 
Mohammad Saeedi
Join Date: Apr 2011
Location: Canada
Posts: 76
Rep Power: 6
saeedi is on a distinguished road
Hi All,
I wanted to make sure about the methodology to calculate the energy spectrum:
According to Pope's book, we should first calculate the autocovariance of the signal and then take the Fourier transform of that.

What I do is first calculate the autocovariance of u V.S. t signal( from a DNS) in MATLAB. when I want to take the FFT of that, there are seceral examples in its help with some subtle differences. Each gives a different answer (specially in the frequency range). Thay are all some how (not perfectly) parallel to -5/3 line. the last segment of the curve which should be vertically going down (pure dissipation) is also not visible.
I am sure about my signal (accuaracy and enough time of probing) but doubt about my methodology.

I'd appreciate if you give me your useful comment.

Thank you
saeedi is offline   Reply With Quote

Old   August 12, 2012, 04:00
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,664
Rep Power: 23
FMDenaro will become famous soon enough
Quote:
Originally Posted by saeedi View Post
Hi All,
I wanted to make sure about the methodology to calculate the energy spectrum:
According to Pope's book, we should first calculate the autocovariance of the signal and then take the Fourier transform of that.

What I do is first calculate the autocovariance of u V.S. t signal( from a DNS) in MATLAB. when I want to take the FFT of that, there are seceral examples in its help with some subtle differences. Each gives a different answer (specially in the frequency range). Thay are all some how (not perfectly) parallel to -5/3 line. the last segment of the curve which should be vertically going down (pure dissipation) is also not visible.
I am sure about my signal (accuaracy and enough time of probing) but doubt about my methodology.

I'd appreciate if you give me your useful comment.

Thank you

- what kind of turbulence are you simulating?
- how did you select the number of temporal-samples?
- are you considering the right range of frequencies until the Nyquist one?
- have you tried the same FFT for the spatial spectrum?
FMDenaro is offline   Reply With Quote

Old   August 12, 2012, 04:25
Default
  #3
Member
 
Mohammad Saeedi
Join Date: Apr 2011
Location: Canada
Posts: 76
Rep Power: 6
saeedi is on a distinguished road
-I am considering homogeneous turbulent channel flow
-I have a signal with 50,000 samples (each time step one sample and the total time of sampling is several times of the time required for flow to sweep the entire streamwise direction)

-I just know Nyquist theorem says sampling frequency should be properly higher than the signal frequency. What do you mean exactly of (considering the right range of frequencies until the Nyquist one) ?

-I did not get what you mean in your last question.
saeedi is offline   Reply With Quote

Old   August 12, 2012, 04:52
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,664
Rep Power: 23
FMDenaro will become famous soon enough
first, the Nyquist frequency is the highest one you solve in your simulation, for a well-done DNS must be of the order of the Kolmogorov one.

Try to do the spectrum of the stream-wise velocity along Kx and Kz at fixed half-height of the channel. You should repeat along the plane and take the average of each FFT. Remember that the k^-5/3 slope is theoretically expected for isotropic/homogeneous turbulence, in channel flow you could have such a range only far from the walls and at high Re_tau number (>200 - 300).
FMDenaro is offline   Reply With Quote

Old   August 12, 2012, 14:40
Default
  #5
Member
 
Mohammad Saeedi
Join Date: Apr 2011
Location: Canada
Posts: 76
Rep Power: 6
saeedi is on a distinguished road
My Re_tau is 300 and I considering half height.
You mean I need to take the fft along the spanwise direction and then take the average? what if we first do the plane averaging for different signals and then take the fft?

BTW, my main question is still remained: methodology for taking fft. There are several examples in Matlab help to take the fft of signals. which one should be used for my case?
should I take the fft of the signal itself or its autocovariance?
saeedi is offline   Reply With Quote

Old   August 12, 2012, 14:49
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,664
Rep Power: 23
FMDenaro will become famous soon enough
Quote:
Originally Posted by saeedi View Post
My Re_tau is 300 and I considering half height.
You mean I need to take the fft along the spanwise direction and then take the average? what if we first do the plane averaging for different signals and then take the fft?

BTW, my main question is still remained: methodology for taking fft. There are several examples in Matlab help to take the fft of signals. which one should be used for my case?
should I take the fft of the signal itself or its autocovariance?

if you take the FFT of a statistically averaged field that makes no sense!! Imagine the case of a RANS solution, you can not compute energy spectrum ....

I don't remember the examples in Matlab, some time ago I read of a sunspot sample. But you have only one way to do what you are looking for, you have u(x,y,z,t), v(x,y,z,t), w(x,y,z,t), so at fixed t=T and y = H/2:

- do the FFT of each velocity component along x for all the z position. The number of frequencies you have to use is fixed and extends up to the Nyquist one.
- do the average of the Fourier coefficient in the plane
- Compute the modulus of the averaged coefficients and plot along kx

- Repeat such procedure for computing the spectra along kz changing the role of x and z


Usually no temporal spectra are computed
FMDenaro is offline   Reply With Quote

Old   August 12, 2012, 16:36
Default
  #7
Member
 
Mohammad Saeedi
Join Date: Apr 2011
Location: Canada
Posts: 76
Rep Power: 6
saeedi is on a distinguished road
Do you mean spectrum at a fixed time (one snapshot)?

You say usually no temporal spectrum is computed? I am confused now. Because in turbulent flows of Pope, chapter 3 (equation 3-134 and round that) he explains how to calculate the spectrum from a temporal signal.
That's where he talks about calculation of autocovariance and then taking the Fourier transform?
saeedi is offline   Reply With Quote

Old   August 12, 2012, 16:42
Default
  #8
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,664
Rep Power: 23
FMDenaro will become famous soon enough
Quote:
Originally Posted by saeedi View Post
Do you mean spectrum at a fixed time (one snapshot)?

You say usually no temporal spectrum is computed? I am confused now. Because in turbulent flows of Pope, chapter 3 (equation 3-134 and round that) he explains how to calculate the spectrum from a temporal signal.
That's where he talks about calculation of autocovariance and then taking the Fourier transform?
Note that for channel flows a well-known database of DNS exists, it contains profiles, RMS, (spatial) spectra.
The spectra are computed at several times and then averaged ...
FMDenaro is offline   Reply With Quote

Old   April 21, 2014, 05:46
Default
  #9
Senior Member
 
Huang Xianbei
Join Date: Sep 2013
Location: CAU,China
Posts: 268
Rep Power: 5
huangxianbei is on a distinguished road
Hi,saeedi:
Have you figure out this problem now? I'm struggling in the problem
can you give me some advice?
huangxianbei 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
turbulent energy spectrum problem cfd.newbie Main CFD Forum 3 March 3, 2015 11:31
turbulent energy spectrum cfd.newbie CD-adapco 1 June 19, 2008 23:48
turbulent energy spectrum cfd.newbie FLUENT 0 June 18, 2008 18:34
LES correlation and turbulent energy spectrum Fabian Main CFD Forum 4 October 18, 2005 02:04
Energy Spectrum Emad Khalifa Main CFD Forum 3 June 30, 2003 16:03


All times are GMT -4. The time now is 13:56.