Agreement between dissipation estimates: Difference between revisions
No edit summary |
No edit summary |
||
| Line 1: | Line 1: | ||
When two or more shear probes, in close proximity, are used to collect simultaneous data, the rate of dissipation derived from | When two or more shear probes, in close proximity, are used to collect simultaneous data, the rate of dissipation derived from such data will not agree exactly. | ||
Even for nearly flawless measurements, there will be disagreement for purely statistical reasons. | Even for nearly flawless measurements, there will be disagreement for purely statistical reasons. | ||
Measuring a turbulent shear is sampling a statistical process. | Measuring a turbulent shear is sampling a statistical process. | ||
| Line 17: | Line 17: | ||
<math> \mathrm{CF_{95}} = \sqrt{\varepsilon_1\,\varepsilon_2} \, \exp\left(\pm1.96\,\sigma_{\ln\varepsilon}/\sqrt{2} \right) \ \ .</math> | <math> \mathrm{CF_{95}} = \sqrt{\varepsilon_1\,\varepsilon_2} \, \exp\left(\pm1.96\,\sigma_{\ln\varepsilon}/\sqrt{2} \right) \ \ .</math> | ||
Thus, there is less than a 5% chance that the ratio of two dissipation estimates is outside of the range of | Thus, there is less than a 5% chance that the ratio of two simultaneous dissipation estimates is outside of the range of | ||
<math> | <math> \exp\left(\pm1.96\,\sigma_{\ln\varepsilon}\,\sqrt{2} \right) \ \ .</math> | ||
Because the rate of dissipation is not identical for a pair of probes, the estimate of the standard deviation, <math>\sigma_{\ln\varepsilon}</math>, will also differ. | |||
But, only slightly because of the quarter-power dependence of the Kolmogorov length on the rate of dissipation. | |||
So, one could use either the smaller of the two standard deviations or their average for the testing of the dissipation ratios. | |||
If there are more than two simultaneous dissipation estimates, then they should be sorted into ascending order. | |||
The ratio test is first applied to the smaller and the largest dissipation estimate. | |||
If this pair pass the test and all other pairs will also pass the ratio test. | |||
If this pair does not pass the test, then the larger of the two should be flagged, and the test applied to the next largest (and the the smallest) until a pair passes the test, or all pairs have been tested. | |||
It is most likely that the larger of a pair of estimates is the erroneous one because signal contamination and other data flaws usually act to increase the variance of the shear=probe signal. | |||
Very large dissipation rates require a fit to the inertial subrange in order to estimate the rate of dissipation because the shear probe will not resolve the shear variance. | |||
In the inertial subrange, the spectrum of shear is | |||
<math> \Psi(k) = A\, \varepsilon^{2/3}\left(2\pi\right)^2 k^{1/3}</math> | |||
where the factor of <math>A</math> is approximately 0.2, depending on the model one wishes to use. | |||
Revision as of 23:07, 23 May 2022
When two or more shear probes, in close proximity, are used to collect simultaneous data, the rate of dissipation derived from such data will not agree exactly. Even for nearly flawless measurements, there will be disagreement for purely statistical reasons. Measuring a turbulent shear is sampling a statistical process. The sample variance will differ from the population variance and this difference will reduce with increasing length of data in the sample. The sampling uncertainty is distributed lognormally with a variance of
[math]\displaystyle{ \sigma^2_{\ln\varepsilon} = \frac{5.5}{1 + \left(\hat{L}_f/4\right)^{7/9}}\ \ ,\ \ \hat{L}_f \equiv \hat{L} V_f^{3/4} = \frac{L}{L_K} V_f^{3/4} }[/math]
where [math]\displaystyle{ L_K=\left(\nu^3/\varepsilon \right)^{1/4} }[/math] is the Kolmogorov length, and [math]\displaystyle{ V_f }[/math] is the fraction of the shear variance that is resolved by terminating the spectral integration at an upper wavenumber of [math]\displaystyle{ k_u }[/math][1] .
The 95% confidence interval on an individual dissipation estimate, [math]\displaystyle{ \varepsilon_1 }[/math] is thus
[math]\displaystyle{ \mathrm{CF_{95}}(\varepsilon_1) = \varepsilon_1\, \exp\left(\pm1.96\,\sigma_{\ln\varepsilon} \right) \ \ . }[/math]
The 95% confidence interval for the geometric mean of a pair of dissipation estimates is
[math]\displaystyle{ \mathrm{CF_{95}} = \sqrt{\varepsilon_1\,\varepsilon_2} \, \exp\left(\pm1.96\,\sigma_{\ln\varepsilon}/\sqrt{2} \right) \ \ . }[/math]
Thus, there is less than a 5% chance that the ratio of two simultaneous dissipation estimates is outside of the range of
[math]\displaystyle{ \exp\left(\pm1.96\,\sigma_{\ln\varepsilon}\,\sqrt{2} \right) \ \ . }[/math]
Because the rate of dissipation is not identical for a pair of probes, the estimate of the standard deviation, [math]\displaystyle{ \sigma_{\ln\varepsilon} }[/math], will also differ.
But, only slightly because of the quarter-power dependence of the Kolmogorov length on the rate of dissipation.
So, one could use either the smaller of the two standard deviations or their average for the testing of the dissipation ratios.
If there are more than two simultaneous dissipation estimates, then they should be sorted into ascending order. The ratio test is first applied to the smaller and the largest dissipation estimate. If this pair pass the test and all other pairs will also pass the ratio test. If this pair does not pass the test, then the larger of the two should be flagged, and the test applied to the next largest (and the the smallest) until a pair passes the test, or all pairs have been tested.
It is most likely that the larger of a pair of estimates is the erroneous one because signal contamination and other data flaws usually act to increase the variance of the shear=probe signal.
Very large dissipation rates require a fit to the inertial subrange in order to estimate the rate of dissipation because the shear probe will not resolve the shear variance. In the inertial subrange, the spectrum of shear is
[math]\displaystyle{ \Psi(k) = A\, \varepsilon^{2/3}\left(2\pi\right)^2 k^{1/3} }[/math]
where the factor of [math]\displaystyle{ A }[/math] is approximately 0.2, depending on the model one wishes to use.
References
- ↑ Lueck, R. G., 2022a: The statistics of oceanic turbulence measurements. Part 1: Shear variance and dissipation rates. J. Atmos. Oceanic Technol., –, in press, doi:--.
