4. 苯-甲苯精馏塔RadFrac严格模型搭建
苯-甲苯精馏塔的严格模型对象与自由度
本例用 RadFrac 搭建苯-甲苯二元精馏塔的严格平衡级模型。与简捷法不同,RadFrac 对每一块理论板同时求解 MESH 方程(物料平衡 M、相平衡 E、摩尔分率加和 S、焓平衡 H),因此必须把塔的“几何与操作”信息补全:理论板数、进料板位置、进料状态、塔顶压力与全塔压降、回流比或馏出率、以及两个产品规格之一。苯-甲苯为理想近似体系,活度系数接近 1,用 Raoult 定律加理想气相即可,物性方法选 IDEAL 或 RK-SOAVE 均可,本例用 IDEAL 便于对照手算。
RadFrac 的自由度关系是搭建时的核心约束:给定进料(流量、组成、温度/压力、气相分率)后,还需指定 理论板数 N、进料板位置 NF、塔顶压力 Ptop、全塔压降 ΔP、以及两个操作/规格变量(如回流比 R 与馏出率 D/F,或 R 与塔顶苯纯度)。少给一个规格,模型欠定;多给一个,模型过定并报自由度错误。
搭建步骤与关键设置
- 新建流程,从
Columns选RadFrac放置,连接进料流FEED、塔顶DIST、塔底BOT,并连接冷凝器回流与再沸器返回线(RadFrac 内部自动生成,无需外接换热器)。 - 在
Setup → Specifications中设Number of stages = 12,Condenser = Total,Reboiler = Kettle,Valid phases = Vapor-Liquid。 - 在
Streams页设进料板Feed stage = 7,Above-Stage Pressure = 101.3 kPa,Stage Pressure Drop = 0。 - 在
Setup → Specifications的 Design Specs 区设两个操作变量:Reflux ratio = 2.0、Distillate rate = 45 kmol/h(对应 D/F=0.45)。 - 进料流设
F = 100 kmol/h、P = 101.3 kPa、Vapor fraction = 0、苯摩尔分率0.50,物性方法选IDEAL。 - 点击
Run,收敛后在Blocks → RadFrac → Profiles查看 T、P、x、y 沿板分布。
示例结果与判读
以下为按上述规格运行后应出现的预期趋势(示例值,用于对照判读,非实测数据):塔顶温度接近纯苯常压沸点 80.1 ℃,塔底接近甲苯沸点 110.6 ℃;温度沿板自上而下单调升高,进料板附近出现斜率变化(进料冷液或饱和液使该板温度略降后再升)。塔顶馏出液中苯摩尔分率约 0.95 以上,塔底釜液中苯摩尔分率约 0.03 以下,具体数值取决于 N、R、N_F 的组合。
| 观察项 | 位置 | 预期现象 | 判读要点 |
|---|---|---|---|
| 温度 T | 第 1 板 → 第 12 板 | 80.1 ℃ → 110.6 ℃ 单调升高 | 若出现温度倒挂,检查压降设置或进料板位置 |
| 苯液相分率 x | 塔顶 / 塔底 | ≈0.95 / ≈0.03 | 与 D/F、R 是否匹配;偏离过大说明规格不合理 |
| 苯气相分率 y | 各板 | y > x(苯为轻组分) | 验证相平衡方向正确 |
| 回流比 R | 冷凝器 | 等于设定值 2.0 | 若被 Design Spec 调整,说明规格冲突 |
| 收敛状态 | Control Panel | Converged,迭代次数 < 20 | 不收敛先查自由度与规格 |
若 RadFrac 报“Degrees of freedom”错误,优先检查是否同时指定了 R、D/F 和产品纯度三个规格;若报“Not converged”,可先把 R 调大(如 3.0)或把 N 增至 15 再试,收敛后再回调到目标规格。
本例用 IDEAL 物性方法仅适用于苯-甲苯这类近似理想体系;换成醇-水等强非理想体系时,必须改用 NRTL/UNIQUAC 等活度系数模型,否则塔顶组成与温度分布会明显失真。