抛物型方程有限差分方法的应用-报告

2015年秋季学期研究生课程考核（读书报告、研究报告）考核科目:偏微分方程数值解法学生所在院（系）:理学院数学系学生所在学科:数学学生姓名:Hiter学号:1XS012000学生类别:考核结果阅卷人第1页共16页抛物型方程有限差分方法的应用摘要抛物型偏微分方程是一类比较重要的偏微分方程。热传导方程是最简单的一种抛物型方程。热传导方程研究的是热传导过程的一个简单数学模型。根据热量守恒定律和傅里叶热传导实验定律可以导出热传导方程。在本篇论文中，将先详述抛物型偏微分方程的有限差分法的相关知识，然后给出抛物型方程的两个具体的应用实例。关键字：抛物型方程，差分格式，应用AbstractParabolicpartialdifferentialequationisakindofimportantpartialdifferentialequation.Theheatconductionequationisoneofthesimplestparabolicequations.Theheatconductionequationisasimplemathematicalmodeloftheheatconductionprocess.HeatconductionequationisderivedbasedonthelawofconservationofheatandFriyege'slawofconduction.Inthisthesis,wefirstgiveadetailedknowledgeofthefinitedifferencemethodforparabolicpartialdifferentialequations,andthengivetwospecificexamplesoftheapplicationoftheparabolicequation.Keywords:parabolicequation,differencescheme,application0前言抛物型方程是偏微分方程中的三大方程（另两种为双曲型方程和椭圆型方程）之一，如何去研究抛物型方程的性质在《偏微分方程数值解法》的课程中占有很大的比例。如果要研究抛物型方程，我们一般从以下几个方面来研究，分别是：定界问题、格林函数、极值原理、解的正则性、抛物方程、拟线性蜕化和反应扩散方程。但是由于所学知识有限，所以我们在此只简单的介绍抛物型方程的有限差分法并给出两个应用实例。1抛物型方程有限差分法1.1简单差分法考虑一维模型热传导方程)(22xfxuatu，Tt0(1.1)其中a为常数。)(xf是给定的连续函数。（1.1）的定解问题分两类：第一，初值问题（Cauchy问题）：求足够光滑的函数txu,，满足方程（1.1）和初始条件：xxu0,，x（1.2）第二，初边值问题（也称混合问题）：求足够光滑的函数txu,，满足方程（1.1）和初第2页共16页始条件：xxu0,,lxl13.1及边值条件0,,0tlutu，Tt023.1假定xf和x在相应的区域光滑，并且于0,0，0,l两点满足相容条件,则上述问题有唯一的充分光滑的解。现在考虑边值问题（1.1），（1.3）的差分逼近。取Nlh为空间步长，MT为时间步长，其中N，M是自然数，jhxxj,Nj,,1,0；kyyk,Mk,,1,0将矩形域GTtlx0;0分割成矩形网格。其中jiyx,表示网格节点；hG表示网格内点（位于开矩形G中的网格节点）的集合；hG表示位于闭矩形G中的网格节点的集合；h表示hG-hG网格边界点的集合；kju表示定义在网点kitx,处的待求近似解，Nj0，Mk0。注意到在节点kitx,处的微商和差商之间的下列关系（(,)kjkjuuxttt）：Otutxutxukjkjkj,,12112,,OtutxutxukjkjkjhOxuhtxutxukjkjkj,,1hOxuhtxutxukjkjkj,,12112,,hOxuhtxutxukjkjkj第3页共16页222211,,2,hOxuhtxutxutxukjkjkjkj可得到以下几种最简差分格式（一）向前差分格式kjkjuu1jkjkjkjfhuuua2112，jjxff14.1jjjxu0，ku0=kNu=024.1其中1,,1,0Nj，1,,1,0Mk。取2har为网比，则进一步有1kju=kjru1+r21kju+kjru1+jf14.1此差分格式是按层计算：首先，令0k，得到1ju=01jru+r210ju+01jru+jf于是，利用初值jjjxu0和边值ku0=kNu=0，可算出第一层的1ju，1,,1,0Nj。再由14.1取1k，可利用1ju和ku0=kNu=0算出2ju，1,,1,0Nj。如此下去，即可逐层算出所有kju（1,,1,0Nj，1,,1,0Mk）。由于第1k层值可以通过第k层值直接得到，如此的格式称为显格式。并视kju为kjtxu,的近似值。若记TkNkkkuuu121,,,u，TNxxx121,,,，TNxfxfxf121,,,f则显格式14.1可写成向量形式011,,1,0,ufAuuMkkk其中rrrrrrrrrr21002100210021A第4页共16页若记22xuatuLukjkjkjhuuuL112112huuuakjkjkj那么截断误差为uRkjkjkjhLutxuL,1=Otxturkj)~,~(2112122=2hO（1.5）其中(,)jkxt是矩形11jjxxx，1jkttt中某一点。事实上，uRkjkjxu222+2Okjxuha442ˆ12=kjxu222+2O22222ˆ112Otuahakj=211212ah222~Otukj=21121r222~Otukj=2hO。这里222244xuxaxuatuaxa122tux22txu2322tutut222xuattxu23故22tu44244xuaxuaa，从而44xu221tua（二）向后差分格式kjkjuu1jkjkjkjfhuuua2111112jjxff16.1jjjxu0，ku0=kNu=026.1其中1,,1,0Nj，1,,1,0Mk。取2har为网比，则进一步有第5页共16页rkju1+r211kjur11kju=kju+jf16.1按层计算：首先，取0k，则利用初值jjjxu0和边值ku0=kNu=0，来确定出第一层的1ju，1,,1,0Nj，即求解方程组：r11ju+r211jur11ju=0ju+jf1,,1,0Nj，ku0=kNu=0。求出1ju，在由14.1取1k，可利用1ju，解出2ju，1,,1,0Nj。如此下去，即可逐层算出所有kju，1,,1,0Mk。如此每层必须解一个三对角线性方程组的格式称为隐格式。并视kju为kjtxu,的近似值。直观地说，采用显式格式进行求解既方便又省工作量。但是，后面我们将看到，有些情况用隐式格式更为便利。（三）Grank-Nicholson法将向前差分格式和向后差分格式做算术平均，得到的差分格式称之为六点对称格式，也称为Grank-Nicholson格式：kjkjuu1jkjkjkjkjkjkjfhuuuhuuua211111211222jjxff18.1jjjxu0，ku0=kNu=028.1进一步，2r11kju+r11kju2r11kju=2rkju1+r1kju2rkju1+jf18.1按层计算：首先，取0k，则利用初值jjjxu0和边值ku0=kNu=0，来确定出第一层的1ju，1,,1,0Nj，即求解方程组：2r11ju+r11ju2r11ju=2r01ju+r10ju2r01ju+jf1,,1,0Nj，ku0=kNu=0。求出1ju，在由18.1,取1k，可利用1ju，解出2ju，1,,1,0Nj。如此下去，即可逐层算出所有kju，1,,1,0Mk。若记22xuatuLu第6页共16页kjkjkjhuuuL13jkjkjkjkjkjkjfhuuuhuuua211111211222那么截断误差为uRkjkjkjhLutxuL,3=22hO（1.9）注意21kjtu=kjkjtxutxu,,12O又122kjxu2122kjxu21232kjtxu2122422kjtxu3Okjxu222122kjxu21232kjtxu2122422kjtxu3O两式相加21kjtukjkjxuxu22122212O211211111,,2,,,2,21htxutxutxuhtxutxutxukjkjkjkjkjkj22hO而21kjtu2122kjxua+jf故有kjkjtxutxu,,1211211111,,2,,,2,2htxutxutxuhtxutxutxuakjkjkjkjkjkj22hO。（四）Richardson格式第7页共16页211kjkjuu2112huuuakjkjkj+jf（1.10）进一步1kju=r2（kju1kju2+kju1）+1kju+2jf110.1这是三层显式差分格式。显然截断误差的阶为22hO。为使计算能够逐层进行，除初值0ju外，还要用到1ju。它可以用其他双层格式提供。Richardson格式的矩阵形式为：