循环流化床总体质量平衡求解算法探讨.pdfVIP

循环流化床总体质量平衡求解算法探讨 白泉 李政 倪维斗 (清华大学1o0084北京) 【摘要】 奉史结合循耶流化床的特点，设计出了用r求解循环流化床总体质量平衡的双重迭代 算法，并用牛顿法和参数牛顿法加以实现，不但简化了编程工作．而且使计算速度大幅度提高．文中 还结合循环流化床流动参数，对参数叶。顿法的参数选择进行了有益的思索和探讨。 关键词循环流化床：总体质量平衡；牛顿迭代法；参数牛顿法 在循环流化床锚炉模型化研究中，气圃两相流动状态的确定至关重要，因为它为燃烧、 传热、物质和能量平衡计算提供必需的边界条件。由于循环流化床内固体呈宽筛分分布，因 此，流动状态的确定最终归结到不同粒径颗粒的质量份额的确定。亦即总体质量平衡方程的 求解。 由于存循环流化床总体数学模型中，总体质量平衡计算需要大量、反复进行，因此尽可 能地提高其计算速度是十分必要的。为此，本文拟探讨牛顿迭代法和参数牛顿法在循环流化 床总体质量平衡计算中的席用。 1 牛顿一拉普森(Newton—Raphson)迭代算法和参数牛顿法简介 牛顿拉普森(Newton—Raphson)迭代算法(简称牛顿法)是求解非线性联立方程组的常用 方法，其表达式为：z。=以_[F’(x)一F[(Ⅳ)】，其中，J为多维求解变量组，，‘一 为求解函数。理沧上可以证明，它足“_二阶收敛”的，所以收敛速度比一般的迭代算法(如 一分法、割线法等)都要快。 参数牛顿法是在牛顿迭代法奇勺基础上在相应的迭代项处增加一个参数∞，其算法表达 ，rr、 式为：』。=x。一∞三二竺。参数牛顿法中参数的选择是根据表达式的不同而不同的，没有 … 4‘ ¨，。(x) 一个数学E的崮定的彩式。当参数选择得台适的时候，能够增加收敛的速度，得到更满意的 计尊效果；如果参数选择得不合适，很有HJ能造成相反的结果。 2循环流化床总体质量平衡简介 由前所述，循环流化床密相区的质量构成雎戒n是确定宽筛分流动特性的关键，它是通 过求解循环流化床燃烧系统总体质量平衡得到的。 在将床料分为N个尺寸组时，任意第t组灰颗粒总体质量平衡可以表达为： ％。A一凡t一一9^+毗^+△睨m^=0 (1) 式中，∥F。女为煤成灰生成的第t组灰颗粒流率，Dk和FLk分别为排渣和飞灰中第女组 灰颗粒流率．形．t为石灰石脱硫反应导致第^组灰颗粒的质量增加率，△％ff女为由于磨耗 356 引起的第t组灰颗粒的净质量变化速率。 上式中，除第一项外，其余各项均可表达为历域女)、密相床高hb以及排渣量D的函数 因此式(1)可以简单表示为： (k，j=1，2…，N) (2) F，(BPd(尼)，南6，Do“f)=O 此外，显然所有颗粒组的质量份额之和应当为1，即有 Ⅳ ∑成d(t)=1．o (3) 々=l 3总体质量平衡需要满足的边界条件 个未知量，因此方程封闭，可以得解。而hb则需要根据运行控制压降确定。记循环流化庆 炉膛实际压降为△P，运行压降控制值为△P。，则总体质量平衡需要满足的边界条件为： △P=△‰ (4) 炉膛压降△P为密相床高度以及炉膛内固体颗粒浓度的函数，可表示为： H △P=。gp，(1一s出。)^6+gp。I(1一s(向))d向 (5) 庙 den为悬浮段及密相区的颗粒空隙率，h为高度， 其中，g为重力加速度，ps为固体密度，￡和8 H为炉膛总高度。 由于方程(5)涉及到流化床内固体浓度，因此hb的计算利上节总体质量平衡的求解耦合 在一起，给整体求解及计算机编程带来困难，因此必须精心设计数值计算方法，以简化编程 难度。 4总体质量平