一维非稳态对流扩散问题有限体积

和一维非稳态扩散算例一样的初始和边界条件 ，采用乘方格式 时间步长为0.001s，初始温度场为200，速度为2m/s，长度为2cm，t=0s时刻东侧温度突然降至0C。时间差分采用全隐式格式

function onedicon

% —————————和一维非稳态扩散算例一样的初始和边界条件 p122页码
% —————————采用乘方格式 时间步长为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;
x