计算流体力学(中科院力学所)第10讲有限体积法2精品PPT课件.pptVIP

计算流体力学讲义 第十讲 有限体积法（2） ；力学所主楼219； ;知识回顾;1） 重构—— 两种不同的重构方案，向左偏及向右偏。 给出两种结果： 及; （称为数值流通量） 的含义; 是控制体内的平均值 ;概念：MUSCL与非MUSC类方法;10.1 近似Riemann解;;;以平均值 代替分布值;;最终界面的HLL通量如下：;10.1.2 HLLC Riemann近似求解器 （Toro）;;问题最终描述：求解超定方程—— 4个未知数，6个方程;近似解： 假设左、右波速已知，简化了计算 （求解线性方程） 精确解： 左、右波速都是压力的函数，方程复杂（非线性），需迭代求解;最终HLLC （方案1）的表达式为： ;HLLC 近似解( 方案2);HLLC 近似解 （方案3）; 10.1.3 波速的估算;Roe 平均：; 10.1.4 “扩展一维Riemann问题”的HLLC近似解;10.1.5 HLLC 方法计算通量的步骤 ;在有限体积法中的具体使用方法：;10.2 中心型有限体积法;;Jameson 人工粘性;10. 3 具体算例讲解——RAE2822翼型绕流;有限差分法的计算结果—— 5阶及7阶迎风差分+Steger-Warming FVS+ 3阶Runge-Kutta BL 湍流模型;1. 控制方程、网格及边界条件 控制方程—— 无量纲N-S方程 ;边界条件的处理方法：; 求解 ， 计算边界值，单边差分（重构）;2. 空间离散;2. 空间离散;;（2） 粘性通量的计算—— 采用中心格式; （1）;方案2： 利用积分关系式（“纯正的”有限体积法） ;3. 时间推进;4. 注意事项： 1）当网格不运动（变形）时，几何量不随时间变化，可事先计算出来并储存，以减小计算量。 2） 不要重复计算通量； 例：交界面(i+1/2,j)处的通量可供节点 （i,j）及（i+1,j）使用。 无需计算两次，否则会浪费计算量。 ;经常不断地学习,你就什么都知道。你知道得越多,你就越有力量 Study Constantly, And You Will Know Everything. The More You Know, The More Powerful You Will Be;结束语 当你尽了自己的最大努力时，失败也是伟大的，所以不要放弃，坚持就是正确的。 When You Do Your Best, Failure Is Great, So DonT Give Up, Stick To The End 演讲人：XXXXXX 时 间：XX年XX月XX日 计算流体力学讲义 第十讲 有限体积法（2） ；力学所主楼219； ;知识回顾;1） 重构—— 两种不同的重构方案，向左偏及向右偏。 给出两种结果： 及; （称为数值流通量） 的含义; 是控制体内的平均值 ;概念：MUSCL与非MUSC类方法;10.1 近似Riemann解;;;以平均值 代替分布值;;最终界面的HLL通量如下：;10.1.2 HLLC Riemann近似求解器 （Toro）;;问题最终描述：求解超定方程—— 4个未知数，6个方程;近似解： 假设左、右波速已知，简化了计算 （求解线性方程） 精确解： 左、右波速都是压力的函数，方程复杂（非线性），需迭代求解;最终HLLC （方案1）的表达式为： ;HLLC 近似解( 方案2);HLLC 近似解 （方案3）; 10.1.3 波速的估算;Roe 平均：; 10.1.4 “扩展一维Riemann问题”的HLLC近似解;10.1.5 HLLC 方法计算通量的步骤 ;在有限体积法中的具体使用方法：;10.2 中心型有限体积法;;Jameson 人工粘性;10. 3 具体算例讲解——RAE2822翼型绕流;有限差分法的计算结果—— 5阶及7阶迎风差分+Steger-Warming FVS+ 3阶Runge-Kutta BL 湍流模型;1. 控制方程