喷管计算说明.txt

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

主程序：NFlow.m
计算采用MacCormack格式，时间步进法，给定初始条件，入口和出口边界条件，计算使用了无量纲处理。
运行后可以生成几个全局变量
        rou 密度
        T   温度
        V   速度
        A   喷管型面
        time 运行时间
举例
>> NFlow

time =

   30.3530

>> plot(rou(:,31))
>> plot(rou(1600,:))


子模块说明：
        function time=MC_TM_1D_NC(N,dt,t)
        %MacCormack非守恒型无激波跨声速喷管主程序
        %N是划分网格点数
        %dt是时间步长
        %t是计算的时间格点数 
        %time返回计算所用时间

        function time=MC_TM_1D_NC(N,dt,t,PN)
        %MacCormack非守恒型无激波亚声速喷管主程序
        %N是划分网格点数
        %dt是时间步长
        %t是计算的时间格点数 
        %PN出口压力
        %time返回计算所用时间

        function time=MC_TM_1D_C(N,dt,t,tU2)
        %MacCormack守恒型无激波跨声速喷管主程序
        %N是划分网格点数
        %dt是时间步长
        %t是计算的时间格点数 
        %tU2计算初始流量
        %time返回计算所用时间

        function time=MC_TM_1D_C(N,dt,t,pN,Cx)
        %MacCormack守恒型有激波跨声速喷管主程序
        %N是划分网格点数
        %dt是时间步长
        %t是计算的时间格点数 
        %pN是出口压力
        %Cx人工粘性系数
        %time返回计算所用时间

%子模块的变量对各个计算程序略有不同

 
        function initial_data(t,N)
        %初始程序
           function y=NF_Shape(N)
           %面积函数
           function init_condition(N)
           %初始条件
        
        function preCal(h,dt)
        %处理计算用到的计算变量

        function predictor_step(N,h,dt);
        %预估

        function BoundaryConditions(h,N)
        %边界条件

        function corrector_step(N,h,dt);
        %校正

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