tespy 0.8.0.post4__tar.gz → 0.8.2__tar.gz

{tespy-0.8.0.post4 → tespy-0.8.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: tespy
-Version: 0.8.0.post4
+Version: 0.8.2
 Summary: Thermal Engineering Systems in Python (TESPy)
 Author-email: Francesco Witte <tespy@witte.sh>
 Requires-Python: >=3.9

{tespy-0.8.0.post4 → tespy-0.8.2}/docs/api/components.rst

@@ -78,6 +78,22 @@ tespy.components.combustion.engine module
     :undoc-members:
     :show-inheritance:
 
+tespy.components.heat_exchangers.base module
+--------------------------------------------
+
+.. automodule:: tespy.components.heat_exchangers.base
+    :members:
+    :undoc-members:
+    :show-inheritance:
+
+tespy.components.heat_exchangers.movingboundary module
+------------------------------------------------------
+
+.. automodule:: tespy.components.heat_exchangers.movingboundary
+    :members:
+    :undoc-members:
+    :show-inheritance:
+
 tespy.components.heat_exchangers.condenser module
 -------------------------------------------------
 
@@ -94,14 +110,6 @@ tespy.components.heat_exchangers.desuperheater module
     :undoc-members:
     :show-inheritance:
 
-tespy.components.heat_exchangers.base module
---------------------------------------------
-
-.. automodule:: tespy.components.heat_exchangers.base
-    :members:
-    :undoc-members:
-    :show-inheritance:
-
 tespy.components.heat_exchangers.simple module
 ----------------------------------------------
 

{tespy-0.8.0.post4 → tespy-0.8.2}/docs/modules/components.rst

@@ -777,12 +777,14 @@ Custom subsystems
 -----------------
 
 Create a :code:`.py` file in your working-directory. This file contains the
-class definition of your subsystem and at minimum two methods:
+class definition of your subsystem and at minimum one method:
 
-- :code:`create_comps`: Method to create the components of your subsystem and
-  save them in the :code:`Subsystem.comps` attribute (dictionary).
-- :code:`create_conns`: Method to create the connections of your subsystem and
-  save them in the :code:`Subsystem.conns` attribute (dictionary).
+- :code:`create_network`: Method to create the network of your subsystem.
+
+On top of that you need to add methods to define the available interfaces of
+your subsystem to the remaining network through specifying the number of inlets
+and outlets in the :code:`__init__` method of your class as seen in the code
+example below.
 
 All other functionalities are inherited by the parent class of the
 :py:class:`subsystem <tespy.components.subsystem.Subsystem>` object.
@@ -826,28 +828,33 @@ Create a file, e.g. :code:`mysubsystems.py` and add the following lines:
 
     >>> class WasteHeatSteamGenerator(Subsystem):
     ...     """Class documentation"""
+    ...     def __init__(self, label):
+    ...         self.num_in = 2
+    ...         self.num_out = 2
+    ...         super().__init__(label)
     ...
-    ...     def create_comps(self):
-    ...         """Create the subsystem's components."""
-    ...         self.comps['eco'] = HeatExchanger(f'{self.label}_economizer')
-    ...         self.comps['eva'] = HeatExchanger(f'{self.label}_evaporator')
-    ...         self.comps['sup'] = HeatExchanger(f'{self.label}_superheater')
-    ...         self.comps['drum'] = Drum(f'{self.label}_drum')
-    ...
-    ...     def create_conns(self):
+    ...     def create_network(self):
     ...         """Define the subsystem's connections."""
-    ...         self.conns['eco_dr'] = Connection(
-    ...             self.comps['eco'], 'out2', self.comps['drum'], 'in1')
-    ...         self.conns['dr_eva'] = Connection(
-    ...             self.comps['drum'], 'out1', self.comps['eva'], 'in2')
-    ...         self.conns['eva_dr'] = Connection(
-    ...             self.comps['eva'], 'out2', self.comps['drum'], 'in2')
-    ...         self.conns['dr_sup'] = Connection(
-    ...             self.comps['drum'], 'out2', self.comps['sup'], 'in2')
-    ...         self.conns['sup_eva'] = Connection(
-    ...             self.comps['sup'], 'out1', self.comps['eva'], 'in1')
-    ...         self.conns['eva_eco'] = Connection(
-    ...             self.comps['eva'], 'out1', self.comps['eco'], 'in1')
+    ...         eco = HeatExchanger('economizer')
+    ...         eva = HeatExchanger('evaporator')
+    ...         sup = HeatExchanger('superheater')
+    ...         drum = Drum('drum')
+    ...
+    ...         inlet_eco = Connection(self.inlet, 'out2', eco, 'in2', label='1')
+    ...         eco_dr = Connection(eco, 'out2', drum, 'in1', label='2')
+    ...         dr_eva = Connection(drum, 'out1', eva, 'in2', label='3')
+    ...         eva_dr = Connection(eva, 'out2', drum, 'in2', label='4')
+    ...         dr_sup = Connection(drum, 'out2', sup, 'in2', label='5')
+    ...         sup_outlet = Connection(sup, 'out2', self.outlet, 'in2', label='6')
+    ...
+    ...         self.add_conns(inlet_eco, eco_dr, dr_eva, eva_dr, dr_sup, sup_outlet)
+    ...
+    ...         inlet_sup = Connection(self.inlet, 'out1', sup, 'in1', label='11')
+    ...         sup_eva = Connection(sup, 'out1', eva, 'in1', label='12')
+    ...         eva_eco = Connection(eva, 'out1', eco, 'in1', label='13')
+    ...         eco_outlet = Connection(eco, 'out1', self.outlet, 'in1', label='14')
+    ...
+    ...         self.add_conns(inlet_sup, sup_eva, eva_eco, eco_outlet)
 
 .. note::
 
@@ -881,13 +888,13 @@ different tespy classes required.
     >>> sg = WasteHeatSteamGenerator('waste heat steam generator')
 
     >>> # %% connection definition
-    >>> fw_sg = Connection(feed_water, 'out1', sg.comps['eco'], 'in2')
-    >>> sg_ls = Connection(sg.comps['sup'], 'out2', steam, 'in1')
-    >>> fg_sg = Connection(waste_heat, 'out1', sg.comps['sup'], 'in1')
-    >>> sg_ch = Connection(sg.comps['eco'], 'out1', chimney, 'in1')
+    >>> fw_sg = Connection(feed_water, 'out1', sg, 'in2')
+    >>> sg_ls = Connection(sg, 'out2', steam, 'in1')
+    >>> fg_sg = Connection(waste_heat, 'out1', sg, 'in1')
+    >>> sg_ch = Connection(sg, 'out1', chimney, 'in1')
 
     >>> nw.add_conns(fw_sg, sg_ls, fg_sg, sg_ch)
-    >>> nw.add_subsys(sg)
+    >>> nw.add_subsystems(sg)
 
     >>> # %% connection parameters
     >>> fw_sg.set_attr(fluid={'water': 1}, T=25)
@@ -895,29 +902,30 @@ different tespy classes required.
     >>> sg_ls.set_attr(p=130)
     >>> sg_ch.set_attr(p=1)
 
-    >>> sg.conns['eva_dr'].set_attr(x=0.6)
+    >>> sg.get_conn('4').set_attr(x=0.6)
 
     >>> # %% component parameters
-    >>> sg.comps['eco'].set_attr(
+    >>> sg.get_comp('economizer').set_attr(
     ...     pr1=0.999,  pr2=0.97, design=['pr1', 'pr2', 'ttd_u'],
     ...     offdesign=['zeta1', 'zeta2', 'kA_char']
     ... )
 
-    >>> sg.comps['eva'].set_attr(
+    >>> sg.get_comp('evaporator').set_attr(
     ...     pr1=0.999, ttd_l=20, design=['pr1', 'ttd_l'],
     ...     offdesign=['zeta1', 'kA_char']
     ... )
 
-    >>> sg.comps['sup'].set_attr(
+    >>> sg.get_comp('superheater').set_attr(
     ...     pr1=0.999,  pr2=0.99, ttd_u=50, design=['pr1', 'pr2', 'ttd_u'],
     ...     offdesign=['zeta1', 'zeta2', 'kA_char']
     ... )
 
-    >>> sg.conns['eco_dr'].set_attr(Td_bp=-5, design=['Td_bp'])
+    >>> sg.get_conn('2').set_attr(Td_bp=-5, design=['Td_bp'])
 
     >>> # %% solve
     >>> # solve design case
     >>> nw.solve('design')
+    >>> nw.assert_convergence()
     >>> nw.save('tmp.json')
 
     >>> # offdesign test

{tespy-0.8.0.post4 → tespy-0.8.2}/docs/modules/connections.rst

@@ -132,7 +132,7 @@ example:
     >>> myconn.set_attr(fluid={'water': 1})  # HEOS back end
     >>> myconn.set_attr(fluid={'INCOMP::water': 1})  # incompressible fluid
     >>> myconn.set_attr(fluid={'BICUBIC::air': 1})  # bicubic back end
-    >>> myconn.set_attr(fluid={'INCOMP::MPG[0.5]': 1})  # binary incompressible mixture
+    >>> myconn.set_attr(fluid={'INCOMP::MPG[0.5]|mass': 1})  # binary incompressible mixture
 
 .. note::
 

{tespy-0.8.0.post4 → tespy-0.8.2}/docs/modules/networks.rst

@@ -95,7 +95,7 @@ or via subsystems using the corresponding methods:
 .. code-block:: python
 
     >>> my_plant.add_conns()
-    >>> my_plant.add_subsys()
+    >>> my_plant.add_subsystems()
 
 .. note::
 
@@ -417,7 +417,7 @@ power :math:`P` to be 1000 W, the set of equations will look like this:
     &0 = \dot{m}_{in} - \dot{m}_{out}\\
     \mathrm{additional:} \, &0 = 1000 - \dot{m}_{in} (\cdot {h_{out} - h_{in}})
 
-.. _tespy_modules_convergence_check_label:
+.. _tespy_module_convergence_label:
 
 Convergence stability
 +++++++++++++++++++++
@@ -462,9 +462,25 @@ after the third iteration, further checks are usually not required.
 .. tip::
 
     To check if the solver successfully found a solution for your model you can
-    check the `.converged` attribute of the Network class after calling the
-    `solve` method. It will be `True` in case no linear dependency was and the
-    residual value of all equations is below the minimum threshold.
+    check the `.status` attribute of the Network class after calling the
+    `solve` method. It will be
+
+    - 0 in case a solution was successfully found
+    - 1 in case a solution was found, but some parameters violate physical
+      limits
+    - 2 in case no convergence was achieved after completion of the iterations
+    - 3 in case a linear dependency in the Jacobian matrix is found
+    - 11 in case the number of specified parameters is too small for the given
+      problem
+    - 12 in case the number of specified parameters is too large for the given
+      problem
+    - 99 in case the simulation crashed due to any other reason
+
+    The :code:`solve` does not exit with an exception in case the status is
+    0, 1, 2 or 3. If you want to raise an error in your script, you can call
+    the `Network.assert_convergence()` method. It will raise an
+    :code:`AssertionError` if the simulation did not find a converged solution
+    status 2 or status 3.
 
 Calculation speed improvement
 +++++++++++++++++++++++++++++

{tespy-0.8.0.post4 → tespy-0.8.2}/docs/whats_new/v0-8-0-001.rst

@@ -9,7 +9,7 @@ Bug Fixes
   (`PR #640 <https://github.com/oemof/tespy/pull/640>`__).
 - In case a heat exchanger heats up the same mass flow (e.g. in a recuperator
   or internal heat exchanger), the derivative to the mass flow was assigned
-  wrong (`PR #640 <https://github.com/oemof/tespy/pull/640>`__).
+  wrong (`PR #646 <https://github.com/oemof/tespy/pull/646>`__).
 
 Other Changes
 #############

{tespy-0.8.0.post4 → tespy-0.8.2}/docs/whats_new/v0-8-0.rst

@@ -5,6 +5,7 @@ API breaking changes
 ####################
 - The :code:`OptimizationProblem` implements the following changes
   (see `PR #610 <https://github.com/oemof/tespy/pull/610>`__):
+
   - The dtype for :code:`objective` has been changed to list. With a single
     element list, a single objective optimization is carried out. With lists
     containing more elements, multi-objective optimization will be carried out.

tespy-0.8.2/docs/whats_new/v0-8-1.rst

@@ -0,0 +1,31 @@
+v0.8.1 - Newton's Nature (May, 29, 2025)
+++++++++++++++++++++++++++++++++++++++++
+
+New Features
+############
+- The :code:`CoolPropWrapper` API can now handle CoolProp based mixtures
+  and allows the user to specify which type of mixture fractions to use
+  (mass, molar, volume). This is mostly for REFPROP support and will require
+  further adaptions for the :code:`CoolPropWrapper` in the future, because the
+  pure fluid functions may not be applicable in the context of mixtures
+  (`PR #655 <https://github.com/oemof/tespy/pull/655>`__).
+
+  In case you are working with incompressible mixtures, this feature is API
+  breaking. You need to append :code:`|mass`, :code:`|volume` or
+  :code:`|molar` to your fluid name string at the very end to specify, which
+  type of mixture is used. This information can be retrieved from the
+  CoolProp online documentation on the incompressible fluids.
+
+Bug Fixes
+#########
+- In case parts of a :code:`Network` were missing a :code:`CycleCloser` or
+  a :code:`Source` and :code:`Sink` combination a misleading error was raised.
+  These cases are checked for explicitly now to make debugging easier
+  (`PR #653 <https://github.com/oemof/tespy/pull/653>`__).
+- The logging of warnings for :code:`HeatExchanger` based components includes
+  the label of the component, when cold or hot side effectiveness cannot be
+  calculated (`PR #670 <https://github.com/oemof/tespy/pull/670>`__).
+
+Contributors
+############
+- Francesco Witte (`@fwitte <https://github.com/fwitte>`__)

tespy-0.8.2/docs/whats_new/v0-8-2.rst

@@ -0,0 +1,50 @@
+v0.8.2 - Newton's Nature (June, 12, 2025)
++++++++++++++++++++++++++++++++++++++++++
+
+API changes
+###########
+- The API of the :code:`Subsystem` class has been revised to make it act more
+  like a :code:`Network`. You can get components and connections from the
+  subsystem with respective methods, and the :code:`SubsystemInterface` class
+  is utilized to target the :code:`subsystem.inlet` and
+  :code:`subsystem.outlet` in a :code:`Connection` object. For the required
+  changes please check the respective section
+  :ref:`in the docs <tespy_subsystems_label>`
+  (`PR #652 <https://github.com/oemof/tespy/pull/652>`__).
+- The attribute :code:`progress` of the class :code:`Network` has been removed.
+  (`PR #684 <https://github.com/oemof/tespy/pull/684>`__).
+- To raise an :code:`AssertionError` on non-convergence of a simulation you
+  must now use :code:`Network.assert_convergence()` instead of
+  :code:`Network._convergence_check()`
+  (`PR #684 <https://github.com/oemof/tespy/pull/684>`__).
+
+New Features
+############
+- The combustion based component classes :code:`CombustionChamber`,
+  :code:`DiabaticCombustionChamber` and :code:`CombustionEngine` can now handle
+  carbonmonoxide as fuel
+  (`PR #674 <https://github.com/oemof/tespy/pull/674>`__).
+- The :code:`Network`'s :code:`solve` method will now assign status values
+  depending of the outcome of the simulation. Check the section on
+  :ref:`convergence <tespy_module_convergence_label>` for more information
+  (`PR #684 <https://github.com/oemof/tespy/pull/684>`__).
+
+Other Changes
+#############
+- The partial derivatives for specified temperature are only calculated towards
+  enthalpy and pressure, not towards the fluid compostion. The reason for this
+  is, that it is not intended, that the composition of a fluid can be
+  determined by specifying temperature. Removing this saves a lot of
+  computational overhead for mixtures
+  (`PR #674 <https://github.com/oemof/tespy/pull/674>`__).
+- The calculation of the logarithmic temperature difference in heat exchanger
+  classes was error-prone due to rounding differences in temperature. The
+  calculation is now consistent through all respective classes and can handle
+  temperature differences to a precision of :code:`1e-6`
+  (`PR #679 <https://github.com/oemof/tespy/pull/679>`__).
+- The export method :code:`to_exerpy` now includes the results for components
+  (`PR #680 <https://github.com/oemof/tespy/pull/680>`__).
+
+Contributors
+############
+- Francesco Witte (`@fwitte <https://github.com/fwitte>`__)

{tespy-0.8.0.post4 → tespy-0.8.2}/docs/whats_new.rst

@@ -3,6 +3,8 @@ What's New
 
 Discover notable new features and improvements in each release
 
+.. include::  whats_new/v0-8-2.rst
+.. include::  whats_new/v0-8-1.rst
 .. include::  whats_new/v0-8-0-001.rst
 .. include::  whats_new/v0-8-0.rst
 .. include::  whats_new/v0-7-9.rst

{tespy-0.8.0.post4 → tespy-0.8.2}/pyproject.toml

@@ -24,7 +24,7 @@ exclude = ["docs/_build"]
 
 [project]
 name = "tespy"
-version = "0.8.0.post4"
+version = "0.8.2"
 description = "Thermal Engineering Systems in Python (TESPy)"
 readme = "README.rst"
 authors = [

{tespy-0.8.0.post4 → tespy-0.8.2}/src/tespy/__init__.py

@@ -1,9 +1,18 @@
 # -*- coding: utf-8
 import importlib.resources
 import os
+import sys
+import warnings
+
+if sys.version_info[1] < 10:
+    msg = (
+        "Supprt for python versions below 3.10 will be dropped with the "
+        "next major release"
+    )
+    warnings.warn(FutureWarning(msg))
 
 __datapath__ = os.path.join(importlib.resources.files("tespy"), "data")
-__version__ = '0.8.0.post4 - Newton\'s Nature'
+__version__ = '0.8.2 - Newton\'s Nature'
 
 # tespy data and connections import
 from . import connections  # noqa: F401

{tespy-0.8.0.post4 → tespy-0.8.2}/src/tespy/components/combustion/base.py

@@ -47,7 +47,7 @@ class CombustionChamber(Component):
 
     Available fuels
 
-    - methane, ethane, propane, butane, hydrogen
+    - methane, ethane, propane, butane, hydrogen, carbon monoxide, nDodecane
 
     Inlets/Outlets
 
@@ -334,21 +334,22 @@ class CombustionChamber(Component):
         hf['propane'] = -103.8
         hf['butane'] = -125.7
         hf['nDodecane'] = -289.4
+        hf['CO'] = -110.5
         hf[self.o2] = 0
         hf[self.co2] = -393.51
         # water (gaseous)
         hf[self.h2o] = -241.826
 
         key = set(list(hf.keys())).intersection(
-                set([a.replace(' ', '')
-                     for a in CP.get_aliases(f)]))
+            set([a.replace(' ', '') for a in CP.get_aliases(f)])
+        )
 
         val = (
             -(
                 self.fuels[f]['H'] / 2 * hf[self.h2o]
                 + self.fuels[f]['C'] * hf[self.co2]
                 - (
-                    (self.fuels[f]['C'] + self.fuels[f]['H'] / 4) * hf[self.o2]
+                    (self.fuels[f]['C'] + self.fuels[f]['H'] / 4 - self.fuels[f]['O'] / 2) * hf[self.o2]
                     + hf[list(key)[0]]
                 )
             ) / inl[0].fluid.wrapper[f]._molar_mass * 1000
@@ -601,6 +602,7 @@ class CombustionChamber(Component):
         n_oxy_stoich = {}
         n_h = 0
         n_c = 0
+        n_o = 0
         for f in self.fuel_list:
             n_fuel[f] = 0
             for i in inl:
@@ -608,10 +610,12 @@ class CombustionChamber(Component):
                 n_fuel[f] += n
                 n_h += n * self.fuels[f]['H']
                 n_c += n * self.fuels[f]['C']
+                n_o += n * self.fuels[f]['O']
 
             # stoichiometric oxygen requirement for each fuel
             n_oxy_stoich[f] = n_fuel[f] * (
                 self.fuels[f]['H'] / 4 + self.fuels[f]['C']
+                - self.fuels[f]['O'] / 2
             )
 
         n_oxygen = 0
@@ -622,7 +626,7 @@ class CombustionChamber(Component):
 
         ###################################################################
         # calculate stoichiometric oxygen
-        n_oxygen_stoich = n_h / 4 + n_c
+        n_oxygen_stoich = n_h / 4 + n_c - n_o / 2
 
         ###################################################################
         # calculate lambda if not set
@@ -662,7 +666,7 @@ class CombustionChamber(Component):
             if self.lamb.val < 1:
                 n_fuel_exc = (
                     -(n_oxygen / n_oxygen_stoich - 1) * n_oxy_stoich[fluid]
-                    / (self.fuels[fluid]['H'] / 4 + self.fuels[fluid]['C'])
+                    / (self.fuels[fluid]['H'] / 4 + self.fuels[fluid]['C'] - self.fuels[fluid]['O'] / 2)
                 )
             else:
                 n_fuel_exc = 0
@@ -994,6 +998,7 @@ class CombustionChamber(Component):
 
         n_h = 0
         n_c = 0
+        n_o = 0
         for f in self.fuel_list:
             for i in inl:
                 n_fuel = (
@@ -1002,6 +1007,7 @@ class CombustionChamber(Component):
                 )
                 n_h += n_fuel * self.fuels[f]['H']
                 n_c += n_fuel * self.fuels[f]['C']
+                n_o += n_fuel * self.fuels[f]['O']
 
         n_oxygen = 0
         for i in inl:
@@ -1009,7 +1015,7 @@ class CombustionChamber(Component):
                 i.m.val_SI * i.fluid.val[self.o2]
                 / inl[0].fluid.wrapper[self.o2]._molar_mass
             )
-        n_oxygen_stoich = n_h / 4 + n_c
+        n_oxygen_stoich = n_h / 4 + n_c - n_o / 2
         return n_oxygen / n_oxygen_stoich
 
 
@@ -1218,16 +1224,12 @@ class CombustionChamber(Component):
                     outl.fluid.val[f] = 0.05
 
             elif f == self.co2:
-                if outl.fluid.val[f] > 0.1:
-                    outl.fluid.val[f] = 0.075
-                if outl.fluid.val[f] < 0.001:
-                    outl.fluid.val[f] = 0.02
+                if outl.fluid.val[f] > 0.15:
+                    outl.fluid.val[f] = 0.1
 
             elif f == self.h2o:
-                if outl.fluid.val[f] > 0.1:
-                    outl.fluid.val[f] = 0.075
-                if outl.fluid.val[f] < 0.001:
-                    outl.fluid.val[f] = 0.02
+                if outl.fluid.val[f] > 0.15:
+                    outl.fluid.val[f] = 0.1
 
             elif f in self.fuel_list:
                 if outl.fluid.val[f] > 0:

{tespy-0.8.0.post4 → tespy-0.8.2}/src/tespy/components/combustion/diabatic.py

@@ -42,7 +42,7 @@ class DiabaticCombustionChamber(CombustionChamber):
 
     Available fuels
 
-    - methane, ethane, propane, butane, hydrogen
+    - methane, ethane, propane, butane, hydrogen, carbon monoxide, nDodecane
 
     Inlets/Outlets
 

{tespy-0.8.0.post4 → tespy-0.8.2}/src/tespy/components/combustion/engine.py

@@ -63,7 +63,7 @@ class CombustionEngine(CombustionChamber):
 
     Available fuels
 
-    - methane, ethane, propane, butane, hydrogen
+    - methane, ethane, propane, butane, hydrogen, carbon monoxide, nDodecane
 
     Inlets/Outlets
 

{tespy-0.8.0.post4 → tespy-0.8.2}/src/tespy/components/component.py

@@ -838,6 +838,7 @@ class Component:
 
     def check_parameter_bounds(self):
         r"""Check parameter value limits."""
+        _no_limit_violated = True
         for p in self.parameters.keys():
             data = self.get_attr(p)
             if isinstance(data, dc_cp):
@@ -848,6 +849,7 @@ class Component:
                         f"{self.label}."
                     )
                     logger.warning(msg)
+                    _no_limit_violated = False
 
                 elif data.val < data.min_val - ERR:
                     msg = (
@@ -856,6 +858,7 @@ class Component:
                         f"{self.label}."
                     )
                     logger.warning(msg)
+                    _no_limit_violated = False
 
             elif isinstance(data, dc_cc) and data.is_set:
                 expr = self.get_char_expr(data.param, **data.char_params)
@@ -868,6 +871,8 @@ class Component:
                         char_data.param, **char_data.char_params)
                     char_data.char_func.get_domain_errors(expr, self.label)
 
+        return _no_limit_violated
+
     def convergence_check(self):
         return
 

{tespy-0.8.0.post4 → tespy-0.8.2}/src/tespy/components/heat_exchangers/base.py

@@ -494,7 +494,7 @@ class HeatExchanger(Component):
         ttd_u = T_i1 - T_o2
         ttd_l = T_o1 - T_i2
 
-        if ttd_u == ttd_l:
+        if round(ttd_u, 6) == round(ttd_l, 6):
             td_log = ttd_l
         else:
             td_log = (ttd_l - ttd_u) / math.log((ttd_l) / (ttd_u))
@@ -719,7 +719,7 @@ class HeatExchanger(Component):
 
     def ttd_l_func(self):
         r"""
-        Equation for upper terminal temperature difference.
+        Equation for lower terminal temperature difference.
 
         Returns
         -------
@@ -738,7 +738,7 @@ class HeatExchanger(Component):
 
     def ttd_l_func_doc(self, label):
         r"""
-        Equation for upper terminal temperature difference.
+        Equation for lower terminal temperature difference.
 
         Parameters
         ----------
@@ -755,7 +755,7 @@ class HeatExchanger(Component):
 
     def ttd_l_deriv(self, increment_filter, k):
         """
-        Calculate partial derivates of upper terminal temperature function.
+        Calculate partial derivates of lower terminal temperature function.
 
         Parameters
         ----------
@@ -774,7 +774,7 @@ class HeatExchanger(Component):
 
     def ttd_min_func(self):
         r"""
-        Equation for upper terminal temperature difference.
+        Equation for minimum terminal temperature difference.
 
         Returns
         -------
@@ -1171,7 +1171,7 @@ class HeatExchanger(Component):
         # kA and logarithmic temperature difference
         if self.ttd_u.val < 0 or self.ttd_l.val < 0:
             self.td_log.val = np.nan
-        elif self.ttd_l.val == self.ttd_u.val:
+        elif round(self.ttd_l.val, 6) == round(self.ttd_u.val, 6):
             self.td_log.val = self.ttd_l.val
         else:
             self.td_log.val = (
@@ -1189,7 +1189,7 @@ class HeatExchanger(Component):
         except ValueError:
             self.eff_hot.val = np.nan
             msg = (
-                "Cannot calculate heat exchanger hot side effectiveness "
+                f"Cannot calculate {self.label} hot side effectiveness "
                 "because cold side inlet temperature is out of bounds for hot "
                 "side fluid."
             )
@@ -1202,7 +1202,7 @@ class HeatExchanger(Component):
         except ValueError:
             self.eff_cold.val = np.nan
             msg = (
-                "Cannot calculate heat exchanger cold side effectiveness "
+                f"Cannot calculate {self.label} cold side effectiveness "
                 "because hot side inlet temperature is out of bounds for cold "
                 "side fluid."
             )

{tespy-0.8.0.post4 → tespy-0.8.2}/src/tespy/components/heat_exchangers/condenser.py

@@ -336,7 +336,7 @@ class Condenser(HeatExchanger):
 
         ttd_u = T_i1 - T_o2
         ttd_l = T_o1 - T_i2
-        if ttd_u == ttd_l:
+        if round(ttd_u, 6) == round(ttd_l, 6):
             td_log = ttd_l
         else:
             td_log = (ttd_l - ttd_u) / math.log((ttd_l) / (ttd_u))
@@ -494,7 +494,7 @@ class Condenser(HeatExchanger):
         # kA and logarithmic temperature difference
         if self.ttd_u.val < 0 or self.ttd_l.val < 0:
             self.td_log.val = np.nan
-        elif self.ttd_l.val == self.ttd_u.val:
+        elif round(self.ttd_l.val, 6) == round(self.ttd_u.val, 6):
             self.td_log.val = self.ttd_l.val
         else:
             self.td_log.val = (
@@ -512,7 +512,7 @@ class Condenser(HeatExchanger):
         except ValueError:
             self.eff_hot.val = np.nan
             msg = (
-                "Cannot calculate heat exchanger hot side effectiveness "
+                f"Cannot calculate {self.label} hot side effectiveness "
                 "because cold side inlet temperature is out of bounds for hot "
                 "side fluid."
             )
@@ -525,7 +525,7 @@ class Condenser(HeatExchanger):
         except ValueError:
             self.eff_cold.val = np.nan
             msg = (
-                "Cannot calculate heat exchanger cold side effectiveness "
+                f"Cannot calculate {self.label} cold side effectiveness "
                 "because hot side inlet temperature is out of bounds for cold "
                 "side fluid."
             )