precal.m

来自「这几个算例都是MacCormack格式的,喷管计算用时间步进法,外突钝角绕流用空」· M 代码 · 共 7 行

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;