应急抢险移动泵站用轻载型轴流泵的设计研发

1、应急抢险移动泵站用轻载型轴流泵的设计研发分类号密级题冃：英文并列题目：研宂生：研究方向：导师姓名：导师姓名：学位授予日期：独创性声明本人声明所呈交的学位论文是本人在导师指导下进行的研究工作及取得的研究成果。尽我所知，除了文中特别加以标注和致谢的地方外，论文中不包含其他人已经发表或撰写过的研究成果，也不包含本人为获得江南大学或其它教育机构的学位或证书而使用过的材料。与我一同工作的同志对木研宂所做的任何贡献均已在论文中作了明确的说明并表示谢意。签名：円期：关于论文使用授权的说明本学位论文作者完全了解江南大学有关保留、使用学位论文的规定:江南大学有权保留并向国家有关部门或机构送交论文的复印件和磁盘，

2、允 许论文被查阅和借阅，可以将学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存、汇编学位论文，并且本人电子文档的内容和纸质论文的a容相一致。保密的学位论文在解密后也遵守此规定。签名：导师签名：日期：江南大学硕士学位论文摘要移动泵站具有灵活方便、随到随用的特点，是抗洪排涝抢险救灾的关键装备。轴流泵是应急抢险移动泵站中的最核心设备，既要保留一般轴流泵大流量的特点，又要提高扬程压力，确保具有远程输送能力(尤其在排除城市内涝时)。本文基于抢险移动泵站对轴流泵的特殊需求，以轴流泵 内部复杂的流场机理入手，利用计算流体力学(cfd)技术对给定的轴流泵进 行了多工况全流

3、道的三维湍流模拟计算，并进行优化设计。本文介绍了轴流泵流场模拟计算分析的理论基础，分析了所采用的数值计算方法，给出了其方程模型和求解方法。本文应用rans方程、sa紊流模型和numeca软件对给定的轴流泵进行了网格划分和流场模拟计算，计算结果表明原模型的效率较高，但在不同工况流道内存在人小不一的涡旋、二次回流等复杂的三维湍流现象，详细给出了其内部的流线图，直观地给出了泵内流场流态；对轴流泵的不同叶轮参数对轴流泵性能进行了研究：研究了叶片数对轴流泵性能的影响并得到了其性能曲线；研究了叶顶间隙对轴流泵性能的影响，得到了其性能曲线和间隙对叶轮压力分布的影响等；计算分析了有无导叶两种情况下的流场。在轴

4、流泵叶轮、导叶叶型的优化中，叶轮的叶片叶型采用独立压力面、 吸力面型线的优化方法，导叶的叶片叶型采用屮弧线和厚度分布型线优化方法，使优化的轴流泵模型性能增加明显，结果表明：运用numeca/design3d可很好地优化轴流泵，为轴流泵优化设计提供了新的方法，具有重要的实际价值。课题在无锡锡泵制造有限公司的实验室平台上对原型泵进行了性能测量，测量的结果证明了数值模拟的准确性，为轴流泵的叶轮、导叶叶片型线的设计提供了理论依据和实验佐证。（英文部分相应修改）关键词：移动泵站;轴流泵，流场模拟；计算流体力学；叶型优化；numeca； design3d iabstractabstractthe mobi

5、le station has the characteristics of flexible and convenient, andit can be used at anytime and everywhere z it compensates for the lack of afixed pump station widely .while the axial flow pump for emergency mobilepumping station is the core components.lt is necessary to further study andimprovement

6、 of its optimization design theory.ln order to understand thecomplex internal flow field of axial flow pump mechanism,use ofcomputational fluid dynamics cfd technique for a given axial flow pump areworking for full 3d turbulent numerical simulation.introduction of the axial flow pump flow field simu

7、lation analysis of theequations of the model, and introduces the method, in this paper theapplication of rans equation and sa model, by using the commercial softwarenumeca for a given axial flow pump are grid division and flow field simulationcalculation, calculation the results show that the model

