连续激光内通道传输的弱可压缩流模型.pdfVIP

第# 卷 第$ 期 强 激 光 与 粒 子 束 %’( # ，)( $ $**+ 年$ 月 ,-., /0123 45623 5)7 /538-942 :25;6 7= ， $**+ 文章编号： ** ?@$$（$**+ ）$ A#@ *? 连续激光内通道传输的弱可压缩流模型! $ ，@ $ @ 柳 建 ， 金 钢 ， 王世庆 ， 李树民 ， 张 翔 （! 成都理工大学工程技术学院，四川乐山+?**B ；$! 中国空气动力研究与发展中心，四川绵阳+$*** ； @! 中国科学院光电技术研究所，成都+*$*A ） 摘 要： 分析了激光在气体中传输时采用等压近似线性方程求解流场密度分布的优缺点，在高低速流场 统一计算模型的基础上提出了基于压力原变量的分步求解的弱可压缩流计算模型，并分析了该模型的特点。 采用该模型结合标量衍射理论对连续激光在封闭充气管道中受到的气体热效应影响进行了数值仿真。仿真结 果与实验结果的对比表明，弱可压缩流计算模型能更精确地反映非自由边界热对流对气体密度分布的影响，进 而反映流场对光束的影响。这说明弱可压缩流计算模型能较好地适应内通道光传输问题的仿真研究。 关键词： 光传输； 气动光学； 热效应； 标量波动方程； 弱可压缩流 中图分类号： 8)*$ ； 0?@+ ； 0@C ； 0$? 文献标识码： 5 ［，$ ］ 在激光大气传输研究中，气体密度分布常由等压近似的流体控制方程获得 。在开放边界环境压强近 似恒等，大气对流速度较大（D E F 量级）时，其仿真结果与实验结果符合很好。但是等压近似忽略了气体粘性、 压差和体积力等作用，难以处理粘性固壁、压差边界、自然对流等问题。而充气管道中的低速流问题，粘性固 壁、流固温度场耦合是影响流场的重要因素，且管内流场可能为湍流运动，进一步使粘性不可忽略；此外，由于 没有大气环流，由于密度变化（体积力作用）导致的流动变化将占据显著地位；而且，在封闭条件下，尽管压力 变化是小量，但压力变化引起的流场密度变化却不一定是小量。以上原因使得采用等压近似方程时所得的光 束波前畸变结果较为粗糙。金钢、刘顺发等人采用了改进的可压缩流模型代替等压近似模型［@ + ］，获得了比较 精细的结果。可压缩流模型能处理任意流场问题，但采用通常的标准方法求解低速弱可压问题时，由于系统出 现刚性而会产生种种困难，而且模型求解复杂，计算量很大［B ，# ］。 本文在1FF’GHI 等人的高低速统一模型基础上，采用音速归一化流速的压力原变量方法，消除了在马赫 数趋于零时压力项奇性和右端源项消失的困难，得到了适合于低速弱可压流问题的计算模型（权且称为弱可 ［A ］ 压缩流模型） 。采用该模型对封闭充气管道中流场对光束的影响进行仿真计算，并与文献［C ］的仿真结果及 实验结果进行了比较。结果表明该模型可以更精细地捕捉流场密度的细微变化，明显反映出密度细微变化引 起的对流效应以及固壁边界的影响，所得远场光斑与实验结果更趋一致，且计算效率较高。 弱可压缩流模型及其求解方法 流场密度分布可以由完全)JKGL6MNF（)6 ）方程组描述，不同的计算模型都是针对特定的情形对该方程 组的变形和简化。对于低速管道流的计算，目前的解决方法大致有两种：条件化预处理方法和高低速问题统一 方法［A * ］。本文对1FF’GHI 的统一格式进行改进，得到了适用于光、流耦合计算的弱可压缩流计算模型。 为表述简便，从2O’L 方程出发（粘性项在最后添加，文中从略）。把状态方程改写为 P ! （#$ ），其中 ，