exerpy 0.0.1__py3-none-any.whl → 0.0.2__py3-none-any.whl

exerpy/components/heat_exchanger/steam_generator.py

@@ -7,117 +7,129 @@ from exerpy.components.component import component_registry
 @component_registry
 class SteamGenerator(Component):
     r"""
-    Class for exergy analysis of steam generators.
+    Class for exergy and exergoeconomic analysis of heat exchangers.
 
-    This class performs exergy analysis calculations for steam generators.
-    The component has several input and output streams as follows:
-
-        Inlet streams:
-            - inl[0]: Feed water inlet (high pressure)
-            - inl[1]: Steam inlet (intermediate pressure)
-            - inl[2]: Heat inlet (providing the heat input Q)
-            - inl[3]: Water injection (high pressure)
-            - inl[4]: Water injection (intermediate pressure)
-
-        Outlet streams:
-            - outl[0]: Superheated steam outlet (high pressure)
-            - outl[1]: Superheated steam outlet (intermediate pressure)
-            - outl[2]: Drain / Blow down outlet
-
-    The exergy product is defined as:
-
-    .. math::
-        \dot{E}_P = \Bigl[ \dot{m}_{\mathrm{out,HP}}\,e_{\mathrm{out,HP}}
-        - \dot{m}_{\mathrm{in,HP}}\,e_{\mathrm{in,HP}} \Bigr]
-        + \Bigl[ \dot{m}_{\mathrm{out,IP}}\,e_{\mathrm{out,IP}}
-        - \dot{m}_{\mathrm{in,IP}}\,e_{\mathrm{in,IP}} \Bigr]
-        - \dot{m}_{\mathrm{w,HP}}\,e_{\mathrm{w,HP}}
-        - \dot{m}_{\mathrm{w,IP}}\,e_{\mathrm{w,IP}}
-
-    where:
-        - \(\dot{m}\) is the mass flow rate and \(e\) is the specific exergy of the stream.
-        - The subscripts HP and IP denote high and intermediate pressure streams,
-          respectively, and 'w' stands for water injection.
-
-    The exergy fuel is computed from the heat input as:
-
-    .. math::
-        \dot{E}_F = E_q = Q \left( 1 - \frac{T_b}{T_0} \right)
-
-    with the thermodynamic temperature difference defined by
-
-    .. math::
-        T_b = \frac{h_{\mathrm{out,HP}} - h_{\mathrm{in,HP}}}{s_{\mathrm{out,HP}} - s_{\mathrm{in,HP}}}
-
-    where:
-        - \(h\) and \(s\) are the specific enthalpy and entropy,
-        - \(T_0\) is the ambient temperature.
-
-    The exergy destruction and efficiency are then given by:
-
-    .. math::
-        \dot{E}_D = \dot{E}_F - \dot{E}_P \quad\mathrm{and}\quad \varepsilon = \frac{\dot{E}_P}{\dot{E}_F}
-
-    Parameters
-    ----------
-    **kwargs : dict
-        Arbitrary keyword arguments passed to the parent class.
+    This class performs exergy and exergoeconomic analysis calculations for heat exchanger components,
+    accounting for two inlet and two outlet streams across various temperature regimes, including
+    above and below ambient temperature, and optional dissipative behavior.
 
     Attributes
     ----------
     E_F : float
-        Exergy fuel of the component :math:`\dot{E}_F` in :math:`\mathrm{W}`.
+        Exergy fuel of the component :math:`\dot{E}_\mathrm{F}` in :math:`\mathrm{W}`.
     E_P : float
-        Exergy product of the component :math:`\dot{E}_P` in :math:`\mathrm{W}`.
+        Exergy product of the component :math:`\dot{E}_\mathrm{P}` in :math:`\mathrm{W}`.
     E_D : float
-        Exergy destruction of the component :math:`\dot{E}_D` in :math:`\mathrm{W}`.
+        Exergy destruction of the component :math:`\dot{E}_\mathrm{D}` in :math:`\mathrm{W}`.
     epsilon : float
-        Exergetic efficiency :math:`\varepsilon`.
+        Exergetic efficiency of the component :math:`\varepsilon` in :math:`-`.
     inl : dict
-        Dictionary containing inlet streams.
+        Dictionary containing inlet stream data with mass flows and specific exergies.
     outl : dict
-        Dictionary containing outlet streams.
+        Dictionary containing outlet stream data with mass flows and specific exergies.
+    Z_costs : float
+        Investment cost rate of the component in currency/h.
+    C_P : float
+        Cost of product stream :math:`\dot{C}_P` in currency/h.
+    C_F : float
+        Cost of fuel stream :math:`\dot{C}_F` in currency/h.
+    C_D : float
+        Cost of exergy destruction :math:`\dot{C}_D` in currency/h.
+    c_P : float
+        Specific cost of product stream (currency per unit exergy).
+    c_F : float
+        Specific cost of fuel stream (currency per unit exergy).
+    r : float
+        Relative cost difference, :math:`(c_P - c_F)/c_F`.
+    f : float
+        Exergoeconomic factor, :math:`\dot{Z}/(\dot{Z} + \dot{C}_D)`.
+    Ex_C_col : dict
+        Custom cost coefficients collection passed via `kwargs`.
+
+    Notes
+    -----
+    The component has several input and output streams as follows.
+
+    Inlet streams:
+
+    - inl[0]: Feed water inlet (high pressure)
+    - inl[1]: Steam inlet (intermediate pressure)
+    - inl[2]: Heat inlet (providing the heat input Q)
+    - inl[3]: Water injection (high pressure)
+    - inl[4]: Water injection (intermediate pressure)
+
+    Outlet streams:
+
+    - outl[0]: Superheated steam outlet (high pressure)
+    - outl[1]: Superheated steam outlet (intermediate pressure)
+    - outl[2]: Drain / Blow down outlet
+
     """
 
     def __init__(self, **kwargs):
