cal_dksi.m

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

function cal_dksi(M,i,dyita,C); 
global u v T dksi ksi iksi;
global gama R theta OE;
  tMa=sqrt(u(:,i).^2+v(:,i).^2)./sqrt(gama.*R.*T(:,i));
  miu=asin(1./tMa);
  if ksi<=OE
      ttheta=atan(v(:,i)./u(:,i));
      temp=dyita*M;
  else
      ttheta=theta+atan(v(:,i)./u(:,i));
      temp=dyita*M+dyita*(ksi-OE)*tan(theta);
  end
  tdksi1=C.*temp./(abs(tan(ttheta+miu)));
  tdksi2=C.*temp./(abs(tan(ttheta-miu)));
  dksi1=min(tdksi1);
  dksi2=min(tdksi2);
  dksi=min(dksi1,dksi2);
  iksi(i)=dksi;
  ksi=ksi+dksi;