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

exerpy/components/heat_exchanger/base.py

@@ -2,24 +2,17 @@ import logging
 
 import numpy as np
 
-from exerpy.components.component import Component
-from exerpy.components.component import component_registry
+from exerpy.components.component import Component, component_registry
 
 
 @component_registry
 class HeatExchanger(Component):
     r"""
-    Class for exergy analysis of heat exchangers.
+    Class for exergy and exergoeconomic analysis of heat exchangers.
 
-    This class performs exergy analysis calculations for heat exchangers, considering
-    different temperature regimes relative to the ambient temperature. The exergy
-    product and fuel definitions vary based on the temperature levels of the streams.
-    Stream 0 represents the hot stream and stream 1 represents the cold stream.
-
-    Parameters
-    ----------
-    **kwargs : dict
-        Arbitrary keyword arguments passed to 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
     ----------
@@ -32,133 +25,228 @@ class HeatExchanger(Component):
     epsilon : float
         Exergetic efficiency of the component :math:`\varepsilon` in :math:`-`.
     inl : dict
-        Dictionary containing inlet streams data (0: hot stream, 1: cold stream) with
-        temperature, mass flows, and specific exergies.
+        Dictionary containing inlet stream data with mass flows and specific exergies.
     outl : dict
-        Dictionary containing outlet streams data (0: hot stream, 1: cold stream) with
-        temperature, mass flows, and specific exergies.
-
-    Notes
-    -----
-    The exergy analysis considers six different cases based on stream temperatures
-    relative to ambient temperature :math:`T_0`. Stream indices refer to hot (0) and
-    cold (1) streams:
-
-    Case 1 - **All streams above ambient temperature**:
-
-    .. math::
-        \dot{E}_\mathrm{P} &= \dot{m}_{\mathrm{out,1}} \cdot e^\mathrm{T}_{\mathrm{out,1}}
-        - \dot{m}_{\mathrm{in,1}} \cdot e^\mathrm{T}_{\mathrm{in,1}}\\
-        \dot{E}_\mathrm{F} &= \dot{m}_{\mathrm{in,0}} \cdot e^\mathrm{PH}_{\mathrm{in,0}}
-        - \dot{m}_{\mathrm{out,0}} \cdot e^\mathrm{PH}_{\mathrm{out,0}} +
-        \dot{m}_{\mathrm{in,1}} \cdot e^\mathrm{M}_{\mathrm{in,1}}
-        - \dot{m}_{\mathrm{out,1}} \cdot e^\mathrm{M}_{\mathrm{out,1}}
-
-    Case 2 - **All streams below or at ambient temperature**:
-
-    .. math::
-        \dot{E}_\mathrm{P} &= \dot{m}_{\mathrm{out,0}} \cdot e^\mathrm{T}_{\mathrm{out,0}}
-        - \dot{m}_{\mathrm{in,0}} \cdot e^\mathrm{T}_{\mathrm{in,0}}\\
-        \dot{E}_\mathrm{F} &= \dot{m}_{\mathrm{in,1}} \cdot e^\mathrm{PH}_{\mathrm{in,1}}
-        - \dot{m}_{\mathrm{out,1}} \cdot e^\mathrm{PH}_{\mathrm{out,1}} +
-        \dot{m}_{\mathrm{in,0}} \cdot e^\mathrm{M}_{\mathrm{in,0}}
-        - \dot{m}_{\mathrm{out,0}} \cdot e^\mathrm{M}_{\mathrm{out,0}}
-
-    Case 3 - **Hot stream inlet/outlet above ambient, cold stream inlet/outlet below ambient**:
-
-    .. math::
-        \dot{E}_\mathrm{P} &= \dot{m}_{\mathrm{out,0}} \cdot e^\mathrm{T}_{\mathrm{out,0}}
-        + \dot{m}_{\mathrm{out,1}} \cdot e^\mathrm{T}_{\mathrm{out,1}}\\
-        \dot{E}_\mathrm{F} &= \dot{m}_{\mathrm{in,0}} \cdot e^\mathrm{PH}_{\mathrm{in,0}}
-        + \dot{m}_{\mathrm{in,1}} \cdot e^\mathrm{PH}_{\mathrm{in,1}}
-        - (\dot{m}_{\mathrm{out,0}} \cdot e^\mathrm{M}_{\mathrm{out,0}}
-        + \dot{m}_{\mathrm{out,1}} \cdot e^\mathrm{M}_{\mathrm{out,1}})
-
-    Case 4 - **First inlet above ambient, all other streams below or at ambient**:
-
-    .. math::
-        \dot{E}_\mathrm{P} &= \dot{m}_{\mathrm{out,0}} \cdot e^\mathrm{T}_{\mathrm{out,0}}\\
-        \dot{E}_\mathrm{F} &= \dot{m}_{\mathrm{in,0}} \cdot e^\mathrm{PH}_{\mathrm{in,0}}
-        + \dot{m}_{\mathrm{in,1}} \cdot e^\mathrm{PH}_{\mathrm{in,1}}
-        - (\dot{m}_{\mathrm{out,1}} \cdot e^\mathrm{PH}_{\mathrm{out,1}}
-        + \dot{m}_{\mathrm{out,0}} \cdot e^\mathrm{M}_{\mathrm{out,0}})
-
-    Case 5 - **Hot stream inlet/outlet above ambient, cold stream inlet/outlet below or at ambient**:
-
-    .. math::
-        \dot{E}_\mathrm{P} &= \mathrm{NaN}\\
-        \dot{E}_\mathrm{F} &= \dot{m}_{\mathrm{in,0}} \cdot e^\mathrm{PH}_{\mathrm{in,0}}
-        - \dot{m}_{\mathrm{out,0}} \cdot e^\mathrm{PH}_{\mathrm{out,0}}
-        + \dot{m}_{\mathrm{in,1}} \cdot e^\mathrm{PH}_{\mathrm{in,1}}
-        - \dot{m}_{\mathrm{out,1}} \cdot e^\mathrm{PH}_{\mathrm{out,1}}
-
-    Case 6 - **Second outlet above ambient, all others below or at ambient**:
-
-    .. math::
-        \dot{E}_\mathrm{P} &= \dot{m}_{\mathrm{out,1}} \cdot e^\mathrm{T}_{\mathrm{out,1}}\\
-        \dot{E}_\mathrm{F} &= \dot{m}_{\mathrm{in,0}} \cdot e^\mathrm{PH}_{\mathrm{in,0}}
-        - \dot{m}_{\mathrm{out,0}} \cdot e^\mathrm{PH}_{\mathrm{out,0}}
-        + \dot{m}_{\mathrm{in,1}} \cdot e^\mathrm{PH}_{\mathrm{in,1}}
-        - \dot{m}_{\mathrm{out,1}} \cdot e^\mathrm{M}_{\mathrm{out,1}}
-
-    Note that in Case 5, the exergy product :math:`\dot{E}_\mathrm{P}` is undefined (NaN),
-    leading to an exergy destruction equal to the exergy fuel:
-    :math:`\dot{E}_\mathrm{D} = \dot{E}_\mathrm{F}`.
-
-    The exergy destruction is calculated as:
-
-    .. math::
-        \dot{E}_\mathrm{D} = \dot{E}_\mathrm{F} - \dot{E}_\mathrm{P}
-
-    And the exergetic efficiency as:
-
-    .. math::
-        \varepsilon = \frac{\dot{E}_\mathrm{P}}{\dot{E}_\mathrm{F}}
-
-    Where:
-        - :math:`e^\mathrm{T}`: Thermal exergy
-        - :math:`e^\mathrm{PH}`: Physical exergy
-        - :math:`e^\mathrm{M}`: Mechanical exergy
+        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`.
     """
 
     def __init__(self, **kwargs):
