packageR290_IIR_P05_30_T263_343_Formula

Refrigerant model for R290 using a hybrid approach with explicit formulas developed by Sangi et al.

Extends from AixLib.Media.Refrigerants.Interfaces.PartialHybridTwoPhaseMediumFormula (Base class for two phase medium using a hybrid approach without records).

Information

This package provides a refrigerant model for R290 using a hybrid approach developed by Sangi et al.. The hybrid approach is implemented in AixLib.Media.Refrigerants.Interfaces.PartialHybridTwoPhaseMediumRecord and the refrigerant model is implemented by complete the template AixLib.Media.Refrigerants.Interfaces.TemplateHybridTwoPhaseMediumRecord .

Assumptions and limitations

The implemented coefficients are fitted to external data by Sangi et al. and are valid within the following range:

"Range of validiry" cellspacing="0" cellpadding="2" border="1" width="30%" style="border-collapse:collapse;">

Parameter

Minimum Value

Maximum Value

Pressure (p) in bar

0.5

30

Temperature (T) in K

263.15

343.15

Validation

Sangi et al. validated their model by comparing it to results obtained from the Helmholtz equation of state. They found out that relative error of the refrigerant model compared to HelmholtzMedia (Thorade and Saadat, 2012) is close to zero.

References

Thorade, Matthis; Saadat, Ali (2012): HelmholtzMedia - A fluid properties library. In: Proceedings of the 9th International Modelica Conference; September 3-5; 2012; Munich; Germany. Linköping University Electronic Press, S. 63–70.

Sangi, Roozbeh; Jahangiri, Pooyan; Klasing, Freerk; Streblow, Rita; Müller, Dirk (2014): A Medium Model for the Refrigerant Propane for Fast and Accurate Dynamic Simulations. In: The 10th International Modelica Conference. Lund, Sweden, March 10-12, 2014: Linköping University Electronic Press (Linköping Electronic Conference Proceedings), S. 1271–1275

Klasing,Freerk: A New Design for Direct Exchange Geothermal Heat Pumps - Modeling, Simulation and Exergy Analysis. Master thesis

Scalabrin, G.; Marchi, P.; Span, R. (2006): A Reference Multiparameter Viscosity Equation for Propane with an Optimized Functional Form. In: J. Phys. Chem. Ref. Data, Vol. 35, No. 3, S. 1415-1442

Parameters

TypeNameDefaultDescription
Modelica.Media.Interfaces.PartialTwoPhaseMedium.FluidConstantsrefrigerantConstantsThermodynamic constants for Propane
Modelica.Media.Interfaces.Choices.IndependentVariablesThermoStates (from PartialMedium)Enumeration type for independent variables
StringmediumName (from PartialMedium)"unusablePartialMedium"Name of the medium
String[:]substanceNames (from PartialMedium){mediumName}Names of the mixture substances. Set substanceNames={mediumName} if only one substance.
String[:]extraPropertiesNames (from PartialMedium)fill("", 0)Names of the additional (extra) transported properties. Set extraPropertiesNames=fill("",0) if unused
BooleansingleState (from PartialMedium)= true, if u and d are not a function of pressure
BooleanreducedX (from PartialMedium)true= true, if medium contains the equation sum(X) = 1.0; set reducedX=true, if only one substance (see docu for details)
BooleanfixedX (from PartialMedium)false= true, if medium contains the equation X = reference_X
AbsolutePressurereference_p (from PartialMedium)101325Reference pressure of Medium: default 1 atmosphere
Temperaturereference_T (from PartialMedium)298.15Reference temperature of Medium: default 25 deg Celsius
MassFraction[nX]reference_X (from PartialMedium)fill(1/nX, nX)Default mass fractions of medium
AbsolutePressurep_default (from PartialMedium)101325Default value for pressure of medium (for initialization)
TemperatureT_default (from PartialMedium)Modelica.Units.Conversions.from_degC(20)Default value for temperature of medium (for initialization)
SpecificEnthalpyh_default (from PartialMedium)specificEnthalpy_pTX(p_default, T_default, X_default)Default value for specific enthalpy of medium (for initialization)
MassFraction[nX]X_default (from PartialMedium)reference_XDefault value for mass fractions of medium (for initialization)
ExtraProperty[nC]C_default (from PartialMedium)fill(0, nC)Default value for trace substances of medium (for initialization)
IntegernS (from PartialMedium)size(substanceNames, 1)Number of substances
IntegernX (from PartialMedium)nSNumber of mass fractions
IntegernXi (from PartialMedium)if fixedX then 0 else if reducedX then nS - 1 else nSNumber of structurally independent mass fractions (see docu for details)
IntegernC (from PartialMedium)size(extraPropertiesNames, 1)Number of extra (outside of standard mass-balance) transported properties
Real[nC]C_nominal (from PartialMedium)1.0e-6*ones(nC)Default for the nominal values for the extra properties
BooleansmoothModel (from PartialTwoPhaseMedium)false= true, if the (derived) model should not generate state events
BooleanonePhase (from PartialTwoPhaseMedium)false= true, if the (derived) model should never be called with two-phase inputs
FluidConstantsfluidConstants (from PartialTwoPhaseMedium)Constant data for the fluid

