换热器综合实验

最新推荐文章于 2024-07-10 22:11:32 发布

wibibaboo

最新推荐文章于 2024-07-10 22:11:32 发布

阅读量854

点赞数

分类专栏：杂货铺文章标签： matlab 编程语言

本文链接：https://blog.csdn.net/dasdkhaksd/article/details/115013218

版权

杂货铺专栏收录该内容

10 篇文章 0 订阅

订阅专栏

(1) 代码

clc;clear;

heatExchangeArea = 0.14; %换热面积（m^2）。
flowArea = (pi / 4) * ((16 - 1.5 * 2) * 0.001)^2; %内管流体横截面面积（m^2）。

%顺流。
parallelFlowTi1 = [60.2 59.4 60.0 60.1 60.0]; %热流体进口温度（℃）。
parallelFlowTo1 = [54.0 52.7 52.4 51.6 50.7]; %热流体出口温度（℃）。
parallelFlowTi2 = [14.5 13.8 13.3 13.3 13.1]; %冷流体进口温度（℃）。
parallelFlowTo2 = [36.4 35.2 34.6 34.0 33.1]; %冷流体出口温度（℃）。
parallelFlowQ1 = [548.4 503.0 453.4 403.8 352.8]; %热流体流量（L/h）。
parallelFlowQ2 = [160 160 160 160 160]; %冷流体流量（L/h）。
parallelFlowTmax = zeros(1,5);
parallelFlowTmin = zeros(1,5);
parallelFlowTm = zeros(1,5); %对数平均温差（℃）。
parallelFlowDensity1 = [985.38 985.86 985.79 985.95 986.18]; %热流体密度（kg/m^3）。
parallelFlowDensity2 = [997.00 997.24 997.37 997.44 997.57]; %冷流体密度（kg/m^3）。
parallelFlowHeatCapacity1 = [4.1776 4.1770 4.1771 4.1769 4.1767]; %热流体比热容（kJ/(kg*K)）。
parallelFlowHeatCapacity2 = [4.1781 4.1790 4.1794 4.1797 4.1802]; %冷流体比热容（kJ/(kg*K)）。
parallelFlowHeat1 = zeros(1,5); %热流体换热量（kW）。
parallelFlowHeat2 = zeros(1,5); %冷流体换热量（kW）。
parallelFlowHeatErr = zeros(1,5); %热平衡偏差（%）。
parallelFlowK = zeros(1,5); %传热系数。
for i = 1:5
    %求对数平均温差。
    parallelFlowTmax(i) = max(abs(parallelFlowTi1(i) - parallelFlowTi2(i)),abs(parallelFlowTo1(i) - parallelFlowTo2(i)));
    parallelFlowTmin(i) = min(abs(parallelFlowTi1(i) - parallelFlowTi2(i)),abs(parallelFlowTo1(i) - parallelFlowTo2(i)));
    parallelFlowTm(i) = (parallelFlowTmax(i) - parallelFlowTmin(i)) / (log(parallelFlowTmax(i) / parallelFlowTmin(i)));
    %求热流体换热量。
    parallelFlowHeat1(i) = parallelFlowQ1(i) * 0.001/3600 * parallelFlowDensity1(i) * parallelFlowHeatCapacity1(i) * ...
    (parallelFlowTi1(i) - parallelFlowTo1(i));
    %求冷流体换热量。
    parallelFlowHeat2(i) = parallelFlowQ2(i) * 0.001/3600 * parallelFlowDensity2(i) * parallelFlowHeatCapacity2(i) * ...
    (parallelFlowTo2(i) - parallelFlowTi2(i));
    %求热平衡偏差。
    parallelFlowHeatErr(i) = 2 * (parallelFlowHeat1(i) - parallelFlowHeat2(i)) / (parallelFlowHeat1(i) + ...
    parallelFlowHeat2(i)) * 100;
    %求传热系数。
    parallelFlowK(i) = (parallelFlowHeat1(i) + parallelFlowHeat2(i)) / 2 / heatExchangeArea / parallelFlowTm(i) * 1000;