-        r"""Initialize heat exchanger component with given parameters."""
+        r"""
+        Initialize the heat exchanger component.
+
+        Parameters
+        ----------
+        **kwargs : dict
+            Arbitrary keyword arguments. Recognized keys:
+            - dissipative (bool): whether component has dissipative behavior, default False
+            - Ex_C_col (dict): custom cost coefficients, default {}
+            - Z_costs (float): investment cost rate in currency/h, default 0.0
+        """
         self.dissipative = False
         super().__init__(**kwargs)
 
     def calc_exergy_balance(self, T0: float, p0: float, split_physical_exergy) -> None:
         r"""
-        Calculate the exergy balance of the heat exchanger.
+        Compute the exergy balance of the heat exchanger.
+
+        Case 1: All streams above ambient temperature
+
+        If `split_physical_exergy=True`:
+
+        .. math::
+
+            \dot{E}_{\mathrm{P}}
+            = \dot{E}^{\mathrm{T}}_{\mathrm{out},2}
+            - \dot{E}^{\mathrm{T}}_{\mathrm{in},2}
+
+        .. math::
+
+            \dot{E}_{\mathrm{F}}
+            = \dot{E}^{\mathrm{PH}}_{\mathrm{in},1}
+            - \dot{E}^{\mathrm{PH}}_{\mathrm{out},1}
+            + \bigl(\dot{E}^{\mathrm{M}}_{\mathrm{in},2}
+                    - \dot{E}^{\mathrm{M}}_{\mathrm{out},2}\bigr)
+
+        Else:
+
+        .. math::
+
+            \dot{E}_{\mathrm{P}}
+            = \dot{E}^{\mathrm{PH}}_{\mathrm{out},2}
+            - \dot{E}^{\mathrm{PH}}_{\mathrm{in},2}
+
+        .. math::
+
+            \dot{E}_{\mathrm{F}}
+            = \dot{E}^{\mathrm{PH}}_{\mathrm{in},1}
+            - \dot{E}^{\mathrm{PH}}_{\mathrm{out},1}
+
+        Case 2: All streams below or equal to ambient temperature
+
+        If `split_physical_exergy=True`:
+
+        .. math::
+            \dot{E}_{\mathrm{P}}
+            = \dot{E}^{\mathrm{T}}_{\mathrm{out},1}
+            - \dot{E}^{\mathrm{T}}_{\mathrm{in},1}
+
+        .. math::
+
+            \dot{E}_{\mathrm{F}}
+            = \dot{E}^{\mathrm{PH}}_{\mathrm{in},2}
+            - \dot{E}^{\mathrm{PH}}_{\mathrm{out},2}
+            + \bigl(\dot{E}^{\mathrm{M}}_{\mathrm{in},1}
+                    - \dot{E}^{\mathrm{M}}_{\mathrm{out},1}\bigr)
+
+        Else
+
+        .. math::
+            \dot{E}_{\mathrm{P}}
+            = \dot{E}^{\mathrm{PH}}_{\mathrm{out},1}
+            - \dot{E}^{\mathrm{PH}}_{\mathrm{in},1}
+
+        .. math::
+            \dot{E}_{\mathrm{F}}
+            = \dot{E}^{\mathrm{PH}}_{\mathrm{in},2}
+            - \dot{E}^{\mathrm{PH}}_{\mathrm{out},2}
+
+        Case 3: Both stream crossing ambient temperature
+
+        If `split_physical_exergy=True`:
+
+        .. math::
+            \dot{E}_{\mathrm{P}}
+            = \dot{E}^{\mathrm{T}}_{\mathrm{out},1}
+            + \dot{E}^{\mathrm{T}}_{\mathrm{out},2}
+
+        .. math::
+            \dot{E}_{\mathrm{F}}
+            = \dot{E}^{\mathrm{PH}}_{\mathrm{in},1}
+            + \dot{E}^{\mathrm{PH}}_{\mathrm{in},2}
+            - \bigl(\dot{E}^{\mathrm{M}}_{\mathrm{out},1}
+                    + \dot{E}^{\mathrm{M}}_{\mathrm{out},2}\bigr)
+
+        Else:
 
-        Performs exergy balance calculations considering different temperature regimes
-        relative to the ambient temperature. The method automatically determines the
-        appropriate case based on stream temperatures.
+        .. math::
+            \dot{E}_{\mathrm{P}}
+            = \dot{E}^{\mathrm{T}}_{\mathrm{out},1}
+            + \dot{E}^{\mathrm{T}}_{\mathrm{out},2}
+
+        .. math::
+            \dot{E}_{\mathrm{F}}
+            = \dot{E}^{\mathrm{PH}}_{\mathrm{in},1}
+            + \dot{E}^{\mathrm{PH}}_{\mathrm{in},2}
+
+        Case 4: Only the hot inlet above ambient temperature
+
+        If `split_physical_exergy=True`:
+
+        .. math::
+            \dot{E}_{\mathrm{P}}
+            = \dot{E}^{\mathrm{T}}_{\mathrm{out},1}
+
+        .. math::
+            \dot{E}_{\mathrm{F}}
+            = \bigl(\dot{E}^{\mathrm{PH}}_{\mathrm{in},1}
+                    + \dot{E}^{\mathrm{PH}}_{\mathrm{in},2}\bigr)
+            - \bigl(\dot{E}^{\mathrm{PH}}_{\mathrm{out},2}
+                    + \dot{E}^{\mathrm{M}}_{\mathrm{out},1}\bigr)
+
+        Else:
+
+        .. math::
+            \dot{E}_{\mathrm{P}}
+            = \dot{E}^{\mathrm{PH}}_{\mathrm{out},1}
+
+        .. math::
+            \dot{E}_{\mathrm{F}}
+            = \bigl(\dot{E}^{\mathrm{PH}}_{\mathrm{in},1}
+                    + \dot{E}^{\mathrm{PH}}_{\mathrm{in},2}\bigr)
+            - \dot{E}^{\mathrm{PH}}_{\mathrm{out},2}
+
+        Case 5: Only the cold inlet below ambient temperature
+
+        If `split_physical_exergy=True`:
+
+        .. math::
+            \dot{E}_{\mathrm{P}}
+            = \dot{E}^{\mathrm{T}}_{\mathrm{out},2}
+
+        .. math::
+            \dot{E}_{\mathrm{F}}
+            = \bigl(\dot{E}^{\mathrm{PH}}_{\mathrm{in},1}
+                    - \dot{E}^{\mathrm{PH}}_{\mathrm{out},1}\bigr)
+            + \bigl(\dot{E}^{\mathrm{PH}}_{\mathrm{in},2}
+                    - \dot{E}^{\mathrm{M}}_{\mathrm{out},2}\bigr)
+
+        Else:
+
+        .. math::
+            \dot{E}_{\mathrm{P}}
+            = \dot{E}^{\mathrm{PH}}_{\mathrm{out},2}
+
+        .. math::
+            \dot{E}_{\mathrm{F}}
+            = \bigl(\dot{E}^{\mathrm{PH}}_{\mathrm{in},1}
+                    - \dot{E}^{\mathrm{PH}}_{\mathrm{out},1}\bigr)
+            + \dot{E}^{\mathrm{PH}}_{\mathrm{in},2}
+
+        Case 6: Hot stream always above and cold stream always below ambiente temperature (dissipative case):
+
+        .. math::
+            \dot{E}_{\mathrm{P}} = \mathrm{NaN}
+
+        .. math::
+            \dot{E}_{\mathrm{F}}
+            = \bigl(\dot{E}^{\mathrm{PH}}_{\mathrm{in},1}
+                    - \dot{E}^{\mathrm{PH}}_{\mathrm{out},1}\bigr)
+            - \dot{E}^{\mathrm{PH}}_{\mathrm{out},2}
+            + \dot{E}^{\mathrm{PH}}_{\mathrm{in},2}
+
+        If `dissipative` is `True`, the component is treated as dissipative:
+
+        .. math::
+            \dot{E}_{\mathrm{P}} = \mathrm{NaN}
+
+        .. math::
+            \dot{E}_{\mathrm{F}}
+            = \bigl(\dot{E}^{\mathrm{PH}}_{\mathrm{in},1}
+                    - \dot{E}^{\mathrm{PH}}_{\mathrm{out},1}\bigr)
+            - \dot{E}^{\mathrm{PH}}_{\mathrm{out},2}
+            + \dot{E}^{\mathrm{PH}}_{\mathrm{in},2}
 
         Parameters
         ----------
         T0 : float
-            Ambient temperature in :math:`\mathrm{K}`.
+            Ambient temperature (K).
         p0 : float
-            Ambient pressure in :math:`\mathrm{Pa}`.
+            Ambient pressure (Pa).
         split_physical_exergy : bool
-            Flag indicating whether physical exergy is split into thermal and mechanical components.
+            Whether to split thermal and mechanical exergy.
 
         Raises
         ------
         ValueError
-            If the required inlet and outlet streams are not properly defined.
-
-        Notes
-        -----
-        This method updates the following component attributes:
-            - E_P (Exergy product)
-            - E_F (Exergy fuel)
-            - E_D (Exergy destruction)
-            - epsilon (Exergetic efficiency)
-
-        The calculation requires two inlet and two outlet streams, and their
-        temperature relationships with ambient temperature determine which case
-        of exergy analysis is applied.
+            If required inlets or outlets are missing.
         """
         # Ensure that the component has both inlet and outlet streams
         if len(self.inl) < 2 or len(self.outl) < 2:
