alfa_convective.m

ART1算法,此算法还存在缺陷,有此产生art2

function a=alfa_convective(b,c)
global t p lamta_s rou_s cp_s t_a0 n t_an;
g=9.8;
b=1;c=1;%调试
if b==1 
    temp_t=(t(b,c)+t_a0)/2;
elseif b==n
    temp_t=(t(b,c)+t_an)/2;
else
    temp_t=0;
end
rou_f=rou_air(temp_t);%kg/m^3
cp_f=cp_air(temp_t);%j/(kg*k)
lamta_f=lamta_air(temp_t);%w/mk
lamta_m=p.*lamta_f+(1-p).*lamta_s.*1e3;%w/mk
miu_f=air_miu(temp_t);%m^2/s
a_m=lamta_m./(p.*rou_f.*cp_f+(1-p).*rou_s.*cp_s);%含水织物的有效扩散系数 %m^2/s
%第一层
if(b==1)&(t(b,c)<t_a0)
    temp_derta_t=t_a0-t(b,c);%第一层织物与内层空气的温差
    a=1.26.*(2.*rou_f.*g.*temp_derta_t.*lamta_f^3/(273+temp_t)./a_m./miu_f)^0.5;%内边界上空气对流换热系数，w/m^2*C（4-39）
elseif(b==1)&(t(b,c)>t_a0)
    temp_derta_t=t(b,c)-t_a0;%第一层织物与内层空气的温差
    a=1.26.*(2.*rou_f.*g.*temp_derta_t.*lamta_f^3./(273+temp_t)./a_m./miu_f)^0.5;%内边界上空气对流换热系数，w/m^2*C（4-39）
elseif(b==n)&(t(b,c)<t_an)
    temp_derta_t=t_an-t(b,c);%第n层织物与内层空气的温差
    a=1.26.*(2.*rou_f.*g.*temp_derta_t.*lamta_f^3./(273+temp_t)./a_m./miu_f)^0.5;%内边界上空气对流换热系数，w/m^2*C（4-39）
elseif(b==n)&(t(b,c)>t_an)
    temp_derta_t=t(b,c)-t_an;%第一层织物与内层空气的温差
    a=1.26.*(2.*rou_f.*g.*temp_derta_t.*lamta_f^3./(273+temp_t)./a_m./miu_f)^0.5;%内边界上空气对流换热系数，w/m^2*C（4-39）
else 
    a=0;%即2-9层的对流换热系数为０
end