CFD软件在对旋式轴流风机数值模拟中的应用.docVIP

/ CFD软件在对旋式轴流风机数值模拟中的应用 摘要：综合分析介绍了CFD软件在叶轮机械数值模拟中的两个关键问题，即数值算法和湍流模型及其最近的一些研究进展，并以对旋式轴流风机的数值模拟为例对其应用进行了说明，最后对CFD在叶轮机械数值模拟中的应用及发展趋势进行了展望。 关键词：计算流体动力学（CFD）；叶轮机械；数值模拟；数值算法；湍流模型 引言 随着科学技术的进步和社会经济的发展，人们节能意识逐步增强，许多领域对叶轮机械性能的要求越来越高。传统的设计方法需要进行试制和测量大量试验参数等工作，在对叶轮机械、喷管及管道等内部流动进行试验测量时，要求试验装置复杂庞大，从而消耗大量的人力和财力，成本较高，对试验的依赖性较强，研制周期较长。而数值模拟的方法将理论分析与试验研究联系在一起，以其独特的优势逐渐成为研究流体流动的重要手段。随着计算机内存和并行技术的发展，计算流体动力学（CFD）已经广泛应用于叶轮机械的研发过程中，并成为一门独特的学科。它能够描述复杂几何体内部的三维流动现象，可以在设计的初期快速地评价设计并做出修改，而不需要费原型生产和反复测试；在设计的中期，用来研究设计变化对流动的影响，减少未预料到的负面影响；设计完成后，CFD提供各种数据和图像，证实设计目的。近年来，CFD越来越多地应用于叶轮机械的设计和流场的分析中，大大减少了研发费用、时间及新设计带来的风险，成为一种重要的设计和计算方法[1]。 按照流体动力学解决问题的思路，一次成功完整的数值计算必须包括建立计算模型、生成网格、确定控制方程、选择湍流模型、确定初始和边界条件、确定数值算法、离散方法及求解方程等步骤，而CFD应用中的关键问题就是数值算法的应用及湍流模型的选取。 1数值算法 20世纪90年代之前，受到计算机技术的限制，叶轮机械内部流动的数值模拟从无粘势流和线性化处理阶段逐步向综合考虑内流粘性和回流方向的准粘性模拟阶段发展[2]。早期叶轮机械内流计算通常简化为二维不可压势流或三维势流，以势函数、流函数或Euler方程为控制方程进行求解。50年代将三维Euler方程简化为两个相互关联的二维方程求解，在当时技术条件下为叶轮机械内部流场数值计算建立了理论基础。从1952年吴仲华教授提出S1、S2流面理论以来，人们普遍采用S1、S2流面相互迭代的方法来计算叶轮内部流动，并由此产生了流线曲率法和准正交面法等一些数值方法。由于非粘性假设在一定程度上可以反映实际流动情况，与物理过程接近而且从认识上易被接受，因此目前仍有不少学者采用两类流面理论研究叶轮机械内部流动，尤其在叶轮机械水力设计反问题中应用很广泛。之后随着计算机技术的发展又出现了更为复杂的计算方法，使得叶轮机械的数值模拟不再停留在无粘阶段，开始综合考虑粘性、回流以及漩涡对内流的影响，出现了势流—边界层迭代解法、射流—尾流模型及涡量—流函数法等。其中势流—边界层迭代解法将叶轮机械内的流场分为无粘性的势流区和有粘性的湍流边界层区，分别进行计算并迭代。该方法对于叶轮机械内部存在的漩涡、二次流、脱流及叶尖间隙损失等复杂流动来讲是一种较好的解决方法。对于该方法在叶轮机械内流计算中的应用仍是当今的研究热点。 20世纪90年代以后，CFD技术随着大容量、高速度计算机的出现得到迅速发展，进入了一个全三维粘性数值模拟时期，即粘性、时间平均化处理阶段[2]。通过直接求解雷诺时均方程，结合湍流模型来计算叶轮机械内部的三维粘性流动成为叶轮机械数值模拟的主要方法。在离散方法上出现了有限差分法、有限体积法、有限元法、有限分析法、边界元法和谱方法等；在参数解耦方式上出现了压力修正法、时间相关法、拟可压缩法（人工压缩性法）和松弛迭代法等一系列研究成果。1972年问世的SIMPLE算法是压力修正法的典型代表，之后又出现了SIMPLER（Patankar,1979）、SIMPLEST（Spalding,1981）、SIMPLEC（DoormalRaithby,1983）等一系列改进方案[3]。SIMPLE算法及其改进算法仍是目前求解叶轮机械内部不可压流动的重要算法。许多学者采用压力修正法来求解叶轮机械内部流场。现今对相关算法的改进算法仍是许多学者的研究课题。时间相关法（时间推进法）也是同一时期出现的数值模拟方法，除了适用于低亚声速、亚声速、跨声速和超声速流动外，还可同时用于内流和外流、定常和非定常流场的计算。该方 / 法可分为显示和隐式两类。叶轮机械内计算方面应用较广的显示格式有：Lax-Wendroff（L-W）格式、MacCormack预估修正格式及Runge-Kutta格式。对隐式方法的求解有近似因式分解法（AF法）和迎风格式。 2湍流模型 所谓湍流模型，就是建立湍流脉动附加项与时均量之间的关系，从而使控制流动的方程组能够封闭。一个良好的湍流模