在Crank-Nicolson格式下解对流扩散问题的移动网格方法.pdfVIP

数学物理学报 2011，31A(4)：1063—1070 http：／／actams．wipm．ac．cn 在 Crank．Nicolson格式下解对流扩散问题 的移动网格方法 杨银 周琴 (湘潭大学数学与计算科学学院 湖南湘潭411105；。湖南涉外经济学院理学部 长沙410205) 摘要：对一类含源项非定常奇异摄动对流扩散问题，作者构造了移动网格下的Crank—Nicolson 差分格式，介绍了移动网格的网格迭代方法和算法，并通过数值实验实现了算法，将求解的误 差与均匀网格下求解的误差相比较，证明了在移动网格下运用 Crank—Nicolson差分格式求解 这类奇异摄动问题能有效的克服数值振荡，得到较好的结果． 关键词：奇异摄动；对流扩散；移动网格； Crank—Nicolson差分格式． MR(2000)主题分类：65Mll；65M12 中图分类号：O211．4 文献标识码：A 文章编号：1003—3998(2011)04—1063—08 1 引言 f{+0()=E象+，( u(x，0)=西()， 0 L， (1) If上(0，t)= 1@)，r“(L，t)= 2@)， 0 t T． 初值条件函数 ()， 1()， 2()已知．其中0E《1是扩散系数，未知量 u(x，t)是关于空 间X和时间t的函数，简记为饥这是含源项非定常对流扩散模型问题，在具有污染源的大 气、地下水环境及化工流动领域均有应用．函数 a(x，t)和 f(x，t)为已知的对流速度和源汇 函数，设它们满足一定的正则性条件．当a=f“， 为运动黏性系数时，模型方程为Burgers 方程，通常视为流体流动的Navier—Stokes方程组的模型方程． 对于含源项非定常对流扩散模型问题 (1)，文献 [4，6]分别提出了数值求解的高精度紧致 隐式差分方法．文献 f5，71研究了求解对流扩散问题的差分格式．但这些方法在求解方程时 空间和时间均采用均匀网格，不能有效的解决当 ￡很小时的奇异摄动情况．移动网格方法是 求解奇异摄动问题的有效方法之一，它能将网格点集中到解的变化比较剧烈的区域．对于奇 收稿 日期：2009—12—22；修订 日期：2010—10-21 E—m ail：yangyinxtu@xtu．edu．CYl 基金项 目：国家自然科学基金 、湖南省教育厅项目(10Al17，09C943)和湘潭大学校级项目 (10QDZlO)资助 1064 数 学 物 理 学 报 VOl1．3lA 异摄动两点边值问题，文献 [12]在移动网格下进行了理论分析，在半离散差分格式下，利 用离散格林函数等得出数值解一致收敛的结论．而我们将移动网格方法应用到含源项非定 常对流扩散问题，将其与 Crank—Nicolson差分格式相结合来求解，使得求得的数值解能更好 的逼近实际问题．并且，我们通过数值实验证明，在移动网格下运用Crank—Nicolson格式求 解，能有效的克服数值振荡，并能解决均匀网格下很难有效解决的奇异摄动问题． 2 差分格式 D U7 我们将模型方程 (1)第一式变形为 】【 =一￡ ) +Ou 圳 ， 。 。