第一章 区域离散及建立高散.pptVIP

区域离散及建立高散方程的方法 §1-1传热问题控制方程的类型 常见形式：二元二阶线性偏微分方程。 数学分类：P点 §1-2 空间区域的离散化方法 §1-3 控制方程的离散方法 * * 第一章 椭圆 过P点无实特征线 抛物 过P点有一条实特征线 双曲 过P点有两条实特征线 依赖区：P点的依赖区是指要唯一确定P点值，必须完全 给定条件的区。 影响区：P点值的变化会影响该区各点的值。 物理分类: 椭圆方程。 a. 平衡、稳态问题于时间无关。 b. 依赖区：包围求解域的边界的封闭曲线，故 称边值问题。 c. 影响区：整个求解域。 d. 计算：控制方程联立求解。 e. 例如：稳定导热、流动。 P.影响 区 依赖区 1. a. 步进、非稳态问题。 b. 依赖区：P点上游（时间或流动上游） 区域的边界。 依赖 P 影响 c. 影响区：P点下游。 d. 计算：逐步推进，一个时刻一个时刻或一个截面一个 截面的计算 e. 例如：一维非稳态导热，二维稳态边界层。 a. 步进问题 b. 依赖区：上游 c. 影响区：下游 在两条特征线之间 依赖 影响 抛物型方程： 1. 双曲型方程 2. 双向坐标→椭圆问题→空间坐标→回流问题 单向坐标→抛物问题→时间坐标→边界层问题 （边界层主流向几何坐标） 方程的守恒型与非守恒型 在直角坐标常物理问题中，对流项为散度形式，则该方程是守恒型。 从微分角度现有方程都是质量、动量、能量守恒的。 1. 守恒型方程，由扩散和对流作用造成的总转移量在任 何有限大小的容积中都守恒。 2. 非守恒型方程 3. 离散：将计算域划分为许多子区域。 四种几何要素： 两类离散方法： 界面 子区域 控制容积 节点 网格线 节点 控制容积 界面 网格线 1. 2. 3. 4. 外节点法：节点在子域的四 角，先定节点 位置而计算相应的界面 1. 两种标记法： W N P E S w e n s 内节点法：节点在子域中心，子域与控制容积重合。 计算时先定界面后算出节点位置。 2. 在边界上完整的控制容积 或者最外边是一个厚度为零的控制容积 内节点 外节点 在边界上是半个控制容积 两种离散方法的区别： 网络均匀时，内部控制容积有半格错位，但边界节点代表的控制容积不同。 1. 网格不均匀 外节点法的节点不在控制容积的中心 内节点法的节点总是在控制容积的中心 2. 网格不均匀 外：界面永远在相邻节点的中间 界面不在相邻点的中间，导致界面上导数计算截差下降 3. 内节法在计算求解域物性发生阶跃变化时，可以该界面 做为界面，而使控制容积内的物性一致。 4. 在讨论了离散的网格后，再研究如何把控制方程离散。 一类从纯数学方法推导；另一类从物理概念推导，各有特点。 首先介绍一种通用的一维非稳态对流—扩散方程。（直角坐标形式） 非守恒型： 守恒型： 非稳态项 对流项 扩散项 源项 Φ：广义变量（速度、温度、浓度┅┅） Γ：广义扩散函数（粘性函数、导热函数┅┅） S：源项 线性化处理：u.物性为已知的常数。 拟线化处理：每次计算中为常数，两次迭代间进行修正。 一、台劳线数展开法（有限差分法） 　把函数中　　　 在（i+1,n）点 对（i，n）点展开。 1 △x 1 1 i–1 i i+1 向前差分 截断误差(一阶精度） 向后差分 中心差分 二阶精度 向前 向后 中心 向前 二、多项式拟合法 主要用于处理边界条件，获得壁面节点的离散方程。 例：已知壁面热流 壁面附近温度按线性分布 x 1. 时间向前差，空间中心差，显示格式.（FTCS） 隐式格式：主要计算 的值。 *