第三十九篇：换热器设计

摘要

换热器是工业生产中实现热量交换的核心设备，广泛应用于电力、化工、制冷、暖通等领域。本文系统分析了换热器设计的基本理论，推导了对数平均温差法和效能-传热单元数法。详细讨论了管壳式、板式、翅片管等常见换热器类型的结构特点和设计方法。采用热平衡方程和传热系数关联式建立设计模型，计算了换热器的传热面积和压降。通过Python仿真，实现了换热器的热力设计计算，分析了流动布置、传热系数、污垢热阻对换热器性能的影响，为工程换热器设计提供理论指导。

关键词

换热器，对数平均温差，效能-传热单元数，管壳式换热器，板式换热器，翅片管，污垢热阻

---

1. 引言

1.1 换热器的分类

• 间壁式：管壳式、板式、翅片管式
• 蓄热式：回转式、阀门切换式
• 混合式：冷却塔、喷射式

1.2 设计考虑因素

• 传热效率
• 压降损失
• 材料成本
• 维护清洁

---

2. 理论分析

2.1 基本传热方程

$$Q=UA\Delta T_m$$

其中：

• $U$：总传热系数
• $A$：传热面积
• $\Delta T_m$：平均温差

2.2 总传热系数

对于平壁：
$$\frac{1}{U}=\frac{1}{h_i}+\frac{\delta }{k}+\frac{1}{h_o}+R_f$$

对于圆管：
$$\frac{1}{U_o}=\frac{r_o}{r_i{h}_i}+\frac{r_o\ln (r_o/r_i)}{k}+\frac{1}{h_o}+R_{f,o}$$

2.3 对数平均温差

逆流：
$$\Delta T_{lm}=\frac{(T_{h,i}-T_{c,o})-(T_{h,o}-T_{c,i})}{\ln \frac{T_{h,i}-T_{c,o}}{T_{h,o}-T_{c,i}}}$$

并流：
$$\Delta T_{lm}=\frac{(T_{h,i}-T_{c,i})-(T_{h,o}-T_{c,o})}{\ln \frac{T_{h,i}-T_{c,i}}{T_{h,o}-T_{c,o}}}$$

2.4 效能-传热单元数法

效能：
$$\epsilon =\frac{Q_{actual}}{Q_{max}}$$

传热单元数：
$$NTU=\frac{UA}{C_{min}}$$

热容比：
$$C_r=\frac{C_{min}}{C_{max}}$$

1-2型管壳式换热器：
$$\epsilon =2{\left[1+C_r+\sqrt{1+C_r^2}\frac{1+\exp (-NTU\sqrt{1+C_r^2})}{1-\exp (-NTU\sqrt{1+C_r^2})}\right]}^{-1}$$

2.5 压降计算

管程压降：
$$\Delta P_t=4f\frac{L{N}_p}{D_i}\frac{\rho u_m^2}{2}$$

壳程压降（Kern方法）：
$$\Delta P_s=\frac{f{G}_s^2{D}_s(N_b+1)}{2\rho D_e}$$

---

3. Python仿真实现

import numpy as np
import matplotlib.pyplot as plt
import os

output_dir = r'd:\文档\500仿真领域\工程仿真\传热学仿真\主题039'
os.makedirs(output_dir, exist_ok=True)

print("="*60)
print("仿真1：对数平均温差计算")
print("="*60)

# 换热器参数
T_hot_in = 150  # 热流体入口温度，°C
T_hot_out = 90  # 热流体出口温度，°C
T_cold_in = 30  # 冷流体入口温度，°C
T_cold_out = 70  # 冷流体出口温度，°C

# 逆流对数平均温差
delta_T1 = T_hot_in - T_cold_out
delta_T2 = T_hot_out - T_cold_in
LMTD_counter = (delta_T1 - delta_T2) / np.log(delta_T1 / delta_T2)

print(f"逆流对数平均温差: {LMTD_counter:.2f} °C")

