基于电子束电离的磁流体力学进气道流动控制数值模拟

采用基于电子束电离的磁流体力学(MHD)控制系统,对高超声速流场附面层,以及非设计状态下的高超声速进气道流场的磁流体控制进行了深入研究.控制方程为低磁雷诺数Navier-Stokes方程,采用等离子体动力学模型与电子束模型模拟空气电离过程.研究结果表明:①电子束电离能有效提高流场的电导率,增强磁场对流场的控制效率;②基于电子束诱导电离的MHD控制系统能有效地控制高超声速流场的附面层,但其控制效率跟电子束能量大小相关;③基于电子束诱导电离的MHD控制系统能有效地改变非设计状态下高超声速飞行器的斜激波结构,使进气道重新满足Shock-on-lip(SOL)条件,但进气道的总压恢复系数以及流量将会降低.      

超声速氩气流磁流体发电初步实验研究(英文)

利用激波风洞,采用氦气驱动氩气,在平衡接触面运行方式下得到高温气体,通过在低压段注入电离种子K2CO3粉末,实现高温条件下导电流体的产生,开展了超声速氩气流磁流体功率提取初步实验研究。在喷管入口总压0.32MPa、总温6504K,磁场强度约0.5T、喷管出口气流速度1959m/s的条件下,对分段磁流体功率提取通道电极的感应电压和短路电流进行了测量,实验测量结果与理论计算相吻合,并由电压电流计算得出了平均电导率约20S/m左右,在负载系数为0.5的情况下,磁流体功率提取通道最大的功率密度可达4.7971MW/m3,最大焓提取率为0.34%。最后分析并给出了气体状态参数T1,M1,T2,M2的测试原理与方法。     

高超声速等截面磁流体发电机性能研究

采用三维低磁雷诺数磁流体动力学五方程模型,对等截面管道法拉第型分段电极磁流体发电机内的流动进行数值模拟,研究了负载系数、磁场强度、电极-绝缘壁宽度比对发电通道性能的影响.结果表明:当负载系数为0.50时电功率密度达到最大,并且随着负载系数的增大发生器的减速能力降低;当磁场强度为8T时,500mm的管道长度即可将来流马赫数从6降低到3以下,温度保持在800K以下,并且每立方米的管道发电量能够达到兆瓦级的水平.      

法拉第型磁流体发电机试验和数值仿真

针对未来航空航天任务对大功率空间电源的迫切需求，开展了国内首次高温惰性气体法拉第型磁流体发电机试验研究。试验采用电弧加热器作为模拟热源，以氩气作为工质，添加铯作为电离种子以提高工质电导率，成功实现了对法拉第型磁流体发电机的原理性验证，在1 T磁场环境的试验条件下取得了最高194 W的发电功率，功率密度为866 kW/m³。根据试验条件对发电过程进行了三维数值模拟，分析结果表明：发电机输出性能受电极压降和工质速度的影响较大，需要在后续研究中改进发电机工艺以降低电极压降，并对加速喷管重新进行设计。     

电导率模拟对高超声速MHD控制影响

高温气体电导率是高超声速电磁流动控制数值模拟最重要的参数之一。针对电导率模拟准确性及其对高超声速磁流体控制影响的问题，考虑高超声速飞行器流场中高温气体热化学非平衡效应，采用三维低磁雷诺数磁流体动力学（MHD）数值模拟方法及程序，结合国内外常见电导率处理方法开展典型状态高超声速MHD控制数值模拟，分析电导率模拟对高超声速磁流体流场分布、气动力/热特性的影响。研究表明：磁控热流减缓效果与电导率呈非线性关系，电导率较大时将出现电导率的磁控热饱和现象，其产生的原因可能与化学反应趋向于平衡态存在一定关系；采用定电导率方法，会人为放大磁场洛仑兹力的磁阻力效果，使阻力系数的预测值偏大；不同电导率模型计算得到的电导率分布差异很大，甚至存在数量级的差别，显著影响了磁流体的控制效果，这与电导率模型的适用范围、参数选取原则存在很大关联；对于含多种离解、电离组分的高温气体流动来说，采用基于多电离组分迁移碰撞的电导率模型（本文模型M8），计算与试验一致性最好。     