8、has higher efficiency,but in different working flow exists in the size of the scroll, the two returncomplex three-dimensional turbulence, introduced in detail the internal flowchart, intuitive gives pump inner flow field of axial flow pump; impellerparameters on the performance of axial flow pump fo

9、r research and study:blade number on performance of axial flow pump and the impact theperformance curve of the tip clearance, the performance of axial flow pump,the performance curve and clearance of impeller pressure, calculated andanalysed without guide vanes under two conditions of flow field.det

10、ails of the axial flow pump impeller, vane blade optimization method,impeller vane using independent pressure surface, suction surface line ofoptimization methods, vane blade adopts the camber line and the thicknessdistribution lines optimization method. optimal axial flow pump modelperformance incr

11、eased significantly, the result shows that: the use ofnumeca/design3d can be a very good optimization of axial flow pump axialflow pump, it presents a new method for optimal design, has importantpractical value.ii江南大学硕士学位论文目录摘abstractii录1第一章绪论."11.1论文研宄的背> 11.1.1移动泵站的优点11.1.2轴流泵简11.1.3 cfd 的

12、应ffl 21.2国内外研宂现31.2.1国外研究现31.2.2国内研宄现护：31.3本文的研宄内容和方法4第二章流场数值模拟的基本理论和方法52.1数值计算的基本方矛呈52.1.1连续方f呈52.1.2动量方程62.1.3紊流时均方程62.2揣流数值计算的方法62.2.1非直接数值模拟与直接数值模m 72.2.2火涡模拟错误未定义书签。2.2.3 reynolds 平均法72.3三维紊流模型82.3.1标准k-e模型错误!未定义书签。2.3.2 rng k- e 模型错误!未定义书签。2.3.3 reazalizable k- e模型错误!未定义书签。2.3.4 spalart-allmar

13、a (sa)揣流模82.4控制方程组的求解方法92.5计算软件介会召9第三章轴流泵的全流道流场数值模m123.1轴流泵模型的三维建模123.2轴流泵内部计算网格的生13 i江南大学硕士学位论文3.3计算参数的设163.3.1边界条件的给定163.3.2定义求解控制（计算参数设 ®183.3.3计算求解183.4计算结果分丰斤183.4.1设计工况下的内流场分183.4.2不同工况流线速度的比较243.4.3叶片数对轴流泵性能的影u|nj 283.4.4不同叶顶间隙对轴流泵性能的影313.4.5有无导叶对轴流泵性能的影ho 333.5本章小结34量绍法37析结第四章移动泵站用轴流泵的实

14、验测4.1试验台与试验装置介4.2实验方35354.3实验结果及误差分394.4本章小42第五章轴流泵叶型参数化及其优435.1轴流泵叶轮叶型的优435.1.1轴流泵叶轮的参数化构成(fitting) 435.1.2轴流泵叶轮叶片数据库的生成(datebase)455.1.3轴流泵叶轮叶型的最终优化(optimization) 475.2轴流泵导叶叶型的优495.2.1轴流泵导叶的参数化构成(fitting) 495.2.2轴流泵导叶叶片数据库的生成(datebase) 515.2.3优化方法和过程525.2.4.导叶叶片优化的结525.3结果分析535.4本章小53第六章全文总结与展6.1全

15、文总结546.2展望.54致i射56参考文献57附录：作者在攻读硕士学位期间发表的论文602第一章绪论第一章绪论1.1论文研宂的背景1.1.1移动泵站的优点常规的抗旱排涝都采用固定式泵站来解决这些问题，但在超出固定泵站设计能力或特定条件下（如突发急雨造成局部内涝），固定式泵站就显出局限性1;在水利、城市给排水的固定泵站中一般采用低扬程、大流量 的轴流泵。移动式泵站成为目前政府弥补固定泵站局限性，解决突发抗旱 排涝应急预案的关键技术手段，发挥了不可替代的作用。在较高压力、较低流量的移动泵站（如消防车、道路清洗车）屮一般采用离心泵，扬程较高，但流量偏小，很少能达到500m3/h以上，无法适应大流量

