删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

湍流等离子体鞘套中高斯光束的传播特性分析

本站小编 Free考研考试/2021-12-29

摘要:为了研究高斯光束在湍流等离子体鞘套中的传输特性, 根据广义惠更斯-菲涅耳原理, 采用基于快速傅里叶变换的功率谱反演法, 用多随机相位屏来模拟湍流带来的影响. 根据超声速飞行器绕流等离子体流场厚度在厘米级别的特点, 光束在两个相位屏之间的传输过程中采用菲涅耳衍射积分的两次快速傅里叶变换算法(double fast Fourier transform algorithm), 利用多随机相位屏模拟等离子体鞘套湍流对光束传输产生的影响, 解决了多随机相位屏模拟湍流研究中的超短距离传输问题. 当飞行高度为45 km, 飞行速度为18马赫时, 通过对超声速飞行器绕流等离子体流场的统计分析, 发现在此飞行条件下折射率起伏方差的强度范围10–11—10–14. 对高斯光束在湍流等离子体流场中的传输特性进行了数值仿真. 结果表明: 在等离子体鞘套湍流中折射率起伏强度、波长、传输距离等都是影响高斯光束质量的重要因素. 折射率方差越大, 传输距离越长, 光斑弥散越严重, 光强起伏越大, 光强减弱也越明显. 光束的波长越长, 高斯光束抑制湍流的能力越强, 光斑弥散程度越小, 光强起伏也越小.
关键词: 光传输/
随机介质/
随机相位屏

English Abstract


--> --> -->
当飞行器以高超声速在大气中飞行时, 其动能会迅速转化为大量的热能导致飞行器周围的空气产生离解和电离, 最终形成等离子体包覆的流场[1-4]. 等离子体的存在对于电磁波信号有很大的衰减作用, 甚至导致信号中断形成“黑障效应”[5,6]. 目前已经提出几种有希望解决这一难题的几种方法[7-11], 例如通过提高电磁波频率来减小信号衰减, 或者在鞘套层中加入亲电材料, 这两种方法的目的都是为了提高电磁波频率与等离子体频率的比值. 同时有研究表明等离子体鞘套层是不稳定的[12,13], 所以湍流也是影响“黑障效应”的重要因素. 近些年, 国内外也早对湍流等离子体流场进行了较多的研究, 例如Shi等[14]为了研究等离子体鞘套在无线信道上的动态效应, 分析飞行条件变化(例如攻角、小尺度湍流)对介质传输性能的影响, 提出了一种新的自适应马尔可夫信道建模方法, 利用该方法建立的信道的一阶统计特性与无线电波的可变衰减具有良好的一致性, 而且该方法具有良好的自适应能力. He等[15]从通信角度采用透射矩阵法研究了再入等离子体鞘套对电磁波传播的动态影响, 所建立的模型中包含空间湍流和等离子体电子密度的变化, 并对由等离子体的非线性特性引起的幅度和相位变化进行了数值仿真, 发现幅度变化服从对数正态分布, 相位变化服从瑞利分布. Yao等[16]针对等离子体流场的不稳定性, 首先在时-空-频域揭示了电子密度的波动的规律, 并基于这些规律提出了时变电子密度分层模型, 利用Monte Carlo准静态电磁波数值方法计算了动态等离子鞘套层对Ka波段电磁波的影响. Lin和Sproul[17]重点研究了电子密度波动对电磁波传播的影响, 利用Navier-Stokes模型、三维抛物化Navier-Stokes模型、耦合边界层和非黏性流动代码估算了截头锥体的电子密度的均值和方差. Li等[18]基于在高超声速湍流实验中测量的分维数, 建立了分形模型和等离子鞘套湍流功率谱, 为研究等离子体鞘套湍流中的光束传输特性提供了重要参考.
在采用随机相位屏研究湍流问题的文献中, 两个随机相位屏之间的传输过程大多采用单次快速傅里叶变换(single fast Fourier transform, S-FFT)算法[19-21], 但是这种菲涅耳衍射积分的S-FFT算法适用于长距离传输, 例如大气湍流[22-24]和海洋湍流[25-27]. 飞行器周围等离子体包覆的流场的厚度很小, 在厘米或毫米级别[28,29]. 对于厘米和毫米级别的短距离传输问题, S-FFT算法会出现欠采样的问题, 最终将会导致信号失真. 针对等离子体鞘套厚度较小这一特点, 采用菲涅耳衍射积分的两次快速傅里叶变换(double fast Fourier transform, D-FFT)算法研究光束在两个相位屏之间的传输过程. 本文第2部分首先在高度为45 km、飞行速度为18马赫和0攻角的飞行条件, 数值模拟了钝头体飞行目标的高超音速等离子体鞘套湍流流场. 然后利用Gladstone-Dale公式根据等离子体密度计算折射率, 并对折射率的方差进行了估算; 第3部分对比研究了菲涅耳衍射积分的D-FFT算法和S-FFT算法, 说明D-FFT算法在厘米级别的短距离传输问题中的显著优势, 并利用多随机相位屏模拟了高斯光束在等离子体鞘套湍流中的传输特性.
2
2.1.钝头体飞行目标绕流流场数值模拟
-->利用三维Navier-Stokes (N-S)方程和湍流两方程$k $-$\varepsilon $模型数值模拟湍流流场. 在笛卡尔坐标系中, 考虑1989年提出的双温热力学模型[30]和7组分空气化学反应模型(七组分化学反应模型如表1所列), 并对三维N-S方程进行适当封闭, 相应模型如下[31-33]:
No.Reaction
1${\rm {N_2} + M \rightleftharpoons N + N + M}$
2${\rm {O_2} + M \rightleftharpoons O + O + M}$
3${\rm NO + M \rightleftharpoons N + O + M}$
4${\rm O + {N_2} \rightleftharpoons NO + N}$
5${\rm NO + O \rightleftharpoons {O_2} + N}$
6${\rm N + O \rightleftharpoons N{O^ + } + e}$


