processpi 0.2.0__tar.gz → 0.2.0.1__tar.gz

{processpi-0.2.0 → processpi-0.2.0.1}/PKG-INFO

@@ -1,8 +1,8 @@
 Metadata-Version: 2.4
 Name: processpi
-Version: 0.2.0
+Version: 0.2.0.1
 Summary: Python toolkit for chemical engineering simulations, equipment design, and unit conversions
-Home-page: https://github.com/varma666/ProcessPi
+Home-page: https://processpi.org
 Author: Raviteja Varma Nadimpalli
 Author-email: Raviteja Varma Nadimpalli <processpi.dev@gmail.com>
 License: MIT License
@@ -27,9 +27,10 @@ License: MIT License
         OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
         SOFTWARE.
         
-Project-URL: Homepage, https://github.com/varma666/ProcessPi
-Project-URL: BugTracker, https://github.com/varma666/ProcessPi/issues
-Project-URL: SourceCode, https://github.com/varma666/ProcessPi
+Project-URL: Homepage, https://processpi.org
+Project-URL: Repository, https://github.com/varma666/ProcessPi
+Project-URL: Documentation, https://processpi.org/user-guide/introduction/
+Project-URL: IssueTracker, https://github.com/varma666/ProcessPi/issues
 Requires-Python: >=3.8
 Description-Content-Type: text/markdown
 License-File: LICENSE

processpi-0.2.0.1/processpi/calculations/heat_transfer/__init__.py

@@ -0,0 +1,59 @@
+from .conduction_heat_loss import ConductionHeatLoss
+from .convection_heat_loss import ConvectionHeatLoss
+from .heat_exchanger_area import HeatExchangerArea
+from .biot import BiotNumber
+from .nusselt import NusseltNumber
+from .combined_modes import ConductionConvectionCombined
+from .condensation_dropwise import DropwiseCondensation
+from .condensation_nusselt import NusseltCondensation
+from .condensation import CondensingVapourFilm
+from .conduction_heat_loss import ConductionHeatLoss
+from .convection_heat_loss import ConvectionHeatLoss
+from .crossflow_tube import CrossFlowSingleTube
+from .fourier import FourierNumber
+from .fourierlaw import FourierLaw
+from .heat_exchanger_area import HeatExchangerArea
+from .lmtd import LMTD
+from .newton_cooling import NewtonCooling
+from .ntu import NTUHeatExchanger
+from .nusselt import NusseltNumber
+from .overall_u import OverallHeatTransferCoefficient
+from .peclet import PecletNumber
+from .prandtl import PrandtlNumber
+from .radial_cylinder import RadialHeatFlowCylinder
+from .radiation_blackbody import BlackbodyRadiation
+from .radiation_exchange import RadiationExchange
+from .radiation_greybody import GreybodyRadiation
+from .radiation_viewfactor import RadiationWithViewFactor
+from .resistance import ThermalResistanceSeries, ThermalResistanceParallel
+from .reyleigh import RayleighNumber
+from .rohsenow_boiling import RohsenowBoiling
+from .stefan_boltzmann import StefanBoltzmann
+__all__ = ["ConductionHeatLoss", 
+           "ConvectionHeatLoss", 
+           "HeatExchangerArea", 
+           "BiotNumber", 
+           "NusseltNumber", 
+           "ConductionConvectionCombined", 
+           "DropwiseCondensation", 
+           "NusseltCondensation", 
+           "CondensingVapourFilm", 
+           "FourierNumber",
+           "FourierLaw",
+           "LMTD",
+           "NewtonCooling",
+           "NTUHeatExchanger",
+           "OverallHeatTransferCoefficient",
+           "PecletNumber",
+           "PrandtlNumber",
+           "RadialHeatFlowCylinder",
+           "BlackbodyRadiation",
+           "RadiationExchange",
+           "GreybodyRadiation",
+           "RadiationWithViewFactor",
+           "ThermalResistanceSeries",
+           "ThermalResistanceParallel",
+           "RayleighNumber",
+           "RohsenowBoiling",
+           "StefanBoltzmann"
+    ]

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/calculations/heat_transfer/heat_exchanger_area.py

