initial_data.m

这几个算例都是MacCormack格式的,喷管计算用时间步进法,外突钝角绕流用空间步进法.都写成模块形式,便于查错.

function initial_data(M,N)
global rou p T u v F1 F2 F3 F4 G1 G2 G3 G4 SF1 SF2 SF3 SF4 iksi ksi dksi pF1_pksi pF2_pksi pF3_pksi pF4_pksi ;
global gama R OE OL theta;

gama=1.4;
R=287.06;
OE=10;
OL=65;
theta=15.352/180*pi;

ksi=0;
iksi=zeros(1,N);
rou=zeros(M,N);
p=zeros(M,N);
T=zeros(M,N);
V=zeros(M,N);
u=zeros(M,N);
v=zeros(M,N);

F1=zeros(M,N);
F2=zeros(M,N);
F3=zeros(M,N);
F4=zeros(M,N);

G1=zeros(M,N);
G2=zeros(M,N);
G3=zeros(M,N);
G4=zeros(M,N);

pF1_pksi=zeros(M,N);
pF2_pksi=zeros(M,N);
pF3_pksi=zeros(M,N);
pF4_pksi=zeros(M,N);

SF1=zeros(M,N);
SF2=zeros(M,N);
SF3=zeros(M,N);
SF4=zeros(M,N);

irou=1.23;
iT=286.1;
iMa=2;
init_condition(irou,iT,iMa);


function init_condition(irou,iT,iMa)
global rou T p u v F1 F2 F3 F4;
global gama R;
 rou(:,1)=irou;
 T(:,1)=iT;
 p(:,1)=R*irou*iT;
 u(:,1)=iMa*sqrt(gama*R*iT);
 % v(:,1)=0;
 F1(:,1)=rou(:,1).*u(:,1);
 F2(:,1)=rou(:,1).*u(:,1).^2+p(:,1);
 F3(:,1)=rou(:,1).*u(:,1).*v(:,1);
 F4(:,1)=(gama/(gama-1)).*p(:,1).*u(:,1)+rou(:,1).*u(:,1).*(u(:,1).^2+v(:,1).^2)./2;