几类发展方程的计算格式-计算数学专业论文.docxVIP

山东大学硕士学位论文几类发展方程的计算格式 山东大学硕士学位论文 几类发展方程的计算格式 李 广 瑞 (山东大学数学与系统科学学院，济南250100) ● 中文摘要 全文共分三章． 第一章对热篮昱直程提出了有限差分非匹配区域分解算法，并给出了最大模 ’ 误差估计．许多数学物理及工程问题，如油、气藏的勘探开发，天气预报，反应堆 计算等，都归结为对大型偏微分方程的求解．区域分解算法对解决这些问题有很 大鳆帮助作用． f前人已经在区域分解方面作了大量的工作[1，2，3，4】．本文的不同之处在于使用 了非交叠域非匹配的网格，这样就可以对不同的区域，使用相应的或粗或细的网 格，以达到节省计算量又能满足精度的双重目的． 考虑如下方程： 一 = (z，Y)∈Q，t∈(0，列 |l Z y (o，Y)∈Q ，Jll_J(1l-【 丝巩Ⅱ“ 曩％ 舭删㈣ lI 0护O (正，Y)∈aQ 【2为二维区域，aQ为Q的边界．为了讨论简单，Q取为单位正方形区域，而且 全文把主要的精力放在对二分子域情况的处理，多分的情况类似．设左右两个子 - 域在Y方向上的步长取为^牡和hya，并且存在mL，mR满足mLh姐=mnh北， 构造差分格式： ％=札：，在边界点 (3．1) 镰】+1一略1[yt·n。-1_“-“n-镨1】 · 乩％：反埘吧一。。{堕塑些型笔-竺竖卜(1-e1)四凡。％ F 一万it一一。一一一一。“Z1，f(【等】+1)址一 jh叫U1．日-l【等1+(地露一【器】‰)吃tc等1+· 山东大学硕士学位论文===============================：： 山东大学硕士学位论文 ===============================：： ：= ：：：： ：： 型坠： ：=：： K一 等一卜一 兰兰篙竺L ! 一c7孑1+职n+-H1，J)=0 A类内界点 (3．2) 。 L。。％：a出％一。。{堕塑些竺等}兰坐)一(1一。。)鳄山。吧“¨．-燮I】+l_Uwg?t q7。《奎】 2 1，[(【等】+1)hua- jh一1UM“-t“舒】+(地L一[3mm∥R w，。Mn--1．【晋】+1 。 日2。huR 【(【等】+1)h，R—jh 1U讲-晋1)+(‰一曙‰)叼丧m 一呀1+D．，n--H1，，)=0 B类内界点 (3．3) Lsn％=Ot，△t％一磋，n％一骘,hyL uZ“,3在左内点 (3 4) Ltn％=Ot，△t％一磋，n矿z“,l一《，h。％在右内点 (3．5) 用极值原理证明了最大模误差估计： J。乙I≤Q[；(△抖忽。+^乳)+；(At+h：彳知)，；日(△¨晋h2+H2+hyR+huL) =K(At+h2+日3+Hh讧+hy北)少 。 第二章讨论的方程是双曲型热传导方程，即在热传导方程中加入了一个带有 小系数的时间二阶导数． 怀为普遍的自然现象，热传导满足热力学第一定律，即能量守恒定律．然而 利麻同的热流密度本构关系(热流密度矢量与温度梯度之间的关系)，将得到不 同类型的热传导方程．本章所用到的方程是由Cattaneo-Vernotte本构关系得到 一 的，它不仅含有经典热传导方程的某些特性，还含有波动方程的某些特性．双曲型 热传导方程的表达式为： t‘￡+丁bt‘tt=o“$￡(z，t)∈(0，f)×(o，T】 II 山东大学硕士学位论文陛蕊。 山东大学硕士学位论文 陛蕊。 - 的混合问题，给出了全隐的和加权平均的差分格式，利用极值原理证明得到最 ／ 大模误差估计： ，．ma，x， 1 e州i l=0(r+h2) 0{J’1nN一1。 - 第三章对波动方程提出了变网格质量集中有限元．由于波的峰值随时间变化， 在峰值附近网格剖分要局部加才能保证精度而不影响计算量．而质量集中有限元 的数值稳定性．括虑如下问题； 、 f铬=△u+／(x，t)(茁，t)∈Q×R {然篙曷l 象(z，o)=咖(z) ：：三 ㈦，， 。 7 z∈n f o_巩q=△t‘+y(x，t)(z，t)∈Q×R l碧o，o)=Q o，t)∈Q×R {u(z，0)=U0(z) ￡∈n， (1．2) · I u(州)：o z∈aQ×【o，卅． (攀，u)^+口(咚l，)=(“}，u)V∈sn+l (2·4) ◆ (粤≯，”)=(啦芦，”) 怕∈岛+- (2．5) (VUh—v％，)=0 Yv∈So (2．6) (磷一毋，t，)=0 Vv∈岛 (2·7) (v嘏一V硼h，u)=0 V∈+ln≥0(2-8) (识一硪，13)=0 铷∈R+l，n≥0(2．9) III 山东大学硕士学位论文经过理论分析，在u(z，t)和，缸)具有讨论要求的光滑性的条件下，可得到如下的 山东大学硕士学位论文 经过理论分