上盖板拖动流的大涡模拟.pdfVIP

核动力工程 2006年10月 NuclearPower Oct．2OO6 Engineering 上盖板拖动流的大涡模拟 魏泳涛，于建华，刘玉君 (四川大学建筑与环境学院工程力学系．成都．610065) 摘要：论文将SU／PG迎风修正和大涡模拟方法相结合，建立了湍流数值模拟的速度．压力有限元格式，实现了上 盖板拖动流的在高Reynolds数下的流动分析，并完成了流场流线和压力等值线的绘制。数值算例验证了有限元大涡模 拟在湍流分析中的有效性。 ． 关键词：有限单元法；SUNG迎风修正；大涡模拟；方腔上盖板拖动流 1前言 大涡模拟是近年来得到迅速发展的湍流数值 盯=一F，+2乒力 (1-3) 式(1)中，有上波浪线的变量是滤波后的大尺度 模拟方法，该方法哇ISmagodnsky于1963年提出。 大涡模拟来源于这样的事实，即在湍流运动中的 量。P是流体密度，kg·m。3；露为速度向量， 大尺度涡(大涡)对于平均流动有明显的影响，其 m．S一；石是体力向量，m．S．2；孑和r分别是 运动可通过直接求解Navier-Stokes方程获得；而Couchy应力张量和亚格子应力张量， 小尺度涡(／J、涡)主要通过大涡之间的非线性相互 kg．m～．s。2；V是向量微分算子，!tiq；芦表示 作用间接产生，并通过粘性耗散并对大涡产生非 压力，kg．m～．S。2；∥为流体动力粘性系数， 线性作用；由于小涡近似存在各向同性性质，因 kg．m～．s一：J是二阶单位张量； 此可通过建立模型来处理。现在，大涡模拟已被 西=0．5[V石+历v】为形变速率张量，S～。 认为是介于直接数值模拟和湍流封闭模型理论之 r体现了小尺度运动对大尺度运动的影响， 间的一种湍流高级数值模拟方法。 对于这项，通过建立涡粘性模型来模拟 本文将大涡模拟方法与SU／PG迎风修正相结 (2) r=2鸬西=2万(e△)2【2西：西r2西 合，利用有限单元法实现了对方腔中的上盖板拖 式中，A=‰咄出，弦滤波宽度，m以；岛为亚 动流在高Reynolds数下的数值模拟。引入SU／PG 迎风修正的目的在于抑制因对流项而引发的数值 格子涡粘模式OPSmagorinsky模式常数，对各向同 振荡。 性湍流取Cs=0．1。 关于大涡模拟更进一步的知识见文献[1】。 2湍流的大涡模拟技术 2．3有限元离散 2．1 Navier-Stokes方程的滤波处理 利用Galerkin法可导出大涡运动方程(1)有限 元离散形式 大涡模拟通过对Navier-Stokes进行滤波处理 以将湍流运动分解成大尺度运动和小尺度运动。 掣l=姒对udf2+l驱pwl勖硷+ 而对张量函数的滤波处理，就是在空间中通过滤 n n (3) 波函数对其进行卷积积分。常用的滤波函数有盒