有限差分法

sbstoneface

有一很长的方形模具，外部尺寸为3030 mm, 内部为1010 mm，现需进行渗碳处理，以提高内表面耐磨性。现假定其内部充满渗碳剂，内表面碳浓度维持在1.4%，外表面为空气，碳浓度为0。求稳态时的碳浓度分布。
编程如下:
clear
space=10;
h=0.1;
time=2000000;
t=250;
n=(space/h)+1;
m=(time/t)+1;
c=zeros(n,m);
c(1,1:m)=1.4;
c(n,1:m)=0;
c(1:n,1)=0;
for j=1:m
      for i=2:n-1
      c(i,j)=0.5*(c(i-1,j)+c(i+1,j));
   end
end
w=linspace(0,space,n);
w1=c(:,m);
plot(w,w1,'m')
可运行结果是一条下降的曲线,应该是一斜直线吧,请高手指点!

happy

方程给一下

sbstoneface

方程为Laplace方程.X为距离,T为扩散时间,C为浓度