MATLAB在复杂化学反应模拟中的应用_黄允中.pdfVIP

第 卷 第 期 ! ’! 乐山师范学院学报 1’234562! #$$% ’! #$%27889 年 月 !#$%’ ( )*+,% -*.,*$+ /’’*0* +,-.,/ 在复杂化学反应模拟中的应用 黄允中 张元勤 （乐山师范学院 环境与生命科学系，四川 乐山 0’#） 摘 要：解刚性常微分方程组已成为研究复杂化学反应的重要手段，本文讨论使用’()’* 的+,$-. 函数对复杂化 学反应进行数值模拟的方法，并对硫化学反应非线性动力学行为进行了模拟研究。 关键词：复杂化学反应，数值模拟，刚性常微分方程、 ’()’* /01-23 ’ !453:6111 557 !:5515:57 中图分类号： 文献标识码： 文章编号： （ ） 一个反应机理模型的正确性可以通过实验方 通的显示解法（如欧拉法）在控制步长达到每一步 法或数值模拟来验证。应用数学理论和数值方法 的精确度时不能保持数值运算的稳定性，这在微 定性或定量地预测或模拟化学物种在反应过程中 分方程求解中归结为刚性（23455）方程。数值求解刚 的浓度1 时间关系是反应动力学研究的主要内容 性微分方程组的问题当前有隐式方法（隐式欧拉、 1 梯形、中点、龙格$ 6789 之一。描述反应过程中化学物种的浓度 时间关系 库塔）、吉尔（ ）法和半隐式 :’; 的数值方法有确定性方法和随机性方法两类。确 的龙格库塔法 。 定性方法即是将反应组分与时间的关系问题归于 对于化学工作者来说，运用专业化的程序解 微分方程的求解。对于不太复杂的反应可以通过 决复杂化学反应的数据研究是很好的选择。大型 求解微分方程得到解析解，但对于很多的复杂反 科学计算软件+,-.,/ 提供了基本的求解刚性 应只能得到数值解。计算机技术的发展使得很多 常微分方程组的 函数，但将其用于研究化学振 :!; 复杂化学反应的研究通过数值计算模拟得到很好 荡问题鲜见报道 。 的进展。 ’ 求解刚性常微分方程组的 函数 基元反应遵守质量作用定律，非基元反应虽 然不符合质量作用定理，但非基元反应是多个基 在+,-.,/ 中求解微分方程的函数有多个， 元反应的复合，反应速率也能与反应物浓度发生 可以用于求解刚性常微分方程组的有以下几个： ———采用数值差分的方法（ 法）， 关联。动力学模拟也就归结为求解非线性微分方 =7’%2 6789 程组及初值问题。然而，像化学振荡反应这样的复 可变阶方法。 杂反应，是一个包含多个化学物种和多步基元反 =7!?2———采用修正的@=27AB9=CD 方法，解 应的复杂