传热 -- 化工原理实验Python代码

实验代码针对天津大学、重庆理工大学化工原理实验教材设计，若有其他需求，请自行编辑

无强化丝原始数据

最小压差 (Δp_min) = 0.77 kPa, 最大压差 (Δp_max) = 4.35 kPa

序号	Δp (kPa)	t_入 (°C)	t_出 (°C)
1	0.77	33.6	64.8
2	1.37	32.7	65.0
3	1.98	33.9	64.6
4	2.58	35.7	64.8
5	3.18	37.8	65.8
6	3.77	40.5	65.8

用Mathmatica处理后数据

序号	1	2	3	4	5	6
Δp	0.77	1.37	1.98	2.58	3.18	3.77
t_入	33.6	32.7	33.9	35.7	37.8	40.5
t_出	64.8	65.0	64.6	64.8	65.8	65.8
t_平	49.2	48.85	49.25	50.25	51.8	53.15
ρ	1.09441	1.09564	1.09423	1.09073	1.08533	1.08067
λ	0.0282876	0.0282603	0.0282915	0.0283695	0.0284903	0.0285955
μ	0.00001958	0.00001956	0.00001958	0.00001963	0.00001971	0.00001977
Pr	0.695658	0.695733	0.695648	0.695435	0.695108	0.694824
Pr^0.4	0.864885	0.864922	0.86488	0.864774	0.864611	0.86447
V_t	19.9257	26.5634	31.9547	36.535	40.662	44.369
V_修	20.939	27.966	33.5522	38.2561	42.4928	46.1585
u_m	18.5141	24.7274	29.6666	33.8259	37.5718	40.8131
W_c	0.00636549	0.00851129	0.0101983	0.0115908	0.0128108	0.0138562
Q	199.596	276.289	314.653	338.979	360.495	352.315
α_i	53.8055	75.0133	84.7353	89.4697	92.3014	87.9157
Nu_i	38.0418	53.0875	59.9016	63.0746	64.7949	61.4892
Re_i	20696	27696.4	33153.4	37588.2	41387.7	44618.8
Nu/Pr^0.4	43.9848	61.3783	69.26	72.9376	74.9411	71.1293
Δt_1	64.8	65.7	64.5	62.7	60.6	57.9
Δt_2	33.6	33.4	33.8	33.6	32.6	32.6
Δt_m	47.5045	47.7427	47.5082	46.6469	45.1626	44.0456
K_o	50.6598	69.7755	79.8563	87.6186	96.2425	96.4439

无强化丝双对数拟合

import numpy as np
import matplotlib.pyplot as plt

# 原始数据（这里没有用mathmatica处理的数据，用的手算的数据，因个人选取参数不同，结果有略微差异）
re = np.array([20700, 38100, 33078.9, 36800, 40525, 45200])
nu_pr = np.array([44.46, 61.9, 70.165, 74.57, 81.65, 102.99])

# 对数转换
log_nu_pr = np.log10(nu_pr)
log_re = np.log10(re)

# 执行线性拟合
fit = np.polyfit(log_re, log_nu_pr, 1)
fit_fn = np.poly1d(fit)

# 设置字体
plt.rcParams['font.family'] = 'SimHei'
plt.rcParams['axes.unicode_minus'] = False

# 绘制双对数图
plt.figure(figsize=(8, 6))
plt.scatter(log_re, log_nu_pr, color='blue', label='Data')
plt.plot(log_re, fit_fn(log_re), color='red', label='Fit')
plt.xlabel(r'$\log_{10}(Re)$')
plt.ylabel(r'$\log_{10}(\frac{Nu}{Pr^{0.4}})$')
plt.title('无强化丝双对数拟合')
plt.text(4.33,1.8, f'无强化丝拟合:y = {fit[0]:.2f}x + {fit[1]:.2f}', fontsize=12, color='red')
plt.legend()
plt.grid(True)
plt.show()

有强化丝原始数据

最小压差 (Δp_min) = 0.37 kPa, 最大压差 (Δp_max) = 2.35 kPa

序号	Δp (kPa)	t_入 (°C)	t_出 (°C)
1	0.38	33.3	71.6
2	0.7	32.5	74.8
3	1.05	33.9	75.4
4	1.36	36.1	75.4
5	1.73	39.2	75.9
6	2.09	43.2	76.5

用Mathmatica处理后数据