16、抢险救灾的需求，尤其是在城市内涝抢险中直接利用汽车发动机动力、大流量（快速排水）和较高扬程（适应地下道、涵洞等至下水道或河距离较长）为此类移动泵站用泵的主要技术特点。因此，对应急抢险移动泵站用轻载型轴流泵与常规轴流泵相比，必须具有效率高、扬程高的特点，研究轴流泵内叶轮导叶的流场机理，找出其中影响效率、扬程等性能参数的相关因素并优 化，具有很高的理论意义和实际价值2。1.1.2轴流泵简介在叶片泵中，轴流泵的比转数较高，一般情况下在4501550之间3。随着对轴流泵的不断研究，近年来轴流泵的比转速呈现不断提高的趋势，国内己经研制出比转数在2000左右的模型泵，而国外己研制高达2500比转数的泵。轴

17、流泵扬程在110m之间，流量一般介于0.5060m3/s之间,正由于这些特点，轴流泵在国民生产中得到了普遍的应用。在不冋的场合大量使用丫不同型号的轴流泵：如中小型轴流泵一般用在应急抢险和农田浇灌；大型轴流泵用在大型水利工程（如南水北调）和城市给排水等；国民经济各产业中均能见到轴流泵发挥着巨大的贡献，如：造船公司、钢铁 厂等生产过程中的给排水、电厂冷却汽轮机的冷却水4;还有在如造纸、 环境工程、医学工程、石油化工等领域无处不见各种各样轴流泵的身影。因为叶轮的叶片的截面和机翼是类似的，根据空气动力学的机翼升力理论，在叶轮高速旋转的过程中，叶片切削水体，使水体相对于叶轮叶片产生了一个高速绕流，它使水

18、体源源不断的从底部抽吸上来，从而使轴流泵具有输送液体的作用。轴流泵具有g己的性能特点：在小流量的工况下，轴流泵的扬程h流量q曲线有一个马鞍形凹下的部分，此凹点为关死点， 在该点处扬程值最高，差不多是最高效率点处扬程的2倍5,其功率也相 对很高，但其效率不高。造成轴流泵这个特点的主耍原因是：在轴流泵叶片角度不变时，流量变小时，在速度三角形中，水体的圆周速度和相对速度的夹角也随之减小，这就造成轴流泵的冲角变大，流量越小，冲击角度越大，当冲角达到一定值时6,水体不能继续吸附在叶栅的表妞，从而产 生了脱流。由于脱流产生，在轴流泵叶片的不同截面上出现了规模不等的回流，这种回流过程称为二次回流。在回流的过

19、程屮，轴流泵的出门处便产生了大小不等的回旋漩涡，这些漩涡使水体在轴流泵中的圆周速度部分转化成了沿着轴流泵轴心方向的速度，这使得轴流泵的扬程增加。但这些回流、漩涡等会使水体大部分的机械能转化成内能，造成很大的能量损1江南大学硕士学位论文失，从而导致轴流泵在小流量处的效率明显的降低。闵此在轴流泵的工作区域中，流量大时效率高而扬程较低，相反在流量小时效率低而扬程却高，整体而言效率高的工作区域相对不是很大。轴流泵虽然结构简单，属于叶轮旋转机械的一种，随着市场需求的不断发展，对轴流泵性能的要求也不断在提髙，迫切需要对轴流泵进行更深入的研究。轴流泵虽然具有流量大，扬程低的特点，但轴流泵a部复杂的流场流态对

20、轴流泵效率和流量等性能参数的影响规律还是需要更深层次 的研宂，因为轴流泵内部的三维湍流流动相当复杂，具有随机性和不稳定性，而且人们对轴流泵内部的流场流动机理还不是十分明朗，为了掌握轴流泵的动力性能的影响因素，积累关键的实验数椐和为轴流泵的设计原理的提供理论支持和为最终实现轴流泵的优化设计方法成为当今设计研究轴流泵的重要内容。随着社会的发展，研究轴流泵的力*法也越來越多，其屮最主要的方法有两种：实验测量和数值模拟，实验测量技术作为研究轴流泵最基本的技术，在近年来得到了不断的发展，其测量方法也是越来越先进，例如，图像测速仪等大量新的测量工具的产生；数值模拟技术作为近年来出现的一种新方法，在泵领域的

