资源简介

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

资源截图

代码片段和文件信息

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


评论

共有 条评论