Modelica.Media.Common.ThermoFluidSpecial

Package Content

Name	Description
FixedIGProperties	constant properties for ideal gases
ThermoBaseVars
ThermoProperties	Thermodynamic base property data for all state models
ThermoProperties_ph	Thermodynamic property data for pressure p and specific enthalpy h as dynamic states
ThermoProperties_pT	Thermodynamic property data for pressure p and temperature T as dynamic states
ThermoProperties_dT	Thermodynamic property data for density d and temperature T as dynamic states
TransportProps	record with transport properties
gibbsToProps_ph	calulate property record for pressure and specific enthalpy as states from dimensionless Gibbs function
gibbsToBoundaryProps	calulate phase boundary property record from dimensionless Gibbs function
gibbsToProps_dT	calulate property record for density and temperature as states from dimensionless Gibbs function
gibbsToProps_pT	calulate property record for pressure and temperature as states from dimensionless Gibbs function
helmholtzToProps_ph	calulate property record for pressure and specific enthalpy as states from dimensionless Helmholtz function
helmholtzToProps_pT	calulate property record for pressure and temperature as states from dimensionless Helmholtz function
helmholtzToProps_dT	calulate property record for density and temperature as states from dimensionless Helmholtz function
TwoPhaseToProps_ph	compute property record for pressure and specific enthalpy as states from saturation properties
TwoPhaseToProps_dT	compute property record for density and temperature as states from saturation properties

Modelica.Media.Common.ThermoFluidSpecial.FixedIGProperties

Information

Parameters

Type	Name	Default	Description
Integer	nspecies		number of components

Modelica definition

record FixedIGProperties "constant properties for ideal gases"
  extends Modelica.Icons.Record;
  parameter Integer nspecies(min=1) "number of components";
  SI.MolarMass[nspecies] MM "molar mass of components";
  Real[nspecies] invMM "inverse of molar mass of components";
  SI.SpecificHeatCapacity[nspecies] R "gas constant";
  SI.SpecificEnthalpy[nspecies] Hf "enthalpy of formation at 298.15K";
  SI.SpecificEnthalpy[nspecies] H0 "H0(298.15K) - H0(0K)";
end FixedIGProperties;

Modelica.Media.Common.ThermoFluidSpecial.ThermoBaseVars

Information

Model description

ThermoBaseVars is inherited by all medium property models and by all models defining the dynamic states for the conservation of mass and energy. Thus it is a good choice as a restricting class for any medium model or dynamic state model.

Extends from Modelica.Icons.Record (Icon for records).

Parameters

Type	Name	Default	Description
Integer	n		discretization number
Integer	nspecies		number of species

Modelica definition

record ThermoBaseVars
  extends Modelica.Icons.Record;
  parameter Integer n(min=1) "discretization number";
  parameter Integer nspecies(min=1) "number of species";
  SI.Pressure[n] p(
    min=PMIN,
    max=PMAX,
    nominal=PNOM,
    start=fill(1.0e5, n)) "Pressure";
  SI.Temperature[n] T(
    min=TMIN,
    max=TMAX,
    nominal=TNOM) "temperature";
  SI.Density[n] d(
    min=DMIN,
    max=DMAX,
    nominal=DNOM) "density";
  SI.SpecificEnthalpy[n] h(
    min=SHMIN,
    max=SHMAX,
    nominal=SHNOM) "specific enthalpy";
  SI.SpecificEntropy[n] s(
    min=SSMIN,
    max=SSMAX,
    nominal=SSNOM) "specific entropy";
  SI.RatioOfSpecificHeatCapacities[n] kappa "ratio of cp/cv";
  SI.Mass[n] M(
    min=MMIN,
    max=MMAX,
    nominal=MNOM) "Total mass";
  SI.Energy[n] U(
    min=EMIN,
    max=EMAX,
    nominal=ENOM) "Inner energy";
  SI.MassFlowRate[n] dM(
    min=MDOTMIN,
    max=MDOTMAX,
    nominal=MDOTNOM) "Change in total mass";
  SI.Power[n] dU(
    min=POWMIN,
    max=POWMAX,
    nominal=POWNOM) "Change in inner energy";
  SI.Volume[n] V(
    min=VMIN,
    max=VMAX,
    nominal=VNOM) "Volume";
  SI.MassFraction[n,nspecies] mass_x(
    min=MASSXMIN,
    max=MASSXMAX,
    nominal=MASSXNOM) "mass fraction";
  SI.MoleFraction[n,nspecies] mole_y(
    min=MOLEYMIN,
    max=MOLEYMAX,
    nominal=MOLEYNOM) "mole fraction";
  SI.Mass[n,nspecies] M_x(
    min=MMIN,
    max=MMAX,
    nominal=MNOM) "component mass";
  SI.MassFlowRate[n,nspecies] dM_x(
    min=MDOTMIN,
    max=MDOTMAX,
    nominal=MDOTNOM) "rate of change in component mass";
  MolarFlowRate[n, nspecies] dZ(
    min=-1.0e6,
    max=1.0e6,
    nominal=0.0) "rate of change in component moles";
  MolarFlowRate[n, nspecies] rZ(
    min=-1.0e6,
    max=1.0e6,
    nominal=0.0) "Reaction(source) mole rates";
  SI.MolarMass[n] MM(
    min=MMMIN,
    max=MMMAX,
    nominal=MMNOM) "molar mass of mixture";
  SI.AmountOfSubstance[n] Moles(
    min=MOLMIN,
    max=MOLMAX,
    nominal=MOLNOM) "total moles";
  SI.AmountOfSubstance[n,nspecies] Moles_z(
    min=MOLMIN,
    max=MOLMAX,
    nominal=MOLNOM) "mole vector";
