完整版大连理工大学高等数值分析抛物型方程有限差分法.doc

抛物型方程有限差分法 抛物方程差分法的构造在空间方向上与椭圆方程类似，在时间方向上用一阶差商代替代替一阶微商。然后在时间方向上逐层求解。特别当空间维数较高时，可以使用局部一维格式大大降低计算量。 1. 简单差分法 考虑一维模型热传导方程 2uu??， (1.1) )xa?f(?Tt?0? 2?t?x其中为常数。是给定的连续函数。（1.1）的定解问题分a)(xf两类： 第一，初值问题（Cauchy 问题）：求足够光滑的函数，满足方程（1.1）和初始条件： ??tx,u（1.2） , ?????x0u?x,???x?? 第二，初边值问题（也称混合问题）：求足够光滑的函数，满足方程（1.1）和初始条件： ??txu, , ???????31.x?x,u0l??xl?1及边值条件 ?? ， ????3.10ul,t?u0,t?Tt?0?2假定和在相应的区域光滑，并且于，两?????????0x00,xf,l 则上述问题有唯一的充分光滑的解。,点满足相容条件． 现在考虑边值问题（1.1），（1.3）的差分逼近 lT为时间步长，其中 ，为空间步长，是取 ??h?NM MN自然数， , ； , ?????jh??xxk?y?yM,1,?k?1,?,N0,j?0,jk??表 将矩形域分割成矩形网格。其中?? y,xT?t??0?x?l;0Gji 示网格节点；表示网格内点（位于开矩形中的网格节点）的集合； GGh表示位于闭矩形中的网格节点的集合； GGh表示-网格边界点的集合。 G?Ghhh 表示定义在网点处的待求近似解，，。??kt,xuN?j?0M?0?kkij注意到在节点处的微商和差商之间的下列关系??tx,kiku??u?? （：）)x,t?(?? kjt??t??j????kt?uu,x,txu????? kjj?1k?O???? ?t???j????ktx?uux,t,u????? 1jj?k1k?2?O???? ?2t???j????ktxx,t,?uuu????? kk?1jjhO???? h?x??j????ktx?u,ux,tu????? k1kjj?hO???? x?h??j????ktx,,xt?uuu????? kjk?1?1j2h??O?? x2h???j??????ktx,xt,?uu,xt?2u2??u??? kjk?j1k?1j2???Oh??? 22xh???j 可得到以下几种最简差分格式 （一） 向前差分格式 ?1kkkkk?2uuu??uu???? ??1??1jjjjjf?a?x?ff41. jjj12?h??， ==0 ??0kk??41.uuxu??N02jjj?a为网比，则进一步，其中。取?r1M?,?11,?,N?,k?0,1j?0, 2h有 ? =+++ ????kkkk?1?ruuurufr1?24.11j1jj?1j?j此差分格式是按层计算：首先，令，得到 0k?=+++ ??0010?fr2?1ruruuuj1j?1?jjj??和边值==0于是，利用初值，可算出第一层的0kk??uux?u?N0jjj?==0，。再由取，可利用和算出??11kkuuuu1N?j?0,1,?,1k?41.10Njj，。如此下去，即可逐层算出所有2kuu1??,Nj?0,1,jj（，）。 1?1,?,M1,?,N?1,k?0j?0,??层值直接得到，如此的格由于第层值可以通过第??1?kk??的近似值。为 式称为显格式。并视ktu,xukjj若记 ??T，， ????????????????kkkkTTu?,,u?u,u????x?,,x,?x???xx,f?f?,fxf,1N1?211?2N1?12N?可写成向量形式 则显格式??41.1k?1k??f,k?0,1,?,uAu?M?1 ?0??u? 其中 r0?01?2r????1?2rr??r???? ?A0??0????r1?2rr?????r21?0?0r??若记 2u??u a?Lu? 2t?x?kk?1kkkuu??2u?uu ??jjj?1jj?1k1?Lu?a jh2?h那末截断误差 2??1?1u??????~~??????= ）=（1.5??kk12?????O?(x?,t?)?uRhO?LuLu?x,t???? kjjkhjj2t12r2?????其中是矩形，中某一点。 %%t?t?tx?x?xt(x,)1k?j1