数值分析方法在暖通专业的应用.docx

数值计算方法作业背景介绍：我是土木工程学院，供热供燃气通风及空调专业的学生，简称暖通空调专业， 即本科阶段所称的建筑环境与能源应用工程专业本专业有四大基础课程：工程 热力学、工程传热学、工程流体力学和建筑环境学在工程流体力学中有这样的 几类问题：工业管道内的流动计算，包括（1）计算水头损失；（2）计算管道内 的流量；（3）计算工业管道管径第一类问题最简单，一般给定了管径和流量或者流速由管径可以得到相对 粗糙度，由流量可以得到流速，再得到雷诺数，由公式或莫迪图确定摩擦系数， 再得到水头损失第二类问题需要迭代计算，一般给定管径和水头损失先假定摩擦系数，计 算得到流速，再得到雷诺数，由管径得到相对粗糙度，再用公式或者查莫迪图确 定摩擦系数，比较摩擦系数，重新选取摩擦系数再重复计算，直到摩擦系数基本 一致，再由流速确定流量第三类问题也需要迭代计算，一般给定流量和水头损失假设摩擦系数，由 水头损失和流量可以得到管径，再计算雷诺数和相对粗糙度，由公式或者莫迪图 重新确定摩擦系数，比较摩擦系数，再进行计算，直到摩擦系数基本一致，最终 得到管径以上第二类和第三类问题需要用到数值分析中的迭代法计算 下面给出并简单介绍莫迪图。

莫迪图又称为摩擦系数计算曲线图工业管道内摩擦系数的计算非常重要， 在莫迪图出现之前，计算摩擦系数有非常多的计算公式，这些计算式在不同场合 下应用的规定给摩擦系数计算带来了极大的不方便1944年莫迪（Moody）根 据这些计算式画出了摩擦系数曲线图，称为莫迪（Moody ）图，如图1莫迪图 上将不同的流动条件分成了四个区域：层流流动区、临界区（其值不确定，因为 流动可能是层流也可能是紊流）、过渡区（摩擦系数同时受到了雷诺数和相对粗 糙度的影响）和完全湍流区（充分粗糙流动，摩擦系数与雷诺数无关，仅由相对 粗糙度确定，阻力损失与速度的二次方成正比，又称为阻力平方区）莫迪图中，二OTJEJe.&T>uLJ_0010.0000.0060 to0.09O.OS0.070.060.05*TiiriarCnticai |^. low :f£呼 TransitionZonetii-44 ConlpiF?Q.Q40€30.0250.020.0151,1 , t 1 1 1 -』1\ H ...2 3 J 6^ 810* 2 3 4 56 2 3 456 2 3xlCr5 - - x1『- + - xtO6 - + - xlO^Reynolds nutnbGf Ft心二弓0.050.G40 030.020.0150.01O.OOB0.0060.0040.WJ20.4JC10.0C080QOT6O.(XX}40.00020.00010.OOO.O50 000,01「4 5$ B-irja过渡区与完全湍流区没有明确的界限，虚线用于区分过渡区与完全湍流区。

图中， 右侧纵坐标是相对粗糙度，其值与曲线对应而不是与网格对应过渡区最下面的 曲线代表水力光滑区的流动曲线，图中许多曲线在低雷诺数时与光滑管流动曲线 重合莫迪图确定的摩擦系数值与实验误差一般不超过5%图 1 MOODY 图问题提出：上面介绍的3 类问题中，有第二类和第三类问题需要用到数值分析中介绍的 迭代法求解下面就第二类和第三类问题背景各提出一种问题第二类问题：用直径为 50mm 的镀锌管从供水总管上接水，总管上压力为 450kPa,支管长度为40m，管出口水龙头位置比总管高1.2m，计算水龙头全开 时的流量第三类问题：灭火过程中用镀锌管接水，输送水流量为200L/S,如果管道长 度为35m，管内水头总损失不超过50m，试确定最小的管径解决问题：第二类问题的解决过程如下：设水温为20°C,水的运动粘度为1.003X10-6m2/s,密度为998.2kg/m3，支管 上总水头损失为p p 450 x 103h 二(厶 + z )—(红 + z )二( - + 0) - (0 +1.2)二 44.8mf p 1 p 2 998.2 x 9.807gg镀锌管当量粗糙度为0.15mm，相对粗糙度为0.003查莫迪图得，摩擦系数九最小值为0.0265计算流速 v， v — 2gdh = 2-9・81X °.°5 -必8 = 41.46,v = 6.44m/s2 九 l 0.0265 x 406.44 x 0.051.003 x 10-6=3.21 x 107再查莫迪图，流动处于充分粗糙区，九二九=0.0265，九值与初始假设一max兀 兀致，则管内流量q = vA = d2v = x 0.052 x 6.44 = 0.0126m3 /s。

v 4 4第三类问题解答过程如下：镀锌管粗糙度e=0.15mm，e/d （相对粗糙度）=0.00015/d4qv = v兀d 24 x 0.2nd 2=0.254 d 2Re _ vd _ (0.254； d2)d ^^e ———V 1.003x10—6—0.2538 x106；d、35 (0.2546 d2)2/V x d 2gVd(m)e/dRe查莫迪图九结论0.030.14730.001021.723X1060.02以九=0.02计算0.020.13580.00111.87X1060.02一致d — 0.298V 15迭代计算计算得到最小管径为 0.1358m分析问题：在数值计算方法中介绍的迭代法有简单迭代法、牛顿迭代法及其简单变形后 的方法数学上的这些方法是从理论上更深入的分析和解决问题的方法，虽然数 学上给出的方法来源与实际问题，但是它是高于实际问题的在本文给出的两个问题中，都需要进行迭代计算，这两类问题都经过这样一 个过程：先确定摩擦系数、比较摩擦系数、重新选取摩擦系数再重复计算，直到 摩擦系数基本一致这两个问题，一个是求解确定管道流量，一个是求解确定管 径，但是从数学方法上来说，本文用的迭代方法只是简单的迭代法，即先给出一 种数据，得另一个数据。

文中所用的迭代法和数值分析方法中给出的迭代法看似 并不是很相同，但是它们的本质是一样的另外一点，本文给出的计算例子没有提及精度问题，而我们在数值分析中所 接触到的迭代法需要确定精度，这种区别的原因在于：本文给出的例子是工程实 例，本身具有实际意义，即里面的数据本身具有一定的限制条件，而数值分析方 法中给出的理论分析迭代方法是纯理论的，和实际工程是具有比较大的差异的。