两方程湍流模型的可压缩性修正及其应用

本文在不可压两方程湍流模型的基础上进行可压缩性修正,给出了ｋ－ε、ｋ－ε两方程模型的可压缩性修正方程,并对亚声平板湍流流动和跨声喷管流动进行了数值模拟,得到了较好的结果。从而证实了两方程湍流模型经过可压缩性修正,能较好地模拟可压缩流动。     

(高)超声速流动试验技术及研究进展

近年来,与高速飞行器相关的(高)超声速流动受到了极大的关注。这类流动所具有的非定常性、强梯度和可压缩性对试验方法和风洞设计技术提出了挑战。超声速纳米示踪平面激光散射(NPLS)技术是由作者所在团队研发的非接触光学测试技术。它能够以较高的空间分辨率来揭示超声速三维流场的一个瞬态剖面的时间解析的流动结构。介绍了NPLS技术以及基于NPLS开发的密度场测量、雷诺应力测量和气动光学波前测量等方法,并回顾了这些技术在超声速边界层、超声速混合层、超声速压缩拐角、激波/边界层相互作用和光学头罩绕流等流动中的应用,清晰地再现了边界层、混合层、激波等典型流场结构及其时空演化特性。另外,为了模拟和研究高空大气条件下边界层自然转捩和超声速混合层的转捩特性,介绍了高超声速静风洞、超-超混合层风洞的设计技术以及层流化喷管的设计方法。     

一类新的速度入流条件在管道流模拟的应用

基于经典边界层理论,发展了一类简化的、按边界层厚度描述的速度入流边界条件。基于入流条件,对低速后向台阶流动和超声速压缩拐角流动进行了数值模拟。通过细致比较空间流场、边界层内速度型、表面摩擦阻力和压力分布等特性,对该类型边界条件进行了计算确认。结果表明,所发展的简化边界层厚度速度入流条件提法正确、合理,具有描述方法简单、鲁棒性好、适合工程计算的优点,可显著简化不同流速的管道流动的数值模拟。     

亚/超声速楔状流层流边界层速度与温度相似解及拟合解

利用相似变换获得了楔状流层流边界层无量纲流函数的3阶非线性常微分方程,用Runge-Kutta法求解微分方程获得了不同楔形角楔状流层流边界层无量纲速度随相似变量的变化曲线;推导了亚声速和超声速楔状流层流边界层无量纲温度关于相似变量的2阶线性齐次和非齐次微分方程,获得了温度分布的通解,恒壁温条件下亚声速楔状流和绝热壁面条件下超声速楔状流层流边界层无量纲温度解析解及指数函数形式的拟合解.以楔形角为0为例利用相似变换研究了超声速条件下气体压缩性及黏度随温度变化等因素对层流边界层速度与温度的影响,得出不可压缩常物性与可压缩变物性条件下无量纲速度相对误差绝对值小于9.8%的结论.研究表明:Pr越大贴近壁面处无量纲温度变化越剧烈;超声速条件下壁温低于绝热壁温时黏性耗散作用可以使层流边界层气体温度从壁面到主流间出现先升高后降低的变化.      

基于边界层相似性解的放大因子输运模型

为了使二维线性稳定性理论能够适应现代CFD求解技术,通过求解Falkner-Skan边界层相似性方程获得各个形状因子下的相似速度型。基于线性稳定性理论,对各个速度型进行稳定性分析从而获得对应不同速度型的扰动放大因子包络线。最后使用标量输运方程的形式实现包络近似方法中放大因子的当地化求解,并结合原始γ-Re_θt转捩模型中的间歇因子输运方程,实现了自然转捩和分离泡转捩的建模。使用该输运模型对S&K平板、S809翼型、NLR7301翼型和DLR-F5机翼进行转捩预测,结果均与试验结果吻合较好,验证了该模型构建的合理性和可行性。     

翼身融合运输机分布式电推进系统设计及油耗评估

针对翼身融合运输机开展了分布式电推进系统的总体设计与油耗评估。通过数值计算完成了70t载质量翼身融合飞机的气动设计与优化。在巡航马赫数为0.80和10km高度的设计点,最大升阻比达到了24。通过求解积分边界层方程组,完成了电推进系统的总体设计。电推进系统包含10个推进风扇,风扇直径为1.45m,压比为1.35,巡航功率为2.94MW。建立了考虑燃烧过程的发动机一维性能模型,对发动机油耗进行了评估,获得了不同发动机循环参数下燃油消耗。建模结果表明,基于翼身融合布局和分布式电推进技术,可使运输机的油耗较C-17节省近50%。      

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

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

基于温敏漆的边界层转捩测量技术研究

以自然层流翼型RAE5243模型为研究对象,在0.6m跨超声速风洞进行温敏漆（Temperature Sensitive Paint,TSP）转捩测量技术研究,在Ma0.73和Ma0.75条件下开展了模型基本外形和鼓包外形的转捩测量试验。针对缺乏定量分析手段的问题,提出基于温度梯度的转捩位置自动判定算法,包括图像预处理、转捩点定位与筛选和转捩位置计算3个步骤。模型温度分布及转捩测量结果表明:重复性试验结果偏差较小,验证了转捩测量结果的可靠性;相同马赫数条件下,鼓包外形转捩位置相对基本外形向后缘移动;相同外形条件下,Ma0.75的转捩位置相对Ma0.73向后缘移动。TSP试验结果与CFD计算结果吻合较好,变化趋势一致,检验了数值模拟方法的有效性。     

采用长细管法进行脉动压力转捩探测的实验研究

为了简便地使用测量模型表面脉动压力特征的方法探测边界层转捩位置,需要研究脉动压力传感器接在传统测压模型外的适用性,即通过长细管将模型表面的脉动压力信号传递到脉动压力传感器上的方式是否可得到转捩的特征信号。首先采用信号发生器驱动扬声器,在无风条件下,测量了长细管对不同频率声压信号的传递损失情况。证明了所采用的长细管系统具有合适的工作频带。然后在西北工业大学NF-3低速风洞二元实验段、实验风速为30m/s的条件下,对弦长为800mm、展长为1.6m的翼型模型沿弦向进行了脉动压力信号测量,并通过改进的数据处理技术判断了模型表面的转捩位置。研究结果表明,采用长细管系统进行脉动压力方法转捩探测具有一定应用价值,值得进一步深入研究。     

一种修正的Osher通量在高阶WCNS中的特性

将一种基于多维修正的Osher通量运用于高阶加权紧致非线性格式(WCNS)中,该修正通量主要在垂直于激波面的界面上增加耗散,能够改善Osher通量的激波捕捉稳定性,同时对边界层和接触间断的分辨率影响非常小。对修正的Osher通量在高阶WCNS中的特性进行研究,通过数值模拟测试了基于Osher通量的WCNS的激波稳定性、热流预测精度、边界层模拟能力、激波边界层干扰模拟能力,并与Steger-Warming通量和Roe通量进行了对比。结果表明修正后的Osher通量比Harten修正的Roe通量具有更好的激波捕捉鲁棒性,而边界层、驻点热流值和激波边界层干扰的模拟则明显优于Steger-Warming通量。上述结果说明了基于修正的Osher通量的高阶WCNS具有较好的激波捕捉特性、热流预测精度和边界层计算能力。