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

turbulent energy spectrum problem

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 1 Post By Patrick
  • 2 Post By romant

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 18, 2008, 19:32
Default turbulent energy spectrum problem
  #1
cfd.newbie
Guest
 
Posts: n/a
Hi

I want to calculate energy spectrum for my LES simulation. I put a point in domain and monitored u, v and w velocities. I calculated the kinetic energy from my velocities and imported the kinetic data and time into matlab. I, then, converted my time domain into frequency domain by using fft.

How do I draw energy spectrum then? In book there is E(k) on y-axis and k (wave number) on x-axis. Is E(k) same as kinetic energy which is 0.5(u^2+v^2+w^2)? How do I calculate k (wave number) on x-axis? Thanks in advance for your answer.

  Reply With Quote

Old   June 20, 2008, 13:00
Default Re: turbulent energy spectrum problem
  #2
Patrick
Guest
 
Posts: n/a
In order to calculate the energy spectrum you need to have the velocities monitores on all your grid points. You need to compute E as the sum of all the kinetic energies in all the 3 dimension (as you wrote) for all your domain points at only one time t! THen you can either average E in each dimension or just keep E as a function of the entire domain. Now you should have this done for many time t, and averaged over the entire period T=t1+t+2+..tn. ONce you have the time-averaged kinetic energy E as a function of space you can carry out FFT in each dimension (either at a specific location or again averaging over surfaces).

Say in a one-dimensional problem you have

E=rho*v_x^2 at a time t.

You compute E(t) for all t from t=0 to t=T and you average E over T

E = sum_i E(t_i)*(t_i+1 t_i)/T

THen you are left with E(x) at discrete location x_i.

Then you carry out an FFT on E(x_i) (i=1,N) and you obtain E(k).

If your entire domain is D and you have N points there, then delta x= D/N. k=1 correspond to a wavelenght lambda=D k=2 correspoinds to a wavelength D/2

etc...

k=N/2 corresponds to a wavelength 2D/N

In two dimension you will have different k's: k_x and k_y. To compute E(k_x) you need to average E(x,y) over t again but also over y so that you are left with E(x) only and again carry out an FFT. For the y dimension you need to average over x and get E(y). An FFT of E(y) brings you E(k_y).

In 3D you have k_x, k_y, k_z. In order to obtain E(x) you need to averages (discrete integration) over y and z, and of course t, then you can do the FFT for x. similarly for y and z.

The Energy spectrum is not as a function of frequency (time) but it is about wavelength with k=2*pi/lambda where lambda is the wavelength (spatial dimension, NOT time).

Look into Numerical Recepices. Check CFD books and see what they have on this CFD site.

It is not the time domain that you convert into frequencies but the space domain (after averaging over time).
cfdonline2mohsen likes this.
  Reply With Quote

Old   January 25, 2015, 10:36
Default
  #3
New Member
 
wang
Join Date: Dec 2014
Posts: 2
Rep Power: 0
wgq521 is on a distinguished road
Quote:
Originally Posted by Patrick
;55543
In order to calculate the energy spectrum you need to have the velocities monitores on all your grid points. You need to compute E as the sum of all the kinetic energies in all the 3 dimension (as you wrote) for all your domain points at only one time t! THen you can either average E in each dimension or just keep E as a function of the entire domain. Now you should have this done for many time t, and averaged over the entire period T=t1+t+2+..tn. ONce you have the time-averaged kinetic energy E as a function of space you can carry out FFT in each dimension (either at a specific location or again averaging over surfaces).

Say in a one-dimensional problem you have

E=rho*v_x^2 at a time t.

You compute E(t) for all t from t=0 to t=T and you average E over T

E = sum_i E(t_i)*(t_i+1 t_i)/T

THen you are left with E(x) at discrete location x_i.

Then you carry out an FFT on E(x_i) (i=1,N) and you obtain E(k).

If your entire domain is D and you have N points there, then delta x= D/N. k=1 correspond to a wavelenght lambda=D k=2 correspoinds to a wavelength D/2

etc...

k=N/2 corresponds to a wavelength 2D/N

In two dimension you will have different k's: k_x and k_y. To compute E(k_x) you need to average E(x,y) over t again but also over y so that you are left with E(x) only and again carry out an FFT. For the y dimension you need to average over x and get E(y). An FFT of E(y) brings you E(k_y).

In 3D you have k_x, k_y, k_z. In order to obtain E(x) you need to averages (discrete integration) over y and z, and of course t, then you can do the FFT for x. similarly for y and z.

The Energy spectrum is not as a function of frequency (time) but it is about wavelength with k=2*pi/lambda where lambda is the wavelength (spatial dimension, NOT time).

Look into Numerical Recepices. Check CFD books and see what they have on this CFD site.

It is not the time domain that you convert into frequencies but the space domain (after averaging over time).
This is totally wrong method.
wgq521 is offline   Reply With Quote

Old   March 3, 2015, 11:31
Default
  #4
Senior Member
 
romant's Avatar
 
Roman Thiele
Join Date: Aug 2009
Location: Eindhoven, NL
Posts: 374
Rep Power: 21
romant is on a distinguished road
If this is a totally wrong method, please supply the right method. More constructive criticism please.
nskelly and mezavalle like this.
__________________
~roman
romant is offline   Reply With Quote

Old   September 24, 2015, 16:22
Default
  #5
Senior Member
 
Julio Mendez
Join Date: Apr 2009
Location: Fairburn, GA. USA
Posts: 290
Rep Power: 18
juliom is on a distinguished road
Send a message via Skype™ to juliom
Dear Patrick;
I read your reply for the turbulent Energy Spectrum. I need to write my own script to do this in FORTRAN because I have the LES solution. From your explanation, do you have a textbook or paper which describes the procedure?? I am reading Wilcox book but I am still pretty loss!!!

Thanks!!!
juliom is offline   Reply With Quote

Old   September 24, 2015, 17:03
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,882
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
search for similar threads, this issue was already addressed
FMDenaro is offline   Reply With Quote

Old   September 24, 2015, 17:13
Default
  #7
Senior Member
 
Julio Mendez
Join Date: Apr 2009
Location: Fairburn, GA. USA
Posts: 290
Rep Power: 18
juliom is on a distinguished road
Send a message via Skype™ to juliom
Than you very much Mr. Filippo. I just wanted to have textbook or paper as a reference!!!
juliom is offline   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
how to solve a problem involving both laminar and turbulent flow seefd FLUENT 1 June 3, 2011 04:20
sgs turbulent kinetic energy blons Main CFD Forum 0 February 25, 2010 13:16
energy loss due to turbulent flow Al Mazdeh CFX 1 March 24, 2008 17:33
Energy Spectrum Emad Khalifa Main CFD Forum 3 June 30, 2003 17:03
the turbulent problem Daejin, Kang Main CFD Forum 0 November 16, 1998 11:52


All times are GMT -4. The time now is 18:08.