偏微分方程数值解法教学提纲.docxVIP

〈〈偏微分方程数值解法? 课程设计 姓 名： 姓晓霜 学 院： 理学院 专 业： 信息与计算科学 班 级： 0911012 学 号： 091101218 指导老师：翟方曼 2012年12月14日 一、题目 用六点对称差分格式计算如下热传导方程的初边值问题 c 2 即 C 8 U c 』c L 』 —=2— ,0 x1,0t 菱1 ct ex x ，u(x,0) =e ,0 Mx 1 u(0,t) =e2t,u(1,t) =4^,0 宜由 i 已知其精确解为 u(x,t) = ex * 二、理论 考虑的问题 考虑一维模型热传导方程 TOC \o 1-5 \h \z - -2 u - u (1.1) 一 = a— + f (x) , 0 t 玄T ；=t :x 其中a为常数。f (x)是给定的连续函数。(1.1)的定解问题分两类： 第一，初值问题(Cauchy问题)：求足够光滑的函数u(x, t),满足方程 (1.1 )和初始条件： (1.2 ) u(x, 0)=甲(x), —°ox8 第二，初边值问题(也称混合问题)：求足够光滑的函数u(x, t),满足 方程(1.1)和初始条件： 1.3 1 u x, 0[=峪[x , -1 ：x ：：l 及边值条件 (1.3》 u(0, t )=u(1,t )=0 , 0 玄 t 玄 T 假定f (x )和中(x件相应的区域光滑，并且丁(0, 0 ), (1,0 )两点满足相容条件， 则上述问题有唯一的充分光滑的解。 现在考虑边值问题(1.1 ), (1.3)的差分逼近 取 h=L为空间步长， J 工为时间步长，其中N , M是自然数， N M x=xj=jh, (j=0,1，…，N )； y = yk=kj (k=0,1，…，M) 将矩形域G = fo^l; ^^T）分割成矩形网格。其中（x“yj）表小网格节点; Gh表示网格内点（位丁开矩形 G中的网格节点）的集合； Gh表示位丁闭矩形G中的网格节点的集合； L表示Gh- Gh网格边界点的集合。 k _三it（Xj，）：u：表示定义在网点（x,tk） k _ 三it（Xj，）： 注意到在节点（Xi,tk）处的微商和差商之间的下列关系（ ：t u Xj ,tk 1 -u Xj, tk U(Xj, tk* Au(Xj, tkn ) (cu V 2 , =—— O j U Xj i U Xj i,tk -U Xj, tk _ 包 i十°(h) EX /j u(Xj ,tk )—U(Xj_L, tk ) (cu \ Li 、 =——I + O(h )h ex ,j u(Xj+,tk J-uj」, tk ) (du 华工-j = ——i十Olh 2h l 欲 J +O(h2)u Xj 1, tk -2u Xj, tk +O(h2) h2 区域网格剖分 取空间步长h = M 和时间步长丁=%，其中N,M都是正整数。用两族平 行直线 x = jX = (j h OjL 1 ,和 Nt =tk =kw(k =0,1,L ,M )将 矩形域 ￡ =｛。4x l;0 t《t｝分割成矩形网格，网格节点为（Xj,tk）。以Gh表示网格内 点集合，即位丁开矩形的网点集合； Gh表示所有位丁闭矩形 G的网点集合； 「h =Gh -Gh是网格界点集合。 其次，用u：表示定义在网点（Xj,tk）的函数，0勺苴N,0苴k壬M 建立相应差分格式 数值分析中，Crank-Nicolson方法是有限差分方法中的一种，用丁数值求 解热方程以及形式类似的偏微分方程。它在时间方向上是隐式的二阶方法，数值 稳定。该方法诞生丁 20世纪,John Crank 与 稳定。该方法诞生丁 20世纪, ①向前差分格式 k 1 kU k 1 k Uj -Uj T u j 1 - 2u j u j i -J^-2j~~fj fj= f Xj h 1.4 2U0 1.4 2 U0=?j=%Xj ), k k U0 =Un =0 ②向后差分格式 k 1 k k 1 k 理1 k 1 1.5 1u j — u j u j 1 _ 2Uj u j 1.5 1 =a 2 fj fj = f Xj h 1.5 2U0=%* 1.5 2 U0=%* )， k k u 0 = Un =0 将向前差分格式和向后差分格式做算术平均， 得到的差分格式称之为六点对 称格式,也称为 Grank-Nicholson 格式:1.6 1k 1 kUj —Ujk - ka Uj + -2Uj +u2 称格式, 也称为 Grank-Nicholson 格式: 1.6 1 k 1 k Uj —Uj k - k a Uj + -2Uj +u 2 h2 fj f j = f Xj 1.6 2 0 Uj k