21、应用fi益增强，随着流体动力学(cfd)的发展，数值模拟技术己经广泛应用于模拟轴流泵内部的流场流态，预测轴流泵的性能，对不同参数的轴流泵进行比较，轴流泵叶型优化设计，叶轮/导叶相互干扰设计等方面。随着计算机技术和cfd技术的不断发展，人们的终极目标是实现数值模拟完全替代实体模型的实验，能够模拟一切的实际工况，但限于计算机的性能和cfd理论的支持，这一步还要很长的路要走。但当前，数值的模拟方法还是对轴流泵内部流场的研宄产生了巨大的作用。1.1.3 cfd的应用当前研宄轴流泵内流场的方法有很多，其中经验公式、数值模拟计算和实验测试分析三种最为广泛7。经验公式经过多年的发展，取得了巨大进步，但由于其

22、局限性(如预测信息量少，预测准确性差，随机性强，结果偏差大等缺点)，仍然被限制采用。而实验测试搭建实验平台要耗费大量的资金、占用空间大、研究周期长，因此实验测试分析这种方法也只能有条件采用。数值模拟计算随着计算机技术的发展得到了长足的进步，其 仅仅需要一台高性能的计算机便可以进行模拟分析，与经验公式和实验测量分析方法相比有无可比拟的优越性，并与经验公式理论和实验测量结果相互依存、相互补充，成为研究轴流泵流体流场最主耍的手段而被广泛的采用。表1-1为三种水泵性能预测方法的情况对比。表1-1为三种水泵性能预测方法的情况对比：2第一章绪论1.2国内外研宂现状从上世纪80年代以来，随着计算机性能的飞速

23、发展，使得用计算机模拟分析叶轮机械的内流场变成了可能，计算流体力学（cfd）的基本理论也得到了快速的发展，基于cfd模拟分析成为研究叶轮机械a流场流态的最主要手段。1.2.1国外研宂现状lee8等利用cfd软件对带有前置导叶的推进式轴流泵的内流场进行了定常分析，经过模拟分析结果和实验测量数据的对比得到如下结论：其流场结果和实际实验的吻合度很高，验证了 cfd流场分析的实用价值，而且其声学特性与实验结果相对应，但由于局部设计的不合理导致流体岀现分离流动，因此对泵进行了部分改型。lan9对喷水推进式水泵进行了研究，提出了其叶轮和导叶的设计方法，并对设计好的水泵进行了流场分析， 介绍了流体分析所用的

24、软件和cfd的基本理论。hu和zangeneh10对喷水推进泵进行了流场分析，采用不同的cfd软件进行了模拟计算并对其进行了分析比较，然后又对三种不同形式的推进泵进行了模拟计算，分别是只 带有叶轮的推进泵，导叶分别处于叶轮前端和后端的两种推进泵，计算的 结果表明：进口处流体的不均匀性对其扭矩的计算结果基本没有什么影响; kunikiyo, yukll等人采用三维cfd方法，对轴流泉进行了流场分析，研究了轴流泵的叶轮对其性能参数的影响，分析了导叶与轴流泵效率之间的 关系，并最终用实验证明了模拟分析的准确性；muggli12等人对一个带有导叶的轴流泵进行了单流道的流场分析，根椐流场分析的结果绘制出

25、了轴流泵的性能曲线，并对轴流泵进行了实验测试，经过对比测试的结果和模拟的结果相当吻合;miner13使用了 cfd软件flotran分析了无导叶轴流泵内的流场流动，并进行了实验研究，由于实验的轴流泵为带有导叶的轴 流泵，因此计算的结果有些误差；yu14等人对一种高雷诺数的模型泵进 行了计算，计算的结果和实验测量的结果偏差不大，仅泵出口处的速度曲线误差较大，验证了 cfd模拟计算的可靠性极其实用性。1.2.2国内研究现状目前国内对轴流泵的研宄工作取得了比较大的进展，积累了较多的实际经验，丰富了轴流泵的优化设计理论和轴流泵流场分析理论，为提高轴流泵的设计水平提供了依据。如唐宏芬15等计算了一个双极

