偏微分方程的数值解方法及源程序

1、 -240- 第二十章第二十章 偏微分方程的数值解偏微分方程的数值解 自然科学与工程技术中种种运动发展过程与平衡现象各自遵守一定的规律。 这些规 律的定量表述一般地呈现为关于含有未知函数及其导数的方程。 我们将只含有未知多元 函数及其偏导数的方程，称之为偏微分方程。 方程中出现的未知函数偏导数的最高阶数称为偏微分方程的阶。 如果方程中对于未 知函数和它的所有偏导数都是线性的， 这样的方程称为线性偏微分方程， 否则称它为非 线性偏微分方程。 初始条件和边界条件称为定解条件，未附加定解条件的偏微分方程称为泛定方程。 对于一个具体的问题， 定解条件与泛定方程总是同时提出。 定解条件与泛定方程作为一 个整体，称为定解问题。 1 偏微分方程的定解问题 各种物理性质的定常（即不随时间变化）过程，都可用椭圆型方程来描述。其最典 型、最简单的形式是泊松(Poisson)方程 ),( 2 2 2 2 yxf y u x u u= + = （1） 特别地，当0),(yxf时，即为拉普拉斯(Laplace)方程，又称为调和方程 0 2 2 2 2 = + = y u x u u （2） 带有稳定热源或内部无

2、热源的稳定温度场的温度分布， 不可压缩流体的稳定无旋流动及 静电场的电势等均满足这类方程。 Poisson 方程的第一边值问题为 = = + ),(| ),( ),(),( ),( 2 2 2 2 yxyxu yxyxf y u x u yx （3） 其中为以为边界的有界区域，为分段光滑曲线，U称为定解区域， ),(),(yxyxf分别为,上的已知连续函数。 第二类和第三类边界条件可统一表示成 ),( ),( yxu n u yx = + （4） 其中n为边界的外法线方向。当0=时为第二类边界条件，0时为第三类边界 条件。 在研究热传导过程，气体扩散现象及电磁场的传播等随时间变化的非定常物理问 题时，常常会遇到抛物型方程。其最简单的形式为一维热传导方程 )0(0 2 2 = a x u a t u （5） 方程（5）可以有两种不同类型的定解问题： 初值问题（也称为 Cauchy 问题） -241- += xxxu xt x u a t u )()0 ,( , 00 2 2 （6） 初边值问题 = = = = xx t u xxxu xt x u a t u t )( )()0 ,( ,

3、 0 0 2 2 2 2 2 （11） 边界条件一般也有三类，最简单的初边值问题为 -242- = = = = = Tttgtlutgtu lxx t u xxu lxt x u a t u t 0)(),(),(), 0( 0)(),()0 ,( 0, 0 21 0 2 2 2 2 2 如果偏微分方程定解问题的解存在，唯一且连续依赖于定解数据（即出现在方程 和定解条件中的已知函数） ，则此定解问题是适定的。可以证明，上面所举各种定解问 题都是适定的。 2 偏微分方程的差分解法 差分方法又称为有限差分方法或网格法， 是求偏微分方程定解问题的数值解中应用 最广泛的方法之一。它的基本思想是：先对求解区域作网格剖分，将自变量的连续变化 区域用有限离散点（网格点）集代替；将问题中出现的连续变量的函数用定义在网格点 上离散变量的函数代替； 通过用网格点上函数的差商代替导数， 将含连续变量的偏微分 方程定解问题化成只含有限个未知数的代数方程组（称为差分格式） 。如果差分格式有 解， 且当网格无限变小时其解收敛于原微分方程定解问题的解， 则差分格式的解就作为 原问题的近似解（数值解） 。因此，用差分

4、方法求偏微分方程定解问题一般需要解决以 下问题： （i）选取网格； （ii）对微分方程及定解条件选择差分近似，列出差分格式； （iii）求解差分格式； （iv）讨论差分格式解对于微分方程解的收敛性及误差估计。 下面我们只对偏微分方程的差分解法作一简要的介绍。 2.1 椭圆型方程第一边值问题的差分解法 以 Poisson 方程（1）为基本模型讨论第一边值问题的差分方法。 考虑 Poisson 方程的第一边值问题（3） = = + ),(| ),( ),(),( ),( 2 2 2 2 yxyxu yxyxf y u x u yx 取, h分别为x方向和y方向的步长，以两族平行线jyykhxx jk =, ), 2, 1, 0,(L=jk将 定 解 区 域 剖 分 成 矩 形 网 格 。 节 点 的 全 体 记 为 ,| ),(为整数jijykhxyxR jkjk =。定解区域内部的节点称为内点，记内点 集IR为 h 。边界与网格线的交点称为边界点，边界点全体记为 h 。与节点 ),( jk yx沿x方向或y方向只差一个步长的点),( 1jk yx 和),( 1jk yx称为节点 ),(

5、jk yx的相邻节点。如果一个内点的四个相邻节点均属于U，称为正则内点，正 则内点的全体记为 )1( ，至少有一个相邻节点不属于U的内点称为非正则内点， 非正则内点的全体记为 )2( 。我们的问题是要求出问题（3）在全体内点上的数值解。 为简便记，记),(),(),(),(),( ,jkjkjkjk yxffyxujkuyxjk=。对正则内点 -243- )1( ),(jk，由二阶中心差商公式 )( ), 1(),(2), 1( 2 2 ),( 2 2 hO h jkujkujku x u jk + + = )( ) 1,(),(2) 1,( 2 2 ),( 2 2 O jkujkujku y u jk + + = Poisson 方程（1）在点),(jk处可表示为 )( ) 1,(),(2) 1,(), 1(),(2), 1( 22 , 22 += + + + hOf jkujkujku h jkujkujku jk （12） 在式（12）中略去)( 22 +hO，即得与方程（1）相近似的差分方程 jk jkjkjkjkjkjk f uuu h uuu , 2 1,1, 2 , 1,

6、 1 22 = + + + + (13) 式（13）中方程的个数等于正则内点的个数，而未知数 jk u , 则除了包含正则内点处 解u的近似值，还包含一些非正则内点处u的近似值，因而方程个数少于未知数个数。 在非正则内点处 Poisson 方程的差分近似不能按式（13）给出，需要利用边界条件得到。 边界条件的处理可以有各种方案，下面介绍较简单的两种。 (i) 直接转移 (ii) 线性插值 由式（13）所给出的差分格式称为五点菱形格式，实际计算时经常取=h，此时 五点菱形格式可化为 jkjkjkjkjkjk fuuuuu h ,1,1, 1, 1 2 )4( 1 =+ + (14) 简记为 jkjk fu h , 2 1 = （15） 其中 jkjkjkjkjkjk uuuuuu ,1,1, 1, 1, 4+= + 。 求解差分方程组最常用的方法是同步迭代法， 同步迭代法是最简单的迭代方式。 除 边界节点外，区域内节点的初始值是任意取定的。 例 1 用五点菱形格式求解 Laplace 方程第一边值问题 =+= = + )1lg(| ),( ),( 0 22 ),( 2 2 2 2 yxy

7、xu yx y u x u yx 其中1,0| ),(=yxyx。取 3 1 =h。 当=h时，利用点) 1, 1(),1, 1(),(+jkjkjk构造的差分格式 jkjkjkjkjkjk fuuuuu h ,1, 11, 11, 11, 1 2 )4( 2 1 =+ + （16） -244- 称为五点矩形格式，简记为 2 2 1 h jkjk fu , = （17） 其中 jkjkjkjkjkjk uuuuuu ,1, 11, 11, 11, 1, 4+= + 。 2.2 抛物型方程的差分解法 以一维热传导方程（5） )0(0 2 2 = a t u a t u 为基本模型讨论适用于抛物型方程定解问题的几种差分格式。 首先对xt平面进行网格剖分。分别取, h为x方向与t方向的步长，用两族平行直 线), 2, 1, 0(L=kkhxx k ，), 2 , 1 , 0(L=jjtt j ，将xt平面剖分成矩形网 格，节点为), 2 , 1 , 0, 2, 1, 0)(,(LL=jktx jk 。为简便起见，记),(),( jk yxjk=， ),(),( jk yxujku=，)( k

8、k x=，)( 11jj tgg=，)( 22jj tgg=，)( 11jj t=， )( 22jj t=。 2.2.1 微分方程的差分近似 在网格内点),(jk处，对 t u 分别采用向前、向后及中心差商公式，对 2 2 x u 采用二 阶中心差商公式，一维热传导方程（5）可分别表示为 )( ), 1(),(2), 1(),() 1,( 2 2 hO h jkujkujku a jkujku += + + )( ), 1(),(2), 1() 1,(),( 2 2 hO h jkujkujku a jkujku += + )( ), 1(),(2), 1( 2 ) 1,() 1,( 2 2 hO h jkujkujku a jkujku += + + 由此得到一维热传导方程的不同的差分近似 0 2 2 , 1, 1,1, = + + h uuu a uu jkjkjkjkjk （18） 0 2 2 , 1, 11, = + + h uuu a uu jkjkjkjkjk （19） 0 2 2 2 , 1, 11,1, = + + h uuu a uu jkjkjkjkjk （20） 2.2.2 初、边值条件的处理 为用差分方程求解定解问题（6） ， （7）等，还需对定解条件进行离散化。 对初始条件及第一类边界条件，可直接得到 kkk xuu=)0 ,( 0 , ), 1 , 0, 1, 0(nkkLL=或 （21） jjjn ijjj gtluu gtuu 2, , 0 ),( ), 0( = = ) 1, 1 , 0(=mjL （22） -245- 其中 T m h l n=,。 对第二、三类边界条件则需用差商近似。下面介绍两种较简单的处理方法。 （i）在左边界)0( =x处用向前差商近似偏导数 x u ，在右边界)(lx =处用向后差 商近似偏导数 x u ，即 )( ), 1(),( )( ), 0(), 1 ( ),( ), 0( hO h jnujnu x u hO h juju x u jn j + =

《偏微分方程的数值解方法及源程序》由会员n****分享，可在线阅读，更多相关《偏微分方程的数值解方法及源程序》请在金锄头文库上搜索。