偏微分方程数值解课件-2.ppt

1、偏微分方程数值解偏微分方程数值解Email:JPhone:02985583997本章要求本章要求 教学目的教学目的 讲解：偏微分方程离散格式及求解的一般过程 教学要求教学要求E熟记 一阶及二阶偏微分方程的离散格式；C 精通 用EXCEL迭代对偏微分方程求解；探索 用两数组交替更新的办法进行编程求解；F 延伸 对化学反应工程中物理场的模拟进行尝试。教学重点教学重点各种偏微分方程的离散与求解EXCEL 循环迭代问题 教学难点教学难点 特殊边界条件的引入与应用Email:JPhone:029855839975.1 偏微分方程简介&偏微分方程如果一个微分方程中出现多元函数的偏导数，或者说如果未知函数和

2、几个变量有关，而且方程中出现未知函数对几个变量的导数，那么这种微分方程就是偏微分方程。在化工或化学动态模拟方程中，常常有一个自变量是时间，其它的自变量为空间位置。如果只考虑一维空间，则只有两个自变量；如果考虑两维空间，则有3个自变量。许多化工过程均是通过对偏微分方程的求解进行工艺参数的确定或数值模拟。Email:JPhone:029855839975.1 偏微分方程简介&偏微分方程的分类22222()()()()()()()0uuuuuabcdefugx yxyxy 线性微分方程 Linear partial differencial equation拟线性微分方程 Quasilinear p

3、artial differencial equation非线性微分方程 Nonlinear partial differencial equation,x y 1,/,nx yuxy ,/,nx yuxy Email:JPhone:029855839975.1 偏微分方程简介&数学上的分类：椭圆方程 Elliptic抛物线方程 Parabolic双曲线方程 Hyperbolic&物理实际问题的归类：波动方程(双曲型）一维弦振动模型：热传导方程（抛物线型）一维线性热传导方程拉普拉斯方程（椭圆型）稳态静电场或稳态温度分布场）240bac 240bac 240bac 22222uutx 22uutx

4、 22220uuxyEmail:JPhone:029855839975.1 微分方程的求解思路&求微分方程数值解的一般步骤：Step1区域剖分：首先按一定规则将整个定义域分成若干小块Step2微分方程离散：构造离散点或片的函数值递推公式或方程Step3初始、边界条件离散：根据递推公式，将初值或边界值离散化，补充方程，启动递推运算Step4 数值解计算：求解离散系统问题&微分方程的定解问题 离散系统的求解问题Email:JPhone:029855839975.2 离散化公式&将自变量在时间和空间上以一定的间隔进行离散化，则应变量就变成了这些离散变量的函数。&一阶偏导的离散化公式一般采用欧拉公式表

5、示&有时为了保证系统和稳定性，对时间的差分往往采用向后公式,(,)ni j kt n t x i x yj y z k zuu t x y z 1,1,1,1,nni j ki j kt n t x i x yj y z k znnij ki j kt n t x i x yj y z k znni jki j kt n t x i x yj y z k znni j ki j kt n t x i x yj y z k zuuuttuuuxxuuuyyuuuzx 1,(1),nni j ki j ktnt x i x yj y z k zuuutt Email:JPhone:02985583

6、9975.2 离散化公式P对于二阶偏导，我们可以通过对泰勒展开式处理技术得到下面离散化计算公式：112,22,21,1,22,2,1,1,22,222()2()2()nnni j ki j ki j kt n t x i x yj y z k znnnij ki j kij kt n t x i x yj y z k znnni jki j kij kt n t x i x yj y z k zt n tuuuuttuuuuxxuuuuyyuz ,1,12,2()nnni j ki j ki j kx i x yj y z k zuuuz Email:JPhone:029855839975.2

7、 离散化公式推导&将uk+1在uk处按二阶泰勒式展开：&将uk-1在uk处按二阶泰勒式展开：&二式相加得：22312()2!kkkkuuhuuhO hxx 22312()2!kkkkuuhuuhO hxx 211222()kkkuuuuxx Email:JPhone:029855839975.3几种常见偏微分方程的离散化计算&1 1、波动方程波动方程 其中：为初值条件 为边值条件 当该波动方程只提供初值条件时，称此方程为波动方程的初值问题，二者均提供时称为波动方程的混合问题。2222200012(,)(),()(),()ttxx luuaf x ttxuuxxtut ut 00012(),()

8、(),()ttxx tuuxxtutut Email:JPhone:029855839975.3.1 波动方程求解对于初值问题，是已知t=0时，u与 依赖于x的函数形式，求解不同位置，不同时刻的u值。而 u是定义在 的二元函数，即上半平面的函数。对于混合问题除初值外，还有边值。是已知初值及x=0及x=l 时u依赖于t的函数，求解不同位置x，不同时刻的u值。此时u是定义在 的带形区域上的二元函数。ut 0,tx xt 0a)初值问题初值问题tx0lb)混合问题混合问题Email:JPhone:029855839975.3.1 波动方程求解2222200012(,)(),()(),()ttxx l