26、轴流泵模型，对该泵模型的流场流态，压力、速度等参量进行了研宄，研宂结果和实际实验结果相对比，证明了全流道的三维紊流模拟的可行性和准确性；徐朝辉16等对一个高速泵模型进行了三维非定长的计算，得到了非稳态条件下的动叶和静叶之间的干扰流动对轴流泵性能的影响，并得到了三维紊流的流场流态；朱亮17等对轴流泵的叶片厚度进行了研宄，对几种具有 不同叶片厚度的轴流泵进行了模拟计算，计算结果表明：叶片厚度变薄，相同流量下，水泵的效率和轴功率增大，是由于叶片减薄，叶片切削水体时受到的阻力减小，通过水体的表面积变大，也就使叶片冲角变大，致使轴流泵的扬程和轴功率最终增加；朱红耕18利用cfd软件对轴流泵内壁的粗糙度进

27、行了研宄，结果表明：叶轮、导叶、轮毂轮盖等零件的粗糙度对轴流泵的性能有着极大的负面影响，导致轴流泵的效率和扬程等性能参数不同程度的下降，因此应尽量采用新工艺来降低零部件的粗糙度，应重视粗糙度对轴流泵性能影响的严重性；黄欢明19对一个轴流泵模型的内流场进行了数值计算并侧重采用新技术piv实验对轴流泵进行了全血的实验研究，实验研宄一方面证明了数值模拟的有效3江南大学硕士学位论文性和准确性，同时也对轴流泵叶轮的汽蚀性和轴流泵内流场流态、部分区域的二次流、漩涡等问题进行了详细的研究，并对叶顶间隙泄漏进行了详细的阐述。从国内外的研究现状可以看出，s前国内外关于轴流泵叶轮叶型优化的文献虽然有所进步，但关于

28、轴流泵优化设计的部分还是相对较少，主要集中在泵内流场流态的模拟分析。因此提出轴流泵叶轮、导叶叶型优化方法，以及轴流泵轮毂轮盖优化设计的方法，不仅具有理论意义，更具有实际价值。1.3本文的研究内容和方法本文研究的轴流泵模型是基于无锡锡泵制造有限公司己有模型，重新设计来满足移动泵站用泵的需要，通过对其内部流场进行了三维湍流数值模拟计算，分析流场流态并绘制其性能曲线，并在锡泵公司实验室的支持下，对轴流泵原模型进行实验研究，验证数值模拟的准确性和可行性，获得较为准确的速度、压力分布，并分析其性能损失的可能原因，找出改善的方案和最终优化设计。本文研究的主耍内容方法有：1）针对具体研宄对象的实际过流通道，

29、在建模软件中pro/e中建立出轴流泵的三维实体模型。2）应用numeca软件的前处理子模块igg对叶轮、导叶进行模型导入，使numeca可以识别，导叶叶片压力面、吸力面的线型数据和轮毂hub,轮盖shroud线型的dat数据文件，然后分别导出，并在igg屮生成叶片压 力面、吸力面的dat文件，然后导入autogrid4中生成模型的几何文件，再 将模型的几何文件导入autogrids中细致划分叶轮导叶部分的网格，并对轴流泵前动叶和后静叶的网格进行衔接等设置。3）应用numeca里的子模块fine对轴流聚全流道的流场进行求解计算。4）应用numeca里的后处理子模块cfview对求解收敛的结果进行

