一类非线性发展方程特征中心差分法.pdfVIP

维普资讯 计 算 物 理 第 24卷 第6期 CHINESEJOURNALOFCOMPUTATIONALPHYSICS Vo1．24．NO．6 2007年 11月 Nov．．2007 [文章编号】 1001—246X(2007J06—0637—10 一 类非线性发展方程的特征中心差分法 郭双冰， 张志跃 (南京师范大学数学与计算机科学学 院，江苏 南京 210097) [摘 要] 给 出一类非线性发展方程 的特征 中心差分法 ，分别得到非规则网格上的位移 “、速度 “及其对空间 变量 的一阶导数项的差分解和误差估计 ．所讨论方法的计算量与基于线性插值 的特征差分法相当，其近似解 与基于二次插值 的特征差分法的近似解具有相同阶的误差估计 ， “， “对空间变量 的一阶导数近似均具有超收 敛误差估计 ．数值试验说 明了该方法的可行性和有效性 ． [关键词] 非线性发展方程；特征中心差分法；误差估计 [中图分类号] O241 [文献标识码] A O 引言 本文首先考虑非线性发展方程的初边值问题 u+6()一 一等： +G( ∈Q：n×EO,T3， M：M ：O， ∈an ×Eo，]， (1) “(，O)= ()，“(，O)= ()， ∈ ． 其中n=[O，z]，an表示n的边界 ， ()， ()为已知光滑函数 ． 非线性发展方程 (1)在生物学和力学等许多数学物理领域有着广泛的应用，如非线性弹性杆中纵向形变 波传播 ，弱非线性作用下空间变换离子声波传播及神经传播等 问题 ，具有深刻的物理背景 ¨。．文 [3～11]在 抛物型方程的有限差分方法的数值分析方面均取得 了一定的结果 ． 本文讨论上述非线性发展方程的特征中心差分法 ．文[12]提出求解对流扩散方程的特征差分法；文[13] 提出解线性椭圆型和线性抛物型方程的块中心差分法；文 [14]讨论 了解拟线性双 曲型积微分方程的块中心 差分法；文[15]综合特征法和块中心差分法，讨论了具有混合边值条件的非线性对流 一扩散方程的特征块中 心一差分法 ．文 [16]讨论了一维非线性对流 占优扩散方程的变网格特征差分方法 ．在此基础上，本文在非规 则网格上构造了上述非线性发展方程的特征中心差分法 ，同时恰 当使用一种变换 ，同时得到了 u和u 的数 值解 ．由于 u在实际中也是重要的物理量 ，这样可以避免采用其他格式先计算 u，然后再计算 u而造成二 次误差的积累和产生 ．此外，分别得到 u，u 的近似梯度，从而使得整个工作的计算量小于混合有限元方法 ． 本方法的计算量与基于线性插值的特征差分法相 当，其近似解的误差阶与基于二次插值的特征差分法的近 似解具有相同阶误差估计 ，u，u对空间变量 的一阶导数近似具有超收敛误差估计 ． ． ． 设 =u，用 r表示与算式 +b(，t) 有关的特征方 向 0 a 1 fIV ，， 、a＼ t DL，， J， 其中 ≯(，t)=(1十b(，t))“，则问题 (1)能写成 [收稿 日期】2006—07—03；[修 回日