# how to compute the tranfer spectrum of TKE

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

 October 26, 2012, 15:24 how to compute the tranfer spectrum of TKE #1 New Member   Join Date: Oct 2012 Posts: 1 Rep Power: 0 Hello and good day to everyone, I'm new here on the forum and would like to say hello to everyone on board! I have a question for the turbulence experts here: I'm computing homogeneous isotropic turbulence (HIT) in a box with a code in physical space (FV method) and I am using an external FFT tool to compute the energy spectra from E(k) = sum over all shell contributions of u*conj(u). The spectra and dissipation rate of the (HIT) comes out allright, but now I'm stuck at how to compute the transfer spectra T... For HIT, we have dE/dt = T - 2nu k^2 E(k), where T is the transfer between the modes and the integral of T over all k is zero.... Now how do I compute the T starting from my solution in physical space? Could I do the following: 1) my code is set up in conservation form, so U_t + F_x+G_y+H_z+viscousFluxes=0. 2) to get the equation for E, we multiply by u=(u,v,w) as a scalar product 3) that would give me u*F_x+u*G_y+u*H_z, where the first term would be something like u*F(2)_x+v*F(3)_x+w*F(4)_x (where 2:4 is the entry corresponding to the u,v,w momentum). 4) Then, I would take an FFT of this whole term, to get T(k1,k2,k3) and then sum up the T according to k=sqrt(k1^2 + k2^2 +k3^2) 5) then plot T(k) over k This is what I have done so far, but something seems to be wrong... the T(k) integral does not sum up to zero, and I'm getting some strong positive values in the low wavenumbers.... is there anything wrong with my approach? Is there a better way to do it? How is this computation normally done? I'm sorry if this might be a stupid question, but I'm new and not so experienced in turbulence! thank you very much! Kind regards!

 Thread Tools Display Modes Linear Mode

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

 Similar Threads Thread Thread Starter Forum Replies Last Post rg12g09 STAR-CCM+ 0 December 6, 2011 10:51 liu OpenFOAM Running, Solving & CFD 6 November 1, 2011 02:00 Walid Main CFD Forum 2 December 2, 2005 09:29 Emad Khalifa Main CFD Forum 3 June 30, 2003 16:03 Christian Paukert Main CFD Forum 2 March 14, 2001 20:34

All times are GMT -4. The time now is 06:05.