壁湍流传热机理研究.pdfVIP

壁湍流传热机理的研究 张 哲 (天津大学力学系，天津300672) 摘要 本文主要研究了壁湍流的传热机理．其基本思想是确定壁湍流的相干结 构所引起的流场，小尺度湍流用统一的模型来表示。求解了壁湍流的三维温度 场，得到了表征热通量的努塞特无量纲量，与实验吻合．所解温度场具有条纹结 构，说明了数值模拟中所发现的条纹结构产生机理。 关键词相干结构，Nusself敖，标量输运 l 引言 剪切湍流中存在着相干结构，而相干结构对湍流的发生和维持起着关键作 用，这已是众所周知的事实．Nagano等人的实验证实了相干结构在边界层近壁 区起着决定性的作用Ⅱ5B．本文在壁湍流标璧输运问题中(壁湍流传热)采用相 干结构的理论模型，得到比通常的湍流计算更合理的摸型． 2壁湍流的相干结构模型 在1m。中采用共振三波模型来解释壁湍流的相干结构产生机理，得到了相 干结构附近流场的法向削面中沿某一事先给定的路径计算的净环量的概率密度 分布函数，与由分析直接数值模拟结果所得的概率密度分布函数很相近，间接证 明了这一模型是合理的． 这一模型的出发点是一个包含三个不稳定波的共振三波．其基本波的速度可 以表示为： dl“e’拙l‘+卢Ip。I‘’+口，￡^∥(a2x+p2㈨±’)I-口t“，e’m]Hp】=一。3‘’(1) 其中虬、“：、沁是由特征值问题解出的特征函数，已按一定条件正规化，因 而a，、a：、吒分别代表三个波的幅值，x、J，、g-分别为流向、法向和展向坐标， ％、声。(≈=l，2，3)是波数，珊々是频率．采用时间模式，口。、∥；、∞。满足特征 关系．共振条件是： 口l=口2+d3，∥l=∥2+岛，COl，=CO”+∞3r 在Ⅱ1目采用了弱非线性理论求a。、d：，a；随时间的演化规律，并按它们的 增长快慢以确定是否对应于～种可能的相干结构．详细论述请看＆1日，这里就 不重复了． 由于弱非线性理论本身的限制，我们不能指望用这一模型来全面反映相干结 构的流场． 但是，Blackwetder早已指出，湍流的相干结构和转捩区的流动结构很似 Ⅱ6Ⅱ．而Kachanov则给出了转捩初期的扰动波的各次谐波的凑度幅值分布 《7m．在Ⅱ7 D的图13中可以清楚地看到，当总扰动幅值达0．14时，其基本 波的幅值已达O，12，且基本波与总扰动波仅在峰值处稍有差别，丽在其它地方 则几乎可以完垒代表总扰动，而且基本波的形状和T．S波的形状基本一致．因此 在我们的计算中可以取(1)式所给出的速度分布形式来近似地代表相干结构的速 度分布形式，但它们的幅值则不能用弱非线性理论给出的演化规律来表达，丽必 须另外设法给定． ·8d8· 3湍流传热方程及数值弦算方法 不考虑辐射、热源及籼性引起的摩擦生热项，则不可压流体的付立叶能量方 程为： (3) 鹏【≥+(＆-v)-r】=V一(EVT) 其中P是密度，G是定容比热，丁是温度，K为分子热扩散系数，V是梯度 算子． 将速度分解为： “=U+《+“’ (4) 其中U是平均流速度，玎是由(1)式所给的速度(其中的参数如何逸择将在下 面给出)，“’则是小尺度湍流脉动速度。这里的一个基本假定就是小尺度湍流脉 动规律在不同的流动中具有一定的蔷适性，因而其作用可以用基本相同的公式描 述． 出于U及玎是已知的，用z』。=U+玎代表，于是(3)式可以改写为： (4) 心。(≥+(‰·v)r】+，)C。(弘v)r=V·(107趵 由于小尺度湍流的随机性，其作用类于分子混合的作用，所以可以将与其有 关的项写成似于分子传热的梯度形式，即设 pc。(“·v)ra-V·(KV乃(5) 其中K是湍流涡流热系数． 将方程无量纲化，令：