-        r"""Initialize steam generator component with given parameters."""
-        self.dissipative = False
+        r"""
+        Initialize the steam generator component.
+
+        Parameters
+        ----------
+        **kwargs : dict
+            Arbitrary keyword arguments. Recognized keys:
+            - Ex_C_col (dict): custom cost coefficients, default {}
+            - Z_costs (float): investment cost rate in currency/h, default 0.0
+        """        
         super().__init__(**kwargs)
 
     def calc_exergy_balance(self, T0: float, p0: float, split_physical_exergy) -> None:
         r"""
-        Calculate the exergy balance of the steam generator.
+        Compute the exergy balance of the steam generator.
 
-        This method computes the exergy fuel from the heat inlet using the relation
+        The exergy fuel is defined as follows. 
+
+        If `split_physical_exergy` is `True`:
 
         .. math::
-            E_F = Q \left(1 - \frac{T_b}{T_0}\right),
-            \quad T_b = \frac{h_{\mathrm{out,HP}} - h_{\mathrm{in,HP}}}{s_{\mathrm{out,HP}} - s_{\mathrm{in,HP}}}
+            \dot{E}_{\mathrm{F}}
+            = \bigl[\dot{E}^{\mathrm{T}}_{\mathrm{out,HP}} - \dot{E}^{\mathrm{T}}_{\mathrm{in,HP}}\bigr]
+            + \bigl[\dot{E}^{\mathrm{T}}_{\mathrm{out,IP}} - \dot{E}^{\mathrm{T}}_{\mathrm{in,IP}}\bigr]
+            - \dot{E}^{\mathrm{T}}_{\mathrm{w,HP}}
+            - \dot{E}^{\mathrm{T}}_{\mathrm{w,IP}}
 
-        and the exergy product as
+        If `split_physical_exergy` is `False`:
 
         .. math::
-            E_P = \Bigl[ \dot{m}_{\mathrm{out,HP}}\,e_{\mathrm{out,HP}} -
-            \dot{m}_{\mathrm{in,HP}}\,e_{\mathrm{in,HP}} \Bigr]
-            + \Bigl[ \dot{m}_{\mathrm{out,IP}}\,e_{\mathrm{out,IP}} -
-            \dot{m}_{\mathrm{in,IP}}\,e_{\mathrm{in,IP}} \Bigr]
-            - \dot{m}_{\mathrm{w,HP}}\,e_{\mathrm{w,HP}}
-            - \dot{m}_{\mathrm{w,IP}}\,e_{\mathrm{w,IP}}
+            \dot{E}_{\mathrm{F}}
+            = \bigl[\dot{E}^{\mathrm{PH}}_{\mathrm{out,HP}} - \dot{E}^{\mathrm{PH}}_{\mathrm{in,HP}}\bigr]
+            + \bigl[\dot{E}^{\mathrm{PH}}_{\mathrm{out,IP}} - \dot{E}^{\mathrm{PH}}_{\mathrm{in,IP}}\bigr]
+            - \dot{E}^{\mathrm{PH}}_{\mathrm{w,HP}}
+            - \dot{E}^{\mathrm{PH}}_{\mathrm{w,IP}}
 
-        The exergy destruction is given by
+        The exergy product is defined as:
 
         .. math::
-            E_D = E_F - E_P
 
-        and the exergetic efficiency is
+            \dot{E}_\mathrm{P} = \Bigl[ \dot E^{\mathrm{PH}}_{\mathrm{out,HP}}
+            - \dot E^{\mathrm{PH}}_{\mathrm{in,HP}} \Bigr]
+            + \Bigl[ \dot E^{\mathrm{PH}}_{\mathrm{out,IP}}
+            - \dot E^{\mathrm{PH}}_{\mathrm{in,IP}} \Bigr]
+            - \dot E^{\mathrm{PH}}_{\mathrm{w,HP}}
+            - \dot E^{\mathrm{PH}}_{\mathrm{w,IP}}
 
-        .. math::
-            \varepsilon = \frac{E_P}{E_F}
+        where the subscripts HP and IP denote high and intermediate pressure streams,
+        respectively, and w stands for water injection.
 
         Parameters
         ----------
         T0 : float
-            Ambient temperature in Kelvin.
+            Ambient temperature (K).
+        p0 : float
+            Ambient pressure (Pa).
+        split_physical_exergy : bool
+            Whether to split thermal and mechanical exergy.
+
+        Raises
+        ------
+        ValueError
+            If required inlets or outlets are missing.
         """
         # Ensure that all necessary streams exist
         required_inlets = [0]
