Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Equivalent refractive-index structure constant of non-Kolmogorov turbulence: comment

Open Access Open Access

Abstract

The paper by Y. Li, W. Zhu, X. Wu, and R. Rao entitled “Equivalent refractive-index structure constant of non-Kolmogorov turbulence,” Opt. Express 23(18), 23004 (2015). [CrossRef]   relates the non-Kolmogorov turbulence structure constant to the classical structure constant for Kolmogorov turbulence by imposing equality of their respective structure functions at large separation distances, higher than the outer scale. As opposed to previous attempts to relate the two structure constants, the approach of Li et al. is anchored on a measurable meteorological parameter, the outer scale. The error lies in the fact that the authors have used a default Kolmogorov structure function with an infinite outer scale. A subsequent assumption of a finite outer scale is not compatible with the initial assumption. In this paper we show the correct procedure to obtain the relationship between the non-Kolmogorov and Kolmogorov structure constants which is based on an explicitly finite outer scale used throughout all calculations.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Li et al. [1] use the Kolmogorov and non-Kolmogorov structure functions of the refractive index evaluated at the separation set to ${L_0}$ to calculate the variance of the index of refraction fluctuations, $\sigma _n^2$ (see Eq. (2) in this paper). This quantity is then kept fixed when comparing quantities resulting from predictions of non-Kolmogorov theories. The motivation for the approach of Ref. [1] is appealing, but the main deficiency of it is that it starts with the (implicit) assumption of infinite outer scale (Eq. (15) in Ref. [1]), only to introduce a finite outer scale in the subsequent equations. The approach requires starting from von Kármán model, where ${L_0}$ is explicitly present, in the first place. In this paper we rectify this inaccuracy.

2. Anchoring non-Kolmogorov turbulence on the outer-scale in the von Kármán spectrum

Let us focus on the following non-Kolmogorov power spectrum valid only for isotropic turbulence with zero inner scale [2,3]:

$${\Phi _n}({\kappa ,\alpha } )= A(\alpha )\cdot \tilde{C}_n^2 \cdot {\left[ {{\kappa^2} + {{\left( {\frac{{2\pi }}{{{L_0}}}} \right)}^2}} \right]^{ - \frac{\alpha }{2}}},3 < \alpha < 5$$
where: $A(\alpha )= \frac{{\Gamma ({\alpha - 1} )}}{{4{\pi ^2}}}\cos \left( {\alpha \cdot \frac{\pi }{2}} \right)$, $\kappa = |{\vec{\kappa }} |\equiv ({{\kappa_x},{\kappa_y},{\kappa_z}} )$ is the spatial wavenumber (note that for isotropic turbulence ${\kappa _x} = {\kappa _y} = {\kappa _z}$), α is the power law exponent, $\tilde{C}_n^2 = \beta (\alpha )\cdot C_n^2$ is the generalized structure parameter with units $[{{m^{3 - \alpha }}} ]$, $\beta (\alpha )$ is a constant depending on $\alpha$ and has units $[{{m^{11/3 - \alpha }}} ]$, symbol $\Gamma (x )$ denotes the Gamma function and ${L_0}$ is the outer scale of turbulence. When the power law assumes value $\alpha = 11/3$ the generalized structure parameter reduces to the structure parameter $C_n^2$ with units $[{{m^{ - 2/3}}} ]$. Note also that when $\alpha = 11/3$ Eq. (1) reduces to the von Kármán power spectrum [3].

It is shown in Ref. [3] that for Kolmogorov power law the structure function of the fluctuations of the index of refraction associated with the von Kármán power spectrum, Eq. (1) with $\alpha = 11/3$, and at very large separation distance, $\rho \gg {L_0}$ takes the form

$${D_n}({\rho \gg {L_0}} )\cong 1.0468 \cdot {\left( {\frac{{2\pi }}{{{L_0}}}} \right)^{ - \frac{2}{3}}} \cdot C_n^2 = 2\sigma _n^2$$