@@ -169,80 +257,125 @@ class HeatExchanger(Component):
 
         if not self.dissipative:
             # Case 1: All streams are above the ambient temperature
-            if all([stream['T'] >= T0 for stream in all_streams]):
+            if all([stream["T"] >= T0 for stream in all_streams]):
                 if split_physical_exergy:
-                    self.E_P = self.outl[1]['m'] * self.outl[1]['e_T'] - self.inl[1]['m'] * self.inl[1]['e_T']
-                    self.E_F = self.inl[0]['m'] * self.inl[0]['e_PH'] - self.outl[0]['m'] * self.outl[0]['e_PH'] + (
-                        self.inl[1]['m'] * self.inl[1]['e_M'] - self.outl[1]['m'] * self.outl[1]['e_M'])
+                    self.E_P = self.outl[1]["m"] * self.outl[1]["e_T"] - self.inl[1]["m"] * self.inl[1]["e_T"]
+                    self.E_F = (
+                        self.inl[0]["m"] * self.inl[0]["e_PH"]
+                        - self.outl[0]["m"] * self.outl[0]["e_PH"]
+                        + (self.inl[1]["m"] * self.inl[1]["e_M"] - self.outl[1]["m"] * self.outl[1]["e_M"])
+                    )
                 else:
-                    self.E_P = self.outl[1]['m'] * self.outl[1]['e_PH'] - self.inl[1]['m'] * self.inl[1]['e_PH']
-                    self.E_F = self.inl[0]['m'] * self.inl[0]['e_PH'] - self.outl[0]['m'] * self.outl[0]['e_PH']
+                    self.E_P = self.outl[1]["m"] * self.outl[1]["e_PH"] - self.inl[1]["m"] * self.inl[1]["e_PH"]
+                    self.E_F = self.inl[0]["m"] * self.inl[0]["e_PH"] - self.outl[0]["m"] * self.outl[0]["e_PH"]
 
             # Case 2: All streams are below or equal to the ambient temperature
-            elif all([stream['T'] <= T0 for stream in all_streams]):
+            elif all([stream["T"] <= T0 for stream in all_streams]):
                 if split_physical_exergy:
-                    self.E_P = self.outl[0]['m'] * self.outl[0]['e_T'] - self.inl[0]['m'] * self.inl[0]['e_T']
-                    self.E_F = self.inl[1]['m'] * self.inl[1]['e_PH'] - self.outl[1]['m'] * self.outl[1]['e_PH'] + (
-                        self.inl[0]['m'] * self.inl[0]['e_M'] - self.outl[0]['m'] * self.outl[0]['e_M'])
+                    self.E_P = self.outl[0]["m"] * self.outl[0]["e_T"] - self.inl[0]["m"] * self.inl[0]["e_T"]
+                    self.E_F = (
+                        self.inl[1]["m"] * self.inl[1]["e_PH"]
+                        - self.outl[1]["m"] * self.outl[1]["e_PH"]
+                        + (self.inl[0]["m"] * self.inl[0]["e_M"] - self.outl[0]["m"] * self.outl[0]["e_M"])
+                    )
                 else:
-                    logging.warning("While dealing with heat exchnager below ambient temperautre, "
-                    "physical exergy should be split into thermal and mechanical components!")
-                    self.E_P = self.outl[0]['m'] * self.outl[0]['e_PH'] - self.inl[0]['m'] * self.inl[0]['e_PH']
-                    self.E_F = self.inl[1]['m'] * self.inl[1]['e_PH'] - self.outl[1]['m'] * self.outl[1]['e_PH']
-
-            # Case 3: Hot stream from above to lower ambient, cold stream from lower to above ambient
-            elif (self.inl[0]['T'] > T0 and self.outl[1]['T'] > T0 and
-                self.outl[0]['T'] <= T0 and self.inl[1]['T'] <= T0):
+                    logging.warning(
+                        "While dealing with heat exchnager below ambient temperautre, "
+                        "physical exergy should be split into thermal and mechanical components!"
+                    )
+                    self.E_P = self.outl[0]["m"] * self.outl[0]["e_PH"] - self.inl[0]["m"] * self.inl[0]["e_PH"]
+                    self.E_F = self.inl[1]["m"] * self.inl[1]["e_PH"] - self.outl[1]["m"] * self.outl[1]["e_PH"]
+
+            # Case 3: Both stream crossing T0 (hot inlet and cold outlet > T0, hot outlet and cold inlet <= T0)
+            elif (
+                self.inl[0]["T"] > T0 and self.outl[1]["T"] > T0 and self.outl[0]["T"] <= T0 and self.inl[1]["T"] <= T0
+            ):
                 if split_physical_exergy:
-                    self.E_P = self.outl[0]['m'] * self.outl[0]['e_T'] + self.outl[1]['m'] * self.outl[1]['e_T']
-                    self.E_F = self.inl[0]['m'] * self.inl[0]['e_PH'] + self.inl[1]['m'] * self.inl[1]['e_PH'] - (
-                        self.outl[0]['m'] * self.outl[0]['e_M'] + self.outl[1]['m'] * self.outl[1]['e_M'])
+                    self.E_P = self.outl[0]["m"] * self.outl[0]["e_T"] + self.outl[1]["m"] * self.outl[1]["e_T"]
+                    self.E_F = (
+                        self.inl[0]["m"] * self.inl[0]["e_PH"]
+                        + self.inl[1]["m"] * self.inl[1]["e_PH"]
+                        - (self.outl[0]["m"] * self.outl[0]["e_M"] + self.outl[1]["m"] * self.outl[1]["e_M"])
+                    )
                 else:
-                    logging.warning("While dealing with heat exchnager below ambient temperautre, "
-                    "physical exergy should be split into thermal and mechanical components!")
-                    self.E_P = self.outl[0]['m'] * self.outl[0]['e_PH'] + self.outl[1]['m'] * self.outl[1]['e_PH']
-                    self.E_F = self.inl[0]['m'] * self.inl[0]['e_PH'] + self.inl[1]['m'] * self.inl[1]['e_PH']
-
-            # Case 4: Hot stream inlet above ambient, all others below or equal to ambient
-            elif (self.inl[0]['T'] > T0 and self.inl[1]['T'] <= T0 and
-                self.outl[0]['T'] <= T0 and self.outl[1]['T'] <= T0):
+                    logging.warning(
+                        "While dealing with heat exchnager below ambient temperautre, "
+                        "physical exergy should be split into thermal and mechanical components!"
+                    )
+                    self.E_P = self.outl[0]["m"] * self.outl[0]["e_PH"] + self.outl[1]["m"] * self.outl[1]["e_PH"]
+                    self.E_F = self.inl[0]["m"] * self.inl[0]["e_PH"] + self.inl[1]["m"] * self.inl[1]["e_PH"]
+
+            # Case 4: Only hot inlet > T0
+            elif (
+                self.inl[0]["T"] > T0 and self.inl[1]["T"] <= T0 and self.outl[0]["T"] <= T0 and self.outl[1]["T"] <= T0
+            ):
                 if split_physical_exergy:
-                    self.E_P = self.outl[0]['m'] * self.outl[0]['e_T']
-                    self.E_F = self.inl[0]['m'] * self.inl[0]['e_PH'] + self.inl[1]['m'] * self.inl[1]['e_PH'] - (
-                        self.outl[1]['m'] * self.outl[1]['e_PH'] + self.outl[0]['m'] * self.outl[0]['e_M'])
+                    self.E_P = self.outl[0]["m"] * self.outl[0]["e_T"]
+                    self.E_F = (
+                        self.inl[0]["m"] * self.inl[0]["e_PH"]
+                        + self.inl[1]["m"] * self.inl[1]["e_PH"]
+                        - (self.outl[1]["m"] * self.outl[1]["e_PH"] + self.outl[0]["m"] * self.outl[0]["e_M"])
+                    )
                 else:
-                    logging.warning("While dealing with heat exchnager below ambient temperautre, "
-                    "physical exergy should be split into thermal and mechanical components!")
-                    self.E_P = self.outl[0]['m'] * self.outl[0]['e_PH']
-                    self.E_F = self.inl[0]['m'] * self.inl[0]['e_PH'] + (
-                        self.inl[1]['m'] * self.inl[1]['e_PH'] - self.outl[1]['m'] * self.outl[1]['e_PH'])
-
-            # Case 5: Inlets are higher but outlets are below or equal to ambient
-            elif (self.inl[0]['T'] > T0 and self.outl[0]['T'] > T0 and
-                self.inl[1]['T'] <= T0 and self.outl[1]['T'] <= T0):
-                self.E_P = np.nan
-                self.E_F = self.inl[0]['m'] * self.inl[0]['e_PH'] - self.outl[0]['m'] * self.outl[0]['e_PH'] + (
-                    self.inl[1]['m'] * self.inl[1]['e_PH'] - self.outl[1]['m'] * self.outl[1]['e_PH'])
-
-            # Case 6: Cold inlet is lower ambient, others higher
-            else:
+                    logging.warning(
+                        "While dealing with heat exchnager below ambient temperautre, "
+                        "physical exergy should be split into thermal and mechanical components!"
+                    )
+                    self.E_P = self.outl[0]["m"] * self.outl[0]["e_PH"]
+                    self.E_F = self.inl[0]["m"] * self.inl[0]["e_PH"] + (
+                        self.inl[1]["m"] * self.inl[1]["e_PH"] - self.outl[1]["m"] * self.outl[1]["e_PH"]
+                    )
+
+            # Case 5: Only cold inlet <= T0
+            elif self.inl[0]["T"] > T0 and self.inl[1]["T"] <= T0 and self.outl[0]["T"] > T0 and self.outl[1]["T"] > T0:
                 if split_physical_exergy:
