基于流固耦合的激波/边界层干扰作用下壁板颤振特性

为了解激波/边界层干扰作用下壁板气动弹性及其对流动分离的影响，采用自主开发的双向流固耦合求解器，对不同激波冲击位置下壁板的振动响应和流动特性进行了数值模拟研究。壁板几何非线性运动方程采用有限差分法求解，基于有限体积法求解Navier-Stokes方程组，对流通量采用MUSCL和AUSMPW+格式离散，双向流固耦合采用交错迭代算法。研究结果表明：激波/边界层干扰作用下壁板振动位移先增大后减小，经若干振荡周期后达到稳定颤振状态，呈现二阶振动模态，壁板变形相对于激波冲击位置呈现非对称性，壁板前部分的振幅始终小于壁板后部分；激波冲击位置可显著改变壁板的颤振振幅、频率及分离区长度，当激波冲击位置靠近壁板两端时，壁板振动最终收敛达到静稳定状态；壁板振动响应与流场特征不随激波冲击位置的改变而单调变化，对于激波冲击位置x/a=0.35工况，壁板颤振可有效抑制激波/边界层干扰流动分离。     

喷水射流预冷对加力燃烧室特性影响分析

为了探究采用射流预冷技术之后加力燃烧室性能，开展了不同喷嘴布置方案、喷水量和来流温度对预冷效果的影响研究。对射流预冷发动机工作过程进行了简化，建立了加力燃烧室进口前段射流预冷喷水特性计算的数学模型。同时搭建了小型试验台，通过与试验结果的比对验证了该模型的准确性，并利用该模型对射流预冷效果进行了仿真预测。结果表明：提高喷嘴数量与布置均匀性能够小幅度改善预冷效果；当来流温度不变时，射流预冷喷射腔室出口处的液态水蒸发量随着喷水量的增加而提高，但蒸发率却处于下降的趋势；当喷水量达到2%时，加力燃烧室燃烧效率对比不喷水工况会有一定的提升；喷水量达到4%以后，加力燃烧室出口温度及燃烧效率随着喷水量的提高而降低；喷水量大于8%以后，恶化了加力燃烧室（V型火焰稳定器）贫油熄火极限与燃烧效率；喷水量达到最大10%时，油气比需从原来设计工况的0.052上升到0.064才能保持稳定点火且对比不喷水时工况，加力燃烧室出口温度由1860K下降到1373K，燃烧效率由80.2%下降到69.2%。     

基于双滤波器的高精度目标跟踪估计方法

针对连续推力的合作航天器，采用双重无迹卡尔曼滤波(DUKF)算法估计其状态和加速度。通过状态滤波器和参数滤波器的配合，提升滤波精度，完成运动状态和参数的估计，从而实现合作目标的运动轨迹跟踪。与合作航天器相比，非合作航天器存在大小未知、发生时刻未知的机动，无法获得加速度，且信息获取和运动状态的估计难度大。针对非合作航天器，基于简化的相对运动方程，结合天基平台获得目标的观测信息，采用两个扩展卡尔曼滤波(EKF)及基于半正焦弦的机动检测策略实现多未知脉冲机动的运动状态的估计。仿真结果表明：相比于无迹卡尔曼滤波（UKF），DUKF在对合作航天器的状态和加速度估计方面具有更快的收敛速度和更高的滤波精度；对于存在未知机动的非合作航天器，通过对比验证机动检测策略与滤波器切换策略相结合的方法的有效性，该方法能够检测到多次机动并且减少误判。     

前体尾流对降落伞工作性能的影响

为研究前体尾流对降落伞工作性能的非定常影响，基于Realizablek-ε湍流模型采用PISO算法开展了物伞系统的非定常绕流数值计算，获得了精细的流场旋涡结构。在此基础上，研究了不同拖曳比下物伞系统的尾涡演变规律、流场分布规律以及伞衣气动特性变化。结果表明:前体尾涡导致伞衣入口处的涡量大小和方向时刻变化，随拖曳比增加，涡量黏性耗散增强，进入伞衣的旋涡强度逐渐减弱，伞衣入口形成稳定的负涡量区，伞衣尾涡脱离周期随之延长；拖曳比对尾涡区后端（伞衣入口处）流场压力的影响远大于前端，随拖曳比增加，流动形式逐渐由闭式转变为开式，流场的速度分布和压力分布更为对称，伞衣入口形成稳定的正压区，内外压差增加；当拖曳比大于9时，前体尾流对降落伞阻力系数和表面压强系数的影响减小。      

MIMO窄带加宽带非高斯随机振动试验

提出了一种用于多输入多输出窄带加宽带非高斯随机振动试验方法。首先，详细分析了随机相位与非高斯随机信号偏度和峭度的关系，提出了一种可用于快速生成具有指定偏度和峭度的迭代相位调节法。其次，根据控制通道之间的耦合关系，利用各通道之间随机相位的不相关特性，实现了多通道非高斯随机信号生成的解耦。然后，将生成的非高斯随机信号作为参考输出响应，采用时域逆系统方法生成驱动信号，并利用控制算法分别对响应信号的功率谱、偏度和峭度进行控制。最后，通过三轴向振动台试验验证了本文提出的多输入多输出窄带加宽带非高斯随机振动试验方法的可行性。结果表明，控制点上加速度的功率谱被控制在了±3 dB的容差限内，满足预设的容差要求。同时，时域响应信号的峭度和偏度也满足给定的参考要求。     

