可压缩粘性流动笛卡尔网格虚拟单元方法研究

以二维高雷诺数可压缩粘性流动问题为背景,提出了一种全新的笛卡尔网格虚拟单元方法。基于壁面函数基本假设,构造了壁面函数-虚拟单元方法(WF-GCM),用于定义湍流壁面边界条件。引入参考点的概念计算虚拟单元上的基本变量与湍流变量值,定义了"非贴体"笛卡尔网格下的湍流壁面边界条件,并通过壁面函数模型修正近壁面单元与界面单元。基于自适应笛卡尔网格体系,采用发展的具有二阶精度的格心格式有限体积求解器,数值模拟了跨音速RAE2822翼型绕流问题与超音速圆柱绕流问题,计算结果与实验值吻合良好,显示了WF-GCM对高雷诺数可压缩粘性问题是有效的。     

对称交叉激波和湍流边界层相互作用的数值研究

针对两种双尖鳍外形的对称交叉激波与湍流边界层相互作用,采用N-S方程和两种湍流模型进行了计算。研究了网格收敛性、鳍的角度和湍流模型对壁面压强、Stanton数和壁面摩擦力线的影响。弱相互作用的计算结果较好,强相互作用的壁面压强和摩擦力线的计算结果与试验吻合较好,而Stanton数的计算结果较差,峰值高达试验的2.5倍左右。随着鳍的角度的增加,壁面压强和Stanton数的分布从单调分布发展为M型分布,两者的峰值不在相同的位置。湍流模型对壁面压强和壁面摩擦力线影响很小,对Stanton数计算的影响很大,SST模型比BSL模型表现好一些。     

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

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

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

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

非对称交叉激波和湍流边界层相互作用的数值研究

针对7°×11°双尖鳍外形的非对称交叉激波与湍流边界层相互作用,采用Navier-Stokes方程和5种湍流模型进行了计算.主要考察对壁面压强、热传导、绝热壁面温度和壁面摩擦力线分布的计算精度.计算结果表明:两道斜激波相交后的区域的壁面压强和热传导都比较高;计算的压强和壁面摩擦力线与试验吻合很好,绝热壁面温度次之,热传导最差,峰值高达试验的3倍左右.湍流模型对壁面压强和壁面摩擦力线影响很小,对绝热壁面温度和热传导影响很大.在5种湍流模型中,TNT(turbulent/non-turbulent)和SST(shear-stress transport)模型表现较好.      

二维沟槽粗糙槽道内的湍流特性研究

采用激光多普勒测速技术对光滑和粗糙槽道湍流特性进行了实验研究。粗糙元为二维横向V型沟槽,沟槽深度为0.8mm,沟槽间距为6.4mm,对应的槽道半高度与沟槽深度比为12.5。基于中线时均速度和槽道半高度的流动雷诺数范围为2740～17400。实验测量了包括时均速度、湍流强度、雷诺切应力和速度脉动偏斜因子和平坦因子在内的湍流统计量,结果表明沟槽型粗糙度对湍流的影响不仅局限于边界层内区,而是延伸到整个边界层范围。粗糙壁面上的粗糙度函数随雷诺数的增大而增大,时均亏损速度也较光滑壁面高。沟槽抑制了内区的流向湍流强度,同时增大了外区的湍流强度。粗糙壁面上的雷诺切应力高于光滑壁面,且与湍流强度一样表现出对雷诺数的依赖性。尽管沟槽型粗糙度对流向平坦因子影响不大,但对流向偏斜因子有显著影响。     

气膜冷却恢复区的湍流模型

在 L am- Bremhorst低雷诺数 k- ε模型的基础上 ,发展了气膜冷却恢复区的湍流模型。原模型中的常数 cμ 和 c1代之以沿 x、y方向变化的参数 ,低雷诺数函数 fμ、f1和 f2 沿 x方向有所变化。修正后的模型可以在 0 .4～ 1.0的吹风比范围内 ,较好地描写恢复区内时均速度、湍动能和湍应力分布 ,并能正确模拟向不吹风情况过渡的恢复过程。     

SSG模型在压缩拐角激波/湍流边界层干扰模拟中的应用

