计算流体力学第2讲差分方法.pptxVIP

9/295/10第二讲有限差分法（2）李明军湘潭大学数学与计算科学学院数学楼315；计算流体力学讲义第3节差分格式的进一步分析第4节从模型方程推广到NS方程（Euler方程）

第3节差分格式的进一步分析9/29/2022半离散分析与全离散分析----如何理论计算修正波数02耗散与色散误差01

1.耗散与色散误差9/295/103CopyrightbyLiMingjun数值实验时间推进：3步TVD型Runge-Kutta,且时间步长足够小（误差忽略）空间离散：1阶及2阶迎风格式（20个网格点）实验观察到的现象——两类误差：振幅误差相位误差（波速误差）

精确解1阶迎风2阶迎风

对以上“实验现象”进行理论分析考查问题：精确解：差分格式：（1）其他格式……假设对于：有隐含假设:线性差分格式，非线性系统作用于单波，会产生多个谐波（2）差分没有误差1阶迎风格式：2阶迎风格式：

精确解：2.半离散分析与全离散分析----如何理论计算修正波数离散带来的误差（难度小、常用）半离散分析：假设时间推进是精确的，仅分析空间注意：实际上就是普通三角函数，采用复数形式仅仅是为了理论推导方便。若用实数形式sin(kx),cos(kx)推导形式上略显繁琐。全离散分析：同时分析时、空离散的误差（难度大）

如果，无误差！4设：5令：1式化为：2“半离散化”：空间导数差分计算，时间方程（常微）精确计算3精确解：6

分析（修正波数）与误差的关系：0504020301的误差导致解的幅值误差——耗散误差的误差导致解传播速度的误差——色散误差(理想情况：)设则的实部：耗散误差的虚部：色散误差

定义：设，含义：反应波数（谱）空间内差分的误差01任意函数：02求导数，精确解03差分解04修正波数定义为：05

Fourier分析的任务：9/295/1011CopyrightbyLiMingjun反映了一个波内的点数。PPW（波内的点数）=计算出，并考差其与的逼近程度。考察格式分辨率（resolution）的重要指标优秀的差分格式，1个波长里面6个点即可精度分辨率精度：反映时的情况分辨率：网格点数很少（例如波里面只有6个点）时的性能，对于多尺度问题，分辨率更重要。（牺牲精度，提高分辨率）

如何理论计算修正波数？01方法：根据差分具体表达式及定义计算02例1：03令04则：05其中。于是：061阶迎风格式：07

2阶迎风格式：9/295/1013CopyrightbyLiMingjun1例2：2令3则：4其中。

第4节从模型方程推广到NS方程

（Euler方程）014.1Jacobian系数矩阵及其性质02严格特征分裂——基架点上冻结系数034.2对流项的分裂两类典型的逐点分裂04

4.1Jacobian系数矩阵及其性质格式F+1格式F-2对流项：信息（波）从上游传至下游，上游更重要，迎风差分3扩散项：信息从中心向周围扩散，不区分上、下游，中心差分4迎风差分优点：有效利用信息传播的方向，增强稳定性5S方程：6单波方程：7单波方程：一个波，容易判断波传播方向8N-S对流项（Euler）：方程组，多波问题，复杂9

何为迎风格式？9/295/1016CopyrightbyLiMingjun单波方程：对于正数a,说明波沿x正方向传播，如下图所示：采用空间向后差分格式：

两类典型系数矩阵A：9/295/1017CopyrightbyLiMingjun常系数矩阵A变系数矩阵A的情况：局部讨论常系数矩阵A的情况：完全解耦，独立求解双曲方程组的原则：特征