Question for modeling Compressor in heat pump cycle.(PQ flash)(CoolProp)(ExternalMedia)(Dymola) - modelica

Question

I have probelm using ExternalMedia(Coolprop) in Dymola(modelica) modeling the compressor in heat pump system.
The error messages are like below
dymosim started
... "Project_WD.Test.Test_Compressor" simulating
... "dsin.txt" loading (dymosim input file)
Error: The following error was detected at time: 0
Pressure to PQ_flash [-509444 Pa] must be in range [0.0228908 Pa, 3.629e+06 Pa]
The stack of functions is:
Project_WD.Media.R600a_CP.setSat_p_Unique15
Project_WD.Media.R600a_CP.setSat_p_Unique15(compressor.p_su)
Project_WD.Media.R600a_CP.saturationTemperature_Unique13(compressor.p_su)
Non-linear solver will silently attempt to handle this problem.
Non-linear solver gave up after attempt: 305
----------------------------------------------------------------------------------------------------------------------
I think this is a CoolProp error, I did not used -50944 Pa in my script, and I wonder how to fix it.
I tried to fix the initial value of the p.su in the script. But It did not work. and I tried to make proper input value in Source and Sink But nothing changes.
The Main Code is like below
model Test_Compressor
ThermoCycle.Components.FluidFlow.Reservoirs.SinkP sinkPFluid(redeclare
package Medium = Project_WD.Media.R600a_CP, p0=2000000)
annotation (Placement(transformation(extent={{42,48},{62,68}})));
ThermoCycle.Components.FluidFlow.Reservoirs.SourceMdot sourceWF(
h_0=503925,
UseT=true,
Mdot_0=0.05,
redeclare package Medium = Project_WD.Media.R600a_CP,
T_0=323.15)
annotation (Placement(transformation(extent={{-94,-60},{-74,-40}})));
Component.compressor_BCA030NAMV compressor(redeclare package Medium =
Project_WD.Media.R600a_CP, redeclare function CPmodel =
Project_WD.Component_Parameter.BCA030NAMV)
annotation (Placement(transformation(extent={{-38,-14},{4,26}})));
equation
connect(sourceWF.flangeB, compressor.InFlow) annotation (Line(
points={{-75,-50},{-17,-50},{-17,-12}},
color={0,0,255},
smooth=Smooth.None));
connect(compressor.OutFlow, sinkPFluid.flangeB) annotation (Line(
points={{-17,24},{-18,24},{-18,58},{43.6,58}},
color={0,0,255},
smooth=Smooth.None));
annotation (Diagram(coordinateSystem(extent={{-120,-100},{80,100}},
preserveAspectRatio=false),
graphics), Icon(coordinateSystem(extent={{-120,-100},
{80,100}})),
experiment(StopTime=1000),
__Dymola_experimentSetupOutput);
end Test_Compressor;
The Sub Code is like below
model compressor_BCA030NAMV
/****************************************** FLUID ******************************************/
replaceable package Medium =
ThermoCycle.Media.R600a_CP
constrainedby Modelica.Media.Interfaces.PartialMedium "Medium model" annotation (choicesAllMatching = true);
/*Ports */
ThermoCycle.Interfaces.Fluid.FlangeA InFlow(redeclare package Medium =
Medium)
annotation (Placement(transformation(extent={{-10,-100},{10,-80}}),
iconTransformation(extent={{-10,-100},{10,-80}})));
ThermoCycle.Interfaces.Fluid.FlangeB OutFlow(redeclare package Medium =
Medium)
annotation (Placement(transformation(extent={{-10,80},{10,100}}),
iconTransformation(extent={{-10,80},{10,100}})));
replaceable function CPmodel =
Project_WD.Component_Parameter.BCA030NAMV constrainedby
Project_WD.calculate.compressor_para
"Compressor Model - LG compressor Map data model!" annotation (
choicesAllMatching=true);
/****************************************** PARAMETERES ******************************************/
parameter Modelica.Units.SI.Frequency N_rot=46 "Compressor frequency";
parameter Modelica.Units.SI.MassFlowRate M_dot_start=1.38e-05
"Nominal Mass flow rate" annotation (Dialog(tab="Initialization"));
parameter Modelica.Units.SI.Pressure p_su_start=100000
"Inlet pressure start value" annotation (Dialog(tab="Initialization"));
parameter Modelica.Units.SI.Pressure p_ex_start=1500000
"Outlet pressure start value" annotation (Dialog(tab="Initialization"));
parameter Modelica.Units.SI.Temperature T_su_start=253.15
"Inlet temperature start value" annotation (Dialog(tab="Initialization"));
parameter Medium.SpecificEnthalpy h_su_start = Medium.specificEnthalpy_pT(p_su_start, T_su_start)
"Inlet enthalpy start value" annotation (Dialog(tab="Initialization"));
parameter Medium.SpecificEnthalpy h_ex_start= Medium.specificEnthalpy_pT(p_ex_start, Medium.saturationTemperature(p_ex_start)+50)
"Outlet enthalpy start value" annotation (Dialog(tab="Initialization"));
/****************************************** VARIABLES ******************************************/
Medium.ThermodynamicState vaporIn
"Thermodynamic state of the fluid at the inlet";
Medium.ThermodynamicState vaporOut
"Thermodynamic state of the fluid at the outlet - isentropic";
Real epsilon_s "Isentropic Efficiency";
Real epsilon_v "Volumetric efficiency";
Modelica.Units.SI.Volume Vs "Swept Volume";
Real rpm;
Modelica.Units.SI.Power W_dot;
Modelica.Units.SI.VolumeFlowRate V_dot_su;
Modelica.Units.SI.MassFlowRate M_dot(start=M_dot_start);
Medium.Density rho_su(start=Medium.density_pT(p_su_start,T_su_start));
Medium.SpecificEntropy s_su;
Medium.SpecificEnthalpy h_su(start=h_su_start);
Medium.SpecificEnthalpy h_ex(start=h_ex_start);
Medium.AbsolutePressure p_su(start=p_su_start);
Medium.AbsolutePressure p_ex(start=p_ex_start);
Medium.SpecificEnthalpy h_ex_s;
Medium.Temperature T_cd(start=Medium.saturationTemperature(p_ex_start));
Medium.Temperature T_ev(start=Medium.saturationTemperature(p_su_start));
Real results[3];
equation
/* Fluid Properties */
vaporIn = Medium.setState_ph(p_su,h_su);
rho_su = Medium.density(vaporIn);
s_su = Medium.specificEntropy(vaporIn);
vaporOut = Medium.setState_ps(p_ex,s_su);
h_ex_s = Medium.specificEnthalpy(vaporOut);
T_cd = Medium.saturationTemperature(p_ex);
T_ev = Medium.saturationTemperature(p_su);
results=CPmodel(T_ev,T_cd);
M_dot = results[1];
W_dot = results[2];
Vs = results[3];
/*EQUATIONS */
rpm = N_rot*60;
V_dot_su = epsilon_v*Vs*N_rot;
V_dot_su = M_dot/rho_su;
h_ex = h_su + (h_ex_s - h_su)/epsilon_s;
W_dot = M_dot*(h_ex - h_su) "Consumed Power";
//BOUNDARY CONDITIONS //
/* Enthalpies */
h_su = if noEvent(InFlow.m_flow <= 0) then h_ex else inStream(InFlow.h_outflow);
h_su = InFlow.h_outflow;
OutFlow.h_outflow = if noEvent(OutFlow.m_flow <= 0) then h_ex else inStream(
OutFlow.h_outflow);
/*Mass flows */
M_dot = InFlow.m_flow;
OutFlow.m_flow = -M_dot;
/*pressures */
InFlow.p = p_su;
OutFlow.p = p_ex;
annotation (Icon(graphics={Bitmap(extent={{-80,-80},{80,80}}, fileName="modelica://LG_compressor/comp.png")}),
conversion(noneFromVersion=""));
end compressor_BCA030NAMV;
function compressor_para
"Base model for the computation of the compressor performance. Must be redeclared with the coefficients and the swept volume"
input Modelica.Units.SI.Temperature T_ev "Condensing temperature";/* unit : K */
input Modelica.Units.SI.Temperature T_cd "Evaporation temperature";/* unit : K */
output Real result[3]
"Vector with the output flow rate, cp power and swept volume";
protected
parameter Real coef[11,4]=
[0,0,0,0;
0,0,0,0;
0,0,0,0;
0,0,0,0;
0,0,0,0;
0,0,0,0;
0,0,0,0;
0,0,0,0;
0,0,0,0;
0,0,0,0;
0,0,0,0]
"Coefficient matrix. Column1: Capacity. Column2: Power. Column3: Current. Column4: Flow Rate";
parameter Modelica.Units.SI.Volume Vs=0 "Compressor Swept Volume";
protected
Modelica.Units.SI.MassFlowRate Mdot "Mass Flow Rate";/* unit : kg/s */
Modelica.Units.SI.Power Wdot "Compressor Power";/* unit : W */
algorithm
Wdot := Project_WD.calculate.Comp_mapdata_cal(
(T_ev - 273.15),
(T_cd - 273.15),
coef[:, 2]); // Temperature input is in Celcius
Mdot := Project_WD.calculate.Comp_mapdata_cal(
(T_ev - 273.15),
(T_cd - 273.15),
coef[:, 4]); // Temperature input is in Celcius. Output is in kg/hr
result:={Mdot,Wdot,Vs};
end compressor_para;
function BCA030NAMV
"BCA030NAMV BLDC/Compact/A/Displacement 03.0cc/rev = 3E-06/rev / R600a / Cu Wire / 220-240V, 50/60Hz / inverter Compressor by LG elec"
extends Project_WD.calculate.compressor_para(Vs=3.0e-6, coef=[0,87.0312493,0,2.73823813;
0,-23.5090917,0,-0.314445833; 0,-0.576825,0,-0.0085375; 0,-0.001413333,0,
-2.66667E-05; 0,-7.740443,0,-0.116178854; 0,0.234250512,0,0.003251526; 0,
-0.001864852,0,-2.81759E-05; 0,0.80516,0,0.01332; 0,-0.006655,0,-0.0001235;
0,0.01724,0,0.00028; 0,-0.000143,0,-0.0000025]);
end BCA030NAMV;
The Fluid Property I used is below
package R600a_CP "CoolProp model of Isobutane"
extends ExternalMedia.Media.CoolPropMedium(
mediumName = "Isobutane",
substanceNames = {"Isobutane"},
ThermoStates = Modelica.Media.Interfaces.Choices.IndependentVariables.ph,
SpecificEnthalpy(start=2e5));
end R600a_CP;