boundaryconditions.m

来自「这几个算例都是MacCormack格式的,喷管计算用时间步进法,外突钝角绕流用空」· M 代码 · 共 10 行

function BoundaryConditions(h,N)
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);
 U3(h,N)=2*U3(h,N-1)-U3(h,N-2);