predictor_step.m

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

function predictor_step(M,i,dyita);
global F1 F2 F3 F4 G1 G2 G3 G4 SF1 SF2 SF3 SF4 pF1_pksi pF2_pksi pF3_pksi pF4_pksi ksi dksi;
global theta OE;
 th=hx(ksi);
 for j=1:(M-1)
     if ksi<=OE 
         temp=0;
     else
         temp=(1-(j-1)*dyita)*tan(theta)/th;
     end
     pF1_pksi(j,i)=temp*(F1(j,i)-F1(j+1,i))/dyita+(G1(j,i)-G1(j+1,i))/(th*dyita);
     pF2_pksi(j,i)=temp*(F2(j,i)-F2(j+1,i))/dyita+(G2(j,i)-G2(j+1,i))/(th*dyita);
     pF3_pksi(j,i)=temp*(F3(j,i)-F3(j+1,i))/dyita+(G3(j,i)-G3(j+1,i))/(th*dyita);
     pF4_pksi(j,i)=temp*(F4(j,i)-F4(j+1,i))/dyita+(G4(j,i)-G4(j+1,i))/(th*dyita);
 end

 F1(:,i+1)=F1(:,i)+pF1_pksi(:,i)*dksi+SF1(:,i);
 F2(:,i+1)=F2(:,i)+pF2_pksi(:,i)*dksi+SF2(:,i);
 F3(:,i+1)=F3(:,i)+pF3_pksi(:,i)*dksi+SF3(:,i);
 F4(:,i+1)=F4(:,i)+pF4_pksi(:,i)*dksi+SF4(:,i);