end ThermoBaseVars;

Modelica.Media.Common.ThermoFluidSpecial.ThermoProperties

Information

Model description

A base class for medium property models which work with most of the versions of dynamic states that are available in the ThermoFluid library. Currently used by all ideal gas models.

Extends from Modelica.Icons.Record (Icon for records).

Parameters

Type	Name	Default	Description
Integer	nspecies		number of species

Modelica definition

record ThermoProperties 
  "Thermodynamic base property data for all state models"
  extends Modelica.Icons.Record;
  parameter Integer nspecies(min=1) "number of species";
  SI.Temperature T(
    min=TMIN,
    max=TMAX,
    nominal=TNOM) "temperature";
  SI.Density d(
    min=DMIN,
    max=DMAX,
    nominal=DNOM) "density";
  SI.Pressure p(
    min=PMIN,
    max=PMAX,
    nominal=PNOM) "pressure";
  SI.Volume V(
    min=VMIN,
    max=VMAX,
    nominal=VNOM) "Volume";
  SI.SpecificEnthalpy h(
    min=SHMIN,
    max=SHMAX,
    nominal=SHNOM) "specific enthalpy";
  SI.SpecificEnergy u(
    min=SEMIN,
    max=SEMAX,
    nominal=SENOM) "specific inner energy";
  SI.SpecificEntropy s(
    min=SSMIN,
    max=SSMAX,
    nominal=SSNOM) "specific entropy";
  SI.SpecificGibbsFreeEnergy g(
    min=SHMIN,
    max=SHMAX,
    nominal=SHNOM) "specific Gibbs free energy";
  SI.SpecificHeatCapacity cp(
    min=CPMIN,
    max=CPMAX,
    nominal=CPNOM) "heat capacity at constant pressure";
  SI.SpecificHeatCapacity cv(
    min=CPMIN,
    max=CPMAX,
    nominal=CPNOM) "heat capacity at constant volume";
  SI.SpecificHeatCapacity R(
    min=CPMIN,
    max=CPMAX,
    nominal=CPNOM) "gas constant";
  SI.MolarMass MM(
    min=MMMIN,
    max=MMMAX,
    nominal=MMNOM) "molar mass of mixture";
  SI.MassFraction[nspecies] mass_x(
    min=MASSXMIN,
    max=MASSXMAX,
    nominal=MASSXNOM) "mass fraction";
  SI.MoleFraction[nspecies] mole_y(
    min=MOLEYMIN,
    max=MOLEYMAX,
    nominal=MOLEYNOM) "mole fraction";
  SI.RatioOfSpecificHeatCapacities kappa "ratio of cp/cv";
  SI.DerDensityByTemperature ddTp 
    "derivative of density by temperature at constant pressure";
  SI.DerDensityByPressure ddpT 
    "derivative of density by pressure at constant temperature";
  Real dupT(unit="m3.kg-1") 
    "derivative of inner energy by pressure at constant T";
  Real dudT(unit="(J.m3)/(kg2)") 
    "derivative of inner energy by density at constant T";
  SI.SpecificHeatCapacity duTp 
    "derivative of inner energy by temperature at constant p";
  SI.SpecificEnergy ddx[nspecies] 
    "derivative vector of density by change in mass composition";
  SI.SpecificEnergy[nspecies] compu(
    min=SEMIN,
    max=SEMAX,
    nominal=SENOM) "inner energy of the components";
  SI.Pressure[nspecies] compp(
    min=COMPPMIN,
    max=COMPPMAX,
    nominal=COMPPNOM) "partial pressures of the components";
  SI.Velocity a(
    min=VELMIN,
    max=VELMAX,
    nominal=VELNOM) "speed of sound";
  SI.HeatCapacity dUTZ 
    "derivative of inner energy by temperature at constant moles";
  SI.MolarInternalEnergy[nspecies] dUZT 
    "derivative of inner energy by moles at constant temperature";
  SI.SpecificEnthalpy[nspecies] dHMxT(
    min=SEMIN,
    max=SEMAX,
    nominal=SENOM) 
    "derivative of total enthalpy w.r.t. component mass at constant T";
  Real dpT "derivative of pressure w.r.t. temperature";
  Real dpZ[nspecies] "derivative of pressure w.r.t. moles";
