R32在光滑管内冷凝的倾角敏感性：传热、液膜再分布与压力梯度

《Next Energy》：Inclination sensitivity of R32 condensation in smooth tubes: Heat transfer, film redistribution, and pressure gradient

　　

管方位是紧凑式换热器设计中的重要考虑因素，因为倾斜管内冷凝性能取决于重力、界面剪切（interfacial shear）与液膜再分布（film redistribution）的相互作用。尽管已有关于倾角对冷凝传热影响的报道，但鲜有研究将R32在中间倾角下冷凝时重力—剪切转变（gravity-to-shear transition）与液膜分布、传热系数（HTC）及压力梯度（pressure gradient）的变化联系起来。本研究采用流体体积法（VOF）结合SST k?ω湍流模型，对R32在内径5 mm光滑管中于11个倾角（β=?90°至+90°，其中负β表示向上流动、正β表示向下流动）下的冷凝进行数值研究。模拟在饱和温度50 °C、质量流率200和400 kg/(m².s)下进行，并以修正Froude数（Fr_so）作为重力—剪切流型指示量解释模拟结果的倾角敏感性。低蒸气干度（vapor quality, χ）和低质量流率下倾角效应最显著：在G=200 kg/(m².s)、χ=0.3时，HTC从β=?90°至+90°降低约37%。提高质量流率会降低HTC对倾角的敏感性，并削弱液膜再分布和压力梯度变化，但在低蒸气干度下强向下倾角仍可引起压力恢复。无量纲Froude数分析结果表明，倾角敏感性随质量流率和蒸气干度增加而降低，这与现有Dobson–Chato判据（Fr_so≈20）所描述的重力影响向剪切主导冷凝转变一致。

研究背景：冷凝器是制冷、化工、储存、电站等过程中的关键部件，其性能受制冷剂质量流率、管径、入口蒸气干度、管方位等因素影响。受空间、泵送或性能增强等约束，冷凝管可能采用不同方位。已有研究报道倾角会影响传热系数（heat transfer coefficient, HTC）、压降和液膜分布，但针对R32在中间倾角下冷凝，并将重力—剪切转变同液膜再分布、HTC和压力梯度同步关联的研究仍有限。R32具有较低全球变暖潜势（GWP=677）、零臭氧消耗潜势（ODP=0），热物性接近R22，被视为R410A等制冷剂的替代选择。现有倾斜管冷凝研究多针对R134a、R1234ze(E)、R245fa等，且对重力促成液膜周向再分布与界面剪切促成环状流之间的竞争解释不足，尤其缺少宽中间倾角范围内R32冷凝的同步分析。因此，研究人员在5 mm光滑圆管中对R32冷凝开展数值研究，考察11个倾角（β=?90°至+90°）、G=200和400 kg/(m².s)及不同蒸气干度下的HTC、液膜再分布和轴向压力梯度，并以修正Froude数（modified Froude number, Fr_so）作为重力—剪切流型指示量。研究结论表明，倾角敏感性在低蒸气干度和低质量流率下最显著；在G=200 kg/(m².s)、χ=0.3时，HTC从β=?90°至+90°降低约37%；提高质量流率可削弱液膜再分布和压力梯度变化。Fr_so分析表明倾角敏感性随质量流率和蒸气干度增加而降低，并与Dobson–Chato判据（Fr_so≈20）一致。其意义在于为R32冷凝器管方位设计及热工水力性能评估提供依据，但原文限定于单一管径、单一饱和温度、两个质量流率以及常壁温和常物性条件。

关键技术方法：采用ANSYS Fluent 2024 R1有限体积法进行瞬态模拟；用流体体积法（VOF）追踪液—气界面，连续表面力（CSF）模型处理表面张力，Lee相变模型（r_v=10 000 s^?1）计算相变源项；采用SST k?ω湍流模型并启用低雷诺数修正；压力—速度耦合采用PISO，界面重构采用Geo-Reconstruct，动量和能量等方程采用二阶迎风。计算域为5 mm圆管，绝热入口200 mm加1000 mm冷凝段，T_sat=323.15 K、P_sat=3.1412 MPa、壁温308.15 K；β取?90°至+90°共11个角度，G=200和400 kg/(m².s)。六面体网格首层高度1 μm，y⁺约5×10^?3–4×10^?1；网格独立性分析选定645 320 cells，GCI（Grid Convergence Index，网格收敛指数）给出的HTC不确定度为0.0435%；与Cao et al.数据验证的最大偏差为3.06%、平均偏差为1.48%。样本队列来源不适用。

