mcpField
T
thermodynamics-workflow

thermodynamics-workflow

github.com/soljourner其他9551 次安装51 周前更新
代码仓库

关于 thermodynamics-workflow

面向泵系统的热力学分析工作流，涵盖能量与熵平衡、热传递、效率和循环计算。适合分析泵及其在更大热力循环中作用的工程人员。

主要功能

  1. 01定义控制体与状态点
  2. 02评估流体热力学性质
  3. 03计算能量与熵平衡
  4. 04分析效率、损失与火用
  5. 05开展相关循环计算

适用场景

  1. 01核算泵的轴功和温升
  2. 02评估泵效率及不可逆损失
  3. 03分析朗肯循环中的给水泵

安装

复制下面的命令，在终端中运行。

$npx skills add https://github.com/soljourner/claude-engineering-skills --skill thermodynamics-workflow

安装后在 Claude Code 中输入/thermodynamics-workflow确认已加载。

附带文件equations-reference.mdexamples.md

详细介绍

热力学分析工作流

概述

此技能提供了一套系统的泵系统热力学分析方法，包括能量平衡、熵分析、传热和循环计算。该方法遵循经典热力学原理，并专门应用于泵系统和流体机械。

热力学分析工作流

1. 系统定义与边界

定义控制体：

任何热力学分析的第一步都是明确定义系统边界。

对于泵系统：

  • 控制体： 从吸入法兰到排出法兰的泵壳
  • 边界： 流体／能量穿过的物理表面
  • 接口：
    • 入口：吸入法兰（测点 1）
    • 出口：排出法兰（测点 2）
    • 轴：机械能输入
    • 壁面：传热表面（如适用）

系统分类：

  • 开放系统（控制体）： 大多数泵分析
  • 稳态流动： 属性不随时间变化
  • 瞬态流动： 启动、停机、喘振工况

需记录的假设：

  • 稳态运行
  • 入口／出口的属性均匀
  • 动能变化可忽略（如适用）
  • 壁面绝热（或说明传热情况）
  • 无化学反应
  • 不可压缩流动（液体）或可压缩流动（气体）

坐标系：

  • 定义势能基准面（z = 0）
  • 指定流动方向（正向流动）
  • 确定测量平面
2. 状态点识别

热力学状态定义：

对于纯物质，两个相互独立的强度性质即可确定其状态。

状态点 1（入口／吸入）：

识别并测量：

Properties at inlet:
- Pressure: P₁ (Pa, bar, psi)
- Temperature: T₁ (K, °C)
- Velocity: c₁ (m/s)
- Elevation: z₁ (m)
- Phase: Liquid, vapor, two-phase

状态点 2（出口／排出）：

识别并测量：

Properties at outlet:
- Pressure: P₂ (Pa, bar, psi)
- Temperature: T₂ (K, °C)
- Velocity: c₂ (m/s)
- Elevation: z₂ (m)
- Phase: Liquid, vapor, two-phase

其他状态点：

对于复杂系统：

  • 级间工况（多级泵）
  • 叶轮出口（蜗壳损失之前）
  • 中间换热器
  • 再循环回路

性质图：

在以下图表上标出状态点：

  • P-h 图： 压力与比焓的关系
  • T-s 图： 温度与比熵的关系
  • P-v 图： 压力与比容的关系
  • 莫里尔图： h-s 图
3. 性质评估

流体性质数据来源：

对于纯物质：

  • 水／蒸汽： IAPWS-IF97 公式（通过 CoolProp）
  • 制冷剂： NIST REFPROP 数据库
  • 烃类： Peng-Robinson 或 Soave-Redlich-Kwong 状态方程
  • 低温流体： NIST 关联式

对于不可压缩液体：

Specific volume: v ≈ constant
Specific heat: cp ≈ constant
Internal energy: du = cp·dT
Enthalpy: dh = cp·dT + v·dP

对于理想气体：

Equation of state: P·v = R·T
Specific heats: cp, cv = f(T)
Internal energy: du = cv·dT
Enthalpy: dh = cp·dT
Entropy: ds = cp·dT/T - R·dP/P