表面磁流体气动激励控制楔面激波规律数值研究

针对表面磁流体(MHD)气动激励对高超声速进气道在非设计状态下激波控制问题,从唯象学的角度出发,基于低磁雷诺数假设,将电磁作用简化为Navier-Stokes方程组中的源项处理,同时考虑到低气压、低磁场环境下电子回旋效应引起的Hall效应,并联立Ohm定律,建立磁流体动力学模型,通过与实验纹影对比验证了模型的合理性,并利用该模型研究了表面MHD加/减速激励作用位置与宽度、磁场强度、电导率和能量转化率等参数对楔面激波的影响规律。结果表明:表面MHD气动激励包括焦耳热与洛伦兹力作用,当放电功率密度为3.8×10~9k W/m~3,磁场强度为0.34T时,MHD加/减速激励分别使激波位置前移6mm与10mm;而磁场强度较低时,由于焦耳热的主导作用,将会出现激波前后压力比增大,激波强度增加等负面效应;根据激励参数的影响规律,激励器电极应靠近尖端布置,增大磁场强度,并改善等离子体源以提高气体电导率,同时适当增大激励区域宽度。     

横截面形状对超声速磁流体发生器性能的影响

通过三维数值模拟研究了不同横截面形状的磁流体(MHD)发生器的热、电磁流动特性,发生器入口采用超声速来流条件,旨在应用于磁流体旁路超燃冲压发动机推进系统。数值模拟采用低磁雷诺数下磁流体五方程模型,通过熵条件格式和超松弛迭代(SOR)法联合算法分别求解Navier-Stokes方程组和电势方程。研究表明在电极壁面上电场、电流以及洛伦兹力等电磁场参数呈周期性变化,由于电极端点效应,电磁场参数在电极端点出现周期性的极值,而在管道中心电极端点效应对电流影响并不大,沿着流向电流保持连续变化。对不同管道横截面形状的发生器的数值模拟表明:由于涡电流和二次流效应的影响,形状参数较小的发生器焦耳热较为严重,二次流较弱;而形状参数较大的发生器焦耳热较小,二次流增强。因此,当形状参数α=0.8时可获得较好的焓提取率和电效率。     

高温磁流体动力技术实验系统设计与调试

为了开展高温磁流体功率提取、磁流体加速、磁流体流动控制等磁流体动力技术实验,研制了基于高温燃气发生器的磁流体动力技术实验系统。介绍了高温磁流体动力技术实验系统的基本组成,设计思想及运行情况。燃气发生器采用航空煤油和气氧燃烧,产生最高温度可达3256K的高温燃气,通过添加辅助电离种子K2CO3,可以得到电导率为11.6S/m的导电流体;高温燃气经喷管的加速,在出口得到马赫数1.5的超声速导电流体;实验系统调试实现了15s以上稳定运行,表明具备了开展较长时间高温磁流体动力实验的能力。     

磁场对霍尔推力器放电特性及性能参数影响的粒子模拟研究

磁场强度及位形对霍尔推力器放电过程有显著影响。根据霍尔推力器通道尺寸和等离子体放电过程建立二维物理模型,采用粒子模拟方法,研究了不同磁场强度及位形等离子体放电特性,讨论了推力、推功比及放电电流的变化规律。模拟表明:当中轴线磁场强度峰值小于200G时,磁场对电子轴向传导约束减弱;当磁场强度峰值在200G～420G时,电子温度、电离率及电子与壁面碰撞频率降低,出口处离子径向速度增大,壁面腐蚀增加;当磁场强度峰值为280G时,加速区最短,放电电流最小。不同零磁点磁场位形会改变通道电离区和加速区位置,影响推力器放电性能。     

