板翅式及管翅式换热器气流湍流特性研究

2020-01-10徐大川廖达雄裴海涛秦红岗

实验流体力学 2019年6期

符澄, 赵波, 徐大川, 廖达雄, 裴海涛, 朱博, 秦红岗

(1. 中国空气动力研究与发展中心空气动力学国家重点实验室, 四川绵阳 621000; 2. 中国空气动力研究与发展中心设备设计及测试技术研究所, 四川绵阳 621000)

0 引言

随着换热器设计技术的进步，圆管、椭圆管翅片式等新型换热器在风洞中不断成功应用。作为连续式风洞的关键部段之一，换热器的设计除了要兼顾换热效率及压力损失之间的关系之外，其对风洞气流的品质的影响也逐渐被人们所认识和了解[1-3]。

风洞中通常使用的换热器如管翅式及板翅式换热器，在结构形式上具备将上游大尺寸旋涡分割成小尺寸旋涡的能力，同时，由于翅片间距小、沿流向长度较长等特点，大大增加了气流流过时的粘性耗散作用，具备降低风洞气流湍流度、提高换热器下游截面流动均匀性及稳定性的能力，如图1所示。赵波等[4]采用数值模拟方法，对错排椭圆管翅片式换热器下游的速度及温度分布进行了优化，并在连续式跨声速风洞中成功应用。楚攀等[5]对椭圆管翅片式换热器的增强混合及换热效果进行了研究，提出了一种改善速度场和温度场之间协同关系的流动控制方法。刘晓波等[6]综合比较了管翅式、板翅式等换热器在风洞中应用的优势与劣势，给出了换热器在大型风洞中布局的优化方案。魏进家等[7]对板翅式换热器的流动及传热特性进行了研究，为板翅式换热器的设计提供了指导。关欣等[8]则对板翅式换热器的动态特性通用分析解进行了研究，获得了较好的计算精度。上述研究成果直接推动了大型连续式风洞换热器的技术进步并扩大了其应用范围。

图1 换热器试验件

板翅式换热器中的长扁管及管翅式换热器中的圆管/椭圆管、翅片均会产生较强的尾流，该尾流影响范围大，会产生较强的再生湍流，其对换热器气流湍流特性的影响具有较大的不确定性。由于换热器下游的湍流度及其分布特性对风洞试验段流场的动态品质具有较大影响，因此通过研究了解和掌握其特性，对风洞换热器的设计具有重要意义。本文在不考虑不同形式热交换器换热效率的前提下，通过试验研究和数值模拟的手段对不同来流湍流度条件下，板翅式及管翅式换热器下游的湍流特性进行了研究，分析了2种换热器在湍流流场整流能力方面的差异，可为高气流品质要求的风洞换热器设计提供参考。

1 试验设备及方法

1.1 换热器模型

换热器试验件沿气流方向入口截面尺寸为550 mm×400 mm。板翅式换热器的扁管长208 mm，厚8 mm；采用矩形翅片，高19 mm，管间距46 mm，翅片厚度0.25 mm，翅片间距2.3 mm，沿气流方向单排布置。

椭圆管翅片式热交换器的椭圆管外长轴36 mm，外短轴14 mm；横向管间距54 mm，排间距28 mm；翅片长55 mm，宽26 mm，厚0.3 mm，翅片间距3.5 mm，沿气流方向交错4排布局。2种热交换试验件的外形及结构见图1。

1.2 试验设备

在0.55 m×0.40 m航空声学风洞中对换热器的湍流特性进行研究。该航空声学风洞开口试验段最大风速100 m/s，空风洞闭口试验段入口气流湍流度小于0.05%，背景噪声78 dB(A)(v=80 m/s)[9-10]。风洞的气动轮廓如图2所示。

图2 风洞气动轮廓

换热器试验件在声学引导风洞试验段内的安装如图3所示。从收缩段出口开始，依次为扰流体组件、换热器试验件、平直段及收集器。其中，扰流体距离换热器试验件入口的距离为250 mm，换热器出口平直段的长度为300 mm。图4给出了单片扰流体的外形。

