thermofluids 2.0.0__py3-none-any.whl

thermofluids/ProcessEquipment/Turbine.py

@@ -0,0 +1,53 @@
+from thermofluids import Statepoint, Quantity
+from thermofluids.Base import Quantity
+import types
+from CoolProp.CoolProp import PropsSI
+
+class Turbine():
+    def __init__(self, inlet=None, eta=None, Pout=None, output_name=None, work=None):
+        """
+        Models a turbine
+        :param inlet: [Statepoint] turbine inlet stream
+        :param eta: [Quantity] or [float] icentropic efficiency of turbo equipment
+        :param Pout: [Quantity] or [float] pressure at outlet of turbine
+        :return: [Statepoint] turbine outlet stream
+        """
+        self.inlet = inlet
+        self.eta = eta
+        self.Pout = Pout
+        self.output_name = output_name
+        self.work = work
+
+    def set_inlet(self, inlet):
+        self.inlet = inlet
+
+    def set_eta(self, eta):
+        self.eta = eta
+
+    def set_Pout(self, Pout):
+        self.Pout = Pout
+
+    def set_output(self, output_name):
+        self.output_name = output_name
+
+    def get_outlet(self):
+        assert isinstance(self.inlet, Statepoint), '\'inlet\' must be set to a Statepoint'
+        # assert isinstance(self.eta, Quantity or types.FloatType), 'eta must be set to a Quantity or Float'
+
+        self.outlet = Statepoint(name=self.output_name, fluid=self.inlet.fluid)
+        self.outlet.set_source(self.inlet)
+        self.outlet.set_m(self.inlet.m)
+        self.outlet.set_p(self.Pout)
+        try:
+            self.outlet.set_c(self.inlet.c)
+        except:
+            pass
+
+        outlet_s = Statepoint(name='outlet_s', fluid=self.inlet.fluid)
+        outlet_s.set_s(self.inlet.s)
+        outlet_s.set_p(self.Pout)
+
+        self.outlet.set_h(self.inlet.h - self.eta*(self.inlet.h - outlet_s.h))
+        self.work = self.outlet.m*(self.inlet.h - self.outlet.h)
+
+        return self.outlet

thermofluids/ProcessEquipment/Valve.py

@@ -0,0 +1,35 @@
+"""
+Manipulations of definition of flow coefficient, Cv=Q*sqrt(SG/dP), to calculate valve performance.
+https://en.wikipedia.org/wiki/Flow_coefficient
+"""
+from thermofluids import Quantity
+
+
+def get_valve_flow_rate(statepoint, dP, C_v, units='kg*s^-1', debug=False):
+    flow_gpm = Quantity((C_v * (dP.to('psi') / (statepoint.rho / Quantity(1000., 'kg*m^-3')))**0.5).magnitude, 'gal*min^-1')
+    if debug:
+        print(('\tdP [psi]: %f' % dP.to('psi')))
+        print(('\tDensity: %f' % (statepoint.rho / Quantity(1000., 'kg*m^-3'))))
+        print(('\tFlow: %f' % flow_gpm))
+    return (flow_gpm * Quantity(0.00379, 'm^3*gal^-1') / Quantity(60., 's*min^-1') * statepoint.rho).to(units)
+
+
+def get_valve_dP(statepoint, C_v, units='bar'):
+    dP = (statepoint.m / statepoint.rho).to('gal/min')**2. * statepoint.rho / Quantity(1000., 'kg*m^-3') / C_v ** 2.
+    return Quantity(dP.magnitude, 'psi').to(units)
+
+
+def get_valve_C_v(statepoint, dP):
+    return ((statepoint.m / statepoint.rho).to('gal/min') * (statepoint.rho / Quantity(1000., 'kg*m^-3') / dP.to('psi'))**0.5).magnitude
+
+if __name__ == '__main__':
+    from thermofluids.Base import Statepoint
+
+    C_v = 430.
+    sp = Statepoint(fluid='ammonia', p=15., t=80., m=2.)
+    dP = get_valve_dP(statepoint=sp, C_v=C_v)
+    print(dP)
+    mdot = get_valve_flow_rate(statepoint=sp, dP=dP, C_v=C_v)
+    print(mdot)
+    C_v_check = get_valve_C_v(statepoint=sp, dP=dP)
+    print(C_v_check)

