IG-541混合气体灭火系统设计理论和基本计算方法教学内容

1、Good is good, but better carries it.精益求精，善益求善。IG-541混合气体灭火系统设计理论和基本计算方法-IG-541混合气体灭火系统设计理论和基本计算方法一概述IG-541混合气体灭火系统作为一种新型洁净气体灭火系统，由于它兼备有效灭火、绿色环保以及对人体无伤害等特性，目前已在国内外消防领域得到广泛应用。然而，人们在大量应用它的同时，对系统性质、性能、原理等方面的量化研究却是十分不足的。国内至今尚无完整的系统设计规范，尤其缺乏完整的系统设计计算理论和方法，甚至于连基本的单元计算方法也不齐全，现有的一些计算公式基本上照搬了国外的书本，并且缺乏完整性和系统性

2、。这种理论研究远远落后于实际应用的反常现象是消防工程界特有的，也是消防系统建设与使用远远相脱节这一客观情况所造成的。国外公司虽有系统设计软件可以代客计算，但并不提供计算方法，我们只能是知其然而不知其所以然。为了解决我国已有IG-541灭火系统的设备和大量实际应用，却还没有设计计算方法的突出矛盾，确保IG-541灭火系统设计的科学先进性、安全可靠性和经济合理性，达到优化设计的目的，我们在努力学习和吸收国外先进技术的同时，还必须建立自己的理论研究体系和设计计算方法。本文探讨了IG-541气体灭火系统设计计算的理论依据，在此基础上推导了和建立了IG-541灭火系统的基本计算方法，为科学地建立具有自主

3、知识产权的IG-541灭火系统计算机设计软件奠定了基础。二系统特征IG-541灭火系统和其他固定气体灭火系统比较既有共性又具有鲜明的个性。IG-541在储存条件下呈气态，比其他灭火系统需要更大的储存容积;在高压下储存和运行，管道的承压能力要求亦较高，设备投资费用大，精确计算和优化设计可以带来明显的经济效益。IG-541灭火的有效浓度为37.5%而对人体安全的浓度为23秒及40秒，并且又要求60秒钟内达到灭火浓度。这也是一个相当严格的的设计约束条件。IG-541灭火系统和其他灭火系统相比,灭火剂设计浓度以及喷射速率的容差范围小得多，且与平常容易发生的误解不同,宽裕的设计不仅浪费投资，设计结果也未

4、必安全。因此，系统设计应采用精确的、全过程动态模拟的分时计算方法。IG-541设计计算的有利条件是：物系临界温度低，整个过程在单一气相下发生，可以通过严格的方法，借助电子计算机进行精确的计算。IG-541灭火系统设计的主要目标是要保证在装置启动后的指定时间内，防护区中的灭火剂达到设计浓度，其中计算IG-541气体在系统各单元中的流动推动力和阻力是关键，二者又取决于系统的物性和单元的设备特征，兹在下文逐一讨论。三纯组份性质物质的物理化学性质甚多,这里只讨论和IG-541系统设计有关的P-V-T性质、热力学性质和迁移性质。1基础物性氮气、氩气和二氧化碳皆为常见气体，其有关性质可以从手册中查到。兹参

5、照ChemicalPropertiesHandbook（1999）一书，将相关数据罗列如下。其他资料上除分子量以外的数据并不完全相同，但对本过程的设计计算结果影响甚微。表一.基础物性和热力学性质名称符号单位ArCO2N2英文名ArgonCarbondioxideNitrogen分子量Mg/mol39.94844.01028.013常冰点TfK83.80216.5863.15常沸点TbK87.28194.7077.35临界温度TcK150.86304.19126.10临界压力PcBar48.9873.8233.94临界比容Vccm3/mol74.694.090.1临界密度cG/cm30.5356

6、0.46820.3109临界压缩因子Zc无因次0.2910.2740.292偏心因子无因次0.0000.2280.040偶极距Debye0002热性质流体在理想气体状态下的热性质是计算热力学性质的基础。通常表达为热容或焓的多项式。如:hi0=C0,i+C1,iT+C2,iT2+C3,iT3+（式1）这些多项式系数通常是用于相当宽的温度范围，而IG541的工作温度范围较窄，约在200至320K之间，可以专门回归较为简短和精确的多项式。兹将APIProject44SelectedVolumeofPropertiesofChemicalCompounds所列CO2、N2的恒压热容CP0文献值列于表二