研究结果：
3.1 传热系数随蒸气干度和倾角的变化（Variation of HTC with vapor quality and inclination angle）：通过G=200和400 kg/(m².s)下HTC随χ和β变化的曲线发现，低χ时向下倾角通常降低HTC，高χ时响应呈非单调；在G=200 kg/(m².s)、χ=0.3时，β=90°相对水平工况最大降低25%；向上倾角在低χ下可提高HTC，随χ增加界面剪切增强并削弱重力驱动的液膜再分布。G=400 kg/(m².s)时HTC总体更高，倾角引起的偏差更小。
3.2 倾角对HTC的影响（Effect of inclination angle on HTC）：通过HTC随β变化的曲线发现，低χ和低G下倾角影响最大；G=200 kg/(m².s)、χ=0.3时，HTC从β=?90°约4320 W/(m²·K)降至β=90°约2711 W/(m²·K)，降幅近37%；G=400 kg/(m².s)时变化较小；β=?30°在χ=0.15时较水平工况提高10.19%。
3.3 倾角对液膜厚度的影响（Effect of inclination angle on liquid film thickness）：通过周向局部液膜厚度曲线和横截面分布发现，β改变液膜周向分布；G=200 kg/(m².s)、χ=0.3时，β超过60°后膜厚显著降低，β=90°时周向膜厚均匀；G=400 kg/(m².s)、χ=0.5时，175°周向位置处β=45°、60°、90°和?30°的膜厚分别下降约25.84%、34.082%、72.84%和8.25%。液膜厚度变化与HTC变化一致。
3.4 不同倾角下蒸气体积分数云图（Vapor volume fraction contours at different inclination angles）：通过蒸气体积分数云图发现，向上倾角β=?15°至?60°时液膜分布接近水平工况且主要在管底；β=?90°时呈近环状对称分布；向下倾角至60°时变化不显著，β=90°时液膜沿周向均匀分布。
3.5 压力梯度随蒸气干度和倾角的变化（Variation of pressure gradient with vapor quality and inclination angle）：通过总轴向压力梯度随χ和β变化发现，压力梯度随χ增加而增大；向下β降低净压力梯度，向上β增加净压力梯度；G=200 kg/(m².s)时，β=60°和90°在χ=0.3附近出现负压力梯度，表明重力驱动压力恢复；G=400 kg/(m².s)时，β=15°在0.25≤χ≤0.35范围内高于水平工况。
3.6 与已有研究的比较（Comparison with previous study）：通过将G=400 kg/(m².s)、χ=0.6、0.7、0.8下HTC相对水平工况的变化与Gu et al.的R1234ze(E)、40 °C、G=300 kg/(m².s)结果进行定性比较发现，χ=0.6时趋势相反，本研究β=90°最大约增15%，Gu et al.约降10%；原文强调质量流率、饱和温度、管径和物性差异限制制冷剂单独效应判断。
3.7 倾角相关冷凝行为的Froude数判据（Froude number criterion for inclination dependent condensation behavior）：通过Bond数和Confinement数计算得Bo≈59、Co≈0.13，表明重力相对表面张力重要；计算Re_lo、X_tt、Ga和Fr_so发现，G=200 kg/(m².s)时Fr_so在χ=0.4–0.5跨过20，G=400 kg/(m².s)时在χ=0.3–0.4跨过20；Fr_so≤20为重力控制，Fr_so>20为剪切主导，与HTC和膜厚趋势一致。

讨论部分总结：原文将倾角敏感性归因于重力与界面剪切的竞争：低Fr_so时重力驱动周向液膜再分布，高Fr_so时界面剪切增强并降低方位敏感性；Fr_so≈20阈值与Dobson–Chato判据一致，也与Thome et al.流型模型、Xing et al.、Lips and Meyer、Ewim et al.观察到的趋势相符。原文指出该阈值不是通用设计极限，依赖制冷剂物性、管几何和运行条件。局限性包括仅研究一个管径、一个饱和温度、两个质量流率及常壁温和常物性；未系统开展接触角和相变系数的敏感性分析；压力梯度预测未用其他数据独立验证；建议未来实验验证并扩展管径、变物性、壁面润湿性和热阻条件。

研究结论翻译：本研究考察了50 °C饱和R32在5 mm水平圆管不同β下冷凝的HTC、液膜分布和压力梯度。结果表明：（1）HTC随向上β增至?90°而增加、随向下β而降低，低χ和低G下最显著；G=200 kg/(m².s)、χ=0.3时，β=?90°至+90° HTC下降近37%。（2）G由200增至400 kg/(m².s)时，向上和向下流的HTC倾角敏感性均降低，因流动更受剪切主导；蒸气分数云图显示较高质量流率下液膜再分布更小。（3）周向膜厚变化依赖β：重力辅助排液使管底膜变薄、侧部增厚；低χ时向上倾角与水平相近，高χ时偏离。（4）净压力梯度随向下β降低、随向上β增加，低χ最显著；G=200 kg/(m².s)时β=60°和90°在χ=0.3附近出现负压力梯度，表明重力驱动压力恢复；G=400 kg/(m².s)时β=15°在0.25≤χ≤0.35高于水平。（5）HTC和液膜分布的倾角敏感性与Dobson–Chato修正Froude数判据Fr_so=20一致：G=200 kg/(m².s)在χ=0.4–0.5、G=400 kg/(m².s)在χ=0.3–0.4发生转变，倾角效应减弱；这给出物理定量指导，表明在所研究R32冷凝条件下管方位影响有限，并与Del Col et al.对R134a和R32方形微通道的经验判据一致。因此，管倾角显著影响传热、液膜分布和压力梯度，尤其低χ和低G下重力效应超过剪切效应时；结果有助于理解R32冷凝随倾角的行为，并评估冷凝器管方位对热工水力性能的影响。