thermofluids/ProcessEquipment/__init__.py

@@ -0,0 +1,3 @@
+from thermofluids.ProcessEquipment.Compressor import Compressor
+from thermofluids.ProcessEquipment.Turbine import Turbine
+from thermofluids.ProcessEquipment.Valve import get_valve_flow_rate, get_valve_C_v, get_valve_dP

thermofluids/Utilities/FlowDiagram.py

@@ -0,0 +1,111 @@
+from thermofluids import Statepoint, Quantity
+
+def WriteFlow(template, SP, filename, SC):
+    text = open(template).read()
+    for i in range(0,len(SP),1):
+        try:
+            try:
+                # text = text.replace('T%s<'%(i), '%.3G'  %float(SP[i].t) + u'\u2103' + '<')
+                text = text.replace(f'T{i}<', f'{SP[i].t.to("degC").magnitude: .3g}\u2103<')
+            except AttributeError:
+                pass
+            try:
+                # text = text.replace('M%s<'%(i), '%.4G kg/s<' %float(SP[i].m))
+                text = text.replace(f'M{i}<', f'{SP[i].m.to("kg/s").magnitude: .4g} kg/s<')
+            except AttributeError:
+                pass
+            try:
+                # text = text.replace('C%s<'%(i), '%.2G wtCon<' %float(SP[i].c))
+                text = text.replace(f'C{i}<', f'{SP[i].c: .2g} wtCon<')
+            except AttributeError:
+                pass
+            try:
+                # text = text.replace('P%s<'%(i), '%.4G bar<' %float(SP[i].p))
+                text = text.replace(f'P{i}<', f'{SP[i].p.to("bar").magnitude: .4g} bar<')
+            except AttributeError:
+                pass
+            try:
+                # text = text.replace('H%s<'%(i), '%.2G kJ/kg<' %float(SP[i].h))
+                text = text.replace(f'H{i}<', f'{SP[i].h.to("kJ/kg").magnitude: .2g} kJ/kg<')
+            except AttributeError:
+                pass
+            try:
+                # text = text.replace('S%s<'%(i), '%.2G kJ/kg-C<' %float(SP[i].s))
+                text = text.replace(f'S{i}<', f'{SP[i].s.to("kJ/(kg*delta_degC)").magnitude: .2g} kJ/kg-C<')
+            except AttributeError:
+                pass
+            try:
+                if 'humid' not in SP[i].fluid: #humid air has no quality term
+                    # text = text.replace('X%s<'%(i), '%.2G<' %float(SP[i].x))
+                    text = text.replace(f'X{i}<', f'{SP[i].x.to("").magnitude: .2g}<')
+            except AttributeError:
+                pass
+            except ValueError:
+                pass
+        except ValueError as e:
+            raise ValueError(f'Error diagramming statepoint {i}', e)
+
+    keys = SC.keys()
+    for key in keys:
+        try:
+            text = text.replace(str(key), str(SC[key]))
+        except AttributeError:
+            pass
+
+    flow = open(filename, 'w', encoding='utf-8')
+    flow.write(text)
+    flow.close()
+
+def WriteFlowEnglish(template, SP, filename, SC):
+    text = open(template).read()
+    for i in range(0,len(SP),1):
+        try:
+            # text = text.replace('T%s<'%(i), '%.3G' %float(SP[i].t.toEnglish()) + u'\u2109' + '<')
+            text = text.replace(f'T{i}<', f'{SP[i].t.to("degF")}<')
+        except AttributeError:
+            pass
+        try:
+            # text = text.replace('M%s<'%(i), '%.3G lb/s<' %float(SP[i].m.toEnglish()))
+            text = text.replace(f'M{i}<', f'{SP[i].m.to("lb/s")}<')
+        except AttributeError:
+            pass
+        try:
+            # text = text.replace('C%s<'%(i), '%.2G <' %float(SP[i].c))
+            text = text.replace(f'C{i}<', f'{SP[i].c: .2g} wtCon<')
+        except AttributeError:
+            pass
+        try:
+            # text = text.replace('P%s<'%(i), '%.3G psi<' %float(SP[i].p.toEnglish()))
+            text = text.replace(f'P{i}<', f'{SP[i].p.to("psi")}<')
+        except AttributeError:
+            pass
+        try:
+            # text = text.replace('H%s<'%(i), '%.2G Btu/lb<' %float(SP[i].h.toEnglish()))
+            text = text.replace(f'H{i}<', f'{SP[i].h.to("BTU/lb")}<')
+        except AttributeError:
+            pass
+        try:
+            # text = text.replace('S%s<'%(i), '%.2G Btu/lb-F<' %float(SP[i].s.toEnglish()))
+            text = text.replace(f'S{i}<', f'{SP[i].s.to("BTU/(lb*delta_degF)")}<')
+        except AttributeError:
+            pass
+        try:
+            if 'humid' not in SP[i].fluid:  # humid air has no quality term
+                # text = text.replace('X%s<'%(i), '%.2G<' %float(SP[i].x))
+                text = text.replace(f'X{i}<', f'{SP[i].x.to("").magnitude: .2g}<')
+
+        except AttributeError:
+            pass
+        except ValueError:
+            pass
+
+    keys = SC.keys()
+    for key in keys:
+        try:
+            text = text.replace(str(key), str(SC[key]))
+        except AttributeError:
+            pass
+
+    flow = open(filename, 'w')
+    flow.write(text)
+    flow.close()

