processpi 0.0.0__py3-none-any.whl

processpi/components/chlorine.py

@@ -0,0 +1,17 @@
+from .base import Component
+from processpi.units import *
+class Chlorine(Component):
+    name = "Chlorine"
+    formula = "Cl2"
+    molecular_weight = 70.906
+    _critical_temperature = Temperature(417.15, "K")
+    _critical_pressure = Pressure(7.71, "MPa") 
+    _critical_volume = Volume(0.124, "m3")  # Placeholder for critical volume per Kmole
+    _critical_zc = 0.276  # Placeholder for critical compressibility factor
+    _critical_acentric_factor = 0.0688  # Placeholder for critical acentric factor
+    _density_constants = [2.23,0.27645,417.15,0.2926]
+    _specific_heat_constants = [63936,46.35,-0.1623,0,0]
+    _viscosity_constants = [-9.5412,456.62,0,0,0] 
+    _thermal_conductivity_constants = [0.2246, -0.000064,-0.000000788,0,0]
+    _vapor_pressure_constants = [71.334,-3855,-8.5171,1.24E-02, 1] 
+    _enthalpy_constants = [3.068E-7, 0.8458, -0.9001,0.453, 0]  # Placeholder for enthalpy constants

processpi/components/chlorobenzene.py

@@ -0,0 +1,17 @@
+from .base import Component
+from processpi.units import *
+class ChloroBenzene(Component):
+    name = "Chloro Benzene"
+    formula = "C6​H5​Cl"
+    molecular_weight = 112.557
+    _critical_temperature = Temperature(632.35, "K")
+    _critical_pressure = Pressure(4.519, "MPa") 
+    _critical_volume = Volume(0.308, "m3")  # Placeholder for critical volume per Kmole
+    _critical_zc = 0.265  # Placeholder for critical compressibility factor
+    _critical_acentric_factor = 0.2499  # Placeholder for critical acentric factor
+    _density_constants = [0.8711,0.26805,632.35,0.2799]
+    _specific_heat_constants = [-1307500,15338,-53.974,0.063483,0]
+    _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

processpi/components/chloroform.py

@@ -0,0 +1,17 @@
+from .base import Component
+from processpi.units import *
+class Chloroform(Component):
+    name = "Chloroform"
+    formula = "CHCl3​"
+    molecular_weight = 119.378
+    _critical_temperature = Temperature(536.4, "K")
+    _critical_pressure = Pressure(5.472, "MPa") 
+    _critical_volume = Volume(0.239, "m3")  # Placeholder for critical volume per Kmole
+    _critical_zc = 0.293  # Placeholder for critical compressibility factor
+    _critical_acentric_factor = 0.2219  # Placeholder for critical acentric factor
+    _density_constants = [1.0841,0.2581,536.4,0.2741]
+    _specific_heat_constants = [124850,-166.34,0.43209,0,0]
+    _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

processpi/components/chloromethane.py

@@ -0,0 +1,17 @@
+from .base import Component
+from processpi.units import *
+class ChloroMethane(Component):
+    name = "Chloro Methane"
+    formula = "CH3Cl​"
+    molecular_weight = 50.488
+    _critical_temperature = Temperature(416.25, "K")
+    _critical_pressure = Pressure(6.68, "MPa") 
+    _critical_volume = Volume(0.143, "m3")  # Placeholder for critical volume per Kmole
+    _critical_zc = 0.276  # Placeholder for critical compressibility factor
+    _critical_acentric_factor = 0.1531  # Placeholder for critical acentric factor
+    _density_constants = [1.817,0.25877,416.25,0.2833]
+    _specific_heat_constants = [96910,-207.9,0.37456,0.000488,0]
+    _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

processpi/components/cyanogen.py

@@ -0,0 +1,17 @@
+from .base import Component
+from processpi.units import *
+class Cyanogen(Component):
+    name = "Cyanogen"
+    formula = "C2N2​"
+    molecular_weight = 52.035
+    _critical_temperature = Temperature(400.15, "K")
+    _critical_pressure = Pressure(5.98, "MPa") 
+    _critical_volume = Volume(0.195, "m3")  # Placeholder for critical volume per Kmole
+    _critical_zc = 0.351  # Placeholder for critical compressibility factor
+    _critical_acentric_factor = 0.279  # Placeholder for critical acentric factor
+    _density_constants = [1.0743,0.20948,400.15,0.20724]
+    _specific_heat_constants = [65516,-144.7,0.063229,0,0]
+    _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

