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

Offset when calculating TKE using autocorrelation

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 22, 2021, 06:32
Default Offset when calculating TKE using autocorrelation
  #1
New Member
 
Fab
Join Date: Nov 2020
Posts: 25
Rep Power: 5
Falu is on a distinguished road
Hey CFD Community,

I simulated a free jet stream in Ansys CFX and want to plot the Turbulence Kinetic Energy of my Simulation. In my Simulation I computed a variable which tracks the unresolved Turbulent Kinetic Energy. Now Im using Matlab to calculate the TKE with the components of the velocity vector. The two graphs have a similar trend, but you can see an offset in the y-values. I attached the two graphs below.

I guess there is something wrong with my calculation, am I missing something? The relevant part is the following:

for Wort=1:words
Bereich=Liste(Wort);
for d=1:1:datapoints
spalte=d+(Wort-1)*datapoints;
dstring=int2str(d);

energy_vel(:,spalte)=(P2Y_fluc_win_han(:,spalte).* conj(P2Y_fluc_win_han(:,spalte)))/(N^2);

end
end

for d=1:1:datapoints
energy_tot(:,d)=energy_vel(:,d)+energy_vel(:,d+dat apoints)+energy_vel(:,d+datapoints*2);
end
Attached Images
File Type: png Total_EnergySpectrum_Frequency10.png (60.2 KB, 8 views)
File Type: png resTKEmon_10.png (59.1 KB, 4 views)
Falu 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 set periodic boundary conditions Ganesh FLUENT 15 November 18, 2020 06:09
Long output in terminal. ssa_cfd OpenFOAM Running, Solving & CFD 1 March 18, 2019 05:25
Calculating TKE in URANS for comparison with experiment pbachant OpenFOAM Post-Processing 0 April 27, 2015 12:30
''unknown radialModelType type Gidaspow'' PROBLEM WITH THE BED TUTORIAL AndoniBM OpenFOAM Running, Solving & CFD 2 March 25, 2015 18:44
How to update polyPatchbs localPoints liu OpenFOAM Running, Solving & CFD 6 December 30, 2005 17:27


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