非均匀湍流路径光传播数值模拟中相位屏间C2n的选取

针对非均匀湍流路径光传播的数值模拟,为减小相位屏间湍流强度均匀的假设所引起的模拟误差,数值分析了在种选择相位屏间折射率结构常数C2n的可能方案.通过计算平面波和球面波的Rytov指数和Fried参数,对比五种方案得到的Rytov指数和Fried参数与理论结果的相对误差.结果表明,当相位屏间C2n选择为相位屏间湍流大气实际的折射率结构常数的路径平均值时,模拟误差最小,该选择方案为相位屏间C2n选取的最伟方案;当根据相位屏间湍流大气的Rytov指数相等的原则来设置相位屏分布时,选择相位屏间C2n为两相邻相位屏处折射率结构常数的儿何平均值时,在一定的误差允许范围内,该选择方案具有可取性.     

非均匀湍流路径上光传播数值模拟的相位屏分布

针对非均匀湍流路径上光传播的数值模拟,研究了相位屏分布的三种不同设置方案,即均匀分布相位屏（UDPS）、等Rytov指数间隔的相位屏（ERPS）和等Fried参数间隔的相位屏（EFPS）.结果表明,ERPS方案根据折射率起伏的强弱设置密疏相间的相位屏,实现了对不同折射率起伏区的充分采样,是非均匀路径上相位屏分布的首选方案.在满足相位屏间距基本要求的前提下,决定相位屏间隔的Rytov指数间隔的选取具有较大的任意性,统计结果显示,尽管相位屏间Rytov指数变化超过一个数量级,但光束有效半径、能量Strehl比的变化均小于2.1%,模拟结果与理论结果符合很好,表明了ERPS方案的可靠性和稳定性. 关键词： 光传播 非均匀湍流路径 数值模拟 相位屏分布     

内尺度对光强闪烁激光雷达测量结果的影响

利用Kolmogorov谱、修正Hill谱和Rytov改进模型三种大气湍流功率谱模型,得到了光强闪烁激光雷达探测路径上闪烁指数与大气折射率结构常数之间的关系。分析了不同内尺度下大气湍流强度的变化情况,并与不考虑内尺度的情况进行了比较。结合实验数据对比分析了内尺度对光强闪烁激光雷达在探测大气湍流时的影响程度,结果表明在内尺度的取值变化范围内,采用修正Hill谱时,理论上有限内尺度的折射率结构常数与不考虑内尺度时折射率结构常数的比值可达9,实验中传输距离为1020 m和2040 m的传输路径上两者最大偏差为0.4和0.1个量级;Rytov改进模型下,理论上有限内尺度的折射率结构常数与不考虑内尺度时折射率结构常数的比值可达6,实验中同样传输路径上两者最大偏差为0.6和0.3个量级。理论和实验结果表明:有限内尺度的折射率结构常数测量结果在一定程度上偏离不考虑内尺度的折射率结构常数,且影响程度与激光传输距离和内尺度的大小有关。因此,在光强闪烁激光雷达的大气湍流探测过程中,必须考虑内尺度效应。     

激光在湍流大气中的干涉仿真

通过研究激光在湍流大气中的传输情况，旨在为光学高分辨率干涉成像提供理论依据。采用相位屏的近似处理方法对激光通过大气湍流传输进行模拟，研究激光通过湍流大气的干涉情况。试验结果表明：孔径间距越大，干涉条纹越密；在传输距离为2km，大气湍流折射率结构常数C2n为10-14，即接近强湍流的情况下，可以形成比较好的干涉条纹。孔径间距较小时，干涉的明暗条纹清晰可见。在实际干涉成像中，可以选择适当的孔径间距来获得较好的干涉条纹，从而获得目标的高分辨率像。     

非Kolmogorov湍流相位屏仿真及对光束传输模拟的影响

采用正态分布模型来描述传播路径上功率谱幂值的非均匀性,利用功率谱反演法构建了基于等效结构常数的非Kolmogorov湍流相位屏,并利用多重相位屏法进行了高斯光束在均匀各向同性湍流与非Kolmogorov湍流两种模型下的传输模拟。通过观察模拟光束的光强均匀度、光束漂移以及闪烁指数等,研究了非Kolmogorov湍流对光束传输模拟的影响。结果发现,当光束通过单个非Kolmogorov相位屏时,光强最大值与光强均匀度随着幂值先增大后减小,光束漂移随幂值在一定范围内单调变化。当光束通过多重相位屏时,发现模拟光束的光强闪烁指数会受到相位屏数量以及湍流模型的影响,当相位屏数量较多时,均匀各向同性湍流模型下模拟得到的光强闪烁指数会大于非K湍流模型下的结果,且非K湍流模拟的光束漂移与均匀各项同性湍流模拟得到的光束漂移的相对误差会随着相位屏数量增多而趋于0。     

测距光束光斑漂移的数值仿真

