改进余弦微分求积法与紧致差分方法在隐式和类隐式格式中的运用.docVIP

第一章绪论 计算物理学的背景 纵观物理学的发展史，大致可以分为古代、近代与现代三个阶段。在古代时期， 物理学的发展主要是基于一些哲学思考，没有形成特别成熟的体系，发展比较缓慢， 主要原因是欧洲的黑暗教皇统治，直至文艺复兴时期才有所改善。 16世纪至19世纪，这一时期称为近代物理学时期，或是经典物理学时期，这一 时期是经典物理学的建立时期，从哥白尼的日心说，到牛顿的三大力学定律，再到麦 克斯韦的电磁学方程组，物理学有了长足的发展。在这一时期，物理学的体系大致形 成，主要分为理论物理与实验物理两大分支，物理学的发展大致可以看做是实验物理 与理论物理相互推进的过程，人们通过物理实验现象来构造基础理论，再通过理论分 析计算来预测实验现象。在理论分析计算过程中往往会有难以求解的非线性方程，然 而在那一时期对此并没有很好的解决方法。 自19世纪开始，量子力学与相对论的诞生与发展使得近代物理学理论有了长足 的发展，然而不同于牛顿第二定律简洁的形式，量子力学中与之对应的薛定谒方程本 身就以偏微分方程的形式给出，这也就使得理论计算越发地困难，在物理学的工程运 用方面对理论分析计算的要求也越来越高。比如二战时期，美国研制核武器的过程， 就涉及了流体动力学过程、核反应过程、中子输运过程等等，这些都涉及到对十分复 杂的偏微分方程的求解，而精确求解这些偏微分方程儿乎是不可能的。因而数值求解 偏微分方程就成为一种相对比较容易实现的方法，而这一时期计算机技术的飞速发 展，也使得这些大型计算变得可能，所以计算物理学应运而生，它主要研究数值求解 各种物理体系或者物理模型的方法。然而随着计算物理学的发展与实际应用的要求， 数值计算变得不仅仅只是理论计算的一种手段，同时也变成了实验物理的一种手段， 由于某些试验要求的规模太大（如大气流动），或是数据难以测量（如液体流动时各 点液体的速度），或是试验比较危险不能直接尝试（如核聚变或裂变实验），数值模拟 实验就成为一个很好的选择，由于数值模拟实验从根本上说也是数值求解某些物理模 型，所以自然也就归入了计算物理学的范畴。 总体来说，计算物理学是一门新兴的交叉学科，其主要研究领域是运用现代计算 机的高速运算功能来实现物理学以及工程物理学中的复杂模型模拟及运算，是当代物 理学科一个重要分支。计算物理学发展至今已经诞牛了很多实用可靠的数值计算方 法，下面就常用的数值方法做一简要介绍。 有限差分方法其主要思想是将连续的空间离散成分立的网格点，这些网格点称 作节点，然后将定义在连续空间上的函数离散成分立的变量函数。把待求解的方程中 的微分用变量函数的差来表示，积分用变量函数的求和来表示，这样待求解的方程就 可简化为线性代数方程组，该方程组称为有限差分方程组，解此方程组就可获得原方 程在节点处的近似解。有限差分方法是数值计算中最肓接简单的方法，其基木思想是 空间离散化，也是绝大多数数值方法的基木思想。有限差分方法格式简单，运用灵活, 实用性强，在计算物理学中占有很重要的位置。 有限元方法此方法诞生于20世纪中期，最早用于结构力学，后来由于计算机 技术的发展，开始用于流体力学方面的模拟。该方法的基木思想是将方程的求解区域 划分成互不重叠而且相连的单元，在每个单元屮选择基函数，并用基函数的线性组合 来构造单元中方程的近似解，并将求解区域中所有单元近似解的集合看作是方程的近 似解。该方法单元划分灵活，基函数选择多样，且对求解不规则边界问题很有优越性。 但其数学原理复杂，编程复杂，计算量大。 谱方法该方法的主要思想是，使用一组光滑的基函数的线性组合来构造某一方 程的近似解，其线性组合的系数由方程木身决定，类似于傅里叶分析的基木思想。该 方法主要用于流体力学的计算与模拟。该方法在求解规则区域时有很好的精度，但灵 活性差，对于一些不规则区域求解结果不好。 上述方法是数值计算中比较常用的三种，其他还有很多实用的方法，如微分求积 法、蒙特卡洛方法、分子动力学方法和大涡模拟法等等。 「2微分求积方法简介 微分求积方法(Differential Quadrature Method, DQM)是一种新的求解偏微 分方程的数值方法，于1971年由Bellman和Casti提出⑴,至今已经被成功地运 用于许多实际的工程物理问题中。微分求积方法使用方便，适用范围广，计算精度 高，计算速度快，边界条件处理简单，且不依赖于泛函与变分原理，是一种十分有 效的方法。 近年来微分求积方法的理论研究与应用取得了很人的进展。美国堪萨斯州立大 学的Mingle使用该方法求解了非线性扩散方程；美国Oklahoma大学的小组应用广 义DQ方法处理输运过程、泊松方程以及多维问题等等35］；随着微分求积法的发展， 又先后被运用于求解纳维斯托克斯方程⑹、薛定铐方程⑺等。微分求积法在工程