end ThermoProperties;

Modelica.Media.Common.ThermoFluidSpecial.ThermoProperties_ph

Information

Model description

A base class for medium property models which use pressure and enthalpy as dynamic states. This is the preferred model for fluids that can also be in the two phase and liquid regions.

Extends from Modelica.Icons.Record (Icon for records).

Modelica definition

record ThermoProperties_ph 
  "Thermodynamic property data for pressure p and specific enthalpy h as dynamic states"

  extends Modelica.Icons.Record;
  SI.Temperature T(
    min=1.0e-9,
    max=10000.0,
    nominal=298.15) "temperature";
  SI.Density d(
    min=1.0e-9,
    max=10000.0,
    nominal=10.0) "density";
  SI.SpecificEnergy u(
    min=-1.0e8,
    max=1.0e8,
    nominal=1.0e6) "specific inner energy";
  SI.SpecificEntropy s(
    min=-1.0e6,
    max=1.0e6,
    nominal=1.0e3) "specific entropy";
  SI.SpecificHeatCapacity cp(
    min=1.0,
    max=1.0e6,
    nominal=1000.0) "heat capacity at constant pressure";
  SI.SpecificHeatCapacity cv(
    min=1.0,
    max=1.0e6,
    nominal=1000.0) "heat capacity at constant volume";
  SI.SpecificHeatCapacity R(
    min=1.0,
    max=1.0e6,
    nominal=1000.0) "gas constant";
  SI.RatioOfSpecificHeatCapacities kappa "ratio of cp/cv";
  SI.Velocity a(
    min=1.0,
    max=10000.0,
    nominal=300.0) "speed of sound";
  SI.DerDensityByEnthalpy ddhp 
    "derivative of density by enthalpy at constant pressure";
  SI.DerDensityByPressure ddph 
    "derivative of density by pressure at constant enthalpy";
  Real duph(unit="m3/kg") 
    "derivative of inner energy by pressure at constant enthalpy";
  Real duhp(unit="1") 
    "derivative of inner energy by enthalpy at constant pressure";
end ThermoProperties_ph;

Modelica.Media.Common.ThermoFluidSpecial.ThermoProperties_pT

Information

Model description

A base class for medium property models which use pressure and temperature as dynamic states. This is a reasonable model for fluids that can also be in the gas and liquid regions, but never in the two-phase region.

Extends from Modelica.Icons.Record (Icon for records).

Modelica definition

record ThermoProperties_pT 
  "Thermodynamic property data for pressure p and temperature T as dynamic states"

  extends Modelica.Icons.Record;
  SI.Density d(
    min=1.0e-9,
    max=10000.0,
    nominal=10.0) "density";
  SI.SpecificEnthalpy h(
    min=-1.0e8,
    max=1.0e8,
    nominal=1.0e6) "specific enthalpy";
  SI.SpecificEnergy u(
    min=-1.0e8,
    max=1.0e8,
    nominal=1.0e6) "specific inner energy";
  SI.SpecificEntropy s(
    min=-1.0e6,
    max=1.0e6,
    nominal=1.0e3) "specific entropy";
  SI.SpecificHeatCapacity cp(
    min=1.0,
    max=1.0e6,
    nominal=1000.0) "heat capacity at constant pressure";
  SI.SpecificHeatCapacity cv(
    min=1.0,
    max=1.0e6,
    nominal=1000.0) "heat capacity at constant volume";
  SI.SpecificHeatCapacity R(
    min=1.0,
    max=1.0e6,
    nominal=1000.0) "gas constant";
  SI.RatioOfSpecificHeatCapacities kappa "ratio of cp/cv";
  SI.Velocity a(
    min=1.0,
    max=10000.0,
    nominal=300.0) "speed of sound";
  SI.DerDensityByTemperature ddTp 
    "derivative of density by temperature at constant pressure";
  SI.DerDensityByPressure ddpT 
    "derivative of density by pressure at constant temperature";
  Real dupT(unit="m3.kg-1") 
    "derivative of inner energy by pressure at constant T";
  SI.SpecificHeatCapacity duTp 
    "derivative of inner energy by temperature at constant p";