-                    self.E_P = self.outl[1]['m'] * self.outl[1]['e_T']
-                    self.E_F = self.inl[0]['m'] * self.inl[0]['e_PH'] - self.outl[0]['m'] * self.outl[0]['e_PH'] + (
-                        self.inl[1]['m'] * self.inl[1]['e_PH'] - self.outl[1]['m'] * self.outl[1]['e_M'])
+                    self.E_P = self.outl[1]["m"] * self.outl[1]["e_T"]
+                    self.E_F = (
+                        self.inl[0]["m"] * self.inl[0]["e_PH"]
+                        - self.outl[0]["m"] * self.outl[0]["e_PH"]
+                        + (self.inl[1]["m"] * self.inl[1]["e_PH"] - self.outl[1]["m"] * self.outl[1]["e_M"])
+                    )
                 else:
-                    logging.warning("While dealing with heat exchnager below ambient temperautre, "
-                    "physical exergy should be split into thermal and mechanical components!")
-                    self.E_P = self.outl[1]['m'] * self.outl[1]['e_PH']
-                    self.E_F = self.inl[0]['m'] * self.inl[0]['e_PH'] - self.outl[0]['m'] * self.outl[0]['e_PH'] + (
-                        self.inl[1]['m'] * self.inl[1]['e_PH'])
+                    logging.warning(
+                        "While dealing with heat exchnager below ambient temperautre, "
+                        "physical exergy should be split into thermal and mechanical components!"
+                    )
+                    self.E_P = self.outl[1]["m"] * self.outl[1]["e_PH"]
+                    self.E_F = (
+                        self.inl[0]["m"] * self.inl[0]["e_PH"]
+                        - self.outl[0]["m"] * self.outl[0]["e_PH"]
+                        + (self.inl[1]["m"] * self.inl[1]["e_PH"])
+                    )
+
+            # Case 6: hot stream always above T0, cold stream always below T0 (dissipative case)
+            elif (
+                self.inl[0]["T"] > T0 and self.inl[1]["T"] <= T0 and self.outl[0]["T"] > T0 and self.outl[1]["T"] <= T0
+            ):
+                self.E_P = np.nan
+                self.E_F = (
+                    self.inl[0]["m"] * self.inl[0]["e_PH"]
+                    - self.outl[0]["m"] * self.outl[0]["e_PH"]
+                    + (self.inl[1]["m"] * self.inl[1]["e_PH"] - self.outl[1]["m"] * self.outl[1]["e_PH"])
+                )
+
+                logging.warning(
+                    f"Component {self.name} is dissipative. This component should be "
+                    "handled with the `dissipative` flag set to True."
+                )
+
+            # Case 7: Not implemented case
+            else:
+                logging.error(
+                    f"The heat exchanger {self.name} has an unexpected temperature configuration. "
+                    "Please check the inlet and outlet temperatures."
+                )
 
         else:
             self.E_F = (
-                self.inl[0]['m'] * self.inl[0]['e_PH']
-                - self.outl[0]['m'] * self.outl[0]['e_PH']
-                - self.outl[1]['m'] * self.outl[1]['e_PH']
-                + self.inl[1]['m'] * self.inl[1]['e_PH']
+                self.inl[0]["m"] * self.inl[0]["e_PH"]
+                - self.outl[0]["m"] * self.outl[0]["e_PH"]
+                - self.outl[1]["m"] * self.outl[1]["e_PH"]
+                + self.inl[1]["m"] * self.inl[1]["e_PH"]
             )
             self.E_P = np.nan
         # Calculate exergy destruction and efficiency
@@ -254,42 +387,114 @@ class HeatExchanger(Component):
 
         # Log the results
         logging.info(
-            f"HeatExchanger exergy balance calculated: "
+            f"Exergy balance of HeatExchanger {self.name} calculated: "
             f"E_P={self.E_P:.2f}, E_F={self.E_F:.2f}, E_D={self.E_D:.2f}, "
             f"Efficiency={self.epsilon:.2%}"
         )
 
-
     def aux_eqs(self, A, b, counter, T0, equations, chemical_exergy_enabled):
-        """Set up auxiliary equations for heat exchanger exergy cost analysis.
-        This method constructs auxiliary equations for thermal, mechanical, and chemical
-        exergy costs in a heat exchanger. It considers different thermal cases based on the
-        temperatures of the streams relative to the reference temperature T0.
+        r"""
+        Add auxiliary cost equations for the heat exchanger.
+
+        This method appends rows to the cost matrix to enforce:
+
+        Case 1: All streams above ambient temperature
+
+        F rule for thermal exergy of the hot stream:
+
+        .. math::
+            -\frac{1}{\dot{E}^{\mathrm{T}}_{\mathrm{out},1}}\,\dot{C}^{\mathrm{T}}_{\mathrm{out},1}
+            + \frac{1}{\dot{E}^{\mathrm{T}}_{\mathrm{in},1}}\,\dot{C}^{\mathrm{T}}_{\mathrm{in},1}
+            = 0
+
+        Case 2: All streams below or equal to ambient temperature
+
+        F rule for thermal exergy of the cold stream:
+
+        .. math::
+            -\frac{1}{\dot{E}^{\mathrm{T}}_{\mathrm{out},2}}\,\dot{C}^{\mathrm{T}}_{\mathrm{out},2}
+            + \frac{1}{\dot{E}^{\mathrm{T}}_{\mathrm{in},2}}\,\dot{C}^{\mathrm{T}}_{\mathrm{in},2}
+            = 0
+
+        Case 3: Both stream crossing ambient temperature
+
+        P rule for thermal exergy of both outlets:
+
+        .. math::
+            -\frac{1}{\dot{E}^{\mathrm{T}}_{\mathrm{out},1}}\,\dot{C}^{\mathrm{T}}_{\mathrm{out},1}
+            + \frac{1}{\dot{E}^{\mathrm{T}}_{\mathrm{out},2}}\,\dot{C}^{\mathrm{T}}_{\mathrm{out},2}
+            = 0
+
+        Case 4: Only the hot inlet above ambient temperature
+
+        F rule for thermal exergy of the cold stream:
+
+        .. math::
+            -\frac{1}{\dot{E}^{\mathrm{T}}_{\mathrm{out},2}}\,\dot{C}^{\mathrm{T}}_{\mathrm{out},2}
+            + \frac{1}{\dot{E}^{\mathrm{T}}_{\mathrm{in},2}}\,\dot{C}^{\mathrm{T}}_{\mathrm{in},2}
+            = 0
+
+        Case 5: Only the cold inlet below ambient temperature
+
+        F rule for thermal exergy of the hot stream:
+
+        .. math::
+            -\frac{1}{\dot{E}^{\mathrm{T}}_{\mathrm{out},1}}\,\dot{C}^{\mathrm{T}}_{\mathrm{out},1}
+            + \frac{1}{\dot{E}^{\mathrm{T}}_{\mathrm{in},1}}\,\dot{C}^{\mathrm{T}}_{\mathrm{in},1}
+            = 0
+
+        Case 6: Hot stream always above and cold stream always below ambiente temperature (dissipative case):
+
+        The dissipative is not handeld here!
+
+        For all cases, the mechanical and chemical exergy costs are handled as follows:
+
+        F rule for mechanical exergy of the hot stream:
+
+        .. math::
+            -\frac{1}{\dot{E}^{\mathrm{M}}_{\mathrm{out},i}}\,\dot{C}^{\mathrm{M}}_{\mathrm{out},i}
+            + \frac{1}{\dot{E}^{\mathrm{M}}_{\mathrm{in},i}}\,\dot{C}^{\mathrm{M}}_{\mathrm{in},i}
+            = 0
+
+        F rule for chemical exergy on hot branch:
+
+        .. math::
+            -\frac{1}{\dot{E}^{\mathrm{CH}}_{\mathrm{out},i}}\,\dot{C}^{\mathrm{CH}}_{\mathrm{out},i}
+            + \frac{1}{\dot{E}^{\mathrm{CH}}_{\mathrm{in},i}}\,\dot{C}^{\mathrm{CH}}_{\mathrm{in},i}
+            = 0
+
         Parameters
         ----------
         A : numpy.ndarray
