圆截面u形管流有限元模拟-空气动力学学报.pdfVIP

第27卷第1期 空气动力学学报 V01．27。No．1 2009年02月 ACTAAERoDYNAMICASINICA Feb．，2009 文章编号：0258—1825{2009)Ol一0083—05 圆截面U形管流有限元模拟 樊洪明1，何开远1，李先庭2 (1．北京工业大学，北京100022；2．清华大学，北京100084) 摘耍：弯曲圆管内流动具有很强的工业应用背景，管内二次涡的大小、位置及发生的频率会对流动系统的阻力及 热质交换产生重大作用，研究弯曲圆管内的流动具有重要的理论意义。本文应用湍流大涡模拟方法对弯道流动进 行了数值模拟。用时问推进法，给出了隐含流线迎风耗散作用的二阶完全展开Taylor-Galerkin有限元离散格式， 数值计算的结果与实验结果符合较好。表明大涡模拟方法适用于边界形状复杂，存在各向异性的大尺度涡的流动 仿真。 关键词：湍流；大涡模拟；流场；有限元法；弯管流 中图分类号：0353 文献标识码：A 统计理论和拟序结构认识的基础上的一种新的数值 0 引 言 预测方法，它克服了湍流模式理论的时均处理和普适 性方面存在的缺陷。目前，高密度时空网格的LES 弯曲圆管内的流动在工业管道输运系统中大量 已被列入湍流高级数值模拟范畴[9]。本文拟以U形 存在。自1927年Dean[13首次提出圆形截面弯曲管 圆管流为例，采用湍流大涡模拟方法进行数值仿真。 道内存在二次流以来，众多研究者对此给予了格外关 Azzola等对U形圆管内流动进行了系统的实验研 注[2~7]，这不仅因为研究弯曲圆管内流动具有很强的 究[2]，这些实验结果可以作为检验数值仿真可靠性与 工业应用背景，而且在一定条件下这种流动的二次流 精度的依据。 结构会发生变化，产生旋涡结构分叉现象，因而研究 弯曲圆管内的流动还具有重要的理论意义。由于管 1大涡模拟基本方程 内二次涡的大小、位置及发生的频率会对流动系统的 大涡模拟的基本思想是运用空间滤波的方法将 阻力及热质交换产生重大作用，详细了解这些信息就 流动变量划分为大尺度量和小尺度量两部分。大尺 有可能通过调整区段的几何形状或其他参数来控制 度量可通过数值求解运动方程直接计算出来；小尺度 二次流，达到满意的工作状态。过去对于某些重要局 运动对大尺度运动的影响将在运动方程中表现为类 部障碍区段或部件经常是通过模型或原型实验来测 似于雷诺应力的应力项，称之为亚格子应力，通过建 定总流的特征参数，例如阻力系数、总传热系数和传 立亚格子应力模型来对其进行模化。 质系数等等，但是，实测不同工况下的三维流场，投入 大涡模拟中常用的滤波器主要有富式截断滤波 非常之大，使得研究单位或厂家望而却步。近年来·， 器、盒式滤波器和高斯滤波器三种，其中高斯滤波器 计算机技术