OpenModelica中H2+H2O介质建模仿真及H2与饱和水分离问询

一、正确构建水-氢气两相介质模型 #

直接替换湿空气模型的参数会引发数值异常，因为湿空气的热力学框架是针对空气-水蒸气体系设计的，与氢气-水体系的物性差异极大。建议从Modelica.Media的基础框架出发构建：

• 选择适配的基类：使用Modelica.Media.Interfaces.PartialTwoPhaseMedium作为介质基类，该类专为两相介质建模设计，内置相平衡基础逻辑。
• 明确相组成定义：
  • 气相：构建氢气+水蒸气的理想气体混合物，采用Modelica.Media.Mixtures.IdealGases.SingleGases中的H2和Modelica.Media.Substances.PureWater的水蒸气模型作为组分，确保各组分物性参数（临界温度、饱和蒸气压方程等）准确。
  • 液相：以纯水作为液相介质，直接调用Modelica.Media.Substances.PureWater.Liquid或两相模型，保证气液两相满足热力学平衡约束（温度相等、水蒸气分压等于当前温度下的饱和蒸气压）。
• 排查数值异常：
  • 检查初始化条件：确保压力、温度、干度等状态变量的初始值符合实际工况，比如氢气临界温度仅33.2K，常温下必然处于气相，避免不合理的相态假设。
  • 单独验证纯物质物性：先测试纯氢气和纯水在目标工况下的物性计算是否正常，再逐步组合成两相介质，定位问题点。
• 介质模型框架示例：

model H2WaterTwoPhaseMedium "Hydrogen-Water two-phase medium"
  extends Modelica.Media.Interfaces.PartialTwoPhaseMedium;
  // 定义气相混合物
  replaceable package GasPhase = Modelica.Media.Mixtures.IdealGases.Mixture(
    substances={Modelica.Media.Substances.PureWater, Modelica.Media.Substances.H2},
    substanceNames={"Water", "Hydrogen"});
  // 定义液相介质
  replaceable package LiquidPhase = Modelica.Media.Substances.PureWater.Liquid;
  // 完善相平衡逻辑（继承基类后补充平衡状态计算）
end H2WaterTwoPhaseMedium;

二、气液分离器模型开发建议 #

Modelica标准库4.0中无现成的通用气液分离器模型，但可基于热力学平衡原理开发简化版：

简化版分离器核心逻辑 #

假设分离器内气液两相达到完全热力学平衡，满足质量守恒、能量守恒，且气液两相温度、压力一致。

实现步骤 #

• 端口定义：设置三个流体端口：入口（两相流）、气相出口（单相气相）、液相出口（单相液相）。
• 守恒方程：
  • 总质量守恒：m_flow_in = m_flow_gas + m_flow_liquid
  • 组分质量守恒：氢气仅存在于气相，故m_flow_H2_in = m_flow_H2_gas；水的质量流量满足m_flow_water_in = m_flow_water_vapor + m_flow_water_liquid
  • 能量守恒：h_in * m_flow_in = h_gas * m_flow_gas + h_liquid * m_flow_liquid
• 平衡约束：气相中水蒸气的分压等于当前温度下纯水的饱和蒸气压，利用Modelica.Media.Substances.PureWater.saturationPressure(T)计算饱和蒸气压，进而确定气相中水蒸气的摩尔分数。
• 初始化设置：给分离器设置初始压力、温度和液位高度，避免仿真初期的数值震荡。

核心代码片段示例 #

model SimpleGasLiquidSeparator
  Modelica.Fluid.Interfaces.FluidPort_a inlet;
  Modelica.Fluid.Interfaces.FluidPort_a gasOutlet;
  Modelica.Fluid.Interfaces.FluidPort_a liquidOutlet;
  
  // 状态变量
  Modelica.SIunits.Pressure p;
  Modelica.SIunits.Temperature T;
  Modelica.SIunits.MassFlowRate m_flow_gas, m_flow_liquid;
  
equation
  // 压力平衡：所有端口压力一致
  inlet.p = p;
  gasOutlet.p = p;
  liquidOutlet.p = p;
  
  // 温度平衡：气液两相温度一致
  inlet.T = T;
  gasOutlet.T = T;
  liquidOutlet.T = T;
  
  // 总质量守恒
  inlet.m_flow = m_flow_gas + m_flow_liquid;
  
  // 氢气组分守恒（仅存在于气相）
  inlet.X[2] * inlet.m_flow = gasOutlet.X[2] * m_flow_gas;
  
  // 水的组分守恒
  inlet.X[1] * inlet.m_flow = gasOutlet.X[1] * m_flow_gas + liquidOutlet.X[1] * m_flow_liquid;
  
  // 能量守恒
  inlet.h * inlet.m_flow = gasOutlet.h * m_flow_gas + liquidOutlet.h * m_flow_liquid;
  
  // 相平衡：气相中水的分压等于饱和蒸气压
  gasOutlet.X[1] * p = Modelica.Media.Substances.PureWater.saturationPressure(T);
end SimpleGasLiquidSeparator;

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