有限体积法介绍.docxVIP

有限体积法 有限体积法基本原理 上一章讲到的有限差分法将数值网格的节点上定义为计算节点，并在网格节点上对微分形式的流体基本方程进行离散，用网格节点上的物理量的代数方程作为原PDE 的近似。 在本章所要学习的有限体积法则采用了不同的离散形式。首先，有限体积法离散的是积 分形式的流体力学基本方程： ? ??v ? nds ? ? ?n ???ds ? ? S S ?  q d?? (1) ? 计算域用数值网格划分成若干小控制体。和有限差分法不同的是，有限体积法的网格定义了控制体的边界，而不是计算节点。有限体积法的计算节点定义在小控制体内部。一般有限体积法的计算节点有两种定义方法，一种是将网格节点定义在控制体的中心，另一种方法中，相邻两个控制体的计算节点到公共边界的距离相等。第一种方法的优点在于用计算节点的值作为控制体上物理量的平均值具有二阶的精度；第二种方法的好处是在控制体边界上的中心差分格式具有较高的精度。 积分形式的守恒方程在小控制体和计算域上都是成立的。为了获得每一个控制体上的代数方程，面积分和体积分需要用求面积公式来近似。 面积分的近似 采用结构化网格，在二维情况下，每一个控制体有4 个面，二维情况，每一个控制体有 6 个表面。计算节点用大写字母表示，控制体边界和节点用小写字母表示。为了保证守恒性， 控制体不能重叠，每一个面都是相邻两个控制体的唯一公共边界。 控制体边界上的积分等于控制体个表面的积分的和： ? fdS ? ?? fds S k Sk  (2) 上式中，f 可以表示 ?? u n 或? ?? 。 ?n NN NN NW nw N n NE ne WW W w sw P s e se E EE SW S SE SS 显然，为了获得边界上的积分，必须知道f 在边界上的详细分布情况，这是不可能实现的， 由于只是计算节点上的函数值，因此必须采用近似的方法来计算积分。 整个近似过程分成两步 第一步：用边界上几个点的近似积分公式 第二步：边界点上的函数值用计算节点函数值的插值函数近似面积分可采用以下不同精度的积分公式： 二阶精度积分： F ? ? fds ? f S e e e Se ? f S e e (3) 上式中 f e 为边界中点出的函数值。近似为方格中心点的值乘以方格的面积。 三阶精度积分： F ? ? fds ? e Se f ? f 2ne se S 2 e  (4) 四阶精度积分： F ? ? fds ? e Se f ? 4 f ? f ne e se S 6 e  (5) 应该注意的是，采用不同精度的积分公式，在相应的边界点的插值时也应采用相应精度的插值函数。积分公式的精度越高，近似公式就越复杂。 3 体积分的近似 和面积分相似，体积分也有不同精度的近似公式 二阶精度积分公式 Q ? ? qds ? qS e Se 采用双二次样条函数  ? q ?? P  (6) q( x, y) ? a ? a x ? a y ? a x 2 ? a y 2 ? a xy ? a x 2 y ? a xy 2 ? a x 2 y 2 0 1 2 3 4 5 6 7 8 (7) 可以得到四阶精度的积分公式： Q ? ? qds ? Se ?? ?16q 36 P 4q s 4q w 4q n 4q s 4q se 4q sw 4q ne 4q nw ? (8) 函数的插值 在上节讲到的积分的近似公式中用到了非计算节点上的函数值，被积函数 f 中包含了多个物理量及其偏微分，如对流项 f c ? ?? v ? n ,扩散项 f d ? ?n ? ?? ，在源项中也有类似情况，这里假定流场和流体的物性参数是已知的，物理量? 及其偏导数在控制面上的值需要通 过计算节点上物理量的插值得到。下面已e 面为例进行讨论。 迎风插值（UDS） ? 用上游计算节点的函数值近似相当于对一阶偏导数采用迎风格式，因此用 UDS 来表示这 e 种近似方法，在UDS 中： ?? if ?v ? n? ? 0 ??e? ? ? P e ? ? e E if ?v ? n? ? 0 e (9) UDS 是唯一无条件满足有界性要求的近似格式，在数值过程中不会产生数值振荡。UDS 存在 数值粘性。根据Taylor 公式，该格式具有一阶精度，并具有数值粘性： ? num ? ?? u ?  ?x / 2 (10) e e 在多维问题中，如果流动方向和网格是斜交的，截断误差会在垂直于流动方向以及流线方向产生扩散，这是一种非常严重的误差，函数的峰值或函数值的快速变化会被抹平，为了得到高精度结果需要采用非常精细的网格。 线性插值（CDS） ? ? ? ? e E E ? (1 ? ? )? E P (11