应用改进的低雷诺数湍流模型预测转捩流动

基于Launder-Sharma(LS)低雷诺数k-ε两方程湍流模型发展了一种改进的具有转捩敏感性的低雷诺数湍流模型.针对LS模型的可实现性(realizability)问题和前缘滞止点湍动能预测过大的不足,模型进行了改进.改进模型只使用当地物理量,不需要求解壁面距离、y+和边界层积分参数.改进模型能够适用于广泛的流动,且容易应用到通用的计算流体动力学(CFD)程序中.对具有详细数据的零压力梯度平板转捩边界层T3A实验的模拟结果显示,改进模型能够预测转捩流动,并能对自由湍流变化给出合理的响应.      

边界层转捩的数值模拟

在数值模拟中将Abu-Ghannam&Shaw(AGS)转捩模型引入到湍流模型中,以考虑流动中的层流和转捩现象,并在此基础上,分别对低雷诺数情况下的平板和涡轮叶栅流动进行了模拟,并与实验结果进行了比较,表明加入转捩模型可以较好地描述边界层的转捩过程。同时,结果还表明AGS转捩模型对压力梯度变化较大情况下的模拟稍有偏差。      

简化三方程转捩模型在高超声速流动中的应用

在雷诺平均方法的基础上,通过耦合k_(-ω)SST湍流模型和间歇因子转捩模型,引入湍流模型和转捩模型的可压缩修正方法,对高超声速平板、双楔、尖锥三类模型边界层转捩流动开展了数值模拟研究。与实验结果的对比分析表明基于压力梯度表征参数T_w=R_(TΩ/ω)的简化三方程转捩模型,能够准确捕捉高超声速平板边界层流动的转捩起始位置、转捩区域长度以及湍流区壁面热流。而对于双楔、尖锥模型,改进前的简化三方程转捩模型由于受到流动可压缩效应的影响,边界层转捩后湍流区的壁面热流模拟预测结果明显高于实验值。在添加模型可压缩修正方法后,转捩区域长度和湍流区壁面热流模拟结果得到有效改善,与实验值吻合较好。可见,简化三方程转捩模型在添加可压缩修正方法后具备准确模拟预测高超声速边界层流动转捩的潜力。     

转捩模式与转捩准则预测高超声速边界层流动

对原始的k-ω-γ转捩模式和"层流+转捩准则"模型进行了改进,在2种方法中分别增加了横流模态时间尺度和横流转捩准则用于预测横流失稳诱导转捩。通过对网格预处理可并行计算获得边界层外缘信息以及边界层内横流速度。采用不同雷诺数条件下的0°攻角尖锥以及HIFiRE-5外形对2种方法预测高超声速边界层转捩的性能进行了比较分析。研究结果表明,2种方法均能正确反映高超声速边界层转捩起始位置和转捩区长度随雷诺数的变化趋势,但不能捕捉转捩区热流峰值;"层流+转捩准则"模型计算得到的传热系数在全湍流区较k-ω-γ转捩模式偏高。对于同时存在流向不稳定和横流不稳定的HIFiRE-5外形,改进的k-ω-γ转捩模式和改进的"层流+转捩准则"模型相比于原始的模型均能更加准确地预测中心线两侧横流失稳诱导形成的转捩;对于中心线附近因速度剖面拐点引起的边界层转捩,"层流+转捩准则"模型由于与边界层厚度相关,预测得到的转捩位置较试验结果靠前,k-ω-γ转捩模式与试验结果吻合很好。     

γ-Re_θ模式应用于高速边界层转捩的研究

应用γ-Reθ转捩模式对超声速、高超声速边界层转捩进行数值模拟。γ-Reθ模式通过求解关于当地雷诺数和间歇因子两个输运方程给出转捩起始位置和转捩区长度等信息。本文对马赫数3.5至7范围内的四种算例进行计算,研究了来流雷诺数、攻角变化和头部钝化半径等关键参数的变化对γ-Reθ模式预测流动转捩性能的影响。给出了壁面摩擦阻力系数、热流值与实验值的对比以及壁面附近间歇因子等值线分布等计算结果。γ-Reθ模式能够正确预测攻角和钝化半径变化时转捩位置和流动参数的变化趋势;在较大雷诺数时计算结果与实验值吻合很好。γ-Reθ模式对于超声速或高超声速边界层转捩的模拟仍需修正和改进。     

