求大神指导程序关于求无粘通量

longming

!内部流场单元---------------------------------------------------------------
      if (kp>0) then         
                wi(1)=0.5*(w(1,kk)+w(1,kp))
                wi(2)=0.5*(w(2,kk)+w(2,kp))
                wi(3)=0.5*(w(3,kk)+w(3,kp))
                wi(4)=0.5*(w(4,kk)+w(4,kp))
                rhoi=wi(1)
                ui=wi(2)/rhoi
                vi=wi(3)/rhoi
        zi=ui*sij(2,i)-vi*sij(1,i)
        pi=(gamma-1.0)*(wi(4)-0.5*rhoi*(ui*ui+vi*vi))
                ci=sqrt(gamma*pi/rhoi)
        flux(1)=zi*rhoi
        flux(2)=zi*rhoi*ui+pi*sij(2,i)
        flux(3)=zi*rhoi*vi-pi*sij(1,i)
        flux(4)=zi*(wi(4)+pi)
        
        do j=1,4
          q(j,kk)=q(j,kk)+flux(j)
          q(j,kp)=q(j,kp)-flux(j)
        end do
其中w守恒量，请问该程序使用的是什么方法求解通量？在grd网格中边界信息中第一二列是空间坐标，第三第四列呢？

longming

这是naca0012翼型GRD:
5306  15688  10382
  0.0000000E+00  0.0000000E+00
  0.8000000E-02 -0.1531500E-01
  0.1600000E-01 -0.2127000E-01
...（以上是节点数据）
      8      7      1     -1
      7    231      1      9
      8    231      3      1
    145    144      2     -1
    144    232      2      4
    145    232     10      2
    231    233      3     11
（以上就是边界信息数据，前两列存的是边界位置，请问后两列是啥？）
其中kk=第三列，kp=第四列，谢谢。

longming

问题描述：用有限体积法欧拉方程，求解可压缩无粘翼型绕流问题，就是关于求无粘通量q(u)。