化工数学模型--JY-5B.docVIP

5.3 解析解 模型求解很重要，是建模思想和学科理论通向应用的桥梁。 主要有两种方法：解析解和数值解。 解析解：应用方便，容易分析模型参数的影响。 形式：解析函数式（项数少）； 级数解（上百项仍容易计算）； 隐式解（难于使用）。 数值解：以列表或图象表示，便于对具体问题进行分析，但分析一般情况略嫌困难。 最重要的优点是，几乎所有数学模型都能求解。 与偏微分方程模型比，常微分方程模型比较简单。 即使是一阶常微分方程，复杂性也各不相同，无一般解法。 一阶、常系数、齐次、线性方程的初始值问题是相对较易求得解析解的， 高阶、变系数（系数为自变量的函数）、非齐次、非线性方程的边界值问题是难解的。 一阶常微分方程可用的特别解法： 全微分方程。 M(x,y)dx+N(x,y)dy=0 若存在函数U(x,y)，满足 dU=M(x,y)dx+N(x,y)dy 则原方程写为 dU=0 解为 U(x,y)=C 其充分必要条件为 可分离变量的方程。 M(x)dx+N(y)dy=0 f1(x)g1(y)dx+ f2(x)g2(y)dy=0 齐次方程。 M(x,y)dx+N(x,y)dy=0 M(?x,??y)=?nM(x,y), N(?x,??y)=?nN(x,y) 则原方程为n阶齐次方程。也可写为 例如， 求解时可利用变换： y=vx 成为可分离变量的 例5－7A 连串反应动力学 产物为A2, 最大收率、最佳转化率？ 第3个方程是多余的 若不关心浓度作为时间的函数, 而仅想知道c1、c2、c3间的比例关系，以及转化率X等指标, 则可将前两个动力学方程式相除，消去dt后得以c1为自变量的一阶线性、变系数常微分方程：一阶齐次方程 变换 v=c2/c1, 当?=k2/k1≠1时， 可分离变量， 代入定义v=c2/c1和转化率X=1- c1/c10，得 最终得到收率Y2与转化率间的关系： 当?=k2/k1=1时， (5－111) 可找到积分因子的方程 M(x,y)dx+N(x,y)dy=0 若存在函数?(x,y)，使 ?(x,y)M(x,y)dx +?(x,y)N(x,y)dy=0 成为全微分方程，即满足 则称?(x,y)为积分因子 一阶线性常微分方程 其积分因子为 解公式为 或代入初始条件(x=0, y=C0)的 初始条件： t=0时, c1= c10, c20=0 应用(5－114)式和初始条件，即得到 二阶常微分方程的分类： 非线性方程（因变量不显式出现系数）。 非线性方程（自变量不显式出现）。 非线性齐次方程 常系数线性方程 变系数线性方程 有时可用积分变换（拉普拉斯变换、Fourier变换）求解。 例5－3 多孔催化剂中的粒内扩散－反应过程的解。球形催化剂内部，一级等温反应常微分方程模型： R=-kc，展开得到二阶、变系数、线性方程： 变量置换，令u=rc则 , 定义参数?=，单位为1/m。 边界条件： 特征方程： 待定常数须从边界条件得到。 是稳态催化反应时的粒内反应物浓度分布。 催化反应速率，可用两种方法来求得。 用局部速率的体积分来计算, 用穿过外表面的扩散通量来计算： 无内扩散阻力, ，, 反应速率,催化剂的有效因子： 催化反应工程中非常重要的性能指标。 Thiele模数: 例5－5 两平板间的传热问题的求解。 流体温度的渐近分布（x→∞时解的极限）: A为正比于q的常数, t(y)仅为坐标y的函数， 积分t(y)的微分方程两次， 在y= ??的右边界的条件： ， 在y=-??的左边界条件： , 待定常数C2与y无关，但可能与x有关。与渐近解中C合并为一个常数，可从第4个边界条件求得： ，右第一项抵消qx， 右第三项＝，故 可知，此合并的C确实是与x无关的常数。 例5－2 套管换热器 逆流流动 动态分布参数微分方程 边界条件：T1=T10， 在 x=0； T2=T20， 在 x=L处 初始条件：T1(x,t)= T10， T2(x,t)= T20， 在 t=0时 本例为边界值问题，求解不是特别容易。 解析求解的概念和方法可以更简单的例子来说明。 管外温度T20为定值，则简化为一阶线性偏微分方程： 边界条件为：T1=T10， 在 x=0 初始条件为：T1(x,t)= T00， 在 t=0 一阶拟线性偏微分方程的一般形式为 若f，g与u无关，且h(t,x,u)=r(t,x)u+s(t,x)，则上式称为一阶线性偏微分方程。