def test_VaporPhase():
m = ConcreteModel()
m.phase = VaporPhase()
assert not m.phase.is_liquid_phase()
assert not m.phase.is_solid_phase()
assert m.phase.is_vapor_phase()
def test_is_phase_valid_liquid(self, m):
m.Liq = LiquidPhase()
m.Sol = SolidPhase()
m.Vap = VaporPhase()
m.Phase = Phase()
assert m.comp._is_phase_valid(m.Liq)
assert not m.comp._is_phase_valid(m.Sol)
assert not m.comp._is_phase_valid(m.Vap)
assert not m.comp._is_phase_valid(m.Phase)
def test_is_phase_valid_liquid(self, m):
m.comp3 = Component(
default={"valid_phase_types": PhaseType.liquidPhase})
m.Liq = LiquidPhase()
m.Sol = SolidPhase()
m.Vap = VaporPhase()
m.Phase = Phase()
assert m.comp3._is_phase_valid(m.Liq)
assert not m.comp3._is_phase_valid(m.Sol)
assert not m.comp3._is_phase_valid(m.Vap)
assert not m.comp3._is_phase_valid(m.Phase)
def test_is_phase_valid_liquid(self, m):
m.comp3 = Apparent(default={
"valid_phase_types": PhaseType.liquidPhase,
"dissociation_species": {"comp": 1},
"_electrolyte": True})
m.Liq = LiquidPhase()
m.Sol = SolidPhase()
m.Vap = VaporPhase()
m.Aqu = AqueousPhase()
m.Phase = Phase()
assert m.comp3._is_phase_valid(m.Liq)
assert not m.comp3._is_phase_valid(m.Sol)
assert not m.comp3._is_phase_valid(m.Vap)
assert not m.comp3._is_phase_valid(m.Aqu)
assert not m.comp3._is_phase_valid(m.Phase)
def build(self):
"""
Callable method for Block construction.
"""
super(WaterParameterData, self).build()
self._state_block_class = WaterStateBlock
# components
self.H2O = Component()
# phases
self.Liq = LiquidPhase()
self.Vap = VaporPhase()
""" References
This package was developed from the following references:
- K.G.Nayar, M.H.Sharqawy, L.D.Banchik, and J.H.Lienhard V, "Thermophysical properties of seawater: A review and
new correlations that include pressure dependence,"Desalination, Vol.390, pp.1 - 24, 2016.
doi: 10.1016/j.desal.2016.02.024(preprint)
- Mostafa H.Sharqawy, John H.Lienhard V, and Syed M.Zubair, "Thermophysical properties of seawater: A review of
existing correlations and data,"Desalination and Water Treatment, Vol.16, pp.354 - 380, April 2010.
(2017 corrections provided at http://web.mit.edu/seawater)
"""
# Parameters
# molecular weight
self.mw_comp = Param(
self.component_list,
mutable=False,
initialize=18.01528e-3,
units=pyunits.kg / pyunits.mol,
doc="Molecular weight",
)
# Liq mass density parameters, eq. 8 in Sharqawy et al. (2010)
dens_units = pyunits.kg / pyunits.m**3
t_inv_units = pyunits.K**-1
s_inv_units = pyunits.kg / pyunits.g
self.dens_mass_param_A1 = Var(
within=Reals,
initialize=9.999e2,
units=dens_units,
doc="Mass density parameter A1",
)
self.dens_mass_param_A2 = Var(
within=Reals,
initialize=2.034e-2,
units=dens_units * t_inv_units,
doc="Mass density parameter A2",
)
self.dens_mass_param_A3 = Var(
within=Reals,
initialize=-6.162e-3,
units=dens_units * t_inv_units**2,
doc="Mass density parameter A3",
)
self.dens_mass_param_A4 = Var(
within=Reals,
initialize=2.261e-5,
units=dens_units * t_inv_units**3,
doc="Mass density parameter A4",
)
self.dens_mass_param_A5 = Var(
within=Reals,
initialize=-4.657e-8,
units=dens_units * t_inv_units**4,
doc="Mass density parameter A5",
)
# Vap mass density parameters (approximating using ideal gas)
self.dens_mass_param_mw = Var(
within=Reals,
initialize=18.01528e-3,
units=pyunits.kg / pyunits.mol,
doc="Mass density parameter molecular weight",
)
self.dens_mass_param_R = Var(
within=Reals,
initialize=8.31462618,
units=pyunits.J / pyunits.mol / pyunits.K,
doc="Mass density parameter universal gas constant",
)
# vapor pressure parameters, eq. 5 and 6 in Nayar et al.(2016)
self.pressure_sat_param_psatw_A1 = Var(
within=Reals,
initialize=-5.8002206e3,
units=pyunits.K,
doc="Vapor pressure of pure water parameter A1",
)
self.pressure_sat_param_psatw_A2 = Var(
within=Reals,
initialize=1.3914993,
units=pyunits.dimensionless,
doc="Vapor pressure of pure water parameter A2",
)
self.pressure_sat_param_psatw_A3 = Var(
within=Reals,
initialize=-4.8640239e-2,
units=t_inv_units,
doc="Vapor pressure of pure water parameter A3",
)
self.pressure_sat_param_psatw_A4 = Var(
within=Reals,
initialize=4.1764768e-5,
units=t_inv_units**2,
doc="Vapor pressure of pure water parameter A4",
)
self.pressure_sat_param_psatw_A5 = Var(
within=Reals,
initialize=-1.4452093e-8,
units=t_inv_units**3,
doc="Vapor pressure of pure water parameter A5",
)
self.pressure_sat_param_psatw_A6 = Var(
within=Reals,
initialize=6.5459673,
units=pyunits.dimensionless,
doc="Vapor pressure of pure water parameter A6",
)
# specific enthalpy parameters, eq. 55 and 43 in Sharqawy et al. (2010)
enth_mass_units = pyunits.J / pyunits.kg
self.enth_mass_param_A1 = Var(
within=Reals,
initialize=141.355,
units=enth_mass_units,
doc="Specific enthalpy parameter A1",
)
self.enth_mass_param_A2 = Var(
within=Reals,
initialize=4202.07,
units=enth_mass_units * t_inv_units,
doc="Specific enthalpy parameter A2",
)
self.enth_mass_param_A3 = Var(
within=Reals,
initialize=-0.535,
units=enth_mass_units * t_inv_units**2,
doc="Specific enthalpy parameter A3",
)
self.enth_mass_param_A4 = Var(
within=Reals,
initialize=0.004,
units=enth_mass_units * t_inv_units**3,
doc="Specific enthalpy parameter A4",
)
# self.enth_mass_param_B1 = Var(
# within=Reals, initialize=-2.348e4, units=enth_mass_units,
# doc='Specific enthalpy parameter B1')
# self.enth_mass_param_B2 = Var(
# within=Reals, initialize=3.152e5, units=enth_mass_units,
# doc='Specific enthalpy parameter B2')
# self.enth_mass_param_B3 = Var(
# within=Reals, initialize=2.803e6, units=enth_mass_units,
# doc='Specific enthalpy parameter B3')
# self.enth_mass_param_B4 = Var(
# within=Reals, initialize=-1.446e7, units=enth_mass_units,
# doc='Specific enthalpy parameter B4')
# self.enth_mass_param_B5 = Var(
# within=Reals, initialize=7.826e3, units=enth_mass_units * t_inv_units,
# doc='Specific enthalpy parameter B5')
# self.enth_mass_param_B6 = Var(
# within=Reals, initialize=-4.417e1, units=enth_mass_units * t_inv_units**2,
# doc='Specific enthalpy parameter B6')
# self.enth_mass_param_B7 = Var(
# within=Reals, initialize=2.139e-1, units=enth_mass_units * t_inv_units**3,
# doc='Specific enthalpy parameter B7')
# self.enth_mass_param_B8 = Var(
# within=Reals, initialize=-1.991e4, units=enth_mass_units * t_inv_units,
# doc='Specific enthalpy parameter B8')
# self.enth_mass_param_B9 = Var(
# within=Reals, initialize=2.778e4, units=enth_mass_units * t_inv_units,
# doc='Specific enthalpy parameter B9')
# self.enth_mass_param_B10 = Var(
# within=Reals, initialize=9.728e1, units=enth_mass_units * t_inv_units**2,
# doc='Specific enthalpy parameter B10')
# specific heat parameters from eq (9) in Sharqawy et al. (2010)
cp_units = pyunits.J / (pyunits.kg * pyunits.K)
self.cp_phase_param_A1 = Var(
within=Reals,
initialize=5.328,
units=cp_units,
doc="Specific heat of seawater parameter A1",
)
# self.cp_phase_param_A2 = Var(
# within=Reals, initialize=-9.76e-2, units=cp_units * s_inv_units,
# doc='Specific heat of seawater parameter A2')
# self.cp_phase_param_A3 = Var(
# within=Reals, initialize=4.04e-4, units=cp_units * s_inv_units**2,
# doc='Specific heat of seawater parameter A3')
self.cp_phase_param_B1 = Var(
within=Reals,
initialize=-6.913e-3,
units=cp_units * t_inv_units,
doc="Specific heat of seawater parameter B1",
)
# self.cp_phase_param_B2 = Var(
# within=Reals, initialize=7.351e-4, units=cp_units * s_inv_units * t_inv_units,
# doc='Specific heat of seawater parameter B2')
# self.cp_phase_param_B3 = Var(
# within=Reals, initialize=-3.15e-6, units=cp_units * s_inv_units**2 * t_inv_units,
# doc='Specific heat of seawater parameter B3')
self.cp_phase_param_C1 = Var(
within=Reals,
initialize=9.6e-6,
units=cp_units * t_inv_units**2,
doc="Specific heat of seawater parameter C1",
)
# self.cp_phase_param_C2 = Var(
# within=Reals, initialize=-1.927e-6, units=cp_units * s_inv_units * t_inv_units**2,
# doc='Specific heat of seawater parameter C2')
# self.cp_phase_param_C3 = Var(
# within=Reals, initialize=8.23e-9, units=cp_units * s_inv_units**2 * t_inv_units**2,
# doc='Specific heat of seawater parameter C3')
self.cp_phase_param_D1 = Var(
within=Reals,
initialize=2.5e-9,
units=cp_units * t_inv_units**3,
doc="Specific heat of seawater parameter D1",
)
# self.cp_phase_param_D2 = Var(
# within=Reals, initialize=1.666e-9, units=cp_units * s_inv_units * t_inv_units**3,
# doc='Specific heat of seawater parameter D2')
# self.cp_phase_param_D3 = Var(
# within=Reals, initialize=-7.125e-12, units=cp_units * s_inv_units**2 * t_inv_units**3,
# doc='Specific heat of seawater parameter D3')
# Specific heat parameters for Cp vapor from NIST Webbook
# Chase, M.W., Jr., NIST-JANAF Themochemical Tables, Fourth Edition, J. Phys. Chem. Ref. Data, Monograph 9, 1998, 1-1951
self.cp_vap_param_A = Var(
within=Reals,
initialize=30.09200 / 18.01528e-3,
units=cp_units,
doc="Specific heat of water vapor parameter A",
)
self.cp_vap_param_B = Var(
within=Reals,
initialize=6.832514 / 18.01528e-3,
units=cp_units * t_inv_units,
doc="Specific heat of water vapor parameter B",
)
self.cp_vap_param_C = Var(
within=Reals,
initialize=6.793435 / 18.01528e-3,
units=cp_units * t_inv_units**2,
doc="Specific heat of water vapor parameter C",
)
self.cp_vap_param_D = Var(
within=Reals,
initialize=-2.534480 / 18.01528e-3,
units=cp_units * t_inv_units**3,
doc="Specific heat of water vapor parameter D",
)
self.cp_vap_param_E = Var(
within=Reals,
initialize=0.082139 / 18.01528e-3,
units=cp_units * t_inv_units**-2,
doc="Specific heat of water vapor parameter E",
)
# latent heat of pure water parameters from eq. 54 in Sharqawy et al. (2010)
self.dh_vap_w_param_0 = Var(
within=Reals,
initialize=2.501e6,
units=enth_mass_units,
doc="Latent heat of pure water parameter 0",
)
self.dh_vap_w_param_1 = Var(
within=Reals,
initialize=-2.369e3,
units=cp_units,
doc="Latent heat of pure water parameter 1",
)
self.dh_vap_w_param_2 = Var(
within=Reals,
initialize=2.678e-1,
units=enth_mass_units * t_inv_units**2,
doc="Latent heat of pure water parameter 2",
)
self.dh_vap_w_param_3 = Var(
within=Reals,
initialize=-8.103e-3,
units=enth_mass_units * t_inv_units**3,
doc="Latent heat of pure water parameter 3",
)
self.dh_vap_w_param_4 = Var(
within=Reals,
initialize=-2.079e-5,
units=enth_mass_units * t_inv_units**4,
doc="Latent heat of pure water parameter 4",
)
# traditional parameters are the only Vars currently on the block and should be fixed
for v in self.component_objects(Var):
v.fix()
# ---default scaling---
self.set_default_scaling("temperature", 1e-2)
self.set_default_scaling("pressure", 1e-5)
self.set_default_scaling("dens_mass_phase", 1e-3, index="Liq")
self.set_default_scaling("dens_mass_phase", 1, index="Vap")
# self.set_default_scaling('dens_mass_solvent', 1e-3)
self.set_default_scaling("enth_mass_phase", 1e-5, index="Liq")
self.set_default_scaling("enth_mass_phase", 1e-6, index="Vap")
self.set_default_scaling("pressure_sat", 1e-5)
self.set_default_scaling("cp_phase", 1e-3, index="Liq")
self.set_default_scaling("cp_phase", 1e-3, index="Vap")
self.set_default_scaling("dh_vap", 1e-6)