高手帮忙解释一下peric的程序啊

Combustor

高手帮忙解释一下peric的程序啊
2dgl文件夹下面的程序考虑了网格的非正交性，但是看不懂他的边界条件的处理。就是这段，感觉壁面切应力有重复算的可能，SU(IJP)项算得是显式的壁面切应力？AP(IJP) =AP(IJP) +COEF*XTW(IW)**2也在算壁面切应力的X 分量？请高手帮忙看看。
C.....WALL BOUNDARIES
C
      DO IW=IWS(K)+1,IWS(K)+NWAL(K)
        IJP=IJPW(IW)
        IJB=IJW(IW)
        COEF=VIS(IJB)*SRDW(IW)
        AP(IJP) =AP(IJP) +COEF*XTW(IW)**2
        APR(IJP)=APR(IJP)+COEF*YTW(IW)**2
        SU(IJP)=SU(IJP)+COEF*(U(IJB)*XTW(IW)**2-
     *         (V(IJP)-V(IJB))*XTW(IW)*YTW(IW))
        SV(IJP)=SV(IJP)+COEF*(V(IJB)*YTW(IW)**2-
     *         (U(IJP)-U(IJB))*XTW(IW)*YTW(IW))
      END DO

Combustor

论坛上原来有个高手啊，怎么不在了呢 呼唤高手啊

liuhuafeifei

仅供参考

liuhuafei