脉冲横向射流的大涡模拟

为了研究脉冲气膜冷却的流动特性,采用基于高阶对称加权本质无振荡(WENO)-Roe格式和隐式牛顿迭代时间推进、动态亚格子模型和预处理技术下的大涡模拟(LES)计算方法,对吹气比为1.0下的斜管横向射流进行了稳态和脉冲两种不同工况下的数值模拟.计算结果显示了射流流场中存在发卡涡结构,研究表明射流孔进口位置的脉冲改变了射流出口拟序结构,即改变了射流与主流的掺混过程.      

可压缩平板边界层转捩的大涡模拟

入口利用层流速度型叠加二维和三维Tollmien-Schliting(T-S)波,应用动态亚格子模型、预处理技术以及高精度算法,采用大涡模拟(Large eddy simulation,LES)数值计算方法,对可压缩平板边界层空间模式的转捩机理及其拟序结构进行了研究。计算结果展示了T-S波经历线性、非线性发展阶段破碎形成发卡涡,最终演化成具有4个涡环的涡链结构的整个过程,并揭示了发卡涡涡腿并非完整对称及低速条带在涡环产生的位置不连续的现象,探讨了主流向涡与壁面相互作用诱导生成次流向涡的机理。     

基于代理模型的制导火箭炮发射诸元计算方法

针对制导火箭炮发射诸元的快速计算问题,提出了一种结合大样本数据和代理模型计算发射诸元的新方法。运用代理模型建立射角、无控弹道侧偏与炮位纬度、炮位高程、射向、射程、目标点高程及药温之间的函数关系,并根据射程和无控弹道侧偏的预测值对射向进行修正。仿真结果表明,高阶多项式响应面、相关函数为高斯函数的Kriging、高阶单项式径向基函数、核函数为高斯函数的最小二乘支持向量机、激活函数为正弦函数的超限学习机以及由上述单一代理模型构建的组合代理模型均具有较高的预测精度,各种单一代理模型对射角和无控弹道侧偏的预测时间均小于1 ms,证明了基于代理模型的射角和无控弹道侧偏预测方法切实可行,且通过对射向进行修正有效减小了由于地球自转引起的无控弹道侧偏。      

结冰模拟软件AERO-ICE中的关键数值方法

自然积冰对于航空飞行安全造成重大隐患,飞机在穿越富含过冷水汽的云层时冰形将按一定的物理规律积聚生长。介绍了一款三维结冰数值模拟软件AERO-ICE,该软件由网格自动生成、空气流场RANS计算、水滴场欧拉方法计算、结冰热力学分析四个模块组成。在空气流场计算方面,采用SPF k-v²-ω湍流模型,该模型引入湍流非平衡特性修正,预测的带冰翼型最大升力系数和失速攻角相对SA和SST模型有显著的提高。水滴场欧拉方程由于源项较大,迭代求解时容易发散,AERO-ICE软件采用流场光顺、二阶MUSCL空间离散以及LU-SGS隐式时间推进方法改善了数值稳定性。在结冰热力学分析模块,AERO-ICE软件同时具有Messinger与Myers模型,并将Messinger模型预估的壁面温度作为Myers模型的边界条件,从而解决了Myers模型温度设置的经验性问题。AERO-ICE软件支持多块网格、多重网格加速技术与大规模并行计算,其冰形计算结果得到了初步的验证。     

过冷大水滴变形与破碎的影响因素

过冷大水滴（SLD）结冰过程中,影响水滴变形与破碎及其导致的效应的因素众多,多种影响因素耦合作用于水滴变形与破碎的规律还未得到清晰认识。采用数值计算的手段,模拟了SLD运动过程中的变形与破碎现象,分别研究了水滴直径、流场速度及绕流物体尺寸变化对水滴破碎与变形及其效应的影响。通过研究,提出了包含SLD的典型流场区域划分方法,得到了不同因素对水滴破碎区域面积、物面的水滴收集特性、破碎子水滴分布特性的影响规律,发现水滴直径、流场速度、翼型前缘半径与变形和破碎影响撞击特性的程度呈正相关关系。在此基础上,提出了一个综合影响参数X,线性回归分析的结果初步说明,X是表征水滴变形与破碎及其效应的合适参数。相关工作可以为SLD结冰的计算、试验以及防/除冰设计提供参考。     

大型工程结构的有限元动力模型优化法

研究以振动试验测得的飞机结构固有频率为目标值的有限元动力模型的优化问题。应用了一种单元组合法,该方法基于动力平衡方程以及正交性。通过修改质量矩阵的元素来优化,使计算出的固有频率与试验测得的频率值基本一致,并成功地应用于无人机水平尾翼模型上,效果良好     

基于非Favre过滤的动力学非线性亚格子模型

采用直接过滤的Navier-Stocks（NS）方程组作为可压缩湍流大涡模拟控制方程组。方程组中因过滤产生的高阶相关项用Taylor级数展开近似表达,并仅保留级数的一阶导数项。因为截断级数产生的误差相当于丢失了模型的部分耗散作用,采用一种动力学非线性亚格子模型来补偿丢失的耗散影响。用模型模拟了超音速湍流平板边界层及绕双椭球高超声速湍流,并将所得结果与经验公式和计算结果进行了对比。对比表明模型可以有效的模拟可压缩湍流流场。     

超声速燃烧气动斜坡喷注器研究

为增强超声速燃烧过程的燃料混合,设计了通过喷嘴阵列形成流向涡强化混合的气动斜坡实验件.在直连式实验台上对气动斜坡实验件进行了纹影及油流谱等冷试实验研究.并采用大涡模拟(LES)方法数值模拟了超声速流场中气动斜坡的流动及燃烧特性.结果表明:气动斜坡喷注器的喷嘴阵列设计能有效加速喷流后流向涡的卷起,提高掺混特性;其燃烧过程也能使喷流形成的剪切层上抬,增进混合;由于喷流方向与主流方向夹角较小,总压损失也较物理斜坡、横向喷流等混合增强方式小.      

轴流风扇自噪声强度特性的数值研究

利用流场/声场混合模型计算分析轴流风扇自噪声的强度特性。使用大涡模拟(LES)数值计算方法获取湍流流场信息,以此为输入,使用声类比方法噪声预测模型获得了不同来流条件下的自噪声频谱。0°攻角条件下,叶片噪声源主要位于叶片尾缘区域,并且叶背上的噪声强度大于叶盆;攻角增大和速度增大都会引起叶片自噪声强度的增大;攻角增大引起的噪声增大主要集中在低频区域。     

涡破裂诱导的垂尾抖振数值模拟

 采用数值模拟方法研究了大后掠三角翼前缘涡破裂诱导的垂尾抖振问题,分析了大迎角条件下的垂尾抖振特性。采用Navier-Stokes方程求解非定常气动力、耦合结构动力学方程,建立了气动弹性方程,在时域内采用松耦合方式推进以得到垂尾结构响应。研究结果表明:涡破裂流的脉动频带覆盖了垂尾扭转模态的固有频率,诱发了垂尾抖振现象;与传统的颤振频域响应特性不同,垂尾抖振响应的各阶位移与加速度响应主频均位于各阶结构模态固有频率附近。此外,弯曲与扭转响应存在耦合效应,且耦合作用的频率与提取的垂尾表面气动载荷脉动频率一致。垂尾的位移响应由一阶弯曲模态主导,振幅不大;加速度响应主要由扭转模态产生,量级较大,使结构持续遭受严重的附加惯性载荷作用。