使用可变角度的扰流体组件模拟产生不同湍流度及速度不均匀的来流条件。使用Dantec公司的热线风速仪对湍流度和速度进行测量，所使用的热线探头为二维探头，可分别对气流的横向和流向速度以及速度脉动进行测量[11-12]。热线探头移测架安装在风洞驻室地面上，通过一个L形的支架伸入换热器下游，并通过移测架的移动实现不同位置的湍流度测量。

图3 试验件及测试设备的安装及尺寸

图4 扰流体组件外形

1.3 测试条件及方法

在换热器入口及出口的风洞试验段壁面上，四壁开静压孔，4个静压孔汇成一路，使用压力扫描阀分别测量获得换热器入口及出口的平均静压。入口、出口的静压差即为换热器的总压损失，损失系数可表示为：

(1)

其中pin_ave、pout_ave分别为换热器入口及出口的平均静压，ρ为空气密度，v为换热器入口气流速度。

在换热器出口下游截面上，测量沿垂直方向及水平方向各一条直线上的湍流度(轴向和横向)和速度(轴向和横向)。测点等间距分布，垂直方向(y方向)间距为5 mm，垂直方向测量总高度340 mm；水平方向(z方向)间距为10 mm，水平方向测量总宽度400 mm。测量区域均沿风洞轴线对称分布。测试过程中，热线风速仪的采样频率为20 kHz；试验段入口风速分别为6，8，10,12和14 m/s；扰流体偏转角度为0°、20°和40°；湍流度及速度测量截面位于距离扰流体组件下游1000、1100和1200 mm处，分别定义为x=0、100和200 mm截面位置。

2 试验结果及分析

对换热器下游的湍流度分布进行测量的目的是比较不同构型的换热器对不同湍流度来流的整流效果，该整流效果主要体现在降低湍流度以及对流动不均匀性的抑制方面。

在没有安装换热器的条件下对不同偏转角度的扰流体组件下游湍流度及速度分布进行了测量。结果表明，来流风速对扰流体下游湍流度及速度分布影响较小。

图5给出了v=14 m/s条件下、x=200 mm截面垂直方向的流向湍流度测试结果。从图中可以看出，扰流体偏转角度为0°时，下游的湍流度分布表现出明显的周期性变化的特征；随着偏转角度的增加，扰流体下游的湍流度逐渐增大。扰流体尾流随着偏转角度的增加而不断向下偏转，垂直方向的湍流度分布的不均匀程度增加，从而呈现出轴线下部的湍流度低、上半部分湍流度高的现象，而水平方向的湍流度分布则比较均匀。试验结果表明，扰流体尾流所产生的湍流度的横向分量略高于其在流向上的分量。

图5 扰流体下游湍流度分布

随着来流湍流度的增加，换热器的阻力损失系数呈减小的趋势。图6给出了2种换热器在不同扰流体偏转角度情况下，损失系数随速度变化的趋势。扰流体40°偏角时的阻力损失比0°时减小了约10%。同时可以看出，相同来流湍流度条件下，随着来流速度的增加，雷诺数变大从而导致翅片附面层厚度减小，换热器的阻力损失系数逐渐降低。速度对压力损失的影响在板翅式换热器上更加明显，与6 m/s的来流速度相比，流速为14 m/s时的损失系数降低了36%。

图7分别给出了v=14 m/s、扰流体偏转角度为40°、x=200 mm截面垂直及水平方向湍流度流向分量及横向分量的测量结果。从图中可以看出，在扰流体偏转角度为40°时，椭圆管翅片式换热器表现出了较好的整流效果，扰流体下游的湍流强度得到了抑制，同时，原来强烈的不对称分布也到了极大的改善。通过比较椭圆管翅片式换热器下游湍流度流向和横向分量分布特性可以看出，2个分量均降至5%～6%的水平，因此可以判断，椭圆管翅片式换热器对流向及横向湍流强度的抑制效果是相当的。

