计算方法偏微分方程数值解精选课件.ppt
《计算方法偏微分方程数值解精选课件.ppt》由会员分享,可在线阅读,更多相关《计算方法偏微分方程数值解精选课件.ppt(32页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第一页,本课件共有32页上一页下一页回目录休息Email:JPhone:02985583997本章要求本章要求教学目的教学目的 讲解:偏微分方程离散格式及求解的一般过程教学要求教学要求E熟记一阶及二阶偏微分方程的离散格式;C精通用EXCEL迭代对偏微分方程求解;探索用两数组交替更新的办法进行编程求解;F延伸对化学反应工程中物理场的模拟进行尝试。教学重点教学重点各种偏微分方程的离散与求解EXCEL循环迭代问题教学难点教学难点特殊边界条件的引入与应用第二页,本课件共有32页上一页下一页回目录休息Email:JPhone:029855839975.1 偏微分方程简介&偏微分方程如果一个微分方程中出现
2、多元函数的偏导数,或者说如果未知函数和几个变量有关,而且方程中出现未知函数对几个变量的导数,那么这种微分方程就是偏微分方程。在化工或化学动态模拟方程中,常常有一个自变量是时间,其它的自变量为空间位置。如果只考虑一维空间,则只有两个自变量;如果考虑两维空间,则有3个自变量。许多化工过程均是通过对偏微分方程的求解进行工艺参数的确定或数值模拟。第三页,本课件共有32页上一页下一页回目录休息Email:JPhone:029855839975.1 偏微分方程简介&偏微分方程的分类线性微分方程Linear partial differencial equation拟线性微分方程Quasilinear pa
3、rtial differencial equation非线性微分方程Nonlinear partial differencial equation第四页,本课件共有32页上一页下一页回目录休息Email:JPhone:029855839975.1 偏微分方程简介&数学上的分类:椭圆方程 Elliptic抛物线方程 Parabolic双曲线方程 Hyperbolic&物理实际问题的归类:波动方程(双曲型)一维弦振动模型:热传导方程(抛物线型)一维线性热传导方程拉普拉斯方程(椭圆型)稳态静电场或稳态温度分布场)第五页,本课件共有32页上一页下一页回目录休息Email:JPhone:02985583
4、9975.1 微分方程的求解思路&求微分方程数值解的一般步骤:Step1区域剖分:首先按一定规则将整个定义域分成若干小块Step2微分方程离散:构造离散点或片的函数值递推公式或方程Step3初始、边界条件离散:根据递推公式,将初值或边界值离散化,补充方程,启动递推运算Step4 数值解计算:求解离散系统问题&微分方程的定解问题 离散系统的求解问题第六页,本课件共有32页上一页下一页回目录休息Email:JPhone:029855839975.2 离散化公式&将自变量在时间和空间上以一定的间隔进行离散化,则应变量就变成了这些离散变量的函数。&一阶偏导的离散化公式一般采用欧拉公式表示&有时为了保证
5、系统和稳定性,对时间的差分往往采用向后公式第七页,本课件共有32页上一页下一页回目录休息Email:JPhone:029855839975.2 离散化公式对于二阶偏导,我们可以通过对泰勒展开式处理技术得到下面离散化计算公式:第八页,本课件共有32页上一页下一页回目录休息Email:JPhone:029855839975.2 离散化公式推导&将uk+1在uk处按二阶泰勒式展开:&将uk-1在uk处按二阶泰勒式展开:&二式相加得:第九页,本课件共有32页上一页下一页回目录休息Email:JPhone:029855839975.3几种常见偏微分方程的离散化计算&1 1、波动方程波动方程其中:为初值条
6、件 为边值条件 当该波动方程只提供初值条件时,称此方程为波动方程的初值问题,二者均提供时称为波动方程的混合问题。第十页,本课件共有32页上一页下一页回目录休息Email:JPhone:029855839975.3.1 波动方程求解对于初值问题,是已知t=0时,u与 依赖于x的函数形式,求解不同位置,不同时刻的u值。而 u是定义在 的二元函数,即上半平面的函数。对于混合问题除初值外,还有边值。是已知初值及x=0及x=l 时u依赖于t的函数,求解不同位置x,不同时刻的u值。此时u是定义在 的带形区域上的二元函数。xt 0a)初值问题初值问题tx0lb)混合问题混合问题第十一页,本课件共有32页上一
7、页下一页回目录休息Email:JPhone:029855839975.3.1 波动方程求解方程离散化整理可得:边界条件初始条件离散化xxin第十二页,本课件共有32页上一页下一页回目录休息Email:JPhone:029855839975.3.1 波动方程求解&例例5.1:5.1:用数值法求解下面偏微分方程。此微分方程,是在不考虑流体本身热传导时的套管传热微分方程.由计算结果可知,当计算的时间序列进行到7272时,传热过程已达到稳态,各点上的温度已不随时间的增加而改变。如果改变套管长度或传热系数,则达到稳态的时间亦会改变。EXCELEXCEL第十三页,本课件共有32页上一页下一页回目录休息Em
8、ail:JPhone:029855839975.3.2 一维流动热传导方程与波动方程的情形类似,用差商近似代替偏商,可以得到一维流动传热传导方程的混合问题的差分方程,以其解作为流动传热传导方程的近似解。2 2、一维流动热传导方程的混合问题、一维流动热传导方程的混合问题离离散散化化第十四页,本课件共有32页上一页下一页回目录休息Email:JPhone:02985583997将上式进行处理得到:该式是显式格式。只要保证式中各项系数大于零,一般情况下是稳定的,可以获得稳定的解。分析上式可以发现,当为了提高数值精度取适当小的x时,最有可能小于零的系数是uin的系数,若要保证此项系数大于零,此时t必须
9、相应地更小,会导致计算量将大大增加,这是显式格式的缺点,为了克服此缺点,下面提出一种隐式格式:偏微分方程在点上进行离散化,且对时间的偏微分采用向后欧拉公式得到原偏微分方程的离散化公式:5.3.2 一维流动热传导方程第十五页,本课件共有32页上一页下一页回目录休息Email:JPhone:02985583997从图5-3中可见要由初值及边界条件一排一排推上去是不行的,需解线性方程组,同时添上二边界条件:正好共有m+2个方程,同时有m+2个变量,就能解出n+1排上各点值。这样,每解一个线性方程组,就可以往上推算一排点的u值,虽然引入了方程组的求解,有可能增加计算量,但由于隐式格式无条件稳定,t的取
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算方法 微分方程 数值 精选 课件
限制150内