precal.m

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

function preCal(h)
global U1 U2 U3 F1 F2 F3 J2 plnA_px gama;
 F1(h,:)=U2(h,:);
 F2(h,:)=U2(h,:).^2./U1(h,:)+(gama-1)./(gama).*(U3(h,:)-(gama./2).*(U2(h,:).^2./U1(h,:)));
 F3(h,:)=gama.*(U2(h,:).*U3(h,:)./U1(h,:))-(gama.*(gama-1)./2).*(U2(h,:).^3./U1(h,:).^2);
 J2(h,:)=((gama-1)./gama).*(U3(h,:)-(gama./2).*(U2(h,:).^2./U1(h,:))).*plnA_px;