计算机应用基-偏微分方程求解.pptVIP

固定床二维反应器模型 Z为反应管轴向距离，r为径向距离，方程系数如下： 6.4 偏微分方程在化工中的应用 通过反应计量关系获得各组分的分压： 温度初始边界条件： 6.4 偏微分方程在化工中的应用 质量初始边界条件 径向边界条件 6.4 偏微分方程在化工中的应用 微分方程组的差分格式离散化 一 隐式差分（Crank-Nicholson) M-1, n M, n+1 M+1, n M, n M+1, n+1 M-1, n+1 M, n+1/2 ?z ?r 采用差分代替微分，离散化进行迭代。 如图，任一节点的径向与轴向位置可表示为： 6.4 偏微分方程在化工中的应用 对于点（M,n+1/2)的各项导数都可以用该处周围六个节点的函数值计算的差商表示，如温度可表示为： 将上三式代入二维模型的物料衡算方程 6.4 偏微分方程在化工中的应用 类似的可求得浓度的差分方程 6.4 偏微分方程在化工中的应用 上两式中(rA)m,n+1/2是指温度为Tm,n+1/2与浓度cAm,n+1/2条件下的反应速率，有： 6.4 偏微分方程在化工中的应用 管中心和管壁处的边值条件也可用差分方程表示。管中心： 6.4 偏微分方程在化工中的应用 管壁处： 6.4 偏微分方程在化工中的应用 迭代过程： 假设所有节点处浓度和温度的初始值C0Am,n和T0m,n 计算出CAm,n＋1、2和Tm,n＋1、2 ，代入动力学方程求得（-rA)m,n+1/2 将（-rA)m,n+1/2代入上述方程，得到各节点处浓度和温度的线性方程组，方程的系数矩阵均为对三角矩阵，用追赶法求解浓度方程和温度方程，得到各节点处浓度和温度的新值CAm,n和Tm,n 若计算的新值与上次迭代的值（或初始值）之偏差小于规定精度，计算结束，否则以新值为假定值，返回步骤（2）。 6.4 偏微分方程在化工中的应用 计算程序”PDEs2DS_Cranknicolson.m” 6.4 偏微分方程在化工中的应用 第六章 偏微分方程求解 6.1 偏微分方程组求解 6.2 二阶偏微分方程的数学描述 6.3 偏微分方程的求解界面应用举例 6.4 偏微分方程在化工中的应用 一 偏微分方程的分类 当A，B，C为常数时，称为拟线性偏微分方程，可分为三类： 椭圆型方程 抛物型方程 双曲型方程 6.1 偏微分方程组求解 二 偏微分方程边界条件： （1）Dirichlet 边界条件 hu=r 也称为第一类边界条件，对于偏微分方程组，Dirichlet边界条件为 （2）Neumann边界条件 6.1 偏微分方程组求解 也称为第三类边界条件；当q=0时，则变为第二类边界条件。对于偏微分方程组，Neumann边界条件为： 其中 n为边界外法向单位向量，g, q, h, r为在边界上定义的函数 （3）混合边界条件 6.1 偏微分方程组求解 6.1 偏微分方程组求解 三 偏微分方程数值解法 有限差分法 正交配置法 MOL法 有限元法 6.1 偏微分方程组求解 四 采用pdepe( )函数求解一维偏微分方程 边界条件的函数描述： 6.1 偏微分方程组求解 【例6-1】 6.1 偏微分方程组求解 6.1 偏微分方程组求解 6.1 偏微分方程组求解 function [c, f, s]=c7mpde(x, t, u, du) c=[1?; 1]?; y=u(1)-u(2)?; F=exp(5.73*y)-exp(-11.46*y)?; s=F*[-1; 1] f=[0.024*du(1); 0.17*du(2)]; 偏微分方程求解程序 “c7mpde” 6.1 偏微分方程组求解 function [pa, qa, pb, qb]=c7mpbc(xa, ua, xb, ub, t) pa=[0; ua(2)]; qa=[1; 0]; pb=[ub(1)-1; 0]; qb=[0; 1]; 边界条件程序”c7mbc.m” function u0=c7mpic(x) u0=[1; 0]; x=0: 0.05: 1; t=0: 0.05: 2; m=0; sol=pdepe(m, @c7mpde, @c7mpic, @c7mpbc, x, t); surf(x, t, sol(?:,?:, 1)) 6.1 偏微分方程组求解 求解函数：ME_5_1.m, 结果如下 6.2 二阶偏微分方程的求解 一 椭圆型偏微分方程 6.2 二阶偏微分方程的求解 adaptmesh 和assempde函数用于求解椭圆型偏微分方程的解，调用格式如下： [u, p, e, t]=adapmesh(g,b,c,a,f) g: 求解几何区域