OpenModelica流体混合模拟咨询：单位转换与组分模拟问题

一、流量单位转换（Kg/s ↔ Nm³/h） #

Nm³指标准状态（0℃、101325Pa）下的体积，可通过理想气体状态方程推导转换关系：

1. 纯组分气体：
   转换系数 = (22.414e-3 m³/mol / 摩尔质量) * 3600 s/h
   例：H₂摩尔质量为2e-3 kg/mol，转换系数 = (22.414e-3 / 2e-3) * 3600 ≈ 40345.2 Nm³/(kg·h)，即1 kg/s H₂ = 40345.2 Nm³/h

2. N₂-H₂混合物：
   先计算混合物平均摩尔质量 M_avg = x_N2*M_N2 + x_H2*M_H2（x为摩尔分数，M_N2=28e-3 kg/mol，M_H2=2e-3 kg/mol）
   转换系数 = (22.414e-3 / M_avg) * 3600

界面实现方式： #

• 测量显示/曲线绘制：在模型中添加计算块，将massFlowRate（Kg/s）乘以转换系数，输出为Nm³/h变量绑定到界面。
  示例代码片段：

  Modelica.Blocks.Math.Product convertMassFlowToNm3h;
  convertMassFlowToNm3h.k = (22.414e-3 / M_avg) * 3600;
  connect(massFlowSensor.y, convertMassFlowToNm3h.u1);
  

• 设定值输入：将界面输入的Nm³/h除以转换系数，得到Kg/s后传入模型流量设定端口。

二、摩尔分数与质量分数的界面适配 #

1. 摩尔分数转质量分数（模型输入/内部计算） #

直接调用内置molartomassfraction函数，传入组分摩尔分数和各组分摩尔质量：

Modelica.Media.Interfaces.Types.MassFraction y = Modelica.Media.Utilities.Types.molartomassfraction(x, {28e-3, 2e-3});

（x为摩尔分数数组，顺序需与N₂、H₂组分一致）

2. 质量分数转摩尔分数（界面显示/设定） #

自定义转换函数实现反向转换：

function massToMolarFraction
  input Modelica.Media.Interfaces.Types.MassFraction y[:];
  input Real M[:]; // 各组分摩尔质量
  output Modelica.Media.Interfaces.Types.MoleFraction x[:];
algorithm
  x := y .* M / sum(y .* M);
end massToMolarFraction;

界面实现方式： #

• 测量显示：在容器质量分数输出端调用上述函数，转换为摩尔分数后绑定到界面控件或曲线。
• 设定值输入：将用户输入的摩尔分数转换为质量分数，再传入模型的组分设定参数（如注入流的质量分数）。

额外提示 #

确认基于NasaGasMixture创建的HNx介质已正确配置N₂和H₂的组分参数；大容器模型需确保质量平衡方程的正确性，避免因容积过大引发的数值稳定性问题。

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