（动力工程及工程热物理专业论文）用于燃烧模拟的动态自适应机理方法及其应用.pdf

implementation and application of an adaptive chemistry scheme for combustion modeling a thesis submitted to chongqing university in partial fulfillment of the requirement for the degree of master of engineering by shi wan-ling supervised by prof. xiao long gou major: power engineering and engineering thermophysics college of power engineering of chongqing university, chongqing, china may, 2011 重庆大学硕士学位论文 中文摘要 i 摘 要 为了掌握燃烧过程的精细流场结构及污染物的排放特性，实现一次能源的高 效清洁利用，需要对碳氢燃料燃烧反应流进行详细的分析。其中将基元燃烧模型 与计算流体力学模型相耦合的方法被证明是对反应流进行精确数值模拟的一种有 效途径。然而碳氢燃料的详细化学反应机理常常包含数目庞大的基元反应与组分， 在每个计算网格中求解构建化学动力学过程的具有强刚性的非线性常微分方程组 将耗费大量的 cpu 时间，为数值模拟带来巨大障碍。为此，对详细化学反应机理 进行简化处理已成为保证机理可用性的关键。针对目前全局机理简化方法难以兼 顾简化程度与计算精度的矛盾与基于多工况的离线最优机理简化方法搜寻过程中 花费时间巨大的问题，本文对传统的动态自适应化学机理方法(dynamic adaptive chemistry，dac) 进行改进，提出了一种新的变步长动态自适应化学机理方法 ( variable step dynamic adaptive chemistry，vsdac)。本文的主要工作包括： 1）分析了机理简化阈值与初始搜索组分集合的选取对动态自适应化学机理方 法性能的影响，初始搜索组分集合选为燃料+co+ho2时，动态自适应机理方法具 有较好的性能，同时阈值存在一个最佳值。 2）在充分比较常见机理简化方法的基础上，基于直接关系图深度搜索算法完 成了即时简化方法的设计。在燃烧过程的不同时刻点，根据当地的热化学状态对 详细机理进行简化，得到即时的当地简化机理以应用于下一时间步的计算。 3）针对不同燃烧阶段基元组分和基元反应活跃程度变化不尽相同的特点，对 dac 方法进行改进，根据当地温度变化的情况对即时机理简化的时间步长进行动 态调整，而得到一种新的变步长动态自适应化学机理方法 vsdac。 4）为了验证 vsdac 方法的机理适应性、精确性和快速性，应用该方法，在 不同初始压力、温度及当量比条件下对正庚烷(558 种组分/2539 个反应)、异辛烷 以及汽油替代燃料(1034 种组分/4236 个反应)的定压绝热大空间自燃着火过程进行 了模拟，并与全过程采用详细机理的计算结果，以及传统 dac 方法的计算结果进 行了比较。 比较结果证明，与全过程采用详细机理的模拟方法相比，vsdac 方法得到了 精确的温度、各主要组分质量份额随时间变化规律以及着火延迟时间与最高温度 等参数。该方法与传统 dac 相比，在保证模拟精度的同时，模拟效率得到进一步 提高。证明了该方法的正确性、高效性。 重庆大学硕士学位论文 中文摘要 ii 本文提出的 vsdac 方法可用于基于详细反应机理的燃烧过程模拟， 为解决燃 烧反应流模拟的精度和计算耗时之间的矛盾、拓展详细机理的应用范围提供了一 种新的思路。 关键词：关键词：机理简化，直接关系图，自适应，燃烧 重庆大学硕士学位论文 英文摘要 iii abstract direct numerical simulation using detailed mechanism is proved to be one of the most efficient ways to study the dynamic process of reactive flow, and it is more and more popularly used in combustion research field to understand the working process since it can accurately predict flow structure and pollutant emissions. however, in general, detailed mechanisms for hydrocarbon fuels contain hundreds of species and reactions which make up a large system of nonlinear stiff ordinary differential equations (odes). it is difficult to solve the odes for each computational fluid dynamics （cfd） cell since the process will cost a lot of cpu time for reactive flows simulation. therefore, to simplify the complexity of the detailed mechanism has become the key to the numerical simulation for reactive flow. now, a lot of traditional mechanism reduction methods are used to get skeletal mechanisms. however, it is difficult to give consideration to both reduced simulation time and high accuracy while using a fixed skeletal mechanism in the whole simulation process. at the same time, optimization-based adaptive chemistry (ac) method takes a lot of time for the mechanism searching process. as an improvement for traditional dac scheme, a variable step dynamic adaptive chemistry (vsadc) scheme for combustion process simulation is presented in this thesis. the main jobs are listed as follows: 1) the effect of the choice of threshold value and the selection of initial searching species were analyzed. the dac scheme produced highest-performance with an appropriate threshold value and fules+co+ho2 as initial searching species set. 2) compared with different other traditional mechanism reduction methods, directed relation graph breadth-first search algorithm was chosen as the instantaneously mechanism reduction method. based on directed relation graph breadth-first search algorithm, a detailed mechanism can be reduced into instantaneously locally accurate sub-mechanisms which are used in the combustion modeling of the next time steps. 3) as to the problem that active species and reactions number vary differently in different sections of ignition process, the reduction time step is dynamically adjusted according to the local temperature accumulation. according to this idea, a variable step dynamic adaptive chemistry (vsadc) scheme was put forward. 4) under various different initial pressures and temperatures, vsdac method is tested by simulating homogeneous bulk space auto-ignition process with n-heptane 重庆大学硕士学位论文 英文摘要 iv (with 558 species and 2539 reactions), iso-octane and surrogate gasoline fuel prf90 (with 1034 species and 4236 reactions) as the detailed mechanisms. compared with using the detailed mechanisms, the process with vsdac can reproduce temperature profiles, species mass fractions, ignition delay times and peak temperatures with high accuracy. at the same time, this method makes the whole calculation process much more efficient than dac scheme. the vsdac method, presented in this thesis, can solve the conflict between simulation time and accuracy. it provides a new idea to widen the use of detailed mechanism for reactive flows simulation. keywords：mechanism reduction; directed relation graph; adaptive; combustion 重庆大学硕士学位论文 目 录 v 目 录 中文中文摘要摘要 . i 英文摘要英文摘要 . iii 1 绪绪 论论 . 1 1.1 燃烧动力学模拟国内外现状燃烧动力学模拟国内外现状 . 1 1.1.1 燃烧动力学模拟的研究意义及现状 . 1 1.1.2 燃烧动力学模拟面临的挑战及解决方法 . 1 1.2 燃烧过程的多尺度模拟及并行运算方法燃烧过程的多尺度模拟及并行运算方法 . 2 1.3 机理简化方法的研究现状机理简化方法的研究现状 . 3 1.4 本文的研究内容本文的研究内容 . 6 2 化学化学反应动力学与动态自适应机理方法反应动力学与动态自适应机理方法 . 7 2.1 化学反应动力学化学反应动力学 . 7 2.1.1 基元反应速率 . 7 2.1.2 阿累尼乌斯定律 . 8 2.1.3 多步反应机理的反应速率 . 8 2.2 基于基于误差误差传播的直接关系图方法传播的直接关系图方法 . 9 2.2.1 drg 方法 . 9 2.2.2 drgep 方法 . 11 2.2.3 drgbs 算法 . 11 2.3 动态自适应机理方法动态自适应机理方法 . 14 2.4 动态变步长自适应机理方法动态变步长自适应机理方法 . 14 2.5 本章小结本章小结 . 15 3 燃烧模型与动态自适应机理计算过程的实现燃烧模型与动态自适应机理计算过程的实现 . 16 3.1 数值模拟控制方程数值模拟控制方程 . 16 3.2 封闭系统中均匀气体混合物封闭系统中均匀气体混合物燃烧模型燃烧模型 . 17 3.3 动态自适应机理计算过程的实现动态自适应机理计算过程的实现 . 18 3.3.1 燃烧过程计算程序的实现方法与整体结构 . 18 3.3.2 动态自适应机理的程序实现 . 20 3.4 本章小结本章小结 . 22 4 动态自适应机理方法参数影响研究动态自适应机理方法参数影响研究 . 23 4.1 封闭均相大空间正庚烷自燃着火模型参数设置封闭均相大空间正庚烷自燃着火模型参数设置 . 23 4.2 阈值的影响阈值的影响 . 23 重庆大学硕士学位论文 目 录 vi 4.3 初始搜索组分的影响初始搜索组分的影响 . 26 4.4 结论结论 . 28 4.5 本章小结本章小结 . 28 5 变步长动态自适应化学机理方法在燃烧模拟过变步长动态自适应化学机理方法在燃烧模拟过程中的应用程中的应用 . 29 5.1 vsdac 方法与方法与 dac 方法的比较方法的比较 . 29 5.2 vsdac 方法的准确性与适用性分析方法的准确性与适用性分析 . 32 5.2.1 正庚烷自燃着火过程的模拟 . 32 5.2.2 异辛烷自燃着火过程的模拟 . 36 5.2.3 汽油替代燃料自燃着火过程的模拟 . 40 5.3 vsdac 方法快速性分析方法快速性分析 . 44 5.4 结论结论 . 46 5.5 本章小结本章小结 . 46 6 结论与展望结论与展望 . 47 6.1 全文总结全文总结 . 47 6.2 主要创新点主要创新点 . 47 6.3 工作展望工作展望 . 47 致致 谢谢 . 48 参考文献参考文献 . 49 附附 录录 . 53 作者在攻读学位期间发表的论文目录作者在攻读学位期间发表的论文目录. 53 重庆大学硕士学位论文 1 绪 论 1 1 绪 论 1.1 燃烧动力学模拟国内外现状 1.1.1 燃烧动力学模拟的研究意义及现状 碳氢燃料等一次能源的燃烧是人类获得二次能源的重要途径之一。现代社会 的动力来源，主要来自于矿物燃料的燃烧，其应用普及至各种领域，如火力发电 领域、化学化工领域、各种动力机械的发动机等，都是依靠固体、液体、气体燃 料的燃烧产生热能作为动力1。为了实现一次能源的高效清洁利用，掌握燃烧过程 的精细流场结构，以及污染物、有毒物质等微量物质的排放特性等，需要对反应 流进行详细的分析。计算机数值模拟技术作为一种有效手段具有价廉、可靠、方 便的优点。近几十年来，实际燃烧装置内反应流的的数值模拟取得了长足的发展。 然而燃烧过程的传统模拟方法由于常采用简单的总包反应，往往难以全面、精确 地描述燃料与流动过程的完整燃烧特征，尤其在污染物排放、分布的预报等方面 存在巨大的缺陷2。鉴于此, 应用详细化学反应机理代替总包反应对动力机械，如 内燃机、锅炉、燃料电池等，在不同环境、不同运行工况下的燃烧、流动特性进 行数值模拟，找出影响着火延迟时间、火焰传播速度、污染物生成、反应完全程 度等燃烧特性的详细机理性原因，以提高燃烧效率、降低污染物排放、实现节能 减排，已经成为目前国际燃烧界的热门研究方向，其中将基元燃烧模型与计算流 体力学(computational fluid dynamics, cfd)模型相耦合的方法被证明是对反应流进 行三维精确数值模拟的一种有效途径3。 国内外众多学者对该方法的应用与性能进 行了研究。国外学者 law 等4应用该方法研究了火焰形状对非预混扰动流烟灰的 形成与传递特性影响情况， 证明了该方法的精确性； smith 等5基于包含 19 个组分 的乙烯简化机理对瞬变的乙烯喷射火焰进行了三维的精确数值模拟；国内学者樊 建人6-8等常年致力于基于详细反应动力学的燃烧流直接数值模拟技术的研究。 1.1.2 燃烧动力学模拟面临的挑战及解决方法 然而碳氢燃料的详细化学反应机理常包含数目庞大的基元反应与组分，基于 详细机理的反应流燃烧动力学模拟遇到了计算时间和计算资源上的巨大挑战。首 先，详细机理的组分和基元反应数目众多，在每个计算流体网格计算庞大的描述 模型化学动力学的常微分方程组给数值模拟带来了难以想象的大计算量与存储 量，对计算的硬性条件提出了要求；其次，在基于详细机理的反应流计算中，实 际化学动力学模型中一般都既包括反应速率快、特征时间尺度小的基元反应，也 包括反应速率慢、特征时间尺度大的基元反应，将不可避免地面对反应速率的多 时间尺度及组分浓度的多尺度分布的问题。反映到微分方程中就是方程的特征值 重庆大学硕士学位论文 1 绪 论 2 差别明显，特征矩阵的条件数很大，形成了计算的刚性问题，如何处理计算中的 刚性问题也是反应流模拟面临的巨大障碍；第三，燃烧反应流过程常常还具有极 强的空间尺度分布特性，三维反应流常常包含成千上万的不同网格，在各个精细 流体网格上求解构建化学动力学过程的具有强刚性的非线性常微分方程组 (ordinary differential equations, odes)，将耗费反应流模拟大量的 cpu 时间，从而 影响到计算效率。 且由于燃烧和流动的复杂耦合作用, 对反应流的研究面临着很大 的挑战，为燃烧流数值模拟带来巨大障碍9。 为了简化反应流的计算， 目前常采用三种解决方案： (1)寻求新的求解刚性 ode 方程的算法，如多尺度算法10；(2)基于目前高速发展的计算机技术，利用高性能 的计算机机群实现并行运行技术11, 12，减少反应流直接数值模拟的用时；(3)对详 细机理进行简化，如准稳态假设方法13, 14、敏感性分析方法(sensitivity analysis, sa)15, 16、主成份分析法(principle component analysis, pac)17、计算奇异摄动方法 (computational singular perturbation, csp)18, 19、直接关系图方法(directed relation graph, drg)20, 21、基于遗传优化思想的简化算法22, 23、自适应简化法24-26、自适 应列表法27 (存储不同工况点的计算结果以再次利用)及动态自适应机理方法等。 1.2 燃烧过程的多尺度模拟及并行运算方法 为了节约基于详细反应机理的燃烧流模拟时间，中国学者苟小龙10等设计完 成了一种动态多时间尺度的模拟方法。该方法通过在每一个计算时间步在线计算 当前各个基元组分的特征时间，并依据其特征时间将组分和对应的反应分为多个 组。在每一个特定的输出时间步内，从快(特征时间小)到慢(特征时间大)采用不同 步长逐一进行计算，如图 1.1 所示，其中 tf，tm和 ts分别为快速组、中速组和 图 1.1 多时间尺度模拟过程示意图 快速 tf 中速 慢速 0 时间 tout tm ts fig.1.1 diagram of multi-time scale (mts) method 重庆大学硕士学位论文 1 绪 论 3 慢速组的时间步长。对于每一个组而言，mts 算法中时间步长均小于该组内每一 个组分的特征时间，从而保证了整个积分过程的收敛性。相比较于采用同一积分 步长的传统显式积分算法，mts 算法中对具有不同特征时间的组分采用了不同的 积分步长，能够较大程度地降低计算时间。该方法有效地解决了组分间特征时间 的巨大差异，以及同一组分在不同时间段内其特征时间的变化给计算稳定性和精 确性带来的问题，大大地提高了计算的效率。 燃烧流的多维数值模拟占用大量内存、耗时巨大，简单的单处理器的 cpu 不 能满足燃烧流数值模拟的硬件要求，所以基于高性能计算机采用并行计算技术是 解决这一难题的重要手段。高性能计算机有向量机、对称多处理器、集群和大规 模并行处理器等，其中集群和对称多处理器较为普遍。对称多处理器高性能计算 机的多个处理器共享同一个内存空间，而集群是一种分布式内存结构，各个节点 之间是通过网络互联，在各个节点间进行频繁的消息传递来交换中间变量结果。 同时在反应流的模拟中，首先在每个化学计算反应时间步长内求解每个网格内的 详细化学反应动力学模型；然后将其计算结果输出给流动控制方程，作为流动控 制方程的源项；最后在每一步流动计算时间步长内求解流动特性。燃烧流的模拟 过程涉及大量的燃烧模型参数与流动模型参数的交互交换，所以集群更适用于燃 烧流模拟的 mpi 的实现。 但文献11, 12的研究表明耦合化学反应动力学的 cfd 并 行计算加速比与 cpu 数量并非成正比，需要优化匹配 cpu 数量、并行模式、网 格规模、反应机理才能获得最大收益。 虽然，随着计算机技术的发展，利用高性能计算机基于多尺度模拟技术10与 并行计算技术11, 12可以大大地缩短反应流模拟的时间，但是目前常规的计算条件 仍然难以完成基于详细化学反应机理的碳氢燃料反应流的数值模拟。对于复杂燃 烧机理的庞大体系，在保证详细机理的全面性以及可靠的预测精度下，对详细反 应动力学机理进行简化以实现反应流的高效、精确模拟，成为国际燃烧学界 c. k. law、 t. tur nyi、lam、t. f. lu 等著名学者与他们的学术团队为之努力的方向之 一。 1.3 机理简化方法的研究现状 目前国际上出现了多种机理简化方法。 1985 年美国 princeton 大学 lam18教授 提出了用于求解刚性常微分方程的 csp 方法。1994 年，lam 等28人将该方法应用 于反应动力学机理的简化过程中，通过对基元反应数目的简化，实现详细化学反 应机理的简化。csp 方法基于基础向量对详细反应机理中每步基元反应进行时间 尺度的分析，将所有基元反应分为三组：快反应、慢反应、反应中不起作用的反 应，然后相较于燃烧特性的研究尺度，剔除详细反应机理中的快反应，将剩余的 重庆大学硕士学位论文 1 绪 论 4 慢反应集合组成简化机理。相比于传统方法，该方法不依赖于人的经验和直觉。 但该方法的时间尺度域值的选取仍然依赖于经验，而且域值的选取与简化机理的 准确程度没有直接的联系，所以很难找到最优域值。近年国内学者对 csp 方法也 做了很多研究，如用 csp 方法对甲烷/空气预混燃烧系统进行分析29；对 co、h2 在空气中燃烧的化学反应、激波管中的甲烷的点火燃烧反应实例进行了简化30； 对 co+ch4+空气的混合燃烧系统进行简化31；研究得到的简化模型具有较高的精 度，能有效提高计算效率，证实了 csp 分析方法在燃烧系统中应用的可靠性和优 越性。 由于基元反应的各组分与各反应之间具有强耦合性，为了对基元反应的机理 进行研究，1990 年，tam s tur nyi 15提出了敏感性分析方法，研究反应机理中反 应参数变化对整体计算结果的影响程度，即计算结果对参数的敏感性，直观地反 映某一参数的误差对系统精确性的影响，去除详细机理中对系统误差影响较小的 基元反应，以达到简化的目的。 敏感性系数定义为(1.1)式，其中 i z表示系统的因变参数，包括各组分浓度、系 统温度等； j k表示系统参数，包括初始温度、压力、化学反应速率常数、指前因 子等。(1.1)式反映了 j k的微变对 i z产生影响的量的度量。 , i i j j z k (1.1) 其中，为了分析的方便，将敏感性矩阵取对数坐标得到无量纲形式，如(1.2) 式。 , ln ln j ii i j jij k zz kzk (1.2) 同时，为了避免产生组分浓度接近零时，敏感性系数却很大的人为数值误差， 用参数区域最大值代替即时值，将(1.2)式修改为(1.3)式。 , max ln ln j ii i j jj k zz kzk (1.3) 该方法已经实现商业化，大型燃烧流模拟软件 chemkin 中具有敏感性分析 模块。但该方法在机理简化过程中要数值求解敏感性系数矩阵，导致计算量过大， 所以该方法适合对详细机理进行全局简化，不适合进行在线简化。 2005 年，美国学者 law 和 lu20提出的直接关系图法(directed relation graph, drg)，直观地描述了各组分之间的耦合情况，方便地实现了对详细机理的快速高 效简化，得到了广泛的应用。2006 和 2008 年, law32和 pitsch33等学者采用基于 drg 方法扩展的 drgep 方法，针对充分搅拌反应器(psr)中碳氢燃料着火等, 进 行了机理简化研究。2007 年，xin 等21基于直接关系图法对甲烷氧化的详细机理 重庆大学硕士学位论文 1 绪 论 5 进行了简化，证明了 drg 方法的精确性与有效性。然而 drg 方存在着一些重大 缺陷34：(1)对机理进行简化时设定的域值 与简化机理准确程度无明确的对应关 系，难以得到最优的简化机理；(2)drg 方法认为基干机理中每个选定的组分的重 要性都相同，对于某些与选定组分强联系的组分，即使其反应线路不重要也将被 选入简化机理中，这就导致了简化机理中容易出现一些冗余组分；(3)drg 方法无 法描述第三体与组分之间的非化学联系(如气相反应中，由于组分总摩尔数变化导 致的定压条件下的膨胀)对机理简化的影响。基于这些缺陷，对 drg 理论进行分 析， 产生了 two-stage drg35、 drgsa36、 drgep33等 drg 的改进算法。 two-stage drg 算法通过对初步简化机理再次运用 drg 算法而得到更优的简化机理； drgsa 方法集合敏感性分析方法与 drg 方法的优点，在对机理进行简化时，先 用 drg 方法在小的设定阈值的条件下，对详细机理进行简化得到相当庞大的预简 化机理，然后再用敏感性分析方法对采用 drg 方法得到的预简化机理进行分析， 减少其中的冗余组分，解决了 drg 方法的前两点缺陷；drgep 方法首先选择重 要组分，再通过深度搜索分析其余组分与重要组分的紧密程度对机理进行简化， 该方法克服了 drg 算法的第二个缺陷。目前 drg 系列的方法因为其计算方法简 化、思路清晰易懂、计算快速的特点，使得其容易实现动态自适应机理，因此得 到了国内外学者的喜欢。本文的研究也将基于该方法，下文将详细介绍该方法的 思想。 上文介绍的三种传统的机理简化方法大多基于一系列特定反应条件(温度、压 力、燃空比等)，通过分析各个组分或反应对整体计算结果的影响程度，删除对结 果影响小的组分和基元反应，以获得适用于全局的统一简化机理。在燃烧反应流 中，温度、压力、组分摩尔分数等反应条件随时间、空间将不断地发生变化，适 用于全局的简化机理必然包含每个状态点的重要反应和重要组分，这将导致简化 机理过于庞大，简化程度受到极大的限值。因此，在反应流的不同区域，针对不 同的条件采用与之相适应的最优简化机理，将是一条高效的、切实可行的途径。 相对于整个模拟过程采用一种机理的方法， green 教授提出的最优机理简化方法24, 25， 根据反应流热化学条件离线获得覆盖所有反应条件的简化机理集， 并在模拟过 程中，根据当地热化学条件在机理库中搜寻相应的当地最优简化机理。国内学者 乔瑜等基于 green 教授的最优机理简化方法37, 38，结合敏感性分析方法，对甲烷、 氢气等机理进行了简化。然而最优机理简化方法最优机理搜索算法的设计遇到了 瓶颈，使得机理的搜寻过程占用了大量反应流模拟的时间，其实用性受到了的限 值。 2009 年美国学者 long liang3, 39提出的基于直接关系图深度搜索算法(directed relation graph breadth-first search, drgbs)的动态自适应化学机理方法， 不需要预先 重庆大学硕士学位论文 1 绪 论 6 提供简化机理集，在对反应流进行模拟时，对详细机理执行动态的机理简化，避 免了反应流模拟过程中机理搜索所用的时间。本文将在下面章节详细介绍该方法 的思想。 综上所述，机理简化方法在近年得到了深远的发展，很多国内外的学者在该 领域做出了重大贡献，并取得了较好的成果。机理简化方法为燃烧模拟在燃料特 性、反应流动力学特性研究工作中的应用创造了条件。目前，各种机理简化方法 仍存在一定的局限与不足，探寻高效、快速、实用的机理简化方法仍然是燃烧科 学领域一项重要的任务。 1.4 本文的研究内容 通过分析燃烧过程模拟的特点和目前传统机理简化方法的优缺点，探索新型 的高效、实用的机理简化新思路，为燃烧反应流的精确模拟提供必要的条件。本 文的研究工作主要包括： (1) 在对传统简化方法进行完整分析的前提下，提出一种新型的自适应机理方 法； (2) 通过对一系列详细机理进行简化，验证该方法的有效性与精确性，分析该 方法的优缺点，并得到完善的措施； (3) 将新型方法与现有方法进行比较，分析新型方法的准确性、高效性和快速 性。 重庆大学硕士学位论文 2 化学反应动力学与动态自适应机理方法 7 2 化学反应动力学与动态自适应机理方法 2.1 化学反应动力学 燃烧过程是人类利用煤、汽油、天然气等一次能源的有效手段，在燃烧过程 中，化学反应速率控制整个燃烧过程。化学动力学(chemical kinetics)是对基元反应 及其反应速率的研究40，是物理化学的一个专门领域。它的研究对象是物质性质 随时间变化的非平衡的动态体系，其中时间是化学动力学的一个重要变量。 目前研究化学动力学的方法主要有两种。一种是唯象动力学研究方法，也称 经典化学动力学研究方法，它的思路为首先从化学动力学的原始实验数据(浓 度与时间)的关系出发，然后经过分析获得某些反应动力学参数(反应速率常数