7、。Ar的CP0受温度影响极小，用ThePropertiesofGasesandLiquids3rdEd一书附录中的多项式求得。表二.理想气体状态下的恒压热容Cal/mol/K温度，K150200273.15298.15300400Ar4.96814.96804.96784.96784.96784.9679CO27.2287.7338.5948.8748.8949.876N26.9566.9576.9596.9616.9616.9913低压下的气相粘度在低压下气体的粘度和压力关系不大，可以视为仅仅是温度的函数。经与HandbookofChemistryandPhysics,80thEd（1999

8、-2000），ChemicalPropertiesHandbook（1999）的资料上的数据比较，发现用多项式拟合低压气体的粘度的精确度不高。而用Lennard-Jones12-6分子势能位计算，25C下的误差降至0.1%左右。下表中Ar的数据取自美国石油学会API手册,其余数据取自ThePropertiesofGasesandLiquids,后者Ar的参数计算结果误差较大。表三.迁移性质-气相粘度L-J分子势能位ArCO2N2势能参数/k，K124.9195.271.4碰撞半径，3.42337.1329.124Lennard-Jones12-6分子势能位计算低压气体的粘度的公式用:gas=2

9、6.69(MT)1/2/2/V（式2）式中:gas气体粘度,PV=1.16145/T*0.14874+0.52487/exp(0.7732T*)+2.16178/exp(2.43787T*)T*=T/(/k)四混合物性质IG-541混合气体的配方是公开的，即52%（mol）的氮气、40%（mol）的氩气和8%（mol）的二氧化碳气体。混合物物性的详尽实验数据很少。混合物的性质主要通过“混合规则”计算求得，而这些规则则是经过若干离散的实验点来验证的。令混合物中i组份的含量为Xi分子分率。下标i,j,k为组份序号，下标m表示混合物。则有：1分子量Mm=XiMi（式3）2混合物临界参数混合物的临界参

10、数并不是混合物的真临界性质，而是用于计算混合物P-V-T性质和热力学函数用的参数。采用Lee-Kesler状态方程的混合规则：Vci=ZciRTci/Pci（式4）Zci=0.29050.085i（式5）Vcm=0.125XjXk(Vcj1/3+Vck1/3)（式6）Tcm=0.125XjXk(Vcj1/3+Vck1/3)3(Tcj+Tck)2/Vcm（式7）m=xii（式8）Pcm=(0.29050.085m)RTcm/Vcm（式9）3IG-541在理想气体状态下的热性质系数理想气体状态下的热性质系数符合加和性，从表二的数据经过计算和单位换算得到IG541在理想气体状态下焓的热性质系数为：C

11、0=8.8(此项系数用来校准150400K的焓值,计算偏差0.03%)C1=25.648C2=0.0005042C3=1.715E-6C4=81.04(此项系数取,1psia压力下的理想气体单质为零熵)用以上热性质系数计算IG541在理想气体状态下的热容、焓和熵的公式是：Cp0=C1+2C2T+3C3T2（式10）h0=C0+C1T+C2T2+C3T3（式11）s0=C1lnT+2C2T+3C3T2/2+C4-RlnP（式12）上述公式中的单位是:mol,K,J,Pa,上角标0表示理想气体状态。4气体混合物的粘度气体混合物的粘度用美国石油学会APIproject44推荐的方法计算：（式13）其

12、中ij为充间参数（式14）将以上公式求得的IG-541粘度再回归成温度的多项式：m0=178.86+0.5123-0.00039（式15）（式15）中m0是IG-541在低压下的粘度，摄氏温度。在-10至50范围内回归误差小于0.1%。5低压下的气相导热系数纯组分的导热系数由LangsChemicalHandbook15thEd（1999）查得,如下表所示:表四.气相导热系数J/s/m/K温度,02040Ar0.01660.01760.0186CO20.01440.01600.0176N20.02410.02560.0270IG-5410.02020.02160.0228气体混合物的气相导热系

13、数按Wassiljewa方程计算:m=(yii/yjAij)（式16）式中:m气相混合物的导热系数I组分i的气相导热系数Aij组分之间的充间作用参数,用LindsayandBromley方法Aij=0.251+(i/j)(Mj/Mi)0.75(T+Si)/(T+Sj)0.52(T+SiJ)/(T+Si)（式16-1）其中:Sutherland常数Si=1.5Tbi;Sij=Cs(SiSj)0.5非极性气体Cs=1用MasonandSaxena方法以及忽略充间作用参数取得的结果和上述方法偏差不大于1%。按以上方法求得IG-541在0、20、40低压下的气相导热系数，回归成多项式：0=0.0202

