def test_LiquidPhase():
m = ConcreteModel()
m.phase = LiquidPhase()
assert m.phase.is_liquid_phase()
assert not m.phase.is_solid_phase()
assert not 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 build(self):
'''
Callable method for Block construction.
'''
super().build()
self._state_block_class = WaterStateBlock
self.Liq = LiquidPhase()
self.H2O = Solvent()
# Get component set from database if provided
comp_set = None
if self.config.database is not None:
comp_set = self.config.database.get_solute_set(
self.config.water_source)
# Check definition of solute list
solute_list = self.config.solute_list
if solute_list is None:
# No user-provided solute list, look up list from database
if comp_set is None:
# No solute list in database and none provided.
raise ConfigurationError(
f"{self.name} no solute_list or database was defined. "
f"Users must provide at least one of these arguments.")
else:
solute_list = comp_set
elif self.config.database is not None:
# User provided custom list and database - check that all
# components are supported
for j in solute_list:
if j not in comp_set:
_log.info(f"{self.name} component {j} is not defined in "
f"the water_sources database file.")
else:
# User provided list but no database - assume they know what they
# are doing
pass
for j in solute_list:
self.add_component(str(j), Solute())
# Define default value for mass density of solution
self.dens_mass_default = 1000 * pyunits.kg / pyunits.m**3
# Define default value for dynamic viscosity of solution
self.visc_d_default = 0.001 * pyunits.kg / pyunits.m / pyunits.s
# ---------------------------------------------------------------------
# Set default scaling factors
self.default_scaling_factor = {
("flow_vol"): 1e3,
("conc_mass_comp"): 1e2
}
def build(self):
'''
Callable method for Block construction.
'''
super(PhysicalParameterData, self).build()
self._state_block_class = ThermalOilStateBlock
# Add Phase objects
self.Liq = LiquidPhase()
# Add Component objects
self.therminol66 = Component()
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 build(self):
"""
Callable method for Block construction.
"""
super(PropParameterData, self).build()
self._state_block_class = PropStateBlock
# phases
self.Liq = LiquidPhase()
# components
self.H2O = Solvent()
self.Na = Solute()
self.Ca = Solute()
self.Mg = Solute()
self.SO4 = Solute()
self.Cl = Solute()
# molecular weight
mw_comp_data = {
"H2O": 18.015e-3,
"Na": 22.990e-3,
"Ca": 40.078e-3,
"Mg": 24.305e-3,
"SO4": 96.06e-3,
"Cl": 35.453e-3,
}
self.mw_comp = Param(
self.component_list,
mutable=False,
initialize=extract_data(mw_comp_data),
units=pyunits.kg / pyunits.mol,
doc="Molecular weight",
)
self.dens_mass = Param(
mutable=False,
initialize=1000,
units=pyunits.kg / pyunits.m**3,
doc="Density",
)
self.cp = Param(mutable=False,
initialize=4.2e3,
units=pyunits.J / (pyunits.kg * pyunits.K))
# ---default scaling---
self.set_default_scaling("temperature", 1e-2)
self.set_default_scaling("pressure", 1e-6)
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(PropParameterData, self).build()
self._state_block_class = PropStateBlock
# phases
self.Liq = LiquidPhase()
# components
self.H2O = Solvent()
self.NaCl = Solute()
self.CaSO4 = Solute()
self.MgSO4 = Solute()
self.MgCl2 = Solute()
# molecular weight
mw_comp_data = {
'H2O': 18.015e-3,
'NaCl': 58.44e-3,
'CaSO4': 136.14e-3,
'MgSO4': 120.37e-3,
'MgCl2': 95.21e-3
}
self.mw_comp = Param(self.component_list,
mutable=False,
initialize=extract_data(mw_comp_data),
units=pyunits.kg / pyunits.mol,
doc="Molecular weight")
self.dens_mass = Param(mutable=False,
initialize=1000,
units=pyunits.kg / pyunits.m**3,
doc="Density")
self.cp = Param(mutable=False,
initialize=4.2e3,
units=pyunits.J / (pyunits.kg * pyunits.K))
# ---default scaling---
self.set_default_scaling('temperature', 1e-2)
self.set_default_scaling('pressure', 1e-6)
def build(self):
'''
Callable method for Block construction.
'''
super(PhysicalParameterData, self).build()
self._state_block_class = SaponificationStateBlock
# Add Phase objects
self.Liq = LiquidPhase()
# Add Component objects
self.H2O = Component()
self.NaOH = Component()
self.EthylAcetate = Component()
self.SodiumAcetate = Component()
self.Ethanol = Component()
# Heat capacity of water
self.cp_mol = Param(mutable=False,
initialize=75.327,
doc="Molar heat capacity of water [J/mol.K]",
units=units.J / units.mol / units.K)
# Density of water
self.dens_mol = Param(mutable=False,
initialize=55388.0,
doc="Molar density of water [mol/m^3]",
units=units.mol / units.m**3)
# Thermodynamic reference state
self.pressure_ref = Param(within=PositiveReals,
mutable=True,
default=101325.0,
doc='Reference pressure [Pa]',
units=units.Pa)
self.temperature_ref = Param(within=PositiveReals,
mutable=True,
default=298.15,
doc='Reference temperature [K]',
units=units.K)
def build(self):
'''
Callable method for Block construction.
'''
super(NaClParameterData, self).build()
self._state_block_class = NaClStateBlock
# components
self.H2O = Solvent()
self.NaCl = Solute()
# phases
self.Liq = LiquidPhase()
# reference
# this package is developed from Bartholomew & Mauter (2019) https://doi.org/10.1016/j.memsci.2018.11.067
# the enthalpy calculations are from Sharqawy et al. (2010) http://dx.doi.org/10.5004/dwt.2010.1079
# molecular weight
mw_comp_data = {'H2O': 18.01528E-3, 'NaCl': 58.44E-3}
self.mw_comp = Param(self.component_list,
mutable=False,
initialize=extract_data(mw_comp_data),
units=pyunits.kg / pyunits.mol,
doc="Molecular weight kg/mol")
# mass density parameters, eq 4 in Bartholomew
dens_mass_param_dict = {'0': 995, '1': 756}
self.dens_mass_param = Var(dens_mass_param_dict.keys(),
domain=Reals,
initialize=dens_mass_param_dict,
units=pyunits.kg / pyunits.m**3,
doc='Mass density parameters')
# dynamic viscosity parameters, eq 5 in Bartholomew
visc_d_param_dict = {'0': 9.80E-4, '1': 2.15E-3}
self.visc_d_param = Var(visc_d_param_dict.keys(),
domain=Reals,
initialize=visc_d_param_dict,
units=pyunits.Pa * pyunits.s,
doc='Dynamic viscosity parameters')
# diffusivity parameters, eq 6 in Bartholomew
diffus_param_dict = {
'0': 1.51e-9,
'1': -2.00e-9,
'2': 3.01e-8,
'3': -1.22e-7,
'4': 1.53e-7
}
self.diffus_param = Var(diffus_param_dict.keys(),
domain=Reals,
initialize=diffus_param_dict,
units=pyunits.m**2 / pyunits.s,
doc='Dynamic viscosity parameters')
# osmotic coefficient parameters, eq. 3b in Bartholomew
osm_coeff_param_dict = {'0': 0.918, '1': 8.89e-2, '2': 4.92}
self.osm_coeff_param = Var(osm_coeff_param_dict.keys(),
domain=Reals,
initialize=osm_coeff_param_dict,
units=pyunits.dimensionless,
doc='Osmotic coefficient parameters')
# TODO: update for NaCl solution, relationship from Sharqawy is for seawater
# specific enthalpy parameters, eq. 55 and 43 in Sharqawy (2010)
self.enth_mass_param_A1 = Var(within=Reals,
initialize=124.790,
units=pyunits.J / pyunits.kg,
doc='Specific enthalpy parameter A1')
self.enth_mass_param_A2 = Var(within=Reals,
initialize=4203.075,
units=(pyunits.J / pyunits.kg) *
pyunits.K**-1,
doc='Specific enthalpy parameter A2')
self.enth_mass_param_A3 = Var(within=Reals,
initialize=-0.552,
units=(pyunits.J / pyunits.kg) *
pyunits.K**-2,
doc='Specific enthalpy parameter A3')
self.enth_mass_param_A4 = Var(within=Reals,
initialize=0.004,
units=(pyunits.J / pyunits.kg) *
pyunits.K**-3,
doc='Specific enthalpy parameter A4')
self.enth_mass_param_B1 = Var(within=Reals,
initialize=27062.623,
units=pyunits.dimensionless,
doc='Specific enthalpy parameter B1')
self.enth_mass_param_B2 = Var(within=Reals,
initialize=4835.675,
units=pyunits.dimensionless,
doc='Specific enthalpy parameter B2')
# 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-6)
self.set_default_scaling('dens_mass_phase', 1e-3, index='Liq')
self.set_default_scaling('visc_d_phase', 1e3, index='Liq')
self.set_default_scaling('diffus_phase', 1e9, index='Liq')
self.set_default_scaling('osm_coeff', 1e0)
self.set_default_scaling('enth_mass_phase', 1e-5, index='Liq')
def build(self):
"""
Callable method for Block construction.
"""
super(CoagulationParameterData, self).build()
self._state_block_class = CoagulationStateBlock
# phases
self.Liq = LiquidPhase()
# components
self.H2O = Component()
self.TSS = Component()
self.TDS = Component()
self.Sludge = Component()
# heat capacity of liquid
self.cp = Param(mutable=False,
initialize=4184,
units=pyunits.J / (pyunits.kg * pyunits.K))
# reference density of liquid
self.ref_dens_liq = Param(
domain=Reals,
initialize=999.26,
mutable=True,
units=pyunits.kg / pyunits.m**3,
doc='Reference water mass density parameter @ 0 oC and no salts')
# change in liquid density with increasing mass fraction of salts/solids
self.dens_slope = Param(
domain=Reals,
initialize=879.04,
mutable=True,
units=pyunits.kg / pyunits.m**3,
doc=
'Relative increase in liquid density with mass fraction of salts')
# adjustment parameters for density change with temperature
# Density calculation as a function of temperature and pressure
# --------------------------------------------------------------
# Engineering Toolbox. Water - Density, Specific Weight, and
# Thermal Expansion Coefficients. (2003) https://www.engineeringtoolbox.com/
# water-density-specific-weight-d_595.html [Accessed 02-01-2022]
self.dens_param_A = Param(
domain=Reals,
initialize=-2.9335E-6,
mutable=True,
units=pyunits.K**-2,
doc='Density correction parameter A for temperature variation')
self.dens_param_B = Param(
domain=Reals,
initialize=0.001529811,
mutable=True,
units=pyunits.K**-1,
doc='Density correction parameter B for temperature variation')
self.dens_param_C = Param(
domain=Reals,
initialize=0.787973,
mutable=True,
units=pyunits.dimensionless,
doc='Density correction parameter C for temperature variation')
# Correction factors for changes in density with changes in pressure
self.ref_pressure_correction = Param(
domain=Reals,
initialize=1.0135,
mutable=True,
units=pyunits.dimensionless,
doc='Density reference correction parameter for changes in pressure'
)
self.ref_pressure_slope = Param(
domain=Reals,
initialize=4.9582E-10,
mutable=True,
units=pyunits.Pa**-1,
doc='Slope of density change as a function of pressure')
# reference density of solids
self.ref_dens_sol = Param(
domain=Reals,
initialize=2600,
mutable=True,
units=pyunits.kg / pyunits.m**3,
doc='Reference dry solid mass density parameter')
# ---default scaling---
self.set_default_scaling('temperature', 1e-2)
self.set_default_scaling('pressure', 1e-6)
def build(self):
'''
Callable method for Block construction.
'''
super(DSPMDEParameterData, self).build()
self._state_block_class = DSPMDEStateBlock
# components
self.H2O = Solvent()
for j in self.config.solute_list:
self.add_component(str(j), Solute())
# phases
self.Liq = LiquidPhase()
# reference
# Todo: enter any relevant references
# TODO: consider turning parameters into variables for future param estimation
# molecular weight
self.mw_comp = Param(
self.component_list,
mutable=True,
default=18e-3,
initialize=self.config.mw_data,
units=pyunits.kg/pyunits.mol,
doc="Molecular weight")
# Stokes radius
self.radius_stokes_comp = Param(
self.solute_set,
mutable=True,
default=1e-10,
initialize=self.config.stokes_radius_data,
units=pyunits.m,
doc="Stokes radius of solute")
self.diffus_phase_comp = Param(
self.phase_list,
self.solute_set,
mutable=True,
default=1e-9,
initialize=self.config.diffusivity_data,
units=pyunits.m ** 2 * pyunits.s ** -1,
doc="Bulk diffusivity of ion")
self.visc_d_phase = Param(
self.phase_list,
mutable=True,
default=1e-3,
initialize=1e-3, #TODO:revisit- assuming ~ 1e-3 Pa*s for pure water
units=pyunits.Pa * pyunits.s,
doc="Fluid viscosity")
# Ion charge
self.charge_comp = Param(
self.solute_set,
mutable=True,
default=1,
initialize=self.config.charge,
units=pyunits.dimensionless,
doc="Ion charge")
# Dielectric constant of water
self.dielectric_constant = Param(
mutable=True,
default=80.4,
initialize=80.4, #todo: make a variable with parameter values for coefficients in the function of temperature
units=pyunits.dimensionless,
doc="Dielectric constant of water")
# ---default scaling---
self.set_default_scaling('temperature', 1e-2)
self.set_default_scaling('pressure', 1e-6)
self.set_default_scaling('dens_mass_phase', 1e-3, index='Liq')
self.set_default_scaling('visc_d_phase', 1e3, index='Liq')
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)
def build(self):
"""
Callable method for Block construction.
"""
super(CoagulationParameterData, self).build()
self._state_block_class = CoagulationStateBlock
# phases
self.Liq = LiquidPhase()
# components
self.H2O = Component()
self.TSS = Component()
self.TDS = Component()
self.Sludge = Component()
# heat capacity of liquid
self.cp = Param(mutable=False,
initialize=4184,
units=pyunits.J / (pyunits.kg * pyunits.K))
# reference density of liquid
self.ref_dens_liq = Param(
domain=Reals,
initialize=999.26,
mutable=True,
units=pyunits.kg / pyunits.m**3,
doc="Reference water mass density parameter @ 0 oC and no salts",
)
# change in liquid density with increasing mass fraction of salts/solids
self.dens_slope = Param(
domain=Reals,
initialize=879.04,
mutable=True,
units=pyunits.kg / pyunits.m**3,
doc=
"Relative increase in liquid density with mass fraction of salts",
)
# adjustment parameters for density change with temperature
# Density calculation as a function of temperature and pressure
# --------------------------------------------------------------
# Engineering Toolbox. Water - Density, Specific Weight, and
# Thermal Expansion Coefficients. (2003) https://www.engineeringtoolbox.com/
# water-density-specific-weight-d_595.html [Accessed 02-01-2022]
self.dens_param_A = Param(
domain=Reals,
initialize=-2.9335e-6,
mutable=True,
units=pyunits.K**-2,
doc="Density correction parameter A for temperature variation",
)
self.dens_param_B = Param(
domain=Reals,
initialize=0.001529811,
mutable=True,
units=pyunits.K**-1,
doc="Density correction parameter B for temperature variation",
)
self.dens_param_C = Param(
domain=Reals,
initialize=0.787973,
mutable=True,
units=pyunits.dimensionless,
doc="Density correction parameter C for temperature variation",
)
# Correction factors for changes in density with changes in pressure
self.ref_pressure_correction = Param(
domain=Reals,
initialize=1.0135,
mutable=True,
units=pyunits.dimensionless,
doc=
"Density reference correction parameter for changes in pressure",
)
self.ref_pressure_slope = Param(
domain=Reals,
initialize=4.9582e-10,
mutable=True,
units=pyunits.Pa**-1,
doc="Slope of density change as a function of pressure",
)
# Adjustment for dynamic viscosity as a function of temperature
# -------------------------------------------------------------
# D.S. Viswananth, G. Natarajan. Data Book on the Viscosity of
# Liquids. Hemisphere Publishing Corp. (1989)
self.mu_A = Param(
domain=Reals,
initialize=2.939e-5,
mutable=True,
units=pyunits.kg / pyunits.m / pyunits.s,
doc="Pre-exponential factor for viscosity calculation",
)
self.mu_B = Param(
domain=Reals,
initialize=507.88,
mutable=True,
units=pyunits.K,
doc="Exponential numerator term for viscosity calculation",
)
self.mu_C = Param(
domain=Reals,
initialize=149.3,
mutable=True,
units=pyunits.K,
doc="Exponential denominator term for viscosity calculation",
)
# ---default scaling---
self.set_default_scaling("temperature", 1e-2)
self.set_default_scaling("pressure", 1e-6)