FLUENT在超音速蒸汽喷射真空泵中的应用

FLUENT在超音速蒸汽喷射真空泵中的应用徐海涛1,桑芝富1,卓震2(1.南京工业大学机械与动力工程学院，江苏南京，210009;2.江苏石油化工学院机械系，江苏常州，213016)摘要：采用大型计算流体动力学软件FLUENT对蒸汽喷射真空泵的超音速混合过程进行数值模拟。计算并分析了第二喉管与工作蒸汽喷嘴喉管面积比、喷嘴出口截面与混合段入口截面间的距离及混合段的锥度等几何参数对真空泵操作性能的影响。数值结果表明，几何参数的改变极大地影响着波系结构，在一定的设计工况下，总存在一个最佳的面积比及一个最优的相对位置对应于最大的喷射系数，其在物理上的表现形式为通过工作蒸汽喷嘴所产生的激波系刚好能够通过第二喉管。混合段的锥度在一定范围内对真空泵的性能无显著影响，等压混合理论较等面积混合理论具更优的操作性能。关键词：蒸汽喷射真空泵；超音速；喷射器；激波；几何参数；操作性能CFDsimulationofthesupersonicsteam-jetvacuumpumpXuHai-tao1,SangZhi-fu1,ZhuoZhen2(1.SchoolofMechanicalEngineering,NanjingUniversityoftechnology,Nanjing210009;2.JiangsuInstituteofPetrochemicalTechnology,Changzhou213016)Abstract:Thesupersonicsteam-jetejectorasafunctionofvacuumpumpwassimulatedusingFLUENT.Theperformanceofsteam-jetvacuumpumpwasinvestigatedbychangingthesecondarytotheprimarythroatarearatio,therelativepositionofthesteamnozzleandthetaperofthemixingsection.Thesenumericalresultsclearlyindicatethatwhengeometricparametersvary,thestructureofshockalsochanges.Thereareoptimumthroatarearatioandidealsteamnozzlepositionforagivenoperatingcondition,whichcorrespondstothestatethatshockwavesgeneratingatthesteamnozzlemovesdownstreamandjustpassesthesecondthroat.Thetaperofthemixingsectionhaslittleeffectontheperformanceoftheejectoroverarangeandtheconstant-pressuremixingtheoryisbetterthanconstant-areaone.Keywords:steam-jetvacuumpump;supersonic;ejector;shockwave;geometricparameters;performance蒸汽喷射真空泵具有结构简单可靠、运转费用低廉、操作维修方便以及对被抽介质无严格要求等优点，因此被广泛应用于各个工业领域，尤其是在某些特殊生产工艺上，如真空蒸馏、真空脱气、真空结晶等。提高这种设备效率的途径在于使喷射器无论是压缩比(PC/PH)还是喷射系数(GH/GP)都达到最优。这就有必要重新审视传统的超音速蒸汽喷射器的设计方法。作为一种真空获得设备，喷射器诞生于19世纪40年代，最初的设计方法都是以实验结果为依据的[1]。图1为典型的喷射器结构简图。工作蒸汽在拉伐尔喷嘴中加速形成超音速射流，而引射流体则由于与工作蒸汽间的剪切作用被卷吸至混合室，而后逐渐形成单一均匀的混合流体，经过一扩散段减速压缩到一定的背压。在混合过程中激波系不断和边界层相互作用形成了极为复杂的流动结构，再加上粘性干扰、分离涡、真实气体效应等物理现象使得喷射器的性能很难用一般的气体动力学理论解释。现有的超音速喷射器的设计理论大都沿用一维空气动力学的分析方法，其主要困难在于对混合过程的动量守恒方程提出一种合理的解析解：Keenan，Neumann[2-3]，Elrod[4]和Fabri[5]等人先后提出了两种比较可行的计算方法，即定常面积混合理论和定压混合理论。前者认为工作流体和引射流体的混合过程是在截面积不变的情况下完成的(A=constant)，而后者则认为两股流体的混合近似为一种等压过程(dP=0)。无论是一维定常面积混合理论还是一维定压混合理论，其出发点都是一样的，即在等压混合或定常面积混合理论的基础上，假设工作流体和引射流体具有相同的分子量和比热容，工作流体和引射流体以及混合后的流体在任意截面上具有均匀的物性分布，喷嘴和扩散段内都是等熵过程，不计壁面摩擦，也不考虑热量损失，在理想气体的基础上，运用质量、动量及能量守恒方程计算工作过程，推导出了最大喷射系数或最大压缩压力的喷射泵结构参数的计算公式。尽管后来有不少学者[6-11]对其进行修正使其更接近于现实，然而现有的一维设计理论尚不能预测混合段所需的最短距离，也不能给出工作蒸汽喷嘴和混合室的最佳相对位置，以及能达到最佳工作性能的喉管面积比等几何参数。除了这些因素，还有诸如工作蒸汽喷嘴出口马赫数、雷诺数、压缩比、膨胀比、比热比等操作参数均不同程度地影响着喷射泵的性能，探讨这些参数对蒸汽喷射泵性能的影响，对完善现有的一维设计理论具有重要的意义。作为一种研究流体流动的新方法，计算流体动力学(ComputationalFluidDynamics,CFD)在各个工业领域得到越来越广泛的应用，目前，CFD方法已开始成为工程装置优化及放大的定量设计工具。它不仅能给出设备内大致的流场结构，而且能提供一些实验不能或者很难测定的信息。本文采用大型计算流体动力学软件FLUENT，用二阶精度的有限体积法(FVM)离散控制方程，标准k—ε二方程湍流模型，近壁面处使用壁面函数修正的方法对整个喷射泵的内部流场进行数值模拟。探讨蒸汽喷射真空泵主要结构参数对其性能的影响，分析在不同喉管面积比及各种喷嘴与混合室入口截面距离下喷射泵内的激波运动规律，揭示在一定操作条件下最佳喷嘴位置及最优喉管面积比的主要物理特征，为蒸汽喷射真空泵的优化设计建立理论基础。1蒸汽喷射泵的数值模拟算法1.1数值算法[12]采用二阶精度的有限体积法(FVM)对控制方程进行空间离散，时间离散则采用多重Runge-Kutta显式格式迭代，并采用多重网格方法(MultigridMethod)加速迭代收敛，自适应网格方法(GridAdaptive)捕捉激波，以期获得与网格无关的解。1.2计算域及网格鉴于在混合段入口截面处引射流体的速度与工作蒸汽速度相比甚小，故可将引射流体的侧向入口简化成轴向环形入口，因而可简化成二维轴对称模型计算，数值结果表明，这种简化的二维模型和三维模型并无明显差异(如图2所示)，其得出的喷射系数相对误差在1%之内。图3为本文所采用的二维计算模型网格(为便于显示，只取实际计算网格密度的20%)及边界层网格的局部放大图。网格采用分块结构化网格，为捕捉激波和边界层，在喷嘴出口及工作流体与引射流体混合面上适当加密。1.3流体物性及边界条件工作蒸汽及引射流体入口均采用压力入口边界，给定滞止压力、滞止温度及适当的湍流条件；混合流体出口采用压力出口边界，给定静压及适当的回流条件；固体壁面采用无滑移、无渗流、绝热边界；中心边界为对称轴，各变量在对称轴法线方向的梯度及垂直于对称轴的速度为零。引射流体为蒸汽，其密度按理想气体计算，粘度则由温度的指数函数形式给出。2模拟结果及讨论2.1混合室流场结构及引射机理探讨气流在蒸汽喷嘴中做正常的降压膨胀，在喷嘴出口形成了比引射流体还低的压力，所以超音速气流流出喷口后，将受到高背压的压缩而形成口外激波系，如图4所示，轴线上方为马赫数等值线，下方为流线图。由于在本算例中引射流体的压力(8932Pa)比喷嘴出口的压力(轴线处为6365Pa)高得不是很多，喷嘴出口处的斜激波在中心处能正常相交，并再次形成压缩波，气流经两次压缩后，马赫数有所降低，而压力则上升到比周围介质还高的水平，从而气流在压缩波与边界层交界的周沿受扰动后产生膨胀波，而后膨胀波相交再次形成膨胀波，气流经两次膨胀后，速度增加而压力下降，通常下降得比周围压力还低。结果在它继续流动过程中必然受到周围介质的压迫，在边界层上反射再次形成压缩波。若无气流的流程损失，它就压缩波—膨胀波—压缩波这样一个波节一个波节的流下去。但由于能量损失以及壁面边界的作用，射流逐渐衰减，经第二喉管时又重新减压增速，形成了极为复杂的波系结构，这将在下文分析。由于压缩波和膨胀波的不断相交、反射，气流在经过这些波线的时候不断折转一定的角度，形成了图4所示的流线结构。从图4可以看出，超音速射流流出喷嘴后，在射流的外缘形成边界层。在这一边界层中，速度由最外沿的低亚音速一直变到与射流核心相接触处的超音速，射流核心轴线上的静压随波系结构呈准周期性变化。图5给出了混合段压力分布侧型图，可见在射流横断面上，压力分布也是极为不均一的，射流里有横断方向的压力梯度。在压缩波之后的各个截面上压力梯度指向射流边界，也即中心的压力要比边界高；而在膨胀波之后的各个截面上则指向轴心。蒸汽喷射器内流体流动的这些特点，使得射流与周围介质不断的进行质量、动量、能量的交换。卷吸进入的流体取得动量而与原来的射流一起向前运动，结果射流断面不断扩大，而流速则不断降低，流量却沿程增加。由此可见，与一般的流体输送机械采用叶片的正应力提高流体压力所不同的是，蒸汽喷射器采用的是两股流体间的剪切力带动引射流体向前运动，激波系和边界层之间的交互作用起着至关重要的作用，其结构参数极大地影响着喷射器的操作性能。2.2喉管面积比的影响蒸汽喷射式热泵工作能力用喷射系数u来表征，它表示在一定工况下，单位质量工作流体通过泵所能抽吸的引射流体(被抽气体)的量，在数值上等于引射流体的质量流量与工作流体的质量流量之比，即u=GH/GP。维持操作参数及喷嘴的几何尺寸不变，仅改变第二喉管的直径所得到的喷射系数u与喉管面积比ψ之间的关系如图6所示。在所有计算工况下，混合段及扩散段的锥度均保持不变。可见，存在一最佳的喉管面积比ψopt对应于最大喷射系数，在ψopt两侧u(ψ)均呈线性率变化，且下降的速度明显快于上升的速度。图7，8分别给出了不同面积比下马赫数及静压的分布情况，图7中的(a)-(e)分别与图6中的a-e点及图8中的5条曲线相对应。可见，在所有计算工况下，混合段及第二喉管段轴线上的静压始终处于准周期性振荡，并逐渐衰减。而当ψ较小时，在扩散段中再次出现激波系，其第一波阵面随ψ的增大向第二喉管处移动，导致轴线上的静压在出口附近仍处于急剧的振荡之中，这是因为在这种情况下第二喉管出口处的马赫数仍处于一个较高的水平，气流尚不能通过正常的等熵膨胀来适应混合流体的出口压力。这点同样反映在壁面压力的分布上，从图8可以看出，ψψopt时，第二喉管段壁面的压力不断波动，可见此时第二喉管对于激波系的扩展起着相当重要的作用；而在扩散段射流分离点以后，壁面压力渐趋平缓并逐渐上升，这是由于边界层沿扩散管逐渐发展从而激波系对壁面压力的影响逐渐减弱的缘故；锥形混合段壁面的压力则随ψ的增大而降低，意味着混合过程中动量损失逐渐减少，而喷射系数则随之增大，喷射器的工作能力日趋增强。增大ψ，扩散段中激波系的第一波阵面不断向第二喉管处移动，在某一ψ=ψopt下，气流刚好在第二喉管出口截面处不产生激波，喷射系数同时达到最大；ψψopt时，锥形混合段壁面的压力随ψ的增大而增大，意味着混合过程中动量损失亦逐渐增加，喷射系数则随之减小，热泵的工作能力日益减弱。而此时第二喉管段轴线及壁面压力都保持在一个稳定的