corrector_step.m

Matlab下实现的CFD算法

function corrector_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+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+S1(h+1,:);
 U2(h+1,:)=U2(h,:)+pU2_pt(h,:).*dt+S2(h+1,:);
 U3(h+1,:)=U3(h,:)+pU3_pt(h,:).*dt+S3(h+1,:);