气动电磁学概念探讨及若干研究进展介绍

于哲峰, 孙良奎, 马 平, 梁世昌, 张志成, 黄 洁

(中国空气动力研究与发展中心, 四川 绵阳 621000)

0 引 言

空气动力学是研究物体和空气之间有相对运动时，即物体在空气中运动或物体不动而空气流过物体时，空气的运动规律及作用力(空气内部的和空气对物体的)所服从规律的学科。这一学科是随着航空、航天事业的发展而发展起来的。空气动力学在向前发展的同时，也在不断地与其它学科相互融合、相互交叉，产生了气动热、气动光学、气动声学等交叉学科。

采用SPSS 20.0统计学软件对数据进行处理。计量资料以“±s”表示，采用t检验。计数资料以百分数（%）表示，采用x2检验。以P＜0.05为差异有统计学意义。

目前国内外航空和航天一些热点研究方向如再入/临近空间高超声速飞行器目标特性、通信中断、气动隐身一体化、等离子体隐身、等离子体流动控制和电磁推进研究等，形成既涉及到空气动力学又涉及到电磁学的交叉领域[1-8]。这些热点问题涉及到共同的基础理论、分析方法、试验设备和测试仪器，目前有一些概念如磁流体力学、电/磁空气动力学、等离子体空气动力学等来定义这些问题，但都不够全面。例如，磁流体力学是主要考察导电流体在磁场作用下运动规律的学科，也即考察磁场如何影响着流体的运动，反之流体的运动又是如何影响着电磁场[9]。然而磁流体力学未能包括再入/临近空间高超声速飞行器目标特性、通信中断、气动隐身一体化等需要考虑空气动力对电磁信号影响的问题。中国空气动力研究与发展中心张志成研究员在《气动物理学》一书中首次提出了气动电磁学的概念：研究高超声速飞行器的等离子体绕流场和尾流场的变化规律及其与电磁波的相互作用规律[10]。但是该概念未能包涵气动隐身一体化、等离子体隐身、等离子体流动控制和电磁推进等非高超声速涉及到的气动和电磁学问题。总之，现有的磁流体力学、气动电磁学等概念只能涵盖部分内容。

本文将原有的气动电磁学定义进行了扩展，将气动电磁学定义为既涉及到空气动力学又涉及到电磁学，研究飞行器流场与电磁场相互影响、相互作用规律的交叉学科，涉及的学科主要包括空气动力学、电磁学以及等离子体学等，研究方向主要包括再入/临近空间高超声速飞行器目标特性、通信中断、气动隐身、等离子体隐身、等离子体流动控制和电磁推进等。

我于是很怀疑王爷表演朱登春那个角色时，是不是也下意识将自己的心渡到了那角色身上。不过，他的痛不是儿失慈娘的悲痛，而是父失孝子——在王爷近七十岁时，他二儿子一次醉酒，再也没有醒来。

图1 气动电磁学主要研究方向Fig.1 Research directions of aero-electromagnetics

1 典型气动电磁学问题

高超声速飞行器再入大气层或在临近空间飞行时，在飞行器表面会形成等离子鞘套，在其下游形成电离尾迹，严重影响高超声速飞行器的电磁散射特性[11]。最初，只在再入段关心等离子体流场对高超声速飞行器电磁散射特性的影响，随着临近空间飞行器概念的提出，在整个飞行过程中都需要考虑这种影响。临近空间飞行器本身速度快、机动范围大、落点无法预测、临近空间气动环境杂等因素，利用现有防御系统实现临近空间飞行器的远程预警、跟踪识别、有效拦截、实时评估还有较大的难度。为了实现对临近空间飞行器的有效预警和探测，临近空间高超声速飞行器雷达散射特性机理研究、相关的雷达探测跟踪理论方法和探测雷达组网方式等成为近期国内外研究的热点。

1.1 高超声速飞行器电磁散射特性

下面介绍一些航空航天领域比较典型的气动电磁学问题。

1.2 通信中断

高超声速飞行器再入大气层和在临近空间飞行时，在飞行器表面会形成等离子鞘套，电磁波通过等离子体鞘套传播过程中被等离子体鞘套反射和吸收，波束强度发生衰减，并出现偏折、延时、相移等现象，甚至可能发生电磁波传输的完全中断，即“通信中断”或“黑障”，严重影响高超声速飞行器的飞行控制和安全。近几年，随着临近空间高速飞行器的发展，国内又掀起了等离子鞘套下信息传输研究的热潮，在等离子体流场精确诊断、试验设备上建立暗室环境消除复杂环境对通信中断试验的电磁干扰、真实信号在等离子体鞘套中传输特性等方面取得了新的研究进展。采用的提高通信频率(采用毫米波和太赫兹通信)、激光通信、外加强电磁场干预、低频磁场波通信、飞行器外形和弹道优化等方法有望部分解决通信中断问题[12]。