Similarly, for the spectrum with non-Kolmogorov power law shown in Eq. (1), the structure function at very large separation distance, $\rho \gg {L_0}$ can be expressed as

$${D_n}({\rho \gg {L_0}} )\cong 8\pi \cdot A(\alpha )\cdot \tilde{C}_n^2 \cdot {\left( {\frac{{2\pi }}{{{L_0}}}} \right)^{ - \alpha }} \cdot \int\limits_0^\infty {{\kappa ^3}} \cdot {\left[ {{{\left( {\kappa \cdot \frac{{{L_0}}}{{2\pi }}} \right)}^2} + 1} \right]^{ - \frac{\alpha }{2}}}\frac{{d\kappa }}{\kappa }$$

To solve integrals of the form as in Eq. (3), Sasiela uses the Mellin transform [3]

$$h(x )\to H(s )\equiv {\rm M}[{h(x )} ]\equiv \int\limits_0^\infty {\frac{{dx}}{x}h(x ){x^s}}$$
and two of its proprieties
$$h({{x^p}} )\to \frac{{H\left( {\frac{s}{p}} \right)}}{{|p |}},p \ne 0\;\;\;\textrm{and}\;\;\;{({1 + x} )^{ - p}} \to \frac{{\Gamma ({s,p - s} )}}{{\Gamma (p )}}.$$

Using Eqs. (4) and (5), Eq. (3) can be expressed in closed form as follows

$${D_n}({\rho \gg {L_0}} )\cong 8\pi \cdot A(\alpha )\cdot \tilde{C}_n^2 \cdot {\left( {\frac{{2\pi }}{{{L_0}}}} \right)^{3 - \alpha }} \cdot \frac{{\Gamma \left[ {\frac{3}{2}} \right] \cdot \Gamma \left[ {\frac{{\alpha - 3}}{2}} \right]}}{{2 \cdot \Gamma \left[ {\frac{\alpha }{2}} \right]}} = 2\sigma _n^2$$

Note that Eq. (6) reduces to Eq. (2) when $\alpha = 11/3$ as it should be for the Kolmogorov turbulence case. Additionally observe that both structure functions in Eqs. (2) and (6), at large separation distances, $\rho \gg {L_0}$ saturate to double the value of the variance of the refractive index fluctuations, $2\sigma _n^2.$ The variance is by definition the covariance ${B_n}(\rho )$ calculated at zero separation distance, $\rho = 0$. For isotropic and homogeneous turbulence the variance or total power is $\sigma _n^2 = {B_n}(0 )= 4\pi \cdot \int\limits_0^\infty {{\kappa ^2}} {\Phi _n}({\kappa ,\alpha } )d\kappa$. To find the relation $\beta (\alpha )$ between the non-Kolmogorov and Kolmogorov refractive index structure constants, we impose the equality of the two structure functions, Eq. (2) and Eq. (6), and we deduce the main result of this paper

$$\beta (\alpha )= \frac{{\tilde{C}_n^2}}{{C_n^2}} = 0.0940 \cdot \frac{{\Gamma \left[ {\frac{\alpha }{2}} \right]}}{{\Gamma \left[ {\frac{{\alpha - 3}}{2}} \right]}}\frac{{{{({2\pi } )}^{\alpha - \frac{{11}}{3}}}}}{{A(\alpha )}}{L_0}^{\frac{{11}}{3} - \alpha }$$

The refractive index structure constant $\tilde{C}_n^2$ can be obtained directly from $C_n^2$ once the outer scale is fixed and the power law exponent $\alpha$ is known. Note that the factor $\beta (\alpha )\cdot A(\alpha )$ represents the turbulence scaling due to non-Kolmogorov turbulence. In Fig. 1 we plot the 3D power spectrum, Eq. (1), multiplied by $4\pi \cdot {\kappa ^2}$ (to obtain the 1D version) for several power laws with outer scale fixed at 5 m and $C_n^2 = 1.7 \cdot {10^{ - 14}}{\textrm{m}^{ - 2/3}}$. It can be easily seen in Fig. 1 that the total power remains constant while the power law is changed because we defined $\beta (\alpha )$ to keep the total power of spectrum equal to the variance, $\sigma _n^2$ for any power law $\alpha$.

 figure: Fig. 1.

