Rosenau-Burgers方程一种新的数值解法.pdfVIP

第28卷 第5期 工 程 数 学 学 报 v01．28No．5 20U年10月 CHINESE JOURNALOFENGINEERINGMATHEMATICS Oct．2011 文章编号：1005-3085(2011)05—0665—06 Rosenau—Burgers方程一种新的数值解法术 刘 艳， 阿布都热西提 ．阿布都外力 (新疆大学数学与系统科学学院，乌鲁木齐 830046) 摘 要：本文讨论Rosenau—Burgers方程初边值问题的数值解法．针对Rosenau．Burgers方程构造了一个 新的差分格式，把网格分为奇、偶两套独立的网格，在偶数网格点采用显式格式，在奇数网格点 采用Crank-Nicolson格式，这样偶、奇、显、隐交替的方法使计算量减少．同时针对非线性项 进行了线性化，使格式的近似解更精确．给出了稳定性和收敛性的严格理论证明，数值实验结果 表明了理论证明的正确性及格式的有效性和可行性，具有推广价值． 关键词：Rosenau．Burgers方程；非线性项；收敛性；稳定性 分类号：AMS(2000)35G25 中图分类号：O241．82 文献标识码：A 1 弓言 在离散动力学系统的研究中，对于波与波相撞或波与墙相撞这一物理特点，用我们所熟知 的KdV方程不能给出很好的描述，为了克服这个缺点，Rosenau[1，2]提出了Rosenau方程 + 。 t+ +UU =0，X∈[0， ，t∈[0，]． (1) 方程 f1)解的存在性和唯一性Park[3】已给出证明．但其解析解的寻找是困难的，另一方面，为 了更进一步考虑动力学系统中的空间耗散，例如水波等，方程 (1)必须加上粘性项一乱 ，故方 程为 t+钆 t一乱 + 。+ =0， ∈[0，]，t∈[0，】． (2) 由于 (2)的耗散影响等同于如下Burgers方程 t—u。+ + =0， ∈[0， ，t∈[0，】， 所以(2)被称为Rosenau—Burgers方程．我们提出在边值 u(0，t)7--u(，t)=0，u。。(0，t)：u。。(L，t)=0，t∈0【，T】， 乱 (0，t)= (Lt)=0， t∈[0， ， (3) 和初值 (z，0)=uo(x)，X∈0【，L】 (4) 条件下的Rosenau—Burgers方程的新的差分格式．方程 (2)柯西问题解的渐进性质，特别是 行波解的稳定性在文献 [4-6]中已经给出．用Crank—Nicolson方法求解方程 (2)，(3)，(4)在文 献 7【】中已经给出．本文的第二部分将在文献 7【]的基础上把网格分为奇偶两套进行差分格式的 构造，第三部分对构造的差分格式进行其稳定性和收敛性分析，第四部分将给出数值例子验证 理论分析的正确性． 收稿 日期：2009-10—30．作者简介：刘艳 (1983年6月生)，女，硕士．研究方向：偏微分方程数值解 基金项目：国家自然科学基金 10961024)；新疆高校科研计划 (xJEDu2007102)． 666 工 程 数 学 学 报 第28卷 2 差分格式的构造 设h和丁分别为空间步长和时间步长，记XJ=jh(0 J J)，t=礼丁(0 7 Ⅳ)，zo= {11,=(uJ)fIt0=uJ=0，J=0，1，2，… ，l，)，定义差分算子如下 ()=( 一uj)／h， (n)虿=(哼一呼 )／， (n)奢=(略 一乱 ~)／2h (u)：(+1一u)／下， (Ⅱ) =(，“