Computer simulations of optical turbulence in the weak- and strong-scattering regime: angle-of-arrival fluctuations obtained from ray optics and wave optics

David G. Voelz; Erandi Wijerathna; Andreas Muschinski; Xifeng Xiao

doi:10.1117/1.OE.57.10.104102

5 October 2018 Computer simulations of optical turbulence in the weak- and strong-scattering regime: angle-of-arrival fluctuations obtained from ray optics and wave optics

David G. Voelz, Erandi Wijerathna, Andreas Muschinski, Xifeng Xiao

Author Affiliations +

Optical Engineering, Vol. 57, Issue 10, 104102 (October 2018). https://doi.org/10.1117/1.OE.57.10.104102

Abstract

It is known that certain geometrical-optics predictions often agree well with optical turbulence field observations even though theoretical constraints for ignoring diffraction may be violated. Geometrical optics assumptions can simplify analyses, and ray optics can significantly reduce simulation computation time. Here, an investigation into angle-of-arrival fluctuations is presented involving wave optics and geometrical (ray) optics computer simulations of a plane wave of visible light propagating through a turbulent refractive-index field. The simulation and Rytov-based theory results for the variances of aperture-filtered angle-of-arrival fluctuations generally agree well for weak scattering (Rytov variance, σR2≲0.2), but for increasing Rytov variance, the simulation results demonstrate a positive slope that can be significantly shallower than that predicted by the theory. For weak-to-moderate scattering regimes ( σR2≲2.67), a comparison of the ray and wave results show they match for aperture diameters greater than about two Fresnel lengths. This result is consistent with a previous theoretical analysis by Cheon and Muschinski. For the strongest scattering case studied ( σR2=26.7), the wave and ray simulations match for aperture diameters greater than about 10 Fresnel lengths. For smaller apertures, we attribute the disparity between the wave and ray simulation results to a Fresnel filtering effect.

1. Introduction

Understanding the physics of optical wave front distortions caused by turbulence in the atmospheric refractive-index field is essential for the design and operation of optical systems used in free-space optical communication, surveillance, navigation, remote sensing, astronomy, and directed-energy technologies.¹^–⁶ If particle scattering, absorption, coupling, and depolarization effects are negligible then there are two phenomena that determine the optical field: refraction and diffraction. The combined effects of refraction and diffraction lead to phase-gradient fluctuations (angle-of-arrival fluctuations) and intensity fluctuations (scintillation), both of which limit the performance of a variety of optical systems.

The analysis of propagation scenarios through both deterministic and stochastic refractive-index fields may be substantially simplified if diffraction effects can be neglected, and considerable theoretical effort has been devoted to establishing criteria for the negligibility of diffraction effects.¹^–³^,⁷^,⁸ These theoretical constraints, however, tend to be excessively conservative, as pointed out by Strohbehn 50 years ago.⁹ In other words, certain geometrical-optics predictions often agree well with field observations, but it remains unclear, from the theoretician’s perspective, why that is so.

In addition to simplifying theoretical analyses, geometrical optics and ray tracing can also substantially reduce computation time for numerical simulations of propagation through turbulence. The utility of geometrical and ray optics is exemplified in a variety of publications from the last several decades. For example, Churnside and Lataitis¹⁰ presented a simple geometrical optics equation for the variance of beam displacements caused by propagation through weak turbulence, and geometrical optics simulation models have been developed by Yuksel et al.¹¹ to investigate effects of aperture averaging on intensity fluctuations. More recently, Lachinova et al.¹² have compared the performance of a brightness function-based numerical simulation method with a wave optics-based method for incoherent imaging through turbulence. The brightness function at the receiving pupil is estimated using ray-like trajectories that are perturbed by the turbulent refractive index gradients. These authors found that the brightness approach can reduce computation times by several orders of magnitude. Other recent examples include ray-tracing methods applied for simulating laser beam-steering effects in turbulent media,¹³ the derivation of turbulence strength profiles from radiometer temperature measurements with a geometrical ray-tracing analysis,¹⁴ and the application of ray tracing through phase screens for the simulation of images from optical survey telescopes.¹⁵

Since the 1980s, powerful phase-screen techniques have been developed¹⁶^–¹⁹ that allow optical propagation through atmospheric turbulence to be accurately simulated well into the “strong-scattering” regime, that is, for scenarios where the negligibility of diffraction effects and the “weak-scattering” assumption become invalid, or at least questionable. Such simulations can be used to test and refine the traditional, often excessively conservative, criteria for the negligibility of diffraction effects.

