流体输运—反应时空动力学质量守恒方程的数值模拟

摘要

模拟输运—反应耦合动力学的常用方法是有限差分或有限元等数值方法,方程离散的结果将问题转化为在某一时间步求解一套非线性代数方程。在多组分反应体系,在某一时间步所需求解方程的数目总不少于N_c·M_o利用准静态近似法,准静态之间的时间步长不受数值稳定性的制约。改进的Newton-Raphson方法是求解这些非线性代数方程最有效的方法。输运—反应方程差分后所得出的Jacobian矩阵具有序结构,适于块迭代方法的使用,如BSOR法。块迭代方法有利于对角线稠密亚矩阵直接和有效地通过LU分解等标准方法求逆阵。