#################################################################################
# The Institute for the Design of Advanced Energy Systems Integrated Platform
# Framework (IDAES IP) was produced under the DOE Institute for the
# Design of Advanced Energy Systems (IDAES), and is copyright (c) 2018-2021
# by the software owners: The Regents of the University of California, through
# Lawrence Berkeley National Laboratory, National Technology & Engineering
# Solutions of Sandia, LLC, Carnegie Mellon University, West Virginia University
# Research Corporation, et al. All rights reserved.
#
# Please see the files COPYRIGHT.md and LICENSE.md for full copyright and
# license information.
#################################################################################
__author__ = "John Eslick"
# Import IDAES
from idaes.core import declare_process_block_class
import idaes.logger as idaeslog
from idaes.models.properties.general_helmholtz import (
helmholtz_available,
HelmholtzParameterBlockData,
HelmholtzParameterBlock,
HelmholtzStateBlockData,
HelmholtzThermoExpressions,
PhaseType,
StateVars,
AmountBasis,
)
from idaes.models.properties.general_helmholtz.helmholtz_state import _StateBlock
# Logger
_log = idaeslog.getLogger(__name__)
def iapws95_available():
"""Make sure the compiled IAPWS-95 functions are available. Yes, in Windows
the .so extension is still used.
"""
return helmholtz_available()
[docs]def htpx(T=None, P=None, x=None):
"""
Convenience function to calculate steam enthalpy from temperature and
either pressure or vapor fraction. This function can be used for inlet
streams and initialization where temperature is known instead of enthalpy.
User must provide values for two of T, P, or x.
Args:
T: Temperature with units (between 200 and 3000 K)
P: Pressure with units (between 1 and 1e9 Pa), None if saturated vapor
x: Vapor fraction [mol vapor/mol total] (between 0 and 1), None if
superheated or subcooled
Returns:
Total molar enthalpy [J/mol].
"""
prop = HelmholtzParameterBlock(pure_component="H2O")
prop.construct()
return prop.htpx(T=T, p=P, x=x)
[docs]@declare_process_block_class("Iapws95ParameterBlock")
class Iapws95ParameterBlockData(HelmholtzParameterBlockData):
CONFIG = HelmholtzParameterBlockData.CONFIG()
CONFIG.pure_component = "H2O"
CONFIG.get("pure_component")._default = "H2O"
[docs]@declare_process_block_class(
"Iapws95StateBlock",
block_class=_StateBlock,
)
class Iapws95StateBlockData(HelmholtzStateBlockData):
pass