计算流体力学课件：有限体积法1.pptxVIP

有限体积法（1）知识点： 1 有限体积法的基本概念 无粘通量及粘性通量的计算 多块网格 2知识回顾： Roe 格式平均斜率线性化，以平均增长率代替瞬时增长率[j,j+1]区间内连续，且可通过相似变换对角化 应当具有的性质常系数方程的Riemann解 3知识回顾2： LU-SGSStep 1: 求解Step 2: 求解 § 7.1 结构网格有限体积法有限体积法主要优势： 处理复杂网格差分法处理复杂外形 —— 坐标变换无粘项需要计算9次导数，计算量大；守恒性被破坏（1）两端同乘无粘项仅需计算三次导数 5差分法有限体积法优点简单、计算量小、易于提高精度本身包含几何信息，易处理复杂网格不足差分离散与几何解耦，难以处理复杂网格复杂、不易提高精度均匀流场数值解（误差），未必为0几何诱导误差 ：对于均匀场实际问题： 外形复杂， 光滑的结构网格生成困难不易保障几何守恒性坐标变换函数必须足够光滑—— 否则损失精度 61. 基本概念1） 控制体节点（中心）型控制体与网格型控制体 72） 积分型控制方程在控制体上积分物理含义： 控制体内总质量/动量/能量的增加 = 穿过控制体边界流入的净质量/动量/能量无粘通量粘性通量 8物理含义： 单位时间内，无粘流动流过垂直于n方向的单位面积的质量、动量和能量法向速度 93）有限体积法中物理量的含义4） 残差残差=净通量=右端项含义： 控制体内的平均量（平均质量密度、平均动量密度、平均能量密度）控制体几何中心处的当地密度、动量密度、能量密度二阶精度近似 102. 无粘通量的计算常用方法 （流过AB边的通量）： a. 利用周围点的值，计算出(I+1/2,J) 点处的物理量； b. 利用该处的物理量，计算出流过AB边的流通量方法1： 中心型有限体积法人工粘性项用中点的值代替AB面上的均值，二阶精度 11方法2： 迎风型有限体积法Step 1: 利用（偏）左侧点及（偏）右侧点的值，计算出I+1/2,J点的值计算方法： 与差分法完全相同各种差分格式，均可直接使用也称为“差分格式”该过程称为“重构”（很多文献中称为“插值”）有限体积与有限差分共通之处, 可直接使用差分格式 12常见的差分格式：2阶NND格式minmod(a,b) : a,b符号相反时取0， 符号相同时取绝对值小的3阶迎风3阶MUSCL格式TVD, WENO, GVC, 保单调格式…… 13重构方式： 原始变量、守恒变量及特征变量以NND格式为例：守恒变量重构原始变量重构特征变量重构先算出UI+1/2 (可用UI和UI+1的算术平均或Roe平均），再利用该值算出SI+1/2 是控制体内的平均值 （称为数值流通量） 的含义14重要概念澄清：重构与插值A. 有限体积法：j+1/2j-1/2确实为f在xj+1/2点的值 ！ 通常做法： 1） 用 计算出 2）u在xj+1/2点的值！关键： 是用 计算 （称为重构） ，而不是用 计算 （是标准的插值）；否则最高也只能达到2阶精度。 15重要概念澄清：重构与插值B. 有限差分法：j+1/2 h(x)为 f(x) 的“重构函数”差分法与有限体积法的数值格式完全相同 16Step 2: 利用 ， 计算出通量 方法1：流通矢量分裂 （FVS）方法2：通量差分分裂 (FDS)通过Riemann解，获得通量可利用近似Riemann阶（Roe, HLL, HLLC)Roe格式：xyx’y’方法3： AUSM方法压力项单独处理，其余项类似Van Leer 分裂 小技巧： 利用局部坐标系，计算通量时简化为1维问题AByx’y’“扩展的”一维问题x,y坐标系下的Riemann问题（二维）x’,y’ 坐标系下的Riemann问题（一维）x’, y’ 坐标系下，切向速度v’表现相当于被动标量 183. 粘性通量的计算关键问题： 计算 对 的导数同样计算19方法1. 利用Jacobian变换方法2. 利用Green积分公式计算最终公式与方法1相同 20 § 7.2 边界条件及网格分块1. 边界条件的处理方法常用方法： 虚网格设置1或2层虚网格（如果采用高精度格式，需设置更多重虚网格）虚网格点上的值根据边界条件给定 绝热固壁边界条件超声速入口边界条件超声速出口边界条件 21亚声速入口边界条件外部计算域内部边界特征分析： 内点提供一个边界条件，外部提供其余边界条件调节作用， 内点速度降低