9、uuaf x ttxuuxxtut ut 1121122(i1,2,m-1)22(,)()()(1,2,)nnnnnniiiiiiuuuuuuaf x ttxn 22212221211222()()()(22)()(,)()()()nnnnniiiiitttuauauauuxf x txxx 100012(),()(i1,2,m)(),()(n1,2,)iiinnmuuuj xixtun tun t 方程离散化整理可得：边界条件初始条件离散化xxinniu1niu 1niu Email:JPhone:029855839975.3.1 波动方程求解&例例5.1:5.1:用数值法求解下面偏微分方程

10、。0012()3150,30,2501,0WnWjxttTtxTtttxx 110.01,0.1nnjjnnjjttttttxxx 110.020.680.3nnnjWjjtTtt 112()3nnnnjjjjjWntttttTx 此微分方程，是在不考虑流体本身热传导时的套管传热微分方程.由计算结果可知，当计算的时间序列进行到7272时，传热过程已达到稳态，各点上的温度已不随时间的增加而改变。如果改变套管长度或传热系数，则达到稳态的时间亦会改变。Email:JPhone:029855839975.3.2 一维流动热传导方程 与波动方程的情形类似，用差商近似代替偏商，可以得到一维流动传热传导方程

11、的混合问题的差分方程，以其解作为流动传热传导方程的近似解。222001(,)00 00 0 tx lxuuuabf u t(xl,t)txxu(x),(xl)u(t)xu(t)0(t)12111201012(,)()()(i1 2)0 (0,1,2,)(nnnnnnniiiiiiiinnmmnuuuuuuuabf i x n ttxxui x,muunxun )(0,1,2,)tn 1,21,1,22,2,1,1,22,22,2()2()nni j ki j kt n t x i x yj y z k znnnij ki j kij kt n t x i x yj y z k znnni jk

12、i j kij kt n t x i x yj y z k zt n t x i x yj y z kuuuttuuuuxxuuuuyyuz ,1,122()nnni j ki j ki j kzuuuz 2 2、一维流动热传导方程的混合问题、一维流动热传导方程的混合问题Email:JPhone:02985583997将上式进行处理得到：该式是显式格式。只要保证式中各项系数大于零，一般情况下是稳定的，可以获得稳定的解。分析上式可以发现，当为了提高数值精度取适当小的x 时，最有可能小于零的系数是 uin的系数，若要保证此项系数大于零，此时t必须相应地更小,会导致计算量将大大增加，这是显式格式的缺