-            Coefficient matrix for the linear equation system.
+            Current cost matrix.
         b : numpy.ndarray
-            Right-hand side vector of the linear equation system.
+            Current RHS vector.
         counter : int
-            Current row counter in the linear equation system.
+            Starting row index for auxiliary equations.
         T0 : float
-            Reference temperature for exergy calculations.
-        equations : dict
-            Dictionary to store equation descriptions.
+            Ambient temperature (K).
+        equations : dict or list
+            Structure for equation labels.
         chemical_exergy_enabled : bool
-            Flag indicating whether chemical exergy analysis is enabled.
+            Must be True to include chemical exergy mixing.
+
         Returns
         -------
         A : numpy.ndarray
-            Updated coefficient matrix.
+            Updated cost matrix.
         b : numpy.ndarray
-            Updated right-hand side vector.
+            Updated RHS vector.
         counter : int
-            Updated row counter.
-        equations : dict
-            Updated dictionary of equation descriptions.
+            Updated row index after adding equations.
+        equations : dict or list
+            Updated labels.
+
+        Raises
+        ------
+        ValueError
+            If required cost variable indices are missing.
         """
+
         # Equality equation for mechanical and chemical exergy costs.
         def set_equal(A, row, in_item, out_item, var):
             if in_item["e_" + var] != 0 and out_item["e_" + var] != 0:
@@ -346,48 +551,85 @@ class HeatExchanger(Component):
         # Case 1: All temperatures > T0.
         if all([c["T"] > T0 for c in list(self.inl.values()) + list(self.outl.values())]):
             set_thermal_f_hot(A, counter + 0)
-            equations[counter] = f"aux_f_rule_hot_{self.name}"
+            equations[counter] = {
+                "kind": "aux_f_rule_hot",
+                "objects": [self.name, self.inl[0]["name"], self.outl[0]["name"]],
+                "property": "c_T",
+            }
         # Case 2: All temperatures <= T0.
         elif all([c["T"] <= T0 for c in list(self.inl.values()) + list(self.outl.values())]):
             set_thermal_f_cold(A, counter + 0)
-            equations[counter] = f"aux_f_rule_cold_{self.name}"
-        # Case 3: Mixed temperatures: inl[0]["T"] > T0 and outl[1]["T"] > T0, while outl[0]["T"] <= T0 and inl[1]["T"] <= T0.
-        elif (self.inl[0]["T"] > T0 and self.outl[1]["T"] > T0 and
-            self.outl[0]["T"] <= T0 and self.inl[1]["T"] <= T0):
+            equations[counter] = {
+                "kind": "aux_f_rule_cold",
+                "objects": [self.name, self.inl[1]["name"], self.outl[1]["name"]],
+                "property": "c_T",
+            }
+        # Case 3: Both stream crossing T0 (hot inlet and cold outlet > T0, hot outlet and cold inlet <= T0)
+        elif self.inl[0]["T"] > T0 and self.outl[1]["T"] > T0 and self.outl[0]["T"] <= T0 and self.inl[1]["T"] <= T0:
             set_thermal_p_rule(A, counter + 0)
-            equations[counter] = f"aux_p_rule_{self.name}"
-        # Case 4: Mixed temperatures: inl[0]["T"] > T0, inl[1]["T"] <= T0, and both outl[0]["T"] and outl[1]["T"] <= T0.
-        elif (self.inl[0]["T"] > T0 and self.inl[1]["T"] <= T0 and
-            self.outl[0]["T"] <= T0 and self.outl[1]["T"] <= T0):
+            equations[counter] = {
+                "kind": "aux_p_rule",
+                "objects": [self.name, self.outl[0]["name"], self.outl[1]["name"]],
+                "property": "c_T",
+            }
+        # Case 4: Only hot inlet > T0
+        elif self.inl[0]["T"] > T0 and self.inl[1]["T"] <= T0 and self.outl[0]["T"] <= T0 and self.outl[1]["T"] <= T0:
             set_thermal_f_cold(A, counter + 0)
-            equations[counter] = f"aux_f_rule_cold_{self.name}"
-        # Case 5: Mixed temperatures: inl[0]["T"] > T0, inl[1]["T"] <= T0, and both outl[0]["T"] and outl[1]["T"] > T0.
-        elif (self.inl[0]["T"] > T0 and self.inl[1]["T"] <= T0 and
-            self.outl[0]["T"] > T0 and self.outl[1]["T"] > T0):
+            equations[counter] = {
+                "kind": "aux_f_rule_cold",
+                "objects": [self.name, self.inl[1]["name"], self.outl[1]["name"]],
+                "property": "c_T",
+            }
+        # Case 5: Only cold inlet <= T0
+        elif self.inl[0]["T"] > T0 and self.inl[1]["T"] <= T0 and self.outl[0]["T"] > T0 and self.outl[1]["T"] > T0:
             set_thermal_f_hot(A, counter + 0)
-            equations[counter] = f"aux_f_rule_hot_{self.name}"
-        # Case 6: Mixed temperatures (dissipative case): inl[0]["T"] > T0, inl[1]["T"] <= T0, outl[0]["T"] > T0, and outl[1]["T"] <= T0.
-        elif (self.inl[0]["T"] > T0 and self.inl[1]["T"] <= T0 and
-            self.outl[0]["T"] > T0 and self.outl[1]["T"] <= T0):
-            print("you shouldn't see this")
+            equations[counter] = {
+                "kind": "aux_f_rule_hot",
+                "objects": [self.name, self.inl[0]["name"], self.outl[0]["name"]],
+                "property": "c_T",
+            }
+        # Case 6: hot stream always above T0, cold stream always below T0 (dissipative case)
+        elif self.inl[0]["T"] > T0 and self.inl[1]["T"] <= T0 and self.outl[0]["T"] > T0 and self.outl[1]["T"] <= T0:
+            logging.warning(
+                f"Component {self.name} is dissipative. This component should be "
+                "handled with the `dissipative` flag set to True."
+            )
             return
-        # Case 7: Default case.
+        # Case 7: Not implemented case
         else:
-            set_thermal_f_hot(A, counter + 0)
-            equations[counter] = f"aux_f_rule_hot_{self.name}"
-        
+            logging.error(
+                f"The heat exchanger {self.name} has an unexpected temperature configuration. "
+                "Please check the inlet and outlet temperatures."
+            )
+
         # Mechanical equations (always added)
         set_equal(A, counter + 1, self.inl[0], self.outl[0], "M")
         set_equal(A, counter + 2, self.inl[1], self.outl[1], "M")
-        equations[counter + 1] = f"aux_equality_mech_{self.outl[0]['name']}"
-        equations[counter + 2] = f"aux_equality_mech_{self.outl[1]['name']}"
-        
+        equations[counter + 1] = {
+            "kind": "aux_equality",
+            "objects": [self.name, self.inl[0]["name"], self.outl[0]["name"]],
+            "property": "c_M",
+        }
+        equations[counter + 2] = {
+            "kind": "aux_equality",
+            "objects": [self.name, self.inl[1]["name"], self.outl[1]["name"]],
+            "property": "c_M",
+        }
+
         # Only add chemical auxiliary equations if chemical exergy is enabled.
         if chemical_exergy_enabled:
             set_equal(A, counter + 3, self.inl[0], self.outl[0], "CH")
             set_equal(A, counter + 4, self.inl[1], self.outl[1], "CH")
-            equations[counter + 3] = f"aux_equality_chem_{self.outl[0]['name']}"
-            equations[counter + 4] = f"aux_equality_chem_{self.outl[1]['name']}"
+            equations[counter + 3] = {
+                "kind": "aux_equality",
+                "objects": [self.name, self.inl[0]["name"], self.outl[0]["name"]],
+                "property": "c_CH",
+            }
+            equations[counter + 4] = {
+                "kind": "aux_equality",
+                "objects": [self.name, self.inl[1]["name"], self.outl[1]["name"]],
+                "property": "c_M",
+            }
             num_aux_eqs = 5
         else:
             # Skip chemical auxiliary equations.
@@ -398,52 +640,147 @@ class HeatExchanger(Component):
 
         return A, b, counter + num_aux_eqs, equations
 
-    def exergoeconomic_balance(self, T0):
-        """
-        Perform exergoeconomic balance calculations for the general heat exchanger.
-        
-        This method calculates various exergoeconomic parameters including:
-        - Cost rates of product (C_P) and fuel (C_F)
-        - Specific cost of product (c_P) and fuel (c_F)
-        - Cost rate of exergy destruction (C_D)
-        - Relative cost difference (r)
-        - Exergoeconomic factor (f)
-        
+    def exergoeconomic_balance(self, T0, chemical_exergy_enabled=False):
+        r"""
+        Perform exergoeconomic cost balance for the heat exchanger.
+
+        .. math::
+            \dot{C}^{\mathrm{T}}_{\mathrm{in},1}
+            + \dot{C}^{\mathrm{M}}_{\mathrm{in},1}
+            + \dot{C}^{\mathrm{T}}_{\mathrm{in},2}
+            + \dot{C}^{\mathrm{M}}_{\mathrm{in},2}
+            - \dot{C}^{\mathrm{T}}_{\mathrm{out},1}
+            - \dot{C}^{\mathrm{M}}_{\mathrm{out},1}
+            - \dot{C}^{\mathrm{T}}_{\mathrm{out},2}
+            - \dot{C}^{\mathrm{M}}_{\mathrm{out},2}
+            + \dot{Z}
+            = 0
+
+        In case the chemical exergy of the streams is know:
+
+        .. math::
+            \dot{C}^{\mathrm{CH}}_{\mathrm{in},1} =
+            \dot{C}^{\mathrm{CH}}_{\mathrm{out},1}
+
+        .. math::
+            \dot{C}^{\mathrm{CH}}_{\mathrm{in},2} =
+            \dot{C}^{\mathrm{CH}}_{\mathrm{out},2}
+
+        This method computes cost coefficients and ratios:
+
+        Case 1: All streams above ambient temperature
+
+        .. math::
+            \dot{C}_P = \dot{C}^{\mathrm{T}}_{\mathrm{out},2}
+                    - \dot{C}^{\mathrm{T}}_{\mathrm{in},2}
+
+        .. math::
+            \dot{C}_F = \dot{C}^{\mathrm{PH}}_{\mathrm{in},1}
+                    - \dot{C}^{\mathrm{PH}}_{\mathrm{out},1}
+                    + \bigl(\dot{C}^{\mathrm{M}}_{\mathrm{in},2}
+                            - \dot{C}^{\mathrm{M}}_{\mathrm{out},2}\bigr)
+
+        Case 2: All streams below or equal to ambient temperature
+
+        .. math::
+            \dot{C}_P = \dot{C}^{\mathrm{T}}_{\mathrm{out},1}
+                    - \dot{C}^{\mathrm{T}}_{\mathrm{in},1}
+
+        .. math::
+            \dot{C}_F = \dot{C}^{\mathrm{PH}}_{\mathrm{in},2}
+                    - \dot{C}^{\mathrm{PH}}_{\mathrm{out},2}
+                    + \bigl(\dot{C}^{\mathrm{M}}_{\mathrm{in},1}
+                            - \dot{C}^{\mathrm{M}}_{\mathrm{out},1}\bigr)
+
+        Case 3: Both stream crossing ambient temperature
+
+        .. math::
+            \dot{C}_P = \dot{C}^{\mathrm{T}}_{\mathrm{out},1}
+                    + \dot{C}^{\mathrm{T}}_{\mathrm{out},2}
+
+        .. math::
+            \dot{C}_F = \dot{C}^{\mathrm{PH}}_{\mathrm{in},1}
+                    + \dot{C}^{\mathrm{PH}}_{\mathrm{in},2}
+                    - \bigl(\dot{C}^{\mathrm{M}}_{\mathrm{out},1}
+                            + \dot{C}^{\mathrm{M}}_{\mathrm{out},2}\bigr)
+
+        Case 4: Only the hot inlet above ambient temperature
+
+        .. math::
+            \dot{C}_P = \dot{C}^{\mathrm{T}}_{\mathrm{out},1}
+
+        .. math::
+            \dot{C}_F = \bigl(\dot{C}^{\mathrm{PH}}_{\mathrm{in},1}
+                            + \dot{C}^{\mathrm{PH}}_{\mathrm{in},2}\bigr)
+                    - \bigl(\dot{C}^{\mathrm{PH}}_{\mathrm{out},2}
+                            + \dot{C}^{\mathrm{M}}_{\mathrm{out},1}\bigr)
+
+        Case 5: Only the cold inlet below ambient temperature
+
+        .. math::
+            \dot{C}_P = \dot{C}^{\mathrm{T}}_{\mathrm{out},2}
+
+        .. math::
+            \dot{C}_F = \dot{C}^{\mathrm{PH}}_{\mathrm{in},1}
+                    - \dot{C}^{\mathrm{PH}}_{\mathrm{out},1}
+                    + \bigl(\dot{C}^{\mathrm{PH}}_{\mathrm{in},2}
+                            - \dot{C}^{\mathrm{M}}_{\mathrm{out},2}\bigr)
+
+        Case 6: Hot stream always above and cold stream always below ambient temperature (dissipative case):
+
+        .. math::
+            \dot{C}_P = \mathrm{NaN}
+
+        .. math::
+            \dot{C}_F = \bigl(\dot{C}^{\mathrm{PH}}_{\mathrm{in},1}
+                    - \dot{C}^{\mathrm{PH}}_{\mathrm{out},1}\bigr)
+            - \dot{C}^{\mathrm{PH}}_{\mathrm{out},2}
+            + \dot{C}^{\mathrm{PH}}_{\mathrm{in},2}
+
         Parameters
         ----------
         T0 : float
-            Ambient temperature
-            
-        Notes
-        -----
-        The exergoeconomic balance considers thermal (T), chemical (CH),
-        and mechanical (M) exergy components for the inlet and outlet streams.
+            Ambient temperature (K).
+        chemical_exergy_enabled : bool, optional
+            If True, chemical exergy is considered in the calculations.
         """
