基于非结构重叠网格的螺旋桨滑流非定常数值模拟

基于各向异性非结构混合网格及多套网格重叠技术,通过求解非定常雷诺平均Navier-Stokes(URANS)方程,分别研究了螺旋桨动力特性和螺旋桨滑流对涡桨飞机气动特性影响。运用物面相交准则实现多套网格间的挖洞,采用距离权和三线性插值技术传递重叠区变量信息,并针对定轴旋转优化了网格重叠边界。采用双时间步求解控制方程,内迭代计算采用LU-SGS(lower-upper symmetric Gauss-Seidel)方法。模拟了单独螺旋桨旋转的非定常流动,计算的拉力系数和扭矩系数与试验结果一致,表明采用的数值方法和网格技术能有效模拟螺旋桨滑流效应。数值模拟了某涡桨飞机,对比了有无滑流下飞机的气动特性,给出了滑流对飞机影响。结果表明:滑流导致其后部气流加速和旋转,改变飞机气动特性,增大飞机的阻力和俯仰力矩,导致气流下洗,影响机翼及平尾压力分布。      

甲烷摩尔分数及初始压力对甲烷/RP-3航空煤油混合燃料燃烧特性的影响

采用定容燃烧实验装置获得初始温度为450K、初始压力为0.1~0.3MPa、当量比为0.7~1.5以及甲烷摩尔分数为0~0.8工况下甲烷/RP-3航空煤油混合燃料火焰发展特性图片、马克斯坦长度和层流燃烧速度等燃烧特性,分析甲烷摩尔分数及初始压力对甲烷/RP-3航空煤油混合燃料燃烧稳定性及层流燃烧速度的影响。结果表明:当量比为1.3时,随着甲烷摩尔分数增加,甲烷/RP-3航空煤油混合燃料燃烧趋于稳定,初始压力对燃烧稳定性影响较大,随着初始压力增加,燃烧稳定性变差。混合燃料马克斯坦长度随当量比增加而减小,当甲烷摩尔分数增加时,混合燃料马克斯坦长度减小趋势变缓,当初始压力增加时,混合燃料马克斯坦长度减小趋势明显变缓。混合燃料层流燃烧速度随当量比增加呈现先增大后减小的变化趋势。当甲烷摩尔分数为0、0.4和0.6时,随着甲烷摩尔分数增加,混合燃料层流燃烧速度逐渐增大,当初始压力为0.1、0.2、0.3MPa时,随着初始压力增加,混合燃料层流燃烧速度显著降低,随着甲烷摩尔分数和初始压力的增加,混合燃料层流燃烧速度峰值有向当量比大的区移动的趋势。      

涂敷硬涂层的失谐整体叶盘减缩建模及振动分析

为了提高整体叶盘的服役寿命,提出一种对叶片涂敷硬涂层的减振方法。同时,针对失谐叶盘-硬涂层复合结构的运算规模庞大的难题,提出了一种双重减缩建模方法。利用一种改进的混合界面子结构(HISCMS)法对复合结构进行第一重的自由度减缩。利用一种模态密集判别准则与经典模态减缩(SNM)法的思想对综合后的模型进行第二重的模态减缩。选取了对叶片双面涂敷NiCoCrAlY+YSZ(yttria-stabilized zirconia,氧化锆) 硬涂层的失谐叶盘进行仿真分析及试验测试,以探究失谐叶盘的振动情况以及硬涂层厚度对失谐叶盘的影响。结果表明:双重减缩建模方法能够在保证计算精度的前提下提高25.86%~42.05%的计算效率;而且硬涂层表现出较强的阻尼作用,能够在共振区显著抑制失谐叶盘的共振响应。      

基于流-固耦合的混合火箭发动机固体燃料表面退移速率计算

基于流-固耦合的方法,在充分考虑混合火箭发动机工作过程中诸多复杂物理过程的基础上,建立了一个可适用于不同工作状况下混合火箭发动机固体燃料表面退移速率预示的计算模型。计算结果与实验数据的对比验证了所建立计算模型的准确性。对模型发动机进行模拟的结果表明,混合火箭发动机中的燃烧、流动及固体燃料表面的退移速率具有明显的不均匀性,发动机中的固体燃料表面的退移速率沿轴向近似地呈“W”形状的曲线变化;在混合发动机中,突扩形状的预燃室和补燃室有利于燃料热解气体和氧化剂气体的扩散混合,可以强化对固体燃料表面的换热,提高固体燃为表面的退移速率。     

