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

fft of unsteady solution

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 14, 2004, 04:54
Default fft of unsteady solution
  #1
K S Chang
Guest
 
Posts: n/a
I have questions to you about FFT of unsteady solution.

I had acquired unsteady solution from 3D DES simulation.

I want to search the Stouhal number at one point from FFT. So I obtained the velocity history about 5000 time step(approximately 200 period). However the frequency is not unique and have so many peaks and the peak frequencies are different according to the range selected.

Q1) In general, how many periods are needed for FFT ?

Q2) If I have 4096(2^12) data, these 4096 data is totally done from a FFT ? Or 4 data set of 1024 data are done from FFT and averaged ?

Q3) If the peak frquency is different according to the range(approximately 50 periods) selected, the solution is not converged ??

Thank you in advance.

  Reply With Quote

Old   January 14, 2004, 05:38
Default Re: fft of unsteady solution
  #2
Steve
Guest
 
Posts: n/a
I would suggest doing something more rigorous like:

1) Work out your minimum required frequency. Use it to choose the sample length. From experience, make your sample long enough to contain at least two cycles of the minimum frequency required. Assuming that 1024 points is long enough...

2) Divide your time history into 8 sets of 1024 points so that there is 50% overlap between sets.

3) Multiply each set by a Hanning window (it's unlikely that your data is cyclic within the window).

4) Take the mean of the autopowers of the time histories (autopower = Fx . conj(Fx) where Fx is the complex spectrum). Square root it to get your spectrum.

5) Scale your data to account for the window: multply the spectrum by 2 to get the amplitudes right or by sqrt(8/3) to get the energy right (if you are summing across bands).

HTH

- Steve
  Reply With Quote

Old   January 14, 2004, 08:55
Default Re: fft of unsteady solution
  #3
Steve
Guest
 
Posts: n/a
Try also:

http://www.dspguide.com/specanal.htm

This is a section from a well and clearly written on-line DSP book.

- Steve
  Reply With Quote

Old   January 14, 2004, 20:48
Default Re: fft of unsteady solution
  #4
K S Chang
Guest
 
Posts: n/a
Thank you for your explanation.

I have one more question.

In DNS or LES simulation, above methods are adopted??

I can't search the general methods about fft of unsteady solution in CFD.

Any good reference ??

Thank you.
  Reply With Quote

Old   January 15, 2004, 05:22
Default Re: fft of unsteady solution
  #5
Steve
Guest
 
Posts: n/a
Sorry, can't really help any more. My knowledge of signal processing isn't really in the CFD field, although I'm sure the same basic principles apply to post-processing simulated data as do to post-processing measured data.
  Reply With Quote

Old   January 15, 2004, 05:43
Default Re: fft of unsteady solution
  #6
K S Chang
Guest
 
Posts: n/a
Thank you very much.

Your site is very helpful to me.

Thanks again.

^^

  Reply With Quote

Old   January 15, 2004, 07:42
Default Re: fft of unsteady solution
  #7
Steve
Guest
 
Posts: n/a
It's not my site. I just share a Christian name with the author!
  Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Convergence Centurion2011 FLUENT 48 June 15, 2022 00:29
conditions for unsteady problems dvdromnu OpenFOAM Running, Solving & CFD 5 November 9, 2009 03:15
How to get data from unsteady solution? asdf0901 FLUENT 1 June 3, 2009 11:34
Unsteady solution Christophe FLUENT 0 August 11, 2006 12:13
Unsteady simulation of flow past wheel Tom FLUENT 8 January 18, 2006 11:54


All times are GMT -4. The time now is 15:20.