表1化学反应模型中的方程
Table1.Reactions considered in chemistry model

$\dfrac{{\partial {{Q}}}}{{\partial t}} + \dfrac{{\partial {{E}}}}{{\partial x}} + \dfrac{{\partial {{F}}}}{{\partial y}} + \dfrac{{\partial {{G}}}}{{\partial z}} = \dfrac{{\partial {{E}}_v }}{{\partial x}} + \dfrac{{\partial {{F}}_v }}{{\partial y}} + \dfrac{{\partial {{G}}_v }}{{\partial z}} + {{S}}, $
${{Q}} = {[ E_{{\rm{int}}},\rho _i,\rho u,\rho v,\rho w,E _t ]}^{\rm{T}}, $
${{E}} = {\left[ {E _{{\rm{int}}} u,\rho _i u,\rho u ^2 + p,\rho uv,\rho uw,\left( {E _t + p} \right)u} \right]}^{\rm{T}}, $

$\begin{split}& {{E}}_v \!=\! \bigg[\sum\limits_{i = 1}^{7}\!{e _{v,i} \rho _i U _{\rm{a}} \!+\! \kappa _{{\rm{int}}} \dfrac{{\partial T _{{\rm{int}}} }}{{\partial x}}}, \rho _i U _{\rm{a}}, \tau _{xx} \!-\! \dfrac{2}{3}\rho k,\tau _{xy}, \tau _{xz}, u\tau _{xx} \!+\! v\tau _{xy} \!+\! w\tau _{xz} \!+\! (\kappa \!+\! \kappa _t )\dfrac{{\partial T _{{\rm{int}}} }}{{\partial x}} \!+\! \kappa _{{\rm{int}}} \dfrac{{\partial T _{{\rm{in}}t} }}{{\partial x}} \!+\!\! \sum\limits_{i = 1}^{7} {h _i \rho _i } U _{\rm{a}} \bigg]^{\rm{T}}, \end{split}$