@@ -129,12 +141,6 @@ class SteamGenerator(Component):
             if idx not in self.outl:
                 raise ValueError(f"Missing outlet stream with index {idx}.")
 
-        # Calculate T_b for high pressure streams: from superheated steam outlet (HP) and feed water inlet (HP)
-        try:
-            T_b = (self.outl[0]['h'] - self.inl[0]['h']) / (self.outl[0]['s'] - self.inl[0]['s'])
-        except ZeroDivisionError:
-            raise ZeroDivisionError("Division by zero encountered in calculating T_b. Check entropy differences.")
-       
         if split_physical_exergy:
             exergy_type = 'e_T'
         else:
@@ -172,6 +178,15 @@ class SteamGenerator(Component):
 
 
     def aux_eqs(self, A, b, counter, T0, equations, chemical_exergy_enabled):
+        r"""
+        This function must be implemented in the future.
+
+        The exergoeconomic analysis of SteamGenerator is not implemented yet.
+        """
+        logging.error(
+            "The exergoeconomic analysis of SteamGenerator is not implemented yet. "
+            "This method will be implemented in a future release."
+        )
         """
         Auxiliary equations for the steam generator.
         
@@ -218,6 +233,16 @@ class SteamGenerator(Component):
         """
 
     def exergoeconomic_balance(self, T0):
