化学反应动力学方程组的通用数值解法.pdfVIP

化学反应动力学方程组 的通用数值解法 李 汝 雄 一兔曲 一 、 前 曰 当 已知 化 学反 应 的机 理 , 包 括 各基元反 应 的速 度常数及 各物种 的初始浓 度时 , 需要求解 动 力学方程组 , 从 而得 到各物种浓 度 与时 间的关系, 只有少数简单 的动 力学方 程 组可 以求 出 解 析解 多数情 况不 行, 这 时, 一 种 办法 是通 过近 似方法 如稳态 法 等 求解, 但近 似方 法 。 只有在一定 条件 下才 比较 准 确, 另一 种方 法 就 是通过数值积 分方法求解 通 常 的数值 积分 方 法 必须写 出具体 的微分 方程 组才 能进 行 。 提 出了转 移矩 。 , , 阵方 法 〔〕 但 未提 出写 出转移矩 阵 的通 用公 式 计 算精确度相 当于欧 拉法 我 们 由质 量 一 、 作用定律, 提 出 了通用 的计 算公式 编制 了通 用程序 , 用 户 只要能 正确输入 已知 数据 , 即 可 。 得 浓 度 时 间数据 , 而不 必写 出具体 的动 力学方程 组 , 使 用 方便 本文提 出的方 法在反 应机 理 的研 究 中也是有用 的, 当我 们用矩 阵法 〔〕等方 法 拟 定 出 若干 套较为 合理 的机理 后 , 对每 套机理