超声速气流中液体横向射流的气液相互作用过程数值研究

采用Euler-Lagrange方法对来流马赫数为1.94的超声速气流中液体横向射流的气液相互作用过程进行数值研究。计算给出的射流穿透深度、液滴Sauter平均直径(SMD)及液滴速度分布均与实验吻合较好。仿真结果较详细地揭示液体射流喷雾与气流之间的强烈相互作用过程。受液雾影响,射流前形成较强激波,气流依次经过激波及液雾区域,气流速度存在两次下降过程。计算结果揭示,超声速来流可以与射流的液滴轨迹相交,气流经液雾前沿进入液雾区域后,流向往壁面偏折。本文首次发现并提出,由于气液相互作用诱导形成两组反向反转漩涡对,这对于理解两相混合过程具有重要意义。气液相对滑移速度的分析表明,液滴在穿透自由来流并开始转向时受气流作用最为显著,完成转向后气液相互作用逐渐减弱。     

Bump进气道中鼓包诱导的激波/边界层干扰特性

为了探索Bump进气道中鼓包诱导的锥形激波和机身发展而来的湍流边界层干扰问题,分析其气动优势,首先选取了半锥和半棱锥这两种与鼓包的流场结构具有一定相似性的构型作为参照,采用数值仿真方法,分别对这三类典型的三维激波/湍流边界层干扰问题进行了流场分析。在此基础之上,设计了三个不同马赫数的鼓包,并研究了设计马赫数对鼓包流场特性的影响。结果表明:当三类构型的无黏激波强度相等时,半锥诱导产生的旋涡强度最强,鼓包次之,半棱锥最弱。尽管鼓包诱导的流场非常复杂,其干扰流场却呈现出准锥形相似的特性。虽然半锥对边界层的排移能力最强,但是综合考虑边界层排移能力及进气道出口流场畸变下,鼓包最具优势,这也是其被选为超声速进气道前缘压缩面的重要原因之一。此外,在设计状态下,适当增加设计马赫数能改善鼓包排移边界层的能力,但设计马赫数太高,边界层排移能力基本不变,反而使得进气道总压损失急剧增加。      

压气机动静叶干扰的非定常流动特性

采用大涡模拟(large eddy simulation,LES)方法计算动叶尾迹对静叶干扰的流场信息。利用涡量分布揭示动叶尾迹在静叶通道内的演化过程,利用压力梯度识别激波结构及波振源,运用动力学模态分解(dynamic mode decomposition,DMD)方法对静叶通道内流场的时空结构进行模态分解。结果表明:流场中存在3处波振源,分别位于动叶尾缘、静叶前缘和静叶尾缘处;发现静叶通道内流场的频谱具有多峰值现象,模态分解的第1阶流动代表动叶尾迹在通道内随时间迁移,对应频率为动叶通过频率(blade passing frequency,BPF)是通道内旋涡非定常波动的主导频率;第2阶流动是动叶通过频率的2倍频流动,旋涡的空间尺度为1阶模态的1/2,为更小尺度的扰动。      

Unsteady Flow Variability Driven by Rotor-stator Interaction at Rotor Exit

Numerical investigation of the unsteady flow variability driven by rotorstator interaction in a transonic axial compressor is performed. Two models with close and far axial gap between rotor and stator rows are studied in the simulation. Particular attention is attached to the analysis of mechanisms involved in driving rotor wake oscillation, rotor wake skewing and flow angle fluctuation at rotor exit. The results show that smaller axial gap is favorable to enhance the interaction in the region between two adjacent rows, and the fluctuation of the static pressure difference between two sides of rotor wake is improved by potential field from down stator, which is the driving force for rotor wake oscillation. The interaction between rotor and stator is weakened by increasing axial distance, rotor wake shifts to suction side of rotor blade with 5%-10% of rotor pitch, the absolute value of flow angle at rotor exit is less than that in the case of close interspace for every time step, and the fluctuation amplitude is also decreased.     

多块协调变形的网格变形技术及其应用