+        # Case 1: All streams are above the ambient temperature
         if all([c["T"] > T0 for c in list(self.inl.values()) + list(self.outl.values())]):
             self.C_P = self.outl[1]["C_T"] - self.inl[1]["C_T"]
-            self.C_F = self.inl[0]["C_PH"] - self.outl[0]["C_PH"] + (
-                self.inl[1]["C_M"] - self.outl[1]["C_M"])
+            self.C_F = self.inl[0]["C_PH"] - self.outl[0]["C_PH"] + (self.inl[1]["C_M"] - self.outl[1]["C_M"])
+        # Case 2: All streams are below or equal to the ambient temperature
         elif all([c["T"] <= T0 for c in list(self.inl.values()) + list(self.outl.values())]):
             self.C_P = self.outl[0]["C_T"] - self.inl[0]["C_T"]
-            self.C_F = self.inl[1]["C_PH"] - self.outl[1]["C_PH"] + (
-                self.inl[0]["C_M"] - self.outl[0]["C_M"])
-        elif (self.inl[0]["T"] > T0 and self.outl[1]["T"] > T0 and
-              self.outl[0]["T"] <= T0 and self.inl[1]["T"] <= T0):
+            self.C_F = self.inl[1]["C_PH"] - self.outl[1]["C_PH"] + (self.inl[0]["C_M"] - self.outl[0]["C_M"])
+        # Case 3: Both stream crossing T0 (hot inlet and cold outlet > T0, hot outlet and cold inlet <= T0)
+        elif self.inl[0]["T"] > T0 and self.outl[1]["T"] > T0 and self.outl[0]["T"] <= T0 and self.inl[1]["T"] <= T0:
             self.C_P = self.outl[0]["C_T"] + self.outl[1]["C_T"]
