-螺纹换热管

研究生课程论文

(2012-2013学年第二学期)

螺纹换热管数值模拟与分析

摘要：利用数值模拟方法，考察单个螺纹换热管不同结构形式时流道内的换热特点，计算得到压力、速度、湍动能、湍流强度、温度等分布云图并做分析。流体流过螺纹管时受到螺纹槽的影响，使流动受到扰动，在提高热交换的同时也会产生较大的能量损失，计算比较模型的能量损失。

关键词：螺纹换热管；数值模拟；换热特点；能量损失

换热管是换热器的元件之一，置于筒体之内，用于两介质之间热量的交换。除光滑管外，换热器还可采用各种各样的强化传热管，如翅片管、螺纹管、螺旋槽管等。当管内直径两侧给热系数相差较大时，翅片管的翅片应布置在给热系数低的一侧。常用材料有碳素钢、低合金钢、不锈钢、铜、铜镍合金、铝合金、钛等。此外还有一些非金属材料，如石墨、陶瓷、聚四氟乙烯等。设计时应该根据工作压力、温度和介质腐蚀性等选用合适的材料。采用小管径时，可使单位体积的传热面积增大、结构紧凑、金属耗量减少、传热系数提高。但小管径流体阻力大，不便清洗，易结构堵塞。一般大直径管子用于粘性大或污浊的流体，小直径管子用于较清洁的流体。

螺纹管又称低肋管，主要是靠管外肋化扩大传热面积，一般用于管内给热系数比管外给热系数大1倍以上的场合。对于管外冷凝及沸腾，由于表面张力作用，也有较好的强化作用，如下图1所示。

图1螺纹换热管

由于进行实验研究具有周期长，耗时费力的缺点，而数值模拟方法具有方便、灵活，研究和开发周期相对较短，费用较低、限制减少等优点。因此本文利用数值模拟方法，考察单个螺纹换热管不同结构形式时流道内的换热特点。流体流过螺旋槽纹管时受到螺旋槽的影响，使流动受到扰动，在贴近壁面处产生较强的涡流，因而在提高热交换的同时也会产生较大的能量损失。

1 理论基础

数值模拟方法遵循质量守恒、动量守恒和能量守恒这三个基本规律，其理论基础可以用以下方程式表示。

连续性方程：

动量方程：

能量方程：

对于这类的流动损失，可将其看作一个较长的管道，则由达西公式得到沿程阻力系数表达式为

λ=D 2Δp ρV in 2D L 式中Δp 为螺纹管两端的压差，ρ为流体的密度。

对于水力光滑管来说。在5000

λ=0.3164

R e 0.25 将上述式子计算的结果对比，可判断螺纹换热管的流动阻力特性。

对于水力光滑管道，传热经验公式为

N u =0.023R e 0.8P r 0.4(20度水P r =7.02)

对于槽纹管，有

N u =165(a D ⁄)13⁄(t D ⁄)−0.5(R e 10⁄4

)(0.8−3.5a D ⁄)P r 0.4

2 物理模型和数学描述

本文采用有限元分析软件FLUENT 软件，以三种不同螺纹结构的换热管作为模拟对象，考察换热管流道中的流体流动和传热特性，以及比较某一类型的螺纹换热管和光滑换热管的。本文的模型尺寸分别采用（1）光滑换热管，管道长Lo=200mm ，管径D=20mm （2）螺纹换热管，管道长Lo=200mm ，螺纹槽段长L=100mm ，管径D=20mm ，螺纹槽螺距T=5mm ，螺旋槽的深度a=2mm ，宽度b=2mm 。（3）螺纹换热管，管道长Lo=200mm ，螺纹槽段长L=100mm ，管径D=20mm ，螺纹槽螺距T=5mm ，螺旋槽的深度a=2mm ，宽度b=1mm 。

换热管的材料选用铜铝合金，流体介质采用稳态不可压缩的水，假设流体进口温度为360K ，计算采用k-epslion 紊流模型，边界条件采用速度进口和压力出口。壁面的恒温温度设置为300k 。在求解的过程中，当连续性方程、动量方

程和能量方程中的变量残差都达到10-6认为收敛。

3网格划分

在进行数值模拟的过程中，网格的划分都是最关键的工作之一，网格的好坏直接影响着计算速度和计算结果的收敛。本文的模型结构比较规则，所以采用网格质量好、计算速度快、内存占用量少的额六面体结构化网格分别对流体和固体部分进行划分，网格的疏密程度对模型计算结果也有很大的影响。计算之前首先对网格进行独立性检测，采用不同的网格密度重复同样的计算，比较所得的结果，当模型进一步细化网格，在工程允许的偏差范围内数值解已经几乎不再发生变化，即认为此解为网格独立的解。以模型（2）为例，在计算网格加密后，重新分析计算结果，通过比较结果发现，网格数量为10左右时能够满足计算的精度要求。三个模型的网格见图2。

图2 光滑管、螺纹宽2mm、1mm的螺纹管

4边界条件设置

本文的模型厚度不计，流体介质采用稳态不可压缩的水，计算采用

k-epslion紊流模型，边界条件采用速度进口和压力出口。假设流体进口温度为360K，选择强度和水利直径，速度为1m/s。我们将这个模型简化，不考虑模型的厚度，在建模的过程中将壁面画成一个面，假设该面为恒定壁温为300K。在

求解的过程中，当连续性方程、动量方程和能量方程中的变量残差都达到10-6认为收敛。

5 求解计算

进行初始化设置后，对模型迭代计算，迭代次数设置为1000次。如下面三张残差监测曲线图。我们可以得出，虽然其结构不同，有细微的差别，但是计算都是收敛的。

图3 光滑管图4 螺纹管螺纹宽2mm

图5螺纹管螺纹宽1mm

6分析与讨论

为了清晰观察三维流动，需要用某些特殊面上的流动状态来显示整个流动的特征，本文创建 z=0为例，展示了整个管子的流动特性（流体从管子右边流入）。

下面分别是模型（1）（2）（3）的各个参数的分布云图，我们可以观察到z=0平面的压力、速度、湍动能、湍流强度、温度的变化趋势以及不同模型之间的差别。

图6 压力分布图7 速度分布

图 8 湍动能分布