对流占优的特征有限元方法-计算数学专业论文.docxVIP

原告j 性声明 本人郑重声明: rfr 呈交的学位法文，是本人在导 yf!? :j:旨导下，进行研究工作听取萍的或 黯文中己经注明引黑的内蒋外，本学位段文不包含任何他人或集体 i兰经发表的作品内 ，也不包含本人为我得其他学位而使用过的材料。对本论文所带及的研究:川牛做出贡献的 其他个人或集体， J句己在文巾以明确方式标明。本学位论文国创性声明的活悻责任由本人承 担。 喻的作者签名:级d在 日期: :J()JO. 1:1.I 7 关于学位论文版权使用授权的说明 本人完全了解离北工业大学关于收集、保存、使用学位论文的以下规定: 影邓、缩萍、扫描、数字化或其它手段保存诠文:学技有权提供本学位论文含义或者部分内 容的 IRJ宽摄务;学校有较将学位论文的全部或部分内容辑人有关数据撑进行检索、交谊:学 校有投向皇家有关部门或者载，持滋交;在文的复印件和电子截。 (保密位学位语文在挥委道培养;授权说明〉 叫作者鉴名:做 导师签名:如泣 巳现: 2t7/0. /:J. 1 日期1 划。 IJ. 7 河北 河北工业大学硕士学位论文 对流 对流占优的特征有限元方法 PAGE PAGE 11 PAGE PAGE 10 第一章 综述 §1-1 问题的来源与研究背景 包含时间变量的很多数学物理偏微分方程统称为发展方程[1](evolution equation)，又可 称作进化方程或演化方程。这类方程往往用来描述在物理、力学以及其它自然科学中随着时 间而变化的状态或者过程，例如热传导方程、流体与气体力学方程组、反应扩散与对流扩散 方程、S chro¨dinger 方程、Kdv方程等等和由这些方程利用恰当方式相合而得到的相合方程 组。各种发展方程的求解问题在科学与技术的发展中得以提出，然而存在绝大多数情况，使 得这些问题不能通过利用解析的公式来表达问题的解，或者解析的公式过于复杂，从而利用 数值求解并描述随时间而变化的状态或过程的偏微分方程的研究具有很强的理论性和应用 性。同时电子计算机以及计算技术的迅猛发展正强有力地推动着该学科以前所未有的速度 蓬勃发展。研究发展方程数值解产生了两大基本方法有限差分方法和有限元方法。 §1-2 研究现状 该部分主要介绍有限元方法的概况及对流占优扩散问题当前的研究现状。 1-2-1 有限元方法 有限元方法在研究偏微分方程数值解时极为有效。在20世纪50年代中期，为了适应航空 工业、水利工程的需要和计算机技术的发展，在工业界将结构体分为有限个单元，再利用变 分原理求解弹性力学和结构力学问题，构成了一类新的算法。在60年代以后，以冯康先生为 代表的中国学者与西方学者独立并行地完成了这种系统的数值方法，并奠定了其数学基础。 从历史上看，Courant在19年己经提出过一种方法在三角形网格上用逐片线性函数去逼 近Dirichlet问题。从数的观点上看，有限元法本质上可看作古典的变分近似法，诸如提出的 一类数值方法即Ritz法和Galerkin[2,3]法的推广。古典的变分近似法很难构造出满足本质值 条件的基函数，在构造容许函数类的有限维线性集合(或空间)时，特别是形状复杂的区域更 难以构造。古典的变分近似法成为实用的方法，是通过有限元法将构造基函数变成一种规范 的过程来得以实现的。由微分方程所对应的泛函，利用变分原理同时结合区域剖分得到离散 算式即代数方程组，作为其重要归化途径。这种方法对于各种形状的边界都能够灵活处理， 克服了有限差分法对区域形状的限制，具有更好的几何适应性，该法有易于编制大型通用计 算软件的优点。并且一系列理论为有限元法的理论研究提供了基础，诸如泛函分析、变分方 法、Sobolev空间的理论，以及椭圆边值问题的现代理论。 有限元法求解数学物理问题的过程如下:(1)利用某种变分原理，将原问题化为变分问题， 由于同一问题可有不同的变分形式，实际中用适合于原问题性质及要求的变分原理，利用变 分近似法求问题的广义解。(2)对求解区域Q作剖分，遵循剖分原则，剖分方式与分片偏差值 相协调。(3)在单元上采用分片括值的方法构造满足本质边值条件的试探函数类和检验函数 空间，使得括值函数在不同单元间保持光滑性连接。(4)作单元分析和总体合成。(5)处理本质 边值条件，使在边界上的网格点近似解满足本质边值条件，得到一个线性代数方程组。(6)求 解所得代数方程组Ku=f。其系数矩阵K既具有稀疏性，又可保持原问题的一些特性(如对称 性，正定性等)便于求解。(7)从Ku=f的解，通过试探函数空间的基函数，得到变分问题的近 似解，即为原问题的有限元解。 1-2-2 对流扩散型方程，特征有限元方法以及先验误差估计 对流扩散型偏微分方程(或方程组)具有共同的特点:对