计算流体力学课程教学课件第六讲差分方法.pptVIP

Copyright by Li Xinliang Copyright by Li Xinliang 计算流体力学讲义2011 第六讲 差分方法（4） 知识点： 通量技术简介——Roe 常用的隐式处理方法——LU-SGS * * 知识回顾： 单调、保单调和TVD 概念： 网格Reynolds数 单调格式、保单调格式及TVD格式 Harten定理： 正系数原则 TVD 保单调 单调 TVD格式= 1阶迎风+ j *（修正项） 二阶精度区 TVD区 二阶精度TVD区（二者交集） * 知识回顾：WENO格式 基本思路 {j-3, j-2,j-1,j,j+1,j+2} {j-3,j-2,j-1,j}； {j-2,j-1,j,j+1}； {j-1,j,j+1,j+2} 五个基架点被分成三个组 1） 若高精度逼近 ， 必然利用多个基架点 2） 如果该基架点内函数有间断，会导致振荡 3） 间断不可能处处存在 4） 把基架点分成多个组（模板）， 每个模板独立计算j点导数的逼近。 —— 得到多个差分 5）根据每个模板的光滑程度，设定权重 6） 对多个差分结果进行加权平均 。光滑度越高，权重越大。如果某模板存在间断，则权重趋于0； 如果都光滑，则组合成更高阶格式。 * 1. 单方程的Roe格式 线性化，用平均变化率代替(j,j+1)之间的变化率a(u) “平均斜率”，不等于“斜率的平均值”，也不等于中点处的斜率 § 6.1 Roe格式 非线性情况 根据Langrage中值定理，[uL,uR]之间必有一点uRoe, 该点处的斜率为平均斜率; 二次函数f(u)=u2中点处的斜率=平均斜率 * 2. 方程组的情况 平均斜率 线性化，以平均增长率代替瞬时增长率 [j,j+1]区间内 连续，且 可通过相似变换对角化 应当具有的性质 常系数方程的Riemann解 * 平均斜率 x j+1/2 常系数单波方程的Riemann解 Roe 格式： 微分型近似Riemann解 * 3. 矩阵 的构造 关键： “向量除以向量” ？ 直接求平均增长率： u f(u) uL uR uRoe Roe点的斜率为平均斜率 （根据拉格朗日中值定理，[UL,UR]区间内肯定存在Roe点） 思路1： 在UL与UR之间寻找一个点URoe, 该点处的增长率为平均增长率 f(u)=u2 u 二次函数—— Roe点与中点重合 标量函数的启示： Roe点肯定存在（Langrage 中值定理） 二次函数的中点即为Roe点 思路2： 进行坐标变换，得到一个二次（齐）函数 引入 如果 是二次（齐）函数，则其中点 即为Roe点 重要启示 更准确地讲，应当是要求 为W的线性函数， 即增长率为线性函数 （中点处的增长率刚好为平均增长率） * 针对Euler方程的具体构造 引入新变量： 则： 目的： 使得F(w)是W二次齐函数 （增长率为线性函数） f(U)不是U的二次齐函数 二次齐函数！ 中点处的斜率即为平均斜率。 Roe点 Roe点为： 增长率为线性函数！ * 最终： 其中 如下计算： 平均增长率（矩阵） 含义： 左、右两个状态点的某种平均 （称为Roe平均，为密度加权平均） 该状态点对应的增长率（矩阵）为平均增长率（矩阵） 实际上是一种“等效平均”。 效果优于简单的算数（或几何）平均。 三维情况下，还有 其他量（如压力、温度、音速等）用这三个量计算 （5） 简单易记： * Roe 格式的计算步骤 （半离散） 已知n时刻所有网格点上的物理量，对于j点： 1） 利用差分格式计算UR,UL 2） 采用Roe平均公式（5）计算Roe平均值 3） 将Jacobian矩阵 进行特征分解: 计算 4) 计算 5）计算 6） 计算空间导数 7）时间推进，计算下一时间步的值。 j-1 j j+1 与前文（第3，4讲）的形式相同，仅需把式中的密度、压力、速度等换成经过Roe平均的密度、压力、速度即可 其中： * 可能出现导数不连续， 可能引起数值振荡 实际使用时 可用如下函数代替 ——所谓“熵修正” 实际上是在特征值0点周围增加了耗散 Roe 格式的优点： 1） 保持守恒性的同时，严格保证了特征方向 2） 便于