spectrum3d
algorithm
- calculates Energy spectrum of a 3D velocity field
- plots Energy spectrum over wave number
- plots Energy spectrum over eddy size
- calculates turbulent kinetic energy
- calculates Integral length scale
- calculates rms values of velocity
syntax
[ENVP1,S1,S,k,L11,RMS]=spectrum3d(x_velocity,y_velocity,z_velocity,L,nx,ny,nz)
example for input
[ENVP1,S1,S,k,L11,RMS]=spectrum3d(x_velocity,y_velocity,z_velocity,0.005,151,151,151);
input arguments
|
x_velocity, |
cube matrices, x-, y- and z-components of velocity |
|
L |
length of domain |
|
nx,ny,nz |
number of points of domain |
|
|
|
output arguments
|
ENVP1 |
vector, values of energy spectrum |
|
k |
vector, wavenumbers |
|
S1 |
turbulent kinetic energy |
|
S |
normalized turbulent kinetic energy |
|
L11 |
Integral length scale |
|
RMS |
standard deviation of velocity values |
|
plots |
energy cascade along wave number |
matlab tools required
- fft3, fftshift
toolbox tools required