1、1有限体积法1 有限体积法基本原理上一章讲到的有限差分法将数值网格的节点上定义为计算节点，并在网格节点上对微分形式的流体基本方程进行离散，用网格节点上的物理量的代数方程作为原 PDE 的近似。在本章所要学习的有限体积法则采用了不同的离散形式。首先，有限体积法离散的是积分形式的流体力学基本方程：(1)dqsdsSS nv计算域用数值网格划分成若干小控制体。和有限差分法不同的是，有限体积法的网格定义了控制体的边界，而不是计算节点。有限体积法的计算节点定义在小控制体内部。一般有限体积法的计算节点有两种定义方法，一种是将网格节点定义在控制体的中心，另一种方法中，相邻两个控制体的计算节点到公共边界的距离
2、相等。第一种方法的优点在于用计算节点的值作为控制体上物理量的平均值具有二阶的精度；第二种方法的好处是在控制体边界上的中心差分格式具有较高的精度。积分形式的守恒方程在小控制体和计算域上都是成立的。为了获得每一个控制体上的代数方程，面积分和体积分需要用求面积公式来近似。2 面积分的近似采用结构化网格，在二维情况下，每一个控制体有 4 个面，二维情况，每一个控制体有 6 个表面。计算节点用大写字母表示，控制体边界和节点用小写字母表示。为了保证守恒性，控制体不能重叠，每一个面都是相邻两个控制体的唯一公共边界。控制体边界上的积分等于控制体个表面的积分的和：(2)kSSfdsf上式中，f 可以表示 或 。
3、nu2P ENWSNNEEWWSSNW NESW SEewnsnenwsesw显然，为了获得边界上的积分，必须知道 f 在边界上的详细分布情况，这是不可能实现的，由于只是计算节点上的函数值，因此必须采用近似的方法来计算积分。整个近似过程分成两步第一步：用边界上几个点的近似积分公式第二步：边界点上的函数值用计算节点函数值的插值函数近似面积分可采用以下不同精度的积分公式：二阶精度积分：(3)eeSe SffdsFe上式中 为边界中点出的函数值。近似为方格中心点的值乘以方格的面积。三阶精度积分：(4)esneSe SffdsFe2四阶精度积分：(5)esenSe ffse64应该注意的是，采用不同精
4、度的积分公式，在相应的边界点的插值时也应采用相应精度的插值函数。积分公式的精度越高，近似公式就越复杂。3 体积分的近似和面积分相似，体积分也有不同精度的近似公式二阶精度积分公式(6)PeSqdsQe采用双二次样条函数 2827265243210),( yxayxayaxyxayxq (7)3可以得到四阶精度的积分公式：(8) nweswsenwsPS qqqqqdsQe 444163 4 函数的插值在上节讲到的积分的近似公式中用到了非计算节点上的函数值，被积函数 f 中包含了多个物理量及其偏微分，如对流项 ,扩散项 ，在源项中也有类nvcf ndf似情况，这里假定流场和流体的物性参数是已知的，
5、物理量 及其偏导数在控制面上的值需要通过计算节点上物理量的插值得到。下面已 e 面为例进行讨论。4.1 迎风插值（UDS）用上游计算节点的函数值近似相当于对一阶偏导数采用迎风格式，因此用 UDS 来表示e这种近似方法，在 UDS 中：(9)0eEPeifnvUDS 是唯一无条件满足有界性要求的近似格式，在数值过程中不会产生数值振荡。UDS 存在数值粘性。根据 Taylor 公式，该格式具有一阶精度，并具有数值粘性：(10)2/xuenme在多维问题中，如果流动方向和网格是斜交的，截断误差会在垂直于流动方向以及流线方向产生扩散，这是一种非常严重的误差，函数的峰值或函数值的快速变化会被抹平，为了得
6、到高精度结果需要采用非常精细的网格。4.2 线性插值（CDS）(11)PEEe)1(12)PEx线性插值具有二阶精度，线性插值相当于 FDM 中的 CDS 格式，因此用 CDS 表示。CDS格式会产生数值振荡。对于扩散项(13)PEexx44.3 三阶迎风格式（QUICK）和 UDS 类似， QUICK 格式也和流动方向有关(14) 0)1(4343 221 eEEPPWe ifggnv其中：； (15a)WePg,21WePe,21； (15b)PeE,23PeEg,244.4 高阶格式（4 阶精度 CDS）采用三次曲线可拟合出四阶精度的中心插值公式，在均匀网格中，四阶公式为：(16)483
7、27EWEPe (17)xxPe5 边界的处理对于对流项，在入口处一般给出了流量或函数值，在边界和对称面上流量为零，在出口处假设和出口的法向坐标无关，因此可采用迎风格式。对于扩散项则可能需要采用偏心格式。6 有限体积法应用举例例：考虑一标量在已知流场中的输运过程（如图 4.4 所示） ，输运方程为：(18)SSddnnv边界条件：；北部入口边界0；西部壁面边界y1对称条件；南部边界5梯度为 0；东部出口条件， ，流线方程xuycxy(y)壁面对称边界入口，0出口 0n流线，xy=c对流项： (19)eScemdFe nv为质量通量。yumxe (20)CDSfor)1U)0,i(0,ax(Ee
8、Pece若采用 UDS 格式，代数方程组中各项系数为：；)0,in(ecEmA)0,in(wcWmA； (21)NsS)(cNcWEcP若采用 CDS 格式，代数方程组中各项系数为；ecEmAwc； (22)nNsS)(cSNcWEcPA根据连续性方程：(23)0snwem相邻 CV 之间的关系：； (24)WeP,WeP,16其余相邻 CV 有类似关系扩散项采用 CDS 格式(25)PEeSce xyxdFe n代数方程组中扩散项系数为：；PEdxyAWPdxyA； (26)PNdSPd)(dNdWEdPAA对于任意控制体(27)PEPSWQ，l 为任意指标 P，E，W ，S ，N。 (28
9、)dcll边界条件的处理：对于西部和北部边界，由于给定了函数值，对流项可直接代入函数值而无需插值，扩散项则采用一侧差分(29)WPxx这里，W 点和 P 的 w 边中点重合。南边和西边的梯度为零，以南边为例，由于梯度为零， ，代数方程变为：SP(30)ENPS QAA)(6 SIMPLE 方法考虑定常不可压流动问题，控制方程为：连续性方程：(31)0SVdnv动量方程：(32) CVSVSVSV dpdbnv不可压缩问题求解的困难在于压力场的求解。主要原因在于压力 p 没有独立的方程组。先考虑一维问题：对于动量方程：7(33) wewewe pxuu若采用 CDS 格式 222 WPEPWPE
10、WPEP pxuxu 简化后得：(34) 2WEPwEWE pxuu根据连续性方程， ，则有 ，由于相邻节点之间的压力没有ciii11 1ii联系方程，容易造成压力交错现象。为了解决这一问题，可采用交错网格技术，即速度场和压力场采用不同的网格。以二维问题为例，交错网格的布置如下图所示：P EWNSueuwvnvs主控制体为压力控制体（黑色实线网格） ，u 的控制体（红色虚线网格）的计算节点在主控制体的 e 边，控制体的 e，w 边界通过主控制体的计算节点，v 控制体（蓝色双点划线网格）的计算节点在主控制体的 n 边，该控制体的 n，s 面经过主控制体的计算节点。在 u 的控制体中，采用有限体积
11、法离散可得 u 的代数方程：(35)eNPnbn Ee ApQvau)(压力场的求解采用压力校正方法。即采用预估的压力场求速度，再用连续性方程校正压力场。当连续性方程得到满足时，压力场就是真实的压力场。具体步骤如下1 预测压力场 p2 将预测压力场代入动量方程，分别求解速度场 vu,8(36)eNPnbn Ee ApQvau)(3 用连续性方程校正压力设方程的精确解为 u，v，p； ； (37)u其中 为校正量。则校正量满足方程：,(38) eNPnbn Ee ApQva)(略去相邻节点速度校正量的影响，可得：(39a)EPeEPedpAu(39b)NPnNPnpav代入连续性方程：(40)0)()()()( ssnnwee AvvAuu整理得： PspnPwe PSnNPenWeEPvA pdpddpd, , (41)PspnPwe SneWePSnNAvupdd, , , 求解压力校正方程(41)可得压力校正量，4 用(39)校正速度场。5 以 p 为猜测值，重复 15 直至收敛6 计算其他物理量如温度场等。
