大涡模拟中亚格子尺度湍动能与亚格子尺度耗散的时均统计特性的研究.pdfVIP

中国T程热物理学会 燃烧学 2003年学术会议论文 编呼：034083 大涡模拟中亚格子尺度湍动能 和亚格子尺度耗散的时均统计特性的研究 王兵，张会强，王希麟 清华大学工程力学系，北京100084 摘要：本文用大涡模拟方法研究了后台阶湍流流动的时均统计特性。大涡模拟给出了准确的流场，井 和此皋础I‘研究rⅡ格予湍动能和亚格予湍耗散的}f寸均特性。删究给Ⅲr后台阶湍流流动的流场结构 以发流动的册格于湍动能和亚格了耗散的时均统计结果的会场分布蚓。对比研究，流动的小J叫区域内 时均统汁速度、籼性和亚格了湍动能和弧格了湍耗散，分r耗散等耋占果。研究表I州，弧格了八度湍动 能和砸格子湍耗散随着流动在空间的发展而减弱，曲l流区内亚格了湍动能和耗散较弱；在台阶截面(y／H =I处)Ⅱ格了湍动能和耗散最大；亚格子湍动能小于脉动动能。牟场的分子耗散占很小部分。 关键词：破格了湍动能，亚格予湍耗散，时均统计．人涡模拟 1引言 人涡模拟(LES)目前是数值研究湍流流动的有效l：具之一，在人满模拟中对人尺 采川适当的亚格子模式的情况F，人涡模拟结果的准确性很高。对于较高Reynolds数的 湍流流动，大涡模拟计算量比直接模拟(DNS)小很多，但模拟结果能够得到真实的瞬 态流场，有助于深入认识湍流流动的物理本质。 在人涡模拟中大部分亚格子模式的假设是基丁涡粘性的，小尺度运动被认为是各向 均匀同性的。而且，目前人部分亚格子模式假发亚格子尺度的作用是从人尺度携带能鼙， 并通过粘性力耗散掉。然而，近年柴的研究也有表明平均地整体地看能量是从人尺度输 运剑小几度，但同样能量的逆过程，tld,尺度向人尺度的能量反射也是发生的。但是广 泛采用的Smagorinsky模式是绝对耗散的，即不能提供小尺度的能簧反射。但对丁预报 小尺度能量的产生和耗散是上E确的。 研究咂格子尺度的湍动能产生和能姑耗散对丁研究湍流流动机理，深入揭示湍流耗 散机理，为统观模型提供更多依据具有重要意义。本文以后台阶流动为对象，州人满模 拟方法一方面研究了亚格子湍动能和Ⅱ格子湍耗散的时均统计特性，另一方面也对比了 时均粘性与统观模型预估的粘性结果。 2数值方法 2．1控制方程 人涡模拟中大尺度的量(“”)由滤波函数定义，本文滤波函数G采J}{j网格滤波。 滤波函数f(x)=If(x。)O(x，x’)dx。 (1) 对方程无量纲化和滤波操作后，连续性方程和NS方程如卜， 冈家自然科学筚会资助项目(N 210 却 (2) 互％ga 巫％ 』p互％ 旦蠢、 盟％ 甄一熟 、厂 (3) 方程(2)、(3)控制大尺度量的运动。小尺度苗的影响出现在滤波产生的勺中，即 亚格子尺度应力。亚格子应力h=i矛瓦百，是大涡模拟中最重要的亚尺度项，需要用 砸格子模式来模拟该应力项。本文采用Smagorinsky模式，定义 叶；(G△)2I季l=(G△)2√2i忑 (4) 那么亚格子尺度应力模拟为： 勺=一c,2△2同写 (5) 其中 瓦=÷c要+争 ㈤ 模型中的常数C。取为0．17。 2．2计算方法 入口条件才用速度的1／7律给山剖面分布。山口条中卜采Hj 以使人尺度涡能够自由通过出流界面，而不产生1F物理反射和不对内流场造成干扰。肇 面上采用无滑移条件。 采川分步投影法对速度场和压力场解耦。对不含压力项的动量方