Modelica密闭压力锅水蒸发建模：抖振与饱和线问题解决方案

抖振问题的核心根源 #

当前模型在水加热至饱和临界点出现的抖振（Chattering），本质是跨饱和线的硬切换逻辑在数值求解时的固有缺陷：

• 模型使用if noEvent(water.h<=hb)这类精确边界的条件判断，将模型划分为“水未饱和无蒸发”和“水超温开始蒸发”两个完全独立的离散模式
• 数值求解过程中，变量在临界点附近会因离散误差在边界两侧小幅浮动，触发模式反复切换，最终表现为仿真曲线的无规则震荡
• 原有蒸发速率的计算没有平滑过渡机制，一旦触发蒸发，蒸发量会直接从0跳变为热量除以汽化潜热，进一步放大数值震荡

原模型中触发震荡的核心代码块为：

if noEvent(water.h<=hb) then
  MflowEvap = 0;
  HflowEvap = 0;
else
  MflowEvap = (Qflow + WflowW)/dhEvap;
  HflowEvap = MflowEvap * hd;
end if;

正确建模的核心思路 #

按照物理过程的连续性建立方程，彻底消除离散模式切换：

• 输入热量全部作用于液态水，液态水未达饱和时，能量全部用于提升水的温度
• 当液态水达到对应压力下的饱和温度后，温度不再升高，额外输入的热量全部用于相变，不存在传质阻力带来的瞬时蒸发跳变
• 容器总容积固定，气液两相压力时刻相等，忽略壁面热容与热损失
• 全程使用连续可微的状态方程，让ODE求解器不需要在离散模式间反复跳转

蒸发速率的合理建模方法 #

放弃“根据焓值是否超过饱和值决定蒸发量”的硬逻辑，采用连续匹配的计算方式：

1. 定义连续的相变程度变量，替代原来的二元判断

   satDegree = max(0, min(1, (water.h - hb) / dhEvap))
   

   变量取值为0时代表水为未饱和液态，取值在0-1之间时代表处于两相共存状态，不需要离散切换
2. 蒸发量计算不再绑定硬边界：

   // 未饱和时蒸发量为0
   // 进入两相区后，蒸发量等于扣除液态水升温所需能量后的剩余热量除以汽化潜热
   MflowEvap_target = if satDegree <= 1e-6 then 0
                      else (Qflow + WflowW - Mw*water.cp*der(water.T)) / dhEvap;
   

3. 给蒸发量增加极小的一阶惯性约束，既不改变物理过程的精度，又能平滑数值跳变：

   parameter SI.Time tau_evap = 1e-3 "极小传质时间常数，避免数值震荡"
   tau_evap * der(MflowEvap) + MflowEvap = MflowEvap_target
   

   这个设置从根源上消除了蒸发量从0到非0的瞬时跳变，求解器不需要反复判断模式切换点。

消除跨饱和线仿真问题的最优方案 #

• 彻底替换noEvent硬判断：所有在饱和点的条件判断都改用连续的饱和程度变量表示，不要让求解器在离散事件点反复触发状态重置
• 保留完整的能量闭环：不要遗漏液态水受热膨胀带来的功项WflowW = -p*der(Vw)，这部分能量项的缺失会导致饱和点附近的能量不平衡，直接触发数值震荡
• 两相区状态约束的连续衔接：
  • 未饱和液态阶段：水的密度随温度变化，体积由质量和密度的关系Vw = Mw / water.d计算
  • 两相共存阶段：水的密度固定为对应压力下的饱和液密度，空气侧的体积为总容积减去剩余液态水的体积
  • 压力与温度的匹配：两相区满足water.T = Water.saturationTemperature(p)的热力学关系，让压力、温度、蒸发量三者形成自洽的闭环

内容的提问来源于stack exchange，提问作者T. Sergi