+        r"""
+        This function must be implemented in the future.
+
+        The exergoeconomic analysis of SteamGenerator is not implemented yet.
+        """
+        
+        logging.error(
+            "The exergoeconomic analysis of SteamGenerator is not implemented yet. "
+            "This method will be implemented in a future release."
+        )
         """
         Perform exergoeconomic balance calculations for the steam generator.
         

exerpy/components/nodes/deaerator.py

@@ -44,53 +44,52 @@ class Deaerator(Component):
     product and fuel are defined based on temperature relationships:
 
     .. math::
-
-        \dot{E}_\mathrm{P} =
-        \begin{cases}
-        \begin{cases}
-        \sum_i \dot{m}_i \cdot (e_\mathrm{out}^\mathrm{PH} -
-        e_{\mathrm{in,}i}^\mathrm{PH})
-        & T_{\mathrm{in,}i} < T_\mathrm{out} \mathrm{ & }
-        T_{\mathrm{in,}i} \geq T_0 \\
-        \sum_i \dot{m}_i \cdot e_\mathrm{out}^\mathrm{PH}
-        & T_{\mathrm{in,}i} < T_\mathrm{out} \mathrm{ & }
-        T_{\mathrm{in,}i} < T_0 \\
-        \end{cases} & T_\mathrm{out} > T_0\\
-        \mathrm{not defined (nan)} & T_\mathrm{out} = T_0\\
+        \displaystyle
+        \dot E_{P} =
         \begin{cases}
-        \sum_i \dot{m}_i \cdot e_\mathrm{out}^\mathrm{PH}
-        & T_{\mathrm{in,}i} > T_\mathrm{out} \mathrm{ & }
-        T_{\mathrm{in,}i} \geq T_0 \\
-        \sum_i \dot{m}_i \cdot (e_\mathrm{out}^\mathrm{PH} -
-        e_{\mathrm{in,}i}^\mathrm{PH})
-        & T_{\mathrm{in,}i} > T_\mathrm{out} \mathrm{ & }
-        T_{\mathrm{in,}i} < T_0 \\
-        \end{cases} & T_\mathrm{out} < T_0\\
+            \displaystyle
+            \sum_{i}\dot m_{i}\,\bigl(e_{\mathrm{out}}^{\mathrm{PH}}
+            -e_{\mathrm{in},i}^{\mathrm{PH}}\bigr),
+            \quad \text{if }T_{\mathrm{in},i}<T_{\mathrm{out}}\text{ and }T_{\mathrm{in},i}\ge T_{0},\\[8pt]
+            \displaystyle
+            \sum_{i}\dot m_{i}\,e_{\mathrm{out}}^{\mathrm{PH}},
+            \quad \text{if }T_{\mathrm{in},i}<T_{\mathrm{out}}\text{ and }T_{\mathrm{in},i}<T_{0},\\[8pt]
+            \displaystyle
+            \text{not defined (nan)},
+            \quad \text{if }T_{\mathrm{out}}=T_{0},\\[8pt]
+            \displaystyle
+            \sum_{i}\dot m_{i}\,e_{\mathrm{out}}^{\mathrm{PH}},
+            \quad \text{if }T_{\mathrm{in},i}>T_{\mathrm{out}}\text{ and }T_{\mathrm{in},i}\ge T_{0},\\[8pt]
+            \displaystyle
+            \sum_{i}\dot m_{i}\,\bigl(e_{\mathrm{out}}^{\mathrm{PH}}
+            -e_{\mathrm{in},i}^{\mathrm{PH}}\bigr),
+            \quad \text{if }T_{\mathrm{in},i}>T_{\mathrm{out}}\text{ and }T_{\mathrm{in},i}<T_{0}.
         \end{cases}
 
-        \dot{E}_\mathrm{F} =
-        \begin{cases}
-        \begin{cases}
-        \sum_i \dot{m}_i \cdot (e_{\mathrm{in,}i}^\mathrm{PH} -
-        e_\mathrm{out}^\mathrm{PH})
-        & T_{\mathrm{in,}i} > T_\mathrm{out} \\
-        \sum_i \dot{m}_i \cdot e_{\mathrm{in,}i}^\mathrm{PH}
-        & T_{\mathrm{in,}i} < T_\mathrm{out} \mathrm{ & }
-        T_{\mathrm{in,}i} < T_0 \\
-        \end{cases} & T_\mathrm{out} > T_0\\
-        \sum_i \dot{m}_i \cdot e_{\mathrm{in,}i}^\mathrm{PH}
-        & T_\mathrm{out} = T_0\\
+    .. math::
+        \displaystyle
+        \dot E_{F} =
         \begin{cases}
-        \sum_i \dot{m}_i \cdot e_{\mathrm{in,}i}^\mathrm{PH}
-        & T_{\mathrm{in,}i} > T_\mathrm{out} \mathrm{ & }
-        T_{\mathrm{in,}i} \geq T_0 \\
-        \sum_i \dot{m}_i \cdot (e_{\mathrm{in,}i}^\mathrm{PH} -
-        e_\mathrm{out}^\mathrm{PH})
-        & T_{\mathrm{in,}i} < T_\mathrm{out} \\
-        \end{cases} & T_\mathrm{out} < T_0\\
+            \displaystyle
+            \sum_{i}\dot m_{i}\,\bigl(e_{\mathrm{in},i}^{\mathrm{PH}}
+            -e_{\mathrm{out}}^{\mathrm{PH}}\bigr),
+            \quad \text{if }T_{\mathrm{out}}>T_{0}\text{ and }T_{\mathrm{in},i}>T_{\mathrm{out}},\\[8pt]
+            \displaystyle
+            \sum_{i}\dot m_{i}\,e_{\mathrm{in},i}^{\mathrm{PH}},
+            \quad \text{if }T_{\mathrm{out}}>T_{0}\text{ and }T_{\mathrm{in},i}<T_{\mathrm{out}}
+            \text{ and }T_{\mathrm{in},i}<T_{0},\\[8pt]
+            \displaystyle
+            \sum_{i}\dot m_{i}\,e_{\mathrm{in},i}^{\mathrm{PH}},
+            \quad \text{if }T_{\mathrm{out}}=T_{0},\\[8pt]
+            \displaystyle
+            \sum_{i}\dot m_{i}\,e_{\mathrm{in},i}^{\mathrm{PH}},
+            \quad \text{if }T_{\mathrm{out}}<T_{0}\text{ and }T_{\mathrm{in},i}>T_{\mathrm{out}},\\[8pt]
+            \displaystyle
+            \sum_{i}\dot m_{i}\,\bigl(e_{\mathrm{in},i}^{\mathrm{PH}}
+            -e_{\mathrm{out}}^{\mathrm{PH}}\bigr),
+            \quad \text{if }T_{\mathrm{out}}<T_{0}\text{ and }T_{\mathrm{in},i}<T_{\mathrm{out}}.
         \end{cases}
 
-        \forall i \in \mathrm{deaerator inlets}
     """
 
     def __init__(self, **kwargs):
