侧向控制喷流干扰流场特性数值研究

侧向喷流是导弹在中高空进行姿态控制的重要手段.本文用数值方法求解N-S方程,模拟了二种带侧向喷流控制的喷流与绕流的干扰流场.通过定性流场结构和物面压力分布与实验结果的比较,数值模拟能较好地模拟侧向喷流与外流干扰产生的复杂流场特性.给出了合理的流场结构和干扰气动力特性,表明数值方法可以作为型号侧向喷流控制设计的研究手段.     

一种结合壁面模型模拟湍流流动的间断界面浸入边界法

提出了一种具有鲁棒性的间断界面浸入边界法，用来模拟不可压湍流的物体绕流问题。开发了一种新的OpenFOAM求解器，该求解器结合了压力隐式分裂算法和k-ωSST湍流模型求解Navier-Stokes方程。为了减少对近壁面网格的要求，在实施边界条件时，采用壁面模型计算得到的壁面剪切力将速度和湍流变量修正联系起来。壁面切应力采用远离壁面的镜像点求得，镜像点上流场变量通过反距离线性插值由周围流场信息得到。使用该求解器模拟了平板绕流、NACA0012表面压强和切应力分布、Buice扩散管、圆柱绕流、方柱绕流5个湍流流动问题，本文数值模拟结果与文献结果有较好的一致性，验证了该方法模拟湍流流动的有效性和准确性。     

三维薄翼非对称绕流的跨音速积分方程及其有限元解法

与有限差分法相比,积分方程的有限元解法具有计算工作量小,对计算机容量要求低而又能达到一定的计算精度等优点,但国外的研究尚限于二维翼型。最近,文[4]发展了积分方程的有限元解法并成功地应用于求解三维薄翼跨音速对称绕流问题。在此基础上,本文进一步发展该方法并用于求解三维薄翼跨音速有升力的非对称绕流问题。     

用非结构动网格方法模拟有相对运动的多体绕流

发展了用于模拟有相对运动的多体非定常绕流，并计算由流场决定的外弹道的数值方法。该方法采用非结构动网格，显式NND有限体 格式，求解非定常Euler方程，并耦合了刚体动力学方程，作为例证，计算了俯仰振动的NACA0012翼型绕流，计算结果与实验及文献结果吻合，最后，整个策略用于计算二维机翼－外挂物分离流场，并确定了外挂物下落轨迹。     

欧拉方程多段翼型绕流计算的嵌套网格方法

本文采用嵌合体埋入技术，将主段翼型和襟翼分成两子域，对多段翼型的复杂构型进行了分区迭代求解，在有限体体积离散的基础上，求解欧拉方程计算各子域流场的解，通过各流场的耦合迭代到多段翼型绕流的解，典型的带３０％襟翼ＧＡ（Ｗ）－１翼型算例表明，该计算方法稳定，能较地模拟多段翼型中主段翼型与襟翼交接处的复杂流动现象计算结果与实验值有较好的一致性。     

基于非结构重叠网格的二维N—S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

亚声速大迎角模型试验洞壁干扰修正方法研究

通过测量洞壁附近的压力分布来模拟透气壁试验段的洞壁边界条件，采用数值求解Ｅｕｌｅｒ方程的方法模拟模型在风洞中的绕流场，然后将洞壁边界条件用远场边界条件替换，用同样方法算出模型在自由流中的绕流场，从而计算出洞壁干扰对模型的绕流场和气动力的影响。针对大迎角模型的洞壁干扰问题的特殊性，采用重叠网格技术，对飞机大迎角标模和ＳＢ－０３模型进行了亚声速大迎角的洞壁干扰计算与修正。     

三维方柱不可压缩绕流的大涡模拟计算

用非定常不可压缩流动求解方法对三维方柱绕流问题作了大涡模拟计算研究.求解程序采用虚拟压缩方法和双时间步的时间推进方案,对流项的差分采用高阶精度的数值差分格式.通过与实验和文献结果比较分析,验证了基于虚拟压缩方法的湍流大涡模拟计算方案对于复杂的三维钝体绕流计算是可行的.     

旋转坐标系N-S 方程无网格/网格混合算法

为了模拟涡桨发动机等绕固定轴旋转的桨叶流场,发展了一种求解旋转体非定常黏性绕流的无网格/网格混合算法。算法基于求解旋转坐标系纳维尔-斯托克斯方程展开,避免了旋转角速度特征物理量的插值运算;计算域采用整体网格和物面附近局部无网格离散,通过引入无网格和网格对偶点,实现了混合算法绕流信息的跨区交换,并给出了一种对偶点调整选取的方法;基于无网格点云空间导数逼近方法离散控制方程,并结合隐式LU-SGS算法,给出了混合算法求解非定常问题的双时间步隐式推进格式。结果表明:所提算法通过2维振荡翼型非定常绕流、旋转圆柱黏性绕流和模拟发动机桨叶旋转运动的3维悬停旋翼绕流进行了考核,所得升力系数等重要气动数据的变化趋势与文献试验或计算值一致,典型截面处所示的激波强度和位置亦与文献值吻合,展现出算法在模拟旋转部件绕流问题方面具有广阔的工程应用前景。     

隐式无网格算法及其应用研究

本文的主要目的在于研究求解Euler方程的隐式无网格算法，并应用于复杂的流场计算。采用无网格算法，计算区域用点云离散代替通常的网格划分；计算点上的空间导数，用当地点云上引入的二次极小曲面逼近。求解的Euler方程用隐式时间后差离散，结合用Roe的近似Riemann解确定通量，并用LU－SGS算法分步计算，数值求解了单翼型或双翼型模拟的复杂绕流。     

ν-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_∞控制设计。仿真结果证明了该方法动态性能满足设计要求，同时对系统参数的不确定性具有较强的鲁棒性。     

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

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

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

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

基于马赫数分布可控曲面外/内锥形基准流场的前体/进气道一体化设计

提出了一种高超声速飞行器乘波前体的外锥形基准流场设计方法,在锥面马赫数分布规律给定的条件下,通过有旋特征线法实现反设计,提高了基准流场设计的灵活性。该基准流场通过锥形"下凹"弯曲激波和波后等熵压缩波系压缩气流,可以在较短的长度内完成高效压缩。基于反正切马赫数分布外锥形基准流场设计的乘波前体具有较高的容积率,乘波特性良好且出口均匀,设计点时有黏升阻比为1.89。另外,基于该乘波前体和马赫数分布可控的内收缩进气道给出了一种双乘波的前体与进气道一体化设计方案,实现了内外流分别独立乘波,充分发挥了乘波前体和内收缩进气道的各自优势。