外突钝角绕流说明.txt

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

主程序：PMEW.m
计算采用MacCormack格式，空间步进法，给定来流初始条件，壁面边界条件采用Abbett`s boundary，采用贴体网格，加入人工粘性，计算未使用无量纲处理。（输出为计算空间）
运行后可以生成几个全局变量
   rou 密度
   p   压力
   T   温度
   u   x方向速度
   v   y方向速度
   iksi 计算的各个步长
举例
>> PMEW

time =

    3.2550

>> mesh(u)


子模块说明：
   function time=MC_SM_2D_C(M,N,dyita,C,Cy)
   %MacCormack守恒型外突钝角绕流主程序
   %M,N是划分网格数
   %dyita是yita所选取步长
   %C是计算ksi步长的Courant Number
   %Cy是人工粘性系数
  
   function initial_data(M,N)
   %初始程序
     function init_condition(irou,iT,iMa)
     %处理初始条件

   function cal_dksi(M,i,dyita,C); 
   %计算ksi步长

   function preCal(M,i,Cy)
   %处理计算用到的计算变量

   function predictor_step(M,i,dyita);
   %预估

   function Cal_rVT(i);
   %计算 rou V T

   function BoundaryConditions(M,i)
   %边界条件

   function corrector_step(M,i,dyita);
   %校正