5.5对流项的高阶迎风格式讲解.ppt

源项线化 要求线化后的斜率Sp为负值 Sp为负值，才能保证最后求解方程组的系数矩阵对角占优，有利于收敛 界面当量导热系数 加权平均方法 调和平均方法 (最符合物理意义) 网格划分和边界条件 网格划分： 点中心 块中心 三类边界条件： 1 函数值(Dirichilet) 2 导数值(Newman) 3 导数值与函数值的关系(Robin) 方程的求解 方程非线性性质的处理(外迭代) 线化代数方程组的解法： TDMA，PDMA 迭代方法(内迭代) 建立控制方程、确定初始与边界条件 流 程 图 解域离散、方程离散 初边条件离散 给出节点初始温度值 计算系数，固定， 线化代数方程 求解离散的 线化的代数方程组 进入下一时层求解 解收敛否？ Yes No 以新温度值 替代老温度值 线化方程的 迭代解法 (内迭代) 外迭代 特别注意 虽然差分离散格式方程为全隐式，但其迭代求解方法可以是显式迭代(点迭代)，隐式迭代(块迭代/线迭代)，或者交替方向隐式迭代(交替方向线迭代)。 差分离散方程为时间层演进，上标为n, n+1 迭代格式代表迭代层演进，上标(n),(n+1) 5.5 对流项的高阶迎风型格式 胡 茂 彬 /~humaobin/ humaobin@ 目标和主要内容 减轻虚假扩散：人工粘性引起的流向扩散、网格取向效应引起的交叉扩散、非常数源项引起的虚假扩散 增加数值计算的准确性 1.二阶迎风格式 2.三阶迎风格式 3. QUICK格式 5.5.1 二阶迎风格式 基于Taylor展开的待定系数方法 与一阶迎风格式的区别 改写为 一阶迎风 修正 一阶迎风估计斜率为割线， 实际上斜率是切线 与一阶迎风格式的区别 改写为 一阶迎风 修正 控制容积积分法离散 对流项的积分 与有限差分离散结果相同，都是二阶精度 5.5.2 三阶迎风型格式 基于泰勒展开待定系数方法 控制容积积分法离散 界面上的函数值 评 述 在离散节点的流动下游取了一个节点，迎风性差 离散精度提高到三阶 格式变为条件稳定 具有守恒性 5.5.3 QUICK格式 对流运动的二次迎风插值(Quadratic Upwind Interpolation of Convective Kinematics) 割线中点值 修正 曲率修正 e界面上 w界面上 界面值修正采用了迎风思想 5.5.4 对流项采用高阶格式时引出的新问题 涉及的节点数目增多 一维问题：P、W、E、WW、EE 二维问题：P、W、E、S、N、WW、EE、SS、NN 问题！ 1. 紧邻边界的第一个内节点的离散方程如何构造？ 2. 离散方程如何求解？ 紧邻边界的第一个内节点的离散方程如何构造？ 1. 开拓一个虚拟节点0，插值确定 2.采用一阶迎风或 混合格式处理边界条件 离散方程如何求解？ 一维情况：可采用五对角阵方程的直接求解方法求解(PDMA) 二维情况：可采用交替方向的五对角阵算法，还可采用延迟修正方法 5.6 对流项离散格式的稳定性问题 三类稳定性问题 1.代数方程迭代求解过程的不稳定性：迭代方法选择不当，迭代不收敛，导致迭代发散的情形 2.初值问题离散格式的不稳定性：由于时间(类时间)步长取得过大或者时间(类时间)与空间步长匹配不当，导致解振荡发散 3.对流格式的不稳定性：在采用某些离散格式求解对流扩散方程时，由于空间步长过大或者流速过高而网格Peclet数过大，也会使解振荡发散 对流项格式的不稳定性 分析方法： 通常将离散格式用于最简单的一维稳态无源项的模型方程来讨论。如果最简单的问题都不能正确解决，那么。。。 正系数法，离散方程精确解分析法，反馈灵敏度分析法，符号不变法 对流项格式的不稳定性 具有迁移性的对流离散格式无条件稳定： 相对离散出发点的单边型一阶、二阶迎风格式具有迁移性，因此恒稳定。 上下游都有节点的离散格式 中心型或偏心型的格式都不具有迁移性，这些格式只能有条件稳定；且格式在下游节点的系数越小，相对稳定性越强，临界网格 Peclet 数越大 临界网格Peclet数：在此Peclet数之下，离散格式是稳定的 上下游都有节点的离散格式 几种常见差分格式的临界Peclet数： 中心差分格式： Quick格式： 三阶迎风格式： 实际计算发生数值解振荡的临界Peclet数要比简化模型分析得到的结果要大 简化模型的五个条件：一维、线性(物性均为常数)、无源、两点边值问题、均匀网格 解除5个苛刻条件中的任何一个，使解发生振荡的 都要增大 实际使用中，要通过数值实验，逐步扩宽临界Peclet数的选择范