平板边界层转捩的数值模拟

对于自由湍流引起的压力梯度为零的平板边界层转捩过程的数值模拟,本文提出了一个单方程湍流模型 (VL M模型),对湍流粘性系数和湍流长度尺度都进行了修正以正确引入低雷诺数效应。数值计算表明,该模型在自由来流湍流度为 0.03～8%的大范围内能正确模拟边界层转捩全过程,包括转捩区起点位置和转捩区长度、边界层内的平均及脉动参数分布等。本文还提出了上述转捩过程的时均化五阶段模型。      

关于两方程湍流模型的考虑

通过对两方程湍流模型的分析,提出了在主流区或远壁区用两方程模型,近壁区用代数湍流模型的思路,即所谓的“双层”模型概念。通过对平板湍流附面层的计算表明,这种分层计算措施获得的速度分布、壁面剪切力较单纯应用低雷诺数两方程湍流模型时的结果准确。对非均匀来流条件下90°弯管内高湍流度流场的计算表明,双层模型可获得更加准确的速度分布。从而表明本文提出的分层计算措施是有效的,可满足工程计算要求。      

[γ-Reθt]转捩模型在跨声速涡轮叶栅中的应用

为了解涡轮叶栅在跨声速条件下的流动特性和准确预测涡轮叶栅外换热情况,对γ-Reθt转捩模型,依据零压力梯度时自由流湍流度衰减实验结果,给出了一种直接估算来流粘性比的方法,以保证叶片前缘附近具有正确的自由流湍流度分布,提高换热预测准确度;同时减少试算次数。对MARK II与VKI两种叶栅通道跨声速工况下的流动换热情况使用CFX软件,选取层流模型、SST k-ω模型以及缺省粘性比和设定合理粘性比的γ-Reθt转捩模型进行了数值模拟验证,计算结果与实验数据的对比表明:转捩模型优于其他模型;而采用本文方法给定进口粘性比,能准确预测转捩位置,同时显著改善γ-Reθt转捩模型对不同来流湍流度下涡轮叶栅表面换热的预测精度;当入口湍流度较高,相比采用缺省粘性比情况,压力面上换热系数的相对误差降低30%以上,控制在7%左右。     

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

用低雷诺数k-ε湍流模型计算涡轮叶片的对流换热

考虑压力面强顺压梯度及吸力面逆压梯度对Ｓｃｈｍｉｄｔ Ｐａｔａｎｋａｒ低雷诺数湍流模型进行改善,使之能用于模拟涡轮叶片上的对流换热情况。计算了６种涡轮叶片的１８个工况。参数范围是：出口雷诺数Ｒｅ２＝０５６×１０６～２７３×１０６;来流湍流度Ｔｕ∞＝０８％～８３％;平均壁温与气流温度比Ｔｗ／Ｔ０＝０６７～０８２。结果表明,在叶片上的传热计算与实验符合得很好     

对 k-ω 湍流模型的改进

为了更适合于较复杂几何外形的数值计算，本文提出仅引入１个参数来改进ｋ－ω湍流模型并采用了一种双网格系统的有限体积格式来求解Ｎ－Ｓ方程与改进的的ｋ－ω模型方程。通过平板附面层的计算结果表明了用本文提出的对ｋ－ω湍流模型的改进对控制转捩是有效的。     

旋涡运动与湍流模型

本文力图将熟知的湍流大涡运动现象与湍流模型理论联系起来以减少模型理论中的经验关系。 在剪切流中以球的流体动力系数,分析确定了涡球的运动,并以此确定剪切流中的湍流剪应力。由此得到了一个二方程模型,方程中仅包含球的流体动力系数。 平板湍流边界层的实例计算表明,计算的速度剖面在整个边界层中都能很好地与实验符合。     

雷诺数对涡轮叶尖流场影响的数值研究

采用基于压力修正的三维计算流体力学程序，结合雷诺应力湍流模型和剪切应力传输湍流模型加壁面函数的方法，对某一轴流涡轮转子叶尖间隙流场进行了数值计算研究，详细计算了不同涡轮叶尖间隙高度和来流湍流度条件下雷诺数对涡轮转子间隙流场的影响，最后计算了转子效率。结果表明：当泄漏流流经叶尖时因为叶尖剪切力做功有块总压增大区；雷诺数带来的影响比湍流度和叶尖间隙高度带来的影响要大，湍流度的变化对流场影响不大；雷诺数对泄漏涡尺寸的影响不大，但低雷诺数会引起流动分离，带来损失，当雷诺数在文中的范围内变化时，效率下降近八个百分点。     

吹吸气流隔断局部空间的基础研究

本文对用吹吸气流隔断局部空间时的气流运动进行基础研究。本研究的气流流动假定为二维、湍流、等温和稳定流动。动量方程式中的从属变量选用流函数和涡度(ψ—ω),湍流动能方程选用一方程式模型(K—L模型),以此作为吹吸气流流动的数学模型。数值解析结果与实验结果比较表明,用上述方法推导出的数学模型和计算程序能较好的反映实际流动情况。本研究对吹吸气流隔断局部空间的研究、设计、现场调试将具有一定的理论价值和工程应用价值。     

压缩性修正对γ-Re_θ转捩模型的影响研究北大核心CSCD

将基于流场当地变量的γ-Reθ转捩模型加入可压缩RANS方程计算程序,对平均流场控制方程、湍流和转捩控制方程进行基于LU-SGS的隐式紧耦合求解。结合三种压缩性修正方法对高超声速平板、双楔和圆锥流动进行了数值模拟,给出了壁面斯坦顿数、热流值和实验值的比较,以及湍动能和间歇因子分布云图。数值计算结果表明,相同来流湍流度下不同压缩性修正方法得到的转捩位置差别较大。相比于未修正的模型,基于当地马赫数压缩性修正后的转捩位置最靠后,其对带有分离的双楔流动预测较为准确;基于来流马赫数的压缩性修正仅使得转捩位置稍有延迟;对湍动能源项的压缩性修正可提高模型在转捩后湍流段的热流预测精度。