平衡流场的再入飞行器电磁散射特性分析.docVIP

平衡流场的再入飞行器电磁散射特性分析 当具有极高速度(如10个马赫数以上)的飞行器再入大气层时,由于目标与空气摩擦将产生高达几千摄氏度的气动热,使周围的气体发生电离,导致飞行器附近空气呈离子状态存在,形成等离子体鞘套和冗长的等离子尾流。尽管作为一种色散介质的等离子体具有“通高频、阻低频”的特性,即大于等离子体频率的电磁波可以在等离子体中传播,而小于等离子体频率的电磁波被等离子体反射,但对再入飞行器来说,不同的再入速度对等离子体尾流会产生何种影响,等离子体尾流内的电子密度会达到何种量级,高电子密度的等离子体尾流对低频电磁波能否表现出强散射特性,从而有利于雷达的探测与识别,这些都是研究再入飞行器电磁散射特性时值得深入探讨的问题。 早在20世纪60年代初国外就已开展了与等离子体尾流相关的研究,鉴于等离子体尾流情况复杂。在理论研究方面,有用Born近似方法计算等离子体尾流的电磁散射特性,建立了再入尾流散射的畸变波Born近似模型。21世纪初期,国内学者也基于Born近似方法开展了大量有关再入段等离子体尾流散射特性的研究,但由于Born近似方法更适合于计算亚密(等离子体频率小于雷达波工作频率)状态下等离子体与雷达波的相互作用。因此,研究的频段主要集中在L和S波段。近年来,也有一些国内外学者利用电磁场数值计算方法研究了等离子体与电磁波的相互作用机理及其电磁特性,但利用该方法研究 再入飞行器等离子体尾流低频电磁散射特性的论文却鲜见发表。因此,本文根据再入飞行器的物理现象,将平衡流场的计算方法与电磁散射数值计算方法相结合,用于再入飞行器低频电磁散射问题的分析。首先借助真实气体效应情况下等离子体流场计算方法,获得锥球形目标再入时接近于真实尾流的非均匀等离子体分布,然后利用移位算子时域有限差分法(finite difference time domain FDTD)计算和分析锥球形目标以零攻角再入时的低频电磁散射特性,最后给出了一些有价值的结论。 一、再入飞行器周围等离子体流场计算方法。 等离子体作为一种特殊的色散介质,若不考虑外加磁场,与其介电常数有关的参数有入射电磁波频率、等离子体振荡频率和等离子体碰撞频率，即 (1) 式中,;N为等离子体中自由电子密度,e为电子电量;m为电子质量;为真空中介电常数。由此可见,在对再入飞行器进行电磁散射建模时,必须首先计算得到飞行器周围等离子体的流场特性,如内部自由电子密度等。 当飞行器以超音速再入时,稳定的空气被排挤开,但由于空气是按音速传播压力,因此被排挤的空气以高于音速的速度被强行挤压而形成激波。空气受到激波的强烈压缩和粘性的剧烈摩擦,激波层内的温度迅速升高,可达数千度,驻点区域甚至达到10 000 K以上。在如此高的温度条件下,空气内部中不仅分子的平动能、转动能被激发,甚至振动能也将被激发,而且还会发生气体分子的离解和电离,气体成分亦随之发生变化,空气呈等离子状态,这种现象称为真实气体效应。真实气体效应不仅直接影响飞行器的气动特性和热负荷,而且还对流场的光电特性产生重要影响。 在常温常压下,若不计微量元素,可以认为空气由和两种成分组成。然而,在高温(如温度为2 000 K~10 000 K的范围内)下, 和将发生化学反应。发生化学反应后,空气成为多组元的混合气体。针对不同的工程应用目的,选择空气混合气体组成成分的方法往往不同。对再入尾流的光、电特性的研究而言,较合理的方案是选择以下组元的化学成分:,,,,,,,,, ,。 在高超声速飞行器的流场计算中,必须考虑流动的物理化学状态。根据流动的物理化学状态的不同,可以将高温流动分为冻结流动、非平衡流动和平衡流动三类。由于本文的研究条件是30 km的高度和3 km/s~5 km/s的速度,其基本处于平衡状态。因此,可以近似用化学平衡模型来对整个流场开展分析。 在热力学平衡状态下,气体的密度、温度、压力、内能、焓、熵这六个特性参数中,任意两个可以视为独立变量,其余特性参数则可以通过实验或理论由这两个独立变量来确定。由气动热力学可知,平衡气体的化学组分仅与两个独立状态变量有关,而与发生反应的过程无关;当气体的状态参数确定时,其组份也相应地被确定。在给定压力P和温度T的条件下,若经过足够长的时间达到化学平衡状态,那么此时组元有各自确定的数量。显然,各组元的数量是压力P和温度T的函数。 对任何化学平衡混合气体,求解其化学组份时,首先需查明此化学系统包含的组元个数及元素个数(电子也当作一个单独的元素)。如果混合气体含有M个组元和N个元素,则应有M-N个独立的化学反应方程和相应的化学反应平衡常数,其余N个方程为Dalton分压定律和元素守恒定律(包括电荷守恒定律)。这样,通过M个独立方程可以求解各组元分压,由各组元分压,可得到混合气体中M个