18.高级算法初步.ppt

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正负， 其中 其迎风格式恒等价于 其中 上式可写为 其中 设 对应的本征矢量为列向量 ，即 定义矩阵 ，其逆矩阵由 组成 左特征矢量 A的Jordan标准型 以1D欧拉方程组为例 对欧拉方程，有 F=AU （代入即可证） 又 定义 和 分别为由A的正负本征值构成（暂对亚音速流，令u0） 形式上A如同常数矩阵 则 F =F++F－，其中 此即flux-vector splitting一例。另一种分裂方式为 对应沿 x 正向的通量，故应采用后向差分 对应沿 x 负向的通量，故应采用前向差分 有些改进的格式，如只适合Euler方程的Van Leer’s flux splitting: 设局地马赫 数Ma=u/c 则 其中 典型通量分裂法 对形如Euler方程的模型方程 其解的一个重要性质是 差分方程中 若 则此格式称为total-variation-diminishing (TVD) 可见一阶迎风格式满足TVD条件，而二阶迎风及中心差分格式则易导致振荡，从而不满足此条件。采用二阶格式且让它满足TVD条件的方法之一是用一非线性函数乘上差分方程中的相关元素，强制该差分方程满足TVD条件，称为Flux limiters。 Harten给出判据，若格式可写为 且对所有i均有 则该格式为TVD格式 其中R(U)和L(U)分别是由A的右、左特征矢量组成的矩阵。 可构造一阶TVD格式： Q的几种形式 对应后面讲的Roe格式 1959年Godunov提出一种与以前差分方法截然不同的方法，他建议将Euler方程的局部精确解组合起来得到全域流场。 什么是局地流场的精确解？ Riemann (1858)首次解出，故称Riemann问题，它即Euler方程的局部精确解。 求满足以上初条且在间断线上满足间断条件的解 Shock tube Riemann 间断解（亦称wave diagram） 真空 为接触间断（discontinuity） 数值解总是分段常数 可算出新值后求平均作为 亦可求通量。 这样做为一阶精度，故耗散强。二阶精度为 但每格点处都需求解黎曼问题很复杂，因此不断有人提出近似solver。 1961年Godunov把所有的间断都用声波代替（定义M为单位时间单位面积上扫过的质量流量）： 只适用于跨格点的物理量的变化仅为网格间距的数量级。 Roe格式与Roe分解 亦为一种迎风格式。对单个守恒型方程 由此构造Roe格式 其中 可推广迎风格式（由于针对守恒型方程，故适合解含间断解的问题） 推广到双曲型方程组 这样 关键是要找到合适的A(U, V)。 其中A为非线性，故复杂。Roe试图通过左右函数的常数态UL和UR构造一个合理的常矩阵A(UL,UR)，即可将迎风格式推广。 可以证明，这时上式的黎曼问题的通量可写为 矩阵A(U, V)满足以下U特性(Uniform) A对U, V是Lipschitz连续，且A(U, U)=A(U)= A(U, V)有N个实特征值及特征矢量的完备集，记其特征值、左特征行矢量和右特征列矢量分别为?k(U, V), lk(U, V), rk(U, V), 使得由其组成的特征矢量矩阵L(U, V)和R(U, V)互