性质评估方法：

方法 1：使用 CoolProp（推荐）

from CoolProp.CoolProp import PropsSI

# Water at 25°C, 2 bar
h = PropsSI('H', 'T', 298.15, 'P', 200000, 'Water')  # J/kg
s = PropsSI('S', 'T', 298.15, 'P', 200000, 'Water')  # J/kg-K
rho = PropsSI('D', 'T', 298.15, 'P', 200000, 'Water')  # kg/m³

方法 2：性质表

  • 蒸汽表（过热蒸汽、饱和蒸汽、压缩液体）
  • 在表值之间进行插值
  • 线性或二次插值

方法 3：经验关联式

  • 密度：ρ(T, P) 关联式
  • 黏度：Andrade 或 Vogel 方程
  • 比热容：多项式拟合

需评估的关键性质：

每个状态点均需评估：

  • 比焓： h (J/kg 或 kJ/kg)
  • 比熵： s (J/kg-K 或 kJ/kg-K)
  • 比容： v (m³/kg) 或密度 ρ (kg/m³)
  • 比内能： u (J/kg 或 kJ/kg)
  • 干度： x（用于两相混合物）
  • 黏度： μ (Pa·s)——用于损失计算
  • 导热系数： k (W/m-K)——用于传热
4. 第一定律分析（能量平衡）

开放系统的一般能量方程：

对于稳态流动的控制体：

Q̇ - Ẇ = ṁ[(h₂ - h₁) + (c₂² - c₁²)/2 + g(z₂ - z₁)]

其中：

  • Q̇ = 热传递速率（W）
  • Ẇ = 功传递速率（W）
  • ṁ = 质量流量（kg/s）
  • h = 比焓（J/kg）
  • c = 速度（m/s）
  • g = 9.81 m/s²
  • z = 高度（m）

泵系统：

按照惯例，功作用于流体（Ẇ < 0），因此：

Ẇshaft = ṁ[(h₂ - h₁) + (c₂² - c₁²)/2 + g(z₂ - z₁)] + Q̇loss

简化形式：

不可压缩液体（密度恒定）：

Ẇshaft/ṁ = (P₂ - P₁)/ρ + (c₂² - c₁²)/2 + g(z₂ - z₁) + q_loss

可用扬程表示为：

H = (P₂ - P₁)/(ρg) + (c₂² - c₁²)/(2g) + (z₂ - z₁) + h_loss

绝热泵（无热损失）：

Ẇshaft = ṁ[(h₂ - h₁) + (c₂² - c₁²)/2 + g(z₂ - z₁)]

动能和势能变化可忽略：

Ẇshaft = ṁ(h₂ - h₁) + Q̇loss

能量平衡的组成部分：

  1. 焓变：

    • 压力升高：Δh_pressure = ∫v·dP
    • 温度升高：Δh_thermal = cp·ΔT

  2. 动能变化：

    • 对泵而言通常较小（< 总能量的 5%）
    • 对高速应用则很重要

  3. 势能变化：

    • 对立式泵而言影响显著
    • 对卧式安装而言可忽略

  4. 热传递：

    • 向环境散热：Q̇loss
    • 可通过温度测量得出
    • 通常为轴功率的 1-5%

泵送流体的温升：

因效率损失导致的温升：

ΔT = (1 - η)·Ẇshaft / (ṁ·cp)

其中 η 为泵效率。

对于不可压缩液体：

ΔT = (1 - η)·g·H / cp

示例： 水泵，H = 100 m，η = 0.80，cp = 4186 J/kg-K

ΔT = (1 - 0.80)·9.81·100 / 4186 = 0.47°C
5. 第二定律分析（熵与效率）

熵平衡：

对于稳态流动的控制体：

ṁ(s₂ - s₁) = Q̇/T_boundary + Ṡgen

其中：

  • s = 比熵（J/kg-K）
  • Ṡgen = 熵产生速率（W/K）
  • T_boundary = 边界温度（K）

熵产生：