针对远距离激光测距,如地球同步卫星测距、激光测月,利用基于傅里叶变换的分步传播方法对测距光束在湍流大气中的上行传播进行数值仿真.根据大气结构常数廓线,采用等Rytov指数间隔法分布相位屏,与理论值比对的结果表明相位屏的分布合理.在此基础上,对光束的光斑漂移进行分析,得到光斑漂移量与准直测距光束束腰的定量关系,为测距激光...     

非均匀采样的功率谱反演大气湍流相位屏的快速模拟

对大气湍流功率谱非均匀采样可以有效改善传统功率谱反演法低频采样严重不足的缺陷, 实现高精度的大气湍流相位屏的模拟. 但采用的直接求和运算计算复杂度高, 相位屏的模拟速度极慢. 将非均匀快速傅里叶变换(NUFFT)引入到大气湍流相位屏的模拟, 可以实现相位屏的快速模拟. 从随机过程的谱分解出发, 将大气湍流相位随机过程表示为有限谐波分量叠加和的均方极限. 通过一个高斯核函数的卷积, 将非均匀分布的谐波复振幅映射到均匀网格空间, 进而利用快速傅里叶变换, 降低计算复杂度, 加快大气湍流相位屏的模拟速度. 以大气湍流的Kolmogorov 谱为例, 利用NUFFT仿真得到大气湍流相位屏, 并对相位屏的模拟精度、模拟速度和误差进行统计分析. 结果表明, NUFFT的引入可以实现快速、高精度的大气湍流相位屏的模拟.     

大气湍流对轨道角动量态复用系统通信性能的影响

通过数值仿真,以多个相位屏模拟大气湍流,研究大气湍流对轨道角动量(Orbital Angular Momentum,OAM)态复用通信系统质量的影响.改变大气湍流折射率结构常量C2n和传输距离z,分析其对OAM态复用系统归一化接收功率、误码率等的作用.研究结果表明当大气湍流强度增加和传输距离增长时,复用系统归一化接收功率下降、误码率增加.当C2n增大到5×10-13 m-2/3时,系统归一化接收功率下降到0.2以下,误码率增大到0.3以上;当传输距离增大到10 000m,系统归一化接收功率下降到0.3以下,误码率增大到0.35以上.同时,对比不同方位角指数l取值的两个OAM态复用系统受大气湍流的影响,发现:l取值越接近的系统,受大气湍流影响越大.在等同光学信噪比下大气湍流导致OAM态复用系统误码率性能降低2至3个数量级.该研究对OAM态复用系统的实用化具有积极的参考作用.     

基于随机数据元扩张的大气湍流相位屏数值模拟

准确模拟大气湍流相位屏是建立大气湍流数值模型的核心问题。从大气湍流的统计学特性入手,利用大气湍流波前相位结构函数建立了一种新的大气湍流相位屏数值模拟方法,通过随机数据元扩张对大气湍流波前畸变相位分布的尺度随机性和高频分量随机性进行了模拟,并以此为基础通过相位屏近似构建了满足Kolmogorov统计规律的大气湍流数值模型。计算结果表明,随机数据元扩张法生成的相位屏在统计特性上与理论值基本吻合,在低频部分相对于功率谱反演法有明显的提升。同时,随着相位屏网格数的增加,计算结果的高频特性逐渐呈现并逼近理论值。对于由相位屏构建的大气湍流数值模型,在此通过光强闪烁率作为判据进行了验证,结果表明对于弱湍流和中等强度的湍流,模拟结果与理论计算基本相同;对于强湍流则误差较为明显,相对误差最大可达40%。     

基于拉盖尔-高斯光束的单光子捕获概率研究

单光子源通常采用基于高斯光束的高度衰减激光脉冲,假设激光束具有初始高斯时域脉冲波形和TEM01模拉盖尔-高斯空域分布.基于折射率起伏的Rytov近似和修正von Karman谱模型,研究了大气湍流对星地量子通信单光子捕获概率的影响;建立了上行信道和下行信道的单光子捕获概率理论模型;针对低轨卫星-地面站间激光链路,对单光子捕获概率进行了分析.结果表明:上行信道的单光子捕获概率强烈依赖于地面折射率结构常数C2n(0),且随着C2n(0)的增加而减小;然而,下行信道的单光子捕获概率并不依赖于C2n(0),即大气湍流对其没有影响.     

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.     

A three-dimensional refractive index model for simulation of optical wave propagation in atmospheric turbulence

Numerical modeling of optical wave propagation in atmospheric turbulence is traditionally performed with using the so-called “split”-operator method, when the influence of the propagation medium’s refractive index inhomogeneities is accounted for only within a system of infinitely narrow layers (phase screens) where phase is distorted. Commonly, under certain assumptions, such phase screens are considered as mutually statistically uncorrelated. However, in several important applications including laser target tracking, remote sensing, and atmospheric imaging, accurate optical field propagation modeling assumes upper limitations on interscreen spacing. The latter situation can be observed, for instance, in the presence of large-scale turbulent inhomogeneities or in deep turbulence conditions, where interscreen distances become comparable with turbulence outer scale and, hence, corresponding phase screens cannot be statistically uncorrelated. In this paper, we discuss correlated phase screens. The statistical characteristics of screens are calculated based on a representation of turbulent fluctuations of three-dimensional (3D) refractive index random field as a set of sequentially correlated 3D layers displaced in the wave propagation direction. The statistical characteristics of refractive index fluctuations are described in terms of the von Karman power spectrum density. In the representation of these 3D layers by corresponding phase screens, the geometrical optics approximation is used.     

