cal_rvt.m

这几个算例都是MacCormack格式的,喷管计算用时间步进法,外突钝角绕流用空间步进法.都写成模块形式,便于查错.

function Cal_rVT(i);
global rou p T u v F1 F2 F3 F4;
global gama R;
 tA=F3(:,i).^2./(2.*F1(:,i))-F4(:,i);
 tB=gama./(gama-1).*F1(:,i).*F2(:,i);
 tC=-(gama+1)./(2.*(gama-1)).*F1(:,i).^3;
 rou(:,i)=(-tB+sqrt(tB.^2-4.*tA.*tC))./(2.*tA);
 u(:,i)=F1(:,i)./rou(:,i);
 v(:,i)=F3(:,i)./F1(:,i);
 p(:,i)=F2(:,i)-F1(:,i).*u(:,i);
 T(:,i)=p(:,i)./rou(:,i)./R;