OPEN_FPE_IFT icon indicating copy to clipboard operation
OPEN_FPE_IFT copied to clipboard

The integral length scales not working for me

Open KoalaMau opened this issue 4 years ago • 1 comments

Hi,

I'm trying to use the code for my data but I encounter some difficulties, when I try to select the lowest and highest frequency of constant range in Hz I receive the error that is seen in the second figure. Can you help me understand what is the problem that I facing? Error

KoalaMau avatar Aug 26 '21 19:08 KoalaMau

Hi Tiago,

sorry for the late reply. I was on vacation and only now read your message on researchgate. From your screenshot I can see that there is NO problem in fitting the constant range. Rather there is a problem in line 241 in the script "length_scales.m". So the problem occurs when determining the dissipation from structure functions. image

Since here the mean value amongst 5 points closest to peak value of dissipation in inertial range is determined, this error arises I suspect in connection with a very small selected value for "increment_bin" in combination with line 203: r = unique(round(logspace(log10((1./Fs).m_data10^6),log10(20IntegralLengthScale_210^6),increment_bin))./10^6,'stable');

You could try to increase the value "increment_bin". If this is not the solution, then you can also just send me the evaluated.mat file. Here all data are saved before the error occurs and I can look at the error again more precisely by using your data.

Best, André