${{S}} = {[ {\dot \omega }_{{\rm{int}}}, {\dot \omega }_i,0,0,0,0]}^{\rm{T}}, $
其中, QS分别是守恒变量和源项矢量; 以$x$方向为例, E${{E}}_v $分别是对流通量矢量和黏性通量矢量; $\kappa $$\kappa_ t$分别是热传导系数和湍流热传导系数; $u$, $v$$w$分别是3个方向的瞬时速度分量; $k$是湍动能; ${\dot \omega }_{{\rm{int}}} $是振动源项, ${\dot \omega }_i $是每个组分所占的质量分数; $\tau _{ij}\; (ij = x,y,z)$为各个方向的剪应力张量, 需要利用湍流模型进行封闭[34]; $e _{v,i} $,$h _i $$\rho _i $分别是每一种组分的振动能量、显焓和密度, 下标$i$代表每一种组分$(i = 1,\;2,\;3,\cdots, 7)$; $\rho $, $p$$U _{\rm{a}} $分别为混合物的密度、压强和平均速度; $E _{{\rm{int}}} $, $E _t $$T _{{\rm{int}}} $分别是混合物的内部能量、湍流内部能量和内部温度.
基于有限体积法, 利用AUSMPW+ (the advection upstream splitting method by pressure-based weight function)方法[35]求解N-S方程, 其中时间推进采用LU-SGS (lower-upper symmetric Gauss-Seidel relaxation)方法[36]. 值得注意的是, 文中采用的湍流$k $-$\varepsilon $模型是利用不同的方程来控制湍流能量传输并适合高雷诺数环境的两方程模型, 其中控制湍动能$k$和湍动能耗散率$\varepsilon $的两方程分别为
$\dfrac{\partial }{{\partial t}}(\rho k) + \dfrac{\partial }{{\partial x _j }}(\rho u _j k) = P - \rho \varepsilon + \dfrac{\partial }{{\partial x _j }}\left[\left({\mu} + \dfrac{{{\mu}_t }}{{\sigma _k }}\right)\dfrac{{\partial k}}{{\partial x _j }}\right], $
$\begin{split} & \dfrac{\partial }{{\partial t}}(\rho \varepsilon ) + \dfrac{\partial }{{\partial x _j }}(\rho u _j \varepsilon ) \\ =\; & C _{\varepsilon _1 } \dfrac{{P\varepsilon }}{k} - C _{\varepsilon _2 } \dfrac{{P\varepsilon ^2 }}{k} + \dfrac{\partial }{{\partial x _j }}\left[ {\left( {{\mu} + \dfrac{{{\mu}_t }}{{\sigma _\varepsilon }}} \right)\dfrac{{\partial \varepsilon }}{{\partial x _j }}} \right], \end{split}$
其中, ${\mu} $, ${\mu}_t $, ${\mu}_i $分别是等离子体的动力学黏性系数、湍流黏性系数和质量平均速度; 常数$C _{\varepsilon _1 } = $$ 1.44,\;C _{\varepsilon _2 } = 1.92,\;\sigma _k = 1.0,\;\sigma _\varepsilon = 1.3$; $P$ 被定义为 $P = {\mu}_t \left( {\dfrac{{\partial u _i }}{{\partial x _j }} + \dfrac{{\partial u _j }}{{\partial x _i }} - \dfrac{2}{3}\dfrac{{\partial u _m }}{{\partial x _m }}\delta _{ij} } \right)\dfrac{{\partial u _i }}{{\partial x _j }} - \dfrac{2}{3}k\dfrac{{\partial u _m }}{{\partial x _m }},$ $u _i $, $u _j $$u _m $为张量表示法中时均速度, $\delta _{ij} $是克罗内克尔符(当$i = j$时, $\delta _{ij} = 1$, 否则$\delta _{ij} = 0$).
2
2.2.对流场数据做统计分析, 估算等离子体折射率起伏方差
-->基于Gladstone-Dale公式[37], 根据2.1节中计算的等离子体密度求解其折射率, Gladstone-Dale公式为
$n\left( {x,y,z} \right) = 1 + G\left( \lambda \right)\rho \left( {x,y,z} \right), $
其中, $\lambda $为光波波长, $G(\lambda )$为Gladstone-Dale常数, 被定义为
$G(\lambda ) = 0.223 \times {10}^{ - 3} \Big( {1 + \dfrac{{7.52 \times {10}^{ - 15} }}{{\lambda ^2 }}} \Big)~\frac{\rm m^2}{ \rm kg}.$
图1(a)图1(b)分别为飞行目标的区域划分图和五个区域的折射率统计方差$\left\langle {n_1^2 } \right\rangle $图. 图2(a)图2(b)分别为飞行目标头部和侧部选取的一条传输路径, 头部路径三维坐标范围分别为: x = [0.0047 m, 0.0170 m], y = [–0.0581 m, –0.0475 m], z = [0.0627 m, 0.0512 m], 侧部路径三维坐标范围分别为: x = [0.0674 m, 0.0788 m], y = [–0.1555 m,–0.1225 m], z = [0.0628 m, 0.0527 m]. 对两条传输路径上的折射率进行方差统计, 统计方法为: 以传输路径中的每个三维流场坐标点为中心, 以超声速湍流中典型的相关长度[38]${L / {10}}$为半径($L$为等离子体鞘套的厚度)建立球模型, 对球内所有数据进行统计, 计算每个网格点处的折射率方差, 流场边界值采用外推法求解. 通过对比图2(a)图2(b)发现, 在接近飞行目标壁面时, 头部折射率方差较大, 比两侧的折射率方差要高出一个数量级, 这是因为等离子鞘套中的湍流状态还会受到激波、高温、高压等因素的影响, 在靠近头部的区域内温度、压强较高且起伏较大. 由图1图2的计算结果显示, 在高度为45 km、飞行速度为18马赫和0攻角的飞行条件下, 飞行目标的绕流流场中折射率起伏方差的强度的为10–11—10–14.
图 1 (a)五个区域; (b)五个区域的折射率方差$\left\langle { n_1^2 } \right\rangle $
Figure1. (a) Five regions; (b) refractive index variance $\left\langle {n_1^2 } \right\rangle $ in five regions.