14、38+6.75910-5t6.46810-8t2（式17）式中:0IG-541在低压下的气相导热系数,J/s/m/Kt摄氏温度为便于比较和利用现将IG541混合气体物理性质及相关计算公式汇总如下：IG541混合气体的物理性质ThePROPERTIESofIG541MIXTRUE名称符号单位氩气二氧化碳氮气IG541英文名SymbolunitARGONCarbonDioxideNitrogenInergen分子式ArCO2N2分子量MWg/mol39.94844.01028.01334.067常沸点TbK87.28194.777.35临界温度TcK150.86304.19126.1147.7临界

15、压力PcBar48.9873.8233.9442.39临界比容Vccc/mol74.694.090.183.2临界压缩因子Zc0.2910.2740.2920.287偏心因子OM0.0000.2230.0400.039碰撞半径A3.4233.9413.798势能参数/KK124.9195.271.4恒压热容Cp,25J/mol/K20.78537.12929.12526.429恒容热容Cv,25J/mol/K12.47028.81420.81018.115绝热指数,251.66671.2891.4001.459低压下气体粘度,25P224.42150.5175.52191.42偶极距Debye

16、0000IG541分子分率YMol%40852100注:正体数值是文献值,斜体数值是计算值。计算方法和依据另详。IG541的密度:标准状态(0,1atm)1.521kg/m3,对空气比重1.176存储状态(20,15Mpa)233.6kg/m3IG541在低压下的粘度多项式:=178.86+0.5123t-0.00039t2为绝对粘度,P;t为摄氏温度。上式适用于-10至50。IG541在低压下的导热系数多项式:=202.38+0.6759t6.46810-4t2为导热系数,W/cm/K;t为摄氏温度。IG541在理想气体状态下的焓多项式:H0=25.648T+5.04210-4T2+1.71

17、510-6T3上式中:H0为焓,J/mol;T为绝对温度,T=t+273.15五压力下的性质1P-V-T关系表征流体压力、比容和温度（P-V-T）关系的方程叫着状态方程。最简单,也是最古老的状态方程是理想气体状态方程:PV=RT式中：P-压力PaV-比容m3/molR-气体通用常数=8.31441T-温度K理想气体状态方程对常压常温或高温的“永久气体”可以使用，但是用来计算低温高压气体，或是非永久气体误差就比较大，对于液体则完全不适用，于是，近百年来许多科学家提出了对理想气体定律的修正，这类方程式称为真实气体状态方程。真实气体状态方程的通式是：PV=zRT式中z的叫做压缩因子，表示真实气体和理

18、想气体的差异。文献上发表的真实气体状态方程有二百多个。其中最准确的是Lee-Kesler方程，它的优点体现在：1.在很宽的温度压力范围内有良好的准确度；2.既能用于气相也能用于液相；3.对于任何非极性流体，只需要临界温度、临界压力和偏心因子三个参数便可以求解，而这三参数都是宏观可测量的物性；4.对于混合物利用混合规则，可以像纯物质一样，用同样的方程式求解；5.对于混合物，不仅能求解总性质，还能解得偏性质，即混合物中各个组份的性质。在Lee-Kesler方程发表前后，还有一些优秀的状态方程。这些方程求解比较简单，在一定的场合可以满足工程计算需要，但就准确度而言，还是及不上Lee-Kesler方程

19、。状态方程不仅可以用来求解流体的密度，还可以和物质的热性质系数结合起来，求得在真实流体（液体和压力下的气体）状态下的热力学性质，传热性质等等，用于流体力学、化工热力学、工程热力学的计算。对于灭火系统而言，不仅管道输送，而且孔板节流、喷嘴喷射、气体膨胀和压缩的精确计算度需要用到状态方程，因此我们选用了Lee-Kesler方程设计计算IG541混合气体灭火系统。Lee-Kesler状态方程为多参数方程，其形式是：z=Pr/Vr/Tr=1+B/Vr+C/Vr2+D/Vr5+c4（+/Vr2）exp（-/Vr2）/Tr5/Vr2（式18）其中：B=b1b2/Trb3/Tr2b4/Tr3（式18-1）C