具有复杂围护结构光学窗口的传热特性分析

针对具有复杂围护结构光学窗口传热特性,建立了高热流作用下的光学窗口及其复杂围护结构的辐射-导热传热模型,分别讨论了窗口的辐射半透明、光谱选择性、外部加热热流以及窗口尺寸对其瞬态温度响应的影响。结果表明,忽略窗口半透明特征会产生最大31.4%误差,忽略光谱选择性会产生最大40.4%误差;窗口外部气动加热热流影响显著,温度响应随热流的增大而急剧增大;窗口厚度对其传热特性的影响较大,随着厚度的增大窗口温度响应减小;而窗口半径产生的影响可忽略不计。     

翼尖铰接复合飞行器动力学特性研究

航空飞行器通过翼尖铰接机构复合飞行时的气动耦合效应，会造成飞机产生不同于其单独飞行时的动力学特性，出现复合飞行安全问题。为研究翼尖铰接复合飞行器的动力学特性，使用Newton-Euler方法和Robberson-Wittenburg方法建立了双机组成的翼尖铰接复合飞行器多刚体系统整体和内部的7自由度非线性动力学和运动学方程组。在气动准定常假设下建立双机复合系统非耦合气动力表达式，基于CFD方法开展复合飞行器系统的三维实体建模和非结构网格划分工作，获取复合飞行器系统的气动力数据。搭建动力学仿真平台，开展准配平方案和全配平方案下的动力学仿真。仿真结果表明:准配平方案下飞行器无法持续稳定飞行，而全配平方案下复合飞行器系统各运动参数在仿真时间内始终趋于稳定。在全配平方案下，使用小扰动假设的非解耦线性化方法重新整理7自由度动力学方程组，研究复合飞行器系统运动模态的特征值中出现的2个发散新模态，根据对应的特征向量发现2个发散模态分别由相对滚转角度和角速度主导，同时也比较分析了其他运动模态相比单机飞行时的特性变化规律。      

R141b在圆形微通道内的沸腾换热实验研究

为提高换热强度、解决设备内部高热流密度散热问题,采用实验方法研究R141b在不同直径(D=0.5mm和1.0mm)水平圆形微通道内的沸腾换热特性,分析了热流密度(q=2.0kW/m~2～47.6kW/m~2)、质量干度(x=0～0.6)、质量流速(G=111.11kg/(m~2·s)～333.33kg/(m~2·s))的变化对平均传热系数h的影响,探究不同情况下影响沸腾换热的主导因素。实验研究表明:平均传热系数h随热流密度q的增加而减小,在不同范围内减小速率有明显差异;热流密度q=2kW/m~2～5kW/m~2时质量流速G对平均传热系数h影响较明显,热流密度较高时质量流速G对换热影响很小;在质量流速G=111.11kg/(m~2·s)～333.33kg/(m~2·s),质量干度x0.3时,平均传热系数h随质量干度x增加而明显下降,在设计微通道换热器时应尽量使R141b处于初始沸腾阶段以获得更好换热效果,并采取一定措施预防干度过高引起的换热恶化。     

非均匀间距的低副瓣宽带微带阵列天线设计

为降低天线副瓣电平（SLL）和展宽带宽，设计了一款谐振频率为14.25 GHz的16阵元非均匀间距的耦合馈电微带阵列天线。天线采用多层设计，通过在接地板开矩形槽进行耦合馈电，并引入空气层，降低天线Q值，增大带宽。区别于均匀间距阵列天线的激励幅值加权，从阵元间距角度入手，利用差分进化算法降低副瓣电平，构建非均匀间距并联线阵天线。用槽面辐射的能量近似代替阵元接收的能量，观察阵元功率分配情况，并建立馈电网络所有馈线段的数学关系，保证非均匀间距条件下所有阵元为等幅同相激励。测试结果显示，天线在14~14.5 GHz范围内电压驻波比小于2，满足了卫星动中通的带宽要求；工作带宽内增益大于16 dB，副瓣电平低于-16 dB，性能优于均匀间距阵列天线。      

偏心撞击对撞击式喷嘴雾化特性的影响

为研究偏心撞击对撞击式喷嘴雾化特性的影响,建立了求解自燃推进剂冷态射流撞击雾化过程的数值模拟方案,计算了不同偏心度条件下的射流撞击雾化过程。采用树形自适应加密算法直接求解不可压Navier-Stokes方程组,由分段线性的流体体积(VOF)方法对流体界面进行捕捉。结果表明偏心撞击会导致雾场发生偏转,当无量纲偏心度E为1/8时,雾场偏转角度约为9.2°,应控制加工偏差小于该值。随着偏心度的增大,液膜的偏转角度增大,理论推导得到的液膜偏转角度要小于数值计算得到的液膜偏转角度。正心撞击时燃料与氧化剂流强峰值接近,雾场的流强分布呈单峰分布。当发生偏心撞击时,由于燃料与氧化剂部分射流未参与撞击导致流强峰值出现交错,雾场的流强分布呈双峰分布,混合比的空间分布发生较大改变。正心撞击时撞击点下游液滴的速度分布近似呈轴对称分布,而偏心撞击之后的速度分布则呈中心对称分布。偏心撞击导致的射流动量损失使得雾化性能变差,当无量纲偏心度E为1/8时,一甲基肼(MMH)的Sauter平均直径增大约4.8%,四氧化二氮(NTO)的Sauter平均直径增大约5.8%。