def fill(self, fase, estado):
fase._bool = True
fase.M = self.M
fase.rho = unidades.Density(estado.rho)
fase.v = unidades.SpecificVolume(estado.v)
fase.Z = unidades.Dimensionless(self.P * estado.v / R / 1000 * self.M /
self.T)
fase.h = unidades.Enthalpy(estado.h)
fase.s = unidades.SpecificHeat(estado.s)
fase.u = unidades.Enthalpy(estado.u)
fase.a = unidades.Enthalpy(fase.u - self.T * fase.s)
fase.g = unidades.Enthalpy(fase.h - self.T * fase.s)
fi = exp((fase.g - self.g0) / 1000 / R * self.M / self.T)
fase.fi = [unidades.Pressure(fi)]
fase.f = [unidades.Pressure(self.P * f) for f in fase.fi]
fase.cv = unidades.SpecificHeat(estado.cv)
fase.cp = unidades.SpecificHeat(estado.cp)
fase.cp_cv = unidades.Dimensionless(fase.cp / fase.cv)
fase.gamma = fase.cp_cv
fase.w = unidades.Speed(estado.w)
fase.rhoM = unidades.MolarDensity(fase.rho / self.M)
fase.hM = unidades.MolarEnthalpy(fase.h * self.M)
fase.sM = unidades.MolarSpecificHeat(fase.s * self.M)
fase.uM = unidades.MolarEnthalpy(fase.u * self.M)
fase.aM = unidades.MolarEnthalpy(fase.a * self.M)
fase.gM = unidades.MolarEnthalpy(fase.g * self.M)
fase.cvM = unidades.MolarSpecificHeat(fase.cv * self.M)
fase.cpM = unidades.MolarSpecificHeat(fase.cp * self.M)
fase.mu = unidades.Viscosity(estado.mu)
fase.nu = unidades.Diffusivity(fase.mu / fase.rho)
fase.k = unidades.ThermalConductivity(estado.k)
fase.alfa = unidades.Diffusivity(fase.k / fase.rho / fase.cp)
fase.epsilon = unidades.Dimensionless(
iapws._Dielectric(estado.rho, self.T))
fase.Prandt = unidades.Dimensionless(estado.mu * estado.cp / estado.k)
fase.n = unidades.Dimensionless(
iapws._Refractive(fase.rho, self.T, self.kwargs["l"]))
fase.alfav = unidades.InvTemperature(estado.deriv("Tpv") / fase.v)
fase.kappa = unidades.InvPressure(-estado.deriv("pTv") / fase.v)
fase.kappas = unidades.InvPressure(
-1 / fase.v * self.derivative("v", "P", "s", fase))
fase.joule = unidades.TemperaturePressure(estado.deriv("Tph"))
fase.deltat = unidades.EnthalpyPressure(estado.deriv("pTh"))
fase.alfap = unidades.Density(fase.alfav / self.P / fase.kappa)
fase.betap = unidades.Density(-1 / self.P * estado.deriv("vTp"))
fase.fraccion = [1]
fase.fraccion_masica = [1]
def fill(self, fase, estado):
fase._bool = True
fase.M = self.M
fase.rho = unidades.Density(estado.rho)
fase.v = unidades.SpecificVolume(estado.v)
fase.Z = unidades.Dimensionless(self.P*estado.v/R/1000*self.M/self.T)
fase.h = unidades.Enthalpy(estado.h)
fase.s = unidades.SpecificHeat(estado.s)
fase.u = unidades.Enthalpy(estado.u)
fase.a = unidades.Enthalpy(fase.u-self.T*fase.s)
fase.g = unidades.Enthalpy(fase.h-self.T*fase.s)
fi = exp((fase.g-self.g0)/1000/R*self.M/self.T)
fase.fi = [unidades.Pressure(fi)]
fase.f = [unidades.Pressure(self.P*f) for f in fase.fi]
fase.cv = unidades.SpecificHeat(estado.cv)
fase.cp = unidades.SpecificHeat(estado.cp)
fase.cp_cv = unidades.Dimensionless(fase.cp/fase.cv)
fase.gamma = fase.cp_cv
fase.w = unidades.Speed(estado.w)
fase.rhoM = unidades.MolarDensity(fase.rho/self.M)
fase.hM = unidades.MolarEnthalpy(fase.h*self.M)
fase.sM = unidades.MolarSpecificHeat(fase.s*self.M)
fase.uM = unidades.MolarEnthalpy(fase.u*self.M)
fase.aM = unidades.MolarEnthalpy(fase.a*self.M)
fase.gM = unidades.MolarEnthalpy(fase.g*self.M)
fase.cvM = unidades.MolarSpecificHeat(fase.cv*self.M)
fase.cpM = unidades.MolarSpecificHeat(fase.cp*self.M)
fase.mu = unidades.Viscosity(estado.mu)
fase.nu = unidades.Diffusivity(fase.mu/fase.rho)
fase.k = unidades.ThermalConductivity(estado.k)
fase.alfa = unidades.Diffusivity(fase.k/fase.rho/fase.cp)
fase.epsilon = unidades.Dimensionless(
iapws._Dielectric(estado.rho, self.T))
fase.Prandt = unidades.Dimensionless(estado.mu*estado.cp/estado.k)
fase.n = unidades.Dimensionless(
iapws._Refractive(fase.rho, self.T, self.kwargs["l"]))
fase.alfav = unidades.InvTemperature(estado.deriv("Tpv")/fase.v)
fase.kappa = unidades.InvPressure(-estado.deriv("pTv")/fase.v)
fase.kappas = unidades.InvPressure(
-1/fase.v*self.derivative("v", "P", "s", fase))
fase.joule = unidades.TemperaturePressure(estado.deriv("Tph"))
fase.deltat = unidades.EnthalpyPressure(estado.deriv("pTh"))
fase.alfap = unidades.Density(fase.alfav/self.P/fase.kappa)
fase.betap = unidades.Density(-1/self.P*estado.deriv("vTp"))
fase.fraccion = [1]
fase.fraccion_masica = [1]
def _Dielectric(self, rho, T):
try:
nu = _Dielectric(rho, T)
except NotImplementedError:
nu = None
return unidades.Dimensionless(nu)
def _Dielectric(self, rho, T):
return unidades.Dimensionless(_Dielectric(rho, T))
def _Dielectric(self, rho, T):
try:
nu = _Dielectric(rho, T)
except NotImplementedError:
nu = None
return unidades.Dimensionless(nu)
