二氯乙烷裂解炉全周期模拟炉膛炉管藕合数值迭代方法与第三章清洁管时期迭代方法类似，不同之处在于在外壁增加了一个焦层计算的循环迭代步骤，以及全周期模拟结束的判断条件。详细算法框架见图4.1。由于工业二氯乙烷裂解炉的一个完整的运行周期较长，本章中结焦过程的采样时间间隔为一周。整个计算由两层循环嵌套而成，第一层与清洁管数值模拟方法一致，通过联合求解炉膛炉管模型，迭代求解烟气温度，从而满足炉膛炉管热量平衡。外层循环用于更新焦层厚度，计算炉管外壁最高温度，判断是否满足停炉清焦条件。具体步骤为:
    步骤1:设定炉膛炉管模型的操作条件、计算收敛指标及精度，包括烟气温度计算精度(10C)，及压力的计算精度(1 kPa)，同时确定好边界条件，并假设一组炉膛烟气初始温度;
    步骤2:以各区烟气温度分布及入口压力作为边界条件，求解反应管模型，得到沿管长方向的热通量分布。
    步骤3:将炉管模型计算所得沿管长方向的热通量分布作为边界条件，计算炉膛辐射室内辐射传热以及热量平衡模型，得到新的烟气温度分布。
    步骤4:将新炉膛模型计算的烟气温度与上一阶段计算的烟气温度对比，若满足精度要求，转步骤5，若未满足精度要求，则从新计算得到炉膛辐射段烟气温度分布作为边界条件，转步骤20
    步骤S:检查炉管出口压力是否满足收敛精度，若未满足条件，则根据出口压力重新估计入口压力，转步骤2，若满足条件，转步骤60
    步骤6:计算炉管外壁最高温度，判断其是否达到炉管金属材质所能承受的最高温度，若未达到极限温度，则计算新的焦层厚度，更新炉管流通内径，转步骤2，如果达到极限温度，则二氯乙烷裂解炉全周期计算结束，输出运行时间及全周期模拟关键性能指标。www.anhuanchem.com