化工计算方法化工计算方法微分方程.pptVIP

6 常微分方程数值解法 例：在一个封闭系统中，三个组分浓度分别为C1、C2、C3。系统在特定频率光照下发生反应，求每种物质浓度随时间的变化规律。 例：一维均匀介质稳态导热问题，设其一端绝热，另一端恒温为T1，此问题的数学模型为: 一阶常微分方程初值问题的一般形式 微分方程数值求解的基本原理和步骤 6.2 欧拉（Euler）法 欧拉（Euler）格式 Euler格式 Euler格式2 Euler格式3 步长 h＝0.1，求 x=0~1间各节点上的数值解。 再用改进欧拉法求解 运行结果 例6-2 已知在管式反应器中进行液相反应 A→R+S，反应为吸热反应，反应管外油浴温度为340℃，假定已知管内温度与转化率的关系为 例6-2 求解 6.3 龙格-库塔（Runge-Kutta）法 四阶龙格-库塔格式 四阶龙格-库塔-吉尔格式 常微分方程组及及高阶方程求解 求常微分方程组的四阶龙格-库塔格式 高阶微分方程求解 对于一阶以上的高阶微分方程，通常可将其化为一阶方程组求解。对以下形式的二阶常微分方程的初值问题 例6-3：对例6-1改用四阶龙格-库塔法计算, 比较改进欧拉法与四阶龙格-库塔法两种方法计算的精度。 解 为编程计算方便，直接调用MATLAB函数库中龙格-库塔函数。 结果中列出了解析解y*以进行精度比较。 例6-4 求解以下微分方程组，区间为[0,0.2]，步长取h = 0.01，要求每隔5步输出一次计算结果。 例题 6-7 甲烷水蒸气催化转化反应制氢，计算：1、从反应管入口到转化率xA大于0.644时，反应管的轴向温度(T)及转化率(xA)分布；2、所需要的反应管长度和催化剂用量。 已知物料衡算和热量衡算式如下 例题 6-7程序 运行后输出结果 例题 6-8 用一维拟均相模型计算以下两种情况的乙苯脱氢制苯反应器。 （2）采用管式换热反应器 解 本例属于反应器设计的问题 令 本章要求 清楚数值方法求解常微分方程的基本原理和方法 了解Euler法、改进Euler法和龙格-库塔法求解常微分方程（组）初值问题的格式 能够运用改进Euler法和龙格-库塔法求解包括化工问题在内的常微分方程（组） 初始条件为： (c) 绝热反应器计算 计算绝热反应器时，反应器的总传热系数为U=0，管径D＝1.211m。 (d) 管式换热反应器计算 计算管式换热反应器时，反应器与管外加热烟气间的总传热系数为U=2.58，采用管径D＝0.1m 的144根反应管。 参考程序及计算结果见教材p.75。 相应的上述两个方程须变换为 dL/dxA，dT/dxA 复习思考题 为什么要研究常微分方程的数值解法?初值问题数值解的基本思想是什么？ 微分方程的解析解与数值解在形式上什么区别？一般来说，龙格-库塔法与改进欧拉法在精度上有什么不同?为什么? * * 6.1 基本概念及求解思路 微分方程：表示未知函数与未知函数的导数以及自变量之间关系的方程。（含有未知函数的导数的方程) 只含一个自变量的导数的方程叫做常微分方程。 方程中出现的导数的最高阶数叫微分方程的阶 一阶常微分方程 二阶常微分方程 根据所给的(已知)定解条件，将常微分方程分为初值问题和边值问题两类 给出自变量两端的函数值或导数值，称边值条件 # 初值问题 给出自变量为 0 或通过变换可将自变量变为0（自变量一端）时的函数值或导数值，称初值条件； 边值问题 已知自变量 t 为0时的变量值（初始条件），为常微分方程组的初值问题。 # 解：此问题的数学方程为 已知： C1 C2 C3 k1C1 k3C22 k2C2C3 绝热 T1 L 常微分方程初值问题的一般形式表示为 已知自变量两端的条件，属常微分方程的边值问题 常微分方程的初值问题和边值问题求解方法不同,本课程只涉及初值问题的求解。 解的形式为 y =g ( x），是在x y平面上通过初始点（x0，y0）的一条曲线，是连续函数。 x1 x2 … xi xn … y =g ( x） 数值解并不是求这个连续函数的形式 y =g ( x） ，而是寻找函数在一系列离散点 xi 上的函数的近似值 yi 用有限个离散点上的近似值的集合代替连续解。 # 将求解区间分成若干离散点 x1 ，x2， …，xi ，…； 在离散点上用差商近似地代替导数，将微分方程化为代数方程（差分方程），差分方程中的未知量就是未知函数在离散点 xi 上的近似值 yi ； 常微分方程初值问题数值求解基本思路是在离散点上用差商代替导数，基本方法则是递推，从初始值开始，沿着节点排列的顺序向前推进。 3. 求解代数方程得到未知函数在离散点上的近似值 y