对流占优发展方程组的特征法及其变网格有限元法.pdfVIP

致 谢 转眼之间我的三年研究生涯就结束了，在这三年的学习中，我要深深 感谢我的导师袁益让教授。是导师的关怀和指导让我的学习有的放矢，掌握 到扎实的基本功，涉猎到科学研究的前沿。在我的硕士论文写作过程中始终 得到导师的悉心指导、热情关怀和鼓励。我取得的每一点进步和成绩，都凝 结着导师的宝贵心血和辛勤劳动。在平时的生活中．导师用自己的言行教导 我们如何做人、做学问。我会牢记导师的教诲“老老实实做人，踏踏实实做 学问”。 我要感谢王申林教授、羊丹平教授、程爱杰教授、鲁统超教授、孙文涛 教授，感谢他们在我的学习和生活中给予的指导和帮助。 我要感谢张志跃博士、刘小华博士、杜宁博士、杨青博士、王同科博士、 感谢他们在学习上给予我的真诚建议和有益探讨。 我要感谢感谢山东大学数学与系统科学学院对我的教育和培养。 最后．我要感谢所有帮助我的老师、同学和朋友们。 对流占优发展方程组的特征法及变罔格有限元法 ～一一～ 一—一————一——一———__———一一…，一 郑洪文 (山东大学数学与系统科学学院，济南250100) 摘要 全文共分两章。 第一章对一般的对流占优的抛物型方程组提出了特征变网格有限元格 式，并给出了最优误差估计。近年来由于实际计算豹需要而发展起来的变网 格有限元法，对不同时刻的空间区域采取不同的有限元嘲格．这种灵话的变动 网格对于许多实际问题有很大的应用价值：如油水两相驱动问题中的油水前 缘，燃烧中的火线，以及变动区域问题等． f7对于一般的抛物型方程组的变网格有限元，国内外学者作了许多工作 [4 1715]Z6]．然而在许多实际问题中，抛物型方程组都带有对流项，为对流占优 的问题．处理这类问题常采用特征法[1][2]． 考虑如下方程组： (刈)∈Q×(0用 {詈+眦”VpV，o(训)V，H)=几∽， lu(x，o)=Uo(x)， x∈Q 其中QEtt“为有界区域．aQ为Q的逐片光滑边界．工为正整数，边界条件 以D为周期． 假定A(x，H)=(口。(x，“))M一致正定有界，即存在正常数C。，C，，使对任一 向量fER‘，有 c。卿s圭％(刈赡孝，-c,It：]12 r．J·I 这里悯c2=圭㈧2 本章中采用特征法结合变网格有限元法解决了一般的带有对流项的抛物 型方程组．利用特征法，可以采用较大的时间步长而不影响逼近精度．且数值 稳定．因此本文所采用的方法的研究有比较重要的理论和实用价值． 令m=(妒l，…，妒￡)，善e中妒．=(I+f6(工，HX2){ 砒为沿警州刈№的特征施其中毒=去睁驰∽-甲) 则原方程组可以改写为 {m罢～，o(训凡“)-，(枷(x，r)∈Qx(0，T】 IⅣ(‘o)=120(z)， 工∈Q 全离散特征变网格有限元格式如下： 虾一：，v=0，Vv∈S舯。．m0 五：一“^0，v=0， Vv∈s1，m=0． 竺王：：；兰￡，V+口(d；；“：“，V)=_，I二：)，V，VV∈s。“，m。 从而可以得到主要结果，设“∈C2(O，r；Q)(即“的各个分量属于 (即“，_au的各个分量属于L2(o，r；H“，。(哪nH；(Q))．网格的变动次数为 讲 M，M．打为有界数，则对于计算格式(4．a)(4．b)(4．C)有如下能量模估计： --lt州+皿∑”一“砒=o(^2‘+(址)2) “ ’。 ，／ 其中，为一确定的正数，h=sup九，九为S。空间的