数学描述1MathematicalFormulationI.PPT

* 冶金数值—— 数学描述—— 相间传输 气－固相间传输的数学描述 气－固间化学反应之铁矿石气体还原 一界面未反应核模型 界面化学反应 一般铁氧化物还原反应中，氧化亚铁的还原是限制性环节 c c1 c0 H2(CO)的 浓度变化 Fe2O3 Fe层 r r0 边界层 ce 式中，k、k’－分别为正逆反应速度常数。 当反应平衡时，kcCO,e=k’cCO2,e，k/k’=K(K为反应平衡常数)。于是按照质量守恒可以得到 还原速度方程 视铁矿石还原过程为准稳态，即传输及反应的各个环节中没有物质积累，总反应速度R 整理得到 反应速度为 * 冶金数值—— 数学描述—— 相间传输 气－固相间传输的数学描述 气－固间化学反应之铁矿石气体还原 一界面未反应核模型 c c1 c0 H2(CO)的 浓度变化 Fe2O3 Fe层 r r0 边界层 ce 边界层传质阻力 Fe层传质阻力 化学反应阻力 按照控制方程中源项单位，将化学反应速度处理成单位体积的生成速度R’，则 式中r无法直接测量，为此我们定义一个还原度f 式中W0为r=r0时的球中的氧重；W为r时球中的氧重。 于是最后得到 * 冶金数值—— 数学描述—— 相间传输 气－固相间传输的数学描述 气－固间化学反应之铁矿石气体还原 三界面未反应核模型 Fe FeO Fe3O4 Fe2O3 r0 r1 r2 r3 De1 De2 De3 k1 k2 k3 kf 一界面模型假设所有化学反应都发生在内核界面，而这显然是不合适的，实际的Fe2O3的还原过程是一个逐步还原过程，这必然涉及到还原程度不同的分层。 三界面未反应核模型就是认为存在三个反应界面层：Fe-FexO、FexO-Fe3O4、Fe3O4-Fe2O3。还原气体CO（或H2）通过边界层扩散到De3层（ Fe-FexO 层）界面，与FexO反应消耗一部分还原剂，产物CO2（或H2O）向外扩散，其余还原剂继续往中心扩散，通过De3到De2、De1…整个过程共有11个环节，其中扩散环节8个（DCO,界面层，DCO,Fe层，DCO,FeO层，DCO,Fe3O4层，DCO2,Fe3O4层， DCO2,FeO层， DCO2,Fe层， DCO2,界面层），界面反应环节3个。 根据一界面未反应核模型，可以很容易地推导各个环节的速率表达，然后按照准稳态假设条件时各环节速率趋于相等的条件推导出用各阶段还原度f1、f2、f3表达的三个还原速度表达式。 需要说明的是，关于各层有效扩散系数的表达需要用到前面的空隙度和迷宫指数的关系式Deff=Dεξ。 作业：参考教材详细推导三界面未反应核模型的完整数学表达。 * 冶金数值—— 数学描述—— 相间传输 流－流相间传输的数学描述 气－液相间传质 气－液相间传质采用的是双膜理论，同样是在准稳态的假设基础上推导而成。 界面的气－液两相分别定为I、II相，三个过程环节的速度式为 δ1 δ2 I II c1* c1 c2* c2 双膜理论模型 I相中反应物向相界面扩散 界面化学反应（一级可逆反应） II相中的扩散 综合反应速度有 式中，A——反应截面积 km1、km2——在I、II两相中反应物的传质系数 对于式中的传质系数km，液相可以采用边界层理论，km=Deff/δ；也可以采用湍流传质理论，km=(2/π)(αD)1/2，α=7.9×10-5Re1.678。气相则采用前面提到的兰兹－马歇尔式计算。另外，两相流体中传质系数的关系存在：(km1/km2)≈(ν1/ν2)0.5(D1/D2)0.7，式中，ν为运动粘度，D为分子扩散系数。 * 冶金数值—— 数学描述—— 相间传输 流－流相间传输的数学描述 渣金反应 渣金反应是钢铁冶金的核心反应之一，属于典型的液－液界面反应，可以借助双重阻力概念和双膜理论来进行描述。反应过程中，熔渣与钢液中组元分别向相界面传输并发生反应，而反应产物或是作为气体排出，或是进入渣相。因为渣金反应属于高温反应，人们认为反应环节一般不会成为限制性环节，即反应阻力可忽略不计，于是总速度可表达为： 渣金反应中km1/km2一般在20~100范围之内，通过分析阻力部分可知，反应平衡常数值的大小对于两个传质环节的速率构成了影响。一般当K远小于20~60时，过程由熔渣侧传质所控制，而K远大于100时，过程则由钢液侧传质所控制。 炼钢的氧化脱硫过程、氧化脱磷过程都是典型的渣金反应，都可以用双膜理论来进行公式推导，并获得所谓双重阻力公式，然后通过分析参事之间的关系获得实际操作过程中的效果判断。 * 下次课再见！ See you! 冶金过程数值模拟 * 冶金数值—— 数学描述—— 控制方程 控制方程 之 小结 回过头来看看前面所述的控制方程，无论是连续性方程、运动方程、能量方程或者溶质方程，他们的形式非常类似，于