corrector_step.m

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

function corrector_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+dksi);
 for j=2:(M-1)
     if ksi<=OE 
         temp=0;
     else
         temp=(1-(j-1)*dyita)*tan(theta)/hx(ksi);
     end
     pF1_pksi(j,i+1)=temp*(F1(j-1,i+1)-F1(j,i+1))/dyita+(G1(j-1,i+1)-G1(j,i+1))/(th*dyita);
     pF2_pksi(j,i+1)=temp*(F2(j-1,i+1)-F2(j,i+1))/dyita+(G2(j-1,i+1)-G2(j,i+1))/(th*dyita);
     pF3_pksi(j,i+1)=temp*(F3(j-1,i+1)-F3(j,i+1))/dyita+(G3(j-1,i+1)-G3(j,i+1))/(th*dyita);
     pF4_pksi(j,i+1)=temp*(F4(j-1,i+1)-F4(j,i+1))/dyita+(G4(j-1,i+1)-G4(j,i+1))/(th*dyita);
 end
     pF1_pksi(1,i+1)=temp*(F1(1,i)-F1(2,i))/dyita+(G1(1,i)-G1(2,i))/(th*dyita);
     pF2_pksi(1,i+1)=temp*(F2(1,i)-F2(2,i))/dyita+(G2(1,i)-G2(2,i))/(th*dyita);
     pF3_pksi(1,i+1)=temp*(F3(1,i)-F3(2,i))/dyita+(G3(1,i)-G3(2,i))/(th*dyita);
     pF4_pksi(1,i+1)=temp*(F4(1,i)-F4(2,i))/dyita+(G4(1,i)-G4(2,i))/(th*dyita);

 pF1_pksi(:,i)=0.5.*(pF1_pksi(:,i+1)+pF1_pksi(:,i));
 pF2_pksi(:,i)=0.5.*(pF2_pksi(:,i+1)+pF2_pksi(:,i));
 pF3_pksi(:,i)=0.5.*(pF3_pksi(:,i+1)+pF3_pksi(:,i));
 pF4_pksi(:,i)=0.5.*(pF4_pksi(:,i+1)+pF4_pksi(:,i));

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