20、=c1c2/Tr+c3/Tr2（式18-2）D=d1+d2/Tr（式18-3）Vr理想对比比容=Pc*V/n/R/Tc，其中V/n称为比容Pr对比压力=P/PcTr对比温度=T/TrLee-Kesler状态方程是三参数普遍化方程，只要知道物质的临界压力、临界温度和偏心因子就能求解。但是方程中的压缩因子z是隐函数，要用试差法求解。表四.Lee-Kesler状态方程系数系数简单流体参比流体系数简单流体参比流体b10.11811930.2026579c30.00.016901b20.2657280.331511c40.0427240.041577b30.1547900.027655d11040.15

21、54880.48736b40.0303230.203488d21040.6236890.0740336c10.02367440.03133850.653921.226c20.01869840.05036180.0601670.03754Lee-Kesler状态方程有多达七个实根，其中只有一个或两个实根有物理意义，要用试差法求解。本文作者曾对此作过研究，具体方法可参考烃类物理化学手册一书。IG541气相压缩因子CompressibilityfactorofIG541(Z)温度-1001020304050压力MPaPrTr1.7821.8491.9171.9852.0522.1202.18810.

22、2360.98960.99120.99260.99380.99490.99580.996620.4720.97960.98280.98570.98810.99020.99200.993630.7080.97010.97500.97920.98280.98590.98870.991040.9440.96110.96760.97320.97800.98210.98570.988951.1790.95270.96080.96780.97370.97880.98320.987161.4150.94500.95470.96290.96990.97600.98120.985771.6510.93800.9

23、4920.95860.96670.97360.97960.984881.8870.93180.94430.95490.96400.97180.97850.984392.1230.92640.94020.95190.96190.97050.97790.9843102.3590.92190.93690.94960.96040.96970.97770.9847112.5950.91830.93430.94790.95950.96950.97810.9855122.8310.91560.93260.94700.95920.96980.97890.9868133.0670.91390.93160.946

24、70.95960.97060.98020.9884143.3020.91310.93150.94710.96050.97200.98190.9905153.5380.91330.93210.94820.96200.97390.98410.9930163.7740.91440.93360.95000.96410.97630.98680.9959Lee-Kesler状态方程要对三个流体求解：一个是需要计算的流体，即IG541，公式中用上角标（i）表示；一个是偏心因子为0的“简单流体”，上角标用（0）表示；还有一个是偏心因子为0.3978的正辛烷，叫做“参比流体”，上角标（r）。先在对象流体的操作条

25、件求得对比压力Pr和对比温度TrPr=P/Pcm（式19）Tr=T/Tcm（式20）从（式18）分别求得简单流体和参比流体的Z（0）和Z（r），然后用以下公式计算对象流体的压缩因子Z（i）=Z（0）+Z（r）-Z（0）（i）/（r）（式21）（式20）中（i）为对象流体即IG541偏心因子m，（r）为参比流体的偏心因子=0.3978。（式21）的物理意义是：流体的压缩因子等于同对比温度压力下简单流体的压缩因子，加上和偏心因子相关的修正值。偏心因子在微观上反映物质分子的大小和形状，小的球形分子如氩气，偏心因子等于0；宏观上偏心因子通过物质的对比饱和蒸汽压来计算：=-lg(PSTr=0.7/Pc)

26、1（式21-1）（式21-1）中lg是十进对数，PSTr=0.7是对比温度为0.7时的饱和蒸汽压。据此，作者曾在70年代进行过验算，由表四的系数，通过Lee-Kesler方程计算出“简单流体”和“参比流体”气液相逸度，求得对比温度0.7下的相平衡和饱和蒸汽压,按照（式21-1）定义计算简单流体和参比流体的偏心因子，分别为：简单流体（0）=0.00529（式21-2）参比流体（r）=0.39547（式21-3）因此，作者将（式21）改成：Z（i）=Z（0）+（i）-（0）/（r）-（0）Z（r）-Z（0）（式21-4）以下是采用上述方法进行的氩气、二氧化碳和氮气的压缩因子验算。此项工作的目的一是