30、处理，得到研究模型的内部流动特性和流场流态，即轴流泵流道内水体的速度矢量分布、静压力分布等，并根据数值模拟的结果得到轴流泵的扬程、效率的水力特性曲线；分析原型泵叶轮在不同叶片数时引起的特性变化,分析原型泵5种不同叶顶间隙对轴流泵性能的影响，分析原型泵有导叶和无导叶两种情况下对轴流泵性能的影响。5）依托实验室的试验平台进行试验，通过对试验数据的分析处理并与数值计算进行比较，验证了通过数值计算模拟轴流泵内部流场的可行性，为后面轴流泵的结构分析，提出叶型优化方法提供理论支持。6）根据模拟计算的最终结果，分析可能存在的问题，应用numeca 的优化模块design 3d对叶片叶型进行优化设计，并给出叶

31、轮叶片和导叶 叶片不同的优化方法，提出叶型优化的思路和理论依据，即采用改变叶轮叶片吸压力面的型线和导叶中弧线的形状来优化轴流泵。4第二章流场数值模拟的基本方法和理论第二章流场数值模拟的基本理论和法湍流的流动十分复杂，它的结构是高度复杂的三维结构，其流动是非 稳态，杂乱无章的，没有什么规律可循的，并且还在高速的旋转21。在湍流中，描述湍流的各种物理参数也在无时无刻不发生着随机的变化，如压力、温度、速度等，它们在不同时刻和不同空间中的变化规律是随机的。如果可以把湍流看成是一种物理形态，那么它就可以看成是由众多大小不等的漩涡叠加而成的流动，并且这些漩涡的状态也是随着时间和空间的不同发生着随机的变化，

32、有的地方在这时涡的尺寸很大，但在下一刻，涡的尺寸、旋转轴等都将发生了巨大的变化，并且毫无规律可循。在湍流中，有许多较大尺寸的涡，这些涡主要是由流动的边界条件所决定，这些涡的大小可以相对很大，可以充满整个流场，其是引起低频脉动最主要的原因；湍流屮还存在着许多较小尺寸的涡旋，其大小可以非常小，小到只占整个流场的千分之一，万分之一，其主要是靠流体的粘性力附着在一起，其是 影响流动中高频脉动最主要的原因。其中大尺寸的涡和小尺寸的涡都可以相互转化，当大尺的涡在流动中破裂后形成较小尺的涡，较小尺、上的 涡又再破裂成更小尺寸的涡，并且较小尺寸的涡也可能合并生成相对较大的涡，但在这随机变化的过程中，虽然毫无规

33、律，但在一定范围内，涡的大小还是连续变化的。在这不断变化的过程中，大尺、r的涡从主流中获得能量，破裂后，将部分能量传递给较小尺寸的涡，通过这不同大小的涡之间传递着能量，最终部分能量传递给了小尺度的漩涡，这些小漩涡消失后，这部分能量也就转化成了流体的内能，造成了流体的发热现象象。由于流体受扰动的作用，并且在流体的边界条件、速度梯度等一些因素的共同制约下，上述过程不断的循环，循环的过程中，能量不断的从大尺度涡到小尺度涡传递，这个复杂的过程便形成了湍流22。2.1基本控制方程轴流泵内部的流场杂乱无章，内部是极其复杂的湍流流动，虽然看似无规律可循，但还要受到大自然的制约，满足大自然运动的规律，其还是遵

34、循质量守恒定律、动量守恒定律和能量守恒定律23。这三大定律构成了流体动力学的基本方程组navier?stokes方程组的2个部分，包括了连续性方程、动量守恒方程。2.1.1质量守恒方程质量守恒方程即是连续方程。质量守恒定律是大自然界的基本规律,任何物体都遵循着基本的质量守恒定律，而连续力*程则是描述质量守恒定律的基本方程24，其方程为：?0 ?t根据分析流体为水体得:d?0 dt带入上式得：5江南大学硕士学位论文?u?v?w?o ?x?y?z式中，？为水体的密度，u、v和w分别为速度矢量在x、y和z方向上的分量，t为时间。2.1.2动量守恒方程动量守恒定律是大自然界另一基本规律，即牛顿第二定律

