tespy 0.10.0.dev0__tar.gz → 0.10.1__tar.gz

{tespy-0.10.0.dev0 → tespy-0.10.1}/PKG-INFO

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

{tespy-0.10.0.dev0 → tespy-0.10.1}/docs/advanced_features/fluid_properties.rst

@@ -94,6 +94,32 @@ recommend using the REFPROP back end instead of HEOS. Also note, the CoolProp
 mixture back end is not tested thoroughly. Please reach out if you would like
 to support us in adopting the TESPy implementation.
 
+CoolProp/REFPROP mixtures (zeotropic fluids)
+++++++++++++++++++++++++++++++++++++++++++++
+Zeotropic mixtures such as refrigerant blends (e.g. isobutane/isopentane for
+ORC applications, R410A, R404A) can be used via the CoolProp/REFPROP backend.
+The fluid name follows the pattern:
+
+.. code-block:: text
+
+    <BACKEND>::<FLUID1>[<x1>]&<FLUID2>[<x2>]|<fraction_type>
+
+where :code:`<BACKEND>` is :code:`REFPROP` or :code:`HEOS`, the fraction
+values :code:`x1`, :code:`x2`, ... sum to 1, and :code:`<fraction_type>` is
+either :code:`mass` (mass fractions) or :code:`molar` (mole fractions).
+REFPROP fluid names must match the REFPROP fluid file names exactly (e.g.
+:code:`ISOBUTAN`, not :code:`ISOBUTANE`), for example:
+
+.. code-block:: text
+
+    REFPROP::ISOBUTAN[0.5]&IPENTANE[0.5]|mass
+
+To use this in a connection, specify the full name as the fluid key:
+
+.. code-block:: python
+
+    c.set_attr(fluid={"REFPROP::ISOBUTAN[0.5]&IPENTANE[0.5]|mass": 1})
+
 .. _incompressible_wrapper_label:
 
 IncompressibleFluidWrapper
@@ -296,6 +322,7 @@ isentropic change of pressure for an ideal gas.
     ...
     ...         self.coefficients = COEF[fluid]
     ...         self.h_ref = self._h_pT(None, reference_temperature)
+    ...         self.s_ref = self._s_pT(None, reference_temperature)
     ...         self._molar_mass = self.coefficients[-1] * 1e-3
     ...         self._T_min = 100
     ...         self._T_max = 2000
@@ -326,9 +353,31 @@ isentropic change of pressure for an ideal gas.
     ...             + self.coefficients[5] / 3 * y ** 3
     ...         ) / self.coefficients[6]
     ...
+    ...     def _s_pT(self, p, T):
+    ...         y = T * 1e-3
+    ...         return 1e3 * (
+    ...             self.coefficients[1]
+    ...             + self.coefficients[2] * np.log(y)
+    ...             + self.coefficients[3] * y
+    ...             - self.coefficients[4] / (2 * y ** 2)
+    ...             + self.coefficients[5] / 2 * y ** 2
+    ...         ) / self.coefficients[6]
+    ...
+    ...     def s_pT(self, p, T):
+    ...         return self._s_pT(p, T) - self.s_ref
+    ...
     ...     def T_ph(self, p, h):
     ...         return newton(self.h_pT, self.cp_pT, h, p)
     ...
+    ...     def T_ps(self, p, s):
+    ...         return newton(self.s_pT, lambda p, T: self.cp_pT(p, T) / T, s, p)
+    ...
+    ...     def s_ph(self, p, h):
+    ...         return self.s_pT(p, self.T_ph(p, h))
+    ...
+    ...     def h_ps(self, p, s):
+    ...         return self.h_pT(p, self.T_ps(p, s))
+    ...
     ...     def isentropic(self, p_1, h_1, p_2):
     ...         T_1 = self.T_ph(p_1, h_1)
     ...         cp = self.cp_pT(p_1, T_1)
@@ -382,6 +431,12 @@ them to CoolProp.
     400.0
     >>> round(h)
     189769
