initial_data.m

Matlab下实现的CFD算法

function initial_data(t,N)
global rou T V x dx A prou_pt pV_pt pT_pt;
 dx=3/(N-1);
 x=0:dx:3;
 A=NF_Shape(N);
 rou=zeros(t,N);
 T=zeros(t,N);
 V=zeros(t,N);
 init_condition(N);
 prou_pt=zeros(t,N);
 pV_pt=zeros(t,N);
 pT_pt=zeros(t,N);


function y=NF_Shape(N)
global x;
  Nd2=floor(N/2);
  y=linspace(0,0,N);
  y(1:Nd2)=1+2.2.*(x(1:Nd2)-1.5).^2;
  y(Nd2+1:N)=1+0.2223.*(x(Nd2+1:N)-1.5).^2;

  
function init_condition(N)
global rou T V x;
  rou(1,:)=1-0.023*x;
  T(1,:)=1-0.009333*x;
  V(1,:)=0.05+0.11*x;