Merge branch 'fix/pygmo-implementation' into dev

docs/whats_new.rst

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

docs/whats_new/v0-7-1.rst

@@ -0,0 +1,14 @@
+v0.7.1 - Newton's Nature (December, 2, 2023)
+++++++++++++++++++++++++++++++++++++++++++++
+
+Bug Fixes
+#########
+- Several bugs introduced by the restructuring of the package in version 0.7.0
+  have been fixed:
+
+  - `PR #451 <https://github.com/oemof/tespy/pull/451>`__
+  - `PR #453 <https://github.com/oemof/tespy/pull/453>`__
+
+Contributors
+############
+- Francesco Witte (`@fwitte <https://github.com/fwitte>`__)

pyproject.toml

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

src/tespy/__init__.py

@@ -3,7 +3,7 @@
 import os
 
 __datapath__ = os.path.join(importlib.resources.files("tespy"), "data")
-__version__ = '0.7.1 - Newton\'s Nature'
+__version__ = '0.7.1.post1 - Newton\'s Nature'
 
 # tespy data and connections import
 from . import connections  # noqa: F401

src/tespy/components/heat_exchangers/desuperheater.py

@@ -239,6 +239,6 @@ def saturated_gas_deriv(self, increment_filter, k):
         """
         o = self.outl[0]
         if self.is_variable(o.p):
-            self.jacobian[k, o.p.J_col] = -dh_mix_dpQ(o, 1, o.fluid_data)
+            self.jacobian[k, o.p.J_col] = -dh_mix_dpQ(o.p.val_SI, 1, o.fluid_data)
         if self.is_variable(o.h):
             self.jacobian[k, o.h.J_col] = 1

src/tespy/networks/network.py

@@ -847,6 +847,21 @@ def initialise(self):
         self.num_conn_vars = 0
         self.variables_dict = {}
 
+        # in multiprocessing copies are made of all connections
+        # the mass flow branches and fluid branches hold references to
+        # connections from the original run (where network.checked is False)
+        # The assignment of variable spaces etc. is however made on the
+        # copies of the connections which do not correspond to the mass flow
+        # branches and fluid branches anymore. So the topology simplification
+        # does not actually apply to the copied network, therefore the
+        # branches have to be recreated for this case. We can detect that by
+        # checking whether a network holds a massflow branch with some
+        # connections and compare that with the connection object actually
+        # present in the network
+        first_conn = self.massflow_branches[0]["connections"][0]
+        if self.conns.loc[first_conn.label, "object"] != first_conn:
+            self.create_massflow_and_fluid_branches()
+            self.create_fluid_wrapper_branches()
         self.propagate_fluid_wrappers()
         self.presolve_massflow_topology()
         self.presolve_fluid_topology()

src/tespy/tools/fluid_properties/wrappers.py

@@ -16,6 +16,16 @@
 from tespy.tools.global_vars import ERR
 
 
+class SerializableAbstractState(CP.AbstractState):
+
+    def __init__(self, back_end, fluid_name):
+        self.back_end = back_end
+        self.fluid_name = fluid_name
+
+    def __reduce__(self):
+        return (self.__class__, (self.back_end, self.fluid_name))
+
+
 class FluidPropertyWrapper:
 
     def __init__(self, fluid, back_end=None) -> None:
@@ -112,7 +122,7 @@ def __init__(self, fluid, back_end=None) -> None:
             back_end = "HEOS"
 
         super().__init__(fluid, back_end)
-        self.AS = CP.CoolProp.AbstractState(self.back_end, self.fluid)
+        self.AS = SerializableAbstractState(self.back_end, self.fluid)
         self._set_constants()
 
     def _set_constants(self):

tests/test_components/test_combustion.py

@@ -12,8 +12,6 @@
 
 import shutil
 
-import numpy as np
-
 from tespy.components import CombustionChamber
 from tespy.components import CombustionEngine
 from tespy.components import DiabaticCombustionChamber

tests/test_components/test_merge.py

@@ -119,17 +119,3 @@ def test_two_fluid_setup(self):
         target = c1.m.val_SI
         msg = f"Target value for mass flow at connection 3 is {target}"
         assert c6.m.val_SI == approx(target), msg
-
-
-test = TestMerge()
-test.setup_method()
-test.test_single_fluid_at_outlet()
-test.setup_method()
-test.test_massflows_from_two_fluid_fractions()
-
-
-test2 = TestCyclicMerging()
-test2.setup_method()
-test2.test_single_fluid_setup()
-test2.setup_method()
-test2.test_two_fluid_setup()

tests/test_components/test_reactors.py

@@ -12,8 +12,6 @@
 
 import shutil
 
-import numpy as np
-
 from tespy.components import Sink
 from tespy.components import Source
 from tespy.components import WaterElectrolyzer

tutorial/advanced/optimization_example.py

@@ -10,7 +10,6 @@
 from tespy.components import SimpleHeatExchanger
 from tespy.components import Merge
 from tespy.components import Splitter
-from tespy.components import Valve
 from tespy.components import Pump
 from tespy.components import Turbine
 from tespy.connections import Bus
@@ -82,10 +81,7 @@ def __init__(self):
         c32 = Connection(fwh1, "out1", pu3, "in1", label="32")
         c33 = Connection(pu3, "out1", me, "in2", label="33")
 
-        self.nw.add_conns(
-            c21, c22, c23, c24,
-            c31, c32, c33
-        )
+        self.nw.add_conns(c21, c22, c23, c24, c31, c32, c33)
 
         # cooling water
         c41 = Connection(cwi, "out1", con, "in2", label="41")
@@ -108,8 +104,6 @@ def __init__(self):
 
         self.nw.add_busses(self.power, self.heat)
 
-        # parametrization
-        # components
         hpt.set_attr(eta_s=0.9)
         mpt.set_attr(eta_s=0.9)
         lpt.set_attr(eta_s=0.9)
@@ -129,13 +123,16 @@ def __init__(self):
         c2.set_attr(p=20)
         c4.set_attr(p=3)
 
-        c41.set_attr(T=20, p=3, fluid={"water": 1})
-        c42.set_attr(T=28)
+        c41.set_attr(T=20, p=3, fluid={"INCOMP::Water": 1})
+        c42.set_attr(T=28, p0=3, h0=100)
 
+        # parametrization
+        # components
         self.nw.solve("design")
         self.stable = "_stable"
         self.nw.save(self.stable)
         self.solved = True
+        self.nw.print_results()
 
     # %%[sec_2]