第八章 流场计算的新方法 90年代后期相继发展了几类新的差....ppt

小结 新的格式 TVD格式 ENO格式 NND格式 新方法 谱方法 自适应网格 多重网格方法 并行计算与向量计算 非结构网格及其应用 无网格方法 格子波二茨曼方法 2、利用中心型差分格式构造显式二阶TVD格式 对称TVD格式：用中心差分构造的TVD 利用L－W格式构造二阶精度TVD格式 设有标量方程 C＝常数0 用L－W格式 改写后 此格式 第一部分为一阶迎风格式 第二部分为反扩散修正项 引入限量因子 控制反扩散量，格式满足TVD特性 选择 的方法 且非负 几种形式的取法 （1） （2） （3） （Roe-Sweby） 范利尔（Van Leer） 查克拉瓦特-奥雪尔（Chakravarthy-Osher） 新方法 湍应力模型 LES DES DNS 第八章 流场计算的新方法 90年代后期相继发展了几类新的差分格式和技术 TVD格式： 总变量衰减格式（Total variation Diminishing） ENO格式： 基本无振震格式（Essentially Non Oscillatory ） NND格式： 无振震、无自由参数、无耗散格式（Non oscillatory containing no free parameters and dissipative scheme） 谱方法： 快速付立叶变换方法 自适应网格技术（Adaptive Meshes） 多重网格方法 并行计算与向量计算 非结构网格及其应用 无网格方法 §8-1 TVD方法 一、 模型精度与网格技术 流动区域：光滑区、边界层、剪切流动区、分离流区、 激波和含有间断面的区域。 流动结构与尺寸网格尺度与流体运动结构尺度相对应 光滑区：参数变化光滑 易于构造差分格式 重点是提高求解效率 对于粘性影响区与无粘区可根据流场特点简化 简化的N-S（抛物化N-S方程） 分离流涡运动 无粘特征和粘性特征都很明显（大涡控制） 2.激波附近的数值解 激波装配法 ：把间断面分开成内外边界，边界采 用激波关系式 不适用于未知激波（悬挂激波） 激波捕捉法 ：采用统一方法进行计算 一阶格式：激波附近无振震，但会抹平激波 二阶格式： 激波附近会产生非物理振震 高阶格式： 需要采用人工粘性抑制 TVD格式：具有较高的激波分辨率 ENO格式：问题是：①.是否具有二阶精度。②.是否满足熵 条件。 3. 粘性流计算中的网格与逼近精度 误差来源：真实粘性被差分截断误差掩盖 例：伯格斯方程（一维粘性流方程） 及 如果坐标变换满足 上式没有一阶空间导数。若 则坐标变换方程 的解为： 上式表明在计算平面内网格点是均匀分布的，物理平 面接近x=0处是加密的 0 x 差分算子 原方程的差分格式 定义网格Re 网格雷诺数在差分方程中起着权衡作用，代表着惯性力 与粘性力各自起作用的大小不同。 单调性(u(x))要求:各系数为非负数 （差分方程粘性稳定性条件） 在粘性计算时的要求 网格雷诺数与流场性质相关 粘性起主导作用的区域网格Re较小 无粘占主导作用的区域Re很大，可以不考虑网格Re影响 二、一阶逼近精度的TVD格式 激波捕捉法的困难 ①流场参数梯度较大区，二阶或高阶精度格式有非物理震荡产生 ②可能出现非线性不稳定 ③可能出现非物理解 已知的结论： ①一阶逼近可得到满足熵不等式的惟一物理解（弱解） ②二阶或高阶格式求解具有间断面的物理问题时需要加人工粘 性，克服伪装解，但引起信息丢失 ③新发展差分格式具有两个特点（TVD） 对非线性方程及常系数双曲性方程总变差衰减守恒律 与熵不等式相容 TVD具有一阶精度和上述新特点 ENO具有二阶精度和上述新特点 1.TVD格式的物理考虑 总变差：一元函数U(X)的定义域被划分成间隔 的多段 TV(u)为 的上界 x O j-1 J j+1 对离散函数函数 （k= ,…, ）