thermofluids/Utilities/__init__.py

@@ -0,0 +1 @@
+

thermofluids/__init__.py

@@ -0,0 +1,3 @@
+from .Base.Statepoint import Statepoint
+from .Base.Quantity import Quantity
+from .Base.Quantity import Quantity, pint_unit_registry

thermofluids-2.0.0.dist-info/METADATA

@@ -0,0 +1,29 @@
+Metadata-Version: 2.5
+Name: thermofluids
+Version: 2.0.0
+Summary: Physically-consistent thermal and fluids modeling tools
+Project-URL: Homepage, https://gitlab.com/ericksonla/thermofluids
+Author: Luke Erickson, Russell Goldfarbmuren
+Author-email: russell@rebound-tech.com
+License-Expression: MIT
+Classifier: Operating System :: OS Independent
+Classifier: Programming Language :: Python :: 3
+Requires-Python: >3.11
+Requires-Dist: coolprop
+Requires-Dist: pint
+Description-Content-Type: text/markdown
+
+# Thermofluids
+
+Thermofluids makes modeling thermofluid systems easier. The essential element is the `statepoint`. `statepoint` objects work like physical statepoints: they have two independent variables, many dependent variables, and sometimes extensive variables. You can use them like this:
+````python
+sp = Statepoint(fluid='water',  # define a fluid
+                t=25.0,         # define one independent variable
+                p=1.0,          # define a second independent variable
+                m=0.1)          # define extensive varible
+
+sp.rho, sp.h, sp.visc
+# (<Quantity(997.047039, 'kilogram / meter ** 3')>, <Quantity(104.918893, 'kilojoule / kilogram')>, <Quantity(0.000890022674, 'kilogram / meter / second')>)
+````
+
+A set of functions are made available to analyze statepoint systems in physically meaningful ways.

thermofluids-2.0.0.dist-info/RECORD