Here, we present and discuss computer simulations of plane waves of visible light ( $λ = 500 nm$ ) propagating through fully turbulent, homogeneous, and isotropic refractive-index fields. We study the optical angle-of-arrival (AOA) fluctuations by means of both a wave-optics propagation simulation and a ray-tracing simulation through a turbulent path that is characterized by a sequence of phase screens with Kolmogorov statistics. The path length assumed is $L = 2 km$ and the refractive-index structure parameter $C_{n}^{2}$ ranges from $10^{- 16} m^{- 2 / 3}$ to $10^{- 13} m^{- 2 / 3}$ , which corresponds to Rytov variances $σ_{R}^{2}$ that span from 0.0267 (weak scattering) to 26.7 (strong scattering), where $σ_{R}^{2} = 1.23 C_{n}^{2} k^{7 / 6} L^{11 / 6}$ and $k = 2 π / λ$ is the optical wave number. For the two simulation methods, we examine the variance and histograms of the aperture-filtered (aperture-averaged) AOA fluctuations for aperture diameters ranging from 0.6 to 13 Fresnel lengths. We compare the simulation results with theoretical predictions²⁰ based on the Rytov theory, which accounts for both refraction and diffraction effects but does not capture strong-scattering effects. Our specific interests are: (1) to better understand under what conditions the Rytov theory for AOA fluctuations becomes invalid, (2) to reveal the behavior of the AOA fluctuations in strong scattering conditions, and (3) to investigate the consequences of ignoring diffraction and the applicability of ray tracing for modeling AOA fluctuations.

2. Approach

A plane wave is propagated through a set of split-step phase screens that model atmospheric turbulence [Fig. 1(a)] and similarly, a bundle of parallel rays is traced through the same screens [Fig. 1(b)]. The wave optics simulation provides a coherent wave result that accounts for refraction, diffraction, and interference. On the other hand, the ray-tracing simulations provide the trajectories of a large number of individual rays, where the phases and amplitudes are left undefined.

Fig. 1

Illustration of the split-step phase screens that model turbulence in (a) wave and (b) ray simulations. $L$ is the propagation path length and $D$ is the diameter of the aperture that is applied in the observation plane for the aperture-averaged AOA.

The simulation parameters of interest are defined in Table 1 and the values used for the simulations are listed. The turbulence along the path is assumed to be statistically homogeneous and isotropic. The path is divided into equal length segments and a phase screen is positioned at the center of each segment. The screens are created with the Fourier filtering method and have Kolmogorov statistics. To model an infinite outer scale, a random tilt component is added to each screen that compensates for turbulent wave-front tilts with length scales greater than the grid width (see Ref. 18 and references therein). Zero inner scale is assumed for the turbulence spectrum; however, the period of the highest spatial frequency of a screen is twice the grid spacing, so the modeled inner scale is no smaller than 2 mm.

Table 1

Simulation parameter values.

Parameter	Notation	Value
Grid physical size	$S \times S$	$1 m \times 1 m$
Grid samples		$1000 \times 1000$
Grid spacing	$Δ x$	1 mm
Wavelength	$λ$	$0.5 μ m$
Path length	$L$	2 km
Fresnel length	$\sqrt{λ L}$	3.2 cm
Aperture diameter	$D$	0.02 to 0.4 m
Turbulence spectrum		Kolmogorov
Structure parameter	$C_{n}^{2}$	$10^{- 16}$ to $10^{- 13} m^{- 2 / 3}$
Rytov variance	$σ_{R}^{2}$	0.0267 to 26.7
Number of phase screens	$N$	6 (for $10^{- 16} m^{- 2 / 3}$ ); 6 (for $10^{- 15} m^{- 2 / 3}$ ); 10 (for $10^{- 14} m^{- 2 / 3}$ ); 22 (for $10^{- 13} m^{- 2 / 3}$ )

The number of turbulence screens $N$ required along the path is an integer and should be greater than ${(10 σ_{R}^{2})}^{6 / 11}$ , which ensures each screen corresponds to weak scattering, and therefore phase-only screens are adequate.¹⁹ The minimum value of $N$ for the strongest scattering case is 22 screens. Fewer screens can be used for the weaker scattering cases, which reduces computational time. However, we implemented more than the required minimum number of screens for the weak scattering cases to provide a better representation of the continuous volume turbulence that is assumed along the path.

