boundaryconditions.m

Matlab下实现的CFD算法

function BoundaryConditions(h,N,pN)
global U1 U2 U3 A gama;
 U1(h,1)=A(1);
 U2(h,1)=2*U2(h,2)-U2(h,3);
 U3(h,1)=U1(h,1)*(1/(gama-1)+(gama/2)*((U2(h,1)/U1(h,1))^2));

 U1(h,N)=2*U1(h,N-1)-U1(h,N-2);
 U2(h,N)=2*U2(h,N-1)-U2(h,N-2);
 VN=U2(h,N)/U1(h,N);
 U3(h,N)=pN*A(N)/(gama-1)+(gama/2)*U2(h,N)*VN;