明渠中跌坎后突扩分离流数值研究力学学报.PDFVIP

第 30 卷 第 1 期 力 　学 　学 　报 Vol. 30 , No , 1 1998 年 1 月 ACTA MECHAN ICA SIN ICA Jan. , 1998 明渠中跌坎后突扩分离流数值研究1) 刘沛清 　　邓学蓥 (北京航空航天大学流体力学所 , 北京 100083) 摘要 　利用大涡模拟技术 , 对明渠中跌坎后的二维突扩分离流结构进行了数值模拟 , 探讨了这类典 型分离流各特征区的流动性质 , 分析了再附区壁面上不同测点处脉动压力的统计特征 , 并与试验结 果进行了比较. 在数值模拟中 , 采用了弱压缩流的控制方程和非均匀网格系统. 关键词 　大涡模拟 , 分离流 , 跌坎 前 　言 跌坎后突扩分离流不仅是流体力学中复杂分离流的一个典型例子 , 同时也是许多工程领域 中一类常见的流动现象 , 如在水利工程中常见于闸堰下游的水流结构中. 对于这类复杂分离流 动 , 长期以来一直被人们视为研究的重点 , 目前在理论、实验和数值模拟等方面均取得较丰富 的成果. 在理论方面 , 文献 [1～3 ] 利用相似性理论分析了剪切层的发展过程 , 导出了速度分 布的相似公式 ; 在实验方面 , 人们利用流场显示技术、热线或热膜测速仪和激光测速仪 (一维 ) 或二维 , 详细地测量了坎后剪切层的分离和再附过程、时均和湍动场结构等 , 这方面资料较 ε 多 , 目前已成为各种湍流模式的验证资料; 在数值模拟方面 , 利用 K - 湍流模式和其它基于 雷诺时均概念的各种湍流各向异性模式 , 较成功地预报了这类分离流动时均场的主要特征. 然 而 , 由于流动结构的复杂性, 精细模拟湍动场结构目前仍处于探索阶段. 为此 , 本文利用基于子网格空间平均法的大涡模拟技术 , 数值研究了这类分离流结构 , 特 别是分析了再附区壁面上脉动压力场的主要特征参数 , 并与有关试验资料进行了比较. 1 　控制方程与大涡模拟技术 为便于显式时步推进压力场 , 用弱压缩流模型[4 ,5 ]作为控制方程 , 即 连续方程 p p 2 uj ρ ( ) t + uj xj + a xj = 0 1 动量方程 2 ui ui 1 p ui υ ( ) t + uj xj = - ρx i + x 2 2