二维嵌套混合MIMO相控阵雷达接收阵列设计

针对二维混合多输入多输出（MIMO）相控阵雷达发射阵列子阵分割带来的自由度损失问题,提出一种基于二维嵌套阵列的二维混合MIMO相控阵雷达接收阵列设计新方法。首先,给出一种由稀疏阵列和密集阵列构成的嵌套阵接收阵列的二维混合MIMO相控阵雷达信号模型;其次,通过对接收信号的协方差矩阵进行Khatri-Rao乘积处理,得到阵元位置差的差异阵列,形成接收阵元数目的虚拟扩展;最后,通过空间平滑处理进行波达方向估计。仿真实验表明,与传统二维混合MIMO相控阵雷达相比,所提方法在不增加实际阵元数目的情况下可以有效扩展虚拟阵元数目,提高雷达阵列自由度,进而提高二维混合MIMO相控阵雷达波达方向估计精度。     

DDES方法在复杂旋翼流场计算中的应用

应用延迟脱体涡模拟（DDES）方法,进行了复杂分离流动旋翼流场的数值模拟研究。为控制计算规模,更好地模拟旋翼桨叶近壁面附近区域流动分离及远场尾迹发展,采用了贴体与背景自适应直角网格嵌套的网格系统;同时为提高计算速度,采用了基于分布式存储的数据并行模式,实现了贡献单元搜索、洞切割和流场计算的并行化。计算工作包含一个共轴双旋翼悬停和单旋翼下降状态,目的是分别讨论DDES方法在旋翼表面附近空间区域是否存在大范围分离时计算结果与雷诺平均Navier-Stokes结果之间的差异。计算表明,对于下降状态,DDES方法预测出了更强的尾迹区分离流动;对于桨叶表面基本为附着流动的双旋翼中等桨距角情况,两种方法预测的流场结构相似,仅在下旋翼桨叶气动力和旋翼桨根下方空间区域的诱导速度计算结果上存在微细差异。     

小分子燃料对RP-3航空煤油燃烧作用的数值研究

为研究小分子燃料对RP-3航空煤油燃烧的影响,选择合理的RP-3航空煤油替代燃料详细燃烧模型开展工作,该模型能够精确预测RP-3航空煤油和小分子燃料的燃烧特性。以六种重要的小分子燃料H2,CH4,C2H4,C2H6,C3H6和C3H8分别与RP-3航空煤油按1:5的比例(摩尔分数)掺混形成的六种混合燃料(Blended Fuel)为研究对象。在当量比为1.0,压力分别为0.1 MPa和1 MPa下系统模拟了RP-3航空煤油及六种混合燃料在高温下的燃烧特性,分析了各种混合燃料的自点火、燃尽时间、绝热火焰温度、熄火温度、组分浓度变化,并结合ROP(Rate of Production)分析方法,分析了小分子燃料对OH自由基生成速率的影响。结果表明,C2H4将RP-3航空煤油的点火延迟时间缩短了近4.6%;C3H6则将RP-3航空煤油的点火延迟时间推后了8.4%;C2H4和H2对RP-3航空煤油的快速点火和稳定燃烧有着积极的作用,其中C2H4的作用最为突出。     

基于造型系统的曲面及二维域有限元网格剖分

介绍了一个基于几何造型系统的有限元分析前置处理系统．其主要功能是对几何造型系统的ＮＵＲＢＳ曲面以及二维任意形体有限元网格划分．该系统可作为几何造型系统与有限元分析程序之间的接口．     

高速太阳风中α粒子分布的形成

以多级波粒相互作用模型和1维电磁混合模拟方法研究了高速太阳风中α粒子分布的形成机制．认为α粒子从其产生区到1AU,经历4个特征区域：在冕洞下方的小尺度磁重联区,α粒子产生并得到加速：在约十几个太阳半径日冕区,α粒子形成球壳分布函数：在几十个太阳半径处,在Alfven波和温度各向异性不稳定性的作用下,α粒子分布函数演化为漂移Maxwell分布：当α粒子进入0．3－1AU时,Alfven波的作用使α粒子速度超过太阳风速度约1个局域Alfven速度．     

机械臂的动态混合控制

本文研究当机械臂的终端受有约束时的控制问题,文中应用了动态混合控制方法。其中心内容是给出关于“任务规范投影算子”的概念。利用它首先将机械臂的动态方程解耦为两组方程,它们分别描述了运动与约束反力。在此基础上给出了机械臂的控制律,使闭环系统跟踪期望的速度与约束反力。