磁流体管流准一维模型的数值验证与性能分析

为了快速分析磁流体管流,采用涡轮喷气发动机模型对磁流体管流进行类比分析,引入变截面管流理论和修改过的多变效率表达式,推导了磁流体发电器和加速器的准一维模型,采用二维数值方法验证了准一维模型的可行性,利用该准一维模型对磁流体管流进行了性能分析,结果表明,多变效率的引入使准一维模型更加完整,当负载系数不变时,磁流体装置中的平均马赫数越低,其等熵效率越高。     

ν-gap度量及其在飞行控制律评估中的应用

 传统控制律评估方法主要用于单输入单输出（SISO）系统，且对模型参数摄动考虑不够全面,针对这些不足，研究了ν-gap度量方法。在介绍系统广义稳定裕度相关概念的基础上，给出了ν-gap度量的定义、特点和性质以及近似摄动模型的计算，提出ν-gap度量评估控制律的步骤。实例结果表明，该方法不仅克服了上述传统评估方法的缺陷，而且还有根据所求的各摄动影响情况忽略影响小的元素，以减少计算量及可以找到最坏情况下的参数摄动组合等优点。     

基于飞机油箱模型形状特征油量测量切片步长选择方法研究

 在分析飞机数字式油量测量过程中目前广泛使用的切片法油量测量原理的基础上,针对现有的定步长切片法无法得到准确、可靠的燃油质量特性数据库的缺陷,结合对飞机油箱模型形状特征的分析,提出了基于飞机油箱模型形状特征的油量测量切片步长选择方法。此方法包括切片步长整体和局部选择两个过程,整体选择以实现相邻两切片平面所夹油箱模型体积近似相等为目的来确定切片步长,以体现油箱模型截面整体变化规律;局部选择以设计切片平面与截面突变平面重合或尽可能接近的方式,突出油箱截面的局部变化特征。实验结果表明：该切片步长选择方法较定步长方法能够建立更为合理、可靠的燃油质量特性数据库,从而提高了油量测量精度。     

Integrated Entry Guidance for Reusable Launch Vehicle

A method for the implementation of integrated three-degree-of-freedom constrained entry guidance for reusable launch vehicle is presented. Given any feasible entry conditions, terminal area energy management interface conditions, and the reference trajectory generated onboard then, the method can generate a longitudinal guidance profile rapidly, featuring linear quadratic regular method and a proportional-integral-derivative tracking law with time-varying gains, which satisfies all the entry corridor constraints and meets the requirements with high precision. Afterwards, by utilizing special features of crossrange parameter, establishing bank-reversal corridor, and determining bank-reversals according to dynamically adjusted method, the algorithm enables the lateral entry guidance system to fly a wide range of missions and provides reliable and good performance in the presence of significant aerodynamic modeling uncertainty. Fast trajectory guidance profiles and simulations with a reusable launch vehicle model for various missions and aerodynamic uncertain-ties are presented to demonstrate the capacity and reliability of this method.     

二维翼段颤振的μ控制

 采用超声电机作为作动器来实现含控制面的翼段颤振鲁棒抑制。针对作者设计的二维翼段颤振主动抑制系统，通过理论与实验相结合的方法，建立了考虑沉浮方向阻尼和作动器模型参数不确定性的控制系统模型，设计了μ控制器，并对控制器做了降阶处理。数值仿真和风洞试验表明，μ控制器可有效地抑制颤振的发生，将颤振临界速度提高23.4%。相对于H_∞控制器，μ控制器的控制效果和鲁棒性更好。     

The effect of inlet conditions on the flow and heat transfer in multiple rotating cavity with axial throughflow

