平面突扩流动非稳定性的大涡模拟ξ.pdfVIP

第 15 卷第 2 期 计 算 力 学 学 报 V o l. 15 N o. 2 1998 年 5 月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS M ay 1998 平面突扩流动非稳定性的大涡模拟 刘有军 (大庆石油学院机械系, 安达, 黑龙江, 151400) 摘　要　用 T G 有限元法求解二维大涡模拟的非定常N S 方程, 数值模拟了平面突扩流动, 计算结 果展示了回流区内流动的多涡结构及其不稳定的周期振荡过程, 而该振荡过程的周期时间历程平 均, 恰好与诸多湍流模型的计算结果及实验结果相吻合, 证明了数值模拟的正确性。 关键词　突扩流动; 大涡模拟; 周期振荡; 时间平均 分类号　 211 3 V 1　引　言 　　传统上对湍流运动的认识是湍流可以分解为平均运动与脉动运动两部分运动之和, 而其 中的脉动是完全不规则的随机运动。对湍流运动进行分析或计算的主要方法就是所谓的雷诺 时间平均方法。平均运算的结果是将脉动运动的全部行为细节一律抹平, 丢失了包含在脉动运 动中的大量有意义的信息, 计算结果反映的是流动在时间历程上的统计平均表现, 而无法反映 流动的瞬时脉动特性。并且当采用涡粘性湍流模型来封闭雷诺时均方程时, 将湍流运动中所有 大小不同尺度的涡同等对待, 不加区分, 且认为都是各向同性的, 从而使模型的应用受到一定 的限制, 如不适用于诸如强旋流之类的各向异性问题。 [ 1 ] 一种非常有希望的湍流数值模拟方法是大涡模拟 , 它源于对湍流的直接数值模拟, 在目 前计算机所允许采用的计算网格尺度比湍流最小涡尺度大得多的情况下, 只能放弃对湍流全 部尺度范围上的涡都进行直接数值模拟的奢望。作为直接数值模拟的一种变通, 改为只将比网 格尺度大的大涡运动通过数值求解N S 方程直接计算出来, 而认为比网格尺度小的小涡是各 向同性的, 它对大尺度涡运动的影响, 可以通过建立模型来模拟, 这就是大涡模拟的基本思想。 为了使小涡的截断不影响能量由大涡向小涡的传输, 不能直接用N S 方程来表示去掉小涡后 的大涡的运动, 而必须经过某种滤波处理, 该过程类似于对N S 方程的时间平均, 但滤波是对 方程在空间上的处理, 所以也称为空间平均。大涡模拟与时间平均不同, 即使对于统计定 N S 常的流动, 滤波后的流动参数也不是常数, 它仍包含大涡的脉动, 且若滤波参数取得足够小, 则 滤波后的流动参数仍保留全部的脉动。所以, 大涡模拟更能反映流动的某些细节, 为我们更加 ( ) 深刻地认识湍流的本质, 特别是认识某些流动的非定常特性 时间平均无法反映 , 提供一些非 常有价值的东西。 原则上讲, 大涡模拟与直接数值模拟一样, 计算必须是三维和非定常的, 这意味着需花费 昂贵的计算代价。因此, 用大涡模拟法模拟实际流动的例子并不多。大涡模拟的开拓者D ear 收稿日期: 修改稿收到日期: 　刘有军: 男, 1965 年生, 博士 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 2 期　　　　　　　　　刘有军: 平面突扩流动非稳定性的大涡模拟　　　　　　　　　　 187 [ 1 ] do rff 于 19