基于湍流模型的转捩预测方法研究

阻力的精确计算,一直是低雷诺数飞行器设计的关键.而与附面层分离和转捩相关的摩擦阻力的精确计算则是此类飞行器阻力精确计算的关键.因此,开展了基于湍流模型的转捩预测方法的研究,应用几种湍流模型对平板边界层进行计算分析,研究各湍流模型预测流动转捩的能力.由于各模型对扰动的过于敏感,引起转捩提前发生,使预测与实验差距较大,由于k-ωSST两方程湍流模型预测能力优于其它模型,在此模型基础上引入了间隙函数进行粘性系数修正,对平板和NACA 0012翼型的转捩重新进行了预测,计算结果表明改进后的模型对转捩位置具有较好预测能力.     

翼型几何模型对边界层转捩预测精度的影响

对飞行器表面的边界层转捩现象进行准确地预测,是设计高性能层流翼型的关键。基于边界层数值模拟和现代CFD技术的特点,选取与SSTk-ω湍流模型相耦合的γ-Reθt转捩模型,对不同类型的翼型边界层转捩进行预测。结果表明:在转捩计算过程中出现的摩阻曲线震荡,是由翼型表面计算网格不光顺所引起的;提高翼型几何模型光顺程度能够在不同程度上提高翼型气动特性的计算精度,并能得到更加准确的转捩预测结果。     

采用γ-Re_(θt)模型的转捩流动计算分析

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

后掠机翼人工转捩最佳粗糙带高度数值预测

介绍了基于当地变量的γ-Re_θ转捩模型,并将该模型应用到后掠机翼的转捩预测和人工转捩最佳粗糙带高度以及人工转捩技术能够模拟的大气飞行雷诺数的确定中。为检验γ-Re_θ转捩模型对后掠机翼转捩的预测能力,对ONERA M6机翼和DLR-F4标模机翼进行了边界层转捩预测,采用结构化网格和有限体积法求解雷诺平均Navier-Stokes(RANS)方程,得到了机翼表面的摩擦阻力系数分布,从而可以得到相应的转捩位置,预测得到的转捩位置与试验结果比较吻合,说明该模型对后掠机翼转捩预测是可信的。最后在DLR-F4标模机翼上表面固定了粗糙带,通过相同的方法得到了转捩位置,从而确定了马赫数为0.785、雷诺数为3.0×10⁶时最佳粗糙带高度为0.11 mm;通过不断增大雷诺数使自由转捩位置不断向前缘移动,验证了人工转捩对大气飞行雷诺数的模拟能力。结果表明,在最佳粗糙带高度为0.11 mm下,可以实现对大气飞行高雷诺数的模拟。     

微通道流动转捩数值模拟

应用CFD计算软件ANSYS CFX对水力直径为0.4mm的矩形截面微尺度通道中的气体流动进行了数值模拟.采用了3种流动模型,即γ-Re_θt转捩模型、层流模型和剪切应力输运（SST）模型,并将模拟结果与实验结果进行了对比.模拟结果表明:层流模型及SST模型不具有预测转捩的能力;γ-Re_θt转捩模型捕捉到了临界雷诺数且摩擦因数相对误差在20%以下.在未达到临界雷诺数之前,通道后部已经出现局部湍流区,且局部湍流区随着雷诺数的增大而增大.壁面摩擦因数的最低点对应的流向位置与通道中前后湍流区的交结位置几乎一致.      

过渡区球锥气动力特性的桥函数估算与DSMC仿真研究

采用DSMC(DirectSimulationMonte-Carlo)方法模拟球锥体过渡区流动。论文给出了不同Knudsen数球锥绕流流场的温度与压力分布图,以及在不同Knudsen数、不同马赫数条件下球锥阻力系数和底部阻力系数的分布,并将阻力系数模拟结果与桥函数一体化估算结果进行了比较分析。     

Transitional Flow and Heat Transfer Characteristics in a Rectangular Duct with Stagger-arrayed Short Pin Fins

