RLW方程的两种新的有限元方法.docxVIP

第一章绪论 复杂的问题，且便于程序标准化，从而适于工程应用．由于有限元方法有上述优越性，它 自上世纪60年代以来巳作为一种独立的数值计算方法获得了迅速发展和广泛应用(见文 献13-5D． 在众多的偏微分方程中，弛w方程是包含时间变数的许多重要的数学物理偏微分方 程之一，它是一个双瞌一抛物型方程．在孤立子理论的研究及非线性波的物理问题中， 它占据了相当重要的地位，在非常广泛的领域里得到了应用．因此，对该方程的研究已受 到人们的很大重视，有许多专家、学者在这方面已做了不少工作．例如文献【6—18|J 过去，求解RLW方程的数值方法主要有有限差分方法、标准有限元法等，虽然我们 能够用这些方法求得它的近似解，但由于RLW方程的双曲性质，我们用这些方法所得到 的近似解的精度并不高，有时稳定性能也比较差，基于这个原因，近30年来，在上述基本 方法的基础上，探寻RLW方程稳定的高效、高精度的数值方法的努力始终没有间断过， 不断地涌现出新的数值方法．例如，著名学者朱江曾在1990年对更一般的非线性RLW 方程的周期边值问题提出了一种特征数值方法(特征有限元方法、特征有限差分方法)，理 论分析和数值结果表明，当方程的扩散系数J较小时，这种方法要比一般的有限元和有限 差分方法精确得多(见文献[150．除此之外，新的数值方法还有；谱方法、混合有限元方 法等(见文献[16]，[18】，【19】)． 本文对Rnv方程(包括线性和非线性情况)提出了两种新的有限元方法．即经济型差 分一流线扩散法(Economical Finite Difference Streamline Diffusion Method，简称EFDSD 法)和时间连续特征有限元方法，这两种方法都具有迎风性，分别以不同方式在标准伽略 金有限元中加进了某种类型的稳定化因子，因此，这两种方法比传统的伽略金有限元法更 加稳定．对于第一种方法，我们分别对线性和非线性RLW方程分析了这种方法的稳定性 和误差，最终证明了这种方法是稳定的，并且还得到了基于线性元空闻的L2模拟丰满误 差估计和日1模丰满误差估计．对于第二种方法，我们也分别对线性和非线性RLW方程 进行了误差分析，最终均得到了基于线性元空间的最佳驴、日1模误差估计． §1．2预备知识 定义1．2．1设区域Q C形是Lebesgue非空可测集，L”(Q)是在n上与一个有界函 2 第一章绪论 数几乎处处相等的可测函数全体，对于任意的，∈L。(Q)，定义 11／11。=esssup I，(z)I． 刷。称为，的本性最大模． 定义1．2．2若记D8为函数V的n阶广义导数，则空间w“，”(n)上的范数定义为 ㈣m。(n)5㈣m—n 2高剌D％lloo 定义1．2．3N空闻胪(n)上的内积和范数分别定义如下· (Ⅱ，归Zu·毗”∈巩Q) II乱0={(u．u)){，¨∈L2(Q)． 定义1．2．4 Sobolev空间W”9(n)上的范数z ㈣M叩)=Ilvll。n={∑liD～IIr}；，1 s P。。 其中，当P=2时，Sobolev空间w”,2(f1)记作H”(Q)，其范数记为I| 定义1．2．51a]定义Sobolev空间H1(Q)的子空间 础(Q)=倒∈H1(n)，口lr=o) 其中r表示区域n的边界． 定义1．2．6120l设胃是—个Banach空间，且它的范数记为¨恬，若：【o，卅一日是勒贝格可测的，则定义如下范数t 8训IL。(oF；日)=(／IIv(·，t)112Hdt)§， lIvllL一(o．见舯=sup II口(·，t)llH． uStSl 其中p(o，正H)=如：【o，TI一日： b(o，r㈣∞}，P=2，(30． 定理1．2．1[31悖一不等式)设a和b为实数，则对任意的￡0，有 兰￡口2+去62- 3 第一章绪论 由此，特别地有 №”)l≤㈣俐ell训12+扣”112． 定理1．2．2[211(插值逼近定理)设s^C硪，厶：zrn砩一昂，若口∈H‘n础则有 I[Ihv一训l+hilhv—vlh≤Chll,ll。， l≤s≤r 定理1．2．3嘲(连续形式的Gronwall不等式)设连续函数q(t)(Ⅱ≤t≤b)满足 ㈤I≤卢+口Z‘M刊札。≤t6’ 其中o，p为非负常数，则 I目0)l≤pe。#一d)，a≤t≤b． 定理1．2．4嘲(离散形式的Gronwall不等式)设Ot，卢≥0是任意常数，序列{珊)满足 ‰I≤声+n^EI啦l， 佗=≈，％+l，?，nhsT m=O 其中h0是步长，则 hh『≤ear(卢+akhMo)， n≥k，nh S正 其中Mo=m8x“伽I，lqlI， ，l啦一11) 定理1．2．5㈣(推广的离散GronwaU不等式)设广，扩，胪为定义在J={1，2， ) 上的非负离散网格函数， 驴非减，如果