@@ -14,7 +14,7 @@ class HeatExchangerArea(CalculationBase):
         for key in required:
             if key not in self.inputs:
                 raise ValueError(f"Missing required input: {key}")
-        if self.inputs["log_mean_temp_diff"] <= 0:
+        if self.inputs["log_mean_temp_diff"].value <= 0:
             raise ValueError("Log mean temperature difference must be positive")
 
     def calculate(self):

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/calculations/heat_transfer/lmtd.py

@@ -2,22 +2,21 @@ import math
 from ..base import CalculationBase
 from ...units import *
 
-class LMTDHeatExchanger(CalculationBase):
+class LMTD(CalculationBase):
     """
     Heat Exchanger using LMTD method.
     Q = U * A * ΔTlm
     """
     def validate_inputs(self):
-        required = ["U", "area", "dT1", "dT2"]
+        required = ["dT1", "dT2"]
         for key in required:
             if key not in self.inputs:
                 raise ValueError(f"Missing required input: {key}")
 
     def calculate(self):
-        U = self._get_value(self.inputs["U"], "heat_transfer_coefficient")
-        A = self._get_value(self.inputs["area"], "area")
+
         dT1 = self._get_value(self.inputs["dT1"], "temperature")
         dT2 = self._get_value(self.inputs["dT2"], "temperature")
 
         dTlm = dT1 if dT1 == dT2 else (dT1 - dT2) / math.log(dT1 / dT2)
-        return HeatFlow(U * A * dTlm)
+        return dTlm

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/calculations/heat_transfer/radiation_greybody.py

@@ -1,9 +1,9 @@
 from ..base import CalculationBase
 from ...units import *
 
