5 计算流体力学基础一.ppt

计算网格单元与节点存储 代数方程组 合适的节点排序 稀疏矩阵 带矩阵 离散误差 微分方程真实解 截断误差 离散方程真实解 离散误差 线性问题： 离散误差估计与 Richardson 外插 ? 网格足够细 ? 单调收敛 Richardson 外插 : R d h h h f f e = + ( ) , r O h r p > : 有限体积法 有限体积法基本思想 nodes centered in CVs CV faces centered between nodes 保证守恒性 数值积分 代数方程 Well- conditioned 面积分计算 面积分计算 面积分计算 ? 两步近似 ? 数值积分 ? 被积函数近似：插值 数值积分方法： 中值公式： 梯形公式： Simpson 公式： 二阶精度 二阶精度 四阶精度 数值积分误差估计 中值公式： ? ? ? ? 2 a b M f b a f ? ? ? ? ? ? ? ? ? 梯形公式： Simpson 公式： ? ? ? ? ? ? ? ? 2 f a f b T f b a ? ? ? ? ? ? ? ? ? 4 6 2 b a a b S f f a f f b ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 体积分计算 二阶精度 四阶精度 九点插值 均匀网格 Fluent 中的积分计算 二阶精度的中值公式 面函数近似：插值（ 离散格式 ） ? 迎风格式（ UDS ） ? 二阶迎风 ? 线性插值（ CDS ，中心差分格式） ? QUICK 格式 ? 乘方格式 ? 三阶 MUSCL 格式 迎风格式（ UDS ） ? 唯一的无振荡格式 ? 数值耗散大 ： 一阶精度 数值（虚假）扩散 线性插值（ CDS ） ? 最简单的二阶精度格式 ? 等价于有限差分方法中一阶导数的中心差分 (CDS) ? 振荡解 二阶精度 解决方法： 1 、推迟更新 2 、 Gradient limiters: 面上的解不超过相邻单元值。 Fluent 中的 CDS ? pressure-based solver only ? LES 抑制振荡 ? Bounded Central Differencing Scheme ? 混合 CDS 和 UDS ? 消除振荡， LES 的默认格式 线性插值（ CDS ）：面导数近似 ? e 在 P 、 E 中点时为二阶精度 ? 网格加密时收敛速度为二阶 一阶精度 Fluen 中网格面上的导数（扩散项）近似 都采用 CDS ，二阶精度 QUICK 格式 ? Quadratic Upwind Interpolation for Convective Kinematics) ? 与中值积分联用，总精度为二阶，比 CDS 稍好 三阶精度 Fluent 中二阶迎风与 QUICK 格式 ? 二阶迎风 ( 上游两点 线性插值 ) ? QUICK 格式 ? CDS 与二阶迎风的加权平均； ? 四边形（ 2D ）、六面体（ 3D ） ? θ 随流场可变 ? pressure-based solver 高阶格式 ? 高阶离散格式只有与积 分精度相匹配才有意义 ? 只有在足够细的网格上 才有相应的高精度 系统加密 Fluent 中的乘方 (Power-law) 格式 ? |Pe|>1, ? = ? u ? |Pe|=0, ? = a ? u+ b ? d ? 0<|Pe|<1, 精确解的乘 方近似 精确解 Fluent 中的三阶 MUSCL 格式 ? Monotone Upstream-Centered Schemes for Conservation Laws ? CDS 与二阶 UDS 的混合 ? 可用于任意网格 ? 对解没有限制，在间断面附近会产生振荡 计算流体力学 CFD 基本要素 ? 数学模型 ? 离散方法 ? 坐标系 ? 计算网格 ? 有限近似 ? 求解方法 ? 收敛准则 CFD 方法性质 ? 方程的相容性 (Consistency) ? 数值算法的稳定性 ? 解的收敛性 ? Lax 等价定理 ? 守恒性