+    >>> s = kkh_water.s_pT(1e5, 400)
+    >>> T = kkh_water.T_ps(1e5, s)
+    >>> float(round(T, 1))
+    400.0
+    >>> round(float(s))
+    546
 
     >>> from tespy.tools.fluid_properties import CoolPropWrapper
 

{tespy-0.10.0.dev0 → tespy-0.10.1}/docs/advanced_tutorials/heat_exchangers.ipynb

@@ -168,8 +168,8 @@
     "})\n",
     "\n",
     "heat = [0, abs(hx_0d.Q.val)]\n",
-    "T_hot = [c2.T.val_SI, c1.T.val_SI]\n",
-    "T_cold = [d1.T.val_SI, d2.T.val_SI]\n",
+    "T_hot = [c2.T.val, c1.T.val]\n",
+    "T_cold = [d1.T.val, d2.T.val]\n",
     "\n",
     "fig, ax = plt.subplots(1, figsize=(10, 6))\n",
     "\n",
@@ -198,7 +198,7 @@
     "ax.plot([heat[0] - offset - bar_len/2, heat[0] - offset + bar_len/2], [T_cold[0], T_cold[0]], color=annotation_color, lw=2)\n",
     "\n",
     "ax.set_xbound([-5.5 * offset, heat[-1] + 4 * offset])\n",
-    "ax.set_ylabel(\"temperature in K\")\n",
+    "ax.set_ylabel(\"temperature in °C\")\n",
     "ax.set_xlabel(\"heat transferred in MW\")\n",
     "plt.show()"
    ]
@@ -233,11 +233,12 @@
     "    \"UA (kW/K)\": f\"{hx_cond.UA.val:.1f}\",\n",
     "})\n",
     "\n",
-    "T_cond = c2.T.val_SI - c2.calc_td_dew()\n",
-    "heat_to_cond = c1.m.val_SI * (h_mix_pQ(c2.p.val_SI, 1, c2.fluid_data) - c2.h.val_SI) / 1e6\n",
+    "T_cond = c2.T.val - c2.calc_td_dew()\n",
+    "h_dew = h_mix_pQ(c2.p.val_SI, 1, c2.fluid_data)\n",
+    "heat_to_cond = abs(hx_cond.Q.val) * (h_dew - c2.h.val_SI) / (c1.h.val_SI - c2.h.val_SI)\n",
     "heat = [0, heat_to_cond, abs(hx_cond.Q.val)]\n",
-    "T_hot = [c2.T.val_SI, T_cond, c1.T.val_SI]\n",
-    "T_cold = [d1.T.val_SI, d2.T.val_SI]\n",
+    "T_hot = [c2.T.val, T_cond, c1.T.val]\n",
+    "T_cold = [d1.T.val, d2.T.val]\n",
     "\n",
     "fig, ax = plt.subplots(1, figsize=(10, 6))\n",
     "ax.plot(heat, T_hot, \"o-\", color=\"red\")\n",
@@ -267,7 +268,7 @@
     "ax.plot([heat[0] - offset - bar_len/2, heat[0] - offset + bar_len/2], [T_cold[0], T_cold[0]], color=annotation_color, lw=2)\n",
     "\n",
     "ax.set_xbound([-5.5 * offset, heat[-1] + 4 * offset])\n",
-    "ax.set_ylabel(\"temperature in K\")\n",
+    "ax.set_ylabel(\"temperature in °C\")\n",
     "ax.set_xlabel(\"heat transferred in MW\")\n",
     "plt.show()"
    ]
@@ -294,8 +295,9 @@
     "nw, c1, c2, d1, d2 = make_network(hx_mb)\n",
     "nw.solve(\"design\")\n",
     "\n",
-    "heat, T_hot, T_cold, _, _ = hx_mb.calc_sections()\n",
-    "heat /= 1e6\n",
+    "heat = hx_mb.Q_sections.val\n",
+    "T_hot = hx_mb.T_hot_sections.val\n",
+    "T_cold = hx_mb.T_cold_sections.val\n",
     "\n",
     "results.append({\n",
     "    \"model\": \"MovingBoundaryHeatExchanger\",\n",
