高等计算流体力学-04学案.pptxVIP

第二部分 可压缩流动的数值方法;第四讲 有限体积方法初步;概述;研究内容：一维线性对流方程 的有限体积方法。 有限体积方法：对积分型的守恒方程进行离散，从而把积分型方程近似为代数方程进行求解的方法。 有限体积方法的出发方程是积分型的守恒方程。 对应的积分型守恒方程：微分型的守恒方程在有 限控制体上积分，利用Gauss公式，就可以得到积分型的守恒 方程。对于 而言，其守恒形式为 其中 。在控制体 上积分 式，可得积分型守恒 方程为;二、空间控制体;三、有限体积方法1：全离散形式;定义数值解在单元内的平均值为;注意到 即为了计算数值通量，我们必须求 得 之间的 。但是，由标准式，已知量只有 时刻单元平均值 。所以，为了求出 ，必须解决如下的两个问题： （1）由 ，求出 时刻数值解沿x方向的分布 ，这一过程称为重构（reconstruction）。 （2）由 时刻数值解沿x方向的分布 ，求出 之间的 以及数值通量 这一过程称为演化（evolution）过程。 有限体积方法中将在这两个过程中引入各种近似，从而把积分型方程化为代数方程。下面分别讨论这两个过程。 ;1. 重构; 现在 仍是未知的，确定 的 方法可以有很多种，其中之一是把 延拓到第 单元，要求其平均值为 ，即 由此得： 另一种方法是把 延拓到第 单 元，要求其平均值为 。 此时 当然， 也可以取上述两个方案的某种平均。更高阶的重构方法（如ENO/WENO格式 ）超出了本书的范围，在此不再介绍。这里我们指出，理论上，可以进行任意阶次的重构。;2 演化过程;（1）当采用零阶重构时;（2）当采用线性重构时，由式， 以及;;四、有限体积方法2：半离散形式;与全离散型的有限体积方法相同，半离散的有限体积方法中数值通量的构造也需要经过重构和演化过程，但是演化过程比全离散方法相对简单。;同理;在实际应用中，对半离散方法的时间导数项进行离散后，即可得到可以进行数 值计算的全离散形式。时间方向的离散可以采用前面介绍过的Runge-Kutta方法、 Crank-Nicolson方法等。;可以看到： 对于一维问题，有限体积方法和有限差分方法在一定意义上 （如 成立时）是等价的。 对于多维问题，当控制体为矩形（二维）或长方体（三维）时，二者在某种意义上也是等价的。 有限体积方法和有限差分方法有密切的联系。 有限差分方法的相容性、收敛性和稳定性概念以及稳定性分析方法，对于有限体积方法也是适用的。 当多维问题中，控制体不是矩形（二维）或长方体（三维）时，有限体积方法和有限差分方法不再等价。 一般认为，有限体积方法对于复杂的求解域形状有较好的适应性，而且也更适于计算包含间断的流场。;Euler方程的有限体积方法;积分型方程;时间积分;重构;重构的基本要求 守恒性 紧致性 由 及其邻近的有限控制体决定. K-exactness 当 在 及附近是次数小于等于k的多项式, 则:;简单重构 零阶重构(一阶精度) 零阶重构认为在每一个控制体内，物理量为常数，即: 一阶(线性)重构(二阶精度) 一阶或者线性重构认为每一个控制体内，物理量的分布为一次多项式,即 如果 位于控制体 的形心, 则;利用 , 则在充分