27、验算Lee-Kesler方程对于IG541组份的准确性；二是比较（式21）、（式21-4）何者更符合实际情况。一般在近临界区或高压下P-V-T计算误差大，为此本文选择IG-541实际运行范围内的温度和较高的压力。验算结果列于下列表格。表五.氩气压缩因子验算操作条件文献值（式21）（式21）（式21-4）（式21-4）压力温度比容压缩因子比容误差比容误差BarKcc/g无因次cc/g%cc/g%1002002.96250.71172.9616-0.03072.9571-061920.88774.63450.33144.62970.22691003005.96330.95

28、515.98150.34095.97640.21971003507.19630.98797.20910.17757.20380.10411502001.81660.65461.82510.46971.82190.28981502502.98340.96012.99790.48702.99350.33831503003.94630.94803.96730.53193.96260.41181503504.81430.91134.82380.38484.82790.2820平均误差%0.33970.2569表六.二氧化碳压缩因子验算操作条件文献值（式21）（式21）（式21-4）（式21-4）温度压

29、力比容压缩因子比容误差比容误差KBarcc/mol无因次cc/mol%cc/mol%2205.9962744.80.89972772.71.01522772.91.02362308.9351882.80.87971883.30.02881883.50.003924012.8301318.80.84791315.8-0.22981315.9-0.216825017.856941.950.8092938.71-0.3437938.86-0.327926024.194682.560.7639679.37-0.4668679.50-0.447927032.034498.020.7107495.05-0

30、.5973495.16-0.574728041.595361.840.6465359.10-0.7564359.20-0.729429053.152256.250.5649253.37-1.1223253.46-1.089530067.095164.000.4411158.60-3.2955158.67-3.2500平均误差%0.87290.8555表七.氮气压缩因子验算操作条件文献值（式21）（式21）（式21-4）（式21-4）压力温度压缩因子压缩因子误差压缩因子误差BarK无因次无因次%无因次%1002000.844990.845010.00190.84392-0.12721002500

31、.958560.960200.17130.959250.07241003001.00511.006210.11001.005380.02751003501.02691.027120.02191.02639-0.05001502000.849860.850190.03920.84897-0.10471502500.971340.973880.26170.97265002561.028240.25751.02712005111.052990.17961.051980.0835平均误差%0.13040.0936表五、六、七的第列,是国际纯和应用化学

32、学会(IUPAC)资料上提供的数据;表五、六的第列数值，由第列换算而来。这些数据表明：除了二氧化碳外，Lee-Kesler方程的计算结果皆与文献值十分符合；而本文作者建议的偏心因子校正式（式21-4），比原作者提出的（式21）更准确。表六所列的温度范围，在二氧化碳的常冰点和临界点之间，压力为相应温度下的饱和蒸汽压，无论对实验还是计算都是高误差区。二氧化碳还有一个特殊情况：在对比温度0.7时已经是固相，（式21-1）的引用发生了问题，因而不同文献上的偏心因子值有差异。尽管物质固相升华压和液相的蒸汽压函数连续，导数却不连续。相比较之下用液相的蒸汽压外推计算偏心因子要合理一点。经验证，TheProp

33、ertiesofGasesandLiquids3rdEd一书上的Harlacher蒸汽压方程,从二氧化碳的冰点到临界点度很准确，由此算得的偏心因子为0.223。2热力学差值函数热力学差值函数的定义是：真实流体的热力学性质和同温度下理想气体的热力学函数之差。热力学差值函数用无因次数的形式表示，通式为：g（i）=g（0）+g（r）-g（0）（i）/（r）（式22）和g相应的恒压热容、熵、焓的差值函数为：（Cp-Cp0）/R，（ss0）/R，（hh0）/RTc。由（式18）求得压缩因子代入Lee-Kesler状态方程的差值函数表达式，求得g（r）、g（0），经（式21）求得g（i），然后结合（式10

34、）至（式12）即可求得真实流体在指定温度压力下的热力学函数。例如h=RTc（hh0）/RTc+h0=RTc（hh0）/RTc+C1T+C2T2+C3T3+C4T4（式23）IG541在不同压力下的焓EnthalpyatpressureofIG541(J/mol)理想气体:温度-1001020304050ARGONi.g.5472.05679.95887.76095.66303.46511.36719.1CO2i.g.8093.68451.18813.09179.39550.19925.410305.4N2i.g.7651.17942.08233.08524.28815.49106.89398.

