旋转圆管湍流大涡模拟数值的研究.pdfVIP

第 34 卷 第 6 期 力 学 学 报 、 1．34．No．6 2002 年 l1 月 ACTA M ECHANICA SINICA NOV．，2002 旋转圆管湍流的大涡模拟数值研究 刘难生 仲峰泉 陆夕云 庄礼贤 (中131科学技术大学力学和机械工程系，合肥 230026) 摘要 采用大涡模拟 (LES)方法，并结合动力学亚格子尺度应力 (SGS)模型，通过数值求解柱 坐标系下的滤波 Navier-Stokes方程，研究了绕管轴旋转 圆管内的湍流流动特性．为验证计算 的可靠性，以及动力学 SGS模型对于旋转湍流的适用性，将大涡模拟计算所得的结果，与相 应 的直接模拟 (DNS)结果和实验数据进行了对比验证，吻合 良好．进一步对旋转圆管湍流的 物理机理进行了探讨，研究 了湍流特性随旋转速率的变化规律．当旋转速率增加时，湍流流动 有层流化的发展趋势．基于湍动能变化的关系，分析了旋转效应对湍流脉动生成的抑制作用． 关键词 大涡模拟 (LES)，亚格子尺度 (SGS)模型，旋转圆管湍流 引 言 旋转湍流在机械工程领域中具有广泛的应用背景，也是湍流基础研究中一个备受关注的重 要课题．众所周知，旋转流动对于燃烧系统中的火焰稳定有非常重要的影响，并且直接关系到 燃烧效率的提高．在大气和海洋科学的研究领域中，与地球 自转效应关联的旋转湍流是普遍存 在的物理现象，并且是季风和潮汐流动形成的重要因素．而绕轴向旋转的圆管湍流存在流 向剪 切和周向旋转效应的耦合作用，常被作为研究旋转湍流物理机制的一个典型问题．因为旋转效 应的影响，这种流动比纯剪切湍流更为复杂． 关于绕轴向旋转圆管内的湍流流动，主控流动特性的无量纲参数有两个：一个是 Reynolds 数，定义为Re=2ubR／ ，其中 为管截面上的平均速度， R为圆管半径， 为流体运动黏 性系数；另一个则是旋转速率 Ⅳ ，定义为 N ： R／ub，其中 表示圆管的旋转角速度．另 外，在地转运动问题中，后者通常采用 Rossby数 Ro(Ro=ub／MR)作为无量纲参数． 在充分发展的旋转圆管湍流中，由于旋转效应的影响，使壁面附近的流动趋于稳定 1【】．这 种效应抑制了近壁区的湍流脉动，在快速旋转时，会使得流向平均速度剖面分布规律趋向于层 流 PoiseuiUe速度分布．对于旋转圆管湍流的机理研究，前人 已经做 了一些工作．为了描述 出旋 转湍流的旋转效应，Re ich等人 0【】和 Bradshaw[3J在他们各 自的实验和理论研究工作中，还用 浮力来比拟旋转运动产生的离心力，将旋转 Richardson数引入到混合长理论中，研究了旋转效 应对湍流的抑制作用． wllite4【J和 Schehukin[5】通过实验也证实随着旋转速率的增大，湍流的 层流化的趋势会变得更加明显，主要表现为流向速度剖面更趋向于层流解，同时湍流强度在壁 面附近减弱．数值研究方面，Orlandi等人 6【J采用直接数值模拟方法，研究了较低雷诺数下的 旋转圆管湍流问题． Orlandi等人认为旋转效应使得在管壁附近的涡结构重新构造，尺寸及形 状都发生变化，从而导致近壁区的湍流喷射和下刮事件发生的几率减小，强度也会随旋转速率 2001-05—31收到第一稿 ， 2001-12-18收到修改稿． 1)国家自然科学基金 、中国科学院 “百人计划”和教育部 “跨世纪优秀人才培养计划”资助项 目 力 学 学 报 2002年 第 34 卷 的增大而减弱． Imao等人 【】在实验中观测到旋转圆管湍流中，周向平均速度剖面并不是刚体 式旋转．此外，传统的湍流模式理论研究表明，采用涡黏性模型的雷诺平均方法无法模拟出非 零的雷诺应力 p，因此无法数值再现出非零周向平均速度 IS】．Hirai等人 。【】的数值研究发现 标准的 —s模型对于旋转圆管湍流也不能给出满意的结果．所 以， Speziale等人 8【J指出只有 利用高阶雷诺应力模型，如显式代数应力模型 加【】和二阶封闭模型 n【】，雷诺平均方法才能合 理地数值预测旋转湍流特性． 近