@@ -331,7 +333,7 @@
     "ax.plot([heat[0] - offset - bar_len/2, heat[0] - offset + bar_len/2], [T_cold[0], T_cold[0]], color=annotation_color, lw=2)\n",
     "\n",
     "ax.set_xbound([-5.5 * offset, heat[-1] + 4 * offset])\n",
-    "ax.set_ylabel(\"temperature in K\")\n",
+    "ax.set_ylabel(\"temperature in °C\")\n",
     "ax.set_xlabel(\"heat transferred in MW\")\n",
     "plt.show()"
    ]
@@ -359,8 +361,9 @@
     "nw_sectioned, c1, c2, d1, d2 = make_network(hx_sectioned)\n",
     "nw_sectioned.solve(\"design\")\n",
     "\n",
-    "heat_50secs, T_hot_50secs, T_cold_50secs, _, _ = hx_sectioned.calc_sections()\n",
-    "heat_50secs /= 1e6\n",
+    "heat_50secs = hx_sectioned.Q_sections.val\n",
+    "T_hot_50secs = hx_sectioned.T_hot_sections.val\n",
+    "T_cold_50secs = hx_sectioned.T_cold_sections.val\n",
     "\n",
     "results.append({\n",
     "    \"model\": \"SectionedHeatExchanger (50 sections)\",\n",
@@ -395,7 +398,7 @@
     "ax.plot([heat_50secs[0] - offset - bar_len/2, heat_50secs[0] - offset + bar_len/2], [T_cold_50secs[0], T_cold_50secs[0]], color=annotation_color, lw=2)\n",
     "\n",
     "ax.set_xbound([-5.5 * offset, heat_50secs[-1] + 4 * offset])\n",
-    "ax.set_ylabel(\"temperature in K\")\n",
+    "ax.set_ylabel(\"temperature in °C\")\n",
     "ax.set_xlabel(\"heat transferred in MW\")\n",
     "plt.show()"
    ]