@@ -177,14 +176,16 @@ class Deaerator(Component):
         the following auxiliary cost relations:
         
         (1) Mixing equation for chemical exergy costs (if enabled):
-            - The outlet's specific chemical exergy cost is calculated as a mass-weighted 
-              average of the inlet streams' specific chemical exergy costs
-            - This enforces proper chemical exergy cost distribution through the deaerator
+
+        - The outlet's specific chemical exergy cost is calculated as a mass-weighted average of the inlet streams' specific chemical exergy costs
+        
+        - This enforces proper chemical exergy cost distribution through the deaerator
         
         (2) Mixing equation for mechanical exergy costs:
-            - The outlet's specific mechanical exergy cost is calculated as a mass-weighted
-              average of the inlet streams' specific mechanical exergy costs
-            - This ensures mechanical exergy costs are properly conserved in the mixing process
+        
+        - The outlet's specific mechanical exergy cost is calculated as a mass-weighted average of the inlet streams' specific mechanical exergy costs
+        
+        - This ensures mechanical exergy costs are properly conserved in the mixing process
             
         Both equations implement the proportionality rule for mixing processes where
         the outlet's specific costs should reflect the contribution of each inlet stream.

exerpy/components/nodes/mixer.py

@@ -44,53 +44,52 @@ class Mixer(Component):
     product and fuel are defined based on temperature relationships:
 
     .. math::
-
-        \dot{E}_\mathrm{P} =
-        \begin{cases}
-        \begin{cases}
-        \sum_i \dot{m}_i \cdot (e_\mathrm{out}^\mathrm{PH} -
-        e_{\mathrm{in,}i}^\mathrm{PH})
-        & T_{\mathrm{in,}i} < T_\mathrm{out} \mathrm{ & }
-        T_{\mathrm{in,}i} \geq T_0 \\
-        \sum_i \dot{m}_i \cdot e_\mathrm{out}^\mathrm{PH}
-        & T_{\mathrm{in,}i} < T_\mathrm{out} \mathrm{ & }
-        T_{\mathrm{in,}i} < T_0 \\
-        \end{cases} & T_\mathrm{out} > T_0\\
-        \mathrm{not defined (nan)} & T_\mathrm{out} = T_0\\
+        \displaystyle
+        \dot E_{P} =
         \begin{cases}