This article describes an experimental study on friction and heat transfer performances of a transitional airflow in a rectangular channel with stagger-arrayed short pin fins. Friction factors, average Nusselt numbers and overall thermal performances of the transitional flow are obtained. The experimental study has showed that the pin fins enhance the heat transfer performance significantly, however increasing the flow frictional resistance considerably. After comparing the experimental results with the published data in references, it can be concluded that, in the transitional flow region, the pin fin channels of the proposed geometrical configuration could lead to a significant improvement of an overall thermal performance; for instance, the convective heat transfer performance is increased by at least 68%. By contrast, in the fully turbulent flow region, the ability of the proposed pin fin channels to increase heat transfer performances decreases as the Reynolds number increases. When Re > 6 000, the overall thermal performance becomes lower than the others.     

过渡状态的倾转旋翼气动特性计算分析(英文)

针对特殊的旋翼倾转运动,建立了一个过渡状态的旋翼非定常气动力数值计算方法.为正确模拟旋翼倾转运动使桨叶受到的附加惯性力及哥氏力作用,重新推导了旋翼倾转时的桨叶气动力模型和挥舞运动方程;为了适合于旋翼倾转时的入流和气动力计算,入流模型中考虑了倾转运动引起的旋翼尾迹弯曲影响.应用建立的方法,首先进行了旋翼配平计算,以验证计算模型,并给出了倾转旋翼的操纵量.然后,着重计算了旋翼在倾转过渡时的拉力.俯仰和滚转力矩随倾转角的变化,分析了倾转飞行时的前飞速度.倾转时间等对旋翼气动力的影响,得出了一些新的结论.     

确定低雷诺数翼型转捩分离泡位置的实验研究

在翼型模型的表面粘贴表面热膜,由其给出脉动电压的均方根值和波形图,可测出层流边界层分离点和湍流边界层再附点,转捩分离泡的位置也就确定了。     

一种基于结构化贴体网格的DSMC方法研究

基于DSMC方法,发展了一种结构化的贴体网格的分子运动轨迹跟踪和定位的方法———网格面法向量判断法。在此基础上建立了一套DSMC数值模拟程序,对球双锥和Apollo飞船在过渡区域流动进行了数值模拟,并分析了其过渡区域超音速、高超音速的气动力、气动热等问题。通过将本文的计算结果与参考文献的计算结果和实验结果进行了比较,验证了所发展的方法的可靠性。     

过渡领域三维绕流直接统计模拟位置元法的一种新方案

本文发展了ＤＳＭＣ方法计算三维绕流位置元方法的新方案，其总体方案和实施细节都是在分析和探索的基础上由我们自己提出的，特别还解决了除气动力、矩总体量以外的气动加热和压力分布等局部分布量的计算所提出的新问题。对于表面位置元的标识、分子的进入及与表面碰撞的判断、网格初始划分与调整、分布量计算对模拟过程及信息存贮所提的特殊要求等进行了分析和讨论，并通过与准确结果的比较证明了采用的方案和技术的有效性。     

过渡流区N-S/DSMC耦合计算研究

在已有CFD和DSMC方法程序研究基础上,采用MPC耦合处理原理,发展适于流场分区信息交换的亚松弛技术,建立了N-S/DSMC耦合算法.通过对近连续过渡流区不同Kn数条件下钝锥绕流计算与低密度风洞试验结果比较分析,一方面证明了本文所建立的N-S/DSMC耦合计算模型和方法在低Kn数过渡流区的有效性,另一方面剖析了N-S/DSMC耦合计算误差,可望发展一个工程适用的近连续过渡流区高超声速飞行器气动特性预测分析工具.     

飞机扭转敏捷性与飞行品质关系研究

借鉴模拟或试飞获得飞机绕飞行速度滚转且截获 90°倾斜角所需时间 TRC90 的技术,设计数字仿真用的飞机扭转敏捷性机动动作及其相应的舵面操纵规律。具体计算某机不同 H,Ma和 ny 下的 TRC90 和 TA尺度,以及定常水平盘旋飞行时的模态特性。讨论和分析飞机扭转敏捷性指标与飞行品质之间的关系     

复合材料桨叶过渡剖面边界区域应力传递边界效应初探

在复合材料桨叶载荷标定前的解耦或应力应变测量等工作中,会碰到一个“奇怪”的现象：在桨叶过渡剖面边界区域测出的应变（或应力）与“平剖面假设”规定不符。文章从某型号尾桨叶载荷标定前解耦调片工作中碰到的现象引出论点,再用有限元对桨叶传力进行建模分析,给出了定性的结论。文章的目的就是从型号中来,再应用到型号工作中去,为今后的复合材料桨叶载荷测量的贴片位置选择起一指导作用。