1.3 气动隐身一体化

日益兴盛的小卫星和微卫星和微型行星探测器要求星上的推进系统质量更轻、体积更小和效率更加高，因此比冲高、结构紧凑、消耗工质少的电推进引起航天器设计师们的注意和青睐[16]。电推进是使用电能加速的工作介质(工质)形成高速射流而产生反作用推力的技术。电推进按其工质的加速方式分为三种类型:电热式、静电式和电磁式。在提高推力方面，近年来提出了大推力螺旋波等离子体发动机概念，它采用螺旋波电离产生高密度等离子体，并通过离子回旋共振加热方式高效吸收射频能量，将电磁波能量转化为离子的动能，形成高温等离子体流并在洛伦兹力、气动力和静电场等多物理场作用下加速，形成高速等离子体射流产生推力。

1.4 等离子体隐身

飞行器的气动性能和电磁性能涉及两个不同学科的目标，需要根据任务的具体要求，建立描述气动与电磁性能的目标函数，对飞行器的气动特性和电磁特性进行综合分析与优化。

1.5 等离子体流动控制

从节能和最大限度地提高飞行器性能这两种观点来看，减阻都具有最大潜力。用等离子体减小飞行器的气动阻力是近年来国外提出的又一新概念。其主要方法有：逆向喷流减阻、边界层控制减阻和局部能量点源减阻。国内外的一些风洞试验和计算结果已证明了等离子体减阻效果的显著性[15]。近年来随着超燃冲压发动机的发展，为了充分利用超燃冲压发动机在高超声速条件下比冲最大的优点，需要保证发动机在高超声速飞行条件变化范围内有效地工作，为此，就要求研制受控的进气道和燃烧室，以便提供稳定的超声速燃烧。以俄罗斯莫斯科大学、西伯利亚分院理论与应用力学所为代表的研究机构使用电磁流体动力系统实现发动机稳定点火、增强燃烧混合，提高推进系统热电力循环效率，并且能提高其比冲和推力。

对于我国来说，制定促进农户适应新型农业食品体系的政策是非常必要的。政策旨在解决农户参与超市供应链时可能遇到的营销壁垒，以及促进集体行动，从而尽可能地降低农户参与过程中的交易成本。

1.6 电推进

在现代战争中为了提高飞行器的突防能力，除要求飞行器具有良好的空气动力学特性外，还要求飞行器具有良好的隐身性能。也就是说，飞行器设计必须降低雷达散射截面(RCS)。飞行器RCS大小不仅受电子设备、材料结构等因素的影响，而且与飞行器外形有密切的关系，因此，低RCS外形是实现飞行器高隐身性能的基础与关键之一。美国SR721侦察机、B2轰炸机和F22战斗机等都采用了气动隐身一体化设计。在理论分析方面，随着计算机性能的不断提升和计算方法的不断改进，基于数值求解N-S方程和Maxwell方程的飞行器气动与隐身特性精确数值模拟方法被广泛使用。对进气道、尾喷管和弹仓等重点部件或位置气动隐身特性的研究不断深入。在测试方面，重点发展了全尺寸飞行器模型雷达散射特性测量、飞行器模型气动隐身特性自由飞模拟等试验技术[13]。

2 气动电磁学涉及的主要机理问题

真空中相对磁导率和相对电导率都是1，而当空气电离形成等离子体时，其相对磁导率仍然为1，相对介电常数εr则变为复数，并且有如下形式:

2）标准溶液浓度的不确定度。根据经验，由于稀释引起的不确定度很小，可忽略不计，故砷元素形态标准溶液配制的不确定度即为砷元素形态标准溶液的不确定度。

西施在男性文化中被塑造成美女的形象，而在女诗人笔下，却凸显其帮助越国实现霸业的不朽功绩。借助对历史的重新评判，诗人也实现了自身的定位。其72岁高龄所作的《琴楼遗稿序》里，还明确提出了女性生存价值的衡量标准：

2.1 高温等离子体流场特性研究