Fig. 1. 1D power spectrum (Eq. (1) multiplied by $4\pi \cdot {\kappa ^2}$) for different power laws $\alpha$ with fixed outer scale, ${L_0} = 5\,\textrm{m}$ and $C_n^2 = 1.7 \cdot {10^{ - 14}}{\textrm{m}^{ - 2/3}}$. The conservation of the total power for each plot is clearly visible.

Download Full Size | PDF

3. Summary

To summarize, in this paper we found a new relation between the structure constant in the von Kármán spectrum and the structure constant in the von Kármán spectrum with non-Kolmogorov power law. The approach we used is based on a measurable physical parameter, the outer scale, and it ensures that the total power of the fluctuations of the refractive index is conserved when the power law changes. Only normalized spectra allow for a meaningful comparison of their effects on various observables. Because a change in the power law brings with itself a change in how the power is distributed across spatial wavenumbers, different observables will react differently to different normalizations. Our approach provides natural and rigorous normalization. To our knowledge, all previous attempts to compare Kolmogorov and non-Kolmogorov turbulence, including Ref. [1] fail to conserve power and as such are not physical.

Acknowledgments

This work was sponsored by WTD 91 (Technical Center of Weapons and Ammunition) of the Federal Defence Forces of Germany – Bundeswehr in the project ABU-SLS.

Disclosures

The authors declare no conflicts of interest.

Data availability

No data were generated or analyzed in the presented research.

References

1. Y. Li, W. Zhu, X. Wu, and R. Rao, “Equivalent refractive-index structure constant of non-Kolmogorov turbulence,” Opt. Express 23(18), 23004–23012 (2015). [CrossRef]  

2. I. Toselli, L. C. Andrews, R. L. Phillips, and V. Ferrero, “Angle of arrival fluctuations for free space laser beam propagation through non Kolmogorov turbulence,” Proc. SPIE 6551, 65510E (2007). [CrossRef]  

3. R. J. Sasiela, Electromagnetic wave propagation in turbulence. Evaluation and application of Mellin transforms, SPIE, 2nd ed. (2007).

Data availability

No data were generated or analyzed in the presented research.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (1)

Fig. 1.
Fig. 1. 1D power spectrum (Eq. (1) multiplied by $4\pi \cdot {\kappa ^2}$) for different power laws $\alpha$ with fixed outer scale, ${L_0} = 5\,\textrm{m}$ and $C_n^2 = 1.7 \cdot {10^{ - 14}}{\textrm{m}^{ - 2/3}}$. The conservation of the total power for each plot is clearly visible.

Equations (7)

Equations on this page are rendered with MathJax. Learn more.

Φ n ( κ , α ) = A ( α ) C ~ n 2 [ κ 2 + ( 2 π L 0 ) 2 ] α 2 , 3 < α < 5
D n ( ρ L 0 ) 1.0468 ( 2 π L 0 ) 2 3 C n 2 = 2 σ n 2
D n ( ρ L 0 ) 8 π A ( α ) C ~ n 2 ( 2 π L 0 ) α 0 κ 3 [ ( κ L 0 2 π ) 2 + 1 ] α 2 d κ κ
h ( x ) H ( s ) M [ h ( x ) ] 0 d x x h ( x ) x s
h ( x p ) H ( s p ) | p | , p 0 and ( 1 + x ) p Γ ( s , p s ) Γ ( p ) .
D n ( ρ L 0 ) 8 π A ( α ) C ~ n 2 ( 2 π L 0 ) 3 α Γ [ 3 2 ] Γ [ α 3 2 ] 2 Γ [ α 2 ] = 2 σ n 2
β ( α ) = C ~ n 2 C n 2 = 0.0940 Γ [ α 2 ] Γ [ α 3 2 ] ( 2 π ) α 11 3 A ( α ) L 0 11 3 α
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.