二阶抛物问题的最小二乘扩展混合元格式.pdfVIP

第8卷 第23期 2008年 l2月 科 学 技 术 与 工 程 Vo1．8 No．23 Dec．2008 1671—1819(2008)23-6318—03 ScienceTechnologyandEngineering ⑥ 2008 Sei．Tech．Engng． 二阶抛物问题的最小二乘扩展混合元格式 丁 胜 (山东师范大学数学科学学院，济南250014) 摘 要 对二阶抛物问题提出了最小二乘扩展混合元数值模拟格式。数值算例验证了格式的有效性。 关键词 二阶抛物问题 最小二乘 扩展混合有限元方法 中图法分类号 0241．82； 文献标志码 A 考虑下列二阶抛物问题 数为 lI }l0,D=(，)T，lI II，n=(I1 jf：．n+ u(，t)一V-(口()Vu(x，t))=／ ，t)， ． 1 ． IlV n)i，I1gI1 n)=(IIg o+}1 ‘q (，t)∈力×J， l1．。)寺，其中(·，·)表示空间( (力)) 中的内 “(，f)=0， (，￡)∈F×J， 积。作时间剖分 7l：T／N，N为正的常数，在时间t “(，0)= 0()， ∈ ， =nJr，n=1，2，…，Ⅳ处逼近方程的解。n：uK是 ． ， (1) ^ ∈r^ 的一个剖分，h=max{diam(K)：K∈ }。设 ， (1)式中J=[0，T]，c 为有界凸域，边界厂 和 分别是 ，A和 的有限维子空间，具有通 是Lipschitz连续的。假定 口()有界，即存在常数 常逼近性质 。 ， 和 的一个标准的选法是 ，r 和O／‘，使得 和s次的分片多项式空间。 0≤ +≤口()≤ ，V ∈ (2) 记A=一7“，= ，并用一阶差分近似 ，则 对上述问题人们提出了许多数值模拟格式，如 式(1)可变成离散形式， 有限元格式_1]，混合有限元格式 J，扩展混合有限 元格式 3J，最小二乘混合有限元格式_4等。现将最 f0t“= 7 · =厂 +R ， ∈ ， 小二乘思想和扩展混合有限元方法相结合，提出并 fA +VH =0， ∈力， 分析了最小二乘扩展混合有限元格式。该格式继 {口A一 =0， ∈ ， (3) 承了最小二乘和扩展混合元方法的优点，即格式能 lM()=0， ∈ ， 同时高精度逼近未知函数，未知函数的梯度 以及流