# 并流对数平均温差
delta_T1_parallel = T_hot_in - T_cold_in
delta_T2_parallel = T_hot_out - T_cold_out
LMTD_parallel = (delta_T1_parallel - delta_T2_parallel) / np.log(delta_T1_parallel / delta_T2_parallel)

print(f"并流对数平均温差: {LMTD_parallel:.2f} °C")

print(f"逆流/并流温差比: {LMTD_counter/LMTD_parallel:.3f}")

print("\n" + "="*60)
print("仿真2：效能-NTU关系")
print("="*60)

# NTU范围
NTU = np.linspace(0.1, 5, 100)

# 不同流动布置的效能公式
def effectiveness_counterflow(NTU, Cr):
    """逆流换热器效能"""
    if Cr < 1:
        return (1 - np.exp(-NTU * (1 - Cr))) / (1 - Cr * np.exp(-NTU * (1 - Cr)))
    else:
        return NTU / (1 + NTU)

def effectiveness_parallelflow(NTU, Cr):
    """并流换热器效能"""
    return (1 - np.exp(-NTU * (1 + Cr))) / (1 + Cr)

def effectiveness_1shell_2tube(NTU, Cr):
    """1-2型管壳式换热器效能"""
    P1 = 2 / (1 + Cr + np.sqrt(1 + Cr**2) * 
              (1 + np.exp(-NTU * np.sqrt(1 + Cr**2))) / 
              (1 - np.exp(-NTU * np.sqrt(1 + Cr**2))))
    return P1

# 不同热容比
Cr_values = [0.2, 0.5, 0.8, 1.0]

fig, axes = plt.subplots(1, 2, figsize=(14, 5))

ax1 = axes[0]
colors = plt.cm.viridis(np.linspace(0, 1, len(Cr_values)))

for Cr, color in zip(Cr_values, colors):
    eff = [effectiveness_counterflow(n, Cr) for n in NTU]
    ax1.plot(NTU, eff, color=color, linewidth=2, label=f'C_r = {Cr}')

ax1.set_xlabel('NTU', fontsize=11)
ax1.set_ylabel('Effectiveness ε', fontsize=11)
ax1.set_title('Effectiveness-NTU (Counterflow)', fontsize=12, fontweight='bold')
ax1.legend(fontsize=10)
ax1.grid(True, alpha=0.3)
ax1.set_xlim([0, 5])
ax1.set_ylim([0, 1])

# 不同流动布置对比
ax2 = axes[1]
Cr = 0.5
eff_counter = [effectiveness_counterflow(n, Cr) for n in NTU]
eff_parallel = [effectiveness_parallelflow(n, Cr) for n in NTU]
eff_1s2t = [effectiveness_1shell_2tube(n, Cr) for n in NTU]

ax2.plot(NTU, eff_counter, 'b-', linewidth=2, label='Counterflow')
ax2.plot(NTU, eff_parallel, 'r--', linewidth=2, label='Parallel flow')
ax2.plot(NTU, eff_1s2t, 'g-.', linewidth=2, label='1-2 Shell & Tube')
ax2.set_xlabel('NTU', fontsize=11)
ax2.set_ylabel('Effectiveness ε', fontsize=11)
ax2.set_title(f'Effectiveness Comparison (C_r = {Cr})', fontsize=12, fontweight='bold')
ax2.legend(fontsize=10)
ax2.grid(True, alpha=0.3)
ax2.set_xlim([0, 5])
ax2.set_ylim([0, 1])

plt.tight_layout()
plt.savefig(f'{output_dir}/effectiveness_ntu.png', dpi=150, bbox_inches='tight')
plt.close()

print("图1：效能-NTU关系已保存")

print("\n" + "="*60)
print("仿真3：管壳式换热器设计")
print("="*60)

# 设计参数
Q_required = 500e3  # W，换热量

# 流体参数（热流体：水）
m_dot_hot = 2.0  # kg/s
cp_hot = 4186  # J/(kg·K)

