-传热学数值计算方法.pptVIP

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