The grid parameters are chosen to provide a grid spacing of $Δ x = λ L / S$ , where $S$ is the grid side length. This sampling choice has been shown to make artifacts that arise from the discrete Fourier transform implementation of the Fresnel diffraction integral to be sufficiently small to be negligible.²¹ The resulting grid spacing $Δ x$ of 1 mm allows about $3 \times$ sampling of the minimum Fried parameter ( $r_{0}$ ) value encountered in our study. The propagation length $L$ , wavelength $λ$ , and four choices of turbulence structure parameter $C_{n}^{2}$ ( $10^{- 16}$ , $10^{- 15}$ , $10^{- 14}$ , and $10^{- 13} m^{- 2 / 3}$ ) correspond to Rytov variances $σ_{R}^{2}$ of 0.0267, 0.267, 2.67, and 26.7, which range in scattering strength from weak ( $σ_{R}^{2} ≪ 1$ ) to strong ( $σ_{R}^{2} > 1$ ). Corresponding Fried diameter [ $r_{0} = {(0.423 k^{2} C_{n}^{2} L)}^{- 3 / 5}$ ] values for the four cases are calculated to be 21.3, 5.3, 1.33, and 0.33 cm.

Ultimately, our interest is to examine the variance and histogram of the fluctuations of the aperture-averaged AOA. This requires running the simulation many times for independent turbulence realizations, obtaining the aperture-averaged AOA for each realization, and then computing the variance.

2.1.

Wave Optics Simulation

The initial plane wave at the entrance plane, $z = 0$ , is a unit-amplitude, zero-phase field that is propagated in steps between the entrance plane and the observation plane using a discretized form of the Fresnel diffraction integral.²¹ At the $i$ ’th intermediate plane, the phase screen $ψ_{i}$ is applied in the usual way as ${u^{'}}_{i} = u_{i} \exp (j ψ_{i})$ , where $u_{i}$ is the incident field and ${u^{'}}_{i}$ is the field that exits the plane. The discrete Fourier transform-based simulation produces a periodicity of the electric field in the $x$ - $y$ plane across the grid boundaries. In addition, the compensating random tilt component in the turbulence phase screens creates some “ringing” artifacts in the receiving plane near the grid boundaries. To minimize unwanted periodicity and edge effects, the aperture is only applied to the center $0.4 \times 0.4$ fraction of the receiving plane grid.

One approach to find the aperture-averaged AOA for a given turbulence realization is to numerically compute the field phase gradients at the observation plane grid points and average the results within the aperture region. However, direct calculation of the phase gradients at the grid points is problematic for strong scattering because of branch points (implying discontinuous phase changes from $2 π$ to 0, or 0 to $2 π$ ) that appear in the wave front and correspond to theoretically infinite phase gradients.

For the results presented here, we use an approach where the wave front within the aperture for a given turbulence realization is numerically focused with a focus transmittance function.²¹ The $x$ - and $y$ -direction centroid positions of the intensity spot at the focal plane are found and the aperture-averaged AOA is the angle defined by the ratio of the centroid displacement (from the optical axis) and the focal length. This AOA is computed for many independent turbulence realizations and the variance is calculated. This approach has the advantage that it is representative of a practical AOA measurement that can be made with a lens and camera system.²² It also effectively weighs the angle contributions by the spatial field amplitude in the pupil, where, for example, the amplitude is vanishing near a branch point so the sometimes extremely large gradients associated with these features are not significant for the AOA variance calculation.

2.2.

Ray Tracing Simulation

For the ray simulation, as shown in Fig. 1(b), the source is comprised of a set of rays where each ray is positioned at a grid point corresponding to the wave simulation and is initially directed parallel to the optical axis. For a given turbulence realization, the rays are traced from one plane to the next where at the $i$ ’th plane, the phase screen changes the $x$ -direction ray angle as

Eq. (1)

