predictor_step.m

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

function predictor_step(N,h,dt);
global U1 U2 U3 F1 F2 F3 J2 dx pU1_pt pU2_pt pU3_pt;
 for i=2:(N-1)
     pU1_pt(h,i)=-((F1(h,i+1)-F1(h,i))/dx);
     pU2_pt(h,i)=-((F2(h,i+1)-F2(h,i))/dx)+J2(h,i);
     pU3_pt(h,i)=-((F3(h,i+1)-F3(h,i))/dx); 
 end
 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;