采用非线性方案模化雷诺应力再分配项的二阶矩雷诺应力输运模型——SSG模型和两方程线性涡粘模型SST模型,对24°高雷诺数二维压缩拐角激波/湍流边界层相互干扰流动进行数值模拟。通过在壁面压力分布、壁面摩擦阻力系数分布、截面速度分布及分离区大小等方面与实验数据比较,综合评估了上述两种模型模拟该问题的能力。结果表明:两种湍流模型在上述各方面都取得了与实验较为一致的结果,但相比较而言,除SST模型在流动再附后对壁面摩擦阻力系数分布的模拟更接近实验值外,SSG模型结果在上述其他方面均更接近实验值。本文工作为压缩拐角激波/湍流边界层干扰这类流动数值模拟的湍流模型选择提供了参考。     

GAO-YONG湍流模型对湍流传热的数值模拟

将GAO-YONG湍流模型应用于湍流传热的研究,分别计算了平板剪切湍流和二维平面冲击射流的湍流传热问题.边界层剪切湍流流动与换热的计算表明:与传统的湍流模型不同,GAO-YONG湍流模型不需要对近壁区域做任何特殊处理(比如壁面函数、低Reynolds数修正等)即可模拟出从壁面到主流区的全部流动与传热情况;另外,对于冲击射流Nusselt数的模拟也得到了与实验符合较好的计算结果,准确地捕捉到了2种冲击高度下流场换热的不同特征,表明了GAO-YONG湍流模型能够较高精度地计算湍流换热.      

壁面函数在超声速湍流模拟中的应用

通过子迭代的方式计算摩擦速度,更新虚拟点湍流黏性系数实现对壁面剪切应力的修正,耦合SSTk-ω两方程模型,在国家数值风洞软件平台上实现了壁面函数程序模块。通过压缩拐角和高速飞行器典型算例进行考核,初步数值实验结果表明：压缩拐角算例壁面函数在无量纲壁面距离y+≤200范围内,均准确预测湍流边界层速度分布,可显著提高粗网格上壁面湍流边界层和壁面摩擦系数的预测精度,且壁面函数使用粗网格最多可节约75%的计算时间;对于复杂外形,附面层网格间距变宽,可使得整体网格减少约38%的网格总量,在相同的计算设置情况下可节约60%的计算时间。唇口附近出现明显的激波边界层干扰现象,使用壁面函数后稀网格和密网格得到流场中分离、激波反射相同,且分离区最高压力系数偏差从15%降低到2%。从粗网格和密网格全机轴向力系数比较来看,使用壁面函数后摩阻预测偏差可从40%下降到4%。整体来看,壁面函数的引入提供了一种高效的飞行器湍流流动气动力预测方法。     

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

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

可压缩湍流边界层壁面函数方法综述

以建立工程实用的高超声速湍流边界层模拟方法为目标,从湍流壁面函数是湍流边界层方程近似解的角度,梳理了相关文献的研究工作,得到如下认识:1)壁面函数与所求定解问题数值解的相容程度决定了壁面第一层网格允许放粗的程度,在流动分离点和再附点附近区域,目前壁面函数尚需进一步完善,而"子网格"壁面函数从理论上解决了相容性问题,尽管...     

扩散界面浸入边界法结合壁面模型在湍流模拟中的应用

提出了一种模拟高雷诺数湍流的扩散界面浸入边界法（IBM）。该方法采用壁面模型来减少壁面附近的网格量。为了实施壁面模型和边界条件，在物面外部设置了2个辅助层（一系列拉格朗日点）：外侧的参考层用于实施壁面模型来得到壁面剪应力，内侧的强制层用于实施离壁的边界条件。在实施壁面模型时，将动量方程沿物面法向积分，从而把物面切向的速度修正量与由壁面模型得到的壁面剪应力联系起来，而速度的法向分量则按二次曲线分布来近似重构。在实施边界条件时，为了严格满足无穿透条件，应用了基于隐式速度修正的IBM。最后，利用绕平板湍流和绕NACA0012翼型湍流来对方法的可行性进行了验证。     

两种湍流模型在跨声速绕流计算的应用研究

湍流模型在对复杂流场的数值计算中起着非常重要的作用.本文采用SA一方程湍流模型和SST k-ω二方程湍流模型,通过解耦求解雷诺平均N-S方程和湍流模型方程实现对亚跨超声速湍流流场的数值模拟.对NACA0012翼型和ONERA-M6机翼跨声速绕流流场进行了计算,对压力分布和激波位置与实验结果进行了细致的比较,并分析了不同离散格式、不同网格疏密及壁面函数对计算结果的影响,在计算过程中这两个模型体现出了较好的简捷性和健壮性.     

