如何在Dymola中创建含固量可变的黑液介质以动态调整饱和特性？

核心思路 #

将物流追踪到的**实时含固量（重量百分比）**作为介质热力学属性计算的输入参数，替代原固定值，通过插值或经验公式关联不同含固量下的饱和沸点/压力特性。

具体步骤 #

1. 重构介质的T2P（温度转饱和压力）计算函数 #

• 打开你自定义的黑液介质模型，找到负责计算饱和压力的T2P函数，将其输入参数扩展为包含**含固量（solidMassFraction）**的变量，示例函数定义：

  function T2P
    input Real T "温度";
    input Real solidMassFraction "实时含固量（重量占比，0-1）";
    output Real P "对应含固量下的饱和压力";
  algorithm
    // 加入基于含固量和温度的饱和压力计算逻辑
  end T2P;
  

• 针对不同含固量下的饱和特性数据（比如15%、30%、50%、70%、75%含固量对应的温度-压力曲线），使用Dymola内置的Modelica.Math.Interpolations.tableLookup工具做二维插值，实现任意含固量下的饱和压力计算。例如整理成二维数据表：

  Real satData[:, :] = [
    [0.15, 373.15, 101325],  // 15%含固量，100℃对应标准大气压
    [0.15, 380, 133000],     // 15%含固量，380K对应压力
    [0.75, 373.15, 85000],   // 75%含固量，100℃对应压力
    // 补充更多实验/设计数据点
  ];
  P = Modelica.Math.Interpolations.tableLookup(satData, solidMassFraction, T);
  

2. 关联物流实时含固量到介质计算 #

• 在蒸发器模型中，利用instream追踪到的物流含固量变量（比如stream.solidMassFraction），将其传递给介质实例的属性计算参数。
• 确保在流单元（Stream）的介质绑定中，把T2P函数的solidMassFraction输入与当前物流的实时含固量变量关联，而非固定数值。

3. 验证与调试 #

• 搭建简化测试模型：让黑液从15%含固量开始，通过蒸发模块逐步提升到75%，实时监测不同含固量下的饱和温度/压力是否符合预期特性。
• 处理边界值：针对含固量超出15%-75%范围的情况，添加外推逻辑或告警，避免计算异常。

实操提示 #

• 如果是基于Modelica标准介质库扩展，尽量遵循原有库的接口规范，避免破坏模型兼容性；
• 若饱和特性有明确的经验公式，也可以直接将含固量作为公式变量代入，替代插值法；
• 注意变量因果性：确保含固量是由物流计算得出的输入变量，而非介质属性计算的输出，避免出现代数环问题。

内容的提问来源于stack exchange，提问作者drfeno