35、3IG541*i.g.6815.47078.37341.67605.27869.38133.88398.7压力下气体：温度-1001020304050压力MPaPrTr1.7821.8491.9171.9852.0522.1202.18810.2366714.66984.77254.47524.07793.58062.98332.320.4726613.36890.97167.47443.17718.17992.68266.630.7086511.76797.27080.77362.77643.47922.98201.740.9446410.26703.86994.57282.97569.37

36、854.28137.651.1796308.96610.96909.07203.97496.27786.38074.561.4156208.16518.76824.37125.97424.17719.58012.571.6516108.26427.56740.77049.07353.17653.87951.681.8876009.56337.66658.46973.47283.47589.47892.092.1235912.56249.36577.76899.27215.27526.47833.7102.3595817.46162.86498.66826.87148.57464.97776.8

37、112.5955724.76078.46421.66756.17083.57404.97721.4122.8315634.75996.56346.76687.47020.37346.77667.6133.0675547.85917.26274.16620.86959.07290.27615.4143.3025464.25840.86204.06556.46899.77235.57564.8153.5385384.25767.36136.66494.36842.57182.67516.0163.7745307.95697.16071.96434.66787.47131.77468.9注：*=25

38、.648*T+0.0005042*T2+1.175e-6*T3IG541在不同压力下的内能InternalEnergypressureofIG541J/mol理想气体:温度-1001020304050ARGONi.g.3284.13408.83533.53658.23782.93907.64032.3CO2i.g.5905.66180.16458.86742.07029.67321.87618.6N2i.g.5463.15670.95878.86086.86294.96503.16711.5IG541i.g.4627.54807.24987.45167.95348.85530.15711.9压

39、力下气体：温度-1001020304050压力MPaPrTr1.7821.8491.9171.9852.0522.1202.18810.2364549.44733.64917.65101.75285.85470.15654.620.4724470.04658.84847.05034.85222.35409.75596.930.7084389.34583.04775.64967.35158.35348.85538.940.9444307.34506.24703.44899.25093.95287.75480.751.1794224.34428.74630.74830.75029.15226.35

40、422.461.4154140.44350.54557.54761.94964.25164.85364.071.6514055.84271.94484.04692.84899.15103.25305.681.8873970.84192.94410.34623.74834.05041.75247.292.1233885.54113.94336.64554.74769.04980.35189.1102.3593800.44035.04263.14485.84704.34919.25131.1112.5953715.63956.54190.04417.44639.94858.45073.6122.8

41、313631.53878.64117.34349.44575.94798.05016.4133.0673548.33801.54045.44282.04512.54738.14959.6143.3023466.43725.43974.44215.34449.84678.94903.5153.5383386.03650.43904.34149.64387.84620.34847.9163.7743307.23576.93835.44084.84326.74562.44793.0IG541在不同压力下的熵EntropyatpressureofIG541(J/mol/K)温度-10010203040

42、50压力MPaPrTr1.7821.8491.9171.9852.0522.1202.1881psia.00158150.93151.91152.85153.77154.65155.51156.340.10.024128.69129.67130.62131.53132.42133.28134.110.20.047122.93123.91124.85125.77126.65127.51128.350.30.071119.56120.54121.48122.40123.28124.14124.970.40.094117.16118.14119.09120.01120.89121.75122.580

43、.50.118115.31116.29117.24118.15119.04119.89120.730.60.142113.50114.50115.47116.41117.31118.19119.030.70.165112.21113.22114.19115.13116.03116.91117.750.80.189111.10112.11113.08114.02114.92115.80116.640.90.212110.13111.13112.10113.04113.94114.82115.6610.236109.25110.26111.23112.16113.07113.94114.7920.

44、472103.19104.22105.22106.17107.10107.99108.8530.70899.51100.58101.60102.58103.52104.43105.3040.94496.8297.9198.9699.96100.92101.84102.7451.17994.6695.7896.8697.8898.8699.80100.7161.41592.8493.9995.0996.1497.1498.1099.0271.65191.2592.4493.5794.6495.6696.6397.5781.88789.8491.0692.2193.3194.3595.3496.2

45、992.12388.5689.8190.9992.1193.1794.1895.15102.35987.3988.6889.8891.0292.1093.1394.11112.59586.3187.6388.8690.0291.1292.1693.16122.83185.3086.6587.9189.0990.2191.2792.28133.06784.3685.7487.0288.2389.3690.4491.46143.30283.4884.8986.1987.4188.5789.6690.69153.53882.6584.0885.4186.6587.8288.9289.97163.77

