有限体积法MATLAB求解程序

有限体积法MATLAB求解程序

function CDstead
% ===================一维稳态问题===================
% ===================采用中心差分并与解析解比较===================
t0=1;t1=0;
r=1;L=1;k=0.1;u=2.5;
n1=100;
dx=L/n1;
x=linspace（0Ln1+1）;
for i=1:n1
    xx（i）=（x（i）+x（i+1））/2;
end
F=r*u;
D=k/dx;
pe=F/D;
A=sparse（n1n1）;
b=zeros（n11）;
% ===================影响系数赋值===================
ae=D-F/2;
aw=D+F/2;
ap=ae+aw;
% ===================源项赋值===================
  sp1=-（2*D+F）;  sc1=2*D+F;   
  sp2=-（2*D-F）;sc2=2*D-F;
  b（1）=sc1*t0;
  b（n1）=sc2*t1;
% ===================对矩阵A主对角线赋值===================
for i=2:n1-1
    A（ii）=ap;
end
% ===================两边界赋值===================
   A（11）=ae-sp1;
   A（n1n1）=aw-sp2;
% ===================对角线赋值===================
for i=1:n1-1
   A（ii+1）=-ae;
   A（i+1i）=-aw;
end
T=A\b;
plot（T‘-*‘）;
hold on
%===================对应该问题参数下的精确解===================
d=（1-exp（F*xx/k））/（exp（F*L/k）-1）+1;
plot（d‘r‘）
title（‘一维稳态对流扩散问题‘）;
legend（‘中心差分解‘‘解析解‘）
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         743  2011-04-27 22:29  CDstead.asv
     文件        1060  2011-04-28 19:59  CDstead.m
     文件        1031  2011-05-09 15:51  CFstead.m
     文件        2125  2011-04-26 21:10  Copy_of_onedicon.m
     文件        4045  2011-04-29 12:19  Copy_of_secdim.asv
     文件        4030  2011-04-29 12:23  Copy_of_secdim.m
     文件          36  2011-04-26 18:28  matrixa.m
     文件      473647  2011-04-27 10:47  mccExcludedFiles.log
     文件        1978  2011-05-10 21:36  odim.m
     文件        2388  2011-04-27 22:39  onedicon.asv
     文件        2124  2011-05-09 15:53  onedicon.m
     文件        3990  2011-04-30 10:18  secdim.m
     文件        2932  2011-04-26 16:05  secdim2.asv
     文件        3952  2011-04-29 13:31  secdim2.m
     文件        3741  2011-05-10 22:16  sedim.m
     文件        1917  2011-05-10 22:22  yiwei.m