This paper discusses experimental results from two different build configurations of a heated multiple rotating cavity test rig.Measurements of heat transfer from the discs and tangential velocities are presented.The test rig is a 70% full scale version of a high pressure compressor stack of an axial gas turbine engine.Of particular interest are the internal cylindrical cavities formed by adjacent discs and the interaction of these with a central axial throughflow of cooling air.Tests were carried out for a range of non-dimensional parameters representative of high pressure compressor internal air system flows(Re up to 5×106 and Rez up to 2×105).Two different builds have been tested.The most significant difference between these two build configurations is the size of the annular gap between the(non-rotating) drive shaft and the bores of the discs.The heat transfer data were obtained from thermocouple measurements of surface temperature and a conduction solution method.The velocity measurements were made using a two component,LDA system.The heat transfer results from the discs show differences between the two builds.This is attributed to the wider annular gap allowing more of the throughflow to penetrate into the cavity.There are also significant differences between the radial distributions of tangential velocity in the two builds of the test rig.For the narrow annular gap,there is an increase of non-dimensional tangential velocity V/Ωr with radial location to solid body rotation V/Ωr=1.For the wider annular gap,the non-dimensional velocities show a decrease with radial location to solid body rotation.      

Abnormal Shape Mould Winding

为解决网格化芯模的缠绕问题,本文提出了复合材料面片缠绕机理;接着详细分析了面片缠绕过程中的芯模凹曲面上纤维滑线和架空现象,应用微分几何曲面理论和空间几何理论,提出判据及其解决方案;最后,针对飞机发动机进气道的缠绕成型,编制缠绕控制程序并进行相应的实验,验证了面片缠绕方法的可行性。     

基于H∞性能指标的质量矩拦截弹鲁棒控制

 以所建立的质量矩拦截弹数学模型为基础，利用微分几何的反馈线性化理论，得到一个解耦线性系统，考虑到拦截弹的鲁棒性要求和3个滑块的协调控制问题，提出采用双回路的设计方法，内回路采用线性二次调节器（LQR），外回路采用考虑混合灵敏度问题的H_∞控制设计。仿真结果证明了该方法动态性能满足设计要求，同时对系统参数的不确定性具有较强的鲁棒性。     

航空发动机管路测量数据分割方法

 为了提高航空发动机管路测量数据的反求建模效率,提出了一种区域增长分割算法。该方法主要是利用管路表面在曲面索引系数映射、主曲率映射、高斯映射上的特性,并基于均值漂移算法和遗传算法提出了区域分割算法中种子区域的选择策略。然后利用区域增长分割方法实现了对管路测量数据的分割。经过仿真和实测数据验证,所提管路分割算法具有较好的分割质量和效率。     

航天器返回地球的气动特性综述

航天器返回地球的飞行过程中,气动特性是实现将宇宙飞行速度减到落地前速度、保证再入飞行得到有效控制以及再入防热安全可靠的关键因素。针对简单旋成体气动外形、半弹道式再入控制、烧蚀防热类返回航天器,综述了返回地球过程中变化的空气流域特性、航天器周围的气体绕流环境、空气与航天器作用产生的动力学与热效应等。系统地给出了该类航天器的再入气动特性参数与飞行性能的共性规律,包括:气动阻力与再入减速、气动升力与再入轨迹控制、配平攻角与飞行稳定性、气动加热与防热,以及再入过程中不同气动特性航天器、气象条件变化等对再入飞行性能的影响规律。为航天器开展返回飞行过程的跨流域气动性能工程研制提供设计参考。     

Interfacial Properties Modification of Carbon Fiber/ Polyarylacetylene Composites

This work was dedicated to performing surface oxidation and coating treatments on carbon fibers (CF) and investigating the changes of fiber surface properties after these treatments, including surface composition, relative volume of functional groups, and sur- face topography with X-ray photoelectron spectroscopy (XPS) and atom force microscopy (AFM) technology. The results show that, after oxidation treatments, interfacial properties between CF and non-polar polyarylacetylene (PAA) resin are remarkably modified by removing weak surface layers and increasing fiber surface roughness. Coating treatment by high char phenolic resin solution after oxida- tion makes interface of CF/PAA composites to be upgraded and the interfacial properties further bettered.