predictor_step.m

来自「Matlab下实现的CFD算法」· M 代码 · 共 12 行

function predictor_step(N,h,dt);
global U1 U2 U3 F1 F2 F3 J2 S1 S2 S3 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+S1(h,:);
 U2(h+1,:)=U2(h,:)+pU2_pt(h,:).*dt+S2(h,:);
 U3(h+1,:)=U3(h,:)+pU3_pt(h,:).*dt+S3(h,:);