46、481.8783.3284.6785.9387.1188.2389.29注：表列数据,参照美国API采用1K，1psia,气体状态下的单质为零熵基准,并且考虑到IG541气体的混合熵效应。psia=每平方英寸磅,绝对压力。3压力下的粘度气体粘度随压力增高而增大，气体在压力下的粘度由下式估算：（m-m0）m=1.08exp(1.439Prm)-exp(-1.111rm1.856)（式24）式中：m高压下气体的粘度，pm0低压下气体的粘度，Prm虚拟对比密度=m/cmm为操作条件下混合气体的密度=n/Vcm为混合气体的临界密度m=Tcm1/6/（Mm1/2Pcm2/3）为物性参数由Lee-Kesl

47、er状态方程求得的Vr计算出m，并由（式15）求得IG-541在低压下的粘度，结合（式24）和IG-541的物性常数即可求得压力下的粘度。IG541在不同压力下的粘度ViscosityatpressureofIG541(P)=0。03245低压气体:温度-1001020304050ARGONi.g.202.17208.53214.80220.96227.04233.02238.92CO2i.g.134.31139.15143.92148.64153.29157.89162.43N2i.g.160.91165.36169.74174.04178.27182.44186.54IG541i.g.17

48、3.68178.85183.94188.95193.88198.73203.51IG541*i.g.173.70178.86183.94188.95193.88198.73203.50注：i.g.=理想气体IdeaGas压力下气体：温度-1001020304050压力MPaPrTr1.7821.8491.9171.9852.0522.1202.18810.236175.68180.76185.77190.71195.57200.36205.0820.472177.95182.93187.84192.69197.47202.19206.8430.708180.52185.36190.15194.

49、88199.57204.19208.7640.944183.39188.06192.69197.30201.86206.38210.8551.179186.57191.03195.48199.92204.35208.74213.0961.415190.06194.27198.51202.76207.02211.27215.5071.651193.86197.78201.77205.81209.89213.97218.0581.887197.98201.55205.27209.07212.94216.84220.7692.123202.39205.59208.99212.52216.16219.

50、86223.61102.359207.10209.88212.93216.17219.56223.04226.59112.595212.09214.41217.08220.00223.11226.36229.71122.831217.35219.17221.43224.01226.83229.82232.96133.067222.86224.14225.96228.17230.68233.42236.32143.302228.59229.31230.67232.50234.68237.13239.79153.538234.53234.66235.53236.96238.80240.96243.

51、36163.774240.65240.18240.55241.55243.04244.90247.04注：*=178.86+.5123*t-.00039*t*t4压力下的导热系数导热系数随压力增高而增大，气体在压力下的导热系数由下式估算：rm0.5（m-m0）Zc5=14.010-8exp(0.535rm)-10.5rm2（m-m0）Zc5=13.110-8exp(0.67rm)-1.069（式24）式中：m高压下气体的导热系数，J/s/m/Km0低压下气体的导热系数，J/s/m/Krm虚拟对比密度=m/cm=Tcm1/6Mm1/2/Pcm2/3为物性参数,其中Pcm的单位是大气压。IG541

52、在不同压力下的导热系数ThermalConductivityatpressureofIG541W/cm/KZc5=0.00227温度-1001020304050压力MPaPrTr1.7821.8491.9171.9852.0522.1202.18800.000195.56202.38209.07215.64222.07228.38234.5610.236200.91207.53214.03220.42226.69232.84238.8820.472206.48212.86219.15225.34231.43237.42243.3030.708212.29218.40224.44230.4123

53、6.30242.10247.8140.944218.32224.13229.90235.62241.29246.89252.4251.179224.59230.05235.52240.97246.40251.79257.1161.415231.08236.16241.29246.46251.63256.78261.9071.651237.80242.45247.22252.07256.96261.86266.7681.887244.73248.92253.30257.80262.40267.03271.7092.123251.87255.55259.51263.65267.93272.2927

54、6.70102.359259.18262.34265.84269.60273.54277.61281.77112.595266.67269.26272.29275.65279.24283.01286.90122.831274.29276.30278.84281.77285.01288.46292.07133.067284.24284.24285.47287.97290.83293.96297.29143.302294.15293.32293.49294.23296.71299.51302.54153.538304.21302.53301.99302.32303.29305.08307.8316