图6 换热器的压力损失系数

图7 换热器下游湍流度分布

与椭圆管翅片式换热器类似，在扰流体偏转角度为40°时，板翅式换热器下游的湍流强度及不对称分布均得到了明显的抑制。但其下游流向与横向湍流度分布有较大差别，横向湍流度明显小于流向湍流度，且横向湍流度均匀性更佳。由此可以判断，单排板翅换热器对横向湍流强度的抑制效果优于对流向湍流强度的抑制效果。总体来说，板翅式换热器的湍流度抑制效果优于椭圆管翅片式换热器，但椭圆管翅片式换热器下游的气流湍流度分布更均匀。

造成上述整流效果差异的主要原因是板翅式换热器和椭圆管翅片换热器在结构上的巨大差异。由于椭圆管翅片式换热器内部独立的椭圆管结构，不断破坏气流在换热器内部的附面层增长，从而让湍流度分布更趋于均匀；而板翅式换热器基管与翅片形成独立气流通道的结构形式则使得不同基管之间没有压力及涡量的传递，对来流不均匀性的抑制要明显差于椭圆管换热器，但得益于基管和翅片所分割出来的独立的气流通道，使得板翅式换热器的整流特性类似于风洞中常用的蜂窝器整流装置[13]，对横向湍流强度有更好的整流效果。

当来流的湍流度较高时，换热器表现出较好的整流效果。但当来流湍流度较低时，由于换热器自身的再生湍流度高于来流湍流度，因此，换热器在低来流湍流条件会表现出增强湍流度的效果。图8给出了扰流体偏转角度分别为0°、20°和40°条件下，椭圆管换热器下游的流向湍流强度分布。可以看出，在这3种条件下，换热器下游的湍流度分布特性基本保持一致，这说明即使来流湍流度有较大幅度的变化，换热器也可以表现出较好的整流效果。

图8 不同来流条件下的换热器下游湍流度分布

Fig.8TurbulentintensitydistributiondownstreamtheHEXunderdifferentconditions

3 数值模拟结果及分析

采用有限体积法分别对板翅式及管翅式换热器的流场进行三维数值模拟，获得不同类型换热器的再生湍流度特性，从而进一步评价这2种换热器的整流效果。

三维坐标系下守恒形式的控制方程为：

(2)

其中t为时间，Q为守恒变量矢量，F、G和H为无粘通矢量，Fv、Gv和Hv为粘性通矢量。压力与速度的耦合方程采用SIMPLE算法求解，湍流模型采用SST模型。

计算网格采用结构化网格，管壁和翅片壁面边界层做加密处理，壁面第一层网格高度0.1 mm，对应的y+在1～5之间，通过比较网格加密程度做网格无关性分析，在此不再赘述。入口采用速度边界条件，出口采用压力边界条件，其余流体边界设为周期性边界条件，壁面采用无滑移边界条件。计算模型见图9。

图9 换热器数值计算模型

图10给出了不同结构型式换热器流向中心线上的总湍流度分布。从图中可以看出，在给定计算域入口湍流度为38%的情况下，椭圆翅片管换热器出口截面下游100 mm处(x=0.2 m处)的湍流度因尾流的影响达到26%左右，随后在粘性耗散的作用下呈现逐渐降低的趋势，在下游600 mm处湍流度降至约为9%；板翅式换热器出口截面下游100 mm处(x=0.3 m处)湍流度约为9%，下游600 mm处的湍流度则降至约3%。因此可以判断，板翅式换热器的再生湍流度仅为椭圆管翅片式换热器的1/3左右。图11给出了板翅式换热器和椭圆翅片管换热器下游600 mm处的湍动能数值模拟结果。可以看出，椭圆翅片管换热器下游截面上的湍动能最大值约为板翅式换热器的3～4倍。