使用HLL或HLLC型通量.PPT

Copyright by Li Xinliang * 边界条件的处理方法： 方案1 ：严格按照特征方向 （精确，但复杂） 外边界1 外边界2 边界 对于每个边界单元上讨论特征方向， 在法方向局部一维无粘化 边界处的物理量 特征值 物理含义 沿特征线守恒量 流特征 熵 流特征 切向速度 右波特征 右Riemann不变量 左波特征 左Riemann不变量 在边界法向的坐标系中，控制方程为近似为“扩展的1维Euler方程” 求解 ， 计算边界值，单边差分（重构） Copyright by Li Xinliang * 计算模型：4个特征波，分别以速度 在该一维直线上传播，根据传播方向独立讨论。 若传播方向指向边界内，则给定流动参数 ；否则，用内部点计算 边界 方案2： 近似处理 外边界1 ： 流场远离翼型，扰动很弱，认为 （无论亚、超均按超音速处理） 外边界2： 按超音速出流处理 外边界1 外边界2 近似 模型： 对于内流，有问题 固壁边界条件 物理条件： 补充计算条件（常用）： Copyright by Li Xinliang * 2. 空间离散 无粘通量 粘性通量 控制体上积分： 注： 选取控制体的方法很多 本计算采用的控制体——节点中心型 控制体的多种选择方式 Copyright by Li Xinliang * 2. 空间离散 无粘通量 粘性通量 （1） 无粘通量的计算 旋转到x’y’坐标系 （迎风）重构控制体边界处的左、右值 线性重构：平均值 =中心点值 求解“扩展”Riemann问题，得到控制体边界m处的通量 建议：使用HLL 或HLLC型通量 （公式见本PPT 12、17页） q 是x’轴与x轴的夹角 Copyright by Li Xinliang * 具体步骤: 无粘通量 讨论以i,j 为中心的控制体 在 所在的控制体表面（记为m）上 1） 重构自变量的左、右值 2）进行坐标旋转，得到坐标系（x’,y’）下的自变量 q 是x’轴与x轴的夹角 3） 求解Riemann问题，获得界面m上的通量 m 使用HLL 或HLLC型通量， 见12，17页 4） 在其他3个控制面（图中虚线所示）上，同样求出通量 5）4个面的通量求和，得到该控制体的总无粘通量 Copyright by Li Xinliang * （2） 粘性通量的计算—— 采用中心格式 可用周围角点的平均 关键问题： 如何计算 ？ 方案1： 借鉴有限差分法，引入坐标变换 具体表达式见: 任玉新等 《计算流体力学基础》 130-131页及127页 转化到计算坐标系 通常，令 计算坐标 实际上就是下标（i,j） （1） Copyright by Li Xinliang * 于是： 坐标变换系数的计算： 控制体的体积 半点上的Jocabian系数（体积加权）插值得到 完全的有限差分法， 坐标变换Jocabian系数的计算是全局的； 本方法中坐标变换Jocabian系数的计算是局部的，只用到周围几个点的坐标。Jocabian系数实际上是一组坐标值。 全部导数均按（1）式计算 Copyright by Li Xinliang * 方案2： 利用积分关系式（“纯正的”有限体积法） 如图示红色虚线控制体, 有： 假设内部 为常数分布，则 控制体的体积（面积） 可用周围四个点插值获得 其中： Copyright by Li Xinliang * 3. 时间推进 方案 1） 显格式， 例如3阶Runge-Kutta方法 方案2） 隐格式， LU-SGS （详见第8讲） 对于定常问题 （或LU-SGS的内迭代）可采用局部时间步长法，加快收敛速度 局部时间步长： 每个单元计算时，时间步长 可根据当地的网格及流动情况选取 局部时间步长法中间结果（未收敛时）是非物理的。 Copyright by Li Xinliang * 4. 注意事项： 1）当网格不运动（变形）时，几何量不随时间变化，可事先计算出来并储存，以减小计算量。 2） 不要重复计算通量； 例：交界面(i+1/2,j)处的通量可供节点 （i,j）及（i+1,j）使用。 无需计算两次，否则会浪费计算量。 Copyright by Li Xinliang * 作业 10.1 计算