35、25，动量方程即是动量守恒定律的数学表达式，方程表达式为：2?uu?l?p?ii(uu)?fiji?t?x?x?x?xjiij当泵处于稳态工作的状态时，其内部的流场处于恒定，由于其介质为不可压缩的水体，则表达式为：?2ui?l?p(uiuj)?fi?xj?xi?xi?xj其中，p表示压强，？是第二粘度，？为常数，fi表示i方向的体积力分量。2.1.3紊流时均方程对上述n?s方程进行时均推算得到如下公式：?u?u?uul?p?j'iii?f) iij?t?u?x?x?x?xjijji 由 boussinesq 对紊动粘性系数的假设，可以得到如下公式：?u?u2j'

36、i ?)k?ijij?x?x3ji根据上式换算出紊流雷诺时均方程:?u?u?u?ul?p?jiit? ?f?it?t?x?xx?x?x?ji?j?ji?上述三式屮，其屮u和p为时均值，？t表示紊动粘性系数，？为运动粘性系数。当泵体内的流场在恒定不可压缩的条件下时，紊流雷诺时均方程则可表示为：?u?u?ul?p?jii ?f?)?it?x?xx?x?xji?j?ji?2.2湍流数值计算的方法人们认识到，在大自然中，湍流在流动中占了很大比例，除了湍流， 大自然界中还有层流，相对来说，湍流是常见的，层流是稀少的26。为 了数值模拟湍流，人们提出了很多方法，如图2-1介绍了揣流数值计算方法的组成，目前

37、最常见的就是两种：直接数值模拟和非直接数值模拟。直接求解方程就是直接数值模拟方法，当前直接数值模拟的难度较大；相反的，非直接就是不直接的求解方程，而是先要对湍流进行一定的处理、简化使其能够满足现冇计算模拟的条件，非直接模拟在当前应用的比较广泛，研究的比较深入，非直接模拟根据其简化和近似处理的方法的区别，又可以将其分成三种不同的6第二章流场数值模拟的基本方法和理论方法：大涡模拟法；统计平均法；和reynolds平均法27。图2-1湍流的数值模拟方法结构图2.2.1非直接数值模拟与直接数值模拟非直接数值模拟方法中，最常用的方法是所谓的雷诺平均法(rnas)。这种方法存在两个明显缺陷29: (1)在

38、模拟计算中，要对湍流进行近似处 理，但这种方法由于采用平均运算，导致脉动运动屮的行为细节等因素完全排除在外，没有考虑其在流动中的作用，但这些细节是在高频脉动中不可或缺的内容。不仅这些大量微小的涡不能够描述，在流体中，那些较大尺寸的涡也不能很好的描述，它们在揣流的运动屮起到举足轻重的作用， 但雷诺平均法还是不能将其描述。(2)各种湍流模型的根据自己的特点都有 其自身的局限性，因为其大多是半经验公式，对经验数据有很大的依赖性，还需要大量实验数据的支持，尽管如此，其分析计算的准确性还是相对较差，不能满足预测信息的要求。直接数值模拟(direct numerical simulation,简称dns)

39、方法是指直技利 用瞬时的navierstokes方程对湍流模型进行计算30。直接模拟相对与非直接模拟有其自身的优点和缺点，其主要好处是不用对揣流模型进行假设，在理论上，其能够很好的模拟湍流。但其自身的特点是直接对湍流进行模拟计算没冇任何假设，就要其对湍流的各种信息都要加以描述，这样对其要求就更高，大到充满整个流场的尺寸很大的涡旋，小到只宥整个流场万分之一的小涡旋都要能够描述，这就要求划分的网格可以很大，也可以很小，但在实际情况中，这种要求是很难达到的，计算起来需要的计算机性能也要很高很高，当前还远远达不到。因此，当今依靠直接数值模拟湍流还不是很现实，需要很长的探索之路要走。2.2.2 reyn

