#################################################################################
# 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).
#
# Copyright (c) 2018-2024 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.
#################################################################################
"""
This module contains the base class for constructing flowsheet models in the
IDAES modeling framework.
"""
# TODO: Missing docstrings
# pylint: disable=missing-function-docstring
from warnings import warn
import pyomo.environ as pe
from pyomo.dae import ContinuousSet
from pyomo.network import Arc
from pyomo.common.config import ConfigValue, ListOf
from pyomo.core.base.units_container import _PyomoUnit
from idaes.core import (
ProcessBlockData,
declare_process_block_class,
UnitModelBlockData,
useDefault,
)
from idaes.core.util.config import (
is_physical_parameter_block,
is_time_domain,
DefaultBool,
)
from idaes.core.util.misc import add_object_reference
from idaes.core.util.exceptions import DynamicError, ConfigurationError
from idaes.core.util.tables import create_stream_table_dataframe
import idaes.logger as idaeslog
# Some more information about this module
__author__ = "John Eslick, Qi Chen, Andrew Lee"
__all__ = ["FlowsheetBlock", "FlowsheetBlockData"]
# Set up logger
_log = idaeslog.getLogger(__name__)
class UI:
"""Encapsulate checks for installed 'idaes_ui' dependency.
Functions exported by this class will run normally if 'idaes_ui' is installed and
just print warnings if it is not.
Functions:
- `visualize(model, model_name, **kwargs)`
Also has an attribute 'installed' to check directly.
"""
def __init__(self):
# pylint: disable=import-outside-toplevel
try:
import idaes_ui
except ImportError:
idaes_ui = None
if idaes_ui is None:
self.visualize = self._visualize_null
self.installed = False
else:
# FIXME the explicit submodule import is needed because the idaes_ui doesn't import its fv submodule
# otherwise, you get "AttributeError: module 'idaes_ui' has no 'fv' attribute"
import idaes_ui.fv
self.visualize = idaes_ui.fv.visualize
self.installed = True
def _visualize_null(self, model, model_name, **kwargs):
self._warn("idaes_ui.fv.visualize")
@staticmethod
def _warn(name):
message = f"Call to {name}() ignored: 'idaes_ui' package is not installed"
# with stacklevel=3, show the caller of this function's caller
warn(message, category=RuntimeWarning, stacklevel=3)
[docs]
@declare_process_block_class(
"FlowsheetBlock",
doc="""
FlowsheetBlock is a specialized Pyomo block for IDAES flowsheet models, and
contains instances of FlowsheetBlockData.""",
)
class FlowsheetBlockData(ProcessBlockData):
"""
The FlowsheetBlockData Class forms the base class for all IDAES process
flowsheet models. The main purpose of this class is to automate the tasks
common to all flowsheet models and ensure that the necessary attributes of
a flowsheet model are present.
The most signfiicant role of the FlowsheetBlockData class is to
automatically create the time domain for the flowsheet.
"""
# Create Class ConfigBlock
CONFIG = ProcessBlockData.CONFIG()
CONFIG.declare(
"dynamic",
ConfigValue(
default=useDefault,
domain=DefaultBool,
description="Dynamic model flag",
doc="""Indicates whether this model will be dynamic,
**default** - useDefault.
**Valid values:** {
**useDefault** - get flag from parent or False,
**True** - set as a dynamic model,
**False** - set as a steady-state model.}""",
),
)
CONFIG.declare(
"time",
ConfigValue(
default=None,
domain=is_time_domain,
description="Flowsheet time domain",
doc="""Pointer to the time domain for the flowsheet. Users may provide
an existing time domain from another flowsheet, otherwise the flowsheet will
search for a parent with a time domain or create a new time domain and
reference it here.""",
),
)
CONFIG.declare(
"time_set",
ConfigValue(
default=[0],
domain=ListOf(float),
description="Set of points for initializing time domain",
doc="""Set of points for initializing time domain. This should be a
list of floating point numbers,
**default** - [0].""",
),
)
CONFIG.declare(
"time_units",
ConfigValue(
description="Units for time domain",
doc="""Pyomo Units object describing the units of the time domain.
This must be defined for dynamic simulations, default = None.""",
),
)
CONFIG.declare(
"default_property_package",
ConfigValue(
default=None,
domain=is_physical_parameter_block,
description="Default property package to use in flowsheet",
doc="""Indicates the default property package to be used by models
within this flowsheet if not otherwise specified,
**default** - None.
**Valid values:** {
**None** - no default property package,
**a ParameterBlock object**.}""",
),
)
[docs]
def build(self):
"""
General build method for FlowsheetBlockData. This method calls a number
of sub-methods which automate the construction of expected attributes
of flowsheets.
Inheriting models should call `super().build`.
Args:
None
Returns:
None
"""
super(FlowsheetBlockData, self).build()
self._time_units = None
# Set up dynamic flag and time domain
self._setup_dynamics()
@property
def time(self):
# _time will be created by the _setup_dynamics method
return self._time
@property
def time_units(self):
return self._time_units
[docs]
def is_flowsheet(self):
"""
Method which returns True to indicate that this component is a
flowsheet.
Args:
None
Returns:
True
"""
return True
[docs]
def model_check(self):
"""
This method runs model checks on all unit models in a flowsheet.
This method searches for objects which inherit from UnitModelBlockData
and executes the model_check method if it exists.
Args:
None
Returns:
None
"""
_log.info("Executing model checks.")
for o in self.component_objects(descend_into=False):
if isinstance(o, UnitModelBlockData):
try:
o.model_check()
except AttributeError:
# This should never happen, but just in case
_log.warning(
"{} Model/block has no model_check method.".format(o.name)
)
[docs]
def stream_table(self, true_state=False, time_point=0, orient="columns"):
"""
Method to generate a stream table by iterating over all Arcs in the
flowsheet.
Args:
true_state : whether the state variables (True) or display
variables (False, default) from the StateBlocks should
be used in the stream table.
time_point : point in the time domain at which to create stream
table (default = 0)
orient : whether stream should be shown by columns ("columns") or
rows ("index")
Returns:
A pandas dataframe containing stream table information
"""
dict_arcs = {}
for a in self.component_objects(ctype=Arc, descend_into=False):
dict_arcs[a.local_name] = a
return create_stream_table_dataframe(
dict_arcs, time_point=time_point, orient=orient, true_state=true_state
)
[docs]
def visualize(self, model_name, **kwargs) -> "VisualizeResult":
"""
Starts up a flask server that serializes the model and pops up a
webpage with the visualization
Args:
model_name : The name of the model
Keyword Args:
**kwargs: Additional keywords for :func:`idaes.core.ui.fv.visualize()`
Returns:
The :class:`idaes_ui.fv.fsvis.VisualizeResult` instance returned by :meth:`UI.visualize`
"""
visualize_result = UI().visualize(self, model_name, **kwargs)
return visualize_result
def _get_stream_table_contents(self, time_point=0):
"""
Calls stream_table method and returns result
"""
return self.stream_table(time_point)
def _setup_dynamics(self):
# Look for parent flowsheet
fs = self.flowsheet()
# Check the dynamic flag, and retrieve if necessary
if self.config.dynamic == useDefault:
if fs is None:
# No parent, so default to steady-state and warn user
_log.warning(
"{} is a top level flowsheet, but dynamic flag "
"set to useDefault. Dynamic "
"flag set to False by default".format(self.name)
)
self.config.dynamic = False
else:
# Get dynamic flag from parent flowsheet
self.config.dynamic = fs.config.dynamic
# Check for case when dynamic=True, but parent dynamic=False
elif self.config.dynamic is True:
if fs is not None and fs.config.dynamic is False:
raise DynamicError(
"{} trying to declare a dynamic model within "
"a steady-state flowsheet. This is not "
"supported by the IDAES framework. Try "
"creating a dynamic flowsheet instead, and "
"declaring some models as steady-state.".format(self.name)
)
# Validate units for time domain
if self.config.time is None and fs is not None:
# We will get units from parent
pass
elif self.config.time_units is None and self.config.dynamic:
raise ConfigurationError(
f"{self.name} - no units were specified for the time domain. "
f"Units must be specified for dynamic models."
)
elif self.config.time_units is None and not self.config.dynamic:
_log.debug("No units specified for steady-state time domain.")
elif not isinstance(self.config.time_units, _PyomoUnit):
raise ConfigurationError(
"{} unrecognised value for time_units argument. This must be "
"a Pyomo Unit object (not a compound unit).".format(self.name)
)
if self.config.time is not None:
# Validate user provided time domain
if self.config.dynamic is True and not isinstance(
self.config.time, ContinuousSet
):
raise DynamicError(
"{} was set as a dynamic flowsheet, but time domain "
"provided was not a ContinuousSet.".format(self.name)
)
add_object_reference(self, "_time", self.config.time)
self._time_units = self.config.time_units
else:
# If no parent flowsheet, set up time domain
if fs is None:
# Create time domain
if self.config.dynamic:
# Check if time_set has at least two points
if len(self.config.time_set) < 2:
# Check if time_set is at default value
if self.config.time_set == [0.0]:
# If default, set default end point to be 1.0
self.config.time_set = [0.0, 1.0]
else:
# Invalid user input, raise Exception
raise DynamicError(
"Flowsheet provided with invalid "
"time_set attribute - must have at "
"least two values (start and end)."
)
# For dynamics, need a ContinuousSet
self._time = ContinuousSet(initialize=self.config.time_set)
else:
# For steady-state, use an ordered Set
self._time = pe.Set(initialize=self.config.time_set, ordered=True)
self._time_units = self.config.time_units
# Set time config argument as reference to time domain
self.config.time = self._time
else:
# Set time config argument to parent time
self.config.time = fs.time
add_object_reference(self, "_time", fs.time)
# We control time units
# pylint: disable-next=protected-access
self._time_units = fs._time_units