18.高级算法初步.pptVIP

18. 高级算法初步 1. Flux-Vector Splitting 2. TVD及Flux limiters 3. Godunov Approach * * 以前算法简单实用。现代新算法得益于应用数学的帮助，效率较高。由于多应用激波捕捉方法，故常用守恒型方程，即使流场不含激波。 其准线性形式为 A, B, C为雅可比(Jacobian)行列式，其本征值决定了方程的数学性质：若为不同实数，则为双曲型，若相同实数则抛物型，若均为虚数则椭圆型。 1D流体力学方程组（Euler方程组） 即 其本征值： 故Euler方程组为双曲型。 很多新算法的差分格式都与特征线方向有关。 简单的例子：迎风格式 第 i -1 点影响第 i 点，若采用中心差分则不稳，而迎风格式 中心差分常包含影响域之外的信息。在平滑流场处这不是问题，但出现激波后，中心差分并不好，会引起振荡等，从而导致迎风格式(upwind scheme)的发展，以真实反映特征线。若处理得好将无振荡。又如单波模型方程 则稳定，缺点是仅一阶精度，且强耗散。 其改进方法如TVD格式、flux splitting、flux limiters、Godunov格式及近似Riemann solvers（它们可统称为迎风格式，因为它们都尽量正确处理信息的传播）。 单波模型方程 无论c