# Turbulent Channel 1D Energy Spectra Calculation

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

 October 23, 2015, 17:13 Turbulent Channel 1D Energy Spectra Calculation #1 New Member   stran1028 Join Date: Oct 2015 Posts: 1 Rep Power: 0 Hello, I have 2 point correlation data for a turbulent channel case and I'm trying to put together a script to compute the 1D energy spectra from this. To verify my script, I'm using both 2 point correlations from DNS runs done by Moser et al (http://turbulence.ices.utexas.edu/data/MKM/chan395/) and trying to recreate their 1D energy spectra. Based on what I've seen in literature, the 1D energy spectra should just be the FFT of the 2 point correlation function. However when I try to do this in matlab using the fft function, I'm not able to match Moser's spectra. A snippet of my matlab code which calculates this is shown below. This was based on the examples shown for the fft function on the matlab website. yhat=fft(Ruu); tmpEuu= abs(yhat/L); Euu = tmpEuu(1:L/2 + 1); Euu(2:end-1) = Euu(2:end-1); Any help would be greatly appreciated! Thank you!

October 24, 2015, 03:38
#2
Senior Member

Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,290
Rep Power: 67
Quote:
 Originally Posted by stran1028 Hello, I have 2 point correlation data for a turbulent channel case and I'm trying to put together a script to compute the 1D energy spectra from this. To verify my script, I'm using both 2 point correlations from DNS runs done by Moser et al (http://turbulence.ices.utexas.edu/data/MKM/chan395/) and trying to recreate their 1D energy spectra. Based on what I've seen in literature, the 1D energy spectra should just be the FFT of the 2 point correlation function. However when I try to do this in matlab using the fft function, I'm not able to match Moser's spectra. A snippet of my matlab code which calculates this is shown below. This was based on the examples shown for the fft function on the matlab website. yhat=fft(Ruu); tmpEuu= abs(yhat/L); Euu = tmpEuu(1:L/2 + 1); Euu(2:end-1) = Euu(2:end-1); Any help would be greatly appreciated! Thank you!

you should search for similar post wherein this issue was already discussed

 Tags correlations, energy spectrum, fft, matlab, turbulent channel flow

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

 Similar Threads Thread Thread Starter Forum Replies Last Post kumar OpenFOAM Post-Processing 2 February 23, 2016 12:25 D_L Main CFD Forum 5 July 17, 2010 02:35 Joseph CFX 14 April 20, 2010 15:45 ulrich bieder Main CFD Forum 8 June 25, 1999 10:27 Zhong Lei Main CFD Forum 23 May 14, 1999 13:22

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