Generation of infinitely long phase screens for modeling of optical wave propagation in atmospheric turbulence

Numerical modeling of optical wave propagation in atmospheric turbulence is traditionally performed by using the so-called 'split'-operator method, where the influence of the propagation medium's refractive index inhomogeneities is accounted for only within a set of infinitely narrow phase distorting layers (phase screens). These phase screens are generated on a numerical grid of finite size, which corresponds to a rather narrow slice (spatial area) of atmospheric turbulence. In several important applications including laser target tracking, remote sensing, adaptive optics, and atmospheric imaging, optical system performance depends on atmospheric turbulence within an extended area that significantly exceeds the area associated with the numerical grid.

In this paper we discuss methods that allow the generation of a family of long (including infinitely long) phase screens representing an extended (in one direction) area of atmospheric turbulence-induced phase distortions. This technique also allows the generation of long phase screens with spatially inhomogeneous statistical characteristics.     

Study on beam propagation through a double-adaptive-optics optical system in turbulent atmosphere

Beam propagation through a double-adaptive-optics optical system in turbulent atmosphere has been investigated. A typical model of optical systems with the double-adaptive-optics configuration is established, principles of the two adaptive optics and theory of beam propagation through the optical system in turbulent atmosphere are analyzed. Power efficiency and beam quality (BPF) of the received beam are introduced to evaluate performance of the optical system. Under the H-V 5/7 turbulent model, influences of the two adaptive optics installations are numerically calculated, respectively. Results show that the adaptive optics at the launcher can improve both power efficiency and beam quality of the received beam, the adaptive optics at the receiver can significantly improve beam quality of the received beam. Evolution of the system performance with the increase of adaptive optics correction order is numerically simulated, and results show that the double-adaptive-optics configuration has great advantages in improving beam quality of the received beam in optical systems with beam long-distance propagation through the turbulent atmosphere.     

The propagation of radially polarized partially coherent beam through an optical system in turbulent atmosphere

The propagation of a radially polarized partially coherent beam through an optical system in turbulent atmosphere is investigated. Analytical propagation formulae for the radially polarized partially coherent beam through an optical system in turbulent atmosphere are derived, and some analyses are illustrated by numerical examples related to propagation properties of the radially polarized partially coherent beam through a double-lens system in turbulent atmosphere. It is shown that the optical system, structure constant of turbulent atmosphere and the initial coherent length have significant influence on intensity, degree of polarization and degree of coherence for radially polarized partially coherent beam.     

次声波在非均匀大气中的超视距传播特性研究

本文基于Nonlinear Progressive Equation (NPE方程)开展了对非均匀大气中次声波超视距传播特性的研究, 通过数值模拟实验对武汉上空四季次声波传播情况以及路径传输损耗进行了模拟, 获取了次声波在非均匀大气中的超视距传播特性. 计算结果表明, 非均匀大气的性质及其中存在的风对次声波传播有明显的影响, 平流层折射与风速和声波传播方向密切相关, 数值模拟结果表明, 当高斯声源主频为0.1 Hz时, 在不同的背景风场传播条件下, 存在着两个反射高度, 其中40 km 高度反射传输损耗约为25 dB, 110 km反射传输损耗约为50 dB. 关键词： 次声波 超视距传播 非均匀大气 平流层折射     

次声波在非均匀运动大气中非线性传播特性的研究

基于有限时域差分方法将大气中近似到二阶的非线性波动方程进行离散化,得到了数值模拟所采用的差分方程. 在此基础上,对线阵列辐射的脉冲声波在非均匀运动大气中的垂直和斜向传播进行了二维数值模拟,模拟了武汉地区(114:20°E, 30:37°N)在夏季和冬季UT=29000 s时开始传播的脉冲声波在不同时刻的声压分布. 模拟时通过采用Msise00和HWM93 两个大气模型,考虑了由于大气温度和密度变化以及大气风场存在所引起的大气不均匀性和运动性. 通过研究上述两季有风与无风条件下的声压差值p_r,可以发现:风场对次声波在传播中声压分布的影响较大;由于不同季节和不同传播距离上"有效声速"的不同,导致了两季p_r分布波形存在差异;风场对声波非线性传播的影响要远大于其对线性传播的影响. 关键词： 次声波传播 非均匀运动大气 有效声速     

高能激光大气传输的仿真实验研究

 对高能激光在大气中传输时产生的热晕效应及其相位补偿进行了室内仿真实验测量,获得了自适应光学系统开闭环时用来衡量光束传输效果参量（如远场的光斑图像、质心漂移、二阶矩半径、及Strehl比等）。结果表明,Bradley-Hermann热畸变参数在300范围内,自适应光学系统相位补偿效果显著,Strehl比基本大于0.4。然而随着热畸变效应的进一步加强,出现相位补偿不稳定性现象,补偿效果变差。另外,对自适应光学系统开闭环时的实验结果分别与薄透镜近似分析理论及几何光学近似理论作了比较。比较结果表明:薄透镜近似理论与未补偿时的实验结果吻合较好,而几何光学近似理论与补偿后的实验结果存在较大的偏差。     

