基于干扰距离控制的直升机桨-涡干扰噪声降噪方法

提出了一种适用于直升机旋翼桨-涡干扰距离计算的方法。针对桨-涡干扰距离计算中必需的旋翼桨盘入流及桨盘倾角,又分别建立了旋翼耦合N-S （Navier-Stokes）/自由尾迹模型和桨盘配平模型进行求解,其特点是兼顾计算精度和计算效率;并从悬停、前飞两方面验证了方法的有效性。应用所建立的桨-涡干扰距离计算方法,着重分析了直升机纵向平面力以及飞行参数对旋翼桨-涡干扰距离的影响机理,并进一步计算了不同纵向平面力作用下的旋翼地面声场特性及影响规律。结果表明:合理施加纵向平面力、加速度等措施后,地面最大噪声降低约4dB。同时,还通过研究桨-涡干扰距离与飞行参数的影响关系,得到了一些对直升机低噪声飞行具有实际指导意义的结论。      

基于振动时滞法的非零叶间相位角叶片颤振

发展了一种利用叶片延迟振动设置叶间相位角的振动时滞法和多通道叶片非同相振动的流固耦合颤振分析模型。模型通道数选取相邻节径线之间通道数的两倍,在循环扇区的不同通道中,令叶片的各阶振动模态位移滞后于前一叶片,结合基于虚拟弹性体的快速动网格算法实现流场及叶片网格的高效更新。针对Rotor 37多通道模型,研究了不同叶间相位角对叶片气弹稳定性及通道流场特性的影响。结果表明:多通道方法与全环叶片颤振分析的计算结果基本一致,而18节径振动下多通道方法的计算时间是全环分析的1486%;节径振动形式对气动阻尼有显著影响,且在2节径时发生气弹失稳;叶间相位角引起流道内激波位置和强度变化和非定常激波脉动异相冲击,是影响颤振的主要原因。     

Parachute dynamics and perturbation analysis of precision airdrop system

To analyze the parachute dynamics and stability characteristics of precision airdrop system,the fluid–structure interaction(FSI) dynamics coupling with the flight trajectory of a parachute–payload system is comprehensively predicted by numerical methods.The inflation behavior of a disk-gap-band parachute is specifically investigated using the arbitrary Lagrangian–Euler(ALE) penalty coupling method.With the available aerodynamic data obtained from the FSI simulation,a nine-degree-of-freedom(9DOF) dynamic model of a parachute–payload system is built and solved to simulate the descent trajectory of the multi-body dynamic system.Finally,a linear five-degree-of-freedom(5DOF) dynamic model is developed,the perturbation characteristics and the motion laws of the parachute and payload under a wind gust are analyzed by the linearization method and verified by a comparison with flight test data.The results of airdrop test demonstrate that our method can be further applied to the guidance and control of precision airdrop systems.     

含激波的湍流流动高精度大涡数值模拟方法

针对采用亚格子模型进行含激波的湍流流动模拟时会面临激波附近的精度损失问题,考虑从通过亚格子模型以及数值模拟方法两方面的改进来实现湍流流动大涡模拟的精度提高.大涡模拟采用了Yee及Sj(o)green (2009)提出的高阶低耗散方法.该方法采用自适应的流场探测器以控制计算中所需区域的数值耗散,并考虑对动力学模型采用在激波位置使用Sagaut和Germano(2005)提出的单边亚格子过滤器和(或)直接禁用亚格子项等方法加以改进.对于标准的马赫数1.5和3条件下的激波-湍流干扰问题,上述新方法相较于全区域采用亚格子模型的方法均表现出了相似的精度提升.同时实现的数值精度改进方案采用了Harten的亚单元分辨过程来定位和锐化激波,并在精确激波位置附近的网格点处采用了单边测试滤波.     

共轴刚性双旋翼/机身干扰流场数值模拟

采用滑移网格技术求解Navier-Stokes (RANS)方程的方法,研究了共轴刚性双旋翼/机身的干扰问题。通过Caradonna-Tung旋翼、Robin直升机、Maryland直升机旋翼/机身干扰和Harrington 2共轴双旋翼等算例,验证了所提出的旋翼流场数值模拟方法的正确性。在此基础上,以Maryland机身为原型,分析了不同桨距的共轴刚性双旋翼与机身之间的干扰特性。结果表明:所提出的数值模拟方法能够很好地模拟共轴刚性双旋翼/机身的气动干扰特性;由于机身对于共轴刚性双旋翼下洗流场的阻滞作用,旋翼的悬停效率增加5%左右,并且随着拉力系数的增大使得悬停效率的增量更加明显;旋翼的悬停效率增加主要来源于下旋翼0°方位角附近的桨叶升力系数的增大,并且拉力系数的增量由桨根向桨尖方向逐渐减小。      

钝体稳定射流火焰的火焰面模型数值研究

基于稳态火焰面(SLFM)和交互式非稳态欧拉颗粒火焰面(EPFM)模型对Sydney大学CH4/H2钝体稳定湍流扩散火焰行了数值研究,采用修正的雷诺应力模型(RSM),同时对两种不同规模的甲烷详细化学反应动力学机理进行研究,比较了燃烧模型和反应机理对湍流火焰结构、活性自由基以及氮氧化物预测精度的影响,与实验数据对比结果表明:两种反应机理得到的温度场和主要组分分布基本相同;SLFM模型能对速度场和标量场的分布进行较为准确预测,采用EPFM模型修正后,部分区域OH预测结果更加靠近实验结果;采用EPFM模型对SLFM模型耦合GRI-Mech 211的计算结果修正后,NO量级降低近2倍,预测精度明显改善,与实验结果实现较好的符合,验证了化学反应动力学机理以及非稳态效应对氮氧化物预测精度的影响。     

轴向间距对涡轮时序效应影响的数值研究

时序效应可以改进多级叶轮机械的效率,为了探讨轴向间距对时序效应的影响,并应用到叶轮机械气动设计中,对1.5级静/动/静布局的涡轮叶片进行了3维非定常数值模拟。结果表明:保持轴向间距不改变,转子居中时涡轮效率整体最高,时均效率最大值为0.866,转子前移时次之,时均效率最大值为0.861,转子后移时效率最低,时均效率最大值仅为0.859,且转子相对位置发生改变并不影响效率-时序位置曲线;总轴向间距变大,时序效应变弱,效率提升幅度仅为0.7个百分点,效率-时序位置曲线随之改变。     

基于响应面法的短距/垂直起降飞机近地面升力损失

建立了短距/垂直起降(S/VTOL)飞机近地面升力损失的流场计算模型.通过数值模拟得出特定升力布局的飞机近地面状态各工况的升力损失.采用响应面法获得了飞机升力损失关于喷管落压比（NPR）、来流速度及飞机高度的2阶响应曲面函数及显著影响飞机升力损失的关键因素.并分析了喷管落压比、来流速度及飞机高度对飞机升力损失的交互影响作用,优化得出给定工况范围内升力损失最小的工作点.研究表明:仅考虑单因素影响时,升力损失随高度、落压比的增大而减小,随来流速度的增大而增大;考虑两因素交互作用时,高度与落压比及来流速度与落压比对升力损失存在交互影响,而高度与来流速度对升力损失无交互影响;优化获得的升力损失最小的工作点是飞机距地面高度为9D（D为喷管直径）、喷飞机高度为3、来流速度为0m/s,此时的升力损失为1.3%.      

液体火箭发动机充液导管流固耦合动力学特性

为了深入研究流固耦合（FSI）作用对液体火箭发动机充液导管频率特性的影响,采用传递矩阵法(TMM)建立了空间导管流固耦合动力学计算模型.针对真实发动机导管开展了传递矩阵模型与传统基于附加质量的有限元(FEM)（非耦合）模型仿真计算以及模态试验验证,比较了管径、壁厚等结构参数对导管流固耦合作用的影响.结果表明:在流固耦合作用下,导管各阶谐振频率减小、而对应的流体振荡与结构振动幅值增大.管径对导管低阶频率特性的影响较壁厚对其影响更大.对于该算例,当管径大于设计值30%后,耦合作用引起的1阶频率误差高于10%,此时流固耦合不能忽略;而壁厚对1阶谐振频率的影响则小于8%.      

适用于倾转旋翼机气动分析的新型嵌套网格方法

为克服围绕复杂倾转旋翼机生成高质量一体化结构网格的困难,将分块网格策略与嵌套网格方法融合,提出并建立了一套适用于倾转旋翼机气动干扰分析的新型嵌套网格生成方法。该方法中,将整体网格按照倾转旋翼机结构特性进行了分块生成,在保证插值精度的基础上降低了网格规模。为高效连接不同部件网格,首先,引入Inverse Map方法来辅助部件网格定位,采用新颖的多向投影方法保证了部件间洞边界的连续性;其次,基于独立挖洞嵌套策略完成机身与背景网格的组合,并使用了相应的辅助计算坐标进行背景边界网格标识。在贴体网格区域采用基于S A（Spalart Allmaras）湍流模型可压缩RANS(雷诺平均Navier Stokes）方程求解流场,背景网格区域采用欧拉方程,并运用了SPMD(single program, multiple data)模式的并行加速技术,建立了适合倾转旋翼机气动干扰特性分析的高效混合CFD方法。采用Robin直升机机身作为数值算例,验证了CFD方法的有效性。在此基础上着重对倾转旋翼机在悬停、过渡和前飞巡航模式下的干扰流场进行了数值分析研究,分别得出了悬停状态下倾转旋翼机典型的“喷泉效应”干扰现象、过渡状态下的非定常气动干扰特性以及巡航状态下的部件气动变化特性,计算结果表明建立的新型嵌套网格方法能够较好地表征倾转旋翼机外形特性,并能够有效地用于倾转旋翼机的气动特性分析,加速比能超过5.0。