二维槽道湍流拟序结构的大涡模拟.pdfVIP

中国}一程热物理学会 燃烧学 第十届年会 论文编号:014073 二维槽道湍流拟序结构的大涡模拟 ‘ 向屏 郭印诚 北京 ·清华人学工程力学系 1‘00089) 联系电话(郭印诚)Email:guoyc画singhua.edu.cn 摘要;本文对不同入口边界条件的二维槽道流动进行了大涡模拟。采用Chorin的分步投 影法求解大尺度涡运动的Navier-Stokes方程，小尺度涡采用三种亚格子(SGS)模式分别 模拟，给出了不同亚格子祸档性模式下的模拟结果。对固壁面采用了壁函数。模拟结果 冉现了二维槽道流动拟序结构的发展演变过程 通过对不同入口速度下的瞬态流场的比 较.揭示了入口速度分布对流场的影响 关键词:大涡模拟:槽道流动:拟序结构:数值模拟 1.引言 槽道湍流流动是典型的壁湍流，是许多燃烧装置 如〔流化床锅炉)和化五设备 (如 催化裂化提升管反应器或下行床反应器)中经常遇到的流动情形。壁面附近的流动状态 对燃烧、传热及化学反应过程有显著的影响。本文以二维槽道流动为对象，采用大涡模 拟方法，探讨槽道流动中壁面附近的流动状态，为以后两相流动的大涡模拟奠定基础 研究表明 卜3]，与人们早期认为湍流是完全随机无序的认识不同，湍流中存在着大 尺度的拟序结构，这些大尺度的拟序结构在湍流的发展中起主导作用。近年来，数值模 拟成为全面了解流体动力学特性的有效途径。日前在实际工程中经常用到的统观数值模 拟方法，虽然能给出流动的宏观现象以及各种统计平均量，但对深入揭示流场的瞬态过 程却无能为力。直接数值模拟(DNS)求解了流动中的所有尺度，无疑可以得到流场的瞬 态过程，但目前的计算机性能使得DNS仍然只能局限于相当低Reynolds数和规则}L何 形状的流动问题。相比之下，大涡模拟以其准确性和可以接受的计算量成为模拟湍流流 动的有效方法。大涡模拟((LES)的基本思想是对流动的大尺度涡运动进行直接求解 Navier-Stokes方程，对小尺度涡通过亚格子模型进行模拟。在采用适当的亚格子模式的 情况F，大涡模拟结果的准确性很高。对于较高Reynolds数的湍流流动，大涡模拟的模 拟结果是真实的瞬态流场，有助于深入认识湍流流动的物理本质。 本文采用大涡模拟的方法对二维槽道湍流流动进行了数值模拟，模拟结果给出了流 场中拟序结构的瞬态发展演变过程。通过对不同入口速度下瞬态流场的比较，揭示了入 口速度分布对流场的影响。通过二种亚格子模式 卜摸拟结果的比较，揭示了亚格子模式 对流场描述准确性的影响。以期深入理解槽道流动的动力学特征和涡结构的作用，并为 ’国家自然科学基金 (编资助项日 253 以后开展两相流动过程的三维LES奠定基础 2.数值模拟的对象和控制方程 2.1控制方程 本文的模拟对象为图1所示的无限大平板间的二维槽道流动 采用比度基准ba=I.OX10气二，主流速度U,=10.0111/s初始时刻入口处的流体物性 参数作为无量纲化的基本参量 见〔表一)，对直角坐标系 卜的_维槽道流动的质量连续 方程和Navier-Stokes方程进行滤波和无量纲化，得到大涡模拟的基本控制方程: 币,oPitcPv_。 连续性方程 代二丁 甲 一甲二丁 ， -二二一 一 U 己 OX 即 一 夕 终 *:cu,)=_鱼一耳G;1{I,!十G!!,I 沙 动量方程 o‘f,/ GT:,at-L}厉 e-,) P - - 韶 _act[PW11 一 厉 标志物输运方程 OX少QY,LScl既/」