Simulation of extended phase screens in problems of optical wave propagation through the atmosphere

Turbulent medium in problems of optical wave propagation through the atmosphere is usually nodeled as a set of statistically independent plane screens with a random two-dimensional field of phase progress. In this paper, we develop methods for the formation of nonperiodic phase screens infinitely extended in a certain direction, which are required in problems of dynamic simulation of wave propagation. __________ Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 49, No. 1, pp. 21–34, January 2006.     

不同方法测量折射率结构常数的比较

对几种测量折射率结构常数（Ｃ＾２ｎ）的方法进行了分析和比较,指出这几种方法的测量结果有时不一致的原因。结果表明,和双点法测量相比,温度谱幂率的变化对闪烁法测量结果造成较大偏差,而对到达角起伏测量和声雷达测量结果差别较小;外悄度发迹对闪烁法测量几乎没有影响,但当外尺度变小时,利用双点法测量,到达角起伏测量和声雷达测量结果都有较大的影响。给出了双点测量Ｃ＾２ｎ和外尺度的关系的实验验证,并给出了两固定点     

Propagation of a Gaussian-beam wave through a random phase screen

Tractable analytic expressions are developed for a variety of basic statistical quantities involving a Gaussian-beam wave propagating through a random medium confined to a portion of the propagation path between input and output planes, the limiting case of which defines a thin random phase screen. For a plane wave incident on a phase screen located midway between input and output planes, it is well known that the statistics in the receiver plane are in close agreement with those associated with a plane wave propagating through an extended random medium between input and output planes. For a similar comparison between a phase screen and extended turbulence in the case of a Gaussian-beam wave at the input plane, the present analysis reveals that the phase screen must be positioned between input and output planes differently from the plane-wave case, the position being dependent upon the Fresnel ratio of the Gaussian beam. The analytic results developed in this paper for the thin phase screen model are based on the Kolmogorov power-law spectrum for refractive-index fluctuations and the Rytov approximation. Extension of these results to multiple phase screens is also discussed.     

Diffraction method for recording optical wave fields

We propose a method for recording an optical wave field that is based on diffraction of the optical wave field by two screens that are complementary to each other. We show that the new method ensures the measurement of the phase of the wave field with a high accuracy and spatial resolution. An optical scheme for implementing the diffraction method and an algorithm for reconstructing the phase of the wave field from the measurement data are developed. We perform a computer simulation of the experiment on recording the optical wave field based on the method proposed.     

大气湍流随机相屏模拟方法研究

对几种模拟随机湍流相位屏的方法进行了总结,利用结构函数对每种方法进行了分析,并从Zernike多项式分解的思路和激光大气传输的角度,通过焦斑扩展和焦斑质心漂移对几种方法进行了完整比较,最后对湍流低频倾斜信息的补偿方式进行了探讨。通过分析发现:傅里叶变换(FFT)法模拟的相屏体现的湍流高频信息较充分,但低频信息不足;Zernike法模拟的相屏体现的湍流低频信息较充分,但高频信息不足。随着采用的Zernike多项式阶数的增加,其体现的高频信息会得到改善,但这些高频信息主要集中在圆域的边缘区域。分形法模拟的相屏体现的湍流高低频信息均较好。针对焦斑扩展和焦斑质心漂移,除分形法外,其他算法均存在不同缺陷。采用FFT法与Zernike多项式展开法相结合的方式或FFT法与湍流自相似理论相结合的方式模拟大气湍流随机相屏,也不失为一种较有效的方式,但从激光大气传输的角度来说,采用分形法模拟湍流随机相屏最好。     

Multicritical behaviors and an induced twist grain boundary phase in a binary liquid crystalline mixture

We report here the composition-temperature phase diagram of a mixture between n=9 and n=16 terms of a chiral tolans series shortly noted nF(2)BTFO1M7. This diagram constitutes an experimental illustration of one of the three theoretical phase diagrams predicted by the Renn and Lubensky's model. The pure compounds n=9 and n=16 exhibit, respectively, the phase sequences Cry-Sm-C*-Sm-A-TGB(A)-N*-BP-I and Cry-Sm-C*-N*-BP-I, where TGB refers to a twist grain boundary phase. Phases identification and transition temperatures, at atmospheric pressure, have been determined by both optical microscopy and photothermal methods. The experimental phase diagram shows the disappearance of Sm-A and TGB(A) mesophases and the appearance of a TGB(C) phase when the n=16 composition increases. Four of the five multicritical points theoretically anticipated by Renn and Lubensky are pointed out.