四高精度差分格式及其数值解的逼近程度分析.ppt

(四) 高精度差分格式及其数值解的逼近程度分析 指大于二阶精度的格式 要求准确模拟小扰动量长时间、远距离传播的速度和幅值 用于计算噪声、DNS等 有传统高精度差分和紧致差分两种 4.1模型方程及半离散方程 模型方程及半离散方程(cont.) 半离散方程的精确解 4.2 高精度差分格式 特别适用于光滑问题 减少网格点数 应用于湍流等多尺度问题 4.2.1 传统型差分格式 传统型差分格式(cont) 传统型差分格式(cont) 传统型差分对模型方程的逼近特性 所有中心型差分：无耗散 迎风偏斜差分：有耗散，但可能为负耗散 可从差分格式的精确解分析色散和耗散 传统型差分格式(cont.) 二阶导数的差分逼近 4.2.2 紧致差分格式 同样的精度比传统差分的基架点少 截断误差的系数较小 二阶导数的(中心型)紧致差分 一阶导数的对称紧致差分 一阶导数的迎风紧致差分 二阶导数的紧致差分 传统型差分的截断误差项的再次离散 四阶紧致： 若边界点S0和SN已知，可用解三对角矩阵方程得到所有网格点上的差分 二阶导数的紧致差分(cont.) 一般形式 获得S需要求解矩阵方程组 二阶导数的紧致差分(cont.) 5点六阶紧致： 一阶导数的对称型紧致差分 一般形式: 获得F需要求解矩阵方程组，左端点最好不超过3 一阶导数的对称型紧致差分(cont.) 四阶精度: 也可从传统的二阶中心差分的截断误差的再次离散 一阶导数的对称型紧致差分(cont.) 六阶精度: 一阶导数的迎风紧致差分 有数值耗散，抑制高频振荡 一般形式： 三阶迎风紧致格式(3点）： 一阶导数的迎风紧致差分(cont.) 五阶迎风紧致格式(5点）： F+按j增加的方向求解, F-按j减少的方向求解 非等距网格上的格式 计算精度比 要高。 三阶迎风紧致： 一阶导数的迎风紧致差分(cont.) 三阶迎风紧致格式的截断误差： 对模型方程的逼近程度： 一阶导数的迎风紧致差分(cont.) 三阶迎风紧致格式对模型方程的逼近程度： 4.3 格式的精度及分辨率 耗散误差和色散误差的行为 耗散：影响振幅 色散：产生振荡 物理上的高波数(K很大） 它使首项误差，甚至更高阶项误差被放大： 模拟高波数问题对差分格式的要求 4.4 数值解中的耗散效应和色散效应 所有中心型格式都是无耗散的 负耗散型格式数值上是不稳定的 强迎风格式（点全在上风区）可能是负耗散 耗散效应 色散效应 群速度 正弦波： 群速度 (cont) 两个正弦波叠加： 叠加后看作一个简单的正弦波，其振幅为 该振幅本身也是一个波（虚线所示），其相速度为 数值解的群速度 4.6 时间离散的耗散和色散效应 第四次课阅读提示 傅《计算流体力学》第四章全部。 第四次课后作业 傅德薰《书》习题 4-2 傅德薰《书》习题 4-4 傅德薰《书》习题 4-5. 一阶隐式：正耗散 二阶：色散和耗散，色散占主导 。 * K=1 K=2 直接由网格点的函数值线性组合 N阶精度至少需要N+1个网格点 Taylor展开确定系数 J J-2 J-3 J-1 J+1 J+2 J+1 J-1 J-2 J J+2 J+3 *