气动电磁学除了会涉及到常规的低速、高速和高超声速空气动力学外，还往往涉及到高温等离子体流场特性研究。空气在达到一定温度时将发生离解和电离，分子和原子的内能模式有不同程度的激发并形成高温气体流场，过程中非常复杂。流场高温化学非平衡过程涉及到各种高温气体效应，包括化学平衡/非平衡/冻结效应和热力学平衡/非平衡/冻结效应等。

2.2 气动对电磁特性的影响

气动电磁学涉及到空气动力学、电磁学和等离子体物理学等多个学科，有许多机理问题需要解决。

电弧风洞、高频等离子体风洞、高温激波管/粉末激波管和弹道靶四种设备产生等离子体的典型参数及特点由表2给出。

2.3 电磁对气动特性的影响

电磁对气动特性的影响主要包含两个方面：一是电磁场可以向流场中注入能量，使飞行器周围的流场结构(波系结构和边界层结构状态)发生变化，致使飞行器的气动特性和物理特性发生改变，从而改变飞行器的气动性能；二是流场中的带电粒子在电磁场中将受到电场力和磁场力的作用， 利用电磁力可以控制边界层，减小边界层的湍流度和/或防止湍流涡系的产生，使流体加速产生推力。

2.4 气动特性和电磁特性的综合分析与优化

等离子体隐身技术作为一种全新的隐身概念，具有吸波频带宽、效率高、使用简便、价格便宜等优点，特别是应用于飞行器隐身时无须改变飞行器的外形，解决了隐身措施与气动性能的矛盾，日益受到国内外国防决策机构和军事专家的关注。据报道，俄罗斯已成功地将等离子体隐身技术应用到新型战斗机T-5上，美国和加拿大等国也积极开展相关研究。近年来，电磁波在等离子体中的衰减已经被理论分析和数值模拟结果广泛证实，研究者们目前更关心如何在工程上实现飞行器的等离子体隐身。由于开放式的等离子体的可控性与稳定性都较差，因而越来越多的研究者开始关注如何利用封闭式的发生器来产生等离子体并实现对飞行器的隐身[14]。

3 主要试验装置和测量技术

开展气动电磁学试验时，除了利用常规风洞和微波暗室等气动和电磁研究设备外，往往还需要等离子体地面模拟装置。从机理上现有的比较实用的等离子体地面模拟装置可以分为两类：一种为热激发式等离子体，利用高温产生的等离子体；另一种为电激发式等离子体，利用气体放电产生出等离子体。前者主要包括激波管、燃烧型等离子体发生器和弹道靶等。后者包括电弧风洞、等离子体风洞和其它各种气体放电型等离子体发生器。不同类型的气体放电等离子体发生器参数如表1所示。

要创设“说”书的时机。教师要给学生创设“说”课外书的时机，让他们来交流讨论自己读过的课外书，可以开课外书交流会、“说”书会等。有时可以围绕一本书来“说”，有时可以不定主题地“说”，最重要的是让学生们能围绕着书大胆自由地说起来就行了。

(1)填补监测空白，提高站网密度。持续开展地下水监测站点建设工作，填补鲁中南中低山丘陵区和鲁东低山丘陵区部分地区裂隙岩溶水和基岩裂隙水的监测空白，加密地下水降落漏斗区、地面沉降区、岩溶塌陷区和大中型水源地的目标含水层监测点密度[12]，达到《区域地下水位监测网设计规范》(DZ/T0271-2014)中对地下水位监测点密度的要求，进一步优化监测网络，使站网布局更趋合理。

表1 不同类型的气体放电等离子体发生器参数Table 1 Parameters of the different plasma generators

表2 四种设备产生等离子体的典型参数Table 2 Plasma parameters of four typical facilities

气动电磁学还涉及到一些共用的测量技术，主要包括流场测量技术、等离子体诊断技术和电磁特性测量技术等。

4 已开展的工作

课题组前期在再入/临近空间高超声速飞行器雷达散射截面(RCS)物理靶试验与数值模拟、磁窗天线增强等离子体鞘套透波特性、喷流等离子体隐身技术和钝锥体气动隐身一体化设计方面开展了一些研究，下面分别加以简单介绍。

4.1 再入/临近空间飞行器RCS试验与模拟

课题组为了研究等离子体对高超声速模型RCS的影响，利用气动物理靶对高超声速模型的RCS进行了测量，给出了雷达接收信号的频谱和模型RCS的一维距离像[17-20]。针对试验状态，对飞行器本体及绕流RCS特性进行了模拟，分析了等离子体鞘套及尾迹对模型电磁散射特性的影响。钝锥模型和尾迹流场纹影照片如图2(a)所示，气动物理靶试验设备如图2(b)，模型本体和尾迹RCS试验与仿真结果对比如图2(c)所示。在此基础上，系统地对升力体外型临近空间高超声速飞行器RCS特性开展了研究[17]。利用临近空间高超声速飞行器后面会形成几百米、甚至上千米的等离子体尾迹流场，对低于等离子体频率的电磁波产生很强散射的特点，探讨了利用天波超视距雷达探测类HTV2临近空间高超声速飞行器的可行性。