end
disp("顺流下各工况=>");
fprintf("\t1.对数平均温差：(℃)\n");
disp(parallelFlowTm);
fprintf("\t2.热流体换热量：（kW）\n");
disp(parallelFlowHeat1);
fprintf("\t3.冷流体换热量：（kW）\n");
disp(parallelFlowHeat2);
fprintf("\t4.热平衡偏差：（%%）\n");
disp(parallelFlowHeatErr);
fprintf("\t5.传热系数：\n");
disp(parallelFlowK);

%逆流。
counterFlowTi1 = [59.6 59.8 60.3 60.2 60.4]; %热流体进口温度（℃）。
counterFlowTo1 = [53.2 52.8 52.5 51.5 50.7]; %热流体出口温度（℃）。
counterFlowTi2 = [13.0 13.0 13.0 13.0 13.0]; %冷流体进口温度（℃）。
counterFlowTo2 = [35.8 35.3 35.1 34.5 33.8]; %冷流体出口温度（℃）。
counterFlowQ1 = [554.3 503.9 453.5 402.6 353.2]; %热流体流量（L/h）。
counterFlowQ2 = [160 160 160 160 160]; %冷流体流量（L/h）。
counterFlowTmax = zeros(1,5);
counterFlowTmin = zeros(1,5);
counterFlowTm = zeros(1,5); %对数平均温差（℃）。
counterFlowDensity1 = [985.70 985.74 985.70 985.95 986.09]; %热流体密度（kg/m^3）。
counterFlowDensity2 = [997.26 997.32 997.35 997.42 997.50]; %冷流体密度（kg/m^3）。
counterFlowHeatCapacity1 = [4.1772 4.1772 4.1772 4.1769 4.1768]; %热流体比热容（kJ/(kg*K)）。
counterFlowHeatCapacity2 = [4.1790 4.1793 4.1794 4.1796 4.1799]; %冷流体比热容（kJ/(kg*K)）。
counterFlowHeat1 = zeros(1,5); %热流体换热量（kW）。
counterFlowHeat2 = zeros(1,5); %冷流体换热量（kW）。
counterFlowHeatErr = zeros(1,5); %热平衡偏差（%）。
counterFlowK = zeros(1,5); %传热系数。
for i = 1:5
    %求对数平均温差。
    counterFlowTmax(i) = max(abs(counterFlowTi1(i) - counterFlowTo2(i)),abs(counterFlowTo1(i) - counterFlowTi2(i)));
    counterFlowTmin(i) = min(abs(counterFlowTi1(i) - counterFlowTo2(i)),abs(counterFlowTo1(i) - counterFlowTi2(i)));
    counterFlowTm(i) = (counterFlowTmax(i) - counterFlowTmin(i)) / (log(counterFlowTmax(i) / counterFlowTmin(i)));
    %求热流体换热量。
    counterFlowHeat1(i) = counterFlowQ1(i) * 0.001/3600 * counterFlowDensity1(i) * counterFlowHeatCapacity1(i) * ...
    (counterFlowTi1(i) - counterFlowTo1(i));
    %求冷流体换热量。
    counterFlowHeat2(i) = counterFlowQ2(i) * 0.001/3600 * counterFlowDensity2(i) * counterFlowHeatCapacity2(i) * ...
    (counterFlowTo2(i) - counterFlowTi2(i));
    %求热平衡偏差。
    counterFlowHeatErr(i) = 2 * (counterFlowHeat1(i) - counterFlowHeat2(i)) / (counterFlowHeat1(i) + ...
    counterFlowHeat2(i)) * 100;
    %求传热系数。
    counterFlowK(i) = (counterFlowHeat1(i) + counterFlowHeat2(i)) / 2 / heatExchangeArea / counterFlowTm(i) * 1000;
end
disp("逆流下各工况=>");
fprintf("\t1.对数平均温差：(℃)\n");
disp(counterFlowTm);
fprintf("\t2.热流体换热量：（kW）\n");
disp(counterFlowHeat1);
fprintf("\t3.冷流体换热量：（kW）\n");
disp(counterFlowHeat2);
fprintf("\t4.热平衡偏差：（%%）\n");
disp(counterFlowHeatErr);
fprintf("\t5.传热系数：\n");
disp(counterFlowK);