40、olds 平均法为了寻找描述瞬时湍流的合理方法，科研学者们不懈努力，最终用瞬 时的n?s方程对其进行描述得到了人部分人们的认可。但其方程的非线性，却使描述三维湍流内部的全部内容变得不是那么容易，因为在湍流中，湍流引起流场的变化是耍研究的主耍内容，它是一个宏观的变化。为了求解时均化的n?s方程，人们寻求各种方法将瞬态的脉动7江南大学硕士学位论文量用一定的模型方程加以表示，最终将其转化成时均方程加以求解, 由这种思想，人们找出了 reynolds平均法，这种方法节省了大量的运算时间，使模拟成为现实。湍流的形态是复杂的，根据上述内容，将湍流瞬时 量分解为时均量和脉动量之和。ui?i?ui&#

41、39;p?p'其屮，ui表示i方向上的速度瞬时值；i表示i方向上的速度平均值； ui'表示速度i方向上的速度脉动值；p表示压力瞬时值；表示压力平 均值；p'表示压力脉动值。连续方程：动量方程：di?0 dxi?i?ui'u'j?x?j?t?ji?i?t?xj?xi?xj上式即为雷诺方程(reynoldsequation)。根据：'?ij?i'j'''?ij为雷诺应力(reynolds stress),其共有9项。把雷诺应 力表示为一个对

42、称的张量?ij:''?ul'ul'?ul'u2?ul'u3?''&# 39;''''?ij?u2ul?u2u2?u2u3?''''''?u3ul?u3u2?u3u3?其屮，？?ui'ui'为正应力项，？?ui'u'j为切应力项。2.3 sp

43、alart-allmara (sa)湍流模型1992年spalart和allmaras提出了 sa模型。该模型是介于代数方程和两方程模型中间的一种湍流模型，其是一种一方程模型。通常情况下，它被认为是连接代数零方程和两方程模型的桥梁33。这种模型具有计算精度高、cpu占用率低、内存消耗低等优点，故近年来得到了广泛的应用。它是根据实际经验和量纲分析出发而不是仅仅对k-e两方程模型直接简化得到的，它是一种先从简单的流动逐渐充分发展，最后适用于带有层流流动的方程模型，总的来说它是一种适用于固壁湍流模型的一方程模型。在模型'于是有：中，通常情况下，应变量选用的是与涡粘性?相关的量?tf，f?t?

44、vlvl?, ? 33?cvl?3上式中下标?表示“粘性”。cvl取值为6.1。？表示粘性系数。上述公式满足从壁面?ky到对数律这一区间内存在着这一条件。可得：8第二章流场数值模拟的基本方法和理论-d?-?l?-?c(?-)2?cf?cblf?-/d?2?f?u2?cbl?l?ft2?-s?blwl wt2tl?dt?k2?其中雷诺应力项：f，f?t?vlvl?, ? 33?cvl?3?ui?uj?/2 ?uiuj?2?tsij, sij?x?j?xi?？、s?u?/kds?s?22fv2, fv2?l?l?fvlkd666fw?g(l?cw, g?r?cw2(r6?r)，3)/(g?cw3)

45、？r?22, ft2?ct3exp(?ct4?2) skd?l/6wt2222d?gd 转捩函数:ftl?ctlgtexp(?ct2tt) 2?u?上式中d表示壁诎到流场中场点的距离;dt表示转捩带到场点的距离;wt是转捩处的旋度；？u表示转捩处到流场中某一场点的速度增量； gt?min(0.1,?u/wt?x)； ?x表示转捩点处壁面第一层网格的大小。各个常数的值取：cbl=0.1355,?=2/3,cb2=0.622,k=0.41,cwl?cbl/k2?(l?cb2)/?,cw2?0.3,cw3=2,cvl=7.1,ctl=l/ct2=2,ct3=l.l,ct4=2o2.4掠制方程组的求解方法在进行湍流模型的数值计算时，虽然在求解域内建立了数学模型，但要直接求解这些微分方程组还是有一定的困难，这就必须转化成微分的思想，将计算区域的模型分割成一个个微小的单元。在进行三维湍流的模拟计算问题中，首先应该对计算域进行网格划分，常见的网格单元有：四面体、六面体、棱锥体和楔形体等35。根据节点的数量，每种单元又可以分成几种类型，例如，六面体可以分为8节点、20节点、27节点等不