-            self.C_F = self.inl[0]["C_PH"] + self.inl[1]["C_PH"] - (
-                self.outl[0]["C_M"] + self.outl[1]["C_M"])
-        elif (self.inl[0]["T"] > T0 and self.inl[1]["T"] <= T0 and
-              self.outl[0]["T"] <= T0 and self.outl[1]["T"] <= T0):
+            self.C_F = self.inl[0]["C_PH"] + self.inl[1]["C_PH"] - (self.outl[0]["C_M"] + self.outl[1]["C_M"])
+        # Case 4: Only hot inlet > T0
+        elif self.inl[0]["T"] > T0 and self.inl[1]["T"] <= T0 and self.outl[0]["T"] <= T0 and self.outl[1]["T"] <= T0:
             self.C_P = self.outl[0]["C_T"]
-            self.C_F = self.inl[0]["C_PH"] + self.inl[1]["C_PH"] - (
-               self.outl[1]["C_PH"] + self.outl[0]["C_M"])
-        else:
+            self.C_F = self.inl[0]["C_PH"] + self.inl[1]["C_PH"] - (self.outl[1]["C_PH"] + self.outl[0]["C_M"])
+        # Case 5: Only cold inlet <= T0
+        elif self.inl[0]["T"] > T0 and self.inl[1]["T"] <= T0 and self.outl[0]["T"] > T0 and self.outl[1]["T"] > T0:
             self.C_P = self.outl[1]["C_T"]
-            self.C_F = self.inl[0]["C_PH"] - self.outl[0]["C_PH"] + (
-                self.inl[1]["C_PH"] - self.outl[1]["C_M"])
+            self.C_F = self.inl[0]["C_PH"] - self.outl[0]["C_PH"] + (self.inl[1]["C_PH"] - self.outl[1]["C_M"])
+        # Case 6: hot stream always above T0, cold stream always below T0 (dissipative case)
+        elif self.inl[0]["T"] > T0 and self.inl[1]["T"] <= T0 and self.outl[0]["T"] > T0 and self.outl[1]["T"] <= T0:
+            logging.warning(
+                f"Component {self.name} is dissipative. This component should be "
+                "handled with the `dissipative` flag set to True."
+            )
+            self.C_P = np.nan
+            self.C_F = self.inl[0]["C_PH"] - self.outl[0]["C_PH"] + (self.inl[1]["C_PH"] - self.outl[1]["C_PH"])
+        # Case 7: Not implemented case
+        else:
+            logging.error(
+                f"The heat exchanger {self.name} has an unexpected temperature configuration. "
+                "Please check the inlet and outlet temperatures."
+            )
 
         self.c_F = self.C_F / self.E_F
         self.c_P = self.C_P / self.E_P
         self.C_D = self.c_F * self.E_D
         self.r = (self.c_P - self.c_F) / self.c_F
-        self.f = self.Z_costs / (self.Z_costs + self.C_D)
+        self.f = self.Z_costs / (self.Z_costs + self.C_D)