end ThermoProperties_pT;

Modelica.Media.Common.ThermoFluidSpecial.ThermoProperties_dT

Information

Model description

A base class for medium property models which use density and temperature as dynamic states. This is a reasonable model for fluids that can be in the gas, liquid and two-phase region. The model is numerically not well suited for liquids except if the pressure is always above approx. 80% of the critical pressure.

Extends from Modelica.Icons.Record (Icon for records).

Modelica definition

record ThermoProperties_dT 
  "Thermodynamic property data for density d and temperature T as dynamic states"

  extends Modelica.Icons.Record;
  SI.Pressure p(
    min=1.0,
    max=1.0e9,
    nominal=1.0e5) "pressure";
  SI.SpecificEnthalpy h(
    min=-1.0e8,
    max=1.0e8,
    nominal=1.0e6) "specific enthalpy";
  SI.SpecificEnergy u(
    min=-1.0e8,
    max=1.0e8,
    nominal=1.0e6) "specific inner energy";
  SI.SpecificEntropy s(
    min=-1.0e6,
    max=1.0e6,
    nominal=1.0e3) "specific entropy";
  SI.SpecificHeatCapacity cp(
    min=1.0,
    max=1.0e6,
    nominal=1000.0) "heat capacity at constant pressure";
  SI.SpecificHeatCapacity cv(
    min=1.0,
    max=1.0e6,
    nominal=1000.0) "heat capacity at constant volume";
  SI.SpecificHeatCapacity R(
    min=1.0,
    max=1.0e6,
    nominal=1000.0) "gas constant";
  SI.RatioOfSpecificHeatCapacities kappa "ratio of cp/cv";
  SI.Velocity a(
    min=1.0,
    max=10000.0,
    nominal=300.0) "speed of sound";
  Real dudT(unit="m5/(kg.s2)") 
    "derivative of inner energy by density at constant T";
end ThermoProperties_dT;

Modelica.Media.Common.ThermoFluidSpecial.TransportProps

Information

Modelica definition

record TransportProps "record with transport properties"
  extends Modelica.Icons.Record;
  SI.DynamicViscosity eta;
  SI.ThermalConductivity lam;
end TransportProps;

Modelica.Media.Common.ThermoFluidSpecial.gibbsToProps_ph

Information

Inputs

Type	Name	Default	Description
GibbsDerivs	g		dimensionless derivatives of Gibbs function

Outputs

Type	Name	Description
ThermoProperties_ph	pro	property record for pressure and specific enthalpy as dynamic states

Modelica definition

function gibbsToProps_ph 
  "calulate property record for pressure and specific enthalpy as states from dimensionless Gibbs function"

  extends Modelica.Icons.Function;
  input GibbsDerivs g "dimensionless derivatives of Gibbs function";
  output ThermoProperties_ph pro 
    "property record for pressure and specific enthalpy as dynamic states";
protected 
  Real vt(unit="m3.kg-1.K-1") 
    "derivative of specific volume w.r.t. temperature";
  Real vp(unit="m4.kg-2.s2") "derivative of specific volume w.r.t. pressure";
algorithm 
  pro.T := g.T;
  pro.R := g.R;
  pro.d := g.p/(pro.R*pro.T*g.pi*g.gpi);
  pro.u := g.T*g.R*(g.tau*g.gtau - g.pi*g.gpi);
  pro.s := pro.R*(g.tau*g.gtau - g.g);
  pro.cp := -pro.R*g.tau*g.tau*g.gtautau;
  pro.cv := pro.R*(-g.tau*g.tau*g.gtautau + (g.gpi - g.tau*g.gtaupi)*(g.gpi
     - g.tau*g.gtaupi)/(g.gpipi));
  pro.a := abs(g.R*g.T*(g.gpi*g.gpi/((g.gpi - g.tau*g.gtaupi)*(g.gpi - g.
    tau*g.gtaupi)/(g.tau*g.tau*g.gtautau) - g.gpipi)))^0.5;
  vt := g.R/g.p*(g.pi*g.gpi - g.tau*g.pi*g.gtaupi);
  vp := g.R*g.T/(g.p*g.p)*g.pi*g.pi*g.gpipi;
  pro.kappa := -1/(pro.d*g.p)*pro.cp/(vp*pro.cp + vt*vt*g.T);
  pro.ddhp := -pro.d*pro.d*vt/(pro.cp);
  pro.ddph := -pro.d*pro.d*(vp*pro.cp - vt/pro.d + g.T*vt*vt)/pro.cp;
  pro.duph := -1/pro.d + g.p/(pro.d*pro.d)*pro.ddph;
  pro.duhp := 1 + g.p/(pro.d*pro.d)*pro.ddhp;