%用最小二乘法拟合数据求实验关联式。
disp("实验关联式拟合=>");
parallelFlowVelocity = parallelFlowQ1 * 0.001 / 3600 / flowArea;
counterFlowVelocity = counterFlowQ1 * 0.001 / 3600 / flowArea;
parallelFlowP = polyfit(log(parallelFlowVelocity),log(parallelFlowK),1); %最小二乘拟合算出顺流的n和lnc。
fprintf("\t1.顺流实验关联式系数为：\n");
fprintf("\t\tc = %6.3f\tn = %6.3f\n",exp(parallelFlowP(2)),parallelFlowP(1));
counterFlowP = polyfit(log(counterFlowVelocity),log(counterFlowK),1); %最小二乘拟合算出逆流的n和lnc。
fprintf("\t2.逆流实验关联式系数为：\n");
fprintf("\t\tc = %6.3f\tn = %6.3f\n",exp(counterFlowP(2)),counterFlowP(1));
%数据可视化。
figure;
hold on;
axis tight;
grid on;
xlabel('u/(m·s^{-1})');
ylabel('k');
title('套管式换热器传热系数与管内流体流速的关系');
scatter(parallelFlowVelocity,parallelFlowK,'rp');
h1 = fplot(@(x)exp(parallelFlowP(2)).*x.^parallelFlowP(1),[0.6,1.2]);
set(h1,'linestyle','-.');
scatter(counterFlowVelocity,counterFlowK,'b*');
h2 = fplot(@(x)exp(counterFlowP(2)).*x.^counterFlowP(1),[0.6,1.2]);
set(h2,'linestyle','--');
legend('顺流工况数据点','顺流工况拟合曲线','逆流工况数据点','逆流工况拟合曲线','location','northwest');
parallelFlowCurveName = strcat('k=',num2str(exp(parallelFlowP(2)),3),'u','^','{',num2str(parallelFlowP(1),3),'}');
counterFlowCurveName = strcat('k=',num2str(exp(counterFlowP(2)),3),'u','^','{',num2str(counterFlowP(1),3),'}');
text(0.73,920,parallelFlowCurveName);
text(0.88,890,counterFlowCurveName);

(2) 运行结果

顺流下各工况=>
	1.对数平均温差：(℃)
   29.4488   29.3409   29.9623   29.8573   29.8943
	2.热流体换热量：（kW）
    3.8879    3.8550    3.9414    3.9264    3.7540
	3.冷流体换热量：（kW）
    4.0545    3.9637    3.9461    3.8355    3.7067
	4.热平衡偏差：（%）
   -4.1943   -2.7821   -0.1187    2.3421    1.2688
	5.传热系数：
  963.2227  951.7040  940.1705  928.4457  891.3228
逆流下各工况=>
	1.对数平均温差：(℃)
   31.2869   31.5338   31.8162   31.6701   31.8281
	2.热流体换热量：（kW）
    4.0574    4.0345    4.0458    4.0068    3.9197
	3.冷流体换热量：（kW）
    4.2231    4.1310    4.0942    3.9835    3.8544
	4.热平衡偏差：（%）
   -4.0015   -2.3653   -1.1909    0.5827    1.6787
	5.传热系数：
  945.2342  924.8049  913.7278  901.0711  872.3324
实验关联式拟合=>
	1.顺流实验关联式系数为：
		c = 945.094	n =  0.166
	2.逆流实验关联式系数为：
		c = 920.978	n =  0.166

在这里插入图片描述

$参考资料来源$

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

$源码作者：Aiden\ Lee$

$博客创作：Aiden\ Lee$

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

wibibaboo

关注

0
点赞
踩
8

收藏

觉得还不错? 一键收藏
1
评论
换热器综合实验

目录(1) 代码(2) 运行结果(1) 代码clc;clear;heatExchangeArea = 0.14; %换热面积（m^2）。flowArea = (pi / 4) * ((16 - 1.5 * 2) * 0.001)^2; %内管流体横截面面积（m^2）。%顺流。parallelFlowTi1 = [60.2 59.4 60.0 60.1 60.0]; %热流体进口温度（℃）。parallelFlowTo1 = [54.0 52.7 52.4 51.6 50.7]; %热流体出
复制链接

扫一扫