首页
登录
注册
-
回收和利用工业余热、地热能和太阳能等低品位热能，是实现全球能源的可持续发展的重要途径。喷射器作为一种能量回收装置，其结构参数是影响其性能的重要因素。经典的喷射器模型[1-3]能较准确地预测喷射器的径向尺寸与其性能的关系，而喷射器的轴向结构，如喷嘴距（NXP）和混合室轴向长度对喷射器性能影响还需更多研究。于文艳等[4]数值研究混合室的轴向结构发现，当混合室的圆锥段和圆柱段长度分别为混合室圆柱段直径的2倍和6倍时，蒸汽喷射器的喷射系数达到最大。Dong 和Wang等[5-6]的数值结果均显示当混合室收敛段长度小于或等于某一长度时，混合室喉部长度对蒸汽喷射器的喷射系数影响不大。Wu等[7]的数值结果显示，混合室收敛段长度存在最佳范围，该最佳范围取决于工作压力。
喷嘴距（NXP）是喷射器的重要几何结构。迄今为止，喷嘴距主要通过实验或CFD方法确定。Yapc等[8]通过实验指出，NXP的精确值是无法准确地定义的，因为它会随着操作条件的变化而变化。此外，最佳喷嘴距随其他几何结构的变化而不同[8-10]。工程科学数据单元（ESDU）[11]建议，最佳喷嘴距与混合室喉部直径的比率，即最佳无量纲喷嘴距应为0.5~1。其中，NXP是指从喷嘴出口到混合室喉部入口的距离。
Poirier[12]的研究结果显示，最佳NXP值受喷嘴几何形状和临界喷射系数影响，首次证明了最佳喷嘴距是临界喷射系数的函数，而不是二次压力的函数，进出口压力的影响可忽略不计。Yan等[13]的实验结果显示，喷射系数随NXP的增加先增加后减小，最佳DNXP为1.88，NXP对临界背压影响不大，得出相同结论的还有文献[14]。Ruangtrakoon等[15]的研究结果指出，随着喷嘴距减小，喷射器的临界压力增加。Rusly等[16]研究结果显示最佳喷嘴距为1.5, 且在±20%的范围内变化不会显著影响喷射器性能。Chen等[17]实验结果显示，喷射器的临界喷射系数随NXP的增加而增加，然后保持最大值不变, 并提出了一个一维模型来确定不同NXP下的临界喷射系数。Valle[18]实验研究三种混合室喷射器在不同喷嘴距的临界条件，结果显示锥形混合室喷射器的喷射系数在喷嘴距−15~−1 mm范围保持恒定，超过该范围后急剧下降。Kumar等[19]研究变面积混合喷射器的喷嘴距影响，发现喷嘴远离或进入混合段时，喷射系数降低，最佳喷嘴距为混合室收敛段入口处。Shi等[20]的实验结果显示，最佳无量纲喷嘴距为0.5，此时喷射系数和升压比均达到最大值。Chen等[21]通过引入环形混合层理论，建立了一种新的最佳喷嘴距预测模型，并对圆柱形混合室的喷射器进行了验证。Cagri等[22]以R134a为工质，对超音速喷射器进行了数值研究，结果显示最优喷嘴距为2.04。Nguyen等[23]以R600a为工质，实验研究了不同操作条件下，喷嘴距对太阳能驱动可变几何喷射器系统性能影响。结果表明，发生压力越高，NXP对COP的影响越大，最佳NXP值不随操作条件而改变。Wang等[24]通过设计了一种自调谐面积比和NXP的喷射器，研究发现最佳NXP随工作压力的增加而减小。Han等[25]数值研究显示最佳喷嘴距为4.0，增加NXP会加剧边界层分离。Ramesh等[26]使用CFD仿真和实验方法，研究了蒸汽喷射器吸入室角度和喷嘴距影响，结果显示喷嘴远离混合室，喷射系数将增加到最大值，然后降低，最佳喷嘴距值随着吸入室角度的增加而减小。
综上所述，针对喷嘴距的研究，主要集中在单一结构优化、操作条件影响、喷嘴距与喷射器的径向结构共同优化这几个方面，对喷嘴距和混合室喉部长度这两个关键轴向长度对喷射器性能的综合影响，以及不可逆性研究尚不充分。此外，喷射器性能常采用喷射系数和临界压力来衡量，而在这两个性能指标的变化规律不一致的情况下，通常难以评估喷射器的实际性能。本文采用实验与数值研究相结合的方法，对混合室喉部长度和喷嘴距对喷射器性能的综合影响进行充分研究，并引入热力学完善度对喷射器整体性能进行评估，建立喷射器熵产模型，确定了发生不可逆性的区域，分析混合室喉部长度和喷嘴距对喷射器熵产影响的内在机制，为实现喷射器结构及综合性能的多目标优化，以减小不可逆损失提供参考。
本研究使用无量纲喷嘴距DNXP，其定义为喷嘴出口到混合室喉部入口的距离，如图1所示。此外，还对混合室喉部长度进行无量纲化，其定义为混合室喉部长度与混合室喉部直径的比值。
式中，NXP、Lm,t、Dm.t为喷嘴距、混合室喉部的长度和直径，mm。
-
图2为可调式喷射器实验系统示意图[27]。实验系统由发生器、蒸发器、可调式喷射器、过热器、真空泵、冷却循环系统、测量和控制系统组成。
实验装置如图3所示。其中，发生器和蒸发器都配置有浸没式电加热器组，通过控制加热功率，从而控制实验压力和温度。喷嘴座和喷射器主体通过螺纹连接，从而实现喷嘴距可调。喷射器几何结构见图1。在理论计算过程中，由于滑移或粘性效应，有效理论流通截面积被认为是实际截面积的0.88倍[21]。冷却循环系统包括冷凝器、水泵和冷却塔。蒸发器、发生器和冷凝器均设石英玻璃管液位计，以便读取液位。过热器采用电加热带，并包裹在工作蒸汽管道外，从而保证工作蒸汽具有一定的过热度。真空泵用于调节出口冷凝压力。测量和控制系统包括测量仪器和数据采集控制单元。主要的测量仪器是压力和温度传感器。数据采集控制单元由数据采集模块、A/D转换器和自动控制系统组成，通过自动调节发生器和蒸发器中电加热器的功率，可以准确地保持系统的预设条件[28]。
-
实验以水蒸气为工质，工况为发生器：pp = 160.65 kPa，Tp = 396.57 K，蒸发器 ：ps = 2.81 kPa，Ts = 296.15 K。喷嘴距和混合室喉部长度参数通过无量纲命名，考虑了12个无量纲喷嘴距DNXP （4.07、4.26、4.54、5.72、5.0、5.19、5.37、5.65、5.83、6.11、6.3）。类似地，考虑了5个无量纲混合室喉部长度DMTL （0.5, 2, 4, 6, 8）。喷射器的固定结构尺寸如表1所示。实验的目的是测试不同无量纲混合室喉部长度下，无量纲喷嘴距发生变化时喷射器的性能。
需测试出喷射器在不同喷嘴距的临界值，具体实验步骤为：
（1）选择并安装目标锥形混合室喷射器，调整无量纲喷嘴距DNXP到目标值。开始前，确保发生器和蒸发器中有足够的水，且全系统具有足够真空度。
（2）确保电脑控制系统与仪器连通，设置发生器目标压力160.65 kPa，温度396.57 K，并启动加热。当发生器达到目标压力和温度后，启动外部循环冷却系统。
（3）依次打开冷凝器的混合蒸汽进口阀，发生器、蒸发器的出口蒸汽阀，设置蒸发器的目标压力2.81 kPa，启动电加热带。
（4）启动真空泵，调节冷凝器内的不凝气体量，保证冷凝压力处于目标值。当蒸发器参数稳定在目标值后，记录发生器和蒸发器的液位。
（5）计算喷射系数并判断是否为临界值。如果是，该DNXP的测试完成，调节下个DNXP；如果不是，继续调节冷凝压力，直至临界值。
-
利用Fluent软件进行数值模拟，建立喷射器二维轴对称模型，采用结构化网格，在壁面及流速变化较大的区域进行网格加密。选择基于密度的稳态求解器和二阶迎风格式，离散化系统采用高斯-塞德尔方法。工质为水蒸气，其密度采用真实气体模型，其他属性均为温度的多项式函数。喷射器壁采用无滑移和绝热壁型，边界条件选择压力入口、压力出口。选择SST k-ω 湍流模型，因为它对全局和局部流动现象的预测精度更高[29-30]。工况与实验相同，工作蒸汽过热度为10 K。计算时，当进出口的质量流量稳定，且计算的每项残差小于10−3时，则计算值已收敛。
为测试网格独立性，当DMTL= 4， DNXP = 6.11时，分别设置网格数为105622、142530、201175（模型的一半），不同网格数的临界喷射系数误差如表2所示。当网格数为142530和201175时，喷射器的临界喷射系数误差几乎相等。不同网格数的马赫数如图4所示，可以看到，三个网格数的喷射器的马赫数波动基本相同，考虑到误差及计算成本，最终选定142530左右的网格数进行模拟。
-
图5为DMTL = 4， DNXP = 6.11时，出口冷凝压力对喷射系数的影响。由图可知，模拟和实验结果的变化趋势相同，存在一临界压力值，当出口压力小于临界压力值时，喷射系数保持不变；超过临界值后喷射系数开始下降；继续增加压力，喷射器将不能正常工作。固定冷凝压力pc = 5 kPa, 表3为各DMTL、DNXP组合结构的模拟和实验结果对比。
由表可知，在测试的所有DMTL和DNXP中，通过实验和模拟得到的临界喷射系数值接近，但模拟值更高。模拟值与实验值相比，临界喷射系数的最大误差在±13.23%以内，而相对偏差绝对值的平均值为9.58%；临界冷凝压力的最大误差在±19.89%以内，而相对偏差绝对值的平均值为11.77%。对差异的存在主要有以下几个原因：首先，SST湍流模型具有过度预测特性；其次，实际环境中无法避免喷射器与周围环境之间的热传导，且光滑无任何摩擦的壁面在目前是难以达到的；最后，可能因为喷射器尺寸小，制造加工不够精确。
-
喷射器性能采用临界喷射系数(ucr)、临界冷凝压力（pcr）以及热力学完善度（η）进行衡量。喷射系数定义为引射蒸汽与工作蒸汽的质量流量之间的比值。
当喷射系数达最大值时，称为临界喷射系数ucr，相应的出口压力为临界压力pcr，pcr越大，说明喷射器的工作能力范围更宽。
热力学完善度反映了设备实际状况与理想状态之间的差距以及不可逆性。Valle等[18]将热力学完善度定义为喷射器的实验喷射系数和理论喷射系数的比值，并用于评估三种混合室喷射器的整体性能，确定了最佳混合室结构。Chen等[31]通过分析不同压缩比和膨胀比时双级喷射器的热力学完善度，为一级和二级引射器优化选型提供了指导。Sierra-Pallares等[32]研究三种混合室喷射器熵产，将理论计算出的熵产与热力学完善度联系起来，结果显示这两者成反比关系，确定了喷射器内与熵产相关的四种不可逆机制。
在本研究中，临界喷射系数和临界冷凝压力由喷嘴距和混合室喉部长度共同决定，其数值无法反映出喷射器的不可逆性；且随DNXP的增加，两个性能参数的变化规律不一致，导致喷射器性能在不同的结构条件下也没有可比性。因此，引入热力学完善度来讨论喷射器偏离理想工作状态的程度。假定喷射器是可逆的，熵产为零，喷射器在理想状态下有：
式中：ui为理想喷射系数；η为热力学完善度；ucr为临界喷射系数；sp、ss、sc为工作蒸汽、引射蒸汽和混合蒸汽的比熵，kJ/(kg·K)； mp、ms、mc为工作蒸汽、引射蒸汽和混合蒸汽的流量，kg/s。
考虑喷射器的有效性，预计喷射器熵产大于零。总熵产率分为热耗散和粘性耗散两部分[33-34]：
由于绝热假设，传热部分忽略不计，喷射器的总熵产率为：
式中，T为温度，K；τ为粘性应力张量，N/m2；ν为速度矢量，m/s。而粘滞耗散分为平均流场耗散和速度波动引起的耗散[35]：
式中，平均流场熵产Smean为：
由速度湍流波动引起的耗散熵产Sflu：
本文采用SST k-ω湍流模型，由文献[36]可知：
得到速度波动引起的耗散熵产：
式中，ρ为密度，kg/m3；ε为湍流耗散率，kg/( m·s4)；ω为比耗散率，kg/ (m3·s2)；Cμ为经验常数,取0.09；k为湍流动能，kg/(m·s3)；μ为动力粘度，N·S/m2；μt/μ为湍流粘度比，取10%[22, 36-37]。
-
在无量纲混合室喉部长度DMTL= 0.5, 2, 4, 6, 8情况下，喷射器的临界喷射系数ucr随DNXP的变化如图6所示，可以看到，在DNXP = 4.07时，临界喷射系数处于最低值，随着DNXP的增加，临界喷射系数逐渐增加，并均在最佳无量纲喷嘴距5.83时达最大值0.45，继续增加喷嘴距，临界喷射系数下降。以DMTL= 0.5为例，在DNXP由4.07增加到5.83过程中，临界喷射系数由0.393增加到0.45，最大偏差达14.5%，可看出DNXP对临界喷射系数的影响是显著的。
以上分析得出，对于DMTL= 0.5, 2, 4, 6, 8 喷射器，均存在一最优喷嘴距使得临界喷射系数最大，且最佳DNXP值相同，最佳无量纲喷嘴距DNXP为5.83时，喷射器的引射性能最佳。
通过模拟喷射器的内部流场来分析喷嘴距影响。模拟工况为pp = 160.65 kPa，Tp = 396.57 K，ps = 2.81 kPa，Ts = 296.15 K ，pc = 5 kPa。以DMTL = 4为例，图7、图8分别为不同DNXP下的速度云图和质量流量曲线。由图可知，开始时DNXP较小，引射蒸汽和工作蒸汽没有充分的时间混合就进入混合室喉部，工作蒸汽继续膨胀，形成较长的激波链，占据混合室喉部大部分通道，同时由于喷嘴的阻挡，减小了引射蒸汽的入口流通面积，引射蒸汽流量较小，如图8所示；随着DNXP的增加，工作蒸汽有足够的时间和空间来引射更多的蒸汽，引射流量逐渐增加并在DNXP = 5.83时达到最大值；继续增加DNXP，由于喷嘴出口到混合室喉部距离较长，喷嘴出口的工作蒸汽射流在混合室收敛段内迅速扩散，引射流体获得动量减小，激波链缩短并集中在混合室收敛段内，引射蒸汽流量减少。
从上述分析可知，不同DMTL的喷射器，喷嘴距改变均使得临界喷射系数出现最大值。然而实际工艺对混合蒸汽的压力有明确需求，如果提高临界喷射系数时带来了临界冷凝压力显著下降，那这种意义上的效率提升将毫无意义。因此，需要进一步对喷射器在变喷嘴距下的临界冷凝压力进行分析。
-
在无量纲混合室喉部长度DMTL= 0.5, 2, 4, 6, 8情况下，DNXP对喷射器的临界冷凝压力影响如图9所示。就整体而言，临界冷凝压力均随DNXP的增加逐渐下降，且下降的局部变化因DMTL有所不同。DMTL= 4时，在DNXP由4.07增加到6.3过程中，pcr由9.5 kPa下降到8.8 kPa，下降了7.36%，下降过程呈现波动性；而DMTL= 8时，临界冷凝压力随DNXP的增加几乎呈稳定的线型下降，pcr由DNXP = 4.07的10.3 kPa下降到DNXP = 6.3的7.6 kPa，降低了26.2%；DMTL= 0.5喷射器的临界冷凝压力下降了17.4%。
实际上，临界冷凝压力的下降与混合流体的总动能有关。通过对DMTL= 4的喷射器进行数值模拟发现，随着DNXP的增加，混合室收敛段内的涡流逐渐增大，如图10所示，原因是DNXP的增加，加剧了流体与混合室壁的碰撞，导致流体边界层分离愈加严重，混合流体总动能减小，临界冷凝压力下降。由此可知，减小DNXP能够增加流体的总动能，使喷射器获得更高的临界冷凝压力，但也要避免临界喷射系数的进一步降低（见图6）。
为解释不同DMTL喷射器的临界冷凝压力下降变化的差异，图11为DMTL = 0.5, 4, 8 喷射器在不同DNXP的速度云图。由图可知，DMTL= 4的激波链长度最长，较强的射流冲量提高了扩压室的压力恢复能力，使得流体动能衰减速度随着DNXP的增加有所减缓，但可能由于激波链在混合室喉部缓慢移动时，流体之间的相互作用力使得临界冷凝压力下降有所波动；当DMTL= 0.5时，过短的喉部长度极易导致激波链由扩压室直接向混合室上游移动，临界冷凝压力急剧降低；当DMTL= 8时，激波强度减弱（相比于DMTL= 4），随着DNXP的增加，受反向压力梯度影响，激波链向混合室移动速度加快，临界冷凝压力下降速度增大。
喷嘴距从4.07增加到6.3的过程中，在所有DMTL中，DMTL= 8的喷射器在DNXP = 4.07~5范围有最高的临界喷射系数，但最高临界冷凝压力在DNXP = 4.07~4.26范围内；DMTL= 0.5的喷射器在DNXP = 5.19~6.3范围有最高的临界喷射系数，但DNXP = 4.07~6.11范围的临界冷凝压力最低，这与临界喷射系数所不一样的。
以上分析可知，一定混合室喉部长度的喷射器，其最高的临界喷射系数和临界冷凝压力对应的喷嘴距可能并不一致，不同DMTL的喷射器的性能更是无法直观比较。因此，有必要对喷射器的不可逆性进行分析，在综合考虑喷射器临界喷射系数和临界冷凝压力的基础上，引入热力学完善度来对喷射器实际性能进行评判。
-
不同DMTL喷射器的热力学完善度随喷嘴距变化情况如图12所示。由图可知。在所有的DMTL中，当DMTL= 4，DNXP = 5.83时，热力学完善度达到最大值0.21，此时增加或减小DNXP都会导致热力学完善度减小。
为确定喷射器不可逆损失增加的原因及来源通过在CFD仿真后处理中建立熵产模型，得到喷射器内部熵产率分布。图13为DMTL= 4时，不同DNXP下喷射器的熵产分布情况。
由图可知，随着DNXP的增加，混合室的熵产逐渐增加，扩压室内的熵产则相反，原因是DNXP的增加，导致混合室内的流体边界层分离愈加严重，流体粘滞损失增加。相反，减小DNXP，流体边界层分离减弱，射流总动量增加，第二冲击波逐渐向亚音速扩压室移动（见图7），提高了喷射器的临界冷凝压力，这在Ruangtrakoon等[15]的研究中得到证实。然而，他们却没有对第二冲击波带来的不可逆损失做更多研究。对此，本研究将做出重要补充。
第二冲击波向扩压室移动，导致扩压室内流体的粘滞损失增大，熵产增加。因此可知，由喷嘴距引起的喷射器熵产主要来源于，两股入口流在混合室收敛段的混合和第二冲击波向扩压室出口的移动。
图14为DMTL= 0.5, 4, 8 喷射器在不同DNXP的熵产分布。由图可知，DMTL= 0.5喷射器的熵产变化最为急剧，如DNXP = 4.07时，熵产主要分布在整个扩压室，混合室的熵产较少，增加DNXP至6.3，扩压室熵产急剧减小，混合室的熵产成为不可逆损失的主要来源，原因是DNXP较小时，工作蒸汽和引射蒸汽没有足够的长度混合均匀就进入扩压室，第二冲击波移动到扩压室中部（见图11），增加扩压室内的粘性耗散；随着DNXP的增加，受反向压力梯度影响，激波链和第二冲击波迅速向混合室移动，加剧流体边界层分离，混合室熵产急剧增加。当 DMTL= 4时，随着DNXP的增加，足够的喉部长度促使流体混合均匀，提高了扩压室的压力恢复能力，激波链及第二冲击波逐渐向上游移动，有效减缓了混合室熵产的增加。当DMTL= 8时，其激波链长度缩短（见图11），DNXP较小时，较长的喉部长度避免了第二冲击波直接冲向扩压室内部，减小了扩压室的熵产，此时热力学完善度最大。随着DNXP增加，多余的长度增加了流动阻力，导致没有足够的射流冲量，加快了第二冲击波向混合室移动，混合室熵产增加，扩压室的熵产最小。
由以上分析可知，改变喷嘴距，喷射器的临界喷射系数或临界冷凝压力达到最大时，其热力学完善度不一定最大。喷嘴距改变时，喷射器熵产主要来源于两股入口流在混合室收敛段的混合和第二冲击波向扩压室出口的移动，而混合室喉部长度直接影响熵产的变化速率。
因此，对喷射器进行结构设计和优化时，在多种性能指标变化不一致的情况下，考虑减小设备的不可逆性，引入新的性能参数进行全面评估，可通过熵分析或CFD仿真进行结构的最佳组合，这对于实现喷射器结构和综合性能的多目标优化具有重要意义。
-
本文以水蒸气为工质，通过实验与数值分析相结合方法，研究不同无量纲混合室喉部长度下，无量纲喷嘴距对喷射器性能的影响，研究结论如下:
（1）在不同的无量纲混合室喉部长度下，喷射器的临界喷射系数随着DNXP的增加先增加后减小，在达到最佳无量纲喷嘴距（DNXP = 5.83）时，临界喷射系数达最大值0.45。合理的喷嘴距应使得工作蒸汽有足够的时间和空间引射更多引射蒸汽，但应避免激波链出现在混合室收缩段。
（2）在不同的无量纲混合室喉部长度下，DNXP的增加使得喷射器的临界冷凝压力呈下降趋势，当无量纲混合室喉部长度DMTL为4时，其 临界冷凝压力下降了7.36%，下降速度最慢。
（3）在无量纲混合室喉部长度DMTL = 0.5, 2, 4, 6, 8 的喷射器中，当DMTL= 4时，喷射器存在一最佳喷嘴距，使得临界喷射系数和热力学完善度同时达到最大值。由喷嘴距改变而引起的喷射器熵产主要来源于，两股入口流在混合室收敛段的混合和第二冲击波向扩压室出口的移动，提高扩压室的压力恢复能力能有效避免喷射器熵产的急剧增加。
（4）对于既定的工作蒸汽、引射蒸汽状态参数及喷射器径向尺寸，合理的喷嘴距和混合室喉部长度组合能使得喷射器的热力学完善度最高。当DMTL= 4，DNXP = 5.83时，热力学完善度达到最大值0.21。对喷射器进行结构设计和优化时，可通过CFD仿真来实现结构的最佳组合。
喷嘴距及混合室喉部长度对喷射器性能综合影响研究
The Comprehensive Influence of the Nozzle Distance and Throat Length of Mixing Chamber on the Performance of Ejector
-
摘要: 在实际环境条件下，采用变几何喷射器是提高喷射器性能的重要途径。基于喷嘴距可调式的喷射器实验测试平台，采用实验和数值分析相结合的方法，研究不同无量纲混合室喉部长度下，无量纲喷嘴距对喷射器的性能及不可逆性影响。结果表明：在测试范围内，喷射器的临界喷射系数随着喷嘴距的增加先增加后减小，当达到最佳无量纲喷嘴距5.83时，临界喷射系数达到最大值0.45，合理的喷嘴距应避免激波链出现在混合室收缩段；喷射器的临界冷凝压力随喷嘴距的增加逐渐下降，当无量纲混合室喉部长度为4时，其临界冷凝压力下降7.36%，下降速度最慢，原因是较大的射流冲量提高了扩压室的压力恢复能力，使流体动能衰减速度有所减缓；当无量纲混合室喉部长度为4，无量纲喷嘴距为5.83时，热力学完善度达到最大值0.21；由喷嘴距引起的喷射器熵产主要来源于两股入口流在混合室收敛段的混合和第二冲击波向扩压室出口的移动，提高扩压室的压力恢复能力能有效避免熵产的急剧增加。研究结果对喷射器的结构及综合性能进行多目标优化，以减小喷射器不可逆性提供重要参考。Abstract: Under the actual environmental conditions, the use of variable geometry ejectors is an important way to improve the performance of ejectors. The effect of dimensionless nozzle distance on the performance and irreversibility of the ejector under different dimensionless throat lengths of the mixing chamber was studied by combining experimental and numerical analysis. The results show that in the test range, the critical entrainment ratio of the ejector first increases and then decreases with the increase of nozzle distance, and when the optimal dimensionless nozzle distance is 5.83, the critical entrainment ratio reaches a maximum value of 0.45. The critical condensation pressure of the ejector decreases gradually with the increase of nozzle distance. When the dimensionless throat length of the mixing chamber is 4, the critical condensation pressure decreases by 7.36%, and the decrease rate is the slowest. When the dimensionless throat length of the mixing chamber is 4 and the dimensionless nozzle distance is 5.83, the maximum thermodynamic perfection reaches 0.21. The entropy production of the ejector caused by the nozzle distance mainly comes from the mixing of the two inlet flows in the convergence section of the mixing chamber and the movement of the second shock wave towards the outlet of the diffuser chamber. Improving the pressure recovery ability of the diffuser can effectively avoid a sharp increase in entropy production. The results provide an important reference for the multi-objective optimization of the ejector structure and comprehensive performance to reduce the irreversibility of the ejector.
-
-
图 3 实验装置图。（a）喷射器性能测试实验平台，（b）喷嘴连接，（c）喷嘴及混合-扩压室
Figure 3. Diagram of the experimental device. (a) Ejector performance test platform, (b) nozzle connection, (c) nozzle and mixer-diffuser chamber
图 5 出口冷凝压力对喷射系数的影响
Figure 5. Effect of outlet condensing pressure on entrainment ratio
图 7 DMTL = 4喷射器在不同喷嘴距的速度云图
Figure 7. Velocity cloud image of ejector at different nozzle distances with DMTL = 4
图 8 DMTL = 4喷射器在不同喷嘴距的质量流量
Figure 8. Mass flow rate of ejector at different nozzle distances with DMTL = 4
图 9 喷嘴距对临界冷凝压力的影响
Figure 9. Effect of nozzle distance on critical condensation pressure
图 10 DMTL = 4喷射器在不同喷嘴距的涡流图
Figure 10. Eddy current diagram of ejector at different nozzle distances with DMTL = 4
图 11 DMTL = 0.5，4，8 喷射器在不同喷嘴距的速度云图
Figure 11. Velocity cloud image of ejector at different nozzle distances with DMTL = 0.5, 4, 8
图 12 临界模式下喷嘴距对热力学完善度的影响
Figure 12. Effect of nozzle distance on thermodynamic perfection under double-chocking mode
图 13 DMTL= 4喷射器在不同喷嘴距的熵产分布
Figure 13. Ejector entropy production distribution at different nozzle distances with DMTL = 4
图 14 DMTL = 0.5, 4 ,8 喷射器在不同喷嘴距的熵产分布
Figure 14. Ejector entropy production distribution at different nozzle distances with DMTL = 0.5, 4, 8
表 1 喷射器固定结构尺寸
Table 1. Fixed structural dimension of ejector
名称 数值 喷嘴长度Lnozzle/mm 83.70 喷嘴入口直径Dnoz,in/mm 13.00 喷嘴喉道直径Dnoz,t/mm 2.00 喷嘴出口直径Dnoz,out/mm 6.20 混合室收敛段长度Lm/mm 56.00 混合室收敛段入口直径Dm/mm 16.67 混合室喉部直径Dm,t/mm 10.80 扩压室长度Ldif/mm 88.00 扩压室出口直径Ddif/mm 23.11 
下载: 导出CSV
表 2 网格独立性检验
Table 2. Grid independence test
网格数 模拟值 实验值 误差% 105622 0.490 0.441 11.1 142530 0.488 0.441 10.6 201175 0.487 0.441 10.4
下载: 导出CSV
表 3 实验和模拟结果对比
Table 3. Comparison of experimental and simulation results
DMTL DNXP 临界喷射系数 临界冷凝压力 实验 模拟 误差% 实验 模拟 误差% 0.5 6.3 0.442 0.480 8.60 7.70 6.60 −14.29 5.19 0.445 0.485 8.99 7.75 7.90 1.94 4.07 0.393 0.445 13.23 9.32 7.75 −16.85 2 6.3 0.440 0.480 9.09 9.15 7.33 −19.89 5.19 0.442 0.485 9.73 9.30 8.00 −13.98 4.07 0.395 0.443 12.15 9.80 8.40 −14.29 4 6.3 0.433 0.475 9.70 8.80 7.60 −13.64 5.19 0.438 0.490 11.87 9.40 8.20 −12.77 4.07 0.400 0.429 7.32 9.50 8.60 −9.47 6 6.3 0.432 0.471 9.03 8.20 7.60 −7.32 5.19 0.441 0.487 10.43 9.00 8.70 −3.33 4.07 0.402 0.430 6.91 9.90 9.50 −4.04 8 6.3 0.430 0.457 6.28 7.60 7.10 −6.58 5.19 0.440 0.488 10.91 9.10 7.30 −19.78 4.07 0.405 0.443 9.49 10.30 8.40 −18.45
下载: 导出CSV
-
[1] 索科洛夫E, 津格尔H M. 喷射器[M]. 北京: 科学出版社, 1977:280(in Chinese) Sokolo E, Zinger H M. Ejector[M]. Beijing: science press, 1977:280 [2] Shestonaloy K O, Huang B J, Peternko V O, et al. Investigation of an experimental ejector refrigeration machine operating with refrigerant R245fa at design and off-design working conditions. Part 1. Theoretical analysis[J]. International Journal of Refrigeration,2015,55:201−211 doi: 10.1016/j.ijrefrig.2015.01.016 [3] Huang B J, Chang J M, Wang C P, et al. A 1-D analysis of ejector performance[J]. International Journal of Refrigeration,1999,22(5):354−364 doi: 10.1016/S0140-7007(99)00004-3 [4] 于文艳, 王海博, 田瑞. 混合室轴向结构参数对蒸汽喷射器性能的影响[J]. 真空科学与技术学报,2018,38(06):455−458(in Chinese) Yu Wenyan, Wang Haibo, Tian Rui. Influence of axial structure parameters of mixing chamber on performance of steam ejector[J]. Journal of Vacuum Science and Technology,2018,38(06):455−458 [5] Dong J M, Hu Q Y, Yu M Q, et al. Numerical investigation on the influence of mixing chamber length on steam ejector performance[J]. Applied Thermal Engineering,2020,174:115204 doi: 10.1016/j.applthermaleng.2020.115204 [6] 王志艳, 胡林静, 席东民, 等. 混合室结构对蒸汽喷射器性能的影响[J]. 真空科学与技术学报,2020,40(02):180−186(in Chinese) Wang Zhiyan, Hu Linjing, Xi Dongmin, et al. Effect of mixing chamber structure on the performance of steam ejector[J]. Journal of Vacuum Science and Technology,2020,40(02):180−186 [7] Wu H Q, Liu Z L, Han B, et al. Numerical investigation of the influences of mixing chamber geometries on steam ejector performance[J]. Desalination:The International Journal on the Science and Technology of Desalting and Water Purification,2014,353:15−20 [8] Yapc R, Ersoy H K, Aktoprako A, et al. Experimental determination of the optimum performance of ejector refrigeration system depending on ejector area ratio[J]. International Journal of Refrigeration,2008,31(7):1183−1189 doi: 10.1016/j.ijrefrig.2008.02.010 [9] Yana J, Caia W, Li Y. Geometry parameters effect for air-cooled ejector cooling systems with R134a refrigerant[J]. Renewable Energy,2012,46:155−163 doi: 10.1016/j.renene.2012.03.031 [10] Tashtoush B, Alshare A, Al-rifai S. Performance study of ejector cooling cycle at critical mode under superheated primary flow[J]. Energy Conversion & Management,2015,94(2):300−310 [11] ESDU. Jet pumps, data item 86030. London: ESDU International Ltd, 1985 [12] Poirier M. Influence of operating conditions on the optimal nozzle exit position for vapor ejector[J]. Applied Thermal Engineering,2022,210:118377 doi: 10.1016/j.applthermaleng.2022.118377 [13] Yan J, Lin C, Cai W J. et al. Experimental study on key geometric parameters of an R134A ejector cooling system[J]. International Journal of Refrigeration,2016,67:102−108 doi: 10.1016/j.ijrefrig.2016.04.001 [14] Chong D T, Hu M Q, Chen W X, et al. Experimental and numerical analysis of supersonic air ejector[J]. Applied Energy,2014,130:679−684 doi: 10.1016/j.apenergy.2014.02.023 [15] Ruangtrakoon N, Thongtip T, Aphornratana S, et al. CFD simulation on the effect of primary nozzle geometries for a steam ejector in refrigeration cycle[J]. International Journal of Thermal Sciences,2013,63:133−145 doi: 10.1016/j.ijthermalsci.2012.07.009 [16] Rusly E, Aye L, Charters W, et al. CFD analysis of ejector in a combined ejector cooling system[J]. International Journal of refrigeration,2005,28(7):1092−1101 doi: 10.1016/j.ijrefrig.2005.02.005 [17] Chen S J, Chen G M, Fang L Y. An experimental study and 1-D analysis of an ejector with a movable primary nozzle that operates with R236fa[J]. International Journal of Refrigeration,2015,60:19−25 doi: 10.1016/j.ijrefrig.2015.08.011 [18] Valle D G J, Jabardo S J, Ruiz C F, et al. An experimental investigation of a R-134a ejector refrigeration system[J]. International Journal of Refrigeration,2014,46:105−113 doi: 10.1016/j.ijrefrig.2014.05.028 [19] Kumar V, Subbarao P M V, Singhal G. Effect of nozzle exit position (NXP) on variable area mixing ejector[J]. SN Applied Sciences,2019,1(11):1−9 [20] 史海路, 刘华东, 魏新利等. 喷嘴距对喷射器及双蒸发压缩/喷射制冷系统性能的影响研究[J]. 高校化学工程学报,2019,33(02):321−328(in Chinese) Shi Hailu, Liu Huadong, Wei Xinli, et al. Study on the influence of nozzle distance on the performance of ejector and double evaporation compression/jet refrigeration system[J]. Journal of Chemical Engineering in Universities,2019,33(02):321−328 [21] Chen H J, Zhu J H, Ge J, et al. A cylindrical mixing chamber ejector analysis model to predict the optimal nozzle exit position[J]. Energy,2020,208:118302 doi: 10.1016/j.energy.2020.118302 [22] Cagri M, Okan G, Ayşe Ugurcan A, et al. Numerical investigation of the flow structures inside mixing section of the ejector[J]. Energy,2018,166:1216−1228 [23] Nguyen V V, Varga S, Soares J, et al. Applying a variable geometry ejector in a solar ejector refrigeration system[J]. International Journal of Refrigeration,2020,113:187−195 doi: 10.1016/j.ijrefrig.2020.01.018 [24] Wang L, Liu J P, Zou T, et al. Auto-tuning ejector for refrigeration system[J]. Energy,2018,161:536−543 doi: 10.1016/j.energy.2018.07.110 [25] Han Y, Wang X D, Sun H, et al. CFD simulation on the boundary layer separation in the steam ejector and its influence on the pumping performance[J]. Energy,2019,167:469−483 doi: 10.1016/j.energy.2018.10.195 [26] Ramesh A, Sekhar J S. Experimental and numerical investigations on the effect of suction chamber angle and nozzle exit position of a steam-jet ejector[J]. Energy,2018,164:1097−1113 doi: 10.1016/j.energy.2018.09.010 [27] Ge J, Chen H J, Jin Y, et al. Conical-cylindrical mixer ejector design model for predicting optimal nozzle exit position[J]. Energy,2023,283:129190 doi: 10.1016/j.energy.2023.129190 [28] Ge J, Chen H J, Li J, et al. Experimental comparison of critical performance for variable geometry ejectors with different mixer structures[J]. Chemical Engineering Journal,2023,478:147487 doi: 10.1016/j.cej.2023.147487 [29] Giorgio B, Nicolò C, Lorenzo C, et al. Computational fluid-dynamics modelling of supersonic ejectors: screening of modelling approaches, comprehensive validation and assessment of ejector component efficiencies[J]. Applied Thermal Engineering,2020,186:116431 [30] Szabolcs V, João S, Rafael L, et al. On the selection of a turbulence model for the simulation of steam ejectors using CFD[J]. International Journal of Low-Carbon Technologies,2017,12(3):233−243 doi: 10.1093/ijlct/ctx007 [31] Chen H J, Zhu J H, Lu W. Optimized selection of one-and two-stage ejectors under design and off-design conditions[J]. Energy Conversion and Management,2018,173:743−752 doi: 10.1016/j.enconman.2018.06.113 [32] Sierra-Pallares J, Valle D G J, Carrascal G P, et al. A computational study about the types of entropy generation in three different R134a ejector mixing chambers[J]. International Journal of Refrigeration,2016,63:199−213 doi: 10.1016/j.ijrefrig.2015.11.007 [33] Herwig H, Kock F. Direct and indirect methods of calculating entropy generation rates in turbulent convective heat transfer problems[J]. Heat and Mass Transfer,2007,43(3):207−215 [34] Sciacovelli A, Verda V, Sciubba E. Entropy generation analysis as a design tool—A review[J]. Renewable and Sustainable Energy Reviews,2015,43:1167−1181 doi: 10.1016/j.rser.2014.11.104 [35] Kock F, Herwig H. Local entropy production in turbulent shear flows: a high-Reynolds number model with wall functions[J]. International Journal of Heat and Mass Transfer,2004,47:2205−2215 doi: 10.1016/j.ijheatmasstransfer.2003.11.025 [36] ANSYS, ANSYS-FLUENT 15.0 User's Guide, 2013 [37] Petrovic A, Svorcan J, Pejcev A, et al. Comparison of novel variable area convergent-divergent nozzle performances obtained by analytic, computational and experimental methods[J]. Applied Mathematical Modelling,2018,57:206−225 doi: 10.1016/j.apm.2018.01.016 -
图( 14) 表( 3)
计量
- 文章访问数: 210
- HTML全文浏览数: 210
- PDF下载数: 2
- 施引文献: 0