Sobolev方程的半有限元方法.pdfVIP

第49卷第3期 华中师范大学学报(自然科学版) V01．49No．3 CHINAN()RMALUNlVERSITY(Nat．Sci．)Jun．2015 OFCENTRAL 2015年6月 JOURNAL 05 文章编号：1000—1190(2015)03—0334 Sobolev方程的半有限元方法 辜继明1，俞 辉H，瞿少成2 (1．三峡大学理学院，湖北宜昌443000；2．华中师范大学信息技术系，武汉430079) 摘 要：对Sobolev方程采用半有限元法进行数值模拟．通过将空间变量和时间变量分离，得到 Sobolev方程的离散格式．首先对空间变量应用有限元方法进行离散化，得到常微分方程组的初值 问题；再对时间变量应用有限差分法进行离散化，得到一系列线性方程组，求解可得到Sobolev方 程的数值解．本文从理论上推导出了本文所讨论的Sobolev方程半有限元算法的矩阵算法格式，分 析了其可行性．在最后给出了数值例子，从数值例子中进一步验证了半有限元方法的可行性． 关键词：半有限元法；Sobolev；方程有限元法；有限差分法 中图分类号：()175 文献标识码：A 讨论以下Sobolev方程 题中的时间变量应用有限差分法再进行离散化，得 象一杀(Ⅱ3u 3u)+弛∽， 到一系列线性代数方程组，这一过程称为全离散化 过程． (工，f)∈力×(o，丁]； 很多的物理或工程问题，都可以用微分方程来 M(z，f)一0， (z，￡)∈an×Eo，￡]； 表示，求解微分方程的数值解最早是用有限差分 ，，，●●●●●●●●●●●●●●●，，．、●●．●●●，，●●●●●●●L “(z，0)一“o(上)，z∈力． 法，有限差分法在求解微分方程时效率很高，尤其 c 其中，n R“(d一1，2，3)为方程的有界区域，an 是在求解一、二维问题时，但是在高维复杂的问题 是方程的边界，Uo(T)，f(x，￡)是已知函数，盘是关 时，因为跳跃误差存在，导致迭代出现误差累积，从 于z，t的函数口(T，f)的连续函数，且＆(z，￡)0， 而会有不收敛现象．有限元法是求解微分方程数值 6(z，t)≥0． 解最重要的方法之一，在高维复杂的问题求解时， Sobolev方程是偏微分方程，主要是研究不定 可以得到较高的精度，但是所需计算的时间往往 常微分方程．对Sobolev方程的研究具有重大的物 较大． 理工程意义．通过对Sobolev方程的研究，可以分 本文结合有限差分法和有限元法，发挥各自的 析流体穿过裂缝岩石的理论，可得出土壤湿气迁移 优势，在空间变量上，应用有限元求解，得到精度较 的过程，同样也可解决不同介质的热传导问题．人 高的关于时间的微分方程，再用有限差分法，可以 们对Sobolev方程很早就开始研究了，如标准混合 高效的求解动态的偏微分方程．这种半有限元法比 有限元法应用于Sobolev方程口j，Galerkin一有限元