-        \sum_i \dot{m}_i \cdot e_\mathrm{out}^\mathrm{PH}
-        & T_{\mathrm{in,}i} > T_\mathrm{out} \mathrm{ & }
-        T_{\mathrm{in,}i} \geq T_0 \\
-        \sum_i \dot{m}_i \cdot (e_\mathrm{out}^\mathrm{PH} -
-        e_{\mathrm{in,}i}^\mathrm{PH})
-        & T_{\mathrm{in,}i} > T_\mathrm{out} \mathrm{ & }
-        T_{\mathrm{in,}i} < T_0 \\
-        \end{cases} & T_\mathrm{out} < T_0\\
+            \displaystyle
+            \sum_{i}\dot m_{i}\,\bigl(e_{\mathrm{out}}^{\mathrm{PH}}
+            -e_{\mathrm{in},i}^{\mathrm{PH}}\bigr),
+            \quad \text{if }T_{\mathrm{in},i}<T_{\mathrm{out}}\text{ and }T_{\mathrm{in},i}\ge T_{0},\\[8pt]
+            \displaystyle
+            \sum_{i}\dot m_{i}\,e_{\mathrm{out}}^{\mathrm{PH}},
+            \quad \text{if }T_{\mathrm{in},i}<T_{\mathrm{out}}\text{ and }T_{\mathrm{in},i}<T_{0},\\[8pt]
+            \displaystyle
+            \text{not defined (nan)},
+            \quad \text{if }T_{\mathrm{out}}=T_{0},\\[8pt]
+            \displaystyle
+            \sum_{i}\dot m_{i}\,e_{\mathrm{out}}^{\mathrm{PH}},
+            \quad \text{if }T_{\mathrm{in},i}>T_{\mathrm{out}}\text{ and }T_{\mathrm{in},i}\ge T_{0},\\[8pt]
+            \displaystyle
+            \sum_{i}\dot m_{i}\,\bigl(e_{\mathrm{out}}^{\mathrm{PH}}
+            -e_{\mathrm{in},i}^{\mathrm{PH}}\bigr),
+            \quad \text{if }T_{\mathrm{in},i}>T_{\mathrm{out}}\text{ and }T_{\mathrm{in},i}<T_{0}.
         \end{cases}
 
-        \dot{E}_\mathrm{F} =
-        \begin{cases}
-        \begin{cases}
-        \sum_i \dot{m}_i \cdot (e_{\mathrm{in,}i}^\mathrm{PH} -
-        e_\mathrm{out}^\mathrm{PH})
-        & T_{\mathrm{in,}i} > T_\mathrm{out} \\
-        \sum_i \dot{m}_i \cdot e_{\mathrm{in,}i}^\mathrm{PH}
-        & T_{\mathrm{in,}i} < T_\mathrm{out} \mathrm{ & }
-        T_{\mathrm{in,}i} < T_0 \\
-        \end{cases} & T_\mathrm{out} > T_0\\
-        \sum_i \dot{m}_i \cdot e_{\mathrm{in,}i}^\mathrm{PH}
-        & T_\mathrm{out} = T_0\\
+    .. math::
+        \displaystyle
+        \dot E_{F} =
         \begin{cases}
-        \sum_i \dot{m}_i \cdot e_{\mathrm{in,}i}^\mathrm{PH}
-        & T_{\mathrm{in,}i} > T_\mathrm{out} \mathrm{ & }
-        T_{\mathrm{in,}i} \geq T_0 \\
-        \sum_i \dot{m}_i \cdot (e_{\mathrm{in,}i}^\mathrm{PH} -
-        e_\mathrm{out}^\mathrm{PH})
-        & T_{\mathrm{in,}i} < T_\mathrm{out} \\
-        \end{cases} & T_\mathrm{out} < T_0\\
+            \displaystyle
+            \sum_{i}\dot m_{i}\,\bigl(e_{\mathrm{in},i}^{\mathrm{PH}}
+            -e_{\mathrm{out}}^{\mathrm{PH}}\bigr),
+            \quad \text{if }T_{\mathrm{out}}>T_{0}\text{ and }T_{\mathrm{in},i}>T_{\mathrm{out}},\\[8pt]
+            \displaystyle
+            \sum_{i}\dot m_{i}\,e_{\mathrm{in},i}^{\mathrm{PH}},
+            \quad \text{if }T_{\mathrm{out}}>T_{0}\text{ and }T_{\mathrm{in},i}<T_{\mathrm{out}}
+            \text{ and }T_{\mathrm{in},i}<T_{0},\\[8pt]
+            \displaystyle
+            \sum_{i}\dot m_{i}\,e_{\mathrm{in},i}^{\mathrm{PH}},
+            \quad \text{if }T_{\mathrm{out}}=T_{0},\\[8pt]
+            \displaystyle
+            \sum_{i}\dot m_{i}\,e_{\mathrm{in},i}^{\mathrm{PH}},
+            \quad \text{if }T_{\mathrm{out}}<T_{0}\text{ and }T_{\mathrm{in},i}>T_{\mathrm{out}},\\[8pt]
+            \displaystyle
+            \sum_{i}\dot m_{i}\,\bigl(e_{\mathrm{in},i}^{\mathrm{PH}}
+            -e_{\mathrm{out}}^{\mathrm{PH}}\bigr),
+            \quad \text{if }T_{\mathrm{out}}<T_{0}\text{ and }T_{\mathrm{in},i}<T_{\mathrm{out}}.
         \end{cases}
 