end gibbsToProps_ph;

Modelica.Media.Common.ThermoFluidSpecial.gibbsToBoundaryProps

Information

Inputs

Type	Name	Default	Description
GibbsDerivs	g		dimensionless derivatives of Gibbs function

Outputs

Type	Name	Description
PhaseBoundaryProperties	sat	phase boundary properties

Modelica definition

function gibbsToBoundaryProps 
  "calulate phase boundary property record from dimensionless Gibbs function"

  extends Modelica.Icons.Function;
  input GibbsDerivs g "dimensionless derivatives of Gibbs function";
  output PhaseBoundaryProperties sat "phase boundary properties";
protected 
  Real vt(unit="m3.kg-1.K-1") 
    "derivative of specific volume w.r.t. temperature";
  Real vp(unit="m4.kg-2.s2") "derivative of specific volume w.r.t. pressure";
algorithm 
  sat.d := g.p/(g.R*g.T*g.pi*g.gpi);
  sat.h := g.R*g.T*g.tau*g.gtau;
  sat.u := g.T*g.R*(g.tau*g.gtau - g.pi*g.gpi);
  sat.s := g.R*(g.tau*g.gtau - g.g);
  sat.cp := -g.R*g.tau*g.tau*g.gtautau;
  sat.cv := g.R*(-g.tau*g.tau*g.gtautau + (g.gpi - g.tau*g.gtaupi)*(g.gpi
     - g.tau*g.gtaupi)/(g.gpipi));
  vt := g.R/g.p*(g.pi*g.gpi - g.tau*g.pi*g.gtaupi);
  vp := g.R*g.T/(g.p*g.p)*g.pi*g.pi*g.gpipi;
  // sat.kappa := -1/(sat.d*g.p)*sat.cp/(vp*sat.cp + vt*vt*g.T);
  sat.pt := -g.p/g.T*(g.gpi - g.tau*g.gtaupi)/(g.gpipi*g.pi);
  sat.pd := -g.R*g.T*g.gpi*g.gpi/(g.gpipi);
end gibbsToBoundaryProps;

Modelica.Media.Common.ThermoFluidSpecial.gibbsToProps_dT

Information

Inputs

Type	Name	Default	Description
GibbsDerivs	g		dimensionless derivatives of Gibbs function

Outputs

Type	Name	Description
ThermoProperties_dT	pro	property record for density and temperature as dynamic states

Modelica definition

function gibbsToProps_dT 
  "calulate property record for density and temperature as states from dimensionless Gibbs function"

  extends Modelica.Icons.Function;
  input GibbsDerivs g "dimensionless derivatives of Gibbs function";
  output ThermoProperties_dT pro 
    "property record for density and temperature as dynamic states";
protected 
  Real vt(unit="m3.kg-1.K-1") 
    "derivative of specific volume w.r.t. temperature";
  Real vp(unit="m4.kg-2.s2") "derivative of specific volume w.r.t. pressure";
  Modelica.SIunits.Density d;
algorithm 
  pro.R := g.R;
  pro.p := g.p;
  pro.u := g.T*g.R*(g.tau*g.gtau - g.pi*g.gpi);
  pro.h := g.R*g.T*g.tau*g.gtau;
  pro.s := pro.R*(g.tau*g.gtau - g.g);
  pro.cp := -pro.R*g.tau*g.tau*g.gtautau;
  pro.cv := pro.R*(-g.tau*g.tau*g.gtautau + (g.gpi - g.tau*g.gtaupi)*(g.gpi
     - g.tau*g.gtaupi)/g.gpipi);
  vt := g.R/g.p*(g.pi*g.gpi - g.tau*g.pi*g.gtaupi);
  vp := g.R*g.T/(g.p*g.p)*g.pi*g.pi*g.gpipi;
  pro.kappa := -1/((g.p/(pro.R*g.T*g.pi*g.gpi))*g.p)*pro.cp/(vp*pro.cp + vt
    *vt*g.T);
  pro.a := abs(g.R*g.T*(g.gpi*g.gpi/((g.gpi - g.tau*g.gtaupi)*(g.gpi - g.
    tau*g.gtaupi)/(g.tau*g.tau*g.gtautau) - g.gpipi)))^0.5;

  d := g.p/(pro.R*g.T*g.pi*g.gpi);
  pro.dudT := (pro.p - g.T*vt/vp)/(d*d);