# 流体参数（冷流体：水）
m_dot_cold = 3.0  # kg/s
cp_cold = 4186  # J/(kg·K)

# 热容流率
C_hot = m_dot_hot * cp_hot
C_cold = m_dot_cold * cp_cold
C_min = min(C_hot, C_cold)
C_max = max(C_hot, C_cold)
Cr = C_min / C_max

print(f"热容流率比 C_r = {Cr:.3f}")

# 假设传热系数
U_assumed = 800  # W/(m²·K)，水-水换热器

# 计算出口温度（假设逆流）
# 迭代求解
T_hot_out_calc = T_hot_in - Q_required / C_hot
T_cold_out_calc = T_cold_in + Q_required / C_cold

print(f"热流体出口温度: {T_hot_out_calc:.1f} °C")
print(f"冷流体出口温度: {T_cold_out_calc:.1f} °C")

# 对数平均温差
delta_T1 = T_hot_in - T_cold_out_calc
delta_T2 = T_hot_out_calc - T_cold_in
LMTD = (delta_T1 - delta_T2) / np.log(delta_T1 / delta_T2)

# 传热面积
A_required = Q_required / (U_assumed * LMTD)
print(f"所需传热面积: {A_required:.2f} m²")

# 使用NTU法验证
Q_max = C_min * (T_hot_in - T_cold_in)
epsilon = Q_required / Q_max
print(f"换热器效能: {epsilon:.3f}")

# 计算NTU（逆流）
if Cr < 1:
    NTU_calc = np.log((1 - epsilon * Cr) / (1 - epsilon)) / (1 - Cr)
else:
    NTU_calc = epsilon / (1 - epsilon)

A_ntu = NTU_calc * C_min / U_assumed
print(f"NTU法计算面积: {A_ntu:.2f} m²")

print("\n" + "="*60)
print("仿真4：污垢热阻影响")
print("="*60)

# 清洁和污垢条件下的总传热系数
h_i = 5000  # 管内换热系数
h_o = 3000  # 管外换热系数
r_o = 0.0125  # 管外半径，m
r_i = 0.0100  # 管内半径，m
k_wall = 50  # 管壁导热系数

# 清洁状态
R_f_clean = 0
U_clean = 1 / (r_o/(r_i*h_i) + r_o*np.log(r_o/r_i)/k_wall + 1/h_o + R_f_clean)

# 不同污垢热阻
R_f_values = np.linspace(0, 0.0002, 50)  # m²·K/W
U_values = []

for R_f in R_f_values:
    U = 1 / (r_o/(r_i*h_i) + r_o*np.log(r_o/r_i)/k_wall + 1/h_o + R_f)
    U_values.append(U)

fig, ax = plt.subplots(figsize=(10, 5))
ax.plot(R_f_values*1e4, U_values, 'b-', linewidth=2)
ax.axhline(y=U_clean, color='r', linestyle='--', alpha=0.7, label='Clean condition')
ax.set_xlabel('Fouling Resistance R_f (×10⁻⁴ m²·K/W)', fontsize=11)
ax.set_ylabel('Overall Heat Transfer U (W/(m²·K))', fontsize=11)
ax.set_title('Effect of Fouling on Heat Transfer Coefficient', fontsize=12, fontweight='bold')
ax.legend(fontsize=10)
ax.grid(True, alpha=0.3)

plt.tight_layout()
plt.savefig(f'{output_dir}/fouling_effect.png', dpi=150, bbox_inches='tight')
plt.close()

print("图2：污垢热阻影响已保存")

print("\n" + "="*60)
print("仿真5：换热器优化")
print("="*60)

# 分析不同设计参数的影响
# 管径选择
D_tube_range = np.linspace(0.01, 0.05, 50)  # m

# 假设流速恒定
u_tube = 1.5  # m/s

# 计算Reynolds数和换热系数
rho = 1000  # kg/m³
mu = 1e-3  # Pa·s
k_water = 0.6  # W/(m·K)
Pr = 7  # 水的Prandtl数

