diff --git a/ThermofluidStream/Media/myMedia/Incompressible.mo b/ThermofluidStream/Media/myMedia/Incompressible.mo index db785c26..18c5394d 100644 --- a/ThermofluidStream/Media/myMedia/Incompressible.mo +++ b/ThermofluidStream/Media/myMedia/Incompressible.mo @@ -392,6 +392,21 @@ which is only exactly true for a fluid with constant density d=d0. annotation(smoothOrder=2); end thermalConductivity; + redeclare function extends isobaricExpansionCoefficient + "Return isobaric expansion coefficient (beta) as a function of the thermodynamic state record" + algorithm + beta := -density_derT_p(state) / density(state); + annotation( + Documentation(info = "

The isobaric expansion coefficient beta is defined as

+
1/v * der(v,T)
+
+

with v = 1/d, at constant pressure p. +Using the chain rule:

+
1/v * der(v,T) = d * (-der(d, T) / d^2)
+= -der(d, T) / d

+")); + end isobaricExpansionCoefficient; + function s_T "Compute specific entropy" extends Modelica.Icons.Function; input Temperature T "Temperature"; @@ -498,6 +513,18 @@ which is only exactly true for a fluid with constant density d=d0. ddph :=0; //incompressable and h = cp*T; h=const; -> T=const end density_derp_h; + redeclare function extends density_derh_p + "Return density derivative w.r.t. specific enthalpy at constant pressure" + algorithm + ddhp := density_derT_p(state) / specificHeatCapacityCp(state); + end density_derh_p; + + redeclare function extends density_derT_p + "Return density derivative w.r.t. temperature at constant pressure" + algorithm + ddTp := Polynomials.derivativeValue(poly_rho, if TinK then state.T else Cv.to_degC(state.T)); + end density_derT_p; + redeclare function extends specificEnthalpy "Return specific enthalpy as a function of the thermodynamic state record" algorithm