茹卡乌斯卡斯横掠顺排管束的数值模拟与公式验证教程.doc

茹卡乌斯卡斯横掠顺排管束的数值模拟与公式验证 周梦 （南京工业大学 化工过程机械 652080706034） Abstract: Zhukauskas experimental correlation has been widely used in heat exchanger design. Fluid for Zhukauskas ACROSS LINE BUNDLE test section flow and heat transfer characteristics of the prototype, after appropriate simplification, to establish a two-dimensional model, using the Fluent CFD software in large model fluid flow and heat transfer characteristics of numerical simulation. The numerical results are compared with the formula, and examine the factors of inlet velocity, inlet temperature, the heat transfer tube wall temperature and the outside diameter of the heat transfer performance of heat transfer tubes (Nusselt number) and use Zhukauskas formula conducting impact error of calculation. Keywords: numerical simulation,ACROSS LINE BUNDLE formula ,Nu number 流体外掠管束的流动比较复杂,其换热和流动特性主要采用实验数据和实验关联式进行预测。在众多应用于换热器设计的实验关联式中,茹卡乌斯卡斯公式有着非常广泛的应用[1]。本文研究Re数在，流体外掠顺排管束为20排的实验模型，根据计算流体力学理论，运用大型通用流体力学计算软件Fluent对实验模型内流体的换热和流动特性进行数值模拟，并将数值模拟结果与公式计算结果进行误差分析。 11.1 模型与参数 模型为干空气横掠顺排管束，示意图见图1几何尺寸及操作条件见表1。 图1 模型正面示意图 几何尺寸 操作条件 干空气50℃物性参数 模型长度L/mm 520 进口速度V/(m/s) 200 密度ρ/(kg/m3) 1.093 模型宽度H/mm 104 进口温度T/℃ 50 热容cp/(J/(kg*K)) 1005 换热管外径D/mm 20 管壁温度Tw/℃ 20 热导率λ/(W/(m*K)) 0.0283 中心距S1、S2/mm 26 流动介质 空气 动力粘度v106/(m2/s) 17.95 1.2 本问题所涉及的流动为二维稳态定常流动, 采用SIMPLE方法求解 N-S方程，湍流模型采用标准k-ε模型。[1-4]方程式如下： 连续性方程 动量方程 k方程 ε方程 能量方程 其中 式中： x为坐标；u为速度分量；i，j为张量指标；t为时间；k为湍流脉动动能；ε为脉动动能耗散率；T 为温度；ρ为流体密度；μ为动力粘度；cｐ为定压比热容；μｔ 为湍流粘度；σｔ为能量普朗特数，一般取0.85 ； σk为脉动动能普朗特数，一般取 1.0；σ ε 为脉动动能耗散率普朗特数，一般取1.3；模型系数cμ ，c1和 c2分别取值0.09，1.44和1.92。 1.3 网格划分 对计算区域采用非结构化壳网格划分，以四边形网格为主，允许一部分三角形网格单元的存在；在为了准确的计算流体在管束壁面边界层内的变化，在管子外划分4层的边界层，底层网格高度为0.5mm，增长率为1.1。全局网格如图2所示，部分网格如图3所示。 图2 全局网格显示 图3 部分网格显示 1.4模型计算边界条件设定如下：管束壁面温度为恒温20℃；流体进口温度为50℃；流体入口采用速度进口；出口采用自由流动；为避免两侧壁面的对体流动产生影响，对轴向中间截面采用对称边界。 1.5在进行模拟时,湍流模型选用标准k-ε模型，壁面处流体无滑移。控制方程采用有限体积法离散，压力速度耦合采用SIMPLE算法,压力插补格式选择standard格式,其他项采用一阶迎风格式。 22.1 网格独立性分析 为了检查数值计算的准确性，分别对网格数目为5380、8804、16061、36541、49499进行独立性检验。模拟计算结果如图4