求大神指导程序关于求无粘通量
!内部流场单元---------------------------------------------------------------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网格中边界信息中第一二列是空间坐标,第三第四列呢?
回复 1# longming 的帖子
这是naca0012翼型GRD:53061568810382
0.0000000E+000.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=第四列,谢谢。
回复 2# longming 的帖子
问题描述:用有限体积法欧拉方程,求解可压缩无粘翼型绕流问题,就是关于求无粘通量q(u)。
页:
[1]