在 Dymola 中声明热容量 Cp

问题描述

我在调用我的工作流体的比热容时遇到了一些问题，在这种情况下是氢，我不能使用压力或温度来调用它，如果有人可以帮助我，请提前致谢。这是我的代码

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;

解决方法

我不确定您到底要达到什么目的，但您正在将 SaturationProperties 对象传递给需要 ThermodynamicState 的函数，该函数无法工作（并且在使用时报告为这样） OpenModelica）。

这是在 300 K 的饱和压力下计算 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;

dymola modelica