序号	1	2	3	4	5	6
Δp	0.38	0.7	1.05	1.36	1.73	2.09
t_入	33.3	32.5	33.9	36.1	39.2	43.2
t_出	71.6	74.8	75.4	75.4	75.9	76.5
t_平	52.45	53.65	54.65	55.75	57.55	59.85
ρ	1.0831	1.079	1.0755	1.0718	1.0657	1.0581
λ	0.02854	0.02863	0.02871	0.0288	0.02894	0.02912
μ	0.00002	0.00002	0.00002	0.00002	0.00002	0.00002
Pr	0.69497	0.69472	0.69451	0.69428	0.69391	0.69344
Pr^0.4	0.86454	0.86442	0.86431	0.8642	0.86401	0.86378
V_t	14.071	19.134	23.471	26.759	30.266	33.386
V_修	14.95	20.458	25.057	28.459	32.044	35.143
u_m	13.219	18.089	22.156	25.163	28.333	31.074
W_c	0.0045	0.0061	0.0075	0.0085	0.0095	0.0103
Q	173.128	260.657	312.231	334.653	349.888	345.682
α_i	43.778	64.438	75.775	79.614	80.635	76.604
Nu_i	30.678	45.007	52.782	55.292	55.73	52.619
Re_i	14508.165	19720.195	24019.392	27113.555	30225.646	32731.759
Nu/Pr^0.4	35.4842	52.0663	61.0683	63.98	64.501	60.917
Δt_1	65.1	65.9	64.5	62.3	59.2	55.2
Δt_2	26.8	23.6	23.0	23.0	22.5	21.9
Δt_m	43.154	41.192	40.246	39.439	37.936	36.02
K_o	48.372	76.296	93.542	102.308	111.203	115.711

有强化丝双对数拟合

import numpy as np
import matplotlib.pyplot as plt

# 原始数据
re = np.array([14500, 20900, 23200, 27280.4, 29498.78, 32700])
nu_pr = np.array([33.31, 52.3, 61.0683, 65.96, 63, 60.87])

# 对数转换
log_re = np.log10(re)
log_nu_pr = np.log10(nu_pr)

# 执行线性拟合
fit = np.polyfit(log_re, log_nu_pr, 1)
fit_fn = np.poly1d(fit)

# 设置字体
plt.rcParams['font.family'] = 'SimHei'
plt.rcParams['axes.unicode_minus'] = False

# 绘制双对数图
plt.figure(figsize=(8, 6))
plt.scatter(log_re, log_nu_pr, color='blue', label='Data')
plt.plot(log_re, fit_fn(log_re), color='red', label='Fit')
plt.xlabel(r'$\log_{10}(Re)$')
plt.ylabel(r'$\log_{10}(\frac{Nu}{Pr^{0.4}})$')
plt.title('有强化丝双对数拟合')
plt.text(4.3,1.8, f'有强化丝拟合:y = {fit[0]:.2f}x + {fit[1]:.2f}', fontsize=12, color='red')

plt.legend()
plt.grid(True)
plt.show()

无强化丝vs.有强化丝

import numpy as np
import matplotlib.pyplot as plt

# 原始数据
re_no_enhancement = np.array([20700, 38100, 33078.9, 36800, 40525, 45200])
nu_pr_no_enhancement = np.array([44.46, 61.9, 70.165, 74.57, 81.65, 102.99])

re_enhancement = np.array([14500, 20900, 23200, 27280.4, 29498.78, 32700])
nu_pr_enhancement = np.array([33.31, 52.3, 61.0683, 65.96, 63, 60.87])

# 对数转换
log_re_no_enhancement = np.log10(re_no_enhancement)
log_nu_pr_no_enhancement = np.log10(nu_pr_no_enhancement)

log_re_enhancement = np.log10(re_enhancement)
log_nu_pr_enhancement = np.log10(nu_pr_enhancement)

# 执行线性拟合
fit_no_enhancement = np.polyfit(log_re_no_enhancement, log_nu_pr_no_enhancement, 1)
fit_fn_no_enhancement = np.poly1d(fit_no_enhancement)

fit_enhancement = np.polyfit(log_re_enhancement, log_nu_pr_enhancement, 1)
fit_fn_enhancement = np.poly1d(fit_enhancement)

# 设置字体
plt.rcParams['font.family'] = 'SimHei'
plt.rcParams['axes.unicode_minus'] = False

# 绘制双对数图
plt.figure(figsize=(8, 6))

# 无强化丝数据及拟合直线
plt.scatter(log_re_no_enhancement, log_nu_pr_no_enhancement, color='blue', label='无强化丝数据')
plt.plot(log_re_no_enhancement, fit_fn_no_enhancement(log_re_no_enhancement), color='red', label='无强化丝拟合线')

# 有强化丝数据及拟合直线
plt.scatter(log_re_enhancement, log_nu_pr_enhancement, color='green', label='有强化丝数据')
plt.plot(log_re_enhancement, fit_fn_enhancement(log_re_enhancement), color='orange', label='有强化丝拟合线')

plt.xlabel(r'$\log_{10}(Re)$')
plt.ylabel(r'$\log_{10}(\frac{Nu}{Pr^{0.4}})$')
plt.title('有强化丝vs.无强化丝')

# 添加拟合方程
plt.text(4.25, 1.80, f'无强化丝拟合方程: y = {fit_no_enhancement[0]:.2f}x + {fit_no_enhancement[1]:.2f}', fontsize=10, color='red')
plt.text(4.10, 1.70, f'有强化丝拟合方程: y = {fit_enhancement[0]:.2f}x + {fit_enhancement[1]:.2f}', fontsize=10, color='orange')

plt.legend()
plt.grid(True)
plt.show()