Contents

NameDescription
SmoothTransitionRecord that contains ranges to calculate a smooth transition between different regions
f_IdgDimensionless Helmholtz energy (Ideal gas contribution alpha_0)
f_ResDimensionless Helmholtz energy (Residual part alpha_r)
t_fIdg_tShort form for tau*(dalpha_0/dtau)_delta=const
tt_fIdg_ttShort form for tau*tau*(ddalpha_0/(dtau*dtau))_delta=const
t_fRes_tShort form for tau*(dalpha_r/dtau)_delta=const
tt_fRes_ttShort form for tau*tau*(ddalpha_r/(dtau*dtau))_delta=const
d_fRes_dShort form for delta*(dalpha_r/(ddelta))_tau=const
dd_fRes_ddShort form for delta*delta(ddalpha_r/(ddelta*delta))_tau=const
td_fRes_tdShort form for tau*delta*(ddalpha_r/(dtau*ddelta))
ttt_fIdg_tttShort form for tau*tau*tau*(dddalpha_0/(dtau*dtau*dtau))_delta=const
ttt_fRes_tttShort form for tau*tau*tau*(dddalpha_r/(dtau*dtau*dtau))_delta=const
ddd_fRes_dddShort form for delta*delta*delta* (dddalpha_r/(ddelta*ddelta*ddelta))_tau=const
tdd_fRes_tddShort form for tau*delta*delta*(dddalpha_r/(dtau*ddelta*ddelta))
ttd_fRes_ttdShort form for tau*tau*delta*(dddalpha_r/(dtau*dtau*ddelta))
saturationPressureSaturation pressure of refrigerant (Ancillary equation)
saturationTemperatureSaturation temperature of refrigerant (Ancillary equation)
bubbleDensityBoiling curve specific density of refrigerant (Ancillary equation)
dewDensityDew curve specific density of refrigerant (Ancillary equation)
bubbleEnthalpyBoiling curve specific enthalpy of refrigerant (Ancillary equation)
dewEnthalpyDew curve specific enthalpy of refrigerant (Ancillary equation)
bubbleEntropyBoiling curve specific entropy of refrigerant (Ancillary equation)
dewEntropyDew curve specific entropy of propane (Ancillary equation)
temperature_phCalculates temperature as function of pressure and specific enthalpy
temperature_psCalculates temperature as function of pressure and specific entroy
density_pTComputes density as a function of pressure and temperature
dynamicViscosityCalculates dynamic viscosity of refrigerant
thermalConductivityCalculates thermal conductivity of refrigerant
surfaceTensionSurface tension in two phase region of refrigerant

Revisions

  • June 12, 2017, by Mirko Engelpracht, Christian Vering:
    First implementation (see issue 408).