对流扩散方程-一致收敛差分格式-计算数学专业论文.docxVIP

万方数据 万方数据 独 创 性 声 明 本人声明所呈交的学位论文是本人在导师指导下进行的研究工作 及取得的研究成果。据我所知，除了文中特别加以标注和致谢的地方 外，论文中不包含其他人已经发表或撰写过的研究成果，也不包含为 获得电子科技大学或其它教育机构的学位或证书而使用过的材料。与 我一同工作的同志对本研究所做的任何贡献均已在论文中作了明确的 说明并表示谢意。 签名： 日期： 年 月 日 关于论文使用授权的说明 本学位论文作者完全了解电子科技大学有关保留、使用学位论文 的规定，有权保留并向国家有关部门或机构送交论文的复印件和磁盘， 允许论文被查阅和借阅。本人授权电子科技大学可以将学位论文的全 部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描 等复制手段保存、汇编学位论文。 （保密的学位论文在解密后应遵守此规定） 签名： 导师签名: 日期： 年 月 日 电子科技大学硕士学位论文 摘 要 对流扩散方程，是用来描述黏性流体的非线性方程的线性化模型方程，是一 类基本的运动方程。对流扩散方程对偏微分方程数值计算方法的研究具有重要的 理论和实际意义。 数值边界问题一直是偏微分方程数值解的难题。求解数值方法有很多种，如 有限差分法、有限元法等。用传统的有限差分法和 Galerkin 有限元法求解往往会 产生数值振荡。格林函数可以很好有效的处理偏微分方程边界问题，避免数值振 荡，而且还提高了算法的计算效率。而且在求解多变量偏微分方程情况下，传统 有限元和有限差分方法很难推出很好的解析格式。格林函数可把求解过程归结为 一个特定的边值问题，对于一些特定的边值可以得到具体的表达式。格林函数可 以很好的处理多变量偏微分方程。??函数的一些特殊性质可以解决一些突出的边界 问题，从而简化推导过程。 在一维对流扩散方程情况下，把对流项和扩散项作为一个算子 M ，推出算子 的局部格林函数，把格林函数和对流扩散方程相互作用积分，通过分部积分法推 出一致收敛差分格式，并对差分格式做误差分析。采用经典算法追赶法进行数值 求解，数值实验验证了格式的有效性。然后对一致收敛差分格式做了稳定性分析， 收敛因子是指数级。 在二维对流扩散方程情况下，根据一维情况推出二维情况的积分方程。通过 格林公式和格林函数的性质推出了五点差分格式和九点差分格式。对源函数应用 二维拉格朗日插值法，推出差分格式的计算模板，对差分格式做误差分析。采用 经典算法 BICG 和 BICGSTAB 进行数值求解，数值实验验证了差分格式的有效性。从 数值实验看出一致收敛差分格式运行时间短，迭代次数少，误差低，而且简单， 从而可以继续推广到 n 维情况或者应用到其它的方面。 关键词：对流扩散方程，局部 Green 函数，Galerkin 原理，分部积分法 IV ABSTRACT ABSTRACT Numerical boundary problems are generally acknowledged to be a hard task for PDE. There are many methods to deal with it, for example, finite difference method, finite element method etc. It is well known that if the standard finite difference method or finite element method is used in regions where layers occur, then unphysical oscillations arise. Green’s function method for the numerical solution of sharp boundary problem has high efficiency to avoid numerical oscillations and improve efficiency of algorithm compute. However, in several variables PDE, used as standard finite difference method and finite element method are hard to derive in an analytical scheme. Green function can reduce the solution procedure to specific boundary problems, and it can get the exact expression. Green’s function can d