cal_rvt.m

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

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);