通过代数方程计算出网格点的位置优点.ppt

Copyright by Li Xinliang Copyright by Li Xinliang 计算流体力学讲义 第十一讲 代数方程的求解及网格生成 李新亮 lixl@ ；力学所主楼219； 知识点： * 讲义、课件上传至 (流体中文网） - “流体论坛” -“ CFD基础理论 ” 讲课录像及讲义上传至网盘 /browse.aspx/.Public Copyright by Li Xinliang 代数方程组的求解 网格生成 Copyright by Li Xinliang * 知识回顾—— 有限体积法 在以某节点为中心的控制体上积分 i,j k 非结构网格的控制体 i+1,j i-1,j i,j+1 i,j-1 k3 k1 k2 k4 k5 结构网格的控制体 x y n 体积平均 控制体边界垂直于节点连线（也可选其他方式） 垂直平分线 n 1) 建立控制体 2) 在控制体上积分 离散方程 重构： 由节点上平均值 给出函数分布，最终给出通量 表示第m个界面上的值 1. 有限体积法的离散过程——重构 1） 重构—— 两种不同的重构方案，向左偏及向右偏。 给出两种结果： 及 Copyright by Li Xinliang * 例如： 0阶重构： 线性重构: 2. 迎风型有限体积法 注： 有激波的情况，需使用限制器 （即激波捕捉格式） j-1 j j+1/2 如不使用限制器，很可能“寸步难行” （第1步就发散） 例： Sod 问题 压力出现负值，计算发散 例如： NND 格式 在 二阶中心 / 二阶迎风/1阶迎风 三个方案中选取最合适的。 1阶迎风 修正项 线性插值， 得到二阶迎风格式 两个修正方案（二阶迎风及二阶中心），给出趋势相同时，选修正量最小的，趋势相反时，不修正 Copyright by Li Xinliang * i,j i+1,j i-1,j i,j+1 i,j-1 n 左重构 右重构 2） 由左右重构得到的自变量： 和 给出通量 方案A: FVS 方案B: 解Riemann问题 （常用） x y 看似二维Riemann问题，其实是一维的，坐标旋转一下就行了 精确Riemann解 Roe-近似Riemann解 （求解常系数线性化的Euler方程） HLL型近似解—— 双激波近似，积分平均解 HLLC型近似解—— 三波近似 3. 粘性通量的计算—— 中心型 Copyright by Li Xinliang * Part 1 代数方程组的求解 微分方程（组） 代数方程组 数值解（离散解） 差分 有限体积 大部分计算量 11.1.1 Gauss 消去法 消元 为了计算稳定，通常使用主元消去法 列主元消去法； 全主元消去法 计算量： 乘法： 加法： 11.1 代数方程组求解的直接法 优点： 简单精确，缺点：计算量大 上三角 矩阵 Copyright by Li Xinliang * 11.1.2 LU分解法 …… 0 0 Step 1 Step 2 Step k 对角线上不能有0， 计算之前先交换矩阵A的元素，将主值交换到对角线上 Copyright by Li Xinliang * 回代过程 计算量： 分解O(n3/3), 回代 O(n2) 优点： 1） 重复求解 ， …… 时，仅需一次LU分解，计算量小； 2） LU分解不破坏带状稀疏矩阵的性质，可大幅减小计算量。 L带宽的带状矩阵： LU 分解: O(nL) 回代： O(nL) Copyright by Li Xinliang * 11.1.3 带状矩阵求解的追赶法 追赶法：等价于带状矩阵的LU分解 例： 三对角矩阵 一般项： 边界项： 追赶法 令： Step 1: Step 2: Step 3: Step 4: 计算量： 9n次 （乘法） A为固定值时， 3n次（乘法） 简单易用，计算量小 Copyright by Li Xinliang * 11.2 代数方程组求解的迭代法 直接法 迭代法 优点 算法简便，准确（未知数少时） 计算量小，误差容易控制 缺点 计算量大O(n3) 舍入误差积累，不易控制 快速收敛的算法设计较为复杂 11.2.1 Jocabi 及Gauss-Seidel迭代 解出对角元素 Jocabi迭代 Gauss-Seidel迭代 “对角占优” Copyri