processpi/components/cyclohexene.py

@@ -0,0 +1,17 @@
+from .base import Component
+from processpi.units import *
+class Cyclohexane(Component):
+    name = "Cyclohexane"
+    formula = "C6H12​"
+    molecular_weight = 84.159
+    _critical_temperature = Temperature(553.8, "K")
+    _critical_pressure = Pressure(4.08, "MPa") 
+    _critical_volume = Volume(0.308, "m3")  # Placeholder for critical volume per Kmole
+    _critical_zc = 0.273  # Placeholder for critical compressibility factor
+    _critical_acentric_factor = 0.2081  # Placeholder for critical acentric factor
+    _density_constants = [0.88998,0.27376,553.8,0.28571]
+    _specific_heat_constants = [-220600,3118.30,-9.4216,0.010687,0]
+    _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

processpi/components/ethane.py

@@ -0,0 +1,17 @@
+from .base import Component
+from processpi.units import *
+class Ethane(Component):
+    name = "Ethane"
+    formula = "C2H6​"
+    molecular_weight = 30.069
+    _critical_temperature = Temperature(305.32, "K")
+    _critical_pressure = Pressure(4.872, "MPa") 
+    _critical_volume = Volume(0.1455, "m3")  # Placeholder for critical volume per Kmole
+    _critical_zc = 0.279  # Placeholder for critical compressibility factor
+    _critical_acentric_factor = 0.0995  # Placeholder for critical acentric factor
+    _density_constants = [1.9122,0.27937,305.32,0.29187]
+    _specific_heat_constants = [44.009,89718,918.77,-1886,0]
+    _viscosity_constants = [-7.0046,276.38,-0.6087,-3.11E-18,7] 
+    _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

processpi/components/ethanol.py

@@ -0,0 +1,20 @@
+
+from .base import Component
+from processpi.units import *
+
+
+class Ethanol(Component):
+    name = "Ethanol"
+    formula = "C2H5OH"
+    molecular_weight = 46.068
+    _critical_temperature = Temperature(514.00, "K")
+    _critical_pressure = Pressure(6.137, "MPa") 
+    _critical_volume = Volume(0.168, "m3")  # Placeholder for critical volume per Kmole
+    _critical_zc = 0.241  # Placeholder for critical compressibility factor
+    _critical_acentric_factor = 0.6436  # Placeholder for critical acentric factor
+    _density_constants = [1.6288, 0.27469, 514, 0.23178]
+    _specific_heat_constants = [102640,-139.63,-0.030341,0.0020386,0]
+    _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

processpi/components/ethylacetate.py

@@ -0,0 +1,17 @@
+from .base import Component
+from processpi.units import *
+class EthlylAcetate(Component):
+    name = "Ethlyl Acetate"
+    formula = "C4​H8​O2​​"
+    molecular_weight = 88.105
+    _critical_temperature = Temperature(523.3, "K")
+    _critical_pressure = Pressure(3.88, "MPa") 
+    _critical_volume = Volume(0.286, "m3")  # Placeholder for critical volume per Kmole
+    _critical_zc = 0.255  # Placeholder for critical compressibility factor
+    _critical_acentric_factor = 0.3664  # Placeholder for critical acentric factor
+    _density_constants = [0.8996,0.25856,523.3,0.278]
+    _specific_heat_constants = [226230,-624.8,1.472,0,0]
+    _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

processpi/components/methanol.py

@@ -0,0 +1,20 @@
+
+from .base import Component
+from processpi.units import *
+
+
+class Methanol(Component):
+    name = "Methanol"
+    formula = "CH3OH"
+    molecular_weight = 32.042
+    _critical_temperature = Temperature(512.6, "K")
+    _critical_pressure = Pressure(8.1, "MPa")
+    _critical_volume = Volume(0.117, "m3")  # Placeholder for critical volume per Kmole
+    _critical_zc = 0.222  # Placeholder for critical compressibility factor
+    _critical_acentric_factor = 0.5658  # Placeholder for critical acentric factor
+    _density_constants = [2.3267, 0.27073, 512.5, 0.24713]
+    _specific_heat_constants = [105800,-362.23,0.9379,0,0]
+    _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

