弯曲管段内流动的大涡模拟ξ.pdfVIP

　 辑第 16 卷第 1 期　　　　　 水 动 力 学 研 究 与 进 展 　　　 . , . 16, . 1 A Ser A V ol N o 　　2001 年 3 月　　　　　　　 　　　　　 . , 2001 JOU RNAL O F H YDROD YNAM ICS M ar 文章编号:(2001) 弯曲管段内流动的大涡模拟 1 2 2 樊洪明 , 　　何钟怡 , 　　王小华 ( 1. 清华大学, 北京 100084; 2. 哈尔滨工业大学, 黑龙江哈尔滨 150006) 　　摘　要: 　本文应用湍流大涡模拟方法对三维空间的弯管流进行了数值模拟。计算中使用了 隐含流线迎风耗散作用的具有三阶部分展开的 T aylor_ Galerk in 离散格式, 导出了经 Gauss 滤波 后的基本方程组的有限元列式。数值计算结果与实验结果符合较好, 表明具有适当计算格式的湍 流大涡模拟方法适用于边界形状复杂和存在各向异性大尺度涡的流动仿真。 关　键　词: 　湍流; 大涡模拟; 流场; 有限元法; 弯管流 中图分类号: 　 351. 3　　　文献标识码: O A 1　引言 在河道以及工业管道输运系统中, 除了长直区段外, 还大量存在几何形状不规则的区段, 这些区段内的流体一般都处于湍流状态, 运动十分复杂, 往往会诱导出二次流, 产生尺度很大 的旋涡, 造成局部障碍。随着对湍流机理认识的不断深入, 发现大涡与平均流动之间有强烈的 相互作用, 涡的发生和发展对边界几何形状以及边界条件的类型有强烈的依赖性, 其形态与强 度因流动的不同而不同, 是非均匀和高度各向异性的。大涡对平均运动也有强烈的影响, 大部 分质量、动量和能量的输运都是由大涡引起的。因此, 二次涡的大小、位置及发生的频率会对流 动系统的阻力及热质交换产生重大作用, 详细了解这些信息就有可能通过调整区段的几何形 状或其他参数来控制二次流, 达到满意的工作状态。过去对于某些重要局部障碍区段或部件经 常是通过模型或原型实验来测定总流的特征参数, 例如阻力系数、总传热系数和传质系数等 等, 但是, 实测不同工况下的三维流场, 投入非常之大, 使得研究单位或厂家望而却步。近年来, 计算机技术的飞速发展以及湍流数值模拟理论的不断改进, 通过数值仿真逐步代替物理实验 来提供场信息开始有了希望。八十年代以来的数值仿真结果表明, 由Smogrinsky, D eardoff [ 1, 2 ] 等针对气象预报问题提出的大涡模拟 (L ES) 理论, 特别适合于边界形状复杂, 存在各向异性 大尺度涡的内流情况, 这不仅包括有压流动系统中的不动部件区段, 也包括转动或移动部件等 等。大涡模拟是建立在湍流统计理论和拟序结构认识的基础上的一种新的数值预测方法。它 克服了湍流模式理论的时均处理和普适性方面存在的缺陷。其基本思想是: 首先, 把包括脉动 在内的湍流瞬时运动通过某种滤波方法分解成大尺度运动和小尺度运动两部分, 大尺度量要 收稿日期: ( ) ( ) 作者简介: 　樊洪明 1964～ , 男, 博士 博士后 。 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 樊洪明等: 弯曲管段内流动的大涡模拟