|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
!内部流场单元---------------------------------------------------------------
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网格中边界信息中第一二列是空间坐标,第三第四列呢? |
|