湍流雷诺应力的本构关系α.pdfVIP

　 辑第 15 卷第 1 期　　　　　 水 动 力 学 研 究 与 进 展 　　　 . , . 15, . 1 A Ser A V o l N o 　　2000 年 3 月　　　　　　　 　　　　　 . , 2000 JOU RNAL O F H YDROD YNAM ICS M ar 文章编号:(2000) 湍流雷诺应力的本构关系 梁在潮, 　　陆　晶 (武汉水利电力大学, 武汉 430072) 　　摘　要: 　本文根据建立力学本构关系的本构公理和湍流运动特征, 分析了雷诺应力与平均 速度梯度之间的本构关系, 得出了该雷诺应力表达式不具备适应各种边界条件的通用性, 但在某 些湍流或某些湍流区域具有存在条件的结论。文中还给出了适用条件的判断准则和计算方法。 关　键　词: 　湍流; 雷诺应力; 本构关系 中图分类号: 　 357. 5　　　文献标识码: O A 1　前言 从 1895 年雷诺给出著名的雷诺方程以后, 雷诺应力的表达式一直是湍流研究的重点课 题之一。在一阶矩湍流模型中, 通常假定雷诺应力与平均速度梯度成比例, 即 U i U j 2 - u u = ( + ) - K ( 1 _ 1) i j t X j X i 3 ij 其中 为涡旋粘性系数; K = u u 2 是湍流动能。 t i j 工程计算中, 对二维流动更多采用 P randtl 公式, 即 2 d U d u = l ( 1 _ 2) dy dy 其中 l = ky 为混合长度; U 为平均速度; u 为脉动速度。一方面由于形式简单, 计算工作量小, 在某些流动中计算结果的精度满足工程的要求, 因而至今在工程计算中仍广泛地采用。另一方 面在理论上存在许多不合理的地方, 严格讲一般不能用平均速度梯度的形式来表示雷诺应力。 因此, 国内外出现了许多其它湍流模式, 而且为了更多的反映湍流的复杂因素, 模式做得越来 越复杂, 造成计算工作量大, 甚至工程界难于理解。有些模式由于引入的假定太多, 使得计算精 度并不比简单的涡旋粘性模式的高。从工程的要求来说, 只要计算精度能满足工程要求, 计算 模式越简单越好。因此, 有必要对雷诺应力表达式进行分析, 找出涡旋粘性模式的使用条件, 进 行必要的修正, 扩大其使用范围。这就是本文的目的。 2　本构关系简介及牛顿流体的本构方程 收稿日期: ( ) 作者简介: 梁在潮 1930～ , 男, 教授, 博导。 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 梁在潮等: 湍流雷诺应力的本构关系 9