添加压力边界时出现奇异不一致标量系统错误求解决方案

我之前在使用ClaRa库的分流器模型搭建多用户供热系统时，遇到过完全相同的初始奇异问题（尤其是0/0的标量系统错误），给你几个亲测有效的排查和解决思路：

1. 排查初始条件的约束冲突 #

当添加第二个用户后，系统的初始约束数量增加，很容易出现过度约束或者约束矛盾，这是最常见的原因：

• 检查supply_pipe.volume的初始压力设置：如果同时给体积模块设置了固定初始压力，又通过压力边界的流量方程推导压力，就会出现冲突。建议手动给体积模块指定一个合理的初始压力（比如系统设计压力，例如1e5 Pa），并设置为fixed=true。
• 调整用户阀门的初始开度：不要将两个用户的阀门初始开度都设为0或全开。尝试给其中一个阀门设置极小的初始开度（比如0.01），避免初始时刻流量为0，导致压损计算的分母出现0。

2. 优化压损计算函数的鲁棒性 #

你提到用函数根据质量流量计算压差，这大概率是触发奇异错误的直接原因。要在函数里处理流量接近0的边界情况：

• 给流量判断添加一个极小的阈值，避免因浮点精度问题触发错误分支：

  // 修改前的判断逻辑
  -(if boundary.m_flow > 0.0 then 0.0 else 1.0)
  // 修改后添加epsilon阈值，规避浮点精度问题
  -(if boundary.m_flow > 1e-6 then 0.0 else 1.0)
  

• 当流量为0时，直接返回默认压差值（比如0），而不是执行除法运算，从根源避免0/0的奇异计算。

3. 调整Dymola初始化的优先级设置 #

Dymola的初始化逻辑会优先处理fixed=true的变量，你可以通过以下方式优化：

• 打开模型的Initialization面板，勾选Fixed initial conditions选项，手动指定关键变量的初始值，比如每个用户的初始流量、体积模块的压力。
• 将压力边界的初始流量设置为fixed=true，并给一个极小的初始值（比如1e-3 kg/s），让系统初始化时有明确的流量参考，避免无约束的奇异计算。

4. 逐步简化模型定位问题 #

既然单用户模型正常，你可以逐步添加第二个用户的组件，每次添加后运行模型，定位具体的问题组件：

• 先添加第二个用户的体积模块，不加阀门，检查是否报错；
• 再添加阀门，调整阀门的初始设置，看是否是阀门的初始约束导致的冲突。

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