processpi/components/toluene.py

@@ -0,0 +1,22 @@
+
+from .base import Component
+from processpi.units import *
+
+
+class Toluene(Component):
+    name = "Toluene"
+    formula = "C7H8"
+    molecular_weight = 92.138
+    _critical_temperature = Temperature(591.8, "K")
+    _critical_pressure = Pressure(4.1, "MPa")
+    _critical_volume = Volume(0.316, "m3")  # Placeholder for critical volume per Kmole
+    _critical_zc = 0.264  # Placeholder for critical compressibility factor
+    _critical_acentric_factor = 0.264  # Placeholder for critical acentric factor
+    _density_constants = [0.8792, 0.27136, 591.75, 0.29241]
+    _specific_heat_constants = [140140,-152.3,0.695,0,0]
+   
+    _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
+    

processpi/components/water.py

@@ -0,0 +1,24 @@
+from .base import Component
+from processpi.units import *
+from math import *
+
+class Water(Component):
+    name = "Water"
+    formula = "H2O"
+    molecular_weight = 18.015
+    _critical_temperature = Temperature(647.096, "K")
+    _critical_pressure = Pressure(22.06, "MPa")
+    _density_constants = [-13.851, 0.654, -0.00191, 0.000002]
+    _specific_heat_constants = [276370,-2090.10,8.125,-0.014116,9.37e-6]
+    _viscosity_constants = [-52.843, 3703.6, 5.866, -5.879e-29, 10] 
+    _thermal_conductivity_constants = [-0.432, 0.0057255, -0.000008078,1.861E-09,0]
+    _vapor_pressure_constants = [73.649, -7258.2, -7.3037, 4.1653E-06, 2] 
+    _enthalpy_constants = [5.2053E-7, 0.3199, -0.212, 0.25795, 0]
+
+    def density(self, temperature: Temperature = None) -> Density:
+        temperature = temperature or self.temperature
+        Tr = temperature.value / self._critical_temperature.value
+        tou = 1 - Tr
+        den = 17.863 + (58.606 * (tou**0.35)) -  (95.396 * (tou**(2/3))) + (213.89 * tou) - (141.26 * (tou**(4/3)))
+        den = den * self.molecular_weight
+        return Density(den, "kg/m3")

processpi/pipelines/base.py

@@ -0,0 +1,76 @@
+# processpi/pipelines/base.py
+
+from abc import ABC, abstractmethod
+from typing import Any, Dict
+
+
+class PipelineBase(ABC):
+    """
+    Abstract base class for all pipeline-related calculations in ProcessPI.
+    
+    Notes:
+    - Keeps parameter handling and validation generic.
+    - Specific pipe attributes (diameter, schedule, roughness, material)
+      are defined in `pipes.py` and should be used there.
+    """
+
+    def __init__(self, **kwargs):
+        """
+        Initialize a pipeline calculation with input parameters.
+
+        Args:
+            **kwargs: Arbitrary keyword arguments for inputs.
+        """
+        self.params = kwargs
+
+    @abstractmethod
+    def calculate(self) -> Dict[str, Any]:
+        """
+        Perform the pipeline calculation.
+        
+        Must be implemented by subclasses.
+
+        Returns:
+            dict: Results of the calculation in key-value format.
+        """
+        pass
+
+    def get_param(self, key: str, default: Any = None) -> Any:
+        """
+        Safely fetch an input parameter.
+
+        Args:
+            key (str): Parameter name.
+            default (Any): Default value if key not found.
+
+        Returns:
+            Any: The value of the parameter or the default.
+        """
+        return self.params.get(key, default)
+
+    def validate_required(self, required_keys: list):
+        """
+        Ensure that required input parameters are present.
+
+        Args:
+            required_keys (list): List of required parameter names.
+
+        Raises:
+            ValueError: If any required parameter is missing.
+        """
+        missing = [k for k in required_keys if k not in self.params]
+        if missing:
+            raise ValueError(f"Missing required parameters: {', '.join(missing)}")
+
+    def _getattr(self, attr_name: str, default: Any = None) -> Any:
+        """
+        Safely fetch an attribute from the pipeline instance.
+
+        Args:
+            attr_name (str): Name of the attribute.
+            default (Any): Default value if attribute not found.
+
+        Returns:
+            Any: The value of the attribute or the default.
+        """
+        return getattr(self, attr_name, default)