(a) 高超声速模型流场气动物理靶试验阴影照片(中国空气动力研究与发展中心谢爱民高工提供)

(b) 气动物理靶

(c) 尾迹RCS试验与仿真结果对比

4.2 磁窗天线增强等离子体鞘套透波特性研究

外加强磁场可以形成“磁窗”，有效地增强右旋圆极化电磁波在等离子体鞘中传播时的透波特性，但要

形成磁窗还存在磁体的重量和复杂性等问题，为解决这一问题，提出了天线和强永磁体一体化综合设计的新方法——磁窗天线，使得磁窗天线既具有天线发射通信信号的功能，又具有产生强磁场改善等离子体透波特性的功能，并且天线和永磁铁之间不会相互影响[21-22]。磁窗天线结构图如图3(a)所示，周围磁场强度分布如图3(b)所示，典型状态下增强等离子体鞘套透波特性效果如图3(c)所示。

(a) 一体化磁窗天线

(b) 磁窗天线周围磁场强度分布云图

(c) 磁窗天线与普通天线的对比

4.3 喷流等离子体隐身技术研究

利用微型固体火箭发动机作为等离子体发生器产生非均匀等离子体流场，分析了等离子体流场对平板RCS特性的影响[23-24]。等离子体发生器和微波反射仪测量示意图如图4所示，接收机输出功率随时间变化曲线如图5所示，利用时域有限差分方法计算得到的RCS随测量角度变化曲线如图6所示。

图4 等离子体发生器和微波反射仪测量示意图Fig.4 Plasma generator and the sketch of microwave test

图5 接收机输出功率随时间变化曲线Fig.5 Curve of the reciever output power vs. test time

图6 RCS随双站角变化曲线Fig.6 RCS vs. bistatic angle

4.4 气动隐身一体化设计

从多目标、多学科综合分析的角度出发，结合实际工程问题，研究了缩比模型和小翼弹头的隐身与气动综合特性，提出了两种综合性能更优的弹头外形[25]。不同底部形状的压力等值线如图7所示，底部形状分别为锥体和椭球体时RCS随椭球体轴长和锥体高变化曲线如图8所示，零升阻力系数随椭球体轴长和锥体高变化曲线如图9所示。

以“2.3.2”项下拟合模型为目标函数，使用Matlab 2014b软件的GADS工具箱，采用遗传算法求解指标成分的最优提取工艺条件。遗传算法工具GUI参数设置见表5；随机搜索10次，结果见表6。由表6可知，10次随机搜索结果的变异范围较小，对目标函数最佳值的逼近程度较好。因此，取10次随机搜索的平均水平即60.872 000%乙醇提取3.105 564 h、液料比15.719 832∶1（mL/g），预测综合评分为2.895 433。

(a) 普通钝锥模型

(b) 底部为锥体，a=8 mm

(c) 底部为椭球，a=8 mm

图8 RCS随椭球体轴长和锥体高变化曲线Fig.8 RCS of the model with ellipsoidal or cone bottoms

图9 零升阻力系数随椭球体轴长和锥体高变化曲线Fig.9 Resistance coefficient of the model at 0°attack angle

5 总 结

气动电磁学涉及的相关研究在航空和航天领域具有巨大的应用潜力，成为近年来世界空气动力学领域研究的新热点。但是应该看到，由于气动电磁学是一门涉及空气动力学、电磁学、等离子体学、热力学和化学动力学等许多学科的交叉学科，涉及的专业面很宽，技术复杂，因此，研究难度也很大。在基础理论、数值模拟、试验设备与测试技术等方面有许多难关需要攻克。

今后需建设可以开展各类气动电磁学问题研究的通用试验和仿真平台。建立复杂外形高超声速飞行器雷达散射特性、气动传输效应、强电磁场流动控制试验方法和仿真手段。形成高超声速飞行器等动目标光电特性分析评估能力、动目标等离子体绕流对通信导航电磁波影响分析评估能力、强电磁场产生和流动控制能力，气动隐身一体化综合设计与评估能力。为新一代航空、航天飞行器气动、隐身和通信等方面综合性能提升提供技术支持。