Double-distance propagation of Gaussian beams passing through a tilted cat-eye optical lens in a turbulent atmosphere

By using the extended Huygens-Fresnel diffraction integral and the method of expanding the aperture function into a finite sum of complex Gaussian functions, an approximate analytical formula of the double-distance propagation for Gaussian beam passing through a tilted cat-eye optical lens and going back along the entrance way in a turbulent atmosphere has been derived. Through numerical calculation, the effects of incidence angle, propagation distance, and structure constant on the propagation properties of a Gaussian beam in a turbulent atmosphere are studied. It is found that the incidence angle creates an unsymmetrical average intensity distribution pattern, while the propagation distance and the structure constant can each create a smooth and symmetrical average intensity distribution pattern. The average intensity peak gradually deviates from the centre, and the central average intensity value decreases quickly with the increase in incidence angle, while a larger structure constant can bring the average intensity peak back to the centre.     

Dynamics of ray propagation beyond a series of space-periodic phase screens

The dynamics of ray propagation through a regularly inhomogeneous medium consisting of a series of identical, deep, sinusoidal phase screens is examined. Discrete mapping is used to analyze the propagation process. A characteristic parameter of the problem is found that determines the transition from regular (nonlinear oscillator and ray capture) to chaotic regimes of ray behavior. Diffuse (stochastic) ray propagation in the given medium is studied—by analogy with a random walk—by a method of fractal analysis based on finding the scaling characteristics of the given process (associated with the fractal dimension).Scientific-Research Radiophysics Institute. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 36, No. 8, pp. 825–831, August, 1993.     

Power coupling of Gaussian vortex beam propagation through an optical system with the Cassegrain-telescope receiver in turbulent atmosphere

Model of Gaussian vortex beam propagation through an optical system with the Cassegrain-telescope receiver in turbulent atmosphere is established. With this model, the analytical formulas of the average intensity distribution at the receiver plane are derived, and the influences of the optical topological charge, the propagation distance and the turbulence strength are numerically analyzed. These studies show that optical power at the receiver plane concentrates in an annular area, which is suitable for power coupling by the Cassegrain-telescope receiver; the optical topological charge of the vortex source need to be optimized to access the most power coupling. Under the H-V 5/7 turbulence model, power coupling efficiencies of the optical system with different parameters are calculated. Results show that in comparison with the Gaussian beams, Gaussian vortex beams have great advantages in power coupling of optical systems with the Cassegrain-telescope receivers in turbulent atmosphere, which can be a new attractive application of the vortex beams.     

Scintillation index of optical wave propagating in turbulent atmosphere

A concise expression of the scintillation index is proposed for a plane optical wave and a spherical optical wave both propagating in a turbulent atmosphere with a zero inner scale and a finite inner scale under an arbitrary fluctuation condition. The expression is based on both the results in the Rytov approximation under a weak fluctuation condition and the numerical results in a strong fluctuation regime. The maximum value of the scintillation index and its corresponding Rytov index are evaluated. These quantities are affected by the ratio of the turbulence inner scale to the Fresnel size.     

Spherical wave propagation through inhomogeneous, anisotropic turbulence: log-amplitude and phase correlations

Inhomogeneity and anisotropy are intrinsic characteristics of daytime and nighttime turbulence in the atmospheric boundary layer. In the present paper, line-of-sight sound propagation through inhomogeneous, anisotropic turbulence with temperature and velocity fluctuations is considered. Starting from a parabolic equation and using the Markov approximation, formulas are derived for the correlation functions and variances of log-amplitude and phase fluctuations of a spherical sound wave. These statistical moments of a sound field are important for many practical applications in atmospheric acoustics. The derived formulas for the correlation functions and variances generalize those already known in the literature for two limiting cases: (a) homogeneous, isotropic turbulence, and (b) inhomogeneous, anisotropic turbulence with temperature fluctuations only. Furthermore, the formulas differ from those for the case of plane wave propagation. Using the derived formulas and Mann's spectral tensor of velocity fluctuations for shear-driven turbulence, the correlation functions and variances of log-amplitude and phase fluctuations are studied numerically. The results obtained clearly show that turbulence inhomogeneity and anisotropy significantly affect sound propagation in the atmosphere.     

Gaussian wave formalism model for propagation of charged-particle beam through a first-order optical systems

An elliptical Gaussian wave formalism model of a charged-particle beam is proposed by analogy with an elliptical Gaussian light beam. In the paraxial approximation, the charged-particle beam can be described as a whole by a complex radius of curvature in the real space domains. Therefore, the propagation and transform of charged-particle beam passing through a first-order optical system is represented by the \textit{ABCD}-like law.As an example of the application of this model, the relation between the beam waist and the minimum beam spot at a fixed target is discussed. The result, well matches that from conventional phase space model, and proves that the Gaussian wave formalism model is highly effective and reasonable.     

Time-dependent propagation of high energy laser beams through the atmosphere

