关于peric程序，高手帮一下

harry

对于peric的动量方程，对于u速度方程，没求用于压力方程的动量方程倒数项，而v速度却求了，即apr＝1/ap
在peric的直角坐标程序内就求了两个，是apu=1/ap和apv＝1/ap
看了半天，没想通，为什么呀？ 下边是对应他的程序部分
      SUBROUTINE CALCUV(K)
.....
C
C.....FINAL COEFFICIENT AND SOURCES MATRIX FOR U-EQUATION
C
      DO I=2,NIM
      DO IJ=LI(I+IST)+2,LI(I+IST)+NJM
        AP(IJ)=(AP(IJ)-AE(IJ)-AW(IJ)-AN(IJ)-AS(IJ))*URFU
        SU(IJ)=SU(IJ)+(1.-URF(IU))*AP(IJ)*U(IJ)
      END DO
      END DO
.....
C
C.....FINAL COEFFICIENT AND SOURCES MATRIX FOR V-EQUATION
C
      DO I=2,NIM
      DO IJ=LI(I+IST)+2,LI(I+IST)+NJM
        AP(IJ)=(APR(IJ)-AE(IJ)-AW(IJ)-AN(IJ)-AS(IJ))*URFV
        SU(IJ)=SV(IJ)+(1.-URF(IV))*AP(IJ)*V(IJ)
        APR(IJ)=1./(AP(IJ)+SMALL)  v方程有这个，为什么u方程就没有呢???????????????????????????
      END DO
      END DO
.......

anyone

没读过PERIC的程序.
让我猜的话, 如果U V处理不一样有可能是2维轴对称情况下, X(U), 和Y(V)方向处理不一样. 查一下那里用APR就会知道了.

wwwwwwwww

你好好把压力修正这块看看就明白了！还是有一点错！但不是你想的那样？wxgjf5110@yahoo.com.cn

xiaofish

一般来说, u和v的系数是相等的,所以只要求一次apr就可以了
当然, 即使u和v的系数不相等, 只用一个apr也是可以得到正确结果的
这是由于apr是求压力修正方程用的,当速度要收敛的时候, 压力修正值将趋近0

sandg

考古下
u,v的系数分别对应AP,APR,两个是相等的,
在压力泊松方程计算时,要用到它们的倒数,所以APR存贮它们的倒数,AP用于存贮压力校正方程的矩阵