偏微分方程数值解2省名师优质课获奖课件市赛课一等奖课件.ppt

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,上一页,下一页,回目录,休 息,Email:,Jansweili,Phone:[phone],本幻灯片资料仅供参考，不能作为科学依据，如有不当之处，请参考专业资料。谢谢,偏微分方程数值解,第1页,本章要求,教学目标,讲解：,偏微分方程离散格式及求解普通过程,教学要求,熟记 一阶及二阶偏微分方程离散格式；,精通 用,EXCEL,迭代对偏微分方程求解；,探索,用两数组交替更新方法进行编程求解,；,延伸 对化学反应工程中物理场模拟进行尝试。,教学重点,各种偏微分方程离散与求解,EXCEL,循环迭代问题,教学难点,特殊边界条件引入与应用,第2页,5.1,偏微分方程介绍,偏微分方程,假如一个微分方程中出现多元函数偏导数，或者说假如未知函数和几个变量相关，而且方程中出现未知函数对几个变量导数，那么这种微分方程就是偏微分方程。,在化工或化学动态模拟方程中，经常有一个自变量是时间，其它自变量为空间位置。假如只考虑一维空间，则只有两个自变量；假如考虑两维空间，则有3个自变量。许多化工过程均是经过对偏微分方程求解进行工艺参数确实定或数值模拟。,第3页,5.1,偏微分方程介绍,偏微分方程分类,线性微分方程,Linear partial differencial equation,拟线性微分方程,Quasilinear partial differencial equation,非线性微分方程,Nonlinear partial differencial equation,第4页,5.1,偏微分方程介绍,数学上分类：,椭圆方程,Elliptic,抛物线方程,Parabolic,双曲线方程,Hyperbolic,物理实际问题归类：,波动方程,(,双曲型）一维弦振动模型：,热传导方程（抛物线型）一维线性热传导方程,拉普拉斯方程（椭圆型）稳态静电场或稳态温度分布场）,第5页,5.1,微分方程求解思绪,求微分方程数值解普通步骤：,Step1,区域剖分,：首先按一定规则将整个定义域分成若干小块,Step2,微分方程离散,：结构离散点或片函数值递推公式或方程,Step3,初始、边界条件离散,：依据递推公式，将初值或边界值离散化，补充方程，开启递推运算,Step4,数值解计算,：求解离散系统问题,微分方程定解问题 离散系统求解问题,第6页,5.2,离散化公式,将自变量在时间和空间上以一定间隔进行离散化，则,应变量就变成了这些离散变量函数,。,一阶偏导离散化公式,普通采取欧拉公式表示,有时为了确保系统和稳定性，,对时间差分往往采取向后公式,第7页,5.2,离散化公式,对于二阶偏导，我们能够经过对,泰勒展开式,处理技术得到下面离散化计算公式：,第8页,5.2,离散化公式推导,将,u,k+1,在,u,k,处按二阶泰勒式展开：,将,u,k-1,在,u,k,处按二阶泰勒式展开：,二式相加得：,第9页,5.3,几个常见偏微分方程离散化计算,1,、,波动方程,其中：为初值条件,为边值条件,当该波动方程只提供初值条件时，称此方程为波动方程初值问题，二者均提供时称为波动方程混合问题,。,第10页,5.3.1,波动方程求解,对于初值问题，是已知,t,=0,时，,u,与 依赖于,x,函数形式，求解不一样位置，不一样时刻,u,值。而,u,是定义在 二元函数，即上半平面函数。,对于混合问题除初值外，还有边值。是已知初值及,x=0,及,x,=,l,时,u,依赖于,t,函数，求解不一样位置,x,，,不一样时刻,u,值。此时,u,是定义在 带形区域上二元函数。,x,t,0,a),初值问题,t,x,0,l,b),混合问题,第11页,5.3.1,波动方程求解,方程离散化,整理可得：,边界条件,初始条件,离散化,x,x,i,n,第12页,5.3.1,波动方程求解,例,5.1:,用数值法求解下面偏微分方程。,此微分方程，是在不考虑流体本身热传导时套管传热微分方程,.,由计算结果可知，当计算时间序列进行到,72,时，传热过程已到达稳态，各点上温度已不随时间增加而改变。假如改变套管长度或传热系数，则到达稳态时间亦会改变。,EXCEL,第13页,5.3.2,一维流动热传导方程,与波动方程情形类似，用差商近似代替偏商，能够得到一维流动传热传导方程混合问题差分方程，以其解作为流动传热传导方程近似解。,2、一维流动热传导方程混合问题,离,散,化,第14页,将上式进行处理得到：,该式是显式格式。只要确保式中各项系数大于零，普通情况下是稳定，能够取得稳定解。,分析上式能够发觉，当为了提升数值精度取适当小,x,时，最有可能小于零系数是,u,i,n,系数，若要确保此项系数大于零，此时,t,必须对应地更小,会造成计算量将大大增加，这是显式格式缺点，为了克服此缺点，下面提出一个隐式格式：,偏微分方程在 点上进行离散化，且对时间偏微分采取向后欧拉公式得到原偏微分方程离散化公式：,5.3.2,一维流动热传导方程,第15页,从图5-3中可见要由初值及边界条件一排一排推上去是不行，需解线性方程组,同时添上二边界条件：,恰好共有,m,+2,个方程，同时有,m,+2,个变量，就能解出,n,+1,排上各点值。这么，每解一个线性方程组，就能够往上推算一排点,u,值，即使引入了方程组求解，有可能增加计算量，,但因为隐式格式无条件稳定,，,t,取法与,x,无关，能够少计算许多排节点上,u,值，对应于显式格式来说，最终反而节约了计算量。,5.3.2,一维流动热传导方程,第16页,例,5.2,考虑纵向导热套管换热器内管各点温度分布微分方程：,解：首先依据前面知识，将所求 方程离散化：,代入微分方程并化简得：,分析上式可知，假如知道了某一时刻各点,t,，（,j,=0,1,2.10,11），,就能够求下一时刻各点温度值,t,(,j,=1,2.10),现在已经知道了零时刻管内各点温度分布及入口处于任何时刻温度，如想求下一时刻温度值，依据上面离散化计算公式，还需知道在,j,=11,处温度，这个温度可利用给定边界条件离散化求得：,有了以上各式，上面微分方程就能够求解了,。,5.3.2,一维流动热传导方程,EXCEL,第17页,5.3.3,稳态导热,/,扩散方程,3,、稳态导热,/,扩散方程,在化工导热及扩散过程中，没有物流流动，仅靠导热及扩散进行热量及质量传递。假如此时系统到达稳定状态，也就是说系统中每一个控制单元各项性质如温度、浓度等不再随时间改变而改变，系统中各种性质只与其所处位置相关，利用化工知识，我们能够得到下面二维、三维稳态导热或扩散偏微分方程：,二维：,三维：,二维稳态导热或扩散偏微分,方程又称,调和方程,。,常见有三种边界条件：,第一类边界条件：,第二类边界条件：,第三类边界条件：,第18页,离散化公式：,取 ，经化简得：,外节点（边界节点）和内节点,求解方法,划分网格,建立节点离散方程,迭代求解（或解稀疏方程组）,x,y,求解区域,N,节点,边界,五点格式示意图,5.3.3,稳态导热/扩散方程求解,第19页,5.3.3,稳态导热/扩散方程求解,惯用,3,种迭代格式：,(1),同时迭代：,(2),异步迭代：,(3),超松弛迭代：,当计算范围,R,为,矩阵区域，,x,方向,m,等分，,y,方向,n,等分，,最正确松弛因子为：,由数学知识可知，用这些迭代法求解上面偏微分方程均收敛。,紧凑迭代,第20页,5.3.3,稳态导热/扩散方程求解,例5.3,：处于传热平衡状态某保温，假设其形状为长方体，在,x，y,两个方向上存在热传导，且导热系数相等，已知边界温度分布以下列图所表示：,解：取某一微元进行能量衡算，,因为已达传热平衡状态，故可得：,传导入热量,-,传导出热量,=0,1,x,y,1,0,(1,1),温度分布,x,y,z,第21页,5.3.3,稳态导热/扩散方程求解,Microsoft Excel,迭代计算公式中,循环引用,在“,工具,”菜单上，单击“,选项,”，再单击“,重新计算,”选项卡。,选中“,迭代,”复选框。,若要设置,Microsoft Excel,进行重新计算最大次数，请在“,最多迭代次数,”框中键入迭代次数。迭代次数越高，,Excel,用于计算工作表时间越多。,若要设置两次迭代结果之间能够接收最大误差，请在“,最大误差,”框中键入所需数值。数值越小，结果越准确，,Excel,用于计算工作表时间也越多。,第22页,5.4,吸附床传热传质模型中偏微分方程求解实例,5.4.1,基本设定及假设,1.,吸附器结构参数设定,上图所表示是套筒式吸附器，该吸附器有效长度为,L,，,其有效内径为,D,环隙宽度为,，,吸附器壁厚为,b,。,导热流体经过环隙将热量传入或传出吸附器，吸附质经过吸附器上端小管进入或离开吸附器。,吸附器结构示意图,D,L,热流体,第23页,5.4.1,基本设定及假设,2.吸附床外流体传热一些基本假设：,1).,忽略流体在环隙宽度,上温度梯度；,2).,忽略热损失；,3).,忽略吸附器壁厚,b,上温度梯度，用集中参数法求取吸附器壁面温度。,.吸附床内传热传质一些基本假设：,1),.,吸附床内吸附质气体处于气滞状态；,2).,忽略蒸发器、冷凝器和吸附床之间压力差；,3).,吸附床内各计算微元内到达吸附平衡。吸附量可利用回归方程计算；,4.,吸附热利用微分吸附热，随吸附量和吸附温度改变而改变；比热采取有效比热，亦随温度改变，但在计算微元内，可认为是常数；,5.,床层活性炭导热系数采取当量导热系数，可由试验测量得到。,第24页,5.4.2流体传热模型建立,在轴方向上取一环隙微元，作能量分析以下：,1.,流体经过流动流入环隙微元能量为,2.流体经过流动流出环隙微元能量,3.流体热传导在,x,处热量导入,7 总能量平衡方程,流体传热微元模型,其中：,f,流体密度,u,f,环隙流体速度，,S,f,环隙横截面积，,C,pf,流体比热。,4.流体热传导在,x+,x,处热量导入,5.微元体传递给吸附床热量,q,t,6.,微元体内能量改变率,为流体横截面积,第25页,5,.4.3,吸附床内吸附剂传热传质模型建立,吸附床内发生着热量和质量传递，但质量传递是建立在热量传递基础上，故只要建立热量传递方程，就能够依据平衡吸附量方程求出各处吸附量。吸附床内热量传递主要以热传导为主，现有经向热传导，也有轴向热传导，为了便于建模分析，选取如图所表示吸附床微元体，进行衡算：,x+,x,x,x,r,r,+,r,吸附床内传热传质微元体,1.轴向导入热量,:,2,.,轴向导出热量,3.径向导入热量,4,.,径向导出热量,5.微元体内能量改变率,其中,为吸附床层内有效比热。,6.总能量平衡方程,第26页,5.4.4,吸附器内,/,外无量纲化方程,吸附器内,/,外无量纲化方程,无量纲化处理,第27页,5.4.4,吸附器内,/,外无量纲化方程,整理可得：,其中：,第28页,5.4.4,吸附器内,/,外无量纲化方程,初始条件：,边界条件,第29页,5.4.6,模型离散化,偏导数差分离散化,边界条件离散化处理,第30页,5.4.6,模型离散化,离散化,方 程,离散化系数,第31页,5.4.7,模型求解,参数及求解区域初始化,偏微分方程系数计算,离散化方程系数计算,收敛性判断,各偏导系数,0,方程组迭代求解,输出结果,Y,N,缩小时间步长,重新划分网格,迭代循环,第32页,