{θ^{'}}_{i} = θ_{i} + \frac{1}{k} \frac{d ψ_{i}}{d x},

where

θ_{i}

and

{θ^{'}}_{i}

are the incoming and outgoing ray angles and

\frac{d ψ_{i}}{d x}

is the phase-screen gradient in the

x

-direction. Ray angle changes in the

y

-direction are computed similarly using the

y

-direction gradients. A three-point calculation is used to compute the screen gradients, and the gradients are interpolated at the position where the ray strikes the plane. The

x

- and

y

-direction angles for the rays that eventually fall within the aperture region at the observation plane are averaged to obtain the aperture-averaged AOA. Our approach for tracing rays through the split-step phase screens is similar to the method presented by Lachinova et al.,¹² however, we do not compute the brightness function that is part of their incoherent image simulation.

3. Theory

The variance of the aperture-filtered AOA fluctuations for a circular aperture and for a plane wave propagated through homogeneous and locally isotropic inertial subrange turbulence is given by Cheon and Muschinski²⁰ as

Eq. (2)

{⟨ \bar{θ^{2}} ⟩}_{P} = γ_{P} (q) C_{n}^{2} L D^{- \frac{1}{3}} = 0.81 γ_{P} (q) σ_{R}^{2} k^{- \frac{7}{6}} L^{- \frac{5}{6}} D^{- \frac{1}{3}},

where

Eq. (3)

γ_{P} (q) = \frac{\sqrt{3}}{16} Γ (\frac{1}{6}) Γ (\frac{8}{3}) {(\frac{β}{2})}^{- \frac{1}{3}} {1 + \frac{6}{5} {(\frac{π}{2})}^{\frac{1}{6}} β^{\frac{1}{3}} q^{\frac{1}{3}} {[1 + \frac{π^{2}}{4} β^{4} q^{4}]}^{\frac{5}{12}} \sin [\frac{5}{6} \arctan (\frac{2}{π β^{2} q^{2}})]},

where

D

is the aperture diameter,

q = D / \sqrt{λ L}

is the aperture diameter normalized by the Fresnel length,

β = 0.5216

is a constant, and

Γ

is the gamma function. Equations (2) and (3) were developed assuming the Rytov approximation for the propagating field. This result is applicable for all values of

q

and assumes infinite outer scale and zero inner scale.

4. Results

4.1.

Observation Plane Illustration

Before presenting the AOA results, we first examine the character of the rays and waves arriving at the observation plane for a single turbulence realization. Figure 2 shows single realization examples within the central $0.4 m \times 0.4 m$ area for the scattering scenarios corresponding to $σ_{R}^{2} = 0.0267$ , 0.267, 2.67, and 26.7. Shown for comparison are: (1) ray-density fields obtained from the ray simulation and (2) intensity and phase fields obtained from the wave optics simulation. The ray density fields were obtained by counting the rays within areas defined by the binning values. The binning is relative to the grid spacing $Δ x$ , for example, $10 \times 10$ binning indicates an area of $10 mm \times 10 mm$ . The bin size was selected to emphasize the most apparent characteristic feature sizes within the frame. A contour plot algorithm was also applied to the ray density plot to reduce the blockiness of the display.

Fig. 2

Observation plane examples. Left-to-right in each row: ray-density, wave-intensity, and wave-phase cross-sections. Turblulent scattering strength increasing from top-to-bottom along columns: (a–c) $σ_{R}^{2} = 0.0267$ , $C_{n}^{2} = 10^{- 16} m^{- 2 / 3}$ , $r_{0} = 21.3 cm$ ; (d–f) $σ_{R}^{2} = 0.267$ , $C_{n}^{2} = 10^{- 15} m^{- 2 / 3}$ , $r_{0} = 5.3 cm$ ; (g–i) $σ_{R}^{2} = 2.67$ , $C_{n}^{2} = 10^{- 14} m^{- 2 / 3}$ , $r_{0} = 1.33 cm$ ; (j–l) $σ_{R}^{2} = 26.7$ , $C_{n}^{2} = 10^{- 13} m^{- 2 / 3}$ , $r_{0} = 0.33 cm$ . Ray density binning: (a) $10 \times 10$ bins, (d) $8 \times 8$ bins, (g) $3 \times 3$ bins, and (j) $1 \times 1$ bins.

Common spatial features are apparent in the intensity and ray-density fields as shown in Fig. 2. The corresponding features are related to the effects of refraction as ray tracing does not account for the interference and diffractive effects provided in the wave simulation. The phase cross-sections, displayed modulo- $2 π$ , illustrate the increasing wave front complexity as a function of the Rytov variance. As noted previously, after the ray and wave fields are generated at the observation plane as shown in Fig. 2, a circular aperture of diameter $D$ is imposed and the methods described in Sec. 2 are applied to find the aperture-averaged AOA field for the particular turbulence realization.

4.2.

Variances of Aperture-Filtered Angle-of-Arrival Fluctuations

The variances of aperture-averaged AOA fluctuations are computed from many independent turbulence realizations. Figure 3 shows results for the theoretical [Eq. (2)] and simulated aperture-averaged AOA variances as a function of the normalized aperture diameter $q$ for the four values of ${C_{n}}^{2}$ listed in Table 1. Each data point was estimated from 1000 turbulence realizations. The general behavior is that the AOA variance decreases as a function of increasing aperture diameter. Both simulation results generally agree well with the AOA theory²⁰ for weak scattering (e.g., $σ_{R}^{2} ≲ 0.2$ ) but decrease significantly below the theory curve as the Rytov variance increases. The results also show the agreement between theory and simulation is better in stronger scattering for larger apertures (e.g., $σ_{R}^{2} = 2.67$ and $q ≳ 10$ ). The over-prediction of the AOA variances by the Rytov theory is likely a result of the well-known failure of the Rytov approximation in the strong-moderate- to strong-scattering regime. We note that the Rytov theory is also known to over-predict intensity scintillation as a function of the Rytov variance in the moderate- and strong-scattering regimes.

Fig. 3

Aperture-averaged AOA simulation results and Rytov-based theory as a function of normalized aperture size $q$ : (a) $σ_{R}^{2} = 0.0267$ , $C_{n}^{2} = 10^{- 16} m^{- 2 / 3}$ , (b) $σ_{R}^{2} = 0.267$ , $C_{n}^{2} = 10^{- 15} m^{- 2 / 3}$ , (c) $σ_{R}^{2} = 2.67$ , $C_{n}^{2} = 10^{- 14} m^{- 2 / 3}$ , and (d) $σ_{R}^{2} = 26.7$ , $C_{n}^{2} = 10^{- 13} m^{- 2 / 3}$ .

In general, the ray-based and wave-based simulation variances compare favorably in Fig. 3. The largest discrepancy is a 15% higher variance for the rays in strong scattering with a small aperture ( $σ_{R}^{2} = 26.7$ and $q = 0.63$ ). For the weak-to-moderate scattering cases [Figs. 3(a)–3(c)], the ray and wave results match well for $q ≳ 2$ . That is, diffraction effects are negligible in the weak- and moderate-scattering regime if the aperture diameter is larger than about twice the Fresnel length. This result is consistent with the theoretical analysis by Cheon and Muschinski²⁰ who found that for weak scattering, the AOA variance predicted by geometrical optics deviates from that predicted by the Rytov theory by $< 1 %$ if $q > 1.65$ . For smaller $q$ , the wave simulations predict a lower AOA variance than the ray simulations. For the strongest scattering case [Fig. 3(d)], the wave and ray variance values converge for $q ≳ 10$ .

Figure 4 shows an alternative presentation of the results where aperture-averaged AOA variances for the small aperture ( $q = 0.63$ ; $D = 0.02 m$ ) and the large aperture ( $q = 12.6$ ; $D = 0.4 m$ ) are plotted as a function of the Rytov variance, $σ_{R}^{2}$ . The Rytov theory curves demonstrate the linear dependence described by Eq. (2). For the small aperture [Fig. 4(a)], the simulations show a positive slope as a function of scattering strength that is significantly shallower than the theory line. For large aperture diameters, the Rytov-based theory is expected to collapse to the geometrical optics solution.²⁰ This behavior is suggested in the large aperture case [Fig. 4(b)], where the simulation results follow the theory line much more closely.

Fig. 4

Aperture-averaged AOA simulation results and Rytov-based theory as a function of $σ_{R}^{2}$ for: (a) $q = 0.63$ ; $D = 0.02 m$ and (b) $q = 12.6$ ; $D = 0.4 m$ .

Figure 4 also shows a small difference between the ray-based and wave-based results for the small aperture ( $q = 0.63$ ), but essentially no difference for the large aperture ( $q = 12.6$ ). As noted previously, geometrical optics is known to be questionable for $q < 1$ , which is consistent with the fact that our wave-optics simulations produce smaller AOA variances than the AOA variances produced by our ray-tracing simulations. We attribute this to “Fresnel filtering.” The Fresnel-filter kernel in the wave-number representation of the AOA variance for plane waves predicted by the Rytov theory has the form $K_{F} (x) = 1 + (2 π / x^{2}) \sin (x^{2} / 2 π)$ ,²⁰ where $x = κ {(λ L)}^{1 / 2}$ is the wave number normalized by the Fresnel wave number $1 / {(λ L)}^{1 / 2}$ . The case $x ≪ 1$ (Fresnel length is small compared with $κ^{- 1}$ , the wave-front tilt length scale under consideration) is the geometrical-optics limit where the Fresnel length has no effect on the AOA variance. For $x$ comparable to and $> 1$ , the Fresnel length (and therefore diffraction) does play a role: $K_{F} (x)$ becomes smaller than $K_{F} (x)$ in the geometrical-optics limit. Therefore, diffraction causes the AOA variance predicted by wave optics to be smaller than the AOA variance predicted by geometrical optics.

That is, a finite Fresnel length acts as a characteristic filter scale similar to the filtering effects due to a finite inner turbulence length scale or a finite aperture diameter. Essentially, diffraction causes wave-front tilt components with length scales comparable to and smaller than the Fresnel length scales to be smoothed out. The Fresnel-filtering effect is significant when the aperture size and the inner scale of turbulence are both comparable to or smaller than the Fresnel length, and it is negligible otherwise. The Fresnel filtering is accounted for in the wave-optics simulations but not in the ray-tracing simulations.

To examine the AOA distribution results in detail, we generated histogram estimates of the probability density functions (pdfs). Figures 5(a)–5(c) show the histogram plots for the weakest scattering case and for three aperture diameters. To reduce noise in the histogram, we used 10,000 realizations for these results. We also generated corresponding quantile–quantile (q–q) plots²³ shown in Figs. 5(d)–5(f). The q–q plots compare the simulation sample distribution to a standard normal (Gaussian) distribution. In all cases, the ray and wave samples create straight lines that indicate Gaussian distributions. The two different slopes in the q–q curves in Fig. 5(d) correspond to a difference in the AOA variance for the wave and ray results, which was discussed previously for small $q$ . Histogram and q–q plots for strong scattering are shown in Fig. 6. These distributions are also Gaussian and q–q plot slope differences in Figs. 6(d)–6(f) are consistent with the variance results in Fig. 3(d). We expect a Gaussian distribution for the aperture-averaged AOA if the propagation through turbulence and aperture averaging represents a large sum of independent tilt components such that the central limit theorem applies.

Fig. 5

Statistics of aperture-filtered AOA fluctuations for weak scattering, $σ_{R}^{2} = 0.0267$ . Left-to-right in each row: $q = 0.63$ ( $D = 0.02 m$ ), $q = 6.3$ ( $D = 0.2 m$ ), and $q = 12.6$ ( $D = 0.4 m$ ). (a–c) pdf estimates and (d–f) corresponding q–q plots.

Fig. 6

Statistics of aperture-filtered AOA fluctuations for strong scattering, $σ_{R}^{2} = 26.7$ . Left-to-right in each row: $q = 0.63$ ( $D = 0.02 m$ ), $q = 6.3$ ( $D = 0.2 m$ ), and $q = 12.6$ ( $D = 0.4 m$ ). (a–c) pdf estimates and (d–f) corresponding q–q plots.

5. Conclusions

The simulations illustrate that after propagation through homogeneous turbulence, wave and ray results share obvious intensity/ray-density and AOA characteristics related to the effects of refraction. We find that both the wave and ray simulation results generally predict variances of the aperture-averaged AOA fluctuations that agree well with the Rytov theory for weak scattering (e.g., $σ_{R}^{2} ≲ 0.2$ ), but as scattering strength (i.e., the Rytov variance) increases, the simulation results demonstrate a positive slope that can be significantly shallower than that predicted by the weak-scattering (Rytov) theory. This over-prediction by the Rytov theory is likely a result of the failure of the Rytov approximation in the strong-scattering regime.

Some insights related to the negligibility of diffraction effects are gained from a comparison of the wave and ray simulation results. For weak-to-moderate scattering regimes ( $σ_{R}^{2} ≲ 2.67$ ), the ray and wave results match well for apertures where $q ≳ 2$ . This result is consistent with the theoretical analysis by Cheon and Muschinski²⁰ who found that for weak scattering, the AOA variance predicted by geometrical optics deviates from that predicted by the Rytov theory by $< 1 %$ if $q > 1.65$ . For the strongest scattering case studied here ( $σ_{R}^{2} = 26.7$ ), the wave and ray variance values agree well for $q ≳ 10$ . For smaller apertures, we find the wave simulation variances are smaller than those produced by the ray simulations, and we attribute this to a “Fresnel filtering” effect associated with diffraction. Both the wave and ray AOA distributions are found to be Gaussian.

Acknowledgments

Earlier results of this work appeared in the conference proceedings paper: D. Voelz, E. Wijerathna, X. Xiao, A. Muschinski, “Angle-of-arrival variance of waves and rays in strong atmospheric scattering: split-step simulation results,” Proc. SPIE 10410, Unconventional and Indirect Imaging, Image Reconstruction, and Wavefront Sensing 2017, 104100Y (6 September 2017). This work was funded by Air Force Office of Science Research (AFOSR) (FA9550-17-C-0021 and FA9550-18-C-0011), Air Force Office of Science Research (AFOSR) Multidisciplinary University Research Initiatives (MURI) Program (FA9550-12-1-0449 and FA9550-18-1-0009), National Science Foundation (NSF) Physical and Dynamics Meteorology Program (AGS-1547476), and U.S. Department of Defense (DOD) High Performance Computing Modernization Program (Frontier Project FP-CFD-FY14-007).

References

1.

V. I. Tatarskii, Wave Propagation in a Turbulent Medium, 285 McGraw-Hill, New York (1961). Google Scholar

2.

V. I. Tatarskii, The Effects of the Turbulent Atmosphere on Wave Propagation, 472 Israel Program for Scientific Translation, Jerusalem (1971). Google Scholar

3.

A. D. Wheelon, Electromagnetic Scintillation I—Geometric Optics, 455 Cambridge University Press, Cambridge, United Kingdom (2001). Google Scholar

4.

A. D. Wheelon, Electromagnetic Scintillation II—Weak Scattering, 440 Cambridge University Press, Cambridge, United Kingdom (2003). Google Scholar

5.

L. C. Andrews and R. L. Phillips, Laser Beam Propagation through Random Media, 783 2nd ed.SPIE Press, Bellingham, Washington (2005). Google Scholar

6.

O. Korotkova, Random Light Beams: Theory and Applications, 361 CRC Press, Boca Raton, Florida (2014). Google Scholar

7.

A. Sommerfeld and I. Runge, “Anwendung der Vektorrechnung auf die Grundlagen der geometrischen Optik,” Ann. Phys., 340 277 –298 (1911). https://doi.org/10.1002/andp.19113400705 Google Scholar

8.

D. L. Fried, “Optical resolution through a randomly inhomogeneous medium for very long and very short exposures,” J. Opt. Soc. Am., 56 1372 –1379 (1966). https://doi.org/10.1364/JOSA.56.001372 JOSAAH 0030-3941 Google Scholar

9.

J. W. Strohbehn, “Line-of-sight wave propagation through the turbulent atmosphere,” Proc. IEEE, 56 1301 –1318 (1968). https://doi.org/10.1109/PROC.1968.6572 IEEPAD 0018-9219 Google Scholar

10.

J. H. Churnside and R. J. Lataitis, “Wander of an optical beam in the turbulent atmosphere,” Appl. Opt., 29 926 –930 (1990). https://doi.org/10.1364/AO.29.000926 APOPAI 0003-6935 Google Scholar

11.

H. Yuksel, W. Atia and C. C. Davis, “A geometrical optics approach for modeling atmospheric turbulence,” Proc. SPIE, 5891 589109 (2005). https://doi.org/10.1117/12.615522 PSISDG 0277-786X Google Scholar

12.

S. L. Lachinova et al., “Comparative analysis of numerical simulation techniques for incoherent imaging of extended objects through atmospheric turbulence,” Opt. Eng., 56 071509 (2017). https://doi.org/10.1117/1.OE.56.7.071509 Google Scholar

13.

Y. Wang and W. D. Kulatilaka, “Optical ray tracing method for simulating beam-steering effects during laser diagnostics in turbulent media,” Appl. Opt., 56 E106 –E115 (2017). https://doi.org/10.1364/AO.56.00E106 APOPAI 0003-6935 Google Scholar

14.

B. E. Vyhnale, “Path profiles of ${C_{n}}^{2}$ derived from radiometer temperature measurements and geometrical ray tracing,” Proc. SPIE, 10096 100961G (2017). https://doi.org/10.1117/12.2252278 PSISDG 0277-786X Google Scholar

15.

J. R. Peterson et al., “Simulation of astronomical images from optical survey telescopes using a comprehensive photon Monte Carlo approach,” Astrophys. J. Suppl. Ser., 218 14 (2015). https://doi.org/10.1088/0067-0049/218/1/14 APJSA2 0067-0049 Google Scholar

16.

D. L. Knepp, “Multiple phase-screen calculation of the temporal behavior of stochastic waves,” Proc. IEEE, 71 722 –737 (1983). https://doi.org/10.1109/PROC.1983.12660 IEEPAD 0018-9219 Google Scholar

17.

J. M. Martin and S. M. Flatté, “Intensity images and statistics from numerical simulation of wave propagation in 3-D random media,” Appl. Opt., 27 2111 –2126 (1988). https://doi.org/10.1364/AO.27.002111 APOPAI 0003-6935 Google Scholar

18.

J. D. Schmidt, Numerical Simulation of Optical Wave Propagation, 149 –182 SPIE Press, Bellingham, Washington (2010). Google Scholar

19.

X. Xiao and D. Voelz, “On-axis probability density function and fade behavior of partially coherent beams propagating through turbulence,” Appl. Opt., 48 167 –175 (2009). https://doi.org/10.1364/AO.48.000167 APOPAI 0003-6935 Google Scholar

20.

Y. Cheon and A. Muschinski, “Closed-form approximations for the angle-of-arrival variance of plane and spherical waves propagating through homogeneous and isotropic turbulence,” J. Opt. Soc. Am. A, 24 415 –422 (2007). https://doi.org/10.1364/JOSAA.24.000415 JOAOD6 0740-3232 Google Scholar

21.

D. G. Voelz, Computational Fourier Optics, 63 –96 SPIE Press, Bellingham, Washington (2011). Google Scholar

22.

Y. Cheon et al., “Angle-of-arrival anemometry by means of a large-aperture Schmidt–Cassegrain telescope equipped with a CCD camera,” J. Opt. Soc. Am. A, 24 3478 –3492 (2007). https://doi.org/10.1364/JOSAA.24.003478 JOAOD6 0740-3232 Google Scholar

23.

W. L. Martinez and A. R. Martinez, Computational Statistics Handbook with MATLAB, 124 –132 3rd ed.CRC Press, Boca Raton, Florida (2016). Google Scholar

Biography

David Voelz is a professor in the Klipsch School of Electrical and Computer Engineering at New Mexico State University. His research interests include spectral and polarization sensing, laser beam propagation through atmospheric turbulence, laser communications, imaging theory, and astronomical instrumentation development. He earned his PhD EE degree from the University of Illinois in 1987. He is a fellow of SPIE.

Erandi Wijerathna is a PhD student in the Klipsch School of Electrical and Computer Engineering at New Mexico State University. Her research interests include laser beam propagation through atmospheric turbulence, imaging theory, and polarimetric lidar. She earned her BS (honors) in physics from the University of Colombo, Sri Lanka, in 2011. She received her MS degree in physics and in electrical engineering from New Mexico State University in 2016.

Andreas Muschinski (MS in physics, TU Braunschweig, Germany, 1990; PhD in meteorology, U Hannover, Germany, 1992; habilitation in meteorology, U Hannover, Germany, 1998) is an atmospheric physicist. He was a CIRES research scientist at the University of Colorado at Boulder, Colorado (1998–2004) and a professor of electrical engineering at UMass Amherst (2004–2011). Since 2011, he has been a senior research scientist at NorthWest Research Associates, Boulder, Colorado. During the last 28 years, he has conducted research on atmospheric turbulence and atmospheric wave propagation.

Xifeng Xiao received her BS and MS degrees in physics from Xiamen University Fujian, China, in 1998 and 2001, respectively. She earned her MS and PhD degrees in electrical engineering from New Mexico State University in 2004 and 2008, respectively. Currently, she is a software engineer at Ball Aerospace, Albuquerque, NM. Her research interests include simulation and modeling of free-space laser communication, liquid-crystal polarization, and demonstration and implementation of AOTF technologies.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation Download Citation

David G. Voelz, Erandi Wijerathna, Andreas Muschinski, and Xifeng Xiao "Computer simulations of optical turbulence in the weak- and strong-scattering regime: angle-of-arrival fluctuations obtained from ray optics and wave optics," Optical Engineering 57(10), 104102 (5 October 2018). https://doi.org/10.1117/1.OE.57.10.104102

Received: 23 April 2018; Accepted: 29 August 2018; Published: 5 October 2018

Access the abstract

JOURNAL ARTICLE
9 PAGES

DOWNLOAD PAPER SAVE TO MY LIBRARY

GET CITATION

CITATIONS

Cited by 20 scholarly publications.

Explore citations on Lens.org

KEYWORDS

Geometrical optics

Scattering

Turbulence

Computer simulations

Diffraction

Optical turbulence

Ray tracing

1.

Introduction

2.

Approach

Fig. 1

Table 1

2.1.

Wave Optics Simulation

2.2.

Ray Tracing Simulation

Eq. (1)

3.

Theory

Eq. (2)

Eq. (3)

4.

Results

4.1.

Observation Plane Illustration

Fig. 2

4.2.

Variances of Aperture-Filtered Angle-of-Arrival Fluctuations

Fig. 3

Fig. 4

Fig. 5

Fig. 6

5.

Conclusions

Acknowledgments

References

Biography

Show All Keywords

Keywords/Phrases

Search In:

Publication Years