有限体积法介绍

1、.有限体积法1 有限体积法基本原理上一章讲到的有限差分法将数值网格的节点上定义为计算节点，并在网格节点上对微分形式的流体基本方程进行离散，用网格节点上的物理量的代数方程作为原PDE 的近似。在本章所要学习的有限体积法则采用了不同的离散形式。 首先，有限体积法离散的是积分形式的流体力学基本方程：v ndsndsq d ?(1)SS计算域用数值网格划分成若干小控制体。和有限差分法不同的是，有限体积法的网格定义了控制体的边界，而不是计算节点。有限体积法的计算节点定义在小控制体内部。一般有限体积法的计算节点有两种定义方法，一种是将网格节点定义在控制体的中心，另一种方法中，相邻两个控制体的计算节点到公共

2、边界的距离相等。第一种方法的优点在于用计算节点的值作为控制体上物理量的平均值具有二阶的精度；第二种方法的好处是在控制体边界上的中心差分格式具有较高的精度。积分形式的守恒方程在小控制体和计算域上都是成立的。 为了获得每一个控制体上的代数方程，面积分和体积分需要用求面积公式来近似。2 面积分的近似采用结构化网格，在二维情况下，每一个控制体有4 个面，二维情况，每一个控制体有6 个表面。计算节点用大写字母表示，控制体边界和节点用小写字母表示。为了保证守恒性，控制体不能重叠，每一个面都是相邻两个控制体的唯一公共边界。控制体边界上的积分等于控制体个表面的积分的和：fdSfds(2)SSkk上式中， f

3、可以表示un 或。n.NNNWNNEnwnneWWWwPeEEEswsseSWSSESS显然，为了获得边界上的积分，必须知道 f 在边界上的详细分布情况，这是不可能实现的，由于只是计算节点上的函数值，因此必须采用近似的方法来计算积分。整个近似过程分成两步第一步：用边界上几个点的近似积分公式第二步：边界点上的函数值用计算节点函数值的插值函数近似面积分可采用以下不同精度的积分公式：二阶精度积分：Fefds f eSe f e Se(3)Se上式中 f e 为边界中点出的函数值。近似为方格中心点的值乘以方格的面积。三阶精度积分：Fefdsf nef seSe(4)Se2四阶精度积分：Fefdsf n

4、e4 f ef se Se(5)Se6应该注意的是， 采用不同精度的积分公式， 在相应的边界点的插值时也应采用相应精度的插值函数。积分公式的精度越高，近似公式就越复杂。3 体积分的近似和面积分相似，体积分也有不同精度的近似公式二阶精度积分公式Qqdsq Se qP(6)Se采用双二次样条函数q( x, y) a0a1x a2 y a3 x2a4 y 2a5 xy a6 x 2 y a7 xy 2a8 x2 y 2(7).可以得到四阶精度的积分公式：Qqds16qP 4qs 4qw 4qn 4qs 4qse 4qsw 4qne 4qnw (8)Se364 函数的插值在上节讲到的积分的近似公式中用

5、到了非计算节点上的函数值，被积函数f 中包含了多个物理量及其偏微分，如对流项f cv n ,扩散项 f dn，在源项中也有类似情况，这里假定流场和流体的物性参数是已知的，物理量及其偏导数在控制面上的值需要通过计算节点上物理量的插值得到。下面已e 面为例进行讨论。4.1 迎风插值（ UDS ）e 用上游计算节点的函数值近似相当于对一阶偏导数采用迎风格式，因此用 UDS 来表示这种近似方法，在UDS 中：PeEifv nifv nee0(9)0UDS 是唯一无条件满足有界性要求的近似格式，在数值过程中不会产生数值振荡。UDS 存在数值粘性。根据Taylor 公式，该格式具有一阶精度，并具有数值粘性

6、：numu ex / 2(10)e在多维问题中， 如果流动方向和网格是斜交的，截断误差会在垂直于流动方向以及流线方向产生扩散， 这是一种非常严重的误差， 函数的峰值或函数值的快速变化会被抹平，为了得到高精度结果需要采用非常精细的网格。4.2 线性插值（ CDS ）eE E(1E ) P(11)xexP(12)ExPxE线性插值具有二阶精度，线性插值相当于FDM 中的 CDS 格式，因此用CDS 表示。 CDS 格式会产生数值振荡。对于扩散项EP(13)xexExP.4.3 三阶迎风格式（QUICK ）和 UDS 类似， QUICK 格式也和流动方向有关g1 Eg2 W(1 g1g2 ) Peg

