板壳式热交换器作为一种新型换热器，集合了管壳式换热器和板式换热器的优点，具有高的传热效率、能够耐高温高压、结构紧凑以及质量轻，并且弥补了板式换热器应力分布不均匀的缺点。优良的性能使得板壳式换热器能够适用于多种复杂的工况，如压缩机级间冷却、浮式海上石油平台、有机朗肯循环发电装置等。
目前对于板壳式换热器圆形波纹通道内气-液两相流型与压降的研究鲜有报道，其流道内的两相流相界面宏观结构和各相宏观分布的状态特征都与管内两相流的流动形态有较大差异，与传统的方形板式换热器的波纹通道也有不同，缺乏对板壳式换热器波纹通道内流动结构的多样性、流型转变机制及系统压力分布的认知。本文主要对板壳式换热器圆形波纹通道内气-水两相垂直向上流动流型和压降进行了实验研究。获得了圆形波纹通道内气-水两相流动的宏观相界面分布特性，两相流动压降预测关联式以及两相流型与压降的相关性规律。实验系统如图l所示，实验系统由储液罐、离心泵、空气压缩机、流量计、气液混合器以及相关的连接管段和多种阀门组成。实验测试段为前后两片透明波纹板片交错叠合构成菱形通道(图2)，其材料为环氧树脂，尺寸参数列于表l，为保证实验段的充分发展以及进水和排水的均匀稳定，板片的进口和出口各设置150 mm的延长段。为了获得高清静态相界面特征以及捕捉流型转变特征，采用背光拍摄的方法。在实验结果的观察与流型识别过程中，根据光学原理，在透明波纹板片流道内气相与液相在光反射强度下的差异来区分液体连续相和气体连续相，可以清晰地观察到气相区域颜色是偏暗的，而液相区域视觉上则偏亮(图3)，气-液相界面可以被明显区分，这一方法在流型识别中被广泛采用，图3中的矩形区域表示图像拍摄区域，该区域是长为140mm、宽为110 mm的矩形。压力传感器分别布置在实验段的进出口，采样率为1000 Hz，时间为l0s，以获得压降变化规律；采用高速摄影机和尼微透镜进行流型捕捉以及流型转变的记录，采样频率和拍摄时间分别设置为500帧／秒和5 s。
- 空气压缩机，2.气体减压阀；3。球阀；4.空气流量计；5.液体流量计；6.调节阀；7.气液混合器；8.可视化实验段；9.储液罐；10.离心泵
根据两相流宏观相界面特征，将波纹通道内垂直向上流动的流型划分为泡状流、弹状流、膜状流和搅混流。如图4所示，在泡状流中，流道内液体为连续相，不同大小和形状各异的气泡分布在连续液相内，并与液体一起流动。流道内的气泡和液体的流动方向分为两种：一种是沿着菱形通道的间隙在竖直方向螺旋上升流动；另一种是沿着内外板片波纹方向的沟槽交替流动。在垂直向上流动中，由于重力作用极大地阻碍了竖直方向的螺旋上升流动，通过动态的图像显示，流动主要是沿着内外板片波纹方向的沟槽交替流动，如图4(a)所示沿着内板的沟槽方向流动，如图4(b)所示沿着外板的沟槽方向流动。
图4：泡状流（USL=0.025m/s, USG=0.5m/s）
图5：弹状流（USL=0.05m/s, USG=1.5m/s）随着气速的增加，流道内的气泡会合并成连续的气相，在接触点附近会出现液体滞留区，气速继续增加，整个流道内出现了连续的气相区域，接触点附近的液体滞留区变得更加明显，流型过渡至弹状流。如图5所示，由于垂直向上流动受重力的阻碍作用，使得高速流动的气体将小液团从接触点位置带走的能力减弱，使得液体滞留区这一现象更加明显。由于流体之间的拖拽作用，使得接触点周围的小液团形状各异，并且显示出向上流动的趋势。随着气速的增加，流道内部分区域的液体滞留区会消失，部分会逐渐拉长成为细长的液团，气速继续增加，流道内的液体滞留区全部消失，在流道内形成连续的气相，液体在内外波纹板的壁面以液膜的形式随气芯快速流动，流型过渡至膜状流。如图6所示，膜状流中气体流量较大，液体流量较小，空气在中问流道流动，液体除了以液膜的形式在波纹板的壁面流动外，还有极少量液体以夹带液滴的形式存在于流道中。在垂直向上流动中，重力作用阻碍液膜沿着壁面向上流动，因此液膜的流动速度较小，厚度较厚。
图6：膜状流（USL=0.075m/s, USG=5.0m/s）随着液体流量的增大，膜状流中的夹带液滴出现得更加频繁，导致液膜剧烈波动并且气体流通截面积减小，当膜状流界面上的剪切力远大于表面张力时，会形成液块并破坏膜状流中间的气芯，流型过渡至搅混流。如图7所示，一部分液体以液膜的形式沿着波纹板的壁面流动，剩余液体与气体在中问流道混合流动，在搅混流中，气液两相剧烈搅拌，很难区分出气液相界面，由于垂直向上流动会产生更剧烈的振荡，在波纹槽道内能看到很多小气泡。
图7：搅浑流（USL=0.175m/s, USG=5.0m/s）降板壳式换热器中的单相压降包括摩擦压降、重力压降和压力损失。总压降表示为
图8显示了单相流动中空气和水的Reynolds数与单相摩擦因子之间的关系，结果表明单相空气和水的摩擦因子都随着Reynolds数的增加逐渐减小。实验结果与其他文献有着类似的趋势，数值的差异可以归因于不同的板片结构和板片的几何参数。基于实验结果，给出了单相空气和水的摩擦因子与Reynolds数的关联式。
图8：单相流动中Reynolds数与摩擦因子的关系
图9显示出了气-液两相摩擦压降与表观气速和表观液速的关系。结果表明，表观气速和表观液速都对两相摩擦压降有较大的影响，随着表观气速和表观液速的增大，垂直向上流动会导致流道内产生较强的振荡，湍流也变得更加剧烈，两相摩擦压降也增大。本文采用均相流模型计算两相摩擦压降。两相流的混合密度由式(8)计算
式中，pm、pl、pg、α和x分别是气液两相的混合密度、液体密度、气体密度、含气率和气体干度。
Lockhart—Maninelli方法通常用于计算两相摩擦压降，计算式如下
式中，C值的大小取决于两种流体的不同流态。Chisholm给出了光滑管内常数C的建议值，液体湍流-气流湍流，C=2l；液体层流-气体湍流，C=12；液体湍流-气体层流，C=10；液体层流-气体层流，C=5。两相摩阻系数与Maninelli参数x之间的关系如图10所示。结果表明，随着Maninelli参数的增大，两相摩阻系数逐渐减小。基于实验结果，提出了一个新的关联式用于预测板壳式换热器波纹通道垂直向上流动的两相摩擦压降，如式(14)所示。大多数预测数据位于±30％偏差范围内。
图10：Martinelli参数随两相摩阻系数的变化图11显示了实验结果与板式换热器流道中内两相摩阻系数与Martinelli参数之间关系的比较。结果表明，板壳式换热器圆形波纹通道内两相摩阻系数与Martinelli参数的关系与板式换热器内具有较好的一致性。
图11：两相摩阻系数与Martinelli参数的文献对比图12描述了板壳式换热器波纹流道内四种典型流型与其压降波动幅值的变化规律。结果表明：泡状流由于表观气速与表观液速较小，压降波动的幅度很小，基本维持在6.4~8.5 kPa范围内波动；随着表观气速的增加，小气泡聚集合并形成弹状流，此时压降波动的幅度有所增加，基本在7.1～13.3 kPa范围内波动；随着表观气速的增加形成膜状流，膜状流流动较为稳定，压降波动幅度不大，压降基本在24.6~35.5 kPa范围内波动；当液体流量增大时形成搅混流时，此时气-液两相在流道内扰动剧烈，在垂直向上流动伴随有较强的振荡，湍流增强，此时的压降波动幅度最大，基本在34.2～56.1 kPa范围内波动。通过以上分析，可以看出不同的流型其压降的大小以及压降的波动幅度都有较大的差异，因此通过不同流型的压降波动幅度的变化可以为流型判别及其转变提供依据。本文对板壳式换热器垂直向上流动的流型和压降进行了实验研究，主要结论如下。
(1)在本文波纹通道几何条件下，垂直向上气-液两相流动的流型被划分为泡状流、弹状流、膜状流和搅混流四种；此外，发现泡状流的流动是沿着内外波纹板的沟槽方向交替流动。(2)在本文波纹通道几何条件下，基于Lockhart—Martinelli理论获得了气一液两相流动的两相摩阻系数与Maninelli参数的关系，并拟合了两相压降预测关联式，发现chisholm参数c的值与Chisholm最初建议的管内层流-层流的值接近。(3)在本文波纹通道几何条件下，通过流型压降波动曲线得出泡状流中的压降波动幅值最小，弹状流与膜状流次之，搅混流中压降波动幅值最大，其压降波动的幅值可以为流型的种类判别及其转变提供一定的依据。
版权声明：本文内容参考《绝热状态下板壳式换热器壳侧的流型与压降》，作者：吴鹏飞等，由制冷空调换热器技术联盟编辑整理，转载请注明来源。
