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

Plot Turbulence Spectrum Using 2D Spatial FFT

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 21, 2021, 08:31
Question Plot Turbulence Spectrum Using 2D Spatial FFT
  #1
New Member
 
Jeffrey Johnston
Join Date: Oct 2020
Location: Belfast, Northern Ireland
Posts: 21
Rep Power: 6
NotDrJeff is on a distinguished road
Hello!

I am trying to replicate the results of a study (here) that performs LES of an atmospheric boundary layer. One of the results produced is a turbulence spectrum. I want to recreate this, but I'm a bit lost.

Here's the description in the paper:

Quote:
Figure 3 shows velocity spectra taken at 90 m (the turbine's hub height) above the surface. The spectra were obtained by taking two-dimensional fast Fourier transforms in space in the horizontal plane at 90 m. We then summed the energy in wavevectors of equal magnitude (i.e. we integrated the energy between concentric rings in the two-dimensional wavespace) to produce the spectra shown in Figure 3. Spectra are shown for both the vertical component and the magnitude of the horizontal components of the velocity vector. A line of −5/3 slope is shown for reference to judge LES's ability to reproduce the cascade of energy from larger to smaller scales.
So here is my first point of confusion. I thought the turbulence energy spectrum had to be based on kinetic energy. That is, I thought I had to calculate \frac{1}{2}U_{i}U_{i} before performing the FFT. But in this paper, the (horizontal) velocity magnitude is calculated instead, \sqrt{U_{i}U_{i}}. So, have I misunderstood the turbulence spectrum? Does the -5/3 slope region exist in a a spectral transform of the magnitude of the velocity? And will it exist for only the horizontal magnitude? Do we have to assume that the vertical component is negligible here (a reasonable assumption in this case)?

My second point of confusion is how to obtain a single 2-axis plot from a 2D fft. Performing a fft of the 2D array of velocity magnitudes produces a 2D array of velocity intensities. I assume that this resultant array is what the we call "wavespace"? But what is a wavevector? How do I find ones of equal magnitude and sum them? And what are the concentric rings all about?

Many thanks,
NotDrJeff is offline   Reply With Quote

Old   October 21, 2021, 13:23
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by NotDrJeff View Post
Hello!

I am trying to replicate the results of a study (here) that performs LES of an atmospheric boundary layer. One of the results produced is a turbulence spectrum. I want to recreate this, but I'm a bit lost.

Here's the description in the paper:



So here is my first point of confusion. I thought the turbulence energy spectrum had to be based on kinetic energy. That is, I thought I had to calculate \frac{1}{2}U_{i}U_{i} before performing the FFT. But in this paper, the (horizontal) velocity magnitude is calculated instead, \sqrt{U_{i}U_{i}}. So, have I misunderstood the turbulence spectrum? Does the -5/3 slope region exist in a a spectral transform of the magnitude of the velocity? And will it exist for only the horizontal magnitude? Do we have to assume that the vertical component is negligible here (a reasonable assumption in this case)?

My second point of confusion is how to obtain a single 2-axis plot from a 2D fft. Performing a fft of the 2D array of velocity magnitudes produces a 2D array of velocity intensities. I assume that this resultant array is what the we call "wavespace"? But what is a wavevector? How do I find ones of equal magnitude and sum them? And what are the concentric rings all about?

Many thanks,

Often one performs the FFT of a velocity component, the corresponding Fourier coefficient is "squared" and you have dimensionally one component of the kinetic energy. There are theoretical reasons to deduce that.

In the paper the 2D plot is reduced to a 1D plot along the unique k wavenumber computed as described, for different rings.
FMDenaro is offline   Reply With Quote

Reply

Tags
fft, les, spectral analysis, turbulence, waves

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 plot time dependent energy spectrum Falu Main CFD Forum 21 January 16, 2021 09:41
Plot the Energy Spectrum Elyas_Mosibat OpenFOAM Post-Processing 14 June 3, 2019 14:46
Amplitude spectrum of temperature fluctuations selvam2487 CFX 0 March 16, 2014 14:59
energy spectrum using FFT kharnabnew FLUENT 0 November 18, 2009 04:16
Discussion: Reason of Turbulence!! Wen Long Main CFD Forum 3 May 15, 2009 10:52


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