-class GraybodyRadiation(CalculationBase):
+class GreybodyRadiation(CalculationBase):
     """
-    Graybody radiation with emissivity factor.
+    Greybody radiation with emissivity factor.
 
     **Formula:**
         q = ε * σ * T^4

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/acetone.py

@@ -51,4 +51,4 @@ class Acetone(Component):
     _viscosity_constants = [-14.918, 1023.4, 0.5961, 0, 0] 
     _thermal_conductivity_constants = [0.2878, -0.000427, 0,0,0]
     _vapor_pressure_constants = [69.006, -5599.6, -7.0985, 0.0000062237, 2] 
-    _enthalpy_constants = [4.215E-7, 0.3397, 178.45, 3.639, 0]  # Placeholder for enthalpy constants
+    _enthalpy_constants = [4.215E-7, 0.3397, 0, 0, 0]  # Placeholder for enthalpy constants

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/air.py

@@ -46,9 +46,9 @@ class Air(Component):
     _critical_acentric_factor = 0  # Placeholder for critical acentric factor
     
     # Constants for various property correlations
-    _density_constants = [0.26733, 132.45, 0.27341,59.15]
+    _density_constants = [2.8963, 0.26733, 132.45,0.2734]
     _specific_heat_constants = [-214460,9185.1,-106.12,0.41616,0]
     _viscosity_constants = [-20.077, 285.15, 1.784,-6.2382E-22, 10] 
     _thermal_conductivity_constants = [0.28472, -0.0017393, 0,0,0]
     _vapor_pressure_constants = [21.662, -692.39, -0.392, 0.0047574, 1] 
-    _enthalpy_constants = [0.3822E-7, 59.15,0.6759,132.45, 0]  # Placeholder for enthalpy constants
+    _enthalpy_constants = [0.8474E-7, 0.3822,0,0,0]  # Placeholder for enthalpy constants

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/ammonia.py

@@ -52,3 +52,16 @@ class Ammonia(Component):
     _thermal_conductivity_constants = [1.169, -0.002314, 0,0,0]
     _vapor_pressure_constants = [90.483, -4669.70, -11.607, 1.72E-02, 1] 
     _enthalpy_constants = [3.1523E-7, 0.3914, -0.2289,0.2309, 0]  # Placeholder for enthalpy constants
+
+
+    #@PropertyMethod
+    def specific_heat(self) -> SpecificHeat:
+        if self._specific_heat is not None:
+            return self._specific_heat
+        T = self.temperature.value
+        Tr = T / self._critical_temperature.value
+        t = 1 - Tr
+        cp = (self._sepcific_heat_constants[0] / t ) + (self._sepcific_heat_constants[1]) - ( 2 * self._sepcific_heat_constants[0] * self._sepcific_heat_constants[2] * t) - (self._sepcific_heat_constants[0] * self._sepcific_heat_constants[3] * (t**2)) - ((self._sepcific_heat_constants[2]**2)*(t**3)/3) - ((self._sepcific_heat_constants[2] * self._sepcific_heat_constants[3] * (t**4))/2) - ((self._sepcific_heat_constants[0]**2) * (t**5)/5) 
+        return SpecificHeat(cp, "J/kgK")
+
+        

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/base.py

@@ -121,9 +121,11 @@ class Component(ABC):
         if self._specific_heat is not None:
             return self._specific_heat
         T = self.temperature.value
-        cp = sum(c * (T ** i) for i, c in enumerate(self._specific_heat_constants))
-        cp *= self.molecular_weight
-        return SpecificHeat(cp, "J/kgK")
+        cp_kmol = sum(c * (T ** i) for i, c in enumerate(self._specific_heat_constants))
+        cp_kg = cp_kmol / self.molecular_weight
+        
+        return SpecificHeat(cp_kg, "J/kgK")
+
 
     # ----------------------------------------------------------------------
     # Viscosity
@@ -198,5 +200,5 @@ class Component(ABC):
                        (self._enthalpy_constants[2] * tau) +
                        (self._enthalpy_constants[3] * tau ** 2) +
                        (self._enthalpy_constants[4] * tau ** 3))))
-        dH *= self.molecular_weight
+        #dH *= self.molecular_weight
         return HeatOfVaporization(dH, "J/kg")

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/benzene.py

@@ -1,4 +1,4 @@
-from .base import Component
+from .base import Component, PropertyMethod
 from processpi.units import *
 
 
@@ -28,9 +28,10 @@ class Benzene(Component):
     _viscosity_constants = [7.5117, 294.68, -2.794, 0, 0]
     _thermal_conductivity_constants = [0.23444, -0.00030572, 0, 0, 0]
     _vapor_pressure_constants = [83.107, -6486.2, -9.2194, 0.0000069844, 2]
-    _enthalpy_constants = [4.5346E-7, 0.39053, 278.68, 3.4705, 0]  # Placeholder for enthalpy constants
-
-    def specific_heat(self, temperature: Temperature):
+    _enthalpy_constants = [4.5346E-7, 0.39053, 0, 0, 0]  # Placeholder for enthalpy constants
+    
+    @PropertyMethod
+    def specific_heat(self):
         """
         Calculates the specific heat of Benzene at a given temperature.
 
