微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

在 Dymola 中声明热容量 Cp

如何解决在 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;

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。