-        \forall i \in \mathrm{mixer inlets}
     """
 
     def __init__(self, **kwargs):
@@ -193,16 +192,18 @@ class Mixer(Component):
         This function adds rows to the cost matrix A and the right-hand-side vector b to enforce
         the following auxiliary cost relations:
         
-        (1) Chemical exergy cost equation (if enabled):
-            - F-principle: The specific chemical exergy cost of the outlet stream is calculated as 
-              the weighted average of the specific chemical exergy costs of the inlet streams.
-            - For inlets with zero chemical exergy, their specific costs are directly transferred.
+        (1) Mixing equation for chemical exergy costs (if enabled):
+
+        - The outlet's specific chemical exergy cost is calculated as a mass-weighted average of the inlet streams' specific chemical exergy costs
+        
+        - This enforces proper chemical exergy cost distribution through the deaerator
+        
+        (2) Mixing equation for mechanical exergy costs:
         
-        (2) Mechanical exergy cost equation:
-            - F-principle: The specific mechanical exergy cost of the outlet stream is calculated as
-              the weighted average of the specific mechanical exergy costs of the inlet streams.
-            - For inlets with zero mechanical exergy, their specific costs are directly transferred.
+        - The outlet's specific mechanical exergy cost is calculated as a mass-weighted average of the inlet streams' specific mechanical exergy costs
         
+        - This ensures mechanical exergy costs are properly conserved in the mixing process
+       
         Parameters
         ----------
         A : numpy.ndarray

exerpy/components/piping/valve.py

@@ -9,17 +9,11 @@ from exerpy.components.component import component_registry
 @component_registry
 class Valve(Component):
     r"""
-    Class for exergy analysis of valves.
+    Class for exergy and exergoeconomic analysis of valves.
 
-    This class performs exergy analysis calculations for isenthalpic valves with 
-    one inlet and one outlet stream. The exergy product and fuel definitions
-    vary based on the temperature relationships between inlet stream, outlet stream,
-    and ambient conditions.
-
-    Parameters
-    ----------
-    **kwargs : dict
-        Arbitrary keyword arguments passed to parent class.
+    This class performs exergy and exergoeconomic analysis calculations for valve components,
+    accounting for one inlet and one outlet streams across various temperature regimes, including
+    above and below ambient temperature.
 
     Attributes
     ----------
@@ -32,11 +26,27 @@ class Valve(Component):
     epsilon : float
         Exergetic efficiency of the component :math:`\varepsilon` in :math:`-`.
     inl : dict
-        Dictionary containing inlet stream data with temperature, mass flows,
-        and specific exergies.
+        Dictionary containing inlet stream data with mass flows and specific exergies.
     outl : dict
-        Dictionary containing outlet stream data with temperature, mass flows,
-        and specific exergies.
+        Dictionary containing outlet stream data with mass flows and specific exergies.
+    Z_costs : float
+        Investment cost rate of the component in currency/h.
+    C_P : float
+        Cost of product stream :math:`\dot{C}_P` in currency/h.
+    C_F : float
+        Cost of fuel stream :math:`\dot{C}_F` in currency/h.
+    C_D : float
+        Cost of exergy destruction :math:`\dot{C}_D` in currency/h.
+    c_P : float
+        Specific cost of product stream (currency per unit exergy).
+    c_F : float
+        Specific cost of fuel stream (currency per unit exergy).
+    r : float
+        Relative cost difference, :math:`(c_P - c_F)/c_F`.
+    f : float
+        Exergoeconomic factor, :math:`\dot{Z}/(\dot{Z} + \dot{C}_D)`.
+    Ex_C_col : dict
+        Custom cost coefficients collection passed via `kwargs`.
 
     Notes
     -----
@@ -179,14 +189,17 @@ class Valve(Component):
             - Valve is treated as dissipative (warning issued)
         
         For T_out <= T0:
-            (1) 1/E_M_in * C_M_in - 1/E_M_out * C_M_out = 0
-                - F-principle: specific mechanical exergy costs equalized between inlet/outlet
-                - If E_M is zero for either stream, appropriate fallback coefficients are used
+        (1) 1/E_M_in * C_M_in - 1/E_M_out * C_M_out = 0
+
+        - F-principle: specific mechanical exergy costs equalized between inlet/outlet
+
+        - If E_M is zero for either stream, appropriate fallback coefficients are used
         
         When chemical_exergy_enabled is True:
-            (2) 1/E_CH_in * C_CH_in - 1/E_CH_out * C_CH_out = 0
-                - F-principle: specific chemical exergy costs equalized between inlet/outlet
-                - If E_CH is zero for either stream, appropriate fallback coefficients are used
+        (2) 1/E_CH_in * C_CH_in - 1/E_CH_out * C_CH_out = 0
+
+        - F-principle: specific chemical exergy costs equalized between inlet/outlet
+        - If E_CH is zero for either stream, appropriate fallback coefficients are used
         
         Parameters
         ----------
@@ -254,8 +267,12 @@ class Valve(Component):
         components (non-dissipative and non-CycleCloser) in proportion to their exergy destruction (E_D)
         and adding an extra overall cost balance row that enforces:
         
-            (C_in,T - C_out,T) + (C_in,M - C_out,M) - C_diff = - Z_costs
-        
+        .. math::
+           (\dot C_{\mathrm{in},T} - \dot C_{\mathrm{out},T})
+           + (\dot C_{\mathrm{in},M} - \dot C_{\mathrm{out},M})
+           - \dot C_{\mathrm{diff}}
+           = -\,\dot Z_{\mathrm{costs}}
+
         In this formulation, the unknown cost variable in the "dissipative" column (i.e. C_diff)
         is solved for, ensuring the valve’s cost balance.
         
@@ -283,8 +300,7 @@ class Valve(Component):
         
         Notes
         -----
-        - It is assumed that each inlet/outlet stream’s CostVar_index dictionary has keys:
-        "T" (thermal), "M" (mechanical), and "dissipative" (the extra unknown).
+        - It is assumed that each inlet/outlet stream's CostVar_index dictionary has keys: "T" (thermal), "M" (mechanical), and "dissipative" (the extra unknown).
         - self.Z_costs is the known cost rate (in currency/s) for the valve.
         """
         # --- Thermal difference row ---

exerpy/parser/from_ebsilon/ebsilon_config.py

@@ -832,7 +832,7 @@ composition_params = [
 # Define the component groups via unique labels
 grouped_components = {
     "Turbine": [6, 23, 56, 57, 58, 68, 122],
-    "HeatExchanger": [10, 15, 16, 25, 26, 27, 43, 51, 55, 61, 62, 70, 71, 124, 126],
+    "HeatExchanger": [10, 25, 26, 27, 43, 51, 55, 61, 62, 70, 71, 124, 126],
     "CombustionChamber": [22, 90],
     "Valve": [2, 13, 14, 39, 42, 59, 68, 133],
     "Pump": [8, 44, 83, 159],

{exerpy-0.0.1.dist-info → exerpy-0.0.2.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: exerpy
-Version: 0.0.1
+Version: 0.0.2
 Summary: Exergy analysis for tabular input data
 Author-email: Sergio Tomasinelli <s.tomasinelli@tu-berlin.de>, Robert Müller <robert.mueller.2@tu-berlin.de>, Francesco Witte <github@witte.sh>
 Requires-Python: >=3.10
@@ -33,10 +33,10 @@ Requires-Dist: sphinx>=7.2.2 ; extra == "dev"
 Requires-Dist: sphinx-copybutton ; extra == "dev"
 Requires-Dist: sphinx-design ; extra == "dev"
 Requires-Dist: sphinxcontrib.bibtex ; extra == "dev"
-Requires-Dist: tespy ; extra == "dev"
+Requires-Dist: tespy>=0.9 ; extra == "dev"
 Requires-Dist: tox ; extra == "dev"
 Requires-Dist: pywin32 ; extra == "ebsilon"
-Requires-Dist: tespy ; extra == "tespy"
+Requires-Dist: tespy>=0.9 ; extra == "tespy"
 Project-URL: Changelog, https://exerpy.readthedocs.io/en/latest/whats_new.html
 Project-URL: Documentation, https://exerpy.readthedocs.io/
 Project-URL: Homepage, https://github.com/oemof/exerpy
@@ -126,7 +126,7 @@ ongoing development. You can cite ExerPy using the following BibTeX entry:
          title = {{ExerPy}: Exergy Analysis in Python},
          note = {Supervision: Prof. Dr.-Ing. Fontina Petrakopoulou}
          url = {https://github.com/oemof/exerpy},
-         version = {0.0.1},
+         version = {0.0.2},
          year = {2025}
     }