13、点，为了克服此缺点，下面提出一种隐式格式：偏微分方程在 点上进行离散化，且对时间的偏微分采用向后欧拉公式得到原偏微分方程的离散化公式：111111211122(,(1)()nnnnnnniiiiiiiuuuuuuuabf i x nttxx 122211222(,)()(12)()()()nnnniiiitttttutf i x n tabuabuauxxxxx (,(1)i x nt 5.3.2 一维流动热传导方程Email:JPhone:02985583997 从图5-3中可见要由初值及边界条件一排一排推上去是不行的，需解线性方程组,同时添上二边界条件：正好共有m+2个方程，同时有m+2个

14、变量，就能解出n+1排上各点值。这样，每解一个线性方程组，就可以往上推算一排点的u值，虽然引入了方程组的求解，有可能增加计算量，但由于隐式格式无条件稳定，t的取法与x 无关，可以少计算许多排节点上的u 值，相应于显式格式来说，最终反而节省了计算量。111011(1),nnnmmuntuu 5.3.2 一维流动热传导方程Email:JPhone:02985583997&例5.2 考虑纵向导热的套管换热器内管各点温度分布微分方程：解：首先根据前面的知识，将所求 的方程离散化：代入微分方程并化简得：分析上式可知，如果知道了某一时刻的各点t，（，（j=0,1,2.10,11），），就可以求下一时刻的各

15、点温度值t(j=1,2.10),现在已经知道了零时刻管内各点的温度分布及入口处在任何时刻的温度，如想求下一时刻的温度值，根据上面的离散化计算公式，还需知道在j=11处的温度，这个温度可利用给定的边界条件离散化求得：有了以上各式，上面的微分方程就可以求解了。220012()0.0013150,30,3001,0WnWjxtttTtxxTtttxx 11211222()0.01,0.1nnjjnnjjnnnjjjttttttxxttttxxx 1110.020.0010.6780.301nnnnjWjjjtTttt 110nnjjxtttxx 5.3.2 一维流动热传导方程Email:JPhone

16、:029855839975.3.3 稳态导热/扩散方程&3 3、稳态导热、稳态导热/扩散方程扩散方程 在化工导热及扩散过程中，没有物流的流动，仅靠导热及扩散进行热量及质量的传递。如果此时系统达到稳定状态，也就是说系统中每一个控制单元的各项性质如温度、浓度等不再随时间的改变而改变，系统中的各种性质只与其所处的位置有关，利用化工知识，我们可以得到下面二维、三维的稳态导热或扩散偏微分方程：二维：三维：二维的稳态导热或扩散偏微分方程又称调和方程。22220uuxy2222220uuuxyzuu un ()uun 常见有三种边界条件：第一类边界条件：第二类边界条件：第三类边界条件：Email:JPhon

17、e:02985583997&离散化公式：&取 ，经化简得：&外节点（边界节点）和内节点&求解方法划分网格建立节点离散方程迭代求解（或解稀疏方程组）1,1,1,122220()()iji jiji ji ji juuuuuuxyxy ,1,11,11()4i jiji jiji juuuuu xy求解区域N 节点边界五点格式示意图5.3.3稳态导热/扩散方程求解Email:JPhone:029855839975.3.3稳态导热/扩散方程求解 常用的3种迭代格式：(1)同步迭代：(2)异步迭代：(3)超松弛迭代：当计算范围R 为 矩阵区域，x方向m等分，y方向n等分，最佳松弛因子为：由数学知识可知

18、，用这些迭代法求解上面的偏微分方程均收敛。(1)()()()(),1,11,11()4kkkkki jiji jiji juuuuu (1)()()(1)(1),1,11,11()4kkkkki jiji jiji juuuuu ()()(1)(1)1,11,1(1)(),1()4(1)kkkkiji jiji jkki ji juuuuuuwuw u (,)axb cyd 22coscos112wmn 紧凑迭代Email:JPhone:029855839975.3.3 稳态导热/扩散方程求解例例5.35.3：处于传热平衡状态的某保温，假设其形状为长方体，在x，y两个方向上存在热传导，且导热系

19、数相等，已知边界温度分布如下图所示：解：取某一微元进行能量衡算，由于已达传热平衡状态，故可得：传导入热量-传导出热量=0 1xy10(1,1)24(,)30 27001,0t x yyyx 1(,)30 27001,0t x yxxy 2(,)30001,0tx yyx3(,)30001,1t x yxy()yttyyx zy ()xttxxy zx ytx zy xty zx xyz22220 xyttxyEmail:JPhone:029855839975.3.3 稳态导热/扩散方程求解&Microsoft Excel 迭代计算公式中的循环引用在“工具工具”菜单上，单击“选项”，再单击“重新

20、计算”选项卡。选中“迭代”复选框。若要设置 Microsoft Excel 进行重新计算的最大次数，请在“最多迭代次数”框中键入迭代次数。迭代次数越高，Excel 用于计算工作表的时间越多。若要设置两次迭代结果之间可以接受的最大误差，请在“最大误差”框中键入所需的数值。数值越小，结果越精确，Excel 用于计算工作表的时间也越多。Email:JPhone:029855839975.4 吸附床传热传质模型中偏微分方程求解实例&5.4.1 5.4.1 基本设定及假设基本设定及假设 1.吸附器结构参数的设定 上图所示的是套筒式吸附器，该吸附器的有效长度为L，其有效内径为D,环隙宽度为，吸附器壁厚为b

21、。导热流体通过环隙将热量传入或传出吸附器，吸附质通过吸附器上端的小管进入或离开吸附器。吸附器结构示意图吸附器结构示意图 DL热流体Email:JPhone:029855839975.4.1 基本设定及假设2.吸附床外流体传热的一些基本假设：1).忽略流体在环隙宽度上的温度梯度；2).忽略热损失；3).忽略吸附器壁厚b上的温度梯度，用集中参数法求取吸附器壁面温度。.吸附床内传热传质的一些基本假设：1).吸附床内的吸附质气体处于气滞状态；2).忽略蒸发器、冷凝器和吸附床之间的压力差；3).吸附床内各计算微元内达到吸附平衡。吸附量可利用回归方程计算；4.吸附热利用微分吸附热，随吸附量和吸附温度的改变

