在Dymola中声明热容Cp

Question

我有一些问题要称为我的工质的比热容量，在这种情况下是氢气，我不能用压力或温度来称它，如果有人能帮我的话，请提前谢谢。以下是我的代码

import Modelica.SIunits;

  package Hyd
    extends ExternalMedia.Media.CoolPropMedium(
      mediumName="hydrogen",
      substanceNames={"hydrogen"},
      inputChoice=ExternalMedia.Common.InputChoice.pT);
  end Hyd;

SIunits.SpecificHeatCapacity cp_in;//[J/kg*K]
Hyd.AbsolutePressure Pb_0;
Hyd.Temperature Tin;
Hyd.SaturationProperties sat9,sat10;

Equation
sat9=Hyd.setSat_T(Tin);
sat10=Hyd.setSat_p(Pb_0);
cp_in=Hyd.specificHeatCapacityCp(sat9);//[J/kg*K]
cp_in=Hyd.specificHeatCapacityCp(sat10);//[J/kg*K]

该函数声明为：

function specificHeatCapacityCp_Unique8 
input ExternalMedia.Media.BaseClasses.ExternalTwoPhaseMedium.ThermodynamicState state ;
output Modelica.Media.Interfaces.Types.SpecificHeatCapacity cp := 1000.0 "Specific heat capacity at constant pressure";
end specificHeatCapacityCp_Unique8;

Answer 1

我不确定您到底想要实现什么，但是您正在将一个SaturationProperties对象传递给一个需要ThermodynamicState的函数，这是不能工作的(并且在使用OpenModelica时报告是这样的)。

这是一个在300K饱和压力下计算cp的工作版本：

model test_SO_68546587
import Modelica.SIunits;

  package Hyd
    extends ExternalMedia.Media.CoolPropMedium(
      mediumName="hydrogen",
      substanceNames={"hydrogen"},
      inputChoice=ExternalMedia.Common.InputChoice.pT);
  end Hyd;

SIunits.SpecificHeatCapacity cp_in;//[J/kg*K]
Hyd.AbsolutePressure Pb_0;
Hyd.Temperature Tin;
Hyd.ThermodynamicState state;
equation
state = Hyd.setState_pT(p=Pb_0, T=Tin);
Tin = 300;
Pb_0 = Hyd.saturationPressure(Tin);
cp_in=Hyd.specificHeatCapacityCp(state);// 14345.2 J/kg*K @ 300 K, 12.951 bar
end test_SO_68546587;