end gibbsToProps_dT;

Modelica.Media.Common.ThermoFluidSpecial.gibbsToProps_pT

Information

Inputs

Type	Name	Default	Description
GibbsDerivs	g		dimensionless derivatives of Gibbs function

Outputs

Type	Name	Description
ThermoProperties_pT	pro	property record for pressure and temperature as dynamic states

Modelica definition

function gibbsToProps_pT 
  "calulate property record for pressure and temperature as states from dimensionless Gibbs function"

  extends Modelica.Icons.Function;
  input GibbsDerivs g "dimensionless derivatives of Gibbs function";
  output ThermoProperties_pT pro 
    "property record for pressure and temperature as dynamic states";
protected 
  Real vt(unit="m3.kg-1.K-1") 
    "derivative of specific volume w.r.t. temperature";
  Real vp(unit="m4.kg-2.s2") "derivative of specific volume w.r.t. pressure";
algorithm 
  pro.R := g.R;
  pro.d := g.p/(pro.R*g.T*g.pi*g.gpi);
  pro.u := g.T*g.R*(g.tau*g.gtau - g.pi*g.gpi);
  pro.h := g.R*g.T*g.tau*g.gtau;
  pro.s := pro.R*(g.tau*g.gtau - g.g);
  pro.cp := -pro.R*g.tau*g.tau*g.gtautau;
  pro.cv := pro.R*(-g.tau*g.tau*g.gtautau + (g.gpi - g.tau*g.gtaupi)*(g.gpi
     - g.tau*g.gtaupi)/g.gpipi);
  vt := g.R/g.p*(g.pi*g.gpi - g.tau*g.pi*g.gtaupi);
  vp := g.R*g.T/(g.p*g.p)*g.pi*g.pi*g.gpipi;
  pro.kappa := -1/(pro.d*g.p)*pro.cp/(vp*pro.cp + vt*vt*g.T);
  pro.a := abs(g.R*g.T*(g.gpi*g.gpi/((g.gpi - g.tau*g.gtaupi)*(g.gpi - g.
    tau*g.gtaupi)/(g.tau*g.tau*g.gtautau) - g.gpipi)))^0.5;
  pro.ddpT := -(pro.d*pro.d)*vp;
  pro.ddTp := -(pro.d*pro.d)*vt;
  pro.duTp := pro.cp - g.p*vt;
  pro.dupT := -g.T*vt - g.p*vp;
end gibbsToProps_pT;

Modelica.Media.Common.ThermoFluidSpecial.helmholtzToProps_ph

Information

Inputs

Type	Name	Default	Description
HelmholtzDerivs	f		dimensionless derivatives of Helmholtz function

Outputs

Type	Name	Description
ThermoProperties_ph	pro	property record for pressure and specific enthalpy as dynamic states

Modelica definition

function helmholtzToProps_ph 
  "calulate property record for pressure and specific enthalpy as states from dimensionless Helmholtz function"

  extends Modelica.Icons.Function;
  input HelmholtzDerivs f "dimensionless derivatives of Helmholtz function";
  output ThermoProperties_ph pro 
    "property record for pressure and specific enthalpy as dynamic states";
protected 
  SI.Pressure p "pressure";
  DerPressureByDensity pd "derivative of pressure w.r.t. density";
  DerPressureByTemperature pt "derivative of pressure w.r.t. temperature";
  DerPressureBySpecificVolume pv 
    "derivative of pressure w.r.t. specific volume";
algorithm 
  pro.d := f.d;
  pro.T := f.T;
  pro.R := f.R;
  pro.s := f.R*(f.tau*f.ftau - f.f);
  pro.u := f.R*f.T*f.tau*f.ftau;
  p := pro.d*pro.R*pro.T*f.delta*f.fdelta;
  pd := f.R*f.T*f.delta*(2.0*f.fdelta + f.delta*f.fdeltadelta);
  pt := f.R*f.d*f.delta*(f.fdelta - f.tau*f.fdeltatau);
  pv := -pd*f.d*f.d;

  // calculating cp near the critical point may be troublesome (cp -> inf).
  pro.cp := f.R*(-f.tau*f.tau*f.ftautau + (f.delta*f.fdelta - f.delta*f.tau
    *f.fdeltatau)^2/(2*f.delta*f.fdelta + f.delta*f.delta*f.fdeltadelta));
  pro.cv := f.R*(-f.tau*f.tau*f.ftautau);
  pro.kappa := 1/(f.d*f.R*f.d*f.T*f.delta*f.fdelta)*((-pv*pro.cv + pt*pt*f.
    T)/(pro.cv));
  pro.a := abs(f.R*f.T*(2*f.delta*f.fdelta + f.delta*f.delta*f.fdeltadelta
     - ((f.delta*f.fdelta - f.delta*f.tau*f.fdeltatau)*(f.delta*f.fdelta -
    f.delta*f.tau*f.fdeltatau))/(f.tau*f.tau*f.ftautau)))^0.5;
  pro.ddph := (f.d*(pro.cv*f.d + pt))/(f.d*f.d*pd*pro.cv + f.T*pt*pt);
  pro.ddhp := -f.d*f.d*pt/(f.d*f.d*pd*pro.cv + f.T*pt*pt);
  pro.duph := -1/pro.d + p/(pro.d*pro.d)*pro.ddph;
  pro.duhp := 1 + p/(pro.d*pro.d)*pro.ddhp;