7、4 EE(1 g3g 4 ) Eg3 P其中：ifv nifv nee0(14)0121e,P 12e,We, Pe,Wg1； g 21(15a)1e,Pe,We, Pe,W1e,W 122e,Pe,Ee, Pg3； g 41(15b)1e,Ee, Pe,Ee, P4.4 高阶格式（ 4 阶精度 CDS ）采用三次曲线可拟合出四阶精度的中心插值公式，在均匀网格中，四阶公式为：27 P27 E3 W3 EE(16)e4827E 27PWEE(17)x e24 x5 边界的处理对于对流项， 在入口处一般给出了流量或函数值， 在边界和对称面上流量为零， 在出口处假设和出口的法向坐标无关， 因此可采用

8、迎风格式。 对于扩散项则可能需要采用偏心格式。6 有限体积法应用举例例：考虑一标量在已知流场中的输运过程（如图4.4 所示），输运方程为：v ndSndS(18)SS边界条件：0 ；北部入口边界1y ；西部壁面边界对称条件；南部边界梯度为 0；东部出口条件.u xx ， u yy ，流线方程xyc入口，0(y)壁面对称边界对流项： Fecv ndSme eSemev ndSux ey 为质量通量。SeFecmax(me ,0) Pmin(me ,0) Efor UDSm（e1e )me E Efor CDSP若采用 UDS 格式，代数方程组中各项系数为：AEcmin( me ,0) ； AWc

9、min( mw ,0)ANcmin( mn ,0) ； AScmin( ms ,0)APc(AEcAWcANcASc )若采用 CDS 格式，代数方程组中各项系数为AEcme e ； AWcmw wANcmn n ； AScms sAPc(AEcAWcANcASc )根据连续性方程：memwm nms0相邻 CV 之间的关系：mw ,Pme,W ；w, P1e,W其余相邻 CV 有类似关系.流 线 ，xy=c出口n0(19)(20)(21)(22)(23)(24).扩散项采用CDS 格式FecSendSyyEP(25)x exE xP代数方程组中扩散项系数为：AEdxEy； AWdyxPxPx

10、WANdxNx； AWdxPx(26)xPxSAPd( AEdAWdANdASd )对于任意控制体AW WAS SAPP ANNAEE QP(27)Al AlcAld， l 为任意指标P，E，W， S，N。(28)边界条件的处理：对于西部和北部边界， 由于给定了函数值， 对流项可直接代入函数值而无需插值， 扩散项则采用一侧差分PW(29)xWxPxW这里， W 点和 P 的 w 边中点重合。南边和西边的梯度为零，以南边为例，由于梯度为零，PS ，代数方程变为：AWW(AS AP)P ANN AEE QP(30)6 SIMPLE方法考虑定常不可压流动问题，控制方程为：连续性方程：vndS0(31

11、)SV动量方程：vvndSn vdSpndSbd(32)SVSVSVCV不可压缩问题求解的困难在于压力场的求解。主要原因在于压力p 没有独立的方程组。先考虑一维问题：对于动量方程：.uu euupe pw(33)uu wxx ew若采用 CDS 格式uu Puu Euu Puu WuE uPuP uWpP pEpP pW22xx22简化后得：uu Euu Wu E uw2u PpEpW(34)2x2根据连续性方程，ui 1uiui 1c ，则有 pi 1pi 1 ，由于相邻节点之间的压力没有联系方程，容易造成压力交错现象。为了解决这一问题，可采用交错网格技术，即速度场和压力场采用不同的网格。以

12、二维问题为例，交错网格的布置如下图所示：NvnWuwPueEvsS主控制体为压力控制体（黑色实线网格）， u的控制体（红色虚线网格）的计算节点在主控制体的 e 边，控制体的 e， w 边界通过主控制体的计算节点，v 控制体（蓝色双点划线网格）的计算节点在主控制体的n 边，该控制体的n， s 面经过主控制体的计算节点。在 u 的控制体中，采用有限体积法离散可得u 的代数方程：ae ueanb unbQ ( pPpE ) Ae(35)anvnanb vnbQ ( p PpN ) Ae压力场的求解采用压力校正方法。即采用预估的压力场求速度，再用连续性方程校正压力场。当连续性方程得到满足时，压力场就是

13、真实的压力场。具体步骤如下1 预测压力场p2 将预测压力场代入动量方程，分别求解速度场u, vae ueanb unbQ ( pPpE ) Ae(36)anvnanb vnbQ ( p PpN ) Ae3 用连续性方程校正压力.设方程的精确解为u， v， puuu ； vvv ； ppp(37)其中 u , v , p 为校正量。则校正量满足方程：ae ueanb unbQ ( pPpE ) Ae(38)anvnanb vnbQ ( p PpN ) Ae略去相邻节点速度校正量的影响，可得：ueAepPpEdepPpE(39a)aevnAnpP pNd npPpN(39b)an代入连续性方程：(ueue ) Ae(uwuw )Aw( vnvn ) An(vs vs ) As0(40)整理得：de, P Ae, Pp PpEde,W Ae,WpWpPdn, P Ae,P pPpNdn, S An ,S pSpPue Ae,Puw Aw ,Pvn An, pvs A