两类抛物型方程的全离散配置解法-计算数学专业论文.docxVIP

山东大学硕士学位论文两类抛物型方程的全离散配置方法 山东大学硕士学位论文 两类抛物型方程的全离散配置方法 王 宣 欣 (山东大学数学与系统科学学院，济南250100) 中文摘要 本文给出了两类抛物型方程一线性和拟线性抛物方程的全离散配置解法．我 们对求解区域进行剖分。采用分片双三次Hermite插值对空间进行离散，对时间采 用一般的差分，在高斯节点上建立求解格式，并对格式进行分析，均得到最优阶的 日1．模估计． 本文共分为两章 第一章处理了矩形域上的线性可变系数抛物方程的初边值问题．在这一章，我 们构造了一种全离散格式来求解该问题．通过分析我们证明了时间步长充分小的情 况下，梧式是稳定的，并且关于空间的H1．模是最优的．本章分为五节． 第一节是引言．给出r方程 甓+㈦+纠u=m，刚)7(训，北Q=Q×(啦】 “(z，Y，t)=0，(z，Y，t)∈aQ Z(0，T]； 仳(z，Y，0)=720(z，y)，(z，Y)∈Q n=(0，1)Z(0，1)，aQ表示Q的边界． 三lu=一n10，可，￡)“。。+bl(x．Y，f)虬+c(x，Y．t)u； ￡2“=-a2(x，Y，t)让w+620．Y，f)“v 且 0om；ⅡS al(z，Y，￡)，a2(z，Y，t)S amax，(z，v，t)∈Q 第二节给出了所用符号的定义，即高斯节点构成的集合．＼．以及在此集合上定 义的离散型内积·，·^和范数‰ 山东大学硕士学位论文 山东大学硕士学位论文 第三节给出分析中要用到的引理和结论，着重证明了引理3．2和引理3．4，这是 本文的两个难点． 引理3．2．对于u∈Mo，成立 C。 1^S fIv[IL：(n)≤C㈣iA C1IlF{I备-(n)s一Av，”．、≤off8备t(Q) 引理3．4．1段定 ￡u=一(nl(z，Y．f)札z)。一(a2(x．Y，t)uⅣ)目+bt(z，Y，t)u。 +如(z，Y、￡)tfg+c(x，Ⅳ，∞牡，(z：Y，￡)∈奄 那么有 Lv．札．A=Ag’(，叫)+』f’(廿，u)，口，W∈^，o，t∈(o，丁l 其中^≯0=0，l，t∈(o，Tt)：Mo x且，。—÷R的双线性形式，满足： -}’(”，w)=-g’(w，口)， 甜，∈^fo，￡∈(o：T】 Q二i。一△u．F．、≤A0’(u，F)s。，max一△，F^ F∈Mo，t∈(o：T】 l-91’(r．u)一．{P2’(口，u)J sⅣl￡l—t2I一△u，uA， F∈510，t∈(o．T】 I-{‘’(F．w)l s e6一△”．i l『』l～钉，叫∈^，o，￡∈(o，Tl 在此 忙。努州器惭川协fOma删2釉啪lie(o)+116。tic(el+IJcIIc㈤ 第四节给出全离散格式：设(1．1)一(1．5)的解“充分光滑，对于n=0，1， ：J． 令u膏∈^，为u“(·)=Ⅱ(·，t。)的分片双三次Hermite插值．定义全离散格式：寻找 {u0)：：o C A，o满足：对于扎=0，1， ，J—l’有 【(u^)?+(三r}+￡；+12)a”h；㈣=，一{(∈)，{∈A 初始条件为： (uh)o(f)=“胃(∈，O)．f∈A II 山东大学硕士学位论文对应于上述格式，我Hj定义截断误差露+÷(f)，∈∈A，乱=0，1．．一，‘，一1 山东大学硕士学位论文 对应于上述格式，我Hj定义截断误差露+÷(f)，∈∈A，乱=0，1．．一，‘，一1 《一{(∈)：【(a“)n+{一(u日)?+(L?+{+￡2n+{)(ttn+；一面芽})】(∈) 然后分析截断误差得到 定理4．1．假定啦，6。(；：1，2)，c满足： 。i，玩．c∈e(国)，i：1，2；露+{由(4．3) 式定义．如果u满足： u∈c(fo，列，H5(Q))n e加，2(国)n Co船(国)r／CO,O,3(国)， i=0．1，2，J=1，2．acu∈c([o，T】．日4(Q))，男B么 At∑I|刀+5|J专≤c[(△t)4+hs] 第五节利用上一节的定理4．1得到稳定性分析 定理5．I(H1一稳定性)．设{”：：o c Mo满足 【?+(￡1n一{+L：一j1)。n+}](∈)：叫n(∈)．∈∈A．礼=0．1，．一，J—l 那么当j￡充分小时，有下式 max e + △ r【 “∑㈣ 然后进一步得到误差估计的结论 定理5．2(H1一收敛性)．设u是问题(1．1)一(1．6)的解，且满足定理4．1的假设 设{u^)：：o c Mo是全离散格式(4 1)一(4．2)的解，那么有如下估计： 。m。。a。x。II““～u：IIH-(n)≤G{(△￡)2+^3) 第二章处理了拟线性抛物方程的初边值问题，我们采用一维问题的处理方式给 出了近似解的求解格式，并且分析得到近似解的存在唯