用于衡量过程中的不可逆性：

Ṡgen = ṁ(s₂ - s₁) - Q̇/T_boundary

对于绝热泵：

Ṡgen = ṁ(s₂ - s₁)

等熵效率：

将实际过程与理想等熵过程进行比较：

η_isentropic = (h₂s - h₁) / (h₂ - h₁)

其中：

  • h₂s = 在 P₂ 处、s₂s = s₁ 时的比焓（等熵）
  • h₂ = 出口处的实际比焓

对于不可压缩液体：

等熵功为：

w_isentropic = ∫v·dP = v·(P₂ - P₁) = (P₂ - P₁)/ρ

火用分析：

火用（可用能）表示最大有用功：

ψ = (h - h₀) - T₀(s - s₀) + c²/2 + gz

其中下标 0 表示死态（环境条件）。

火用损失：

Ẋdestroyed = T₀·Ṡgen

第二定律效率：

η_II = Ẋuseful / Ẋinput = 1 - Ẋdestroyed / Ẋinput

这表示泵对可用能量的利用效率。

按熵产生划分损失：

不同的损失机制会导致熵产生：

  1. 摩擦损失： 叶轮和蜗壳内的流体摩擦
  2. 冲击损失： 流动分离和入射损失
  3. 泄漏损失： 流体通过间隙回流
  4. 圆盘摩擦： 叶轮表面的黏性剪切
  5. 热传递： 与流体之间的不可逆热传递
6. 循环计算

与泵相关的热力学循环：

泵通常是单独的部件，但它们运行于更大的热力学循环中。

兰金循环（蒸汽发电厂）

组成部分：

  1. 锅炉（加入热量）
  2. 涡轮机（输出功）
  3. 冷凝器（排出热量）
  4. 给水泵（提高压力）

兰金循环中的泵分析：

状态点：

  • 状态 3：冷凝器出口（饱和液体或过冷液体）
  • 状态 4：锅炉入口（高压压缩液体）

泵功（等熵）：

w_pump,s = v₃·(P₄ - P₃) = h₄s - h₃

实际泵功：

w_pump = w_pump,s / η_pump

对循环效率的影响：

η_cycle = (w_turbine - w_pump) / q_in

在兰金循环中，泵功通常为涡轮机功的 1-3%。

布雷顿循环（燃气轮机）

不直接适用于液体泵，但与以下设备相关：

  • 燃气轮机中的燃料泵
  • 液体火箭推进剂泵
  • 高压流体系统

与兰金循环的主要区别：

  • 工质始终为气体
  • 采用压缩而非泵送
  • 温度比更高
制冷循环

适用场景：

  • 热泵中的制冷剂泵
  • 冷水机组
  • 低温泵送

组件：

  1. 蒸发器
  2. 压缩机（或液体系统中的泵）
  3. 冷凝器
  4. 膨胀阀

性能系数（COP）：

COP = Q_useful / W_input

其中泵功计入 W_input。

泵的应用

泵送流体的温升

物理机制：

泵的效率损失会将机械能转化为热能，导致温度升高。

计算方法：

根据能量平衡：

ΔT = (h₂ - h₁) / cp - v·(P₂ - P₁) / cp

液体的简化公式：

ΔT = (1 - η)·v·(P₂ - P₁) / cp

或者用扬程表示：

ΔT = (1 - η)·g·H / cp

影响温升的因素：

  1. **泵效率：**η 越低 → ΔT 越高
  2. **压力升高：**ΔP 越高 → ΔT 越高
  3. **比热容：**cp 越低 → ΔT 越高
  4. **运行工况：**偏离设计工况 → η 越低 → ΔT 越高

实际影响：

  • **冷流体：**可能导致周围环境结露
  • **热流体：**可能接近沸点，导致气蚀
  • **黏性流体：**可能显著升温（η 低，摩擦大）
  • **挥发性流体：**存在汽化风险
  • **聚合物：**温度升高可能导致降解

计算示例：

水泵：

  • H = 200 m，η = 0.75，cp = 4186 J/kg-K
  • ΔT = 0.25·9.81·200/4186 = 1.17°C