@@ -409,8 +412,9 @@
    "source": [
     "hx_sectioned.set_attr(num_sections=6)\n",
     "nw_sectioned.solve(\"design\")\n",
-    "heat_6secs, T_hot_6secs, T_cold_6secs, _, _ = hx_sectioned.calc_sections()\n",
-    "heat_6secs /= 1e6\n",
+    "heat_6secs = hx_sectioned.Q_sections.val\n",
+    "T_hot_6secs = hx_sectioned.T_hot_sections.val\n",
+    "T_cold_6secs = hx_sectioned.T_cold_sections.val\n",
     "\n",
     "results.append({\n",
     "    \"model\": \"SectionedHeatExchanger (6 sections)\",\n",
@@ -463,15 +467,16 @@
     "hx_sectioned.set_attr(num_sections=1)\n",
     "nw_sectioned.solve(\"design\")\n",
     "\n",
-    "heat_1sec, T_hot_1sec, T_cold_1sec, _, _ = hx_sectioned.calc_sections()\n",
-    "heat_1sec /= 1e6\n",
+    "heat_1sec = hx_sectioned.Q_sections.val\n",
+    "T_hot_1sec = hx_sectioned.T_hot_sections.val\n",
+    "T_cold_1sec = hx_sectioned.T_cold_sections.val\n",
     "\n",
     "fig, ax = plt.subplots(1, figsize=(10, 6))\n",
     "ax.plot((heat_1sec, heat_1sec), (list(T_hot_1sec), list(T_cold_1sec)), color=annotation_color, linestyle=\"-\", linewidth=0.5)\n",
     "ax.plot(heat_50secs, T_hot_50secs, \"o-\", color=annotation_color, markersize=2, linewidth=0.5, label=\"SectionedHeatExchanger (50 sections)\")\n",
     "ax.plot(heat_1sec, T_hot_1sec, \"o-\", color=\"red\", label=\"0D heat exchanger\")\n",
     "ax.plot(heat_1sec, T_cold_1sec, \"o-\", color=\"blue\")\n",
-    "ax.set_ylabel(\"temperature in K\")\n",
+    "ax.set_ylabel(\"temperature in °C\")\n",
     "ax.set_xlabel(\"heat transferred in MW\")\n",
     "ax.legend()\n",
     "plt.show()"
@@ -498,7 +503,7 @@
     "ax.plot(heat_50secs, T_hot_50secs, \"o-\", color=annotation_color, markersize=2, linewidth=0.5, label=\"SectionedHeatExchanger (50 sections)\")\n",
     "ax.plot(heat_6secs, T_hot_6secs, \"o-\", color=\"red\", label=\"SectionedHeatExchanger (6 sections)\")\n",
     "ax.plot(heat_6secs, T_cold_6secs, \"o-\", color=\"blue\")\n",
-    "ax.set_ylabel(\"temperature in K\")\n",
+    "ax.set_ylabel(\"temperature in °C\")\n",
     "ax.set_xlabel(\"heat transferred in MW\")\n",
     "ax.legend()\n",
     "plt.show()"
@@ -526,22 +531,24 @@
     "nw_sectioned, c1, c2, d1, d2 = make_network(hx_sectioned)\n",
     "hx_sectioned.set_attr(dp1=2, dp2=0.5)\n",
     "nw_sectioned.solve(\"design\")\n",
-    "heat_sect, T_hot_sect, T_cold_sect, _, _ = hx_sectioned.calc_sections()\n",
-    "heat_sect /= 1e6\n",
+    "heat_sect = hx_sectioned.Q_sections.val\n",
+    "T_hot_sect = hx_sectioned.T_hot_sections.val\n",
+    "T_cold_sect = hx_sectioned.T_cold_sections.val\n",
     "\n",
     "hx_mb = MovingBoundaryHeatExchanger(\"heatexchanger\")\n",
     "nw_mb, c1, c2, d1, d2 = make_network(hx_mb)\n",
     "hx_mb.set_attr(dp1=2, dp2=0.5)\n",
     "nw_mb.solve(\"design\")\n",
-    "heat_mb, T_hot_mb, T_cold_mb, _, _ = hx_mb.calc_sections()\n",
-    "heat_mb /= 1e6\n",
+    "heat_mb = hx_mb.Q_sections.val\n",
+    "T_hot_mb = hx_mb.T_hot_sections.val\n",
+    "T_cold_mb = hx_mb.T_cold_sections.val\n",
     "\n",
     "fig, ax = plt.subplots(1, figsize=(10, 6))\n",
     "ax.plot((heat_mb, heat_mb), (list(T_hot_mb), list(T_cold_mb)), color=annotation_color, linestyle=\"-\", linewidth=0.5)\n",
     "ax.plot(heat_sect, T_hot_sect, \"o-\", color=annotation_color, markersize=2, linewidth=0.5, label=\"SectionedHeatExchanger (50 sections)\")\n",
     "ax.plot(heat_mb, T_hot_mb, \"o-\", color=\"red\", label=\"MovingBoundaryHeatExchanger\")\n",
     "ax.plot(heat_mb, T_cold_mb, \"o-\", color=\"blue\")\n",
-    "ax.set_ylabel(\"temperature in K\")\n",
+    "ax.set_ylabel(\"temperature in °C\")\n",
     "ax.set_xlabel(\"heat transferred in MW\")\n",
     "ax.legend()\n",
     "plt.show()"
@@ -594,21 +601,23 @@
     "hx_sectioned = SectionedHeatExchanger(\"heatexchanger\")\n",
     "nw_sectioned, c1, c2, d1, d2 = make_network_preheating(hx_sectioned)\n",
     "nw_sectioned.solve(\"design\")\n",
-    "heat_sect, T_hot_sect, T_cold_sect, _, _ = hx_sectioned.calc_sections()\n",
-    "heat_sect /= 1e6\n",
+    "heat_sect = hx_sectioned.Q_sections.val\n",
+    "T_hot_sect = hx_sectioned.T_hot_sections.val\n",
+    "T_cold_sect = hx_sectioned.T_cold_sections.val\n",
     "\n",
     "hx_mb = MovingBoundaryHeatExchanger(\"heatexchanger\")\n",
     "nw_mb, c1, c2, d1, d2 = make_network_preheating(hx_mb)\n",
     "nw_mb.solve(\"design\")\n",
-    "heat_mb, T_hot_mb, T_cold_mb, _, _ = hx_mb.calc_sections()\n",
-    "heat_mb /= 1e6\n",
+    "heat_mb = hx_mb.Q_sections.val\n",
+    "T_hot_mb = hx_mb.T_hot_sections.val\n",
+    "T_cold_mb = hx_mb.T_cold_sections.val\n",
     "\n",
     "fig, ax = plt.subplots(1, figsize=(10, 6))\n",
     "ax.plot((heat_mb, heat_mb), (list(T_hot_mb), list(T_cold_mb)), color=annotation_color, linestyle=\"-\", linewidth=0.5)\n",
     "ax.plot(heat_sect, T_cold_sect, \"o-\", color=annotation_color, markersize=2, linewidth=0.5, label=\"SectionedHeatExchanger (50 sections)\")\n",
     "ax.plot(heat_mb, T_hot_mb, \"o-\", color=\"red\")\n",
     "ax.plot(heat_mb, T_cold_mb, \"o-\", color=\"blue\", label=\"MovingBoundaryHeatExchanger\")\n",
-    "ax.set_ylabel(\"temperature in K\")\n",
+    "ax.set_ylabel(\"temperature in °C\")\n",
     "ax.set_xlabel(\"heat transferred in MW\")\n",
     "ax.legend()\n",
     "plt.show()"
@@ -662,22 +671,23 @@
     "hx_sectioned = SectionedHeatExchanger(\"heatexchanger\")\n",
     "nw_sectioned, c1, c2, d1, d2 = make_network_supercritical(hx_sectioned)\n",
     "nw_sectioned.solve(\"design\")\n",
-    "heat_sect, T_hot_sect, T_cold_sect, _, _ = hx_sectioned.calc_sections()\n",
-    "heat_sect /= 1e6\n",
+    "heat_sect = hx_sectioned.Q_sections.val\n",
+    "T_hot_sect = hx_sectioned.T_hot_sections.val\n",
+    "T_cold_sect = hx_sectioned.T_cold_sections.val\n",
     "\n",
     "hx_mb = MovingBoundaryHeatExchanger(\"heatexchanger\")\n",
     "nw_mb, c1, c2, d1, d2 = make_network_supercritical(hx_mb)\n",
     "nw_mb.solve(\"design\")\n",
-    "heat_mb, T_hot_mb, T_cold_mb, _, _ = hx_mb.calc_sections()\n",
-    "heat_mb /= 1e6\n",
-    "heatex = SectionedHeatExchanger(\"heatexchanger\")\n",
+    "heat_mb = hx_mb.Q_sections.val\n",
+    "T_hot_mb = hx_mb.T_hot_sections.val\n",
+    "T_cold_mb = hx_mb.T_cold_sections.val\n",
     "\n",
     "fig, ax = plt.subplots(1, figsize=(10, 6))\n",
     "ax.plot((heat_mb, heat_mb), (list(T_hot_mb), list(T_cold_mb)), color=annotation_color, linestyle=\"-\", linewidth=0.5)\n",
     "ax.plot(heat_sect, T_hot_sect, \"o-\", color=annotation_color, markersize=2, linewidth=0.5, label=\"SectionedHeatExchanger (50 sections)\")\n",
     "ax.plot(heat_mb, T_hot_mb, \"o-\", color=\"red\", label=\"MovingBoundaryHeatExchanger\")\n",
     "ax.plot(heat_mb, T_cold_mb, \"o-\", color=\"blue\")\n",
-    "ax.set_ylabel(\"temperature in K\")\n",
+    "ax.set_ylabel(\"temperature in °C\")\n",
     "ax.set_xlabel(\"heat transferred in MW\")\n",
     "ax.legend()\n",
     "plt.show()"