对流换热实验

(1) 空气横掠单管强制对流换热实验

(1.1) 代码

clc;clear;
d = 3.73e-3;
l = 99.6e-3;
tw = [40.57 42.41 44.27 46.48 49.72];
tf = [21.12 21.59 21.77 21.77 21.99];
deltaP = [213.5 180.3 144.0 108.5 74.4];
U = [0.3946 0.3979 0.3985 0.3992 0.4004];
I = [14.601 14.689 14.688 14.687 14.685];
fprintf('1. 定性温度（℃）：');
tm = (tw + tf)/2
fprintf('1.1 以定性温度查的密度（kg/m^3）：');
rou = (1.165-1.128)/(30-40)*(tm - 30)+1.165
fprintf('1.2 以定性温度查的导热系数（W/(m·K)）：');
lamda = (2.76e-2-2.67e-2)/(40-30)*(tm - 30)+2.67e-2
fprintf('1.3 以定性温度查的运动粘度（m^2/s）：');
miu = ((16.96-16.00)/(40-30)*(tm - 30)+16.00).*1e-6
fprintf('1.4 以定性温度查的普朗特数：');
Pr = (0.699-0.701)/(40-30)*(tm - 30)+0.701
fprintf('2. 流体速度（m/s）：');
u = sqrt(2*deltaP./rou)
fprintf('3. 雷诺数：');
Re = u.*d./miu
fprintf('4. 实验测得努塞尔数：');
h = (U.*I)./(pi*d*l)./(tw-tf);
Nu = h.*d./lamda
Nu_cal = zeros(1,5);
for i = 1:5
    if (Re(i)>=40)&&(Re(i)<=4000)
        Nu_cal(i) = 0.683*(Re(i)^0.466)*(Pr(i)^(1/3));
    elseif (Re(i)>4000)&&(Re(i)<=40000)
        Nu_cal(i) = 0.193*(Re(i)^0.618)*(Pr(i)^(1/3));
    end
end
fprintf('5. 努塞尔数理论计算值：');
Nu_cal
fprintf('6. 相对偏差（%%）：');
e = 100*abs(Nu - Nu_cal)./Nu_cal
fprintf('7. 对流换热量计算（单位J）：');
qConvection = U.*I
fprintf('8. 辐射换热量估算（取表面发射率为0.5，单位J）：');
qRadiation = 0.5*5.67e-8*(tw.^4-tf.^4)*(pi*d*l)

(1.2) 运行结果

1. 定性温度（℃）：
tm =
   30.8450   32.0000   33.0200   34.1250   35.8550
1.1 以定性温度查的密度（kg/m^3）：
rou =
    1.1619    1.1576    1.1538    1.1497    1.1433
1.2 以定性温度查的导热系数（W/(m·K)）：
lamda =
    0.0268    0.0269    0.0270    0.0271    0.0272
1.3 以定性温度查的运动粘度（m^2/s）：
miu =
   1.0e-04 *
    0.1608    0.1619    0.1629    0.1640    0.1656
1.4 以定性温度查的普朗特数：
Pr =
    0.7008    0.7006    0.7004    0.7002    0.6998
2. 流体速度（m/s）：
u =
   19.1705   17.6495   15.7989   13.7382   11.4081
3. 雷诺数：
Re =
   1.0e+03 *
    4.4466    4.0658    3.6176    3.1254    2.5693
4. 实验测得努塞尔数：
Nu =
   35.3561   33.3770   30.8240   28.0113   24.8891
5. 努塞尔数理论计算值：
Nu_cal =
   30.8018   29.1405   27.6115   25.7897   23.5354
6. 相对偏差（%）：
e =
   14.7859   14.5381   11.6346    8.6141    5.7518
7. 对流换热量计算（单位J）：
qConvection =
    5.7616    5.8448    5.8532    5.8631    5.8799
8. 辐射换热量估算（取表面发射率为0.5，单位J）：
qRadiation =
   1.0e-03 *
0.0831    0.0999    0.1197    0.1470    0.1945

(2) 水平管外自然对流换热实验

(2.1) 代码

%用最小二乘拟合计算出实验关联式
clc;clear;
syms LnGrPr;
GrPr = [458934.6 128178.6 81973.53 47893.61 1903764 1523845 241316.5 586097.9 811834.9 1291265];
Nu = [12.77 9.53 9.42 8.02 17.51 17.01 11.39 15.16 13.1013 15.54];
ln_GrPr = log(GrPr);
ln_Nu = log(Nu);
p = polyfit(ln_GrPr,ln_Nu,1); %1次多项式拟合（即线性拟合）
LnNu = vpa(poly2sym(p,LnGrPr),5)
n = vpa(p(1),5)
C = vpa(exp(p(2)),5) %求n和C

%双对数坐标图绘制曲线
clc;clear;
GrPr = [458934.6 128178.6 81973.53 47893.61 1903764 1523845 241316.5 586097.9 811834.9 1291265];
Nu = [12.77 9.53 9.42 8.02 17.51 17.01 11.39 15.16 13.1013 15.54];
figure;
loglog(GrPr,Nu,'*');
hold on;
grid on;
p = polyfit(log(GrPr),log(Nu),1);
x = min(GrPr):0.1:max(GrPr);
y = exp(polyval(p,log(x)));
loglog(x,y,'-');
%legend('原始数据点','拟合曲线','location','northwest');
xlabel('GrPr');
ylabel('Nu');
title('水平管外自然对流换热实验曲线');
scatter(GrPr(9),Nu(9),'o');
legend('原始数据点','拟合曲线','3#组数据点','location','northwest');
gtext(['3#','(',num2str(GrPr(9)),',',num2str(Nu(9)),')']);

(2.2) 运行结果

LnNu =
0.20561*LnGrPr - 0.12204 
n = 
0.20561
C =
0.88511

---

$$参考资料来源$$

1. 《传热学(第五版)》.陶文铨.杨世铭.高等教育出版社.
2. 《热与流体实验教程》.王小丹.孟倩.张可.吴青平.唐上朝.西安交通大学出版社.

---

$$源码作者：AidenLee$$

$$博客创作：AidenLee$$

特别声明：文章仅供学习参考，转载请注明出处，严禁盗用！