precal.m

Matlab下实现的CFD算法

function preCal(h,N,Cx)
global U1 U2 U3 F1 F2 F3 J2 S1 S2 S3 A pA_px gama;
 F1(h,:)=U2(h,:);
 F2(h,:)=U2(h,:).^2./U1(h,:)+(gama-1)./(gama).*(U3(h,:)-(gama./2).*(U2(h,:).^2./U1(h,:)));
 F3(h,:)=gama.*(U2(h,:).*U3(h,:)./U1(h,:))-(gama.*(gama-1)./2).*(U2(h,:).^3./U1(h,:).^2);
 p(h,:)=(gama-1).*(U3(h,:)-(gama./2).*(U2(h,:).^2./U1(h,:)))./A;
 J2(h,:)=p(h,:).*pA_px./gama;
 for i=2:(N-1)
     temp=Cx*abs(p(h,i+1)-2*p(h,i)+p(h,i-1))/(p(h,i+1)+2*p(h,i)+p(h,i-1));
     S1(h,i)=temp*(U1(h,i+1)-2*U1(h,i)+U1(h,i-1));
     S2(h,i)=temp*(U2(h,i+1)-2*U2(h,i)+U2(h,i-1));
     S3(h,i)=temp*(U3(h,i+1)-2*U3(h,i)+U3(h,i-1));
 end