precal.m

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

function preCal(M,i,Cy)
global rou T p u v F1 F2 F3 F4 G1 G2 G3 G4 SF1 SF2 SF3 SF4;
global gama R;


  G1(:,i)=rou(:,i).*v(:,i);
  G2(:,i)=rou(:,i).*u(:,i).*v(:,i);
  G3(:,i)=rou(:,i).*v(:,i).^2+p(:,i);
  G4(:,i)=(gama/(gama-1)).*p(:,i).*v(:,i)+rou(:,i).*v(:,i).*(u(:,i).^2+v(:,i).^2)./2;
    
  for j=2:(M-1)
      SF1(j,i)=Cy*abs(p(j+1,i)-2*p(j,i)+p(j-1,i))/(p(j+1,i)+2*p(j,i)+p(j-1,i))*(F1(j+1,i)-2*F1(j,i)+F1(j-1,i));
      SF2(j,i)=Cy*abs(p(j+1,i)-2*p(j,i)+p(j-1,i))/(p(j+1,i)+2*p(j,i)+p(j-1,i))*(F2(j+1,i)-2*F2(j,i)+F2(j-1,i));
      SF3(j,i)=Cy*abs(p(j+1,i)-2*p(j,i)+p(j-1,i))/(p(j+1,i)+2*p(j,i)+p(j-1,i))*(F3(j+1,i)-2*F3(j,i)+F3(j-1,i));
      SF4(j,i)=Cy*abs(p(j+1,i)-2*p(j,i)+p(j-1,i))/(p(j+1,i)+2*p(j,i)+p(j-1,i))*(F4(j+1,i)-2*F4(j,i)+F4(j-1,i));
  end