计算流体力学第2讲差分方法2初步.pptVIP

第3节 差分格式的进一步分析 第4节 从模型方程推广到NS方程（Euler方程） * * 第二讲 有限差分法（2） 李明军 湘潭大学 数学与计算科学学院 数学楼315； email: alimingjun@163.com 计算流体力学讲义 第3节 差分格式的进一步分析 第4节 从模型方程推广到NS方程（Euler方程） 1. 耗散与色散误差 2.半离散分析与全离散分析 -----如何理论计算修正波数 1. 耗散与色散误差 * Copyright by Li Mingjun 数值实验 时间推进： 3步TVD型Runge-Kutta, 且时间步长足够小 （误差忽略） 空间离散： 1阶及2阶迎风格式 （20个网格点） 实验观察到的现象—— 两类误差： 振幅误差 相位误差 （波速误差） * Copyright by Li Mingjun 精确解 1阶迎风 2阶迎风 * Copyright by Li Mingjun * 对以上“实验现象”进行理论分析 考查问题： 精确解： 差分格式： （1） 其他格式 …… 假设对于： 有 隐含假设: 线性差分格式， 非线性 系统作用于单波， 会产 生多个谐波 （2） 差分没有误差 Copyright by Li Mingjun 1阶迎风格式： 2阶迎风格式： * 半离散分析： 假设时间推进是精确的，仅分析空间 离散带来的误差（难度小、常用） 全离散分析： 同时分析时、空离散的误差 （难度大） 注意：实际上就是普通三角函数，采用复数形式仅仅是为了理论推导方便。若用实数形式 sin(kx), cos(kx)推导形式上略显繁琐。 Copyright by Li Mingjun 精确解： 2.半离散分析与全离散分析 -----如何理论计算修正波数 * 令： （1）式化为： “半离散化”： 空间导数差分 计算，时间方程 （常微）精确计算 如果 ， 无误差！ Copyright by Li Mingjun 设： 精确解： * 分析 （修正波数）与误差的关系： 的误差导致解的幅值误差 ——耗散误差 的误差导致解传播速度的误差 —— 色散误差 Copyright by Li Mingjun (理想情况： ) 设 则 的实部：耗散误差 的虚部：色散误差 * 含义： 反应波数（谱）空间内差分的误差 任意函数： 定义：设 ， 求导数，精确解 差分解 Copyright by Li Mingjun 修正波数 定义 为： 反映了一个波内的点数。 PPW （波内的点数）= * 计算出 ，并考差其与 的逼近程度。 考察格式分辨率（resolution）的重要指标 优秀的差分格式，1个波长里面6个点 即可 Copyright by Li Mingjun 精度 分辨率 Fourier 分析的任务： 精度： 反映 时的情况 分辨率：网格点数很少（例如波里面只有6个点）时的 性能， 对于多尺度问题，分辨率更重要。 （牺牲精度， 提高分辨率） * 如何理论计算修正波数？ 方法：根据差分具体表达式及定义计算 例1： 令 则： 其中 。于是： 1阶迎风格式： Copyright by Li Mingjun * 例2： 2阶迎风格式： 令 则： Copyright by Li Mingjun 其中 。 4.1 Jacobian 系数矩阵及其性质 4.2 对流项的分裂 (2) 严格特征分裂—— 基架点上冻结系数 (1) 两类典型的逐点分裂 格式F+