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

第八章 流场计算的新方法 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处是加密的;差分算子;原方程的差分格式;（差分方程粘性稳定性条件）;二、一阶逼近精度的TVD格式 激波捕捉法的困难 ①流场参数梯度较大区，二阶或高阶精度格式有非物理震荡产生 ②可能出现非线性不稳定 ③可能出现非物理解;对非线性方程及常系数双曲性方程总变差衰减守恒律 与熵不等式相容 TVD具有一阶精度和上述新特点 ENO具有二阶精度和上述新特点;1.TVD格式的物理考虑;对离散函数函数 （k= ,…, ） 2.显示TVD格式 对于简单的方程 三点守恒差分格式为;三种典型的数值通量形式 (1) L-W格式;(2) L-F格式（Lax-Friedrichs）; ------TVD格式的充分条件 TVD格式不满足熵不等式，因此可能收敛于非物理解 改进的通量函数（与熵不等式相容） ;3.隐式TVD格式 具有单参数的五点守恒差分格式 设 则; 其中 通量相容条件 当 时，得到显示格式 当 时，得到梯形格式 当 时，得到欧拉时间后差格式 ;用差分算子表示这一差分格式 其中算子 此格式为TVD的充分条件是 ;4、半离散形式的TVD格式 只在空间离散，则离散方程为 若 ，则方程具有TVD特征。 ;总变差; ;三、高阶TVD格式构造 用一阶TVD格式构造二阶TVD格式 用L-W格式构造二阶TVD格式 利用限量因子（limiter)限制函数 或通量函数 的梯度变化 迎风或对称TVD格式利用迎风格式或 中心格式采取限量因子构造的TVD格式;对于修正后的通量列出差分格式;其中，;;2、利用中心型差分格式构造显式二阶TVD格式;用L－W格式 ;引入限量因子 控制反扩散量，格式满足TVD特性 ;几种形式的取法 （1） （2） （3） ;小结;9、春去春又回，新桃换旧符。在那桃花盛开的地方，在这醉人芬芳的季节，愿你生活像春天一