packageR410a_IIR_P1_48_T233_340_Formula
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 R410a using a hybrid approach developed by Sangi et al.. The hybrid approach is implemented in AixLib.Media.Refrigerants.Interfaces.PartialHybridTwoPhaseMedium and the refrigerant model is implemented by complete the template AixLib.Media.Refrigerants.Interfaces.TemplateHybridTwoPhaseMediumRecord . The fitting coefficients required in the template are saved in the package AixLib.DataBase.Media.Refrigerants.R410a .
Assumptions and limitations
The implemented coefficients are fitted to external data by
Engelpracht and are valid within the following range:
|
Parameter |
Minimum Value |
Maximum Value |
|
Pressure (p) in bar |
1 |
48 |
|
Temperature (T) in K |
233.15 |
340.15 |
The reference point is defined as 200 kJ/kg and 1 kJ/kg/K, respectively, for enthalpy and entropy for the saturated liquid at 273.15 K.
Assumptions and limitations
R410a is calculated as pseudo-pure fluid and, hence, only roughly valid within the two-phase region.
Validation
The model is validated by comparing results obtained from the example model AixLib.Media.Refrigerants.Examples.RefrigerantProperties to external data (e.g. obtained from measurements or external media libraries).
References
Lemmon, E. W. (2003): Pseudo-Pure Fluid Equations of State for the Refrigerant Blends R-410A, R-404A, R-507A, and R-407C. In: International Journal of Thermophysics 24 (4), S. 991–1006. DOI: 10.1023/A:1025048800563.
Geller, V. Z.; Bivens, D.; Yokozeki, A. (2000): Viscosity of Mixed Refrigerants, R404A, R407C, R410A, and R507C. In: International refrigeration and air conditioning conference. USA, S. 399–406. Online available at http://docs.lib.purdue.edu/iracc/508.
Nabizadeh, H.; Mayinger, F. (1999): Viscosity of Gaseous R404A, R407C, R410A, and R507. In: International Journal of Thermophysics 20 (3), S. 777–790. DOI: 10.1007/978-1-4615-4777-8_1.
Geller, V. Z.; Nemzer, B. V.; Cheremnykh, U. V. (2001): Thermal Conductivity of the Refrigerant Mixtures R404A, R407C, R410A, and R507A. In: International Journal of Thermophysics 22 (4), 1035–1043. DOI: 10.1023/A:1010691504352.
Fröba, A. P.; Leipertz, A. (2003): Thermophysical Properties of the Refrigerant Mixtures R410A and R407C from Dynamic Light Scattering (DLS). In: International Journal ofThermophysics 24 (5), S. 1185–1206. DOI: 10.1023/A:1026152331710.
Engelpracht, Mirko (2017): Development of modular and scalable simulation models for heat pumps and chillers considering various refrigerants. Master Thesis
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| Modelica.Media.Interfaces.PartialTwoPhaseMedium.FluidConstants | refrigerantConstants | Thermodynamic constants for R410a | |
| Modelica.Media.Interfaces.Choices.IndependentVariables | ThermoStates (from PartialMedium) | Enumeration type for independent variables | |
| String | mediumName (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 |
| Boolean | singleState (from PartialMedium) | = true, if u and d are not a function of pressure | |
| Boolean | reducedX (from PartialMedium) | true | = true, if medium contains the equation sum(X) = 1.0; set reducedX=true, if only one substance (see docu for details) |
| Boolean | fixedX (from PartialMedium) | false | = true, if medium contains the equation X = reference_X |
| AbsolutePressure | reference_p (from PartialMedium) | 101325 | Reference pressure of Medium: default 1 atmosphere |
| Temperature | reference_T (from PartialMedium) | 298.15 | Reference temperature of Medium: default 25 deg Celsius |
| MassFraction[nX] | reference_X (from PartialMedium) | fill(1/nX, nX) | Default mass fractions of medium |
| AbsolutePressure | p_default (from PartialMedium) | 101325 | Default value for pressure of medium (for initialization) |
| Temperature | T_default (from PartialMedium) | Modelica.Units.Conversions.from_degC(20) | Default value for temperature of medium (for initialization) |
| SpecificEnthalpy | h_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_X | Default 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) |
| Integer | nS (from PartialMedium) | size(substanceNames, 1) | Number of substances |
| Integer | nX (from PartialMedium) | nS | Number of mass fractions |
| Integer | nXi (from PartialMedium) | if fixedX then 0 else if reducedX then nS - 1 else nS | Number of structurally independent mass fractions (see docu for details) |
| Integer | nC (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 |
| Boolean | smoothModel (from PartialTwoPhaseMedium) | false | = true, if the (derived) model should not generate state events |
| Boolean | onePhase (from PartialTwoPhaseMedium) | false | = true, if the (derived) model should never be called with two-phase inputs |
| FluidConstants | fluidConstants (from PartialTwoPhaseMedium) | Constant data for the fluid |
Contents
| Name | Description |
|---|---|
| SmoothTransition | Record that contains ranges to calculate a smooth transition between different regions |
| f_Idg | Dimensionless Helmholtz energy (Ideal gas contribution alpha_0) |
| f_Res | Dimensionless Helmholtz energy (Residual part alpha_r) |
| t_fIdg_t | Short form for tau*(dalpha_0/dtau)_delta=const |
| tt_fIdg_tt | Short form for tau*tau*(ddalpha_0/(dtau*dtau))_delta=const |
| t_fRes_t | Short form for tau*(dalpha_r/dtau)_delta=const |
| tt_fRes_tt | Short form for tau*tau*(ddalpha_r/(dtau*dtau))_delta=const |
| d_fRes_d | Short form for delta*(dalpha_r/(ddelta))_tau=const |
| dd_fRes_dd | Short form for delta*delta(ddalpha_r/(ddelta*delta))_tau=const |
| td_fRes_td | Short form for tau*delta*(ddalpha_r/(dtau*ddelta)) |
| ttt_fIdg_ttt | Short form for tau*tau*tau*(dddalpha_0/(dtau*dtau*dtau))_delta=const |
| ttt_fRes_ttt | Short form for tau*tau*tau*(dddalpha_r/(dtau*dtau*dtau))_delta=const |
| ddd_fRes_ddd | Short form for delta*delta*delta* (dddalpha_r/(ddelta*ddelta*ddelta))_tau=const |
| tdd_fRes_tdd | Short form for tau*delta*delta*(dddalpha_r/(dtau*ddelta*ddelta)) |
| ttd_fRes_ttd | Short form for tau*tau*delta*(dddalpha_r/(dtau*dtau*ddelta)) |
 saturationPressure | Saturation pressure of refrigerant (Ancillary equation) |
 saturationTemperature | Saturation temperature of refrigerant (Ancillary equation) |
 bubbleDensity | Boiling curve specific density of refrigerant (Ancillary equation) |
 dewDensity | Dew curve specific density of refrigerant (Ancillary equation) |
 bubbleEnthalpy | Boiling curve specific enthalpy of refrigerant (Ancillary equation) |
 dewEnthalpy | Dew curve specific enthalpy of refrigerant (Ancillary equation) |
 bubbleEntropy | Boiling curve specific entropy of refrigerant (Ancillary equation) |
 dewEntropy | Dew curve specific entropy of propane (Ancillary equation) |
 temperature_ph | Calculates temperature as function of pressure and specific enthalpy |
 temperature_ps | Calculates temperature as function of pressure and specific entroy |
 density_pT | Computes density as a function of pressure and temperature |
 dynamicViscosity | Calculates dynamic viscosity of refrigerant |
 thermalConductivity | Calculates thermal conductivity of refrigerant |
 surfaceTension | Surface tension in two phase region of refrigerant |
Revisions
- June 20, 2017, by Mirko Engelpracht, Christian Vering:
First implementation (see issue 408).