Thomas理论模型及其应用.pptVIP

Thomas理论模型及其应用 Thomas 理论模型 Thomas 理论模型是着眼于实际情况而考虑的一个理论模型，它既考虑体系内的温度分布，又考虑体系与环境的温度突跃。 该模型的特点是不仅体系内的温度分布随空间位置及时间的变化而变化，而且体系与环境的温度突跃也随时间的变化而变化。 体系和环境的温度变化可表示为空间坐标和时间坐标的函数。 在实际应用中，为了解题方便，一般认为Thomas 模型的温度分布具有对称性。 Thomas模型温度分布示意图 三种理论模型的内在关系 上面描述的热自燃模型称为Thomas模型，反应系统为Thomas系统。 实际上，Thomas边界条件包含了Frank-Kamenetskii边界条件和均温系统边界条件。 Semenov边界条件和Frank-Kamenetskii边界条件分别是Thomas条件的两个极限情况。 Thomas模型及其边界条件 Thomas边界条件属于Robin类边界条件，由下式给出： 无量纲化后的热平衡方程： 对于Bi＝0，表示系统内部热阻很小，温差集中在系统边界上，系统内部温度均匀，对应的是均温系统情况，称为Semenov边界条件。 当Bi ，则温差全部表现在系统内部，边界面上温度与环境温度相同，此即Frank-Kamenetskii边界条件。 Thomas系统爆炸判据的分析解 其无量纲化形式为 在系统中心处，边界条件仍为： Thomas系统爆炸判据的分析解 A类形状Thomas系统热自燃临界参数（指数近似） A类形状Thomas系统热自燃临界参数（指数近似） Thomas系统爆炸判据的数值解 Thomas系统中，系统的临界参数是由Biot数和值共同决定的，即: 当考虑无量纲活化能时，Tomas热平衡方程方程无法得到分析结果，这时需要采用数值模拟的方法来确定系统临界参数。 比较准确的计算临界参数的方法是根据分歧理论发展的数值方法。 Thomas系统爆炸判据的数值解 对于一维A类形状的Thomas系统，具体的计算式为： Thomas系统爆炸判据的数值解 上面3个方程式为热平衡方程的“非定常”形式，函数 为 的正值特征函数，C1为一个正的常数，计算中可取为1。三个方程需同步求解，离散格式采用Crank-Nicolson半隐格式，用三对角矩阵算法（TDMA算法[21]）求解每一时间步的变量值。计算时首先设定Bi和 ? 的大小，随着计算时间t的增加，变量?和?都逐渐收敛，?的收敛值即为?cr，?的收敛值为?cr，?cr(0)即为此时的临界中心温度。 平板的?cr随Bi和?的变化情况 平板的 随Bi和的变化情况 圆柱的 随Bi和?的变化情况 圆柱的 随Bi和?的变化情况 球的 随Bi和?的变化情况 球的 随Bi和?的变化情况 总结 通过以上的讨论可知，均温系统是研究放热反应系统发生热自燃的最简单的系统，其处理方法比较简单。在研究化学反应系统的热自燃问题时，经常采用对研究的系统进行简化，即将其简化为均温系统来处理。 总结 首先 详细介绍了均温系统，即Semenov模型下由反应性化学物质组成的体系的热平衡方程，并通过该方程研究了由于反应性化学物质的自身发热反应而引起的热自燃现象，求出了体系在不同条件下发生热自燃的临界条件和临界条件的数学表达式、热自燃发生的着火界限、绝热和非绝热系统的着火延滞期等。 总结 其次 对非均温系统，即Frank-Kamenetskii系统和Thomas系统的热自燃临界参数进行了计算。并就A类形状和非A类形状进行了分别研究。 总结 对A类形状； 为了得到分析解，使用指数近似下得到判断热自燃临界性的参数为Frank-Kamenetskii数 总结 对非A类形状 引入当量球的定义可以对临界参数进行计算。对非均温系统热自燃问题的数值计算方法进行了研究。采用直接迭代法和根据分歧理论发展的数值计算法对Frank-Kamenetskii系统和Thomas系统的热爆炸临界参数进行了计算，并将计算结果与文献值进行了比较，吻合得较好。对Frank-Kamenetskii系统热爆炸的非定常问题进行了数值计算，得到了不同Frank-Kamentskii参数下系统温度随时间变化情况。 * * 中心 边界 温度 Thomas模型温度分布 环境温度T0 体系内温度分布 对A类形状，热平衡方程为： 在指数近似(ε=0)下，可得到方程的分析解。 δcr随Biot数的变化 θ0,cr和θ1,cr随Biot数的变化 平板的 随Bi和的变化情况 3.322 2.000 0.878 ?cr 球 无限圆柱 无限大平板 形状 *