Re_range = rho * u_tube * D_tube_range / mu

# Dittus-Boelter方程
Nu_range = 0.023 * Re_range**0.8 * Pr**0.4
h_range = Nu_range * k_water / D_tube_range

# 单位长度传热面积
A_per_length = np.pi * D_tube_range

# 单位长度换热量（假设温差）
delta_T_assumed = 50  # K
Q_per_length = h_range * A_per_length * delta_T_assumed

fig, ax1 = plt.subplots(figsize=(10, 5))

color = 'tab:blue'
ax1.set_xlabel('Tube Diameter (mm)', fontsize=11)
ax1.set_ylabel('Heat Transfer Coefficient h (W/(m²·K))', color=color, fontsize=11)
ax1.semilogy(D_tube_range*1000, h_range, color=color, linewidth=2)
ax1.tick_params(axis='y', labelcolor=color)
ax1.grid(True, alpha=0.3)

ax2 = ax1.twinx()
color = 'tab:red'
ax2.set_ylabel('Heat Transfer per Unit Length (W/m)', color=color, fontsize=11)
ax2.plot(D_tube_range*1000, Q_per_length/1000, color=color, linewidth=2)
ax2.tick_params(axis='y', labelcolor=color)

plt.title('Tube Diameter Optimization', fontsize=12, fontweight='bold')
plt.tight_layout()
plt.savefig(f'{output_dir}/tube_optimization.png', dpi=150, bbox_inches='tight')
plt.close()

print("图3：管径优化已保存")

print("\n" + "="*60)
print("仿真6：板式换热器设计计算")
print("="*60)

# 板式换热器参数
A_plate = 0.5  # m²，单板面积
N_plates = 50  # 板片数
U_overall = 3000  # W/(m²·K)，总传热系数

# 流体参数
m_dot_hot_plate = 2.0  # kg/s，热流体流量
m_dot_cold_plate = 2.5  # kg/s，冷流体流量
cp_hot_plate = 4200  # J/(kg·K)
cp_cold_plate = 4186  # J/(kg·K)

T_hot_in_plate = 80  # °C
T_cold_in_plate = 20  # °C

# 热容流率
C_hot_plate = m_dot_hot_plate * cp_hot_plate
C_cold_plate = m_dot_cold_plate * cp_cold_plate
C_min_plate = min(C_hot_plate, C_cold_plate)
C_max_plate = max(C_hot_plate, C_cold_plate)
C_r_plate = C_min_plate / C_max_plate

# 总传热面积
A_total_plate = A_plate * (N_plates - 2)  # 减去两端板

# NTU
NTU_plate = U_overall * A_total_plate / C_min_plate

# 效能（逆流）
epsilon_plate = (1 - np.exp(-NTU_plate * (1 - C_r_plate))) / (1 - C_r_plate * np.exp(-NTU_plate * (1 - C_r_plate)))

# 换热量
Q_max_plate = C_min_plate * (T_hot_in_plate - T_cold_in_plate)
Q_actual_plate = epsilon_plate * Q_max_plate

# 出口温度
T_hot_out_plate = T_hot_in_plate - Q_actual_plate / C_hot_plate
T_cold_out_plate = T_cold_in_plate + Q_actual_plate / C_cold_plate

print(f"热容流率比 C_r: {C_r_plate:.3f}")
print(f"传热单元数 NTU: {NTU_plate:.2f}")
print(f"换热器效能 ε: {epsilon_plate:.3f}")
print(f"最大换热量: {Q_max_plate/1000:.1f} kW")
print(f"实际换热量: {Q_actual_plate/1000:.1f} kW")
print(f"热流体出口温度: {T_hot_out_plate:.1f}°C")
print(f"冷流体出口温度: {T_cold_out_plate:.1f}°C")