@@ -0,0 +1,32 @@
+thermofluids/__init__.py,sha256=8CQPY_pPZ_2YR5YtHXNnM6SZKJjnKz-MaiHVlMfubQk,131
+thermofluids/Base/Quantity.py,sha256=BJzbv40aArDHvCToMopktiejjnYJtJSp_6Vi2Y3h8fc,438
+thermofluids/Base/Statepoint.py,sha256=TSX4y01u4A7M_icWOzil4wQnr1S1Sjp4ZSUgEl-sDis,51734
+thermofluids/Base/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+thermofluids/Base/thermofluids_registry.txt,sha256=ngeUgAEP10I9LZU5OnDgjp7SRa2SZsWiEs7Nt9r7ADo,127
+thermofluids/Fluids/Drag.py,sha256=BW8e0GxGP9MyRLNPBfFxoS23NACWoeUU8KE60qztwAo,6798
+thermofluids/Fluids/Moody.py,sha256=bjes0IjehKE-oEi7oGe1IuUz00clPosXzUgbhuhGxZ4,1761
+thermofluids/Fluids/PressureDrop.py,sha256=k1Rgz_WvUo5RKqmT-x87B8aVkJC8CGWqfoFrQOFzIk4,13059
+thermofluids/Fluids/Re.py,sha256=6mbmKVMlOl14k7D1e0fe55KRoa1KTIGBU9IN6ADWfNU,3763
+thermofluids/Fluids/__init__.py,sha256=k0gTYO3okjUznTK9NmtXU4oQUxbMfwYddMZSqADbX9Q,361
+thermofluids/HeatTransfer/Gr.py,sha256=dSztVuazHJ5MlEjfey6lgNoW9QRxcmP1D4dudlEiL-o,365
+thermofluids/HeatTransfer/NTU.py,sha256=OSTXam2T9kwnuyVEHCjZQx-6g8AdEwc4zBPIOw9zL8E,1102
+thermofluids/HeatTransfer/Nusselt.py,sha256=bICchSw-0nWGGzffomi2lljKIJz_KVQBj3CU9aw3Kac,2486
+thermofluids/HeatTransfer/Ra.py,sha256=wlNaVwagmNwO0nHZ6mRzwnNay3hggbHsH9-2W6T8-58,819
+thermofluids/HeatTransfer/TwoPhase.py,sha256=yl5Idq8PvVv5bIw437b-IBKostEuTJlnsWpzfLJNZUs,9537
+thermofluids/HeatTransfer/__init__.py,sha256=BqxMzlHV0p4OiQCu3d-DveVnAgX2DHV9uvO6bVDVzZQ,315
+thermofluids/HeatTransfer/lmtd.py,sha256=afqicephB-Z0QjVyN-Se0V4tzf83m9Ve6SMJ3L5m6Gw,959
+thermofluids/ProcessEquipment/Column.py,sha256=G223CtU4FX_BL-eA5LJosbsaARdLpGw0_-EXOuJ-_Ok,50879
+thermofluids/ProcessEquipment/Compressor.py,sha256=-4DCwryWvB4WRyln2ui3scJz7W9UkTi1AZfpCPh09Lw,6841
+thermofluids/ProcessEquipment/Ejector.py,sha256=05kzO7cg2zm-Nd5qlyb0638qiQQs6WO_ofRWijUrtoA,5716
+thermofluids/ProcessEquipment/HeatExchanger.py,sha256=z2KEcW1k2PJRyjtXZ5JDwIQ5Wy87-y-aFYgldxEBIHE,5968
+thermofluids/ProcessEquipment/IsoCompressor.py,sha256=d2AUENbt9zximAhDfNGacPrxjeX_WzCL3g5iG28oEkY,4880
+thermofluids/ProcessEquipment/Pipe.py,sha256=6EkNssz929K4z0tLuBA1iuCdlyPH52EVq9-vHkJF2Kg,458
+thermofluids/ProcessEquipment/Pump.py,sha256=v5ui1Z7vwO-YAEpt_VDcya3KRxT6hpHzimO2fYUzEwI,6172
+thermofluids/ProcessEquipment/Turbine.py,sha256=GDYRK_VsAnBKfpbNt6YRFSVAxK71kkLoLDxIvlOnFSo,1772
+thermofluids/ProcessEquipment/Valve.py,sha256=SYfo5N1HJTSdyFdXPdEjC2Bhv2r9duSAy8d2Lqdky5E,1436
+thermofluids/ProcessEquipment/__init__.py,sha256=GQ6YyM-WyJJhq-9W-mtFdw2-bxSPcARYh0gTmEIywPg,219
+thermofluids/Utilities/FlowDiagram.py,sha256=Tx3sxX4dZITKHnKqQO2tCunibWzm42lzH7nf-7IT_jw,4377
+thermofluids/Utilities/__init__.py,sha256=AbpHGcgLb-kRsJGnwFEktk7uzpZOCcBY74-YBdrKVGs,1
+thermofluids-2.0.0.dist-info/METADATA,sha256=_fH91j23jgQ983MtuC6miLTJxF6acqafJYk4JsmtVGw,1332
+thermofluids-2.0.0.dist-info/WHEEL,sha256=W3fkpkm7-wf9vBI5Z-7s0eWkeM-spu78I8Neb98DeEg,87
+thermofluids-2.0.0.dist-info/RECORD,,

thermofluids-2.0.0.dist-info/WHEEL

@@ -0,0 +1,4 @@
+Wheel-Version: 1.0
+Generator: hatchling 1.32.4
+Root-Is-Purelib: true
+Tag: py3-none-any