关于“边界条件独立性原理”和“二阶矩湍流模型的改进”.pdfVIP

关于“边界条件独立性原理”和“二阶矩湍流模型的改进” 吕盘明 （中国科学技术大学工程科学学院，合肥，230026，pmlu@） 摘要 这是一个通俗一点的湍流模型发展史背景材料，同时也是作者历经近 20年的研究工作的 一个总结。它包括：1。二阶矩湍流模型的理论体系及其遇到的“高歌疑难”。2。作者提出的 改进二阶矩湍流模型使之能加入湍流结构的方向性特征的边界条件的的一整套理论（共十项内 容）。3。湍流应力张量边界条件的独立性原理及其对“高歌疑难”的解释。4。对后人应用本 模型的可能的工作方向的七条建议。 一．二阶矩湍流模型理论碰到三个难题（ “高歌疑难” ） 1．1895 年，O．Reynolds （雷诺） 推出时间平均的纳维尔-斯托克斯方程（N-S 方程）， ′ ′ 其中出现了一组新的变量--湍流应力张量现今也叫 Reynolds 应力张量,的九个分量 ， u−u i j 而要求解 Reynolds 应力张量，它们是脉动速度的二阶关联量，又需要知道脉动速度的三阶关 联量。如此下去，无有穷尽，因此为了方程组能封闭求解，唯有建立起 Reynolds应力张量与 平均速度场之间近似的关联，这就是所谓的湍流模型。差不多在同时，布辛涅斯克（Boussinesq， 1883）提出，可以仿照分子粘性应力张量和平均速度梯度场之间的关联，提出一个各向同性化 的涡粘性系数（ε 或ν ）和涡粘性模型（EVM，eddy-viscosity-model）的概念。 m t 2．1925 年，普朗特（L. Prandlt）提出混合长度假设（MLH，mixed-length-hypothesis）, 对边界层流动一类的简单剪切流可用经验的混合长度l 算出涡粘性系数,这一方法一直沿用 m 至今。但对内部复杂流动等复杂剪切流就不行了，因为误差太大。 3．1945 年，周培源 （P Y Chou）推出了Reynolds 应力张量的输运方程，但其中的 Reynolds 应力耗散项ε ，扩散项D 和脉动速度压力关联项Φ 这三项因不能精确计算出来需要用模型 ij ij ij 关系式近视表达出来。与此差不多同时，1940 年，柯尔莫果洛夫（Kolmogorov）提出了Reynolds 应力耗散项ε 的各向同性化模型，1951 年，Rotta提出了脉动速度压力关联项Φ 的各向同 ij ij 性化模型。 4．1967 年，英国帝国理工学院航空系的布兰德肖（P.Bradshaw）提出一个一方程（湍流 动能k 的输运方程）的湍流模型，这比零方程的混合长度模型有了很大进步。但还是只适合于 外流边界层流动。 5．1971-1975 年之间，英国帝国理工学院机械系的 D.B.Spalding 的小组发展出两方程涡 粘性模型 k EVM−ε （Launder,B.E.,Spalding,D.B.,1974）,代数 Reynolds 应力模型ASM （Rodi,W.,1976）,微分应力模型 （Launder,B.E.,Reece,G.L.,Rodi,W.,1975）等一套 DSM 通称为二阶矩湍流模型的理论,它们在一定程度上能考虑到历史的影响，因而可以用于复杂剪