计算流体力学2011cfd2011-8讲有限体积法.pptx

计算流体力学讲义2011第八讲有限体积法（2）李新亮；力学所主楼219；知识点：1讲义、课件上传至(流体中文网）-“流体论坛”-“CFD基础理论”讲课录像及讲义上传至网盘CopyrightbyLiXinliang通量分裂技术简述加速收敛技术

CopyrightbyLiXinliang2知识回顾：有限体积法基本流程无粘项常用方法（流过AB边的通量）：a.利用周围点的值，计算出(I+1/2,J)点处的物理量；直接利用“差分格式”b.利用该处的物理量，计算出流过AB边的流通量迎风型方法需利用“通量分裂技术”FVS类：FDS类：利用Riemann解Reimann解：Godunov,Roe,HLL,HLLC利用坐标变换，转化为一维Riemann问题

§8.1通量分裂技术简述多维问题可通过坐标旋转，变为局部扩展一维问题该情况下，v,w表现为被动标量新坐标系下的控制方程：含（一个或多个）被动标量的一维Euler方程被动标量与一维问题相比，增加了（一个或多个）被动标量，对方程性质没有影响。被动标量方程分裂形成简单

CopyrightbyLiXinliang41.流通矢量分裂(FVS)1)Steger-Warming分裂1维：3维：

CopyrightbyLiXinliang52)Lax-Friedrichs(L-F)分裂=+优点：简单，计算量小缺点：耗散偏大足够大L-F分裂=+优点：耗散小缺点：导数间断S-W分裂

CopyrightbyLiXinliang63）VanLeer分裂根据当地Mach数分裂亚声速情况下，均匀过渡-11

CopyrightbyLiXinliang72.通量差分分裂（FDS）利用Riemann解，计算通量1）精确Riemann解（Godunov方法）v,w按照被动标量处理

满足：物理意义为平均增长率CopyrightbyLiXinliang82）Roe近似Riemann解uf(u)uLuRuRoeRoe平均Riemann问题：近似：用平均增长率替代瞬时增长率常系数线性方程组，求解简单（相似变换解耦求解）

CopyrightbyLiXinliang93）HLL近似Riemann解（Harten,LaxvanLeer）Ref.：E.F.Toro:RiemannSolversandNumericalMethodsforFluidDynamics,Springer,2009(ThirdEdition)基本原理：双激波近似t=0t=t0激波1，速度Z1激波2，速度Z2假设间断面产生两道激波，速度分别为Z1,Z2根据质量、动量、能量守恒，容易计算出图中控制体积内的总质量、总动量、总能量t0时刻激波才传到控制体边界，因此0到t0时刻，控制体边界处物理量保持0时刻的值。利用总量，求出图中控制体内的平均值，作为该区域物理量的近似值

CopyrightbyLiXinliang104）HLLC近似Riemann解（Toro）发展了HLL近似解，用三波模型来近似（如图）三波近似，左、右波的速T时刻的流动状态激波激波接触间断模型：左右两道激波，中间有接触间断激波速度已知为:ZL,ZR未知数（4个）：方程（6个）：两道激波的RH关系式方程多了两个？（因为假设激波速度已知）常用方法：去掉两个方程去掉两个能量方程，4个未知数，4个方程，求解求解过程简单，轻易可给出表达式

CopyrightbyLiXinliang11最终，HLLC公式为：三波近似，左、右波的速T时刻的流动状态激波激波接触间断接触间断移动速度Ref.：E.F.Toro:RiemannSolversandNumericalMethodsforFluidDynamics,Springer,2009(ThirdEdition)

CopyrightbyLiXinliang12HLL及HLLC均假设ZL,ZR已知，实际上它们仍需要估算准确计算ZL,ZR，实际是计算Riemann精确解，计算量大方法1：直接估算1a:假设以声速传播（Davis）竟然假设激波以声速传播，太OUT了小常识：激波的传播速度激波相对于波前介质以超声速传播，相对于波后介质以亚声速传播；弱激波（Ma趋近于1）以声速传