22、而改变；比热采用有效比热，亦随温度改变，但在计算微元内，可认为是常数；5.床层活性炭导热系数采用当量导热系数，可由实验测量得到。Email:JPhone:029855839975.4.2流体传热模型的建立在轴方向上取一环隙微元，作能量分析如下：1.流体通过流动流入环隙微元的能量为 2.流体通过流动流出环隙微元的能量3.流体热传导在x 处的热量导入 7 总能量平衡方程fffTTTxxxxx 流体传热微元模型流体传热微元模型 infffpffquSCT ()fou tfffpffTquSCTxx outqinqxqfxffTqSx 其中：f 流体的密度 uf 环隙的流体速度，Sf 环隙的横截面积，

23、Cpf流体的比热。4.流体热传导在x+x处的热量导入 5.微元体传递给吸附床的热量 qt 6.微元体内的能量变化率 为流体的横截面积xxq(/)()fFxxffTTxqxSxx ()tffwqhxDTT cq()fpfffcCSx Tqt 221()ffffffwffffTuThDTTTtxSx fffpfC fSEmail:JPhone:029855839975.4.3 吸附床内吸附剂传热传质模型的建立&吸附床内发生着热量和质量的传递，但质量的传递是建立在热量传递基础上的，故只要建立热量传递方程，就可以根据平衡吸附量方程求出各处的吸附量。吸附床内的热量传递主要以热传导为主，既有经向的热传导，

24、也有轴向的热传导，为了便于建模分析，选取如图所示的吸附床微元体，进行衡算：2bxaTqrrx x+x x xrr+r吸附床内传热传质微元体1.轴向导入热量:2.轴向导出热量3.径向导入热量 4.径向导出热量 5.微元体内的能量变化率 其中 为吸附床层内的有效比热。6.总能量平衡方程 222()bbxxaTTqrrxxx 2braTqrxr 222()()bbdraTTqrrxrrr 2bceffaTqrrx Ct ()effpapbPmCCmCHT 222211bbbbeffTTTTtrrxraeffaeffC Email:JPhone:029855839975.4.4 吸附器内/外无量纲化方

25、程吸附器内/外无量纲化方程2,021()()fwwwfwwbswswshDThDTTTTTxSSx ,sswpsssC 221()ffffffwffffTuThDTTTtxSx fffpfC 222211bbbbeffTTTTtrrxraeffaeffC *2000000,ffwbininintxrxrtLRLTTTTTTTTTTTT 无量纲化处理无量纲化处理Email:JPhone:029855839975.4.4 吸附器内吸附器内/外无量纲化方程外无量纲化方程2*2()PeKA Bitxx 20*2()sswsfswsfDrBiBiBiBitx 22222,fffffbbefffbbfff