图 2 折射率方差$\left\langle {n_1^2 } \right\rangle $ (a)头部区域1; (b)侧部区域2
Figure2. Refractive index variance $\left\langle {n_1^2 } \right\rangle $ in different regions: (a) Head area 1; (b) side area 2.

2
3.1.相邻两屏之间传输算法讨论及仿真实验验证
-->S-FFT算法和D-FFT算法都是从菲涅耳衍射公式出发, 但对于衍射距离$d$的要求不同. S-FFT算法主要适用于衍射距离较大的情况, 而D-FFT算法可适用于衍射距离较小的问题, 下面对两种算法进行理论分析和数值仿真, 说明选用D-FFT算法的原因.
3
3.1.1.菲涅耳衍射积分的S-FFT算法
-->在傍轴近似情况下, 菲涅耳衍射积分可以表示为
$\begin{split}& \tilde U\left( {x,y} \right) \\=\; &\dfrac{{\exp \left( {{\rm{j}}{k_\lambda }d} \right)}}{{{\rm{j}}\lambda d}}\int \!\!\!{\int_{ - \infty }^{ + \infty } {{{\tilde U}_0}\left( {{x_0},{y_0}} \right)} } \\&\times\exp \left\{ {\dfrac{{{\rm{j}}{k_\lambda }}}{{2d}}\left[ {{{\left( {x - {x_0}} \right)}^2} + {{\left( {y - {y_0}} \right)}^2}} \right]} \right\}{\rm{d}}{x_0}{\rm{d}}{y_0}.\end{split} $
其中, $\lambda $为光波波长, $k _\lambda = {{2{\text{π}}} / \lambda }$为波数, $d$是观察面到衍射屏的距离, ${\tilde U}_0 \left( {{x_0},{y_0}} \right)$$\tilde U\left( {x,y} \right)$分别是光场在衍射面上和观察面上某一点的复振幅分布. 对(10)式进行适当变换, 并根据傅里叶变换法则, 上式可改写为
$\begin{split}& \tilde U\left( {x,y} \right)\\ =\; & \dfrac{{\exp \left( {{\rm{j}}{k_\lambda }d} \right)}}{{{\rm{j}}\lambda d}}\exp \left[ {\dfrac{{{\rm{j}}{k_\lambda }}}{{2d}}\left( {{x^2} + {y^2}} \right)} \right]\\& \times {\rm{FFT}}\left\{ {{{\tilde U}_0}\left( {{x_0},{y_0}} \right)} \right.\left. {\exp \left[ {\dfrac{{{\rm{j}}{k_\lambda }}}{{2d}}\left( {x_0^2 + y_0^2} \right)} \right]} \right\}.\end{split}$
由于这种算法只需要进行一次傅里叶变换, 所以称为S-FFT算法. 由(11)式可知其空间频率为$ f_{\mu} = \dfrac{x}{{\lambda d}}, f_v = \dfrac{y}{{\lambda d}}$, 即观察屏上不同的点$\left( {x,y} \right)$其实对应不同的空间频率(方向)的光. 又因为光传播方向与空间频率的关系为$f _{\mu} = {1 / {X = {{\cos \alpha } / \lambda }}}$, 所以$f _{\mu} = {{\sin \left( {{\text{π}}/2 - \alpha } \right)} / {\lambda \approx {1 / {\lambda \cdot {{(x} / d}}}}})$.
为了确定观察面的尺寸, 以$x$方向为例进行讨论说明, $y$方向的情况与之类似. 设衍射面在$x$方向上的宽度为$L _{x _0 } $, 共有${N_x}$个抽样点, 则按照抽样定理, 在用FFT变换计算时得到的空间频率最高为${f}_{{{\mu} _{{\rm{max}}}}} = {{{N_x}} / {2{L_{x _0 }}}}$, 空间频率为$f _{{{\mu} _{{\rm{max}}}}} $的这些衍射光对应的传播方向为$f _{{{\mu} _{{\rm{max}}}}} = {1 / X} = {{\cos {\alpha _{{\rm{max}}}}} / \lambda }$. 这些衍射光出现在衍射观察屏上坐标的位置为$f _{{{\mu} _{{\rm{max}}}}} = {{\sin \left( {{\text{π}}/2 - \alpha } \right)} / \lambda } \approx {1 / \lambda } \cdot ({{{{{L_{x{\rm{max}}}}} / 2}} / d})$, 即衍射屏上的最大坐标为${L_{x{\rm{max}}}} = {{{N_x}\lambda d} / {{L_{x0}}}}$. 由此可见在使用S-FFT算法时, 衍射观察屏的大小不仅是波长的函数, 而且还是再现距离和抽样数的函数, 当$d$很小, 如果保持取样数不变, 则观察屏上的再现结果只是邻近光轴的很小区域, 所以该算法主要适用于衍射距离较大的情况.
3
3.1.2.菲涅耳衍射积分D-FFT算法
-->同S-FFT算法相类似, 在直角坐标系中, 傍轴近似下, 当衍射距离满足菲涅耳衍射条件时, 根据衍射积分表达式(10)可以得到其脉冲响应:
$\tilde h\left( {x,y} \right) = \dfrac{{\exp \left( {{\rm{j}}k _\lambda d} \right)}}{{{\rm{j}}\lambda d}}\exp \left[ {\dfrac{{{\rm{j}}k _\lambda }}{{2d}}\left( {{x^2} + {y^2}} \right)} \right], $
显然, 菲涅耳衍射下的脉冲响应是不变的, 所以可以用卷积运算来完成, 即
$\tilde U\left( {x,y} \right) = {\tilde U}_0 \left( {{x_0},{y_0}} \right) \otimes \tilde h\left( {x,y} \right).$
由于在空域完成卷积运算需要花费很长的机器时间, 所以选择到频域完成. 假设菲涅耳衍射的传递函数为$\tilde H\left( {f _{\mu},f _v } \right)$, 则菲涅耳积分可以分为两次傅里叶变换完成, 所以称其为D-FFT算法. 可以通过傅里叶变换的性质直接得到菲涅耳衍射传递函数$\tilde H\left( {f _{\mu},f _v } \right)$的解析表达式:
$\begin{split} & \tilde H\left( {f _{\mu},f _v } \right) \\ =& \;F\left\{ {\dfrac{{\exp \left( {{\rm{j}}k _\lambda d} \right)}}{{{\rm{j}}\lambda d}}\exp \left[ {\dfrac{{{\rm{j}}k _\lambda }}{{2d}}\left( {{x^2} + {y^2}} \right)} \right]} \right\} \\ =& \;\exp \left\{ {{\rm{j}}k _\lambda d\left[ {1 - \dfrac{{{\lambda ^2}}}{2}\left( {{{f _{\mu} }^2} + {{f _v }^2}} \right)} \right]} \right\}. \end{split} $
综上所述, 菲涅耳衍射公式可以这样计算:
$\begin{split} \tilde U\left( {x,y} \right) =&\; {\rm{FF}}{{\rm{T}}^{ - 1}} \!\bigg\{\! {\rm{FFT}}\!\left\{\!{{{\tilde U}_0}\left( {{x_0},{y_0}} \right)}\!\right\} \\&\times \exp \left\{ {{\rm{j}}{k_\lambda }d\left[ {1\! -\! \frac{{{\lambda ^2}}}{2}\left( {{f_{\text{μ}}}^2 \!+\! {f_v}^2} \right)} \right]}\!\right\}\!\bigg\}\!.\end{split}$
D-FFT算法中根据已经确定的传递函数$\tilde H\left( {f _{\mu},f _v } \right)$的解析表达式, 只需要将频域中$x$方向和$y$方向的频率$f _{\mu} $$f _v $代入解析式, 就可以得到传递函数的数值. 假设实际问题是在宽度为${L_{x0}} \times {L_{y0}}$的空间采样, 并且有${N_x} \times {N_y}$个抽样点, 最高空间频率分别为$f _{{{\mu} _{{\rm{max}}}}} = {{{N_x}} / {2{L_{x0}}}},f _{{\nu _{{\rm{max}}}}} = {{{N_y}} / {2{L_{y0}}}}$. 所以, $f _{\mu} $的取值范围从$ { - f}_{{{\mu} _{{\rm{max}}}}} $$f _{{{\mu} _{{\rm{max}}}}} $${N_x}$个点, 同样, $f _v $的取值范围从$ - f_{{\nu _{{\rm{max}}}}} $$ f_{{\nu _{{\rm{max}}}}} $${N_y}$个点.
2
3.2.S-FFT算法和D-FFT算法的数值模拟结果的比较
-->在衍射计算过程中, 应该特别注意欠采样的问题, 在保证抽样数和光束波长不变的情况下, 改变衍射距离$d$, 比较S-FFT算法和D-FFT算法抑制欠采样发生的能力, 说明D-FFT算法在处理短距离衍射积分问题的优势.
图3(a)图3(b)分别为衍射距离为100 mm时利用S-FFT算法和D-FFT算法得到的衍射光强分布, 图3(c)图3(d)分别为衍射距离为10 mm时利用S-FFT算法和D-FFT算法得到的衍射光强分布. 通过对比图3(a)图3(b)发现: 当衍射距离为100 mm时, 利用S-FFT和D-FFT两种算法得到的衍射光强分布均未出现欠采样现象. 通过对比图3(c)图3(d)可知: 随着衍射距离的减小, 采用S-FFT算法得到的衍射光强出现了明显的欠采样现象, 但是D-FFT算法未出现, 可见D-FFT算法的抗欠采样能力要明显高于S-FFT算法. 结果表明, D-FFT算法不惧怕计算衍射距离近的问题, 即使衍射距离$d$ = 1 mm时, 用D-FFT算法依然可以得到很好的结果. 超高音速飞行器周围形成的等离子体鞘套的厚度属于厘米数量级. 由此可见, 用FFT算法数值模拟光束在相邻两屏之间的传输过程时, 选用D-FFT算法更合适.
图 3 在不同衍射距离下S-FFT和D-FFT两种算法的光强分布 (a) S-FFT, $d$ = 100 mm; (b) D-FFT, $d$ = 100 mm; (c) S-FFT, $d$ = 10 mm; (d) D-FFT, $d$ = 10 mm
Figure3. Intensity distribution of S-FFT and D-FFT algorithm at different diffraction distances: (a) S-FFT, $d$ = 100 mm; (b) D-FFT, $d$ = 100 mm; (c) S-FFT, $d$ = 10 mm; (d) D-FFT, $d$ = 10 mm.

