贯流风扇气动声场数值模拟及分析研究.pdf

中国工程热物理学会 热机气动热力学 学术会议论文 编号：082046 贯流风扇气动声场数值模拟及分析 刘敏 吴克启 (华中科技大学煤燃烧国家重点实验室,武汉,430074) (Tel: 02783,Email:zeliumin110@) 摘要:针对贯流风扇气动噪声传播特性，采用基于非结构网格的 CE/SE 算法对其气动流场进行数值模 拟，并直接模拟了气动噪声的传播。流场模拟采用大涡模型模拟，搭接式滑移网格模型处理动静干涉。 声场的传播采用完全欧拉方程作为控制方程。两相对旋转圆柱之间的流动模拟作为一个标准测试问题 与fluent模拟结果进行对比计算，表明本文算法准确性。随后将使用该方法得到了声场结果与实验数 据进行了对比分析，结果表明与实验测量得到宽带频谱不同，采用大涡模拟湍流模型进行声源模拟时， 会加强部分离散频率上的声级，从而产生误差。 关键词:贯流风扇；CE/SE算法；气动噪声；直接数值模拟；实验分析 0 前言 计算气动声学问题对数值算法提出了高精度，低频散和低耗散的要求。目前采用的数 值算法主要有离散关系保持格式（DRP ），紧致差分格式（HOC ），离散迦辽金方法（galerkin method ）以及最近的格子boltzmman方法(LBM)等。由于算法本身的限制，例如DRP,HOC 等最初建立于结构网格而难于应用于非结构网格，LBM等精度及耗散性欠缺，均难以用 于求解复杂流场的气动声学问题，例如贯流风扇的流场。尽管随后上述算法进行了优化， 例如Foluso将DRP和HOC应用于曲线网格[1]，Liu Min采用HOC及浸入式边界的方法求解 复杂固壁的声反射问题[2] ，Eric采用分区结构化网格及高阶差分格式求解线性欧拉方程 模拟声在复杂外型固壁下的散射现象等[3]，但在本质上这些方法应用于计算气动声学依 然存在一定的不足。计算气动声学的另一个挑战是边界条件的处理。其一是无反射边界 条件的应用，其二是在利用高阶差分格式时如何保证边界条件的高精度。 针对风扇等复杂流场气动噪声问题，yong采用大涡模拟模型求解不可压Navier-Stokes 方程求解出叶片及蜗舌表面处的压力脉动，采用FW-H方程对表面进行积分求解，忽略 了风扇涡壳的影响，即没有考虑声音的反射，衍射等因素[4] 。文献[2]对比了有无蜗壳存 在情况下的声场，结果表明蜗壳的存在对声场的结果存在较大的影响，主要是声反射， 衍射等因素，因而不能忽略蜗壳的影响。 CE/SE方法是NASA Lewis研究中心的S.C.Chang近年来提出的一种全新的高分辨率守恒 型方程计算格式，具有计算精度高和格式构造简单等优点。它在给出网格点物理量值的 同时给出其偏导数，在相同的基点下可以提高格式的精度，这一点类似于紧致差分格式。 它具有局部和整体的守恒特性，在计算区域的每个守恒元内，尤其在具有大梯度和不连 续区域内可以达到真正意义上的时空通量守恒，同时它具有边界处理简单，物理概念清 晰等优点[5]。其特点是精度及分辩率高。X.Y.Wan等在有激波存在的条件下对CE/SE方法 的精度进行了分析，结果表明采用一阶Tayler展开的CE/SE方法其精度与四阶精度的ENO 方法相当[6]。同时能使用少量网格表现一个激波，具有很高的分辨率。 目前CE/SE算法正逐渐应用于计算气动声学领域，并表现出一定的优越性。Ching Y.Loh随后应用该方法求解非定常N-S方程来模拟了圆柱绕流近场离散噪声，计算结果得 到的涡脱落主频与试验结果非常接近，表明了CE/SE算法能相比其它方法得到更精确的 [7] 绕流脱落涡频率 。 数值模拟具有运动边界的问题时，可以直接采用动网格模拟，但其精度不足，并且 需要大量的计算资源，目前文献主要采用滑移网格的方法，从而避免了每时间步上网格 的更新。保证滑移面的精度是保证全局计算精度的关键。由于旋转计算域和固定计算域 的存在，两个计算域之间的界面上的通量及原始变量的计算格外重要。为了精确的模拟 转子和静子之间的干涉，文献中已经发展了很多方法，例如混合平面法，该模型假设流 动在周向