26、fwsswsfsswswuC LhDLSDrPeBiKASShDLhDLLDrBiBiCSSR 22*2*2*bDrCCtxrrr整理可得：其中：Email:JPhone:029855839975.4.4 吸附器内吸附器内/外无量纲化方程外无量纲化方程*(,0)0(,0)0(,0)0 x tx tx rt*0*1(0,)1 0(1,)(1,)0 xxxtxxtxtx *1000 0 0 xxrxxr*0*()/bwrRhRr 初始条件：边界条件Email:JPhone:029855839975.4.6 模型的离散化,1,*,1,*21,1,*2*22i ki ki kikikiki kttxx

27、xx ,1,*1,1,*2*2,1,*1,1,*2*2,1,*,1,1,*2*222i ki kikiki ki j ki j kij kij ki j ki j ki jki jki jki j kttxxttxrrrrr ,00,1,00,1,1,0,0,0,1,1,1,*0,01,0000()ikn knkikkn knki ji jikikn j knj kbi mki m ki ki m kwikhrR 偏导数的差分离散化偏导数的差分离散化边界条件离散化处理边界条件离散化处理Email:JPhone:029855839975.4.6 模型的离散化,111,2,31,4.,111,2,

28、31,4,5,111,2,31,4,1,5,1,i kiki kiki ki kiki kiki m ki ki j kij ki j kij ki jki jkAAAABBBBBCCCCC *12*2*2*34*2*1231*2*2*451*2*2*22,1()2,1,221swsfssssfswssbbtPettAAKA BitxxxPettAAKA BitxxBiBitttBBBBDrDrxDrxBitBittBBCDrDrxtCtCDrxDr *31*2*2*45*2*2*2,bbbbCtCCrDrjrCtCtCtCCDrrDrjrDrr 离散化离散化方方 程程离散化系数离散化系数Em

29、ail:JPhone:029855839975.4.7 模型求解参数及求解区域初始化参数及求解区域初始化偏微分方程系数计算偏微分方程系数计算离散化方程系数计算离散化方程系数计算收敛性判断收敛性判断各偏导系数各偏导系数0方程组迭代求解方程组迭代求解输出结果输出结果YN缩小时间步长缩小时间步长重新划分网格重新划分网格迭代循环迭代循环第二节第二节 现代汉语词汇的构成现代汉语词汇的构成一、词、固定词组一、词、固定词组二、传承词、新造词二、传承词、新造词三、口语词、书面语词、方言词、三、口语词、书面语词、方言词、外来词、行业词外来词、行业词一、词、固定词组一、词、固定词组（一）词（一）词 词是最小的音义

30、结合定型的词是最小的音义结合定型的造句造句单位。单位。如：黑板如：黑板 语素是最小的音义结合定型的语素是最小的音义结合定型的构词构词单位。单位。如：如：技技、术、人、葡萄、术、人、葡萄 没有特定语言环境时，没有特定语言环境时，“人人”既是一个语素，既是一个语素，又是一个词，但在具体的语言环境中，又是一个词，但在具体的语言环境中，“人人”只能是语素（如只能是语素（如“人民、人气人民、人气”），或只能），或只能是词（如是词（如“人是要有一点精神的人是要有一点精神的”）。）。（二）固定词组（二）固定词组 专名一般是由几个词构成的，但跟词一样专名一般是由几个词构成的，但跟词一样具有音义定型性，不能扩展

31、，使用时也只具有音义定型性，不能扩展，使用时也只作为一个造句单位，作句子的一个成分。作为一个造句单位，作句子的一个成分。熟语是由词构成而功能上又相当于一个词熟语是由词构成而功能上又相当于一个词的具有结构凝固性、意义整体性、社会习的具有结构凝固性、意义整体性、社会习用性的固定词组，包括成语、惯用语、歇用性的固定词组，包括成语、惯用语、歇后语等。后语等。1.成语成语 成语是以成语是以“四言四言”为主要形式、具有结构凝固性、为主要形式、具有结构凝固性、意义整体性、社会习用性和书面语色彩的一种熟语。意义整体性、社会习用性和书面语色彩的一种熟语。发号施令发号施令 胡思乱想胡思乱想 切磋琢磨切磋琢磨 孤掌

