傳热与流体流动的数值计算13章.pptVIP

传热与流体流动的数值计算;本课程学习内容;参 考 书 目;第一章 引 论;实验研究 全比例实验 模型实验 结果外推 测量仪表精度 理论计算 理论预测出自于数学模型的结果 数学模型主要由一组微分方程组成 解析解和数值解 ;1-2 预测的方法;2-1 控制微分方程 微分方程的意义 各个微分方程都代表着一定的守恒原理。每一个方程以一定的物理量做为它的因变量，方程本身则代表着那些影响该因变量的各个因素之间必定存在着的某种平衡。 通常以单位质量为基础来表示各因变量。 例如：质量分量、速度（单位质量的动量）、比焓等。;假设J代表一个典型因变量?的流量密度;以单位容积为基础来表达一项——变化速率 代表在单位容积内所包含的相应广延性质的大小 表示单位容积内有关性质的变化率;— 化学组分的守恒;— 能量方程;— 动量方程;— 紊流的时间平均方程;— 紊流的动能方程;— 通用微分方程;? 可以代表无因次的变量 热、质传递，流体流动，紊流以及有关的一些现象的所有有关微分方程都可以看成通用方程的一个特殊情况；可以只编写一个求解通用方程的程序，对不同意义的? 重复使用这个程序； 对不同的? 需要对相应的?和S分别赋以各自合适的表达式，同时给出合适的初始条件和边界条件。;自变量 一般来说，因变量φ是三个空间坐标与时间的函数 φ =φ（x，y，z，t） 其中x，y，z以及t都是自变量。 当有关的物理量只与一个空间坐标有关时，所研究的问题是一维的；当问题与时间无关时，叫做稳态的，否则叫做非稳态或与时间有关的问题。 另一种写法： z=z（x，y，T） z是因变量，代表在位置（x，y）相对于温度T的等温面高度。;恰当明智地选择坐标系统有时可以减少所需要的自变量数。 并非只能使用直角坐标系，任何一种描述空间位置的方式都是可以采用的。 例子： 1. 在一个静止的坐标系上看以恒定速度飞行的飞机周围的流体流动是非稳态的；但是相对于固定在飞机上的移动坐标系而言，流动是稳态的。 2. 在一圆管内的轴对称流动于直角坐标系内是三维的，但在r，θ，z的圆柱极坐标系内则是二维的。 3. 坐标变??可能用来进一步减少自变量数量。 4. 改变因变量可能导致自变量数目的减少。;如果在一个坐标上的一个给定位置处的条件，要受该位置两侧条件变化的影响的话，那么这个坐标就是一个双向的坐标。 如果在一个坐标上的一个给定位置处的条件只受该位置一侧条件变化的影响，这样的坐标就是一个单向的坐标。 空间坐标是双向坐标，但如果在一个坐标方向上有很强的单向流动，也可以近似作为单向坐标。对流——单向，扩散——双向 时间坐标是单向坐标。 抛物型表示一种单向的状态；椭圆型表示双向的概念。 非稳态导热问题实际上是时间坐标上的抛物型和空间坐标上的椭圆型问题；稳态导热对所有的坐标都是椭圆型的。 讨论单向双向坐标的动机：如果可以用一个单向的坐标来规定一个给定的状态，就有可能大大节省计算机的存储量和时间。;3-1 数值方法的本质 一个微分方程的数值解系由一组可以构成因变量?的分布的数所组成 类似于在实验室中进行实验 定义 数值方法就是把计算域内有限数量位置（叫做网格结点）上的因变量值当作为基本的未知量来处理。 任务 任务是提供一组关于这些未知量的代数方程并规定求解这组方程的算法。;离散化的概念;离散化方程的结构;泰勒级数公式 如图中网格结点。结点2为结点1，3的中点，在2周围展开泰勒级数：;- 变分公式;- 加权余数法;控制容积公式;优点：所得到的结果将意味着任何一组的控制容积内，也就是整个计算域内，诸如质量、动量以及能量这样一些物理量的积分守恒都可以精确地得到满足。 无论网格数量多少，这种积分守恒都满足。 求解离散化方程以得到在网格结点上的因变量值时，可以用两种不同方法来表示结果。 有限元法和大多数加权余数法中把假设的由网格结点上的值以及网格结点之间的内插函数（或分布）所构成的φ的变化取作近似解。 有限差分法只考虑把网格结点上的φ值构成解，而不去管在网格结点之间φ是怎样变化。类似于实验。 对微分方程中的不同项可以采用不同的分布假设来进行积分。;其中k是导热系数，T是温度，S是单位容积的发热率。使用下图的网格结点群：; 两种分布假设： 阶梯形分布——斜率 在控制容积面上是不确 定的； 线性分布——网格结 点之间采用线性的内插 函数。;用分段线性分布计算方程（3.11）中的dT/dx，所得方程：;上述有关选择分布函数的自由度，最终导致不同变型的离散化方程形式。当网格结点的数目增加时，所有这些不同形式的方程都会给出相同的解。 附加要求（可大大减少可以接受的公式的数目）： 即便是采用很粗的网格，解也总应该满足： 物理上真实的性状 一个真实的变化应当具有与准确