55、3.774314.36311.84310.59310.31310.77311.81313.13注：*=202.38+.6759*t-6.468E-04*t2六.管道输送压力降IG-541灭火系统输气属于“可压缩流体输送”，其输送阻力计算方式分为：图表、经验公式、理论公式三类。图表和经验公式很多，大多针对特定的系统，并隐含了管壁粗糙度等影响输送阻力的因素，比较简便，适合手工计算。一般说来较新的方法求得的压力降较小，这是因为技术进步以及近代管道的内壁比较为光滑的缘故。从理论推导的公式严谨，适合各种参数变化的场合，但是求解复杂，而且文献上介绍的公式并不完全一致。本课题参考了多份资料，并经过理论推导，

56、确认以下公式适合于水平管道的压力降计算。P12-P22=ZRmTG2fL/D+2ln（P1/P2）（式25）式中：P压力，PaZ气体压缩因子，无因次Rm用重量表达的气体通用常数=1000R/MmT绝对温度，KG重量流率，kg/m2/sf摩擦系数，无因次L从1到2的管段长度，mD管内径，m下标：1管段入口2管段出口（式25）和“EncyclopediaofChemicalProcessingandDesign”Vol22，p241（1985）上的公式一致，只是摩擦系数的单位不一样。式中fL/D部分为摩擦阻力项，2ln（P1/P2）部分为气体动能项。在长距离输送时气体动能项通常予以忽略，可是灭火系

57、统的管道短，气体动能的变化仍需考虑。流体流动的摩擦系数f是雷诺数Re和管壁相对粗糙度的函数,按雷诺数的范围分段确定。IG-541的流动条件在Re4000的紊流区，摩擦系数适用Colebrook式1/f0.5=-2loge/（3.7D）+2.512/Re/f0.5=-0.868589lne/（3.7D）+2.512/Re/f0.5（式26）式中：e管内壁绝对粗糙度，mD管内径，m；e/D称为相对粗糙度Colebrook式中的f为隐函数，求解不便，因此后来有许多学者提出了各种显函数式，经验算以下公式与Colebrook式较为符合：f=-1.325/ln（e/3.7/D）+5.75/Re0.92（式

58、27）Colebrook式经过很多人的检验，已经被公认。因而在计算时采用（式27）求得f初值，然后代入（式26）计算。气体的雷诺数表达为Re=Du/（式28）式中：u气体流动的线速度，m/s气体密度，kg/m3气体粘度，Pas=m/107由于所研究对象是开口系统，且与大气相通，所以如果管道不是水平安装，或者不在设计基准水平面上，则还要加上位能修正项：P21=gc（Z1-Z2）（m-AIR）（式29）式中：P21由位能和空气浮力产生的附加压差，Pagc重力加速度=9.80665m/s2Z1,Z2管段进出口标高，mAIR环境大气密度kg/m3=Pu/RAIR/TRAIR=1000R/28.964其

59、中Pu是防护区环境压力,不是管内压力。该项对于管内压力损失的影响是微乎其微的。七管道热平衡流动系统的能量平衡通式是:(h2-h1)+(u22-u12)/2+gc(H2-H1)=qws（式30）就IG541灭火系统的具体情况而言,轴功ws=0，于是（式30）简化成:(h2-h1)+(u22-u12)/2+g(H2-H1)=q（式31）在等温条件下按求得（式25）求得出口压力即可算出出口流体的h2，u2,再由（式31）计算出系统从环境吸收的热量qJ/kg。消防管道的实际情况是既非等温又非绝热，所以实际的出口状态参数界于等温与绝热之间。消防管道并不保温，其外部传热受风速限制，按王绍周所著管道工程设计

60、施工与维修一书推荐，无风的情况下大气与管壁的传热系数为4250kJ/m2/h，经单位换算和乘以圆周率之后为36.6543.63取下限，便可得到环境向管壁的传热量为：Q1=36.65DOL（TuTw）（式32）式中：Q1环境向管壁的传热量，J/sDO管外径，mL管段长度，mTw管外壁温度，KTu环境温度，KIG-541的喷射属于不稳定过程，起初管壁的温度和环境相等，管内流体首先向管壁吸热,致使管壁温度下降。钢管降温放热量与管壁温度降之间的关系为：Q2=449（DO2D2）/4L7800Tw=275019（DO2D2）LTw（式33）式中：Q2管壁向流体的传热量，JD管内径，mTw管壁温度的变化，