OpenModelica新手求助：实现混合气体摩尔分数传感器的代码问题

问题描述 #

作为OpenModelica新手，希望结合以下两个MSL组件创建输出摩尔分数的传感器：

1. Modelica.Fluid.Sensors.MassFractions
2. Modelica.Media.Interfaces.PartialMixtureMedium.massToMoleFractions

计算摩尔分数需要引用介质中混合气体各组分的分子量，但编码时遇到两个问题：

1. 无法自动获取MMX数组大小
2. 无法自动获取MMX元素的摩尔质量值

补充说明：气体数据基于Modelica.Media.IdealGases及其Ar（氩气）创建，现有模型示例代码：

Modelica.Fluid.Sources.MassFlowSource_T inlet_A(
  redeclare package Medium = example2025.mixABC,
  use_m_flow_in = true,
  use_T_in = false,
  use_X_in = false,
  use_C_in = false,
  nPorts = 1,
  T = 373.15,
  X = {1.0, 0.0, 0.0}
)

Modelica.Fluid.Sensors.MassFractionsTwoPort XC2(
  redeclare package Medium = example2025.mixABC,
  substanceName = "gasC"
)

自行尝试编写的MolarFractionsTwoPort模型代码（存在问题）：

model MolarFractionsTwoPort
  extends Modelica.Fluid.Sensors.BaseClasses.PartialFlowSensor;
  extends Modelica.Icons.RoundSensor;
    
  Modelica.Blocks.Interfaces.RealOutput Xi "Mass fraction in port medium"
    annotation (Placement(transformation(extent={{-10,-10},{10,10}},
          rotation=90,
        origin={0,110}), iconTransformation(
        extent={{-10,-10},{10,10}},
        rotation=90,
        origin={0,110})));
  parameter String substanceName = "water" "Name of mass fraction";
  Real MMX[Medium.nXi]; //***(1)***
  
protected
  parameter Integer ind(fixed=false)
    "Index of species in vector of independent mass fractions";
       
initial algorithm
  ind:= -1;
  for i in 1:Medium.nXi loop
    if ( Modelica.Utilities.Strings.isEqual(Medium.substanceNames[i], substanceName)) then
      ind := i;
    end if;
  //MMX[i]:=Medium.fluidConstants[i].molarMass; //***(2)***
  MMX[i] := Modelica.Media.Interfaces.Types.Basic.FluidConstants[i].molarMass;//***(2)***
  end for;
.......
end MolarFractionsTwoPort;

解决方案 #

核心问题修正 #

1. MMX数组类型修正：组分摩尔质量是介质的常量参数，因此应定义为parameter Real数组，而非普通Real变量，编译阶段即可确定大小和值。
2. 摩尔质量获取方式修正：直接从介质包的fluidConstants中获取，即Medium.fluidConstants[i].molarMass，而非引用通用类型Basic.FluidConstants（不包含当前介质的具体组分数据）。
3. 组分索引初始化优化：用内置函数简化查找逻辑，并增加错误检查，避免未找到组分的运行错误。

完整修正模型代码 #

model MolarFractionsTwoPort
  extends Modelica.Fluid.Sensors.BaseClasses.PartialFlowSensor;
  extends Modelica.Icons.RoundSensor;

  Modelica.Blocks.Interfaces.RealOutput yi "Mole fraction of the specified substance"
    annotation (Placement(transformation(extent={{-10,-10},{10,10}}, rotation=90, origin={0,110}),
                          iconTransformation(extent={{-10,-10},{10,10}}, rotation=90, origin={0,110})));
  parameter String substanceName = "water" "Name of the substance to measure mole fraction for";

protected
  parameter Integer ind(fixed=false) "Index of the substance in the medium's substance list";
  parameter Real molarMasses[Medium.nXi] = [Medium.fluidConstants[i].molarMass for i in 1:Medium.nXi]
    "Molar masses of all components in the medium";
  Real X[Medium.nXi] "Mass fractions from the port medium";
  Real Y[Medium.nXi] "Mole fractions calculated from mass fractions";

initial algorithm
  ind := Modelica.Utilities.Strings.findSubstring(Medium.substanceNames, substanceName);
  if ind < 1 then
    Modelica.Utilities.Streams.error("Substance " + substanceName + " not found in medium " + Medium.name);
  end if;

equation
  // 获取当前介质的质量分数
  X = port_a.X;
  // 调用介质自带的质量分数转摩尔分数函数
  Y = Medium.massToMoleFractions(X, molarMasses);
  // 输出目标组分的摩尔分数
  yi = Y[ind];
end MolarFractionsTwoPort;

使用示例 #

结合你提供的混合介质模型，使用该传感器的代码如下：

model ExampleMolarSensor
  example2025.mixABC medium; // 你的混合介质包
  Modelica.Fluid.Sources.MassFlowSource_T inlet_A(
    redeclare package Medium = medium,
    use_m_flow_in = true,
    use_T_in = false,
    use_X_in = false,
    use_C_in = false,
    nPorts = 1,
    T = 373.15,
    X = {1.0, 0.0, 0.0}
  );
  MolarFractionsTwoPort yi_sensor(
    redeclare package Medium = medium,
    substanceName = "gasC"
  );
equation
  connect(inlet_A.port, yi_sensor.port_a);
end ExampleMolarSensor;

关键说明 #

• 直接复用PartialMixtureMedium内置的massToMoleFractions函数，无需手动实现计算逻辑。
• 通过列表推导式[Medium.fluidConstants[i].molarMass for i in 1:Medium.nXi]自动生成摩尔质量数组，解决数组大小和元素获取问题。
• 增加错误检查，避免未找到组分时的运行异常。

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