计算程序_计算流体力学_对流方程_有限差分法_La格式_迎风格式_FTCS格式.docx

%一维对流方程迎风格式、Lax格式、FTCS格式差分法

计算

%潭花林 清华大学航天航空学院

%FTCS格式对于一维对流方程不稳定，最好不用clc

clearall

%1.参数定义dx=1;

x1=-18;x2=18;

x=x1:dx:x2;L1=length(x);

%dt=0.5*dx; %收敛

dt=2*dx; %不收敛t1=0;

t2=t1+80*dt;t=t1:dt:t2;L2=length(t);alpha=1;lambda=alpha*dt/dx;

geshi=1; %迎风格式

%geshi=2; %Lax格式

%geshi=3; %FTCS格式

%2.显式求解zeta=zeros(L1,L2);

forkk=1:3

geshi=kk;forii=1:L1

ifx(ii)0zeta(ii,1)=1;

else

ifx(ii)==0zeta(ii,1)=1/2;

else

end

end

end

ifx(ii)0zeta(ii,1)=0;

end

ifgeshi==1

forii=2:L1

forjj=1:(L2-1)

zeta(ii,jj+1)=zeta(ii,jj)-lambda*(zeta(ii,jj)-zeta(ii-1,jj));

endzeta(1,jj+1)=zeta(2,jj+1);

endzeta1=zeta;

elseifgeshi==2

forii=2:(L1-1)forjj=1:(L2-1)

zeta(ii,jj+1)=(zeta(ii+1,jj)+zeta(ii-1,jj))/2-.

..

lambda/2*(zeta(ii+1,jj)-zeta(ii-1,jj));end

zeta(1,jj+1)=zeta(2,jj+1);

zeta(L1,jj+1)=zeta(L1,jj)-lambda*(zeta(L1,jj)-zeta(L1-1,jj));

endzeta2=zeta;

elseifgeshi==3

forii=2:(L1-1)forjj=1:(L2-1)

zeta(ii,jj+1)=zeta(ii,jj)-lambda/2*(zeta(ii+1,jj)-zeta(ii-1,jj));

endzeta(1,jj+1)=zeta(2,jj+1);

zeta(L1,jj+1)=zeta(L1,jj)-lambda*(zeta(L1,jj)-zeta(L1-1,jj));

endend

end

zeta3=zeta;end

end

%3.绘图

%3.1 t=0

figure(1)n=1;

plot(x,zeta1(1:L1,n),-k,x,zeta2(1:L1,n),-.k

,x,zeta3(1:L1,n),--k)% 作图

%axisequal %%%

是否要求x、y坐标间距相等

%gridon %%%

是否要求画网格

xlabel(x/m),ylabel(t/s) %

%%x,y轴表示的变量含义

%text(1,2,f(x)) %%%

图中文字标识

legend(迎风格式,Lax格式,FTCS格式

) %%%不同曲线的线

型区分

title(t=0时刻的计算结果) %%%标题axis([-18,18,-0.2,1.2])

%3.2 t=10

figure(2)

n=(10-t(1))/dt;

plot(x,zeta1(1:L1,n),-k,x,zeta2(1:L1,n),-.k

,x,zeta3(1:L1,n),--k)% 作图

%axisequal %%%

是否要求x、y坐标间距相等

%gridon %%%

是否要求画网格

xlabel(x/m),ylabel(t/s) %

%%x,y轴表示的变量含义

%text(1,2,f(x)) %%%

图中文字标识

legend(迎风格式,Lax格式,FTCS格式

) %%%不同曲线的线

型区分

title(t=10s时刻的计算结果) %%%标题

%3.3 t=20

figure(3)

n=(20-t(1))/dt;

plot(x,zeta1(1:L1,n),-k,x,zeta2(1:L1,n),-.k

,x,zeta3(1:L1,n),--k)% 作图

%axisequal %%%

是否要求x、y坐标间距相等

%gridon %%%

是否要求画网格

xlabel(x/m),ylabel(t/s) %

%%x,y轴表示的变量含义

%text(1,2,f(x)) %%%

图中文字标识

legend(迎风格式,Lax格式,FTCS格式

) %%%不同曲线的线

型区分

title(