processpi/pipelines/design_rules.py

@@ -0,0 +1,91 @@
+# processpi/pipelines/design_rules.py
+
+"""
+Pipeline Design Rules Module
+
+This module defines pipeline design rules inspired by industry practices 
+(similar to Aspen). It does not contain calculation logic directly, 
+but instead integrates functions from processpi.calculations to validate 
+pipeline designs against engineering standards.
+
+Each function acts as a design "rule" or "check".
+"""
+
+from processpi.calculations.fluids.reynolds_number import ReynoldsNumber
+from processpi.calculations.fluids.friction_factor import ColebrookWhite
+from processpi.calculations.fluids.flow_type import FlowType
+from processpi.calculations.fluids.optimum_pipe_diameter import OptimumPipeDiameter
+from processpi.calculations.pressure.drop import PressureDrop  # assumed module
+from processpi.calculations.velocity import Velocity           # assumed module
+
+
+def check_flow_regime(velocity: float, diameter: float, density: float, viscosity: float) -> str:
+    """
+    Rule: Determine if flow is laminar, transition, or turbulent.
+    """
+    Re = ReynoldsNumber(velocity, diameter, density, viscosity).calculate()
+    return FlowType(Re).calculate()
+
+
+def check_pressure_drop(flowrate: float, diameter: float, length: float, density: float, viscosity: float, roughness: float) -> float:
+    """
+    Rule: Calculate pressure drop across the pipeline.
+    """
+    return PressureDrop(flowrate, diameter, length, density, viscosity, roughness).calculate()
+
+
+def check_friction_factor(Re: float, diameter: float, roughness: float) -> float:
+    """
+    Rule: Calculate friction factor using Colebrook-White.
+    """
+    return ColebrookWhite(Re, diameter, roughness).calculate()
+
+
+def recommend_optimum_diameter(flowrate: float, density: float) -> float:
+    """
+    Rule: Recommend optimum diameter based on economic sizing correlation.
+    """
+    return OptimumPipeDiameter(flowrate, density).calculate()
+
+
+def check_velocity_limits(velocity: float, fluid_type: str = "liquid") -> bool:
+    """
+    Rule: Check if velocity is within recommended design limits.
+
+    Typical guidelines:
+    - Liquids: 1–3 m/s (avoid erosion/cavitation)
+    - Gases/steam: 10–20 m/s
+    """
+    if fluid_type.lower() == "liquid":
+        return 1.0 <= velocity <= 3.0
+    elif fluid_type.lower() in ["gas", "steam"]:
+        return 10.0 <= velocity <= 20.0
+    return False
+
+
+def check_pressure_drop_limit(pressure_drop: float, limit: float = 10.0) -> bool:
+    """
+    Rule: Ensure pressure drop per 100 m of pipe length is within design limits (default 10 kPa/100 m).
+    """
+    return pressure_drop <= limit
+
+
+def check_design_compliance(flowrate: float, diameter: float, length: float, density: float, viscosity: float, roughness: float, fluid_type: str = "liquid") -> dict:
+    """
+    Perform a complete design compliance check.
+    Returns a dictionary summarizing pass/fail status of all rules.
+    """
+    velocity = Velocity(flowrate, diameter).calculate()
+    Re = ReynoldsNumber(velocity, diameter, density, viscosity).calculate()
+    flow_regime = FlowType(Re).calculate()
+    friction_factor = ColebrookWhite(Re, diameter, roughness).calculate()
+    pressure_drop = PressureDrop(flowrate, diameter, length, density, viscosity, roughness).calculate()
+    optimum_d = OptimumPipeDiameter(flowrate, density).calculate()
+
+    return {
+        "flow_regime": flow_regime,
+        "velocity_ok": check_velocity_limits(velocity, fluid_type),
+        "pressure_drop_ok": check_pressure_drop_limit(pressure_drop),
+        "friction_factor": friction_factor,
+        "optimum_diameter_mm": optimum_d,
+    }