32、难鸣孤掌难鸣 枯木逢春枯木逢春 夜郎自大夜郎自大 墨守成规墨守成规 固执己见固执己见 包罗万象包罗万象 浩然正气浩然正气 一面之词一面之词 摇摇欲坠摇摇欲坠 殃及池鱼殃及池鱼 逍遥法外逍遥法外 重于泰山重于泰山 请君入瓮请君入瓮 调虎离山调虎离山 引狼入室引狼入室 解甲归田解甲归田 画蛇添足画蛇添足 破釜沉舟破釜沉舟 2.惯用语惯用语 惯用语是以惯用语是以“三言三言”为主要形式、结构较为定型、意为主要形式、结构较为定型、意义有所引申、具有浓郁口语色彩的一种熟语。义有所引申、具有浓郁口语色彩的一种熟语。多为动宾结构：多为动宾结构：拍马屁拍马屁 穿小鞋穿小鞋 捅娄子捅娄子 出洋相出洋相 钻空子钻空

33、子 吹牛皮吹牛皮 唱高调唱高调 泼冷水泼冷水 碰钉子碰钉子 和稀泥和稀泥 也有偏正结构：铁饭碗也有偏正结构：铁饭碗 敲门砖敲门砖 马后炮马后炮 闭门羹闭门羹 定型性比成语弱，有时可插入一些词语。定型性比成语弱，有时可插入一些词语。意义的凝固性强。背黑锅意义的凝固性强。背黑锅背着黑色的锅子。背着黑色的锅子。多有比喻义、引申义，不同于一般的三音节词语。多有比喻义、引申义，不同于一般的三音节词语。踢皮球踢皮球踢足球踢足球3.歇后语歇后语 歇后语由两部分组成，前一部分像谜面，后一部分歇后语由两部分组成，前一部分像谜面，后一部分像谜底，两部分之间有明显停顿，整体表示跟后一部分像谜底，两部分之间有明显停顿

34、，整体表示跟后一部分相关的意义。相关的意义。比喻型：周瑜打黄盖比喻型：周瑜打黄盖一个愿打，一个愿挨一个愿打，一个愿挨 骑着毛驴看唱本骑着毛驴看唱本走着瞧走着瞧 徐庶进曹营徐庶进曹营一言不发一言不发 哑巴吃饺子哑巴吃饺子心中有数心中有数 泥菩萨过河泥菩萨过河自身难保自身难保双关型双关型：咸菜缸里的秤砣：咸菜缸里的秤砣一盐（言）难进（尽）一盐（言）难进（尽）猪八戒的脊梁猪八戒的脊梁悟能（无能）之背（辈）悟能（无能）之背（辈）南天门挂镜子南天门挂镜子照（赵）云照（赵）云 孔夫子搬家孔夫子搬家尽是书（输）尽是书（输）外甥打灯笼外甥打灯笼照舅（旧）照舅（旧）二、传承词、新造词二、传承词、新造词 传承词：

35、传承词：在现代汉语形成之前产生并在现代汉语中使用的词。在现代汉语形成之前产生并在现代汉语中使用的词。新造词：新造词：在现代汉语形成之后产生的词。在现代汉语形成之后产生的词。现代汉语词汇现代汉语词汇传承词传承词新造词新造词一、什么是词？它跟语素、一般词组的区别在哪里？一、什么是词？它跟语素、一般词组的区别在哪里？词是最小的音义结合定型的造句单位。词是最小的音义结合定型的造句单位。语素是最小的音义结合定型的构词单位。语素也语素是最小的音义结合定型的构词单位。语素也具有音义定型性，但它不用于造句，不是造句单位。具有音义定型性，但它不用于造句，不是造句单位。语素的功能在于造词。语素的功能在于造词。一般词组也可以用于造句，是造句单位。但一般一般词组也可以用于造句，是造句单位。但一般词组缺少像词那样的音义定型性。一般词组都是由两词组缺少像词那样的音义定型性。一般词组都是由两个或更多个词共同组成的，它们之间的组合是临时性个或更多个词共同组成的，它们之间的组合是临时性的，它们可以在形式上进行扩展（或缩略）而意义保的，它们可以在形式上进行扩展（或缩略）而意义保持不变。持不变。思考和练习二思考和练习二