The computation of time-dependent three-space-dimensional laser beam propagation is described. The methods are applicable to the propagation of high energy laser beams through the atmosphere in the presence of a horizontal wind and turbulence for most situations of interest. Possible cases are propagation of cw beams through stagnation zones, multi-pulse propagation, including the self-consistent treatment of pulse self-blooming, and propagation involving transonic slewing. The solution of the Maxwell wave equation in Fresnel approximation is obtained by means of a discrete Fourier transform method, which, surprisingly, gives excellent results for diffraction problems. The latter provide a stringent test for the accuracy of any solution method. Considerable use is also made of discrete Fourier transform methods in solving the hydrodynamic equations. The treatment of turbulence is based on the generation of random phase screens at each calculation step along the propagation path. In a time-dependent calculation the random phase screens can be either made to move with the wind at a given propagation position or generated anew for each successive time. This work was done under contract to the Advanced Research Projects Agency of the Department of Defense, the Army Missile Command, Huntsville, Alabama, and the U.S. Energy Research and Development Administration.     

湍流折射率谱型对大气闪烁和相位起伏功率谱的影响

基于Taylor湍流冻结假设理论,在不同湍流折射率谱型条件下,推导得出了光波闪烁和相位起伏频谱的表达式;数值计算了湍流谱型中折射率标度指数、内尺度以及外尺度变化时对光波频谱的影响。结果表明:随着折射率起伏标度指数的增大,闪烁频谱的低频段不再仅为常数,高频段下降的幂率逐渐增大,同时相位频谱在整个起伏频率段下降的幂率越来越大;湍流内尺度的增加将引起光波频谱的高频段下降的幂率越来越大;而随外尺度的减小,闪烁频谱低频段的振幅减小,这种影响在大口径接收时较为明显,相位谱的低频段幂率减小。     

Changes in generalized Stokes parameters of stochastic electromagnetic beams on propagation through ABCD optical systems and in the turbulent atmosphere

Using the derived formulas for the transformation of 2 × 2 cross-spectral density matrix of the stochastic electromagnetic beams propagating through ABCD optical systems and in the turbulent atmosphere, the changes in the generalized Stokes parameters of the beams propagating under these conditions can be investigated directly. Some typical numerical calculations are illustrated relating to the electromagnetic Gaussian Schell-model beams passing through free space, focal system, dual-focus system, and the turbulent atmosphere with different structure parameters. Further extensions are also pointed out.     

湍流折射率谱型对大气闪烁和相位起伏功率谱的影响

基于Taylor湍流冻结假设理论,在不同湍流折射率谱型条件下,推导得出了光波闪烁和相位起伏频谱的表达式;数值计算了湍流谱型中折射率标度指数、内尺度以及外尺度变化时对光波频谱的影响。结果表明：随着折射率起伏标度指数的增大,闪烁频谱的低频段不再仅为常数,高频段下降的幂率逐渐增大,同时相位频谱在整个起伏频率段下降的幂率越来越大;湍流内尺度的增加将引起光波频谱的高频段下降的幂率越来越大;而随外尺度的减小,闪烁频谱低频段的振幅减小,这种影响在大口径接收时较为明显,相位谱的低频段幂率减小。     

Computational simulation of wave propagation problems in infinite domains

This paper deals with the computational simulation of both scalar wave and vector wave propagation problems in infinite domains. Due to its advantages in simulating complicated geometry and complex material properties, the finite element method is used to simulate the near field of a wave propagation problem involving an infinite domain. To avoid wave reflection and refraction at the common boundary between the near field and the far field of an infinite domain, we have to use some special treatments to this boundary. For a wave radiation problem, a wave absorbing boundary can be applied to the common boundary between the near field and the far field of an infinite domain, while for a wave scattering problem, the dynamic infinite element can be used to propagate the incident wave from the near field to the far field of the infinite domain. For the sake of illustrating how these two different approaches are used to simulate the effect of the far field, a mathematical expression for a wave absorbing boundary of high-order accuracy is derived from a two-dimensional scalar wave radiation problem in an infinite domain, while the detailed mathematical formulation of the dynamic infinite element is derived from a two-dimensional vector wave scattering problem in an infinite domain. Finally, the coupled method of finite elements and dynamic infinite elements is used to investigate the effects of topographical conditions on the free field motion along the surface of a canyon.     

The mean-field method in the problem of acoustic wave propagation in a turbulent atmosphere

We analyze the problem of acoustic-wave propagation in a turbulent atmosphere using the mean-field method. The equation for the sound pressure is written with accuracy up to terms that are linear with respect to the Mach number of the turbulent air flow. An expression for the attenuation constant of the mean field is obtained. For the turbulence model described by the von Kármán correlation function of fluctuations, the attenuation coefficient of the mean field is numerically studied in detail. It is shown that under typical conditions of the near-ground atmospheric layer, the predominant contribution to scattering of acoustic waves is given by the turbulent motion of the air masses. __________ Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 51, No. 5, pp. 413–424, May 2008.