近壁区理论耗散率κ-ε模型地研究.pdfVIP

中国工程热物理学会 热机气动热力学学术会议 近壁区理论耗散率k一占模型的研究 陆利蓬 陈矛章 (北京航空航天大学404北京100083) ■ 摘要提出在湍流边界层近壁区采用共振三波的理论模型描述湍流相干结构，根据 r 理论模型计算了s的分布。并且在传-统k-E模式基础上依照理论占值计算了平均速度分 卜 布。在粘性作用层理论值与直接数值模拟符台很好。表明该理论方法不仅有益于对湍 l 流机制的了解，而且可能为湍流的近壁模型化开辟一条新的途径。 关键词近壁区：湍流模型：相干结构：湍能耗散率 湍能耗散率输运方程(s方程)是高阶湍流封}iJ模式的重要方程之一。但是，在 各湍流模型中．对s方程的模拟最不能令人满意，因为对方程中的各项知之甚少．特 别是起主要作用的小涡拉伸产生项与粘性破坏项都是由极高波数的小尺度涡决定的． 在现在和可预见的未来都是银难准确测量。因此，只是根据量纲分析．按直觉和类比 进行模拟，没有逻辑推理，缺少物理讥制。如果说现有的湍流模型理论带有一定的随 意性。其中对￡的处理则随意性虽大。湍流边界层的近壁区是许多湍流现象的营发区， 它对整个边界层湍流的发生发展羽l维持都有重要的作用。壁面以及粘性的作用使得其 流动特性区别予充分发展湍流．加之前述的原因导致近壁区成为了湍流楼型理论中的 薄弱环节。 在近壁区，由于强烈的粘性作用流动处于准层流状态，小尺度湍流脉动被抑制， 端流运动是一种火尺度相干结构．研究结果表明这种湍流运动是以波的形式存在fl，2】。 根据流动稳定性理论，采埘不稳定波模型对揣流相干结构的理论研究得到了实验的验 证，取得了许多有意义的成果【3，4，5】，近来采_[}{该理论模型对近壁湍流输运特性的研究 取得了新的进展『6．71。 鉴于上述分析．参照f6、7】的研究，提出在湍流边界层的近擘区采尉共振三波的理 论模型来描述湍流相干结构扰动速度，根据理论模型计算了占的分布，然后在传统k．s 模式基础上放弃求解5方程而代之以本文理论s构成封l=11。所得平均速度分布在近壁 粘性作埘层与Spalding剖面符台很女f。 1、理论分析 设扰动量由一个二维波和一对对称的二维波组成．有如下形式13。61 妒=蛾(y)P“甜一埘’彳+如(y弦‘‘2椰以“’+c．c (1) 对 向扰动速度。a，卢，甜分别为扰动波的流向、法向波数莆I频率。对h”彳=sjn(卢z) 涡域方科、连续方靴秘山，jb体过程可参址【3，51。 J一25 14时， Kachanov[8]的实验结果表明在人工引发的转捩区，当扰动波的总幅值达0 其基本波的幅值已达012，基本波代表了总扰动的特征，而且与T—s波形状一致。 Blackwelder的研究已证实了，湍流边界层近壁区的相干结构与转捩区的结构非常相似 【1】。因此，在我们的计算中可以取(2)式来近似的代表相干结构的扰动速度和压力， 其幅值参照实验结果确定。 文献[6，7】的研究也表明这样作是可以的。 将理论扰动量“，v，W代入5(以“，壁面摩擦速度，l，分子粘性系数无鼍纲化) 的表达式中，可得出以相干结构理论模型为基础的理论s表达式，如F： du6u 6Hau a“6u6uau 。