JessesnowCFD的博客

一、Matlab代码片%亚声速-超声速等熵喷管流动守恒形CFD解法 MacCormack方法clear; %清理内存变量clc; %清理工作窗中的所有显示内容r=1.4; %比热比L=3; %喷管长度i=31; %网格数目C=0.5; %科朗数x=linspace(0,L,i); %网格点横坐标N=1401; %时间步%预分配内存rho=zeros(1401,31); V=zeros(1401,31); T=zeros(1401,31);U1=zeros(1401,31); U2=z

一、Matlab代码片%全亚声速等熵喷管流动 非守恒型麦考马克方法数值求解clear; %清理内存变量clc; %清理工作窗中的所有显示内容r=1.4; %比热比L=3; %喷管长度i=31; %网格数目C=0.5; %科朗数x=linspace(0,L,i); %网格点横坐标N=5001; %时间步dx=L/(i-1); %空间步长dt(1:N)=0; %时间步长Pe=0.93; %喷管两端压力比A(1:31)=0; %预分配内存for j=1:31if(x(j)<1.

一、Matlab代码片%亚声速-超声速等熵喷管流动 非守恒型麦考马克方法数值求解clear; %清理内存变量clc; %清理工作窗中的所有显示内容r=1.4; %比热比L=3; %喷管长度i=31; %网格数目C=0.5; %科朗数x=linspace(0,L,i); %网格点横坐标N=1401; %时间步dx=L/(i-1); %空间步长dt(1:N)=0; %时间步长A=1+2.2*(x-1.5).*(x-1.5); %喷管面积rho(N,i)=0; %流场密度赋值T(N,i