螺纹换热管.docx

研究生课程论文
(2012-2013学年第二学期)

螺纹换热管数值模拟与分析
摘要：利用数值模拟方法，考察单个螺纹换热管不同结构形式时流道内的换热特点，计算得到压力、速度、湍动能、湍流强度、温度等分布云图并做分析。流体流过螺纹管时受到螺纹槽的影响，使流动受到扰动，在提高热交换的同时也会产生较大的能量损失，计算比较模型的能量损失。
关键词：螺纹换热管；数值模拟；换热特点；能量损失
换热管是换热器的元件之一，置于筒体之内，用于两介质之间热量的交换。除光滑管外，换热器还可采用各种各样的强化传热管，如翅片管、螺纹管、螺旋槽管等。当管内直径两侧给热系数相差较大时，翅片管的翅片应布置在给热系数低的一侧。常用材料有碳素钢、低合金钢、不锈钢、铜、铜镍合金、铝合金、钛等。此外还有一些非金属材料，如石墨、陶瓷、聚四氟乙烯等。设计时应该根据工作压力、温度和介质腐蚀性等选用合适的材料。采用小管径时，可使单位体积的传热面积增大、结构紧凑、金属耗量减少、传热系数提高。但小管径流体阻力大，不便清洗，易结构堵塞。一般大直径管子用于粘性大或污浊的流体，小直径管子用于较清洁的流体。
螺纹管又称低肋管，主要是靠管外肋化扩大传热面积，一般用于管内给热系数比管外给热系数大1倍以上的场合。对于管外冷凝及沸腾，由于表面张力作用，也有较好的强化作用，如下图1所示。

图1螺纹换热管
由于进行实验研究具有周期长，耗时费力的缺点，而数值模拟方法具有方便、灵活，研究和开发周期相对较短，费用较低、限制减少等优点。因此本文利用数值模拟方法，考察单个螺纹换热管不同结构形式时流道内的换热特点。流体流过螺旋槽纹管时受到螺旋槽的影响，使流动受到扰动，在贴近壁面处产生较强的涡流，因而在提高热交换的同时也会产生较大的能量损失。
1理论基础
数值模拟方法遵循质量守恒、动量守恒和能量守恒这三个基本规律，其理论基础可以用以下方程式表示。
连续性方程：

动量方程：

能量方程：

对于这类的流动损失，可将其看作一个较长的管道，则由达西公式得到沿程阻力系数表达式为
λ=D2ΔpρVin2DL
式中p为螺纹管两端的压差，为流体的密度。
对于水力光滑管来说。在5000<Re<30000范围内有勃拉休斯公式
λ=0.3164Re0.25
将上述式子计算的结果对比，可判断螺纹换热管的流动阻力特性。
对于水力光滑管道，传热经验公式为
Nu=0.023Re0.8Pr0.4(20度水Pr=7.02)
对于槽纹管，有
Nu=165aD13(tD)-0.5(Re104)(0.8-3.5aD)Pr0.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认为收敛