动态质量平衡通过将非线性汽液平衡关系替换为其局部线性形式,然后合并组分项来进行线性化。 对于第 (n) 级,液体组分平衡为
[ H_n\frac{dx_n}{dt}
Lx_{n-1}-Lx_n+Vy_{n+1}-Vy_n. ]
如果平衡关系在局部近似为 (y_i=m_i x_i),且斜率 (m_i) 为常数,则级方程将变为与相邻级耦合的线性微分方程。这种形式是推导传递函数、研究扰动传播以及为吸收和蒸馏中试装置设计组分控制系统的基础。
核心要点: 线性化通过将持液量、流速和局部平衡斜率视为常数,将级平衡转化为一阶动态模型。由此产生的模型通常是一组耦合的一阶方程组,而不是每个级的独立方程,因为每一级都与相邻级进行物质交换。
从动态级平衡开始
液体组分积累
对于恒定的液体持液量 (H_n),第 (n) 级上被吸收或蒸馏组分的积累为
[ H_n\frac{dx_n}{dt}. ]
这里,(x_n) 是第 (n) 级液体中的摩尔分数。单位必须一致:如果 (H_n) 是摩尔持液量,则 (L) 和 (V) 应为摩尔流速。
流量贡献
进入第 (n) 级的液体流贡献 (Lx_{n-1}),而离开的液体流贡献 (Lx_n)。
从相邻级进入的蒸汽流贡献 (Vy_{n+1}),离开第 (n) 级的蒸汽流移走 (Vy_n)。因此,
[ H_n\frac{dx_n}{dt}
Lx_{n-1}-Lx_n+Vy_{n+1}-Vy_n. ]
该平衡假设 (L) 和 (V) 恒定、持液量恒定、无反应,且气相中的积累可忽略不计。
应用局部平衡近似
线性平衡关系
蒸汽和液体组分通过平衡函数关联
[ y=f(x). ]
在选定的操作点附近,该关系可近似为
[ y_i \approx y_{i,0}+m_i(x_i-x_{i,0}), ]
其中
[ m_i=\left.\frac{dy}{dx}\right|{x{i,0}}. ]
对于具有恒定平衡线斜率且组分相对于操作点测量的简化模型,这简化为
[ y_i=m x_i. ]
斜率 (m) 可能代表局部挥发度或亨利定律关系,具体取决于系统。
代入级平衡方程
使用 (y_n=mx_n) 和 (y_{n+1}=mx_{n+1}),
[ H_n\frac{dx_n}{dt}
Lx_{n-1}-Lx_n+mVx_{n+1}-mVx_n. ]
合并包含 (x_n) 的项得到
[ \boxed{ H_n\frac{dx_n}{dt}
Lx_{n-1}
(L+mV)x_n + mVx_{n+1} } ]
或者,除以持液量后,
[ \boxed{ \frac{dx_n}{dt}
\frac{L}{H_n}x_{n-1}
\frac{L+mV}{H_n}x_n + \frac{mV}{H_n}x_{n+1} }. ]
这是恒定流量假设下理想平衡级的标准线性化动态方程。
该模型中“线性化”的含义
线性化平衡关系
原始平衡关系可能是非线性的,特别是对于相对挥发度变化显著的蒸馏系统。线性化用选定稳态操作点附近的切线替换了该非线性关系。
该近似对于操作点附近相对较小的组分变化是可靠的。大的扰动可能需要完整的非线性汽液平衡关系。
线性化物质平衡
当 (H)、(L) 和 (V) 为常数时,物质平衡本身在 (x) 中已经是线性的。当这些量随组分、温度、压力、液位或总流量变化时,会产生额外的非线性。
对于更通用的模型,可以写成
[ H_n(x_n)\frac{dx_n}{dt}
L(x_{n-1}-x_n) + V(y_{n+1}-y_n), ]
然后在稳态附近展开模型:
[ x_i=x_{i,0}+\delta x_i, \qquad y_i=y_{i,0}+\delta y_i. ]
平衡扰动近似为
[ \delta y_i=m_i\delta x_i. ]
由此产生的扰动模型具有与加框方程相同的结构,但其系数是在操作点处评估的。
重要的索引修正
必须保留相邻蒸汽组分
一个常见的索引错误是将 (y_{n+1}=mx_{n-1}) 代入。在通常的逆流级编号下,进入第 (n) 级的蒸汽是 (y_{n+1}),它与第 (n+1) 级的液体组分相关:
[ y_{n+1}=mx_{n+1}. ]
因此,正确索引的方程是
[ H_n\frac{dx_n}{dt}
Lx_{n-1}
(L+mV)x_n + mVx_{n+1}. ]
表达式
[ H\frac{dx_n}{dt}
Lx_{n-1}-(L+mV)x_n+mVx_{n-1} ]
仅当变量 (x_{n-1}) 被特意定义为进入蒸汽的等效上游液体组分时才有效。这不是通常的级索引平衡替换。
区分的重要性
保留 (x_{n+1}) 会产生一个耦合的塔模型。某一级的组分响应取决于两个相邻级,这对于表示通过塔板或离散化填料塔的扰动传播至关重要。
用相同的上游变量替换相邻级项可能会产生更简单的单级模型,但它改变了物理连接性,并可能给出错误的动态预测。
从级方程到传递函数
拉普拉斯域形式
假设初始扰动为零且系数为常数,进行拉普拉斯变换得到
[ H_n sX_n(s)
LX_{n-1}(s)
(L+mV)X_n(s) + mVX_{n+1}(s). ]
重新整理,
[ \boxed{ \left[H_n s+(L+mV)\right]X_n(s)
LX_{n-1}(s)+mVX_{n+1}(s) }. ]
该方程是塔的动态传递函数或状态空间模型的一行。
单级近似
如果将进入的蒸汽组分视为已知输入,而不是与 (x_{n+1}) 相关的状态,则定义
[ u_n=Lx_{n-1}+Vy_{n+1}. ]
平衡方程变为
[ H_n\frac{dx_n}{dt}
u_n-(L+mV)x_n. ]
相应的传递函数为
[ \boxed{ \frac{X_n(s)}{U_n(s)}
\frac{1}{H_n s+L+mV} } ]
时间常数为
[ \boxed{ \tau_n=\frac{H_n}{L+mV} }. ]
如果输入仅为上游液体组分,且蒸汽贡献被忽略或单独考虑,则
[ \frac{X_n(s)}{X_{n-1}(s)}
\frac{L}{H_n s+L+mV}. ]
在吸收和蒸馏中的应用
气液吸收
在吸收中,当实际气体组分超过其平衡值时,溶质通常从气相转移到液相:
[ y>y^*(x). ]
平衡斜率 (m) 描述了平衡气体组分随液体组分的变化情况。(mV) 项越大,气相耦合的强度就越强,并会改变级的响应。
对于填料吸收塔,相同的概念通常写为分布微分模型,而不是离散塔板平衡。有限体积离散化产生具有持液量和相邻单元耦合的级式方程。
蒸馏
在蒸馏中,平衡关系通常表示为
[ y_i=y^*(x_i), ]
其中 (y^*(x)) 可能是非线性的。在操作组分附近,
[ \delta y_i \approx \left.\frac{dy^*}{dx}\right|{x{i,0}}\delta x_i. ]
局部斜率可能因级而异,因此蒸馏塔通常具有系数 (m_n)、(L_n)、(V_n) 和 (H_n),而不是一个通用的 (m)、(L) 和 (V)。
边界条件
内部级方程本身不足以描述系统。冷凝器、再沸器、进料级、回流流、馏出物流、塔底流和塔端点提供了边界条件和操纵变量关系。
例如,回流改变了进入上级的液体流量,而再沸改变了进入下级的蒸汽流量。这些边界效应通常是控制扰动的主要路径。
将模型扩展到可变持液量和流量
总质量平衡
当液位或持液量变化时,必须包括总平衡。通用的液体持液量平衡为
[ \frac{dH_n}{dt}
L_{n-1}-L_n+V_{n+1}-V_n, ]
符号和流定义根据所选的级约定进行调整。
然后应根据组分持液量编写组分平衡:
[ \frac{d(H_nx_n)}{dt}
\text{组分输入}
\text{组分输出}. ]
展开导数得到
[ H_n\frac{dx_n}{dt} + x_n\frac{dH_n}{dt}
\text{组分输入}
\text{组分输出}. ]
当塔板液位或再沸器和冷凝器持液量发生明显变化时,额外的 (x_n,dH_n/dt) 项很重要。
能量和压力效应
对于精确的中试装置模型,温度会影响平衡和流速。压力会影响蒸汽组分、饱和条件和体积流量。
因此,仅组分的线性化主要适用于固定温度和压力操作点附近的适度扰动。较大的操作变化需要耦合的质量、能量和压力平衡方程。
理解权衡
线性化模型的优点
线性化模型紧凑、计算效率高,适用于传递函数分析。它提供了对时间常数、级间耦合、扰动衰减和可控性的直接洞察。
它特别适用于控制器整定和教学中试实验,在这些实验中,塔在稳定的标称条件下运行。
恒定平衡斜率的局限性
当相对挥发度随组分或温度发生显著变化时,单一的 (m) 值可能不足。该近似在共沸物、强非理想混合物或宽操作范围内也可能失效。
对于这些系统,请使用分级的局部斜率或保留非线性平衡函数。
恒定流量和持液量的局限性
假设 (L)、(V) 和 (H) 恒定排除了液位动力学和水力效应。这可能导致对响应速度的乐观预测,并可能掩盖流量、压力和组分控制之间重要的相互作用。
当中试装置可能出现塔板溢流、漏液、发泡或显著的液位变化时,应使用测量值或独立建模的持液量。
混淆操作线与动态平衡
操作线描述了相组分之间的稳态物质平衡关系。它本身并不描述组分随时间如何演变。
动态模型需要诸如 (H,dx/dt) 之类的持液项。移除这些项会将模型转换为稳态关系,而不是瞬态关系。
为您的目标做出正确的选择
合适的模型取决于目标是控制设计、概念分析还是高保真中试装置预测。
- 如果您的主要重点是控制器设计: 在预期的操作点附近线性化 (y^*(x)),保留相邻级耦合,并以状态空间形式组装方程。
- 如果您的主要重点是简单的单级响应: 将进入的蒸汽组分视为已知输入,并使用一阶传递函数 (1/[Hs+(L+mV)])。
- 如果您的主要重点是全塔扰动传播: 在每个级方程中保留 (x_{n-1}) 和 (x_{n+1}),并应用正确的冷凝器、再沸器、进料、回流和产品边界条件。
- 如果您的主要重点是大型中试装置瞬态: 使用非线性平衡、可变持液量、变化的流速以及耦合的能量和压力平衡。
一个健全的线性化塔模型在保留物理级连接性的同时,简化了明确定义的操作点附近的局部平衡行为。
总结表:
| 关键要素 | 线性化模型见解 |
|---|---|
| 动态平衡 | H_n dx_n/dt = Lx_{n-1} - Lx_n + Vy_{n+1} - Vy_n |
| 平衡近似 | 在稳态操作点附近用 delta y=m delta x 替换 y=f(x) |
| 耦合级方程 | H_n dx_n/dt = Lx_{n-1} - (L+mV)x_n + mVx_{n+1} |
| 主要假设 | 恒定持液量、液体和蒸汽流量,以及可忽略的气相积累 |
| 传递函数形式 | [H_n s + (L+mV)]X_n = LX_{n-1} + mVX_{n+1} |
| 单级时间常数 | tau_n = H_n/(L+mV) |
| 关键索引点 | 进入的蒸汽项与 x_{n+1} 相关,保留了相邻级耦合 |
| 最佳应用 | 控制器整定、扰动分析、传递函数和状态空间建模 |
| 何时需要非线性模型 | 大扰动、可变持液量或流量、强非理想性、压力变化或能量效应 |
需要用于教学、研究或工艺开发的中试装置吗?联系 LABPARK 讨论用于气液吸收、蒸馏、生物工艺和生物技术以及环境和水处理应用的中试装置。LABPARK 可以帮助大学、研究机构和企业选择支持实践培训、工艺演示、控制器整定、动态质量平衡研究和可扩展实验工作的设备。分享您的工艺、产能、操作条件以及学习或研究目标,我们的团队将为您推荐合适的中试装置配置。