计算方法 偏微分方程数值解演示课件.pptVIP

第五章 偏微分方程数值解 Numerical Methods for Partial Differential Equations 5.1 偏微分方程简介 5.2 离散化公式 5.3 几种常见偏微分方程的离散化计算 5.4吸附床传热传质模型中偏微分方程求解 本章要求 教学目的 讲解： 偏微分方程离散格式及求解的一般过程 教学要求 熟记 一阶及二阶偏微分方程的离散格式； 精通 用EXCEL迭代对偏微分方程求解； 探索 用两数组交替更新的办法进行编程求解； 延伸 对化学反应工程中物理场的模拟进行尝试。 教学重点 各种偏微分方程的离散与求解 EXCEL 循环迭代问题 教学难点 特殊边界条件的引入与应用 5. 1 偏微分方程简介 偏微分方程 如果一个微分方程中出现多元函数的偏导数，或者说如果未知函数和几个变量有关，而且方程中出现未知函数对几个变量的导数，那么这种微分方程就是偏微分方程。 在化工或化学动态模拟方程中，常常有一个自变量是时间，其它的自变量为空间位置。如果只考虑一维空间，则只有两个自变量；如果考虑两维空间，则有3个自变量。 许多化工过程均是通过对偏微分方程的求解进行工艺参数的确定或数值模拟。 5.1 偏微分方程简介 偏微分方程的分类 ?线性微分方程 Linear partial differencial equation ?拟线性微分方程 Quasilinear partial differencial equation ?非线性微分方程 Nonlinear partial differencial equation 5.1 偏微分方程简介 数学上的分类： 椭圆方程 Elliptic 抛物线方程 Parabolic 双曲线方程 Hyperbolic 物理实际问题的归类： 波动方程(双曲型）一维弦振动模型： 热传导方程（抛物线型）一维线性热传导方程 拉普拉斯方程（椭圆型）稳态静电场或稳态温度分布场） 5.1 微分方程的求解思路 求微分方程数值解的一般步骤： Step1区域剖分：首先按一定规则将整个定义域分成若干小块 Step2微分方程离散：构造离散点或片的函数值递推公式或方程 Step3初始、边界条件离散：根据递推公式，将初值或边界值离散化，补充方程，启动递推运算 Step4 数值解计算：求解离散系统问题 微分方程的定解问题 离散系统的求解问题 5.2 离散化公式 将自变量在时间和空间上以一定的间隔进行离散化，则应变量就变成了这些离散变量的函数。 一阶偏导的离散化公式 一般采用欧拉公式表示 有时为了保证系统和稳定性， 对时间的差分往往采用向后公式 5.2 离散化公式 对于二阶偏导，我们可以通过对泰勒展开式处理技术得到下面离散化计算公式： 5.2 离散化公式推导 将uk+1在uk处按二阶泰勒式展开： 将uk-1在uk处按二阶泰勒式展开： 二式相加得： 5.3几种常见偏微分方程的离散化计算 1、 波动方程 其中： 为初值条件 为边值条件 当该波动方程只提供初值条件时，称此方程为波动方程的初值问题，二者均提供时称为波动方程的混合问题。 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)混合问题 5.3.1 波动方程求解 方程离散化 整理可得： 边界条件 初始条件 离散化 x xi τn 5.3.1 波动方程求解 例5.1: 用数值法求解下面偏微分方程。 此微分方程，是在不考虑流体本身热传导时的套管传热微分方程.由计算结果可知，当计算的时间序列进行到72时，传热过程已达到稳态，各点上的温度已不随时间的增加而改变。如果改变套管长度或传热系数，则达到稳态的时间亦会改变。 EXCEL 5.3.2 一维流动热传导方程 与波动方程的情形类似，用差商近似代替偏商，可以得到一维流动传热传导方