基于多块网格协调变形的理论,提出了适用于叶轮机内流多物边界的更具一般性的网格变形技术,先获得块边界处的变形量,然后插值得到块内网格各点的变形量。通过二维振荡叶栅和带叶尖间隙的三维跨音风扇叶片的非定常绕流算例,证明网格变形技术在叶轮机气弹稳定性数值模拟中的应用是有效的,能够在保证网格高精度的条件下,提高网格变形速度。     

非定常环境下动叶气膜冷却流场的数值模拟

对非定常环境下动叶气膜冷却的流场进行了数值模拟,研究了动静部件之间非定常效应对气膜冷却效率的影响.计算结果表明,静叶的尾流迁移在冷气喷射位置处时,冷气在喷射位置下游壁面附近的覆盖宽度增大,距喷射位置下游约30%弧长范围内具有较高的气膜冷却效率,但在喷射位置附近的冷却效果略差一些.气膜冷却效率的振荡幅度随着吹风比的增大而增大.由于端壁附近旋涡的影响,冷气很难到达端壁附近,而是被卷吸到喷射位置的下游去,在喷射位置下游靠近端壁附近的冷却效果更好一些.      

月球车轮与月面相互作用的两种仿真模型的比较与应用

文章进行了月球车单个车轮和模拟月壤相互作用的初步研究,并将两种不同轮下分布模型得到的单轮力学性能仿真分析的结果和土槽试验结果进行了对比,从中选出合适的模型对车轮的牵引性能进行仿真计算。根据试验结果,利用动沉陷量-滑转率之间的线性关系获得了动沉陷系数,然后对不同垂直载荷、不同表面形状车轮牵引力的变化规律进行了预测,并进行了相应的试验验证。     

降落伞开伞过程的多结点模型仿真

根据降落伞的结构和力学特征,在轴对称假设下创建了伞衣及回收物系统的多结点结构模型.通过考虑应力,重力和气动力的作用效果,建立了用于无气流攻角平面圆形降落伞充气模拟的多结点结构模型动力学方程组.对开伞过程中的流场变化引入准定常假设,利用simple算法数值模拟求解RNG(Renormalization Group)k-ε湍流模型下的雷诺平均N-S(Navier-Stokes)方程以获得选定时刻的伞衣表面压力分布.结合多结点模型动力学方程组的解算代码和计算流体力学程序,采用流固耦合的方法对选定的平面圆形降落伞模型的开伞过程进行了动态仿真,得到了开伞过程中降落伞外形和特性的变化.通过结果分析和比较,证明了多结点模型的可行性,发展出了一种用于降落伞流固耦合计算的新方法.      

氢氧发动机真空羽流干扰试验研究

为了研究航天器发动机多股羽流产生的羽流干扰对流场结构及参数的影响,在北京航空航天大学真空羽流效应试验系统（PES）中,采用总压皮托管阵列对60 N氢氧模拟发动机单机羽流和双机干扰羽流的压力场进行了测量。试验结果表明:单机羽流流场压力随着轴向距离的增大而迅速降低;受到钟形喷管产生的压缩波束影响,同一轴向距离上的压力最大值逐渐偏离喷管轴线;双机羽流干扰发生于两股发动机羽流主流之间的区域,干扰流的作用范围和强度随着轴向距离的增大而增大,干扰流压力最高可以达到单机羽流相同位置处压力的5倍以上;受到干扰流压缩波束边界的影响,同一轴向距离上的压力最大值逐渐偏离双发动机对称面。      

一种空间目标RCS序列的姿态异常检测方法

准确、快速判断空间目标姿态运动模式异常,对于空间目标监测具有重要意义。针对空间目标雷达散射截面 (Radar Cross Section,RCS)序列,提出一种基于小波包分解（Wavelet Packet Decomposition,WPD）能量谱特征的无监督机器学习异常检测方法,并采用单类支持向量机(One Class Support Vector Machine,OCSVM)验证异常检测效果。设置了几种典型异常场景进行仿真分析,试验结果表明,该方法能有效检测出三轴稳定类空间目标发生失稳旋转的姿态异常。相比于传统统计参数特征、小波变换统计参数特征及能量特征的姿态判别方法,具有检测概率高、鲁棒性好的特点。