cal_rvt.m

Matlab下实现的CFD算法

function Cal_rVT(h);
global rou V T U1 U2 U3 A gama;
 rou(h+1,:)=U1(h+1,:)./A;
 V(h+1,:)=U2(h+1,:)./U1(h+1,:);
 T(h+1,:)=(gama-1).*(U3(h+1,:)./U1(h+1,:)-(gama/2).*V(h+1,:).^2);