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

power spectral density

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 31, 2010, 19:39
Default power spectral density
  #1
Member
 
Naveen
Join Date: Feb 2010
Location: Los Angeles
Posts: 65
Rep Power: 8
vetnav is on a distinguished road
Hello,

I am carrying out a two dimensional simulation in a geometry with Nx X Ny number of grid points in stream wise and cross stream directions. At every time step I write the value of stream wise velocity component at a particular grid point into a datafile, I run for about 2 million time steps with a step size of 5e-04.

At the end of the run, I have a time signal of this stream wise velocity component at that particular grid point. The figure in the attached file shows the variation of u component with time at that particular location, now I want to see this signal on power - frequency plane, after lot of googling I found that if the signal is periodic one should go for power spectrum whereas if the signal is aperiodic one should go for power spectral density (is it true?). So I used the function pwelch in MATLAB to get the PSD which is also shown in attached file, but as per the definition PSD just gives the power between two frequencies but what I want to see is peaks at dominating frequencies and very small values at others, can someone please tell me how can I get it.

Thank you
Vetnav
Attached Files
File Type: doc PulsU.doc (75.0 KB, 46 views)
vetnav is offline   Reply With Quote

Old   September 1, 2010, 20:29
Default
  #2
Member
 
Naveen
Join Date: Feb 2010
Location: Los Angeles
Posts: 65
Rep Power: 8
vetnav is on a distinguished road
Hello

Is there any one who can help me with this?

thanks
vetnav is offline   Reply With Quote

Old   September 2, 2010, 11:22
Default Try...
  #3
New Member
 
Tomer
Join Date: Sep 2010
Posts: 9
Rep Power: 7
Avr.Tomer is on a distinguished road
Try Indirectly to FFT the signal and than extract PSD. Matlab usualy has a straiforward Help Guide but still one should be sure of the math behind the operator used.
Avr.Tomer is offline   Reply With Quote

Old   September 2, 2010, 11:30
Default Thanks
  #4
Member
 
Naveen
Join Date: Feb 2010
Location: Los Angeles
Posts: 65
Rep Power: 8
vetnav is on a distinguished road
Thank u Tomer for the reply, but I am not looking for the PSD, I want to see peaks at some particular frequencies (like a peak at 10Hz on power vs. frequency plot for a 10 Hz sine wave), PSD does not give this, it just gives total amount of power contained between two different frequencies.

Thanks
vetnav is offline   Reply With Quote

Old   September 2, 2010, 12:02
Default
  #5
New Member
 
Tomer
Join Date: Sep 2010
Posts: 9
Rep Power: 7
Avr.Tomer is on a distinguished road
I dont realy understand... If lets say you signal is composed only from a superposition of sine waves you would idealy get a sharp discrete frequency peaks via PSD in the proper fr., on the other hand white noise would give you a band with average frequency content.
now you should ask if the data you aquired is ergodicm otherwise frequancy analysis is somewhat meaningless
Avr.Tomer is offline   Reply With Quote

Old   September 2, 2010, 12:44
Default
  #6
Member
 
Naveen
Join Date: Feb 2010
Location: Los Angeles
Posts: 65
Rep Power: 8
vetnav is on a distinguished road
Thanks for the reply, You seem to be expert in this kind of stuff, but as per me I am doing this for the first time that is why I am using MATLAB. But reading your suggestions I am really lost now .

Let me think more and ask.

Thanks
Naveen
vetnav is offline   Reply With Quote

Old   September 3, 2010, 01:47
Default
  #7
New Member
 
Tomer
Join Date: Sep 2010
Posts: 9
Rep Power: 7
Avr.Tomer is on a distinguished road
Another suggestion which might achieve the goal of watching the frequency "at a particular time" is using wavelet transform methods.
Under some restrictions (ergodicity, sampling rate, etc...) by tuning its parameters it will allow you to get somewhat of "time - frequency" contour map which will enable you to view to some extent the energy content of a specific frequency and time (you might look at it as E=E(f,t)) .
Of course the Idea of E=E(f,t) is somewhat an otopic idea for when applying the transform (which might be considered as somewhat of a moving gaussian windowing) tou get both frequency content and time content of the energy "smudged".
Avr.Tomer is offline   Reply With Quote

Old   September 3, 2010, 08:08
Default Time Serie Analysis
  #8
New Member
 
Nicholas F Camus
Join Date: Sep 2009
Location: London, England.
Posts: 21
Rep Power: 8
Killercam is on a distinguished road
Dear Vetnav,

Tomer has it right. There is a wealth of analysis you can do with such a time-series. The first thing you should try (Tomer has suggested this), is to take the Fourier Transform of your time series; this will provide you with real and imaginary parts in frequency-period space. Take the sum of these parts squared and you have the power spectrum, this will tell you the quasi-periods present in your time-series (the real part should do this alone).

To analyse if there are any dominant noise processes at work you could use a structure function. Structure function analysis provides a method of quantifying time variability without the problem of aliasing, or windowing, that are encountered using the traditional FFT technique. Potentially it is able to provide information on the nature of the process that causes variability. The method is mainly concerned with the categorization of underlying noise processes and the identification of correlation time-scales.

Finally the mother of all analysis tools. Wavelets. Check out Torrence and Compo's website for all the information you need on these. Wavelets essentailly not only tell you if periodicity exists in your time-series but when it exists (temporal nature of the fluctuations).

Please contact me if you require further information and I will send you a PDF of a papar of mine...

Cheers,

Cam
Killercam is offline   Reply With Quote

Old   September 3, 2010, 14:28
Default
  #9
Member
 
Naveen
Join Date: Feb 2010
Location: Los Angeles
Posts: 65
Rep Power: 8
vetnav is on a distinguished road
Dear Tomer and Killercam

Thank you very much for the help, I calculated it again and I have more questions now looking at the result, I will go through the website that killercam suggested, but I want to ask you guys some more questions I prepared a document which contains different figures and its size is bigger than the upload limit here so can you please provide me your email ids? If you don't want to write your email ids here, u can email me at
vnaveen at ucla dot edu

Thank you
Naveen
vetnav 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
power spectrum density hosseini CFX 0 February 21, 2008 05:04
power spectral density cfd_newbie FLUENT 1 December 10, 2007 01:53
Wind Power Frank Grassi FLUENT 0 October 12, 2007 11:20
REAL GAS UDF brian FLUENT 6 September 11, 2006 08:23
Warning 097- AB CD-adapco 6 November 15, 2004 05:41


All times are GMT -4. The time now is 12:52.