离心泵叶轮和蜗壳的非定常数值模拟marcgugau德国达姆施塔特工业.pdfVIP

离心泵叶轮和蜗壳的非定常数值模拟 Marc Gugau 德国达姆施塔特工业大学叶轮机械及流体动力系 1 引言 目前 CFD技术已经在叶轮机械中得到了广泛应用。CFD数值模拟在解决此类工程时显示 出了比实验更快更经济的优势。然而，在模拟叶轮机械非定常流动时需要研究更合适的方法 进行计算。准定常“冻结转子法”和非定常的“滑移网格”技术已经应用于计算旋转叶轮和 静止螺旋型蜗壳。在预测设计点的性能和叶轮和蜗壳的耦合关系基础上比较上述两种方法的 计算结果。由于两种方法产生不同的结果，特别在模拟叶轮和蜗壳的耦合关系时，人们不得 不重新研究瞬态时的计算方法的稳定性和适应性问题。 本文的研究了比转速为 19 的离心泵的流场特征。图 1 为作者在 Darmstadt 大学叶轮机 械和流体动力实验室做的实验离心泵的模型。目的是通过实验测量每部分的效率，从而找到 影响泵效率的主要因素。另外,为了能够进行数值模拟需要采取一些特殊的结构特征。首先， 叶轮和蜗壳壁面非常光滑，使壁面满足水力光滑管的要求，以确保近壁面采用壁面函数法时 对壁面的适应性.考虑到叶轮和蜗壳计算时二者之间没有间隙，而实际泵运行在间隙有流体 运动的状态下。数值模拟的基本结构参数如表 1和图 3。 2 具体模型 计算网格通过 fluent 前处理器 gambit,在近壁面附近用 CAD 软件把流动区域进行了分 割以产生若干内表面。这些措施有助于建立块结构化网格，为了构建内部体积，在 gambit 里添加了许多辅助的面结构。解决了螺旋形蜗壳的环行对称横截面和叶轮进口上游的进口管 处网格很难生成的难点。图 2 为计算区域的一些需要注意的网格区域。计算区域高质量的 hex网格容易进行区域定位，整个计算区域网格背离等体积的最大偏斜为 0.836,这是因为偏 斜在叶轮出口吸力面是不可避免的。整个区域的平均偏斜为 0.183 表明网格质量比较高,可 以不用四面体网格.使用块结构化网格的另一个原因是它有良好的后处理优势,而且可以用 另外的结构化网格求解器计算. 对于使用 fluent 求解器的用户,具有相似边界条件的几个块结构化网格能够合并成一 个整体的流体区域.六个流体区域组成了一个完整的离心泵计算区域,如图 3 分别命名为 F1 至 F6.在惯性坐标系中旋转的离心泵叶轮为流体区域 F2,在使用冻结转子法进行叶轮和蜗壳 耦合时计算采用旋转参考系。同样，非定常计算时采用旋转的动网格（即滑移网格技术）。 除此之外，其他流体区域是定常的，包括吸入管 F1,蜗壳螺旋型部分 F3,出口管 F4 延伸 2 倍长度的出口直径以消除出口边界条件对出口管流场的影响，同样也要考虑叶轮两侧的两个 间隙区域F5和F6.网格界面两侧对应一对有相对运动的流体区域,即F1对应F2, F2对应F3, F3 分别对应 F5 和 F6.表 2 为不同流体区域的网格统计和全部网格的统计数据.KORRELL[2] 曾进行了网格独立性研究,采用旋转参考系和周期性边界条件定义叶轮计算区域.采用三种 不同数量的网格方案,对应的是A网格节点数以及A/2和A/4的节点数,而泵的性能特性没有 发生比较大的变化.对于确定整个泵最佳的网格数,需要考虑网格数和无量纲壁面长度 y+与 壁面函数理论的一致性.选择最常用的标准湍流模型可以不考虑湍流模型的适应性要求,采 用标准k −ε 方程而在光滑近壁面采用标准壁面函数.进口速度边界条件按照给定的工况点 流量计算.湍流强度估计为 5%左右,长度因子为 DS/2mm.出口管出口处设置为压力出口条件, 采用二阶迎风格式进行求解区域离散,压力和速度耦合采用 SIMPLE 算法.按照两次间隔计算 的残差为收敛标准,计算时监视不同位置的质量流量和静压变化.采用冻结转子法时收敛需 要进行 1000 次迭代,然后将结果作为非定常求解时的初始条件.时间离散采用二阶欧拉隐格 式意味着对最大时间步长没有限制. 计算时一个时间步长进行 7 次迭代转子旋转 2 度, 上 述设置接近实际运行特征和计算时间.非定常收敛条件采用监视叶轮流道内的质量流量,进 口压力以及叶轮和出口的动量矩随时间的变化曲线. 从这些残差曲线可以检查一个时间周 期的求解过程. 3 性能结果和分析 从应用的角度来说,最值得注意的结果应该是泵的性能曲线.可以预测泵的性能以及确 定泵的一部分参数的取值等等而不需要知道关于泵内部流动的具体特征. 实验