MATLAB用于振荡反应模拟.pdfVIP

年 月 四川大学学报（自然科学版） (!! ) +N 0(!! S 第 卷第 期 （ ） / M;N:9K;6F5@ON9P5Q8:354 I94N:9KF@58@8E545; T;K0/ I;0 R 文章编号： （ ） !#!$%’%(!!!$!)(#$! *+,-+.用于振荡反应模拟 黄允中，张元勤 （乐山师范学院环境与生命科学系，四川乐山 ） %/!!! 摘要：解刚性常微分方程组已成为研究复杂化学反应的重要手段 作者讨论了使用 0 *+,-+. 对振荡反应进行数值模拟的方法，应用该方法对硫化学反应非线性动力学行为进行了模拟0 关键词：振荡反应；数值模拟；刚性常微分方程； *+,-+. 中图分类号： 文献标识码： 1%2 + —个反应机理模型的正确性可以通过实验方法或数值模拟来验证 应用数学理论和数值方法定性或 0 定量的预测或模拟化学物种在反应过程中的浓度 时间关系是反应动力学研究的主要内容之一 描述反应 $ 0 过程中化学物种的浓度 时间关系的数值方法有确定性方法和随机性方法两类 确定性方法即是将反应组 $ 0 分与时间的关系问题归于微分方程的求解 对于不太复杂的反应可以通过求解微分方程得到解析解，但对 0 于很多复杂反应只能得到数值解 随着计算机技术的发展，使得对很多复杂化学反应的研究可以通过数值 0 计算模拟而得以完成0 在一个包含多个化学物种和多步基元反应的复杂反应体系中，由于各个反应的速率常数和各个物种 的浓度相差可达 个以上的数量级，这样的微分方程组的数值解法存在严重的问题，普通的显示解法（如 /! 欧拉法）在控制步长达到每一步的精确度时不能保持数值运算的稳定性，这在微分方程求解中归结为刚性 （ ）方程 数值求解刚性微分方程组的主要方法有隐式方法（隐式欧拉、梯形、中点、龙格 库塔）、吉尔 34566 0 $ ［］ （ ）法和半隐式的龙格库塔法 / 789: 0 对于化学工作者来说运用商业化的程序来解决复杂化学反应的数据研究是很好的选择，大型科学计 算软件 也为我们提供了基本的求解刚性常微分方程组的 函数，但将其用于研究化学振荡问 *+,-+. ! ［］ 题鲜见报道( 我们讨论用 对振荡反应进行数值模拟的正确方法，并将该方法应用于对硫化学 0 *+,-+