液压油：

  • ΔP = 200 bar，η = 0.70，ρ = 850 kg/m³，cp = 2000 J/kg-K
  • ΔT = 0.30·200e5/(850·2000) = 3.53°C

高黏度流体：

  • H = 50 m，η = 0.50（黏度较高导致效率较低），cp = 2500 J/kg-K
  • ΔT = 0.50·9.81·50/2500 = 0.98°C
效率与损失

总效率：

η = η_h · η_vol · η_mech

各组件效率：

  1. 水力效率（η_h）：

    • 计入流体摩擦损失和冲击损失
    • 典型值：0.85-0.95

  2. 容积效率（η_vol）：

    • 计入内部泄漏
    • 典型值：0.96-0.99

  3. 机械效率（η_mech）：

    • 计入轴承和密封摩擦
    • 典型值：0.95-0.98

热力学损失分析：

每种损失机制都会增加熵：

摩擦损失：

Ṡ_friction = ṁ·g·h_friction / T_avg

泄漏损失：

Ṡ_leakage = ṁ_leak·(s_discharge - s_suction)

热传递损失：

Ṡ_heat = Q̇_loss / T_boundary - Q̇_loss / T_fluid

总熵产：

Ṡ_total = Ṡ_friction + Ṡ_leakage + Ṡ_heat + Ṡ_mechanical

火用效率：

η_exergy = 1 - T₀·Ṡ_total / Ẇshaft

这能更全面地反映热力学性能。

传热考量

泵中的传热方式：

  1. 传导：

    • 通过泵壳壁
    • 通过轴传至轴承
    • 通过叶轮材料传递

  2. 对流：

    • 与泵送流体发生强制对流
    • 与周围空气发生自然对流
    • 内部再循环回路

  3. 辐射：

    • 通常可忽略，除非温度较高
    • 对热油泵（>150°C）可能影响显著

热损失估算：

根据泵壳壁：

Q̇_loss = h·A·(T_fluid - T_ambient)

其中：

  • h = 总传热系数（W/m²-K）
  • A = 表面积（m²）
  • h 的典型值 = 5-20 W/m²-K（自然对流、带保温层）

根据能量平衡：

Q̇_loss = Ẇshaft - ṁ·(h₂ - h₁) - ṁ·(c₂² - c₁²)/2 - ṁ·g·(z₂ - z₁)

对性能的影响：

  • 降低出口温度
  • 影响流体性质（黏度、密度）
  • 可能导致热梯度和变形
  • 对低温泵和高温泵而言至关重要

热管理：

冷却要求：

  • 外部冷却夹套
  • 机械密封的冲洗方案
  • 再循环回路中的换热器

保温：

  • 防止热损失（热流体）
  • 防止结露（冷流体）
  • 提高安全性（高温表面）
多级泵热力学

分级原理：

将较高扬程需求分配给多个叶轮：

H_total = H_stage1 + H_stage2 + ... + H_stage_n

热力学分析：

级间：

State 1 → Impeller 1 → State 2 → Diffuser 1 → State 3 →
Impeller 2 → State 4 → ... → Final discharge

每级温升：

ΔT_stage = (1 - η_stage)·g·H_stage / cp

总温升：

ΔT_total = Σ ΔT_stage = (1 - η_overall)·g·H_total / cp

级间冷却：

对于高扬程或高温应用：

Q̇_cooling = ṁ·cp·(T_beforecooling - T_aftercooling)

优点：

  • 降低最终排出温度
  • 改善流体性质（降低黏度）
  • 防止汽化
  • 延长泵的使用寿命

压力叠加效应：

随着压力逐级升高：

  • 密度可能变化（可压缩流体）
  • 泄漏流量增加（ΔP 更高）
  • 密封载荷增加
  • 结构要求发生变化

能量分配：

并非所有级的贡献都相同：

  • 第一级：NPSH 要求最高
  • 中间级：效率最高
  • 最后一级：承受的压力载荷最高

泵相关的常见循环