利用谱反演法生成等离子体鞘套湍流随机相位屏, 首先需要在频域生成一个均值为1, 方差为0的$N \times N$复随机矩阵$H\left( {{k_x},{k_y}} \right)$, 利用等离子体鞘套湍流功率谱[18]对其进行滤波, 然后再进行傅里叶逆变换便可以得到等离子体鞘套湍流的随机相位屏$\varPhi \left( {{k_x},{k_y}} \right)$. 等离子体鞘套湍流功率谱为
${\varPhi _n}\left( {{k_x},{k_y}} \right) = 0.078\langle n_1^2\rangle L _0^{{{ - 4} / 5}} {\left( {k_x^2 + k_y^2} \right)^{ - 19/10}}, $
其中$\langle n_1^2\rangle $为折射率方差, $L _0 $为湍流外尺度, ${k_x}$${k_y}$分别为不同方向的频率. 由折射率频谱可得到相位频谱为$\varPhi \left( {{k_x},{k_y}} \right) = 2{\text{π}}{k^2}{\varPhi _n}\left( {{k_x},{k_y}} \right)$, 则相位频谱的方差为${\sigma ^2}\left( {{k_x},{k_y}} \right) = {\left( {{{2{\text{π}}d} / {N{\Delta}x}}} \right)}^2 \varPhi \left( {{k_x},{k_y}} \right)$, 其中$d$是每两个随机相位屏的间距, $N$为采样数, $\Delta x$是栅格间距. 通过傅里叶变化便可得到空域中的随机相位屏[39]: $\varPhi \left( {x,y} \right) = {\cal F}\left[ {H\left( {{k_x},{k_y}} \right)\sigma \left( {{k_x},{k_y}} \right)} \right]$.
利用菲涅耳衍射积分的D-FFT算法求解高斯光束在相邻两屏之间的传输过程, 利用功率谱反演法生成的多随机相位屏模拟等离子鞘套湍流带来的影响, 多随机相位屏示意图如图4所示.
图 4 多随机相位屏示意图
Figure4. Schematic diagram of multi-random phase screen.

