对流占优Sobolev方程投影稳定化有限元方法.pdfVIP

2014年 2月 计 算 数 学 第 36卷第 1期 Feb．，2014 MATHEMATICA NUM ERICA SINICA Vo1．36．No．1 对流占优的 Sobolev方程的投影稳定化 有限元方法 ) 周 琴 潘雪琴 冯民富。 (四川大学数学学院，成都610064) 摘 要 对于对流占优的Sobolev方程，提出了一种新的投影稳定化有限元方法，建立了半离散和全 离散的投影稳定化格式，给出了解的稳定性和收敛性分析．该方法能够有效克服对流占优，与内罚 方法相比，投影格式更简单，计算量更小，且得到的C—N格式是无条件稳定的，时间精度达到了 二阶．最后，通过实验证明，数值结果与理论结果完全一致． 关键词：Sobolev方程；对流 占优；投影方法；离散格式 MR (2000)主题分类：65N30，65M30，76M10 1．引 言 本文考虑如下带有狄利克雷边界条件和初值边界条件的Sobolev方程 f札t+ ()V·u— ·(a(x)Vu+b(x)Vut)+ ()=f(x，t)，(，t)∈Q×J， { (，t)=0， (，t)∈aQ×J， (1．1) 【u(z，0)= 0()， z∈Q， 其中 c (d=2，3)是有界凸多边形或多面体区域，aQ为其边界，时间区域 J= (0，]， 0T∞，Ut：=Ou／Ot． ()=(1()，2())为给定的流速，满足 Lipschitz连续条件且在 区域Q上 V ·()=0．f(x，t)为源函数，uo(x)为初始数据，0()j6()，()为已知的非负 函数．假定存在正常数a{，6(=0，1)，对任意的 ∈ 有 0ao a(x) al，0bo b(x) bl， 0()，b(x)《I()I=(()+ ；())／。 Sobolev方程在岩石缝隙的流体力学 ．1J’热力学_2_等许多数学物理问题方面应用十分 广泛．但在实际运用中，Sobolev方程的计算域往往是不规则的，或问题本身比较复杂，要求 出其解析解是比较困难的，有效的方法是求其数值解．目前已有一些采用有限元方法来研究 Sobolev方程的报道．这些报道主要倾向于研究不含对流项的Sobolev方程 [a-sJ．但随着研究 的不断深入，含对流项的Sobolev方程也越来越受到关注． 采用有限元方法研究含对流项的Sobolev方程的困难主要在于：对流项 ()·Vu(x)的 存在，常常会引起数值解在内层或者在边界层的震荡．为克服这一困难，近年来，人们提出了 2013年6月 9日收到． )基金项 目：国家 自然基金 资助项 目 )通讯作者：冯民富，E-Mail：fmf@wtjs．an． 100 计 算 数 学 许多稳定化有限元方法．GU提出了特征有限元方法 【，Sun与Yang提出了两种有限差分流 线扩散格式_101¨】，后来又提出了间断的Galerkin内罚方法 L12_，Gao和 Rui考虑了基于最小 二乘法的特征混合有限元方法 l13】，Sun和Ma给出了基于间断 Galerkin方法的空问与时间 离散格式 |l4]，Zhang和Wang等提出了连续的内罚有限元方法 ．本文的主要 目的在于提 出一种新的投影稳定化有限元方法来研究含对流项的Sobolev方程，据我们所知，目前还没有 关于这一方面的研究报道． 基于投影的稳定化有限元方法最早源于 Guermond的思想，他在 [16】中提出了Subgrid modeling方法，通过 bubble函数来增强有限元空间，使得人工粘性项仅作用在 bubble函数 空间．在此工作基础上，Layto