圆管内CaSO析晶污垢模型与数值模拟..docVIP

圆管内CaSO4析晶污垢模型与数值模拟 徐志明　张进朝 （东北电力大学，吉林　吉林　132012） 摘　要　本文从传热传质的角度建立了圆管内CaSO4析晶污垢形成过程的数学模型，并进行了数值模拟。根据模拟得出的圆管内温度场、速度场、以及CaSO4的浓度场，进而结合污垢模型计算出了CaSO4析晶污垢的沉积率、剥蚀率以及污垢热阻随时间的变化规律。最后将计算结果和实验结果进行了比较，证明了该模型的正确性。 关键词　析晶污垢；数值模拟；CaSO4 中图分类号：TK124　　　文献标识码：A　　　文章编号： The Model and Numerical Simulation of CaSO4 Crystallization Fouling in the Tube XU Zhi-Ming　ZHANG Jin-Chao (Northeast　Dianli　University，Jilin　132012，China) Abstract　 A mathematical model of CaSO4 crystallization fouling process in the tube is proposed based on the theory of heat and mass transfer. A numerical simulation calculation was also carried out. Then the changes of the CaSO4 crystallization fouling deposition rate, the erosion rate and the fouling resistance with time were achieved by the fouling model, with the calculated temperature field, velocity field, and the CaSO4 concentration field in the tube. Finally, the computed results and experimental results are compared. The correctness of this model is proved. Key　words　Crystallization Fouling；Numerical Simulation；CaSO4 前　言 析晶污垢是指在流动条件下呈过饱和的流动溶液中溶解的无机盐淀析在换热面上的结晶体。一般正常溶解度的盐类淀析在冷却面上，而具有反常溶解度的难溶或微溶盐类则淀析在加热面上[1]。本文所研究的CaSO4析晶污垢就是具有反常溶解度的微溶盐类。 对于析晶污垢而言，目前的模型已有多种，但几乎所有的研究者都是以Kern和Seanton提出的微分方程为基础进行研究的[1]。由于污垢的形成过程是涉及质量交换、热量交换和动量交换等复杂的动态过程，影响这一过程的因素很多，如流体性质、壁温、流体流速、壁面状况、流体与壁面的温度梯度、湍流强度、流体与固体壁面的剪切应力等，所以析晶污垢一直没有一个让大家所认同的统一数学模型。对于析晶污垢的形成机理研究，杨传芳等[2]提出了一个预测饱和水系统中碳酸钙析晶污垢沉积率预测模型，并对其诱导期内和以后的沉积机理进行了探讨。Taborek等[3]提出了一个预测冷却水系统的析晶污垢，基本思路仍然是按照沉积与剥蚀分别考虑表达式的。Hasson等[4]在没有剥蚀过程的前提下提出了一个考虑离子组分输运扩散的模型。上述模型都是以宏观参数为基础进行预测的。Brahimd等[5]首次将传热传质基本原理与以Kern和Seanton为基础的污垢模型结合建立了一个矩形截面槽道中两个侧面的污垢形成过程的数学模型，并进行了数值模拟和实验验证。但绝大多数换热设备的基本元件是管。本文以二维圆管道中的析晶污垢为研究对象，在前人研究的基础上，建立了CaSO4析晶污垢形成过程的数学模型，并进行了模拟计算。在不考虑诱导期的情况下，从污垢沉积质量出发计算出了污垢热阻随时间的变化规律，并通过实验证实了本模型的准确性。 物理模型 物理模型采用内直径d=22mm的光管，为避免进口段的影响，根据尼古拉兹试验紊流进口段管长度L*≈（25~40）d的约束，取管长L=1500mm。管内工质为浓度cf=1.2kg/m3的CaSO4溶液。假设流动、污垢等诸特性参数在各个方向上都是相同的，即各向同性且均匀分布，可将模型简化为二维模型。物理平面模型及各区域分布如图1所示： 图1 物理平面模型图 数学模型及边界条件 要得到整个计算区域的速度场、温度场、浓度场以及污垢的沉积率需要将连续性方程、动量方程、能量方程以及传质方程联