mc_sm_2d_c.m

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

function time=MC_SM_2D_C(M,N,dyita,C,Cy)
% MACCORMACK`S TECHNIQUE 
% conservation Form
% Space Marching

tic;
initial_data(M,N);
for i=1:N-1
    cal_dksi(M,i,dyita,C);
    
    preCal(M,i,Cy);
    predictor_step(M,i,dyita);
    Cal_rVT(i+1);
    BoundaryConditions(M,i+1);
    
    preCal(M,i+1,Cy);
    corrector_step(M,i,dyita);
    Cal_rVT(i+1);
    BoundaryConditions(M,i+1);
    
end
time=toc;