2
3.3.等离子体鞘套湍流中影响高斯光束质量的因子分析
-->图5的传输参数为: 束腰半径$ \omega _0$ = 40 mm, 相位屏间隔$d$ = 0.013 m, 相位屏采样数为$1024 \;\times $1024. 分析折射率方差、波长和传输距离对高斯光束质量的影响.
图 5 不同传输条件时的高斯光束光强数值仿真 (a) λ = 1.55 × 10–6 m, ω0 = 40 mm, d = 13 mm; (b) ω0 = 40 mm, d = 13 mm, $\langle n_1^2\rangle$ = 10–12; (c) λ = 1.55 × 10–6 m, ω0 = 40 mm, $\langle n_1^2\rangle$ = 10–12
Figure5. Numerical simulation of intensity of Gaussian beam in different propagation conditions: (a) λ = 1.55 × 10–6 m, ω0 = 40 mm, d = 13 mm; (b) ω0 = 40 mm, d = 13 mm, $\langle n_1^2\rangle$ = 10–12; (c) λ = 1.55 × 10–6 m, ω0 = 40 mm, $\langle n_1^2\rangle$ = 10–12.

图5(a)是高斯光束在折射率起伏方差分别为$\langle n_1^2\rangle = {10^{ - 13}},\,{10^{ - 12}},\;{10^{ - 11}}$时通过10个随机相位屏后系统归一化的三维(3D)和二维(2D)光强数值仿真结果. 由图5(a)可知, 在湍流折射率方差不断增大的过程中, 光斑分裂成形状不规则的若干小光斑, 光斑弥散程度加剧, 而且光强起伏逐渐增大, 强度逐渐减小.
图5(b)是高斯光束波长分别为$\lambda = 0.6328\; \times$ $ {10^{ - 6}},\; 1.05 \times {10^{ - 6}},\; 1.55 \times {10^{ - 6}}$ m时通过10个随机相位屏后系统归一化的3D和2D光强数值仿真结果. 通过对比分析可知, 随着波长的不断增加, 高斯光束抵抗湍流的能力逐渐增强, 具体表现在光斑弥散程度减小、光强起伏减弱上, 所以在等离子体鞘套通信系统中选用波长较长的激光光束作为信息载体比较好.
图5(c)是高斯光束沿z轴传输时距离分别为6.5, 13, 19.5 cm, 即分别通过5, 10, 15个随机相位屏后系统归一化的3D和2D光强数值仿真结果. 结果反映出: 传输距离越大, 光斑弥散越严重, 光强起伏也越明显.
首先利用菲涅耳衍射积分的D-FFT算法解决了湍流问题中的短距离传输问题, 然后利用多随机相位屏数值模拟了高斯光束在等离子体鞘套湍流中的传输特性. 研究表明折射率方差越大, 高斯光束的光斑弥散程度越大, 光强越弱, 光强起伏越明显; 随着波长的不断增加, 高斯光束抵抗湍流的能力逐渐增强, 所以在等离子体鞘套通信系统中选用波长较长的激光光束作为信息载体比较好; 传输距离越大, 光斑弥散越严重, 光强起伏也越明显. 综上所述, 高斯光束在等离子体鞘套湍流中的传输质量与外界因素折射率起伏方差、传输距离及自身因素波长、束腰半径等密切相关. 本文对研究激光在湍流等离子体鞘套中的传输问题具有一定的意义.
相关话题/计算 空间 质量 统计 湍流

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 大质量样品厚度对中子活化定量分析的影响和修正
    摘要:滑石粉的主要化学成分是Mg3[Si4010](OH)2,通过检测面粉中24Mg和28Si原子核的含量,可以计算出面粉中滑石粉的含量.利用中子活化分析方法测量面粉中24Mg和28Si原子核的含量问题时,被测样品内部中子通量和能量随厚度的变化以及γ射线自吸收效应会对测量结果有较大的影响.利用MCN ...
    本站小编 Free考研考试 2021-12-29
  • 基于共心球透镜的多尺度广域高分辨率计算成像系统设计
    摘要:针对实时广域高分辨率成像需求,充分利用具有对称结构的多层共心球透镜视场大且各轴外视场成像效果一致性好的特点,设计基于共心球透镜的多尺度广域高分辨率计算成像系统.该系统基于计算成像原理,通过构建像差优化函数获得光学系统设计参数,结合球形分布的次级相机阵列进行全局性优化,提高系统性能的同时有效简化 ...
    本站小编 Free考研考试 2021-12-29
  • 惯性约束聚变装置中靶面光场特性的统计表征方法
    摘要:在激光驱动的惯性约束聚变装置中,常采用多种束匀滑手段对焦斑的时空特性进行调控.光传输链路中涉及的光学元件众多、传输变换复杂,往往导致光传输模型复杂,且在运用衍射光学方法分析焦斑形态和特征时面临大量的数据处理和计算,致使出现计算量大、计算效率低等问题,亟需寻求快速而简便的新方法来描述焦斑的统计特 ...
    本站小编 Free考研考试 2021-12-29
  • 硅纳米结构晶体管中与杂质量子点相关的量子输运
    摘要:在小于10nm的沟道空间中,杂质数目和杂质波动范围变得十分有限,这对器件性能有很大的影响.局域纳米空间中的电离杂质还能够展现出量子点特性,为电荷输运提供两个分立的杂质能级.利用杂质原子作为量子输运构件的硅纳米结构晶体管有望成为未来量子计算电路的基本组成器件.本文结合安德森定域化理论和Hubba ...
    本站小编 Free考研考试 2021-12-29
  • Ce-La-Th合金高压相变的第一性原理计算
    摘要:采用第一性原理计算对Ce0.8La0.1Th0.1在高压下fcc-bct的结构相变、弹性性质及热力学性质进行了研究讨论.通过对计算结果的分析,发现了合金在压力下的相变规律,压强升高到31.6GPa附近时fcc相开始向bct相转变,到34.9GPa时bct相趋于稳定.对弹性模量的计算结果从另一角 ...
    本站小编 Free考研考试 2021-12-29
  • 螺旋锥束计算机断层成像倾斜扇束反投影滤波局部重建算法
    摘要:螺旋锥束计算机断层成像(CT)作为常用的临床诊断工具,如何尽可能地减少其辐射剂量是热点研究领域之一.局部成像利用准直器减小射线直照区域,能够有效降低CT辐射剂量.然而,局部成像会造成投影数据横向截断,产生局部重建问题.现有螺旋反投影滤波(BPF)算法只能实现局部曲面重建,难以实现局部体区域重建 ...
    本站小编 Free考研考试 2021-12-29
  • 光声光谱检测装置中光声池的数值计算及优化
    摘要:利用光声光谱技术进行痕量气体的检测具有独特的优势,光声池是系统装置中最为重要的核心部件,它决定着整机性能的优劣.以一圆柱形共振型光声池为研究对象,基于声学与吸收光谱学的基本理论,建立了光声池声场激发的数学模型;利用数值模拟方法对光声池空腔结构进行了声学模态仿真,获得了前8阶声学模态值以及声压可 ...
    本站小编 Free考研考试 2021-12-29
  • 正弦波沟槽对湍流边界层相干结构影响的TR-PIV实验研究
    摘要:利用高时间分辨率粒子图像测速(time-resolvedparticleimagevelocimetry,TR-PIV)技术,在不同雷诺数下对光滑壁面和二维顺流向、三维正弦波(two/threedimensional,2D/3D)沟槽壁面湍流边界层流场进行了实验测量,从不同沟槽对湍流边界层相干 ...
    本站小编 Free考研考试 2021-12-29
  • W-In体系溶质晶界偏聚行为的第一性原理计算
    摘要:基于第一性原理构建了钨基合金体系的溶质偏聚模型,以W-In体系为例研究了不同浓度下溶质的晶界偏聚行为和成键特征,从电子结构层面揭示了W-In体系的键合作用,预测了W-In体系界面稳定性随溶质浓度的变化规律.结合键布居、电荷密度、差分电荷密度和态密度等电子结构分析,发现了W-In体系中溶质原子在 ...
    本站小编 Free考研考试 2021-12-29
  • 空间电子辐照聚合物的充电特性和微观机理
    摘要:空间电子辐照聚合物的充电特性和微观机理是研究和防护航天器聚合物充放电特性的基础.采用蒙特卡罗方法模拟空间电子的散射过程,快二次电子模型模拟二次电子的产生,有限差分法求解电荷连续性方程、电流密度方程和泊松方程的电荷输运过程,俘获过程基于Poole-Frenkel效应来实现.基于电子散射/输运同步 ...
    本站小编 Free考研考试 2021-12-29