兰金循环（详解）

蒸汽电厂配置：

Boiler → Turbine → Condenser → Feedwater Pump → Boiler
    ↓        ↓         ↓              ↓
   High P   Power   Low P          High P
   High T  Output   Low T        Moderate T

状态点分析：

状态 1： 进入冷凝器的汽轮机排汽

  • P₁ = 0.05-0.1 bar（真空条件）
  • T₁ = P₁ 下的饱和温度
  • x₁ = 0.85-0.95（湿蒸汽）

状态 2： 冷凝器出口 / 泵入口

  • P₂ = P₁（压降极小）
  • T₂ = T_sat(P₂) - 过冷度（通常为 2-5°C）
  • x₂ = 0（饱和液体或过冷液体）

状态 3： 泵出口 / 锅炉入口

  • P₃ = 100-300 bar（高压）
  • T₃ ≈ T₂ + ΔT_pump（温升较小）
  • 压缩液体

状态 4： 锅炉出口 / 汽轮机入口

  • P₄ = P₃ - ΔP_boiler
  • T₄ = 500-600°C（过热蒸汽）
  • 过热蒸气

泵功计算：

等熵功：

w_pump,s = h₃s - h₂ = v₂·(P₃ - P₂)

对于 40°C 的水（v ≈ 0.001008 m³/kg）：

w_pump,s = 0.001008·(250e5 - 0.08e5) ≈ 2.52 kJ/kg

实际功：

w_pump = w_pump,s / η_pump = 2.52 / 0.85 ≈ 2.96 kJ/kg

温升：

ΔT = (h₃ - h₂) / cp ≈ 0.44 / 4.18 ≈ 0.11°C (negligible)

循环效率：

η_Rankine = (w_turbine - w_pump) / q_in

典型值：

  • 简单循环：30-40%
  • 再热循环：38-45%
  • 再生循环：40-50%

泵功占汽轮机功的比例虽小（约为 1-2%），但对系统运行至关重要。

制冷循环

带液泵的蒸气压缩循环：

部分工业制冷系统会使用液泵代替压缩机，或在压缩机之外再加液泵。

状态点：

  1. 蒸发器出口： 低压、低温、饱和蒸气
  2. 压缩机出口： 高压、高温、过热蒸气
  3. 冷凝器出口： 高压、中温、饱和液体
  4. 泵出口： 极高压、中温、压缩液体
  5. 膨胀阀出口： 低压、低温、气液两相混合物

泵的应用：

  • 液体过量供给系统： 循环制冷剂
  • CO₂ 跨临界系统： 泵送液态 CO₂
  • 复叠系统： 在不同压力等级之间泵送

COP 计算：

COP = Q_evaporator / (W_compressor + W_pump)

泵功通常远小于压缩机功。

闭式传热循环

应用：

  • 太阳能热系统
  • 地热热泵
  • 工业余热回收
  • 区域供热/制冷

热力学分析：

制热模式：

Q̇_delivered = ṁ·cp·(T_supply - T_return)
W_pump = ṁ·g·H / η_pump
COP_heating = Q̇_delivered / W_pump

制冷模式：

Q̇_removed = ṁ·cp·(T_return - T_supply)
COP_cooling = Q̇_removed / W_pump

优化：

  • 尽量降低流量（减少泵功）
  • 尽量增大 ΔT（对于给定负荷，可减少流量）
  • 平衡管道摩擦与泵的尺寸

总结

泵系统的热力学分析需要系统地应用以下内容：

  1. 控制体定义： 明确边界和假设
  2. 状态点识别： 全面评估各项性质
  3. 性质评估： 使用可靠来源提供的准确流体性质
  4. 第一定律： 对各种能量形式进行能量平衡核算
  5. 第二定律： 分析熵生成和效率
  6. 循环分析： 理解泵在更大系统中的作用

主要结果：

  • 功率需求
  • 温升
  • 效率分解
  • 损失机制
  • 优化机会

参考资料

详细方程请参见 equations-reference.md，计算示例请参见 examples.md。

评论

支持 Markdown