高斯热源 直接point flux值 和 face flux 子程序 结果为什么不同

吴宇哲

在热分析中， 给一组平面热源数据 （w/m^2) , 近似高斯分布 ， 首先我计算出我需要的node 平均数值 用ponit flux 方式计算 ；然后用 face flux 方式 用子程序计算 ， 但是两者计算结果有十倍的差异，请问 为什么这样？谢谢大家

子程序如下：

subroutine flux(f,temflu,mibody,time)

real z0,A,xc,w1,yc,w2,xx,yy
        z0=1.16394;        A=11.88473;        xc=25; w1=2.014; yc=10; w2=0.233


      if (bcname.eq.'fflux') then
                xx=-1*0.5*((temflu(4)-xc)/w1)**2
                yy=-1*0.5*((temflu(5)-yc)/w2)**2

                 f(1)=z0+A*exp(xx+yy)            
      end if
      return
      end