混合解析/数值方法及其在湍流数值模拟上的应用

针对带小时间尺度的源项的方程描述的流动问题,提出了混合解析/数值方法。混合解析/数值方法的基本思想是:分裂原始方程组为对流-扩散部分的偏微分方程和源项的常微分方程。偏微分方程采用合适的数值方法求解,而常微分方程采用解析方式积分。模型方程的理论误差研究表明,混合方法提高了源项处理的精度,降低了混合方法的整体数值误差。分析同时表明,基于时间分裂的算法在求解含源项双曲系统的定常类型问题,会存在数值振荡。为此发展了非分裂方式的混合解析/数值方法,在湍流模型数值计算中提高了数值稳定性,而且加快了计算的收敛速度。     

Detached-Eddy Simulation方法模拟不同类型翼型的失速特性

 应用DES(Detached-EddySimulation)方法数值模拟了3种不同失速类型的翼型的升力特性。DES方法结合了RANS(Reynolds-averagedNavier-Stokes)和LES(LargeEddySimulationapproaches)的优点。基于Spalart-Allmaras湍流模型,在近壁面DES体现为RANS模型的特点而在远离物面处又具有LES的亚格子模型的特性。对此模型使用了LU-SGS隐式格式求解。通过和实验结果对比,显示这种方法可以有效地预测翼型的失速特性。     

转静态封严篦齿流场的数值分析和封严特性的试验

采用标准k-ε模型及RNGk-ε模型验证了前人的试验结果,最大误差分别为2.5%,15.9%。对4道、7道环型直通式篦齿,在转静状态下采用非结构化网格和SIMPLE算法,求解了不同转速的速度场和压力场;并试验研究了不同转速、齿数工况的封严效果,探讨了篦齿前后压比与泄漏系数的变化关系。表明转速对篦齿封严的影响不明显,变化幅度在3.6%;低压比情况下,随着压比增大,在相同泰勒数变化范围内,泄漏系数略有增大。      

民用飞机颤振试飞技术研究

为了能够顺利完成ARJ21-700飞机颤振试飞,研究民用运输类飞机颤振试飞的依据与要求,提出民用运输类飞机颤振试飞的测试方法、试飞方法、激励方法和数据处理方法;结合ARJ21-700飞机合格审定试飞进行试飞验证,给出该型飞机颤振特性试飞的符合性验证结果。结果表明:该套方法适合民用运输类飞机合格审定颤振试飞。     

模型的转捩流动计算分析简

 为了在黏性流动数值模拟中实现边界层转捩的自动预测,将γ-Re_θt转捩模型引入到三维非结构混合网格的雷诺平均Navier-Stokes方程求解程序（HUNS3D）。该转捩模型由两个依赖当地变量定义的关于间歇因子和当地化转捩起始动量厚度雷诺数的输运方程组成,其数值求解算法与流场求解程序中湍流模型的求解方法相同。为了考察和验证HUNS3D程序中γ-Re_θt转捩模型对航空工程中的常见附面层自由转捩问题的预测精度,对低速平板流动、Aerospatial-A翼型、NLR 7301超临界翼型和NASA Trap wing 高升力构型等典型外形的自由转捩流动进行了计算,并将计算结果与相关试验结果进行了对比分析。算例结果表明：γ-Re_θt转捩模型对于转捩位置具有很好的敏感性,能比较准确地预测自然转捩和分离转捩,可以有效提高HUNS3D程序对实际流动的模拟能力和预测精度。     

Assessment of Two Turbulence Models and Some Compressibility Corrections for Hypersonic Compression Comers by High-order Difference Schemes

The Spalart-Allmaras (S-A) turbulence model,the shear-stress transport (SST) turbulenee model and their compressibility corrections are revaluated for hypersonic compression corner flows by using high-...     

临近空间螺旋桨低雷诺数高效翼型数值分析

针对临近空间翼型边界层变厚更容易发生分离的特点,在修正了湍流模型后,通过分析比较13种低雷诺数翼型的升阻系数得到了比较适合用于临近空间螺旋桨叶素的高效翼型。结果显示,攻角在-6°～8°、8°～14°、14°～22°三种工况下,分别采用SST k-ω、RNG k-ε、Realizable k-ε湍流模型,可以得到与国外实验比较接近的合理结果;工况条件对翼型气动性能好坏的影响很大,翼型S-1223和FX63-137在所研究的工况内都具有较好的气动性能;因此可以选择这两种翼型作为临近空间螺旋桨用高效翼型。     

用转捩模型预测转捩及确定最佳粗糙带高度

采用SST k-ω二方程湍流模型,附加γ-Reθ转捩模型对不同迎角下的NACA0012翼型绕流进行了转捩预测,计算预测的转捩位置与实验结果符合较好,说明γ-Reθ转捩模型对转捩位置具有较好的预测能力。进一步用模型对带有粗糙带的翼型进行转捩预测,获得了不同雷诺数下刚好激发转捩又不引起过大局部压力扰动和附加阻力的最低粗糙带高度,对固定转捩风洞实验有指导意义。     

抑涡孔气膜冷却的大涡模拟

采用大涡模拟研究抑涡孔气膜冷却的流动和换热机理.通过与相同工况下圆孔气膜冷却流场的湍流结构进行比较分析,得出辅孔射流与主孔射流之间的干涉作用,并探索大尺度湍流结构影响气膜冷却效率的物理机制.结果表明:①辅孔射流抑制主孔射流形成的反转涡对的尺寸和强度,并为主孔射流的卷吸提供冷气补充;②由于上游辅孔出流冷气的保护作用以及对主孔出口尾缘低压区的冷气补充作用,主孔射流在气膜孔出口处温度没有急剧升高;③由于辅孔射流的干涉作用,主孔射流并未形成完整的发卡涡结构,而是无规律的近壁条带结构,减少了冷气和主流的掺混,并使冷气更好覆盖和冷却壁面.      

矩形横截面螺旋管湍流流动与传热特性的数值研究

 采用了剪切应力输运(SST)k-ω两方程湍流模型并考虑近壁低雷诺数的修正对矩形横截面螺旋管内冷却水流动和传热特性进行了数值研究。数值分析了在不同入口雷诺数、曲率半径以及扭距条件下,螺旋管内的温度、速度场以及流线的变化,讨论了螺旋管内、外壁面对流传热系数的差异及产生机理,同时与直通管道传热性能进行了比较。研究发现由于离心力的作用,螺旋管内存在显著的二次流动,管内、外侧壁面对流传热存在差异。旋转一周后,螺旋管即进入了流动稳定状态,入口雷诺数可以显著提升螺旋管整体的对流换热效率,扭矩和曲率对内外壁面传热效果的影响不大,而窄高型的横截面构型可以显著改善螺旋管的传热效果。研究结果对应用矩形横截面螺旋管的冷却设计提供参考。