corrector_step.m

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

function corrector_step(N,h,dt);
global U1 U2 U3 F1 F2 F3 J2 dx pU1_pt pU2_pt pU3_pt gama;
 for i=2:(N-1)
     pU1_pt(h+1,i)=-((F1(h+1,i)-F1(h+1,i-1))/dx);
     pU2_pt(h+1,i)=-((F2(h+1,i)-F2(h+1,i-1))/dx)+J2(h+1,i);
     pU3_pt(h+1,i)=-((F3(h+1,i)-F3(h+1,i-1))/dx);
 end
 
 pU1_pt(h,:)=0.5.*(pU1_pt(h+1,:)+pU1_pt(h,:));
 pU2_pt(h,:)=0.5.*(pU2_pt(h+1,:)+pU2_pt(h,:));
 pU3_pt(h,:)=0.5.*(pU3_pt(h+1,:)+pU3_pt(h,:));
    
 U1(h+1,:)=U1(h,:)+pU1_pt(h,:)*dt;
 U2(h+1,:)=U2(h,:)+pU2_pt(h,:)*dt;
 U3(h+1,:)=U3(h,:)+pU3_pt(h,:)*dt;