一种剪切修正的涡黏性亚格子模型及其在槽道湍流中的应用

提出了一种利用壁面附近强剪切对涡黏系数进行修正的亚格子模型,称为剪切修正的涡黏模型.该模型能使涡黏系数在壁面附近自动满足壁面修正,在脉动为零的地方,涡黏系数自动为零,保证了壁面附近涡黏系数趋于零的特性.用该模型对槽道湍流壁面无吹吸、有吹吸,较低和较高雷诺数的情况进行了数值模拟,与直接数值模拟的结果进行了对比,结果吻合较好.      

基于大涡模拟分析气膜冷却的湍流场

用大涡模拟考察了单孔平板气膜冷却的湍流场,气膜孔沿流向倾斜30°,气膜出流的雷诺数为2600,吹风比为0.5.计算结果表明:①气膜孔内的流动分离增加了湍流动能,导致气膜冷却效率的降低,在设计中要尽量避免或减少流动分离;②与射流侧向扩展有关的涡黏性系数在气膜孔两侧存在峰值,而采用各向同性的湍流模型预测气膜冷却时,涡黏性系数的峰值出现在射流与主流的剪切区,因此会高估射流的垂直穿透而低估射流的侧向扩展;③可以用大涡模拟辅助建立各向异性的湍流模型,以便提高湍流模型的预测精度.      

桥梁断面静风荷载的格子Boltzmann方法数值计算

通过引入反映湍流涡粘性的湍流松弛,得到了模拟高雷诺数湍流的BGK方程.在速度相空间、物理空间和时间上对BGK方程进行离散得到了三维十九速离散速度模型;结合分区计算技术,设计了格子Boltzmann并行算法;根据亚格子Smagorinsky模型,提出了直接从粒子分布函数计算湍流松弛时间的方法.用开发的并行计算程序对分体双箱截面和闭口箱梁截面的静风荷载进行了数值识别,得到的静力三分力系数和流场压力分布与风洞试验结果及CFD宏观方法计算结果吻合,并从表面压力分布入手分析了两种桥梁截面的绕流特点.     

湍流模型对叶片表面换热计算的影响

为了比较湍流模型对涡轮叶片外换热计算结果的影响,采用五种湍流模型及两种壁面函数计算了NASA-MarkII导向叶片表面换热并与实验数据进行了对比,结果表明:在相同边界条件下不同湍流模型的计算结果有很明显的差别,即使是同一种湍流模型,如果采用不同的壁面处理函数其计算结果也是有很大差别的;某些湍流模型的计算值只是在某个区域较为理想,还不能找到在整个叶片表面计算结果与实验数据较为接近的湍流模型.在尚无普适性较好的湍流模型的情况下,研究在不同计算域采用不同湍流模型和壁面处理函数的计算技术,是一种较好的可行的方案.      

槽流的亚临界转捩与局地湍流

在中低雷诺数时,处于线性稳定的槽道流若受到强扰动可发生亚临界转捩变为湍流。近年人们对该类转捩机理的研究取得重要进展。对于平板泊肃叶流,亚临界转捩之初是稀疏湍流态,其特征结构是远间隔的包含小尺度涡和高低速条带的大尺度的湍流带,可倾斜伸长。该阶段的湍流占比有上限但并非雷诺数的单值函数。随着雷诺数的增加转捩进入平衡局地湍流态,即存在统计定常态,其湍流占比是雷诺数单值函数,可由定向逾渗模型描述。进一步增加雷诺数,湍流带的分裂愈发频繁,最终流场会布满湍流带,在更高的雷诺数时变为均匀湍流。论文概述了为比拟亚临界转捩过程所提出的动力学模型,以及为定量表征管流、平板库艾特流和平板泊肃叶流转捩过程的相似性所提出的局地稳定性参数,并在最后对槽道流亚临界转捩研究的发展做了简要展望。     

航空发动机燃烧室流动数值计算中湍流模型的比较

分别采用标准k-ε湍流模型、RNGk-ε湍流模型、Realizablek-ε湍流模型以及雷诺应力模型,对某型航空发动机燃烧室流动进行了数值计算.近壁处理采用标准壁面函数法,计算得到速度矢量分布以及质量流量、湍流粘度比和湍流强度等参数.四种湍流模型计算的总体流动差别较小,但流场的细节有较明显的不同.标准k-ε模型、Realizable k-ε模型和雷诺应力模型的湍流粘性计算结果较为接近,而RNGk-ε模型计算的湍流粘性较小.