end helmholtzToProps_ph;

Modelica.Media.Common.ThermoFluidSpecial.helmholtzToProps_pT

Information

Inputs

Type	Name	Default	Description
HelmholtzDerivs	f		dimensionless derivatives of Helmholtz function

Outputs

Type	Name	Description
ThermoProperties_pT	pro	property record for pressure and temperature as dynamic states

Modelica definition

function helmholtzToProps_pT 
  "calulate property record for pressure and temperature as states from dimensionless Helmholtz function"

  extends Modelica.Icons.Function;
  input HelmholtzDerivs f "dimensionless derivatives of Helmholtz function";
  output ThermoProperties_pT pro 
    "property record for pressure and temperature as dynamic states";
protected 
  DerPressureByDensity pd "derivative of pressure w.r.t. density";
  DerPressureByTemperature pt "derivative of pressure w.r.t. temperature";
  DerPressureBySpecificVolume pv 
    "derivative of pressure w.r.t. specific volume";
  IsobaricVolumeExpansionCoefficient alpha 
    "isobaric volume expansion coefficient";
  // beta in Bejan
  IsothermalCompressibility gamma "isothermal compressibility";
  // kappa in Bejan
 SI.Pressure p "Pressure";
algorithm 
  pro.d := f.d;
  pro.R := f.R;
  pro.s := f.R*(f.tau*f.ftau - f.f);
  pro.h := f.R*f.T*(f.tau*f.ftau + f.delta*f.fdelta);
  pro.u := f.R*f.T*f.tau*f.ftau;
  pd := f.R*f.T*f.delta*(2.0*f.fdelta + f.delta*f.fdeltadelta);
  pt := f.R*f.d*f.delta*(f.fdelta - f.tau*f.fdeltatau);
  pv := -(f.d*f.d)*pd;
  alpha := -f.d*pt/pv;
  gamma := -f.d/pv;
  p := f.R*f.d*f.T*f.delta*f.fdelta;
  // calculating cp near the critical point may be troublesome (cp -> inf).
  pro.cp := f.R*(-f.tau*f.tau*f.ftautau + (f.delta*f.fdelta - f.delta*f.tau
    *f.fdeltatau)^2/(2*f.delta*f.fdelta + f.delta*f.delta*f.fdeltadelta));
  pro.cv := f.R*(-f.tau*f.tau*f.ftautau);
  pro.kappa := 1/(f.d*f.R*f.d*f.T*f.delta*f.fdelta)*((-pv*pro.cv + pt*pt*f.
    T)/(pro.cv));
  pro.a := abs(f.R*f.T*(2*f.delta*f.fdelta + f.delta*f.delta*f.fdeltadelta
     - ((f.delta*f.fdelta - f.delta*f.tau*f.fdeltatau)*(f.delta*f.fdelta -
    f.delta*f.tau*f.fdeltatau))/(f.tau*f.tau*f.ftautau)))^0.5;
  pro.ddTp := -pt/pd;
  pro.ddpT := 1/pd;
  //problem with units in last two lines
  pro.dupT := gamma*p/f.d - alpha*f.T/f.d;
  pro.duTp := pro.cp - alpha*p/f.d;
end helmholtzToProps_pT;

Modelica.Media.Common.ThermoFluidSpecial.helmholtzToProps_dT

Information

Inputs

Type	Name	Default	Description
HelmholtzDerivs	f		dimensionless derivatives of Helmholtz function

Outputs

Type	Name	Description
ThermoProperties_dT	pro	property record for density and temperature as dynamic states

Modelica definition

function helmholtzToProps_dT 
  "calulate property record for density and temperature as states from dimensionless Helmholtz function"

  extends Modelica.Icons.Function;
  input HelmholtzDerivs f "dimensionless derivatives of Helmholtz function";
  output ThermoProperties_dT pro 
    "property record for density and temperature as dynamic states";
