第三章－数值计算2008.08.pptVIP

第三章 离散化方法 主要包含以下内容： §3.1 数值方法的本质 §3.2 推导离散化方程的方法 §3.3 一个说明性的例子 §3.4 四项基本法则 §3.5 结 语 §3.1 数值方法的本质 §3.1-1 任务 1. 什么是微分方程的数值解？ 它是由一组可以构成因变量?分布的数组成的集合，即用一组数字表示待定变量在定义域内的分布。类似于在实验室中进行的实验，仪器的读数构成了所研究区域内被测物理量的分布（有限个离散点的值的集合）。 2. 实验数据的处理—拟合关系式 通常采用的方法是利用实验结果，拟合出遵循的方程，然后利用方程即可计算出任意位置点上因变量的值。 * 太 原 理 工 大 学 * /38 Thermal 本章任务：前两章已阐述了获得物理现象的理论预测有一些明显的好处，且一些有趣的现象是受一些微分方程支配的，这些方程已用通用方程归纳了，接下来的任务就是推导求解这个通用方程的方法—即所谓的离散化方法。 第三章 离散化方法 推导前的假设： 设通用变量? 仅仅是一个变量的函数（一维问题） eg. 设? 随 x 的变化遵循高次多项式的形式 利用实验数据，并采用数值方法求得有限数量的多项式中的各项系数值，进而得到相应的?与x 遵循的方程式，由此即可求得任意点处的? 值了。 若最终的兴趣是得到不同位置上的?值，则此方法有些不便，因为各个系数 a 本身没有什么特别的意义，但要求得 ? 还必须进行代入过程。 3. 数值方法及任务 建立一个把一系列给定点上的? 值作为原始未知量的方程，实际上求解微分方程的多数方法均属此类。该方法的任务是提供一组关于这些未知量的代数方程并规定求解这组方程的算法。 §3.1-2 离散化的概念 1. 离散化方法及基本思想 把注意力集中在网格节点处的值，用离散的值取代包含在微分方程精确解中的连续信息，这样就离散了? 的分布，这类数值方法叫离散化方法。 根据实际研究对象，把定义域分为若干个有限的区域，在定义域内连续变化的待求变量场，由有限区域上的若干个点的待求变量值来表示，这就是离散化的基本思想。 2. 离散化方程：所取网格节点上未知因变量 ? 值的代数方程，此方程由支配 ? 的微分方程推导而得。 在推导过程中，需对网格节点之间? 如何变化作某种假设，变量? 在节点间的分布形式不同，推导离散化方程的方法也就不同；另外可以选择在整个计算域内满足一个简单表达式的分布；更为实际的方法还是采用分段分布，即将计算区域分布一定数量的子域或单元，每个子域可以有一个独立的分布假设。 §3.1-3 离散化方程的结构 1. 离散化方程的结构 一个离散化方程是连接一组网格节点处? 值的代数关系式，由支配? 的微分方程推导而得，并表示与该微分方程相同的物理信息。当节点数很多时，离散方程的解接近于相应微分方程的精确解，相邻点之间? 变化很小，有关? 分段分布的细节就不那么重要了。 相应于一个已知的微分方程，离散化方程的形式决不是唯一的，这起因于分布假设以及推导方法的不同。 网格节点数非常多的极限条件下，所有可能类型的离散化方程将会给出相同的解。 2. 离散化方法 常见的方法主要有：有限差分法和有限元法。 有限差分法(finite difference method, FDM) 基本点：将求解区域用与坐标轴平行的一系列网格线的交点所组成的点的集合来代替，在每个节点上，将控制方程中每一个导数用相应的差分表达式来代替，从而在每个节点上形成一个代数方程，每个方程中包含了本节点及其附近一些节点上的未知值，求解这些代数方程就获得了所需的数值解。 有限容积法(finite volume method, FVM) 基本点：将所求解的区域划分成一系列控制容积，每个控制容积都有一个节点作代表。通过将守恒性的控制方程对控制容积做积分来导出离散方程。 在推导过程中，需对网格节点之间? 如何变化作某种假设，变量? 在节点间的分布形式不同，推导离散化方程的方法也就不同；另外可以选择在整个计算域内满足一个简单表达式的分布；更为实际的方法还是采用分段分布，即将计算区域分布一定数量的子域或单元，每个子域可以有一个独立的分布假设。 用控制容积法导出的离散方程可以保证具有守恒特性，而且离散方程系数的物理意义明确，是目前流动与传热问题数值计算中应用最广的一种方法 两种方法的区别来自于选择分布和推导离散化方程的方法不同。 在此主要关注的方法具有有限差分的外形，但它采用了典型的有限元方法所具有的思想，把此方法叫有限差分法可能在于它坚持遵守习惯的有限差分法做法。 §3.2 推导离散化方程的方法 对于一个已知