# 可视化温度分布
N_channels = N_plates - 1
T_hot_profile = np.linspace(T_hot_in_plate, T_hot_out_plate, N_channels//2)
T_cold_profile = np.linspace(T_cold_out_plate, T_cold_in_plate, N_channels//2)

fig, ax = plt.subplots(figsize=(12, 6))
x_pos = np.arange(len(T_hot_profile))
ax.plot(x_pos, T_hot_profile, 'r-', linewidth=2, marker='o', label='Hot fluid')
ax.plot(x_pos, T_cold_profile, 'b-', linewidth=2, marker='s', label='Cold fluid')
ax.fill_between(x_pos, T_hot_profile, T_cold_profile, alpha=0.2, color='gray', label='Temperature difference')
ax.set_xlabel('Channel Number', fontsize=11)
ax.set_ylabel('Temperature (°C)', fontsize=11)
ax.set_title('Plate Heat Exchanger Temperature Distribution', fontsize=12, fontweight='bold')
ax.legend(fontsize=10)
ax.grid(True, alpha=0.3)

plt.tight_layout()
plt.savefig(f'{output_dir}/plate_hx_temperature.png', dpi=150, bbox_inches='tight')
plt.close()

print("图4：板式换热器温度分布已保存")

print("\n" + "="*60)
print("仿真7：换热器网络优化")
print("="*60)

# 简单换热器网络分析
# 两个换热器串联

# 换热器1参数
U1 = 500  # W/(m²·K)
A1 = 10  # m²
C_min1 = 5000  # W/K

# 换热器2参数
U2 = 400  # W/(m²·K)
A2 = 15  # m²
C_min2 = 6000  # W/K

# 计算单个换热器效能
NTU1 = U1 * A1 / C_min1
NTU2 = U2 * A2 / C_min2

epsilon1 = 1 - np.exp(-NTU1)  # 简化，C_r=0
epsilon2 = 1 - np.exp(-NTU2)

# 串联总效能（简化计算）
epsilon_total = epsilon1 + epsilon2 * (1 - epsilon1)

print(f"换热器1: NTU={NTU1:.2f}, ε={epsilon1:.3f}")
print(f"换热器2: NTU={NTU2:.2f}, ε={epsilon2:.3f}")
print(f"串联总效能: {epsilon_total:.3f}")

# 并联分析
# 流量分配
m_dot_total = 2.0  # kg/s
split_ratio = 0.6  # 60%流量走换热器1

m_dot_1 = m_dot_total * split_ratio
m_dot_2 = m_dot_total * (1 - split_ratio)

print(f"\n并联流量分配:")
print(f"换热器1流量: {m_dot_1:.2f} kg/s")
print(f"换热器2流量: {m_dot_2:.2f} kg/s")

# 可视化网络配置
fig, axes = plt.subplots(1, 2, figsize=(14, 5))

# 串联配置
ax1 = axes[0]
ax1.text(0.2, 0.5, 'HEX1', fontsize=14, ha='center', va='center', 
         bbox=dict(boxstyle='round', facecolor='lightblue'))
ax1.annotate('', xy=(0.4, 0.5), xytext=(0.3, 0.5),
            arrowprops=dict(arrowstyle='->', lw=2))
ax1.text(0.6, 0.5, 'HEX2', fontsize=14, ha='center', va='center',
         bbox=dict(boxstyle='round', facecolor='lightgreen'))
ax1.set_xlim(0, 1)
ax1.set_ylim(0, 1)
ax1.set_title('Series Configuration', fontsize=12, fontweight='bold')
ax1.axis('off')

# 并联配置
ax2 = axes[1]
ax2.text(0.3, 0.7, 'HEX1', fontsize=14, ha='center', va='center',
         bbox=dict(boxstyle='round', facecolor='lightblue'))
ax2.text(0.3, 0.3, 'HEX2', fontsize=14, ha='center', va='center',
         bbox=dict(boxstyle='round', facecolor='lightgreen'))
ax2.annotate('', xy=(0.2, 0.7), xytext=(0.1, 0.5),
            arrowprops=dict(arrowstyle='->', lw=2))
ax2.annotate('', xy=(0.2, 0.3), xytext=(0.1, 0.5),
            arrowprops=dict(arrowstyle='->', lw=2))
ax2.annotate('', xy=(0.5, 0.5), xytext=(0.4, 0.7),
            arrowprops=dict(arrowstyle='->', lw=2))
ax2.annotate('', xy=(0.5, 0.5), xytext=(0.4, 0.3),
            arrowprops=dict(arrowstyle='->', lw=2))
ax2.set_xlim(0, 1)
ax2.set_ylim(0, 1)
ax2.set_title('Parallel Configuration', fontsize=12, fontweight='bold')
ax2.axis('off')

plt.tight_layout()
plt.savefig(f'{output_dir}/hx_network_config.png', dpi=150, bbox_inches='tight')
plt.close()

print("图5：换热器网络配置已保存")
print("\n所有仿真完成！")

---

4. 工程应用

4.1 管壳式换热器

标准TEMA设计：

• TEMA标准分类
• 前封头、壳体、后封头类型
• 设计压力与温度

挡板布置优化：

• 单弓形挡板
• 双弓形挡板
• 螺旋挡板
• 压降与传热平衡

管程数选择：

• 1-2-4-6-8程
• 流速要求
• 清洗维护

4.2 板式换热器

板片波纹设计：

• 人字形波纹
• 水平平直波纹
• 球形波纹
• 传热与压降特性

流道布置：

• 并流、逆流、错流
• 多流程设计
• 温度程序

密封材料选择：

• 橡胶垫片
• 温度限制
• 化学兼容性

4.3 翅片管换热器

翅片类型选择：

• 平直翅片
• 波纹翅片
• 开缝翅片
• 百叶窗翅片

翅片效率计算：
$${\eta }_f=\frac{\tanh (mL)}{mL}$$

其中$m=\sqrt{\frac{2h}{k\delta }}$

结霜除霜设计：

• 结霜对传热的影响
• 除霜方法：热气、电加热、水冲
• 除霜周期优化

---

5. 换热器设计方法

5.1 对数平均温差法（LMTD）

对数平均温差：
$$\Delta T_{lm}=\frac{\Delta T_1-\Delta T_2}{\ln (\Delta T_1/\Delta T_2)}$$

换热量：
$$Q=UA\Delta T_{lm}$$

适用条件：

• 已知进出口温度
• 设计计算
• 校核计算

5.2 效能-NTU法

效能定义：
$$\epsilon =\frac{Q_{actual}}{Q_{max}}$$

NTU定义：
$$NTU=\frac{UA}{C_{min}}$$

关联式：

• 逆流：$\epsilon =\frac{1-e^{-NTU(1-C_r)}}{1-C_r{e}^{-NTU(1-C_r)}}$
• 并流：$\epsilon =\frac{1-e^{-NTU(1+C_r)}}{1+C_r}$

适用条件：

• 已知效能或出口温度
• 优化设计
• 变工况分析

5.3 污垢热阻

常见污垢热阻值：

• 冷却水：0.0001-0.0002 m²·K/W
• 海水：0.0001-0.0003 m²·K/W
• 油类：0.0002-0.001 m²·K/W

设计考虑：

• 预留污垢裕量
• 清洗周期
• 水质处理

---

6. 本章小结

换热器设计需要综合考虑传热效率、压降损失和经济性，对数平均温差法和效能-NTU法是两种主要的设计方法，污垢热阻是长期运行的重要考虑因素。

核心内容

1. 设计方法：

   • LMTD法
   • 效能-NTU法
   • 污垢热阻

2. 换热器类型：

   • 管壳式
   • 板式
   • 翅片管式

3. 优化设计：

   • 传热与压降平衡
   • 经济性分析
   • 可靠性考虑

工程价值

1. 能源系统：提高能源利用效率
2. 化工过程：反应器热管理
3. HVAC系统：建筑环境控制
4. 动力工程：余热回收

换热器是热力系统的核心设备，掌握其设计方法对于工程实践具有重要意义。