protected 
  DerPressureByTemperature pt "derivative of pressure w.r.t. temperature";
  DerPressureBySpecificVolume pv "derivative of pressure w.r.t. pressure";
algorithm 
  pro.p := f.R*f.d*f.T*f.delta*f.fdelta;
  pro.R := f.R;
  pro.s := f.R*(f.tau*f.ftau - f.f);
  pro.h := f.R*f.T*(f.tau*f.ftau + f.delta*f.fdelta);
  pro.u := f.R*f.T*f.tau*f.ftau;
  pv := -(f.d*f.d)*f.R*f.T*f.delta*(2.0*f.fdelta + f.delta*f.fdeltadelta);
  pt := f.R*f.d*f.delta*(f.fdelta - f.tau*f.fdeltatau);

  // calculating cp near the critical point may be troublesome (cp -> inf).
  pro.cp := f.R*(-f.tau*f.tau*f.ftautau + (f.delta*f.fdelta - f.delta*f.tau
    *f.fdeltatau)^2/(2*f.delta*f.fdelta + f.delta*f.delta*f.fdeltadelta));
  pro.cv := f.R*(-f.tau*f.tau*f.ftautau);
  pro.kappa := 1/(f.d*pro.p)*((-pv*pro.cv + pt*pt*f.T)/(pro.cv));
  pro.a := abs(f.R*f.T*(2*f.delta*f.fdelta + f.delta*f.delta*f.fdeltadelta
     - ((f.delta*f.fdelta - f.delta*f.tau*f.fdeltatau)*(f.delta*f.fdelta -
    f.delta*f.tau*f.fdeltatau))/(f.tau*f.tau*f.ftautau)))^0.5;
  pro.dudT := (pro.p - f.T*pt)/(f.d*f.d);
end helmholtzToProps_dT;

Modelica.Media.Common.ThermoFluidSpecial.TwoPhaseToProps_ph

Information

Inputs

Type	Name	Default	Description
SaturationProperties	sat		saturation property record

Outputs

Type	Name	Description
ThermoProperties_ph	pro	property record for pressure and specific enthalpy as dynamic states

Modelica definition

function TwoPhaseToProps_ph 
  "compute property record for pressure and specific enthalpy as states from saturation properties"

  extends Modelica.Icons.Function;
  input SaturationProperties sat "saturation property record";
  output ThermoProperties_ph pro 
    "property record for pressure and specific enthalpy as dynamic states";
protected 
  Real dht(unit="(J/kg)/K") 
    "derivative of specific enthalpy w.r.t. temperature";
  Real dhd(unit="(J/kg)/(kg/m3)") 
    "derivative of specific enthalpy w.r.t. density";
  Real detph(unit="m4.s4/(K.s8)") "thermodynamic determinant";
algorithm 
  pro.d := sat.d;
  pro.T := sat.T;
  pro.u := sat.u;
  pro.s := sat.s;
  pro.cv := sat.cv;
  pro.R := sat.R;
  pro.cp := Modelica.Constants.inf;
  pro.kappa := -1/(sat.d*sat.p)*sat.dpT*sat.dpT*sat.T/sat.cv;
  pro.a := Modelica.Constants.inf;
  dht := sat.cv + sat.dpT/sat.d;
  dhd := -sat.T*sat.dpT/(sat.d*sat.d);
  detph := -sat.dpT*dhd;
  pro.ddph := dht/detph;
  pro.ddhp := -sat.dpT/detph;
end TwoPhaseToProps_ph;

Modelica.Media.Common.ThermoFluidSpecial.TwoPhaseToProps_dT

Information

Inputs

Type	Name	Default	Description
SaturationProperties	sat		saturation properties

Outputs

Type	Name	Description
ThermoProperties_dT	pro	property record for density and temperature as dynamic states

Modelica definition

function TwoPhaseToProps_dT 
  "compute property record for density and temperature as states from saturation properties"

  extends Modelica.Icons.Function;
  input SaturationProperties sat "saturation properties";
  output ThermoProperties_dT pro 
    "property record for density and temperature as dynamic states";
algorithm 
  pro.p := sat.p;
  pro.h := sat.h;
  pro.u := sat.u;
  pro.s := sat.s;
  pro.cv := sat.cv;
  pro.cp := Modelica.Constants.inf;
  pro.R := sat.R;
  pro.kappa := -1/(sat.d*sat.p)*sat.dpT*sat.dpT*sat.T/sat.cv;
  pro.a := Modelica.Constants.inf;
  pro.dudT := (sat.p - sat.T*sat.dpT)/(sat.d*sat.d);
end TwoPhaseToProps_dT;