matlab求解二维导热与对流问题

有限体积法求解二维导热与对流问题

function onedicon

% 和一维非稳态扩散算例一样的初始和边界条件
% 采用乘方格式 时间步长为0.001s
% 初始温度场为200，速度为2m/s，长度为2cm，t=0s时刻东侧温度突然降至0C。
% 时间差分采用全隐式格式
clc
disp（‘一维无内热源非稳态全隐 乘方格式离散，初始场为200C，右侧边界温度变为0，左侧边界条件可选择‘）;
disp（‘------------------‘）
a=input（‘选择左边界条件：1—第一类边界条件，2-第二类边界条件 :   ‘）;
q=0; t1=0;
if（a==0）
    q=input（‘热流为 ： ‘）;
else
    t1=input（‘选择左边界第一类条件温度，取值在0-200内： ‘）;
end
t2=input（‘选择右边界第一类条件温度，取值在0-200内：  ‘）;
u=input（‘please input the velocity :   ‘）;
n=input（‘长度2cm上划分网格数 :   ‘）;
dt=input（‘时间步长 :   ‘）;
%几何条件，物性参数
L=0.02;r=1.0;k=0.03;dx=L/n;
F=r*u;D=k/dx;pe=F/D;
TT=ones（n+11）;
%没有内热源
sp=0;
%初始温度场
tp0=200*ones（n1）;
x=linspace（00.02n+1）;
for i=1:n
     xx（i）=（x（i）+x（i+1））/2;
end
xxx（n+1）=0.02;
xxx（1:n）=xx;
%乘方格式赋值影响系数
ae=D*max（0（1-0.1*abs（pe^5）））+max（-F0）;
aw=D*max（0（1-0.1*a

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1758  2013-06-15 16:02  onedimensional.asv

     文件       1758  2013-06-15 16:05  onedimensional.m

     文件       3125  2013-06-15 16:02  twodimensional.asv

     文件       3125  2013-06-15 16:05  twodimensional.m

----------- ---------  ---------- -----  ----

                 9766                    4