@@ -45,7 +46,7 @@ class Benzene(Component):
         Returns:
             SpecificHeat: The calculated specific heat in units of J/kgK.
         """
-        T = temperature.value
+        T = self.temperature.value
         # Use a conditional to select the correct set of polynomial constants.
         if T <= 353.24:
             constants = self._specific_heat_constants_1
@@ -63,4 +64,14 @@ class Benzene(Component):
         # the molecular weight (g/mol), since 1 kmol is 1000 mol and 1 kg
         # is 1000 g.
         cp_kg = cp_kmol / self.molecular_weight
+        
         return SpecificHeat(cp_kg, "J/kgK")
+
+    def httype(self):
+        P = self.pressure.to("Pa").value
+        Pvap = self.vapor_pressure().to("Pa").value
+        if P < Pvap:
+            httype = "organicvapour"
+        else:
+            httype = "organicliquid"
+        return httype

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/benzoicacid.py

@@ -16,4 +16,4 @@ class BenzoicAcid(Component):
     _viscosity_constants = [-12.947,2557.9,0,0, 0] 
     _thermal_conductivity_constants = [0.2391, -0.0002325, 0,0,0]
     _vapor_pressure_constants = [88.513, -11829, -8.6826, 2.32E-19, 6] 
-    _enthalpy_constants = [10.19E-7, 0.478, 395.45,7.1277, 0]  # Placeholder for enthalpy constants
+    _enthalpy_constants = [10.19E-7, 0.478,0,0,0]  # Placeholder for enthalpy constants

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/bromine.py

@@ -16,4 +16,4 @@ class Bromine(Component):
     _viscosity_constants = [16.775,-314,-3.9763,0, 0] 
     _thermal_conductivity_constants = [-0.2185, 0.0042143,-0.000017753,3.10E-08,-2.01E-11]
     _vapor_pressure_constants = [108.26, -6592, -14.16,1.60E-02, 1] 
-    _enthalpy_constants = [4E-7, 0.351, 265.85,3.2323, 0]  # Placeholder for enthalpy constants
+    _enthalpy_constants = [4E-7, 0.351, 0,0, 0]  # Placeholder for enthalpy constants

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/carbonmonoxide.py

@@ -15,3 +15,14 @@ class CarbonMonoxide(Component):
     _thermal_conductivity_constants = [0.2855, -0.001784,0,0,0]
     _vapor_pressure_constants = [45.698,-1076.60,-4.8814,7.57E-05, 2] 
     _enthalpy_constants = [0.8585E-7, 0.4921, -0.326,0.2231, 0]  # Placeholder for enthalpy constants
+
+
+    #@PropertyMethod
+    def specific_heat(self) -> SpecificHeat:
+        if self._specific_heat is not None:
+            return self._specific_heat
+        T = self.temperature.value
+        Tr = T / self._critical_temperature.value
+        t = 1 - Tr
+        cp = (self._sepcific_heat_constants[0] / t ) + (self._sepcific_heat_constants[1]) - ( 2 * self._sepcific_heat_constants[0] * self._sepcific_heat_constants[2] * t) - (self._sepcific_heat_constants[0] * self._sepcific_heat_constants[3] * (t**2)) - ((self._sepcific_heat_constants[2]**2)*(t**3)/3) - ((self._sepcific_heat_constants[2] * self._sepcific_heat_constants[3] * (t**4))/2) - ((self._sepcific_heat_constants[0]**2) * (t**5)/5) 
+        return SpecificHeat(cp, "J/kgK")

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/carbontetrachloride.py

@@ -14,4 +14,4 @@ class CarbonTetrachloride(Component):
     _viscosity_constants = [-8.0738,1121.1,-0.4726,0, 0] 
     _thermal_conductivity_constants = [0.1589, -0.0001987,0,0,0]
     _vapor_pressure_constants = [78.441,-6128.10,-8.5766,6.85E-06, 2] 
-    _enthalpy_constants = [4.3252E-7, 0.37688, 250.33,3.4528, 0]  # Placeholder for enthalpy constants
+    _enthalpy_constants = [4.3252E-7, 0.37688, 0,0, 0]  # Placeholder for enthalpy constants

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/chlorobenzene.py

@@ -14,4 +14,4 @@ class ChloroBenzene(Component):
     _viscosity_constants = [0.15772,540.5,-1.6075,1.50E+14,-6.00E+00] 
     _thermal_conductivity_constants = [0.1841, -0.0001917,0,0,0]
     _vapor_pressure_constants = [54.144,-6244.40,-4.5343,4.70E-18, 6] 
-    _enthalpy_constants = [5.148E-7,0.36614, 227.95,4.3707, 0]  # Placeholder for enthalpy constants
+    _enthalpy_constants = [5.148E-7,0.36614, 0,0, 0]  # Placeholder for enthalpy constants

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/chloroform.py

@@ -14,4 +14,4 @@ class Chloroform(Component):
     _viscosity_constants = [-14.109,1049.2,0.5377,0,0] 
     _thermal_conductivity_constants = [0.1778, -0.0002023,0,0,0]
     _vapor_pressure_constants = [146.43,-7792.30,-20.614,2.46E-02, 1] 
-    _enthalpy_constants = [4.186E-7,0.3584, 209.63,3.5047, 0]  # Placeholder for enthalpy constants
+    _enthalpy_constants = [4.186E-7,0.3584, 0,0, 0]  # Placeholder for enthalpy constants

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/chloromethane.py

@@ -14,4 +14,4 @@ class ChloroMethane(Component):
     _viscosity_constants = [-25.132,1381.9,2.0811,-4.50E-27,10] 
     _thermal_conductivity_constants = [0.41067, -0.0008478,0,0,0]
     _vapor_pressure_constants = [64.697,-4048.10,-6.8066,1.04E-05, 2] 
-    _enthalpy_constants = [2.9745E-7,0.353, 175.43,2.452, 0]  # Placeholder for enthalpy constants
+    _enthalpy_constants = [2.9745E-7,0.353, 0,0, 0]  # Placeholder for enthalpy constants

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/cyanogen.py

@@ -14,4 +14,4 @@ class Cyanogen(Component):
     _viscosity_constants = [-12.086,994.23,0,0,0] 
     _thermal_conductivity_constants = [0.4685, -0.00086594,0,0,0]
     _vapor_pressure_constants = [81.565,-4808.90,-9.3748,1.39E-05, 2] 
-    _enthalpy_constants = [3.384E-7,0.3707, 245.25,2.3803, 0]  # Placeholder for enthalpy constants
+    _enthalpy_constants = [3.384E-7,0.3707, 0,0, 0]  # Placeholder for enthalpy constants

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/cyclohexene.py

@@ -14,4 +14,4 @@ class Cyclohexane(Component):
     _viscosity_constants = [-33.763,2497.2,3.2236,0,0] 
     _thermal_conductivity_constants = [0.19813, -0.0002505,0,0,0]
     _vapor_pressure_constants = [51.087,-5226.40,-4.2278,19.76E-18, 6] 
-    _enthalpy_constants = [4.4902E-7,0.39881, 279.69,3.392, 0]  # Placeholder for enthalpy constants
+    _enthalpy_constants = [4.4902E-7,0.39881, 0,0, 0]  # Placeholder for enthalpy constants

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/ethane.py

@@ -15,3 +15,13 @@ class Ethane(Component):
     _thermal_conductivity_constants = [0.35758, -0.0011458,6.19E-07,0,0]
     _vapor_pressure_constants = [51.857,-2598.70,-5.1283,1.49E-05, 2] 
     _enthalpy_constants = [2.1091E-7,0.60646, -0.55492,0.32799, 0]  # Placeholder for enthalpy constants
+
+    #@PropertyMethod
+    def specific_heat(self) -> SpecificHeat:
+        if self._specific_heat is not None:
+            return self._specific_heat
+        T = self.temperature.value
+        Tr = T / self._critical_temperature.value
+        t = 1 - Tr
+        cp = (self._sepcific_heat_constants[0] / t ) + (self._sepcific_heat_constants[1]) - ( 2 * self._sepcific_heat_constants[0] * self._sepcific_heat_constants[2] * t) - (self._sepcific_heat_constants[0] * self._sepcific_heat_constants[3] * (t**2)) - ((self._sepcific_heat_constants[2]**2)*(t**3)/3) - ((self._sepcific_heat_constants[2] * self._sepcific_heat_constants[3] * (t**4))/2) - ((self._sepcific_heat_constants[0]**2) * (t**5)/5) 
+        return SpecificHeat(cp, "J/kgK")

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/ethanol.py

@@ -17,4 +17,4 @@ class Ethanol(Component):
     _viscosity_constants = [7.875, 781.98, -3.0418, 0, 0] 
     _thermal_conductivity_constants = [0.2468, -0.000264, 0,0,0]
     _vapor_pressure_constants = [73.304, -7122.3, -7.1424, 0.0000028853, 2] 
-    _enthalpy_constants = [5.5789E-7, 0.31245, 159.05, 4.9694, 0]  # Placeholder for enthalpy constants
+    _enthalpy_constants = [5.5789E-7, 0.31245, 0, 0, 0]  # Placeholder for enthalpy constants

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/ethylacetate.py

@@ -14,4 +14,4 @@ class EthlylAcetate(Component):
     _viscosity_constants = [14.354,-154.6,-3.7887,0,0] 
     _thermal_conductivity_constants = [0.2501, -0.0003563,0,0,0]
     _vapor_pressure_constants = [66.824,-6227.60,-6.41,1.79E-17, 6] 
-    _enthalpy_constants = [4.933E-7,0.3847, 189.6,4.149, 0]  # Placeholder for enthalpy constants
+    _enthalpy_constants = [4.933E-7,0.3847, 0,0, 0]  # Placeholder for enthalpy constants

processpi-0.2.0.1/processpi/components/ethylene.py

@@ -0,0 +1,17 @@
+from .base import Component
+from processpi.units import *
+class Ethlylene(Component):
+    name = "Ethlylene"
+    formula = "C2H4"
+    molecular_weight = 28.053
+    _critical_temperature = Temperature(282.34, "K")
+    _critical_pressure = Pressure(5.41, "MPa") 
+    _critical_volume = Volume(0.131, "m3")  # Placeholder for critical volume per Kmole
+    _critical_zc = 0.281  # Placeholder for critical compressibility factor
+    _critical_acentric_factor = 0.0862  # Placeholder for critical acentric factor
+    _density_constants = [2.0961,0.27657,282.34,0.29174]
+    _specific_heat_constants = [247390,-4428,40.936,-0.1697,0.00026816]
+    _viscosity_constants = [1.8878,78.865,-2.1554,0,0] 
+    _thermal_conductivity_constants = [0.4194,-0.001591,0.000001306,0,0]
+    _vapor_pressure_constants = [53.963,-2443,-5.5643,0.000019079,2] 
+    _enthalpy_constants = [1.8844E-7,0.36485, 0,0, 0]  # Placeholder for enthalpy constants

processpi-0.2.0.1/processpi/components/fluorine.py

@@ -0,0 +1,17 @@
+from .base import Component
+from processpi.units import *
+class Fluorine(Component):
+    name = "Fluorine"
+    formula = "F2"
+    molecular_weight = 37.997
+    _critical_temperature = Temperature(144.12, "K")
+    _critical_pressure = Pressure(5.172, "MPa") 
+    _critical_volume = Volume(0.066547, "m3")  # Placeholder for critical volume per Kmole
+    _critical_zc = 0.287  # Placeholder for critical compressibility factor
+    _critical_acentric_factor = 0.053  # Placeholder for critical acentric factor
+    _density_constants = [4.2895,0.28587,144.12,0.28776]
+    _specific_heat_constants = [-94585,7529.9,-139.6,1.1301,-0.0033241]
+    _viscosity_constants = [8.18,-75.6,-3.5148,0,0] 
+    _thermal_conductivity_constants = [0.2758,-0.0016297,0,0,0]
+    _vapor_pressure_constants = [42.393,-1103.3,-4.1203,0.000057815,2] 
+    _enthalpy_constants = [0.88757E-7,0.34072, 0,0, 0]  # Placeholder for enthalpy constants

processpi-0.2.0.1/processpi/components/fluorobenzene.py

@@ -0,0 +1,17 @@
+from .base import Component
+from processpi.units import *
+class FluoroBenzene(Component):
+    name = "Fluorobenzene"
+    formula = "C6H5F"
+    molecular_weight = 96.102
+    _critical_temperature = Temperature(560.09, "K")
+    _critical_pressure = Pressure(4.551, "MPa") 
+    _critical_volume = Volume(0.269, "m3")  # Placeholder for critical volume per Kmole
+    _critical_zc = 0.263  # Placeholder for critical compressibility factor
+    _critical_acentric_factor = 0.2472  # Placeholder for critical acentric factor
+    _density_constants = [1.0146,0.27277,560.09,0.28291]
+    _specific_heat_constants = [-991200,11734,-40.669,0.047333,0]
+    _viscosity_constants = [-10.064,1058.7,-0.17162,0,0] 
+    _thermal_conductivity_constants = [0.20962,-0.00028034,0,0,0]
+    _vapor_pressure_constants = [51.915,-5439,-4.2896,8.7527E-18,6] 
+    _enthalpy_constants = [4.582E-7,0.3717, 0,0, 0]  # Placeholder for enthalpy constants

processpi-0.2.0.1/processpi/components/formicacid.py

@@ -0,0 +1,17 @@
+from .base import Component
+from processpi.units import *
+class FormicAcid(Component):
+    name = "Formic Acid"
+    formula = "CH2O2"
+    molecular_weight = 46.026
+    _critical_temperature = Temperature(588, "K")
+    _critical_pressure = Pressure(5.81, "MPa") 
+    _critical_volume = Volume(0.125, "m3")  # Placeholder for critical volume per Kmole
+    _critical_zc = 0.149  # Placeholder for critical compressibility factor
+    _critical_acentric_factor = 0.3173  # Placeholder for critical acentric factor
+    _density_constants = [1.938,0.24225,588,0.24435]
+    _specific_heat_constants = [78060,71.54,0,0,0]
+    _viscosity_constants = [-48.529,3394.7,5.3903,0,0] 
+    _thermal_conductivity_constants = [0.302,-0.000108,0,0,0]
+    _vapor_pressure_constants = [50.323,-5378.2,-4.203,0.0000034697,2] 
+    _enthalpy_constants = [2.3195E-7,1.9091,-5.0003,3.2641, 0]  # Placeholder for enthalpy constants

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/methanol.py

@@ -17,4 +17,4 @@ class Methanol(Component):
     _viscosity_constants = [-25.317, 1789.2, 2.069, 0, 0] 
     _thermal_conductivity_constants = [0.2837, -0.000281, 0,0,0]
     _vapor_pressure_constants = [82.718, -6904.5, -8.8622, 0.0000074664, 2] 
-    _enthalpy_constants = [5.0451e-7, 0.33594, 175.47, 4.3825, 0]  # Placeholder for enthalpy constants
+    _enthalpy_constants = [5.0451e-7, 0.33594, 0, 0, 0]  # Placeholder for enthalpy constants

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/oil.py

@@ -16,3 +16,13 @@ class Oil(Component):
     _thermal_conductivity_constants = [0.18, -0.00025, 0, 0, 0]
     _vapor_pressure_constants = [75, -7000, -7.2, 0.000007, 2]
     _enthalpy_constants = [4.5E-7, 0.34, 200, 3.6, 0]
+
+    def httype(self):
+        P = self.pressure.to("Pa").value
+        Pvap = self.vapor_pressure().to("Pa").value
+        if P < Pvap:
+            httype = "heavyoil"
+        else:
+            httype = "lightoil"
+        return httype
+    

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/organic_liquid.py

@@ -5,6 +5,8 @@ class OrganicLiquid(Component):
     name = "Organic Liquid"
     formula = "R-CH"
     molecular_weight = 100.0
+    httype = "organicliquid"
+    p_h = "liquid"
     _critical_temperature = Temperature(550, "K")
     _critical_pressure = Pressure(4.5, "MPa")
     _critical_volume = Volume(0.25, "m3")
@@ -16,3 +18,12 @@ class OrganicLiquid(Component):
     _thermal_conductivity_constants = [0.2, -0.0003, 0, 0, 0]
     _vapor_pressure_constants = [70, -6000, -7.0, 0.000006, 2]
     _enthalpy_constants = [4.2E-7, 0.3, 180, 3.5, 0]
+    
+    def httype(self):
+        P = self.pressure.to("Pa").value
+        Pvap = self.vapor_pressure().to("Pa").value
+        if P < Pvap:
+            httype = "organicvapour"
+        else:
+            httype = "organicliquid"
+        return httype

{processpi-0.2.0 → processpi-0.2.0.1}/processpi/components/toluene.py

@@ -18,5 +18,14 @@ class Toluene(Component):
     _viscosity_constants = [-226.08, 6805.7,37.542,-0.060853, 1] 
     _thermal_conductivity_constants = [0.20463, -0.00024252, 0,0,0]
     _vapor_pressure_constants = [76.945, -6729.8, -8.179, 0.0000053017, 2] 
-    _enthalpy_constants = [4.9507E-7, 0.37742, 178.18, 4.3246, 0]  # Placeholder for enthalpy constants
-    
+    _enthalpy_constants = [4.9507E-7, 0.37742, 0, 0, 0]  # Placeholder for enthalpy constants
+
+    def httype(self):
+        P = self.pressure.to("Pa").value
+        Pvap = self.vapor_pressure().to("Pa").value
+        if P < Pvap:
+            httype = "organicvapour"
+        else:
+            httype = "organicliquid"
+        return httype
+    

processpi-0.2.0.1/processpi/streams/__init__.py

@@ -0,0 +1,2 @@
+from .material import MaterialStream
+from .energy import EnergyStream

processpi-0.2.0.1/processpi/streams/energy.py

@@ -0,0 +1,67 @@
+from ..units import Power
+from typing import Dict, Optional, Union
+
+
+class EnergyStream:
+    """
+    Represents an energy stream (heat or work transfer).
+    Can bind to equipment and log Q/W duties automatically.
+
+    Notes
+    -----
+    * Positive duty = energy input (to equipment)
+    * Negative duty = energy output (from equipment)
+    * Duties are stored with units (`Power`) to ensure consistency.
+    """
+
+    def __init__(self, name: str):
+        self.name = name
+        self.log: list[Dict[str, Union[str, Power]]] = []  # {equipment, tag, duty}
+
+    def bind_equipment(self, equipment) -> None:
+        """Optionally keep reference to equipment if needed later."""
+        self.log.append({
+            "equipment": equipment.name,
+            "tag": "bind",
+            "duty": Power(0.0, "W")
+        })
+
+    def record(self, duty: Union[float, Power], tag: str, equipment: str) -> None:
+        """
+        Log an energy duty.
+
+        Parameters
+        ----------
+        duty : float | Power
+            Energy duty (W). If float is given, assumed in watts.
+        tag : str
+            Label such as 'Q_in', 'Q_out', 'W_in', 'W_out'.
+        equipment : str
+            Name of equipment associated with this duty.
+        """
+        if isinstance(duty, (int, float)):
+            duty = Power(duty, "W")
+
+        self.log.append({
+            "equipment": equipment,
+            "tag": tag,
+            "duty": duty
+        })
+
+    def total_duty(self, tag: Optional[str] = None) -> Power:
+        """
+        Return cumulative duty as a Power object.
+        Optionally filter by tag.
+        """
+        duties = (
+            [event["duty"] for event in self.log if tag is None or event["tag"] == tag]
+        )
+        total_value = sum(d.value for d in duties)
+        return Power(total_value, "W")
+
+    def __repr__(self) -> str:
+        return (
+            f"EnergyStream(name={self.name}, "
+            f"total={self.total_duty()}, "
+            f"entries={len(self.log)})"
+        )