计算区域与控制方程的离散.pptVIP

流体流动与传热的数值计算 第二节 建立离散方程的Taylor展开法 三、一维模型方程的离散 第二节 建立离散方程的Taylor展开法 三、一维模型方程的离散 Crank-Nicolson格式： 流体流动与传热的数值计算 第二章 计算区域与控制方程的离散化 * 复习与总结 所有流动与传热过程都受到最基本的3(4)个物理规律支配： 质量守恒； 动量守恒； 能量守恒； 组分守恒。 守恒定律的数学表达式： 偏微分方程 (控制方程) 质量守恒方程 动量守恒 (Navier-Stokes) 方程 对于粘度系数为常数的不可压缩流体： 能量守恒方程 化简后的能量守恒方程： 组分守恒方程 封闭方程 上述 7 个方程中包含 7 个未知量：ρ，u，v，w，p，T，ml，方程封闭。 控制方程的通用形式： 流动与传热问题的控制方程可以表达为通用形式： 第二章 计算区域与控制方程的离散化 对流与传热问题数值计算的步骤： 网格生成 (区域的离散) 对空间区域进行剖分，将它划分为许多子区域，并确定每个区域中的节点。 控制方程的离散化 (方程的离散) 生成网格后，将描写流动与传热过程的偏微分方程转化为各个节点上的代数方程。 代数方程组的求解 (方程的求解) 建立数学模型(控制方程及定解条件) 第一节 空间区域的离散化 一、区域离散化的内容 一、区域离散化的内容 1. 区域离散化的定义 区域离散化：用一组有限个离散点来代替原来 的连续空间。 第一节 空间区域的离散化 2. 区域离散化的步骤 把计算区域划分为许多互不重叠的子区域； 确定每个子区域的节点位置 确定每个节点所代表的控制容积。 第一节 空间区域的离散化 一、区域离散化的内容 3. 区域离散化的几何要素 节点： 需要求解的未知物理量的几何位置； 控制容积： 应用控制方程的最小的几何单位； 界面： 各控制容积的分界面； 网格线： 沿坐标方向联结相邻节点而成的曲线簇。 第一节 空间区域的离散化 一、区域离散化的内容 4. 说明 一般将节点看成为控制容积的代表； 控制容积 与子区域 并不总是重合的(视节点位置而定)； 子区域：区域离散化过程中，由一系列与坐标轴相应的直线或曲线簇所划分出来的小区域。 在不同的离散方法时，节点在子区域中的位置不同。 第一节 空间区域的离散化 二、设置节点的方法 二、设置节点的方法 1. 外点法(FDM) 节点位于子区域的角顶上，划分子区域的曲线簇就是网格线，但子区域不是控制容积。 在相邻两节点的中间位置做界面线，由这些界面线构成各节点的控制容积。 第一节 空间区域的离散化 二、设置节点的方法 2. 内点法(FVM) 节点位于子区域中心。 子区域即为控制容积 第一节 空间区域的离散化 三、一般规定 三、几个规定 1. 实线表示网格线，虚线表示界面线，黑点表示节点。 第一节 空间区域的离散化 三、一般规定 2. 对节点及有关要素的命名： (1) 对离散方程进行特性分析时： 采用 i-j-n 表示法 所研究的问题的节点位置为 (i, j)； 与该节点相邻的界面分别为： i+1/2, i-1/2, j+1/2, j-1/2 n 表示非稳态问题的时间层。 第一节 空间区域的离散化 三、一般规定 (2) 在进行其他分析时： 采用P, N, E, W, S表示所研究的节点及相邻的4个节点； 用n, e, w, s 表示相应的界面； 用上标0表示非稳态问题中上一时间层的值； 相邻两节点间的距离用δx 表示； 相邻两界面间的距离用Δx 表示。 说明：在均匀的网格系统中，不强调δx与Δx的区别。 第一节 空间区域的离散化 四、两类节点的比较 四、内点法与外点法的比较 相同点： 当网格划分均匀时，在区域内部，内点法和外点法的节点没有太大区别，仅是控制容积有半个控制容积厚度的错位。 第一节 空间区域的离散化