非稳态流体网络方法在发动机空气冷却系统中的应用

将流体网络法与非稳态模拟方法相结合应用到发动机空气冷却系统沿程温度、压力、流量的非稳态特性预测中.本方法将空气系统简化为节点和元件组成的网络,节点和元件处建立守恒方程,针对传统解法的计算稳定性问题引入了压力修正和积分方法相结合的求解方法.程序应用到可参考对照的模型中得到验证.验证后的程序被应用到发动机轴部冷却通道空气系统非稳态特性预测中.      

上游尾迹与涡轮转子泄漏流相互作用数值模拟

叶轮机内部流动本质上是周期性非定常的,研究涡轮转子叶尖区域的非定常相互作用机理,对提高小展弦比高负荷涡轮性能具有重要意义.利用数值模拟方法研究了上游静子尾迹与涡轮转子叶尖泄漏流的非定常相互作用,分析了定常结果、时间平均结果以及瞬时时刻结果的流动图画.结果表明:上游静子尾迹与涡轮转子尖区二次流的相互作用能明显影响泄漏涡和机匣通道涡的时空演化规律,从而改变转子尖区的损失分布.上游尾迹在转子通道中传播时,诱导泄漏涡和通道涡区域出现周期性的扰动涡对,扰动涡对沿着泄漏涡和通道涡的轨迹向下游运动,使得转子尖区二次流结构呈现周期性变化.      

不同轴向间距下涡轮级内非定常流场的数值研究

针对不同轴向间距下涡轮级内非定常流动进行了数值研究.讨论了非定常条件下轴向间距变化对非定常性及流场分布的影响,探讨了转子和静子叶片之间的非定常相干的机理.结果表明,轴向间距L对转子内流动的非定常性有着强烈的影响.随着轴向间距的减小,静叶尾迹的强度及影响范围增加,静叶尾迹对动叶的周期性非定常影响还是主要出现在前缘附近.由于静叶尾迹的恢复效应与来流尾迹强度成正比,50%轴向间距下动叶出口流场的湍动能在动叶通道中央区域最小.      

湍流边界层厚度对三维空腔流动的影响

采用脱体涡模拟(DES)方法开展了不同湍流边界层厚度(TTBL)下的三维空腔非定常流动数值计算。空腔长、宽、深比例为5:1:1,来流马赫数为0.85,雷诺数为13.47×10⁶ m^-1,各工况湍流边界层厚度比值为1:2:4:8。研究结果表明,湍流边界层厚度对自由剪切层的发展、空腔底部静态压力分布、脉动压力及空腔流动类型均有重要影响,且随着边界层厚度的增大,下游剪切层覆盖的范围会增大,但是剪切层增长率降低;空腔前后静态压力压差减小、压力梯度下降;腔内局部测点的脉动压力声压级下降,各阶声压峰值频率向低频方向偏移;空腔流动类型往开式流动方向转换。     

基于阶跃响应的飞行器单独静、动导数数值求解方法

基于阶跃响应方法,结合刚性动网格技术,对飞行器的单独静、动导数的精细化数值计算进行了研究.以纵向为例,通过给物面施加恒定附加攻角,求解得到阶跃响应运动过程的非定常气动力,求导得到静导数.同样给物面施加恒定的俯仰角速度,并同时强迫物面平动以抵消俯仰转动产生的附加攻角影响,可由非定常气动力求导得到动导数值.分别利用NACA0012翼型和三维SACCON飞翼无人机进行了计算验证,各攻角下的静、动导数值与文献、试验结果吻合得很好,最大误差不超过5%.结论表明:基于阶跃响应的单独静、动导数直接模拟方法计算耗时仅为传统强迫振动方法的21%,效率相对较高,且可推广到横航向的动导数计算,为飞行器的稳定性研究提供参考.      

柔性扑翼气动结构耦合特性数值研究

微型扑翼体积小、重量轻,其柔性变形对气动特性有显著的影响。通过求解雷诺平均N-S方程(ReynoldsAveraged Navier-Stokes,RANS)和结构动力学方程,对微型柔性扑翼飞行器的气动结构耦合特性进行了数值模拟研究。针对微型扑翼的大幅运动,发展了适用于扑翼的气动结构耦合数值计算方法,研究了微型扑翼的气动结构耦合特性。通过求解雷诺平均Navier-Stokes(RANS)方程得到微型扑翼的非定常气动特性;利用哈密顿原理(Hamilton Principle)推导了扑翼的结构动力学方程,采用结构有限元方法对该动力学方程进行离散并求解,得到扑翼的动态结构特性;采用松耦合方法进行迭代。计算结果与风洞实验结果相比吻合良好,验证了所发展方法的有效性。在此基础上研究了惯性力和关键运动参数对柔性扑翼气动及结构特性的影响规律,有助于比较详细、全面地了解微型扑翼的气动机理,为柔性扑翼的设计提供了参考依据。     

旋翼非定常气动特性CFD模拟的通用运动嵌套网格方法

针对直升机旋翼非定常气动特性CFD模拟中的网格生成难题,提出了一套高效、通用的运动嵌套网格生成方法.首先,基于Poisson方程求解和翻折法生成旋翼桨叶的正交贴体网格.其次,针对旋翼桨叶的扭转分布及变距、挥舞等复杂运动,建立了一套通用的洞单元识别的扰动衍射法;为保证洞包络面的封闭性,完善了挖洞过程中网格加密策略;在洞边界确立基础上,提出了一种高效、鲁棒的最小距离法贡献单元搜索的改进方法.在此基础上,建立了基于RANS(Reynolds-averaged Navier-Stokes)方程的旋翼非定常流场CFD模拟方法.最后,采用所建立的方法分别对悬停和前飞状态下的C-T(Caradonna-Tung)旋翼和7A(Helishape 7A)旋翼的气动特性、桨尖涡的位置进行了计算,计算结果与试验值误差小于5%,验证了该运动嵌套网格生成方法在旋翼非定常气动特性CFD模拟中的有效性.      

波浪形非均匀间隙封严结构影响涡轮性能的数值模拟

利用数值模拟的方法研究了波浪形非均匀间隙封严结构和均匀轴向间隙封严结构下轮缘封严气流对涡轮性能的影响.研究表明:燃气入侵与出流结构受到静盘、动盘及主流切向速度的影响,以低于动盘转速同向旋转,并改变了转子的进气条件,增强了压力面马蹄涡强度,因此对转子出口流场造成很大影响.封严气流与上游导叶尾迹的相互作用引起转子通道内熵增,造成涡轮效率的下降.与均匀轴向间隙封严结构相比,波浪形非均匀间隙封严结构使大的入侵与出流结构破碎为小的结构,对涡轮性能的负面影响减小,涡轮效率提高了0.9%.结果证明了波浪形非均匀间隙封严结构在具有较好的封严效果的同时提高了涡轮性能.      

穿越微下冲气流的飞翼布局无人机控制方法

微下冲气流是最危险的低空风切变形式,为在起降阶段安全穿越该气流,飞翼布局的无人机控制律应具有快速响应能力和良好的鲁棒性。针对大展弦比飞翼布局无人机舵面附加升力大和低速状态俯仰操纵效能低的特点,提出了舵面附加升力和机体气动力相结合的复合控制方案,改进了以输出误差为参考量的非线性指令分配策略,设计了基于迎角保护的指令分配策略。将风干扰和模型的不确定性视为未知扰动,采用自抗扰控制(ADRC)理论设计飞翼布局无人机非线性控制律,使之对风干扰和模型的不确定性进行估计补偿。仿真结果表明,复合控制与ADRC相结合的方法加速了航迹倾角的单位阶跃响应速度,使上升时间缩短了64%,同时能够实现对风干扰的有效观测和补偿,使高度损失低于2m;能够在风切变中有效保护迎角,使其维持在5.5°以内。因此,该方法能够为飞翼布局无人机安全平稳地穿越微下冲气流提供一种参考方案。     

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

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