Merge pull request #1064 from fernandarossi/dev-lubricant_dict

Updating properties for CFD simulations
petrobras · Oct 17, 2024 · 32c5f1a · 32c5f1a
2 parents 67fbd3a + 88a1e18
commit 32c5f1a
Show file tree

Hide file tree

Showing 4 changed files with 100 additions and 39 deletions.
diff --git a/ross/__init__.py b/ross/__init__.py
@@ -14,6 +14,7 @@
 from .shaft_element import *
 from .units import Q_
 from .utils import get_data_from_figure, visualize_matrix
-from ross.fluid_flow.lubricants import lubricant_dict
+from ross.fluid_flow.lubricants import lubricants_dict
+from ross.fluid_flow.materials import materials_dict
 
 _pio.templates.default = "ross"
diff --git a/ross/fluid_flow/cylindrical.py b/ross/fluid_flow/cylindrical.py
@@ -5,7 +5,7 @@
 from ross.bearing_seal_element import BearingElement
 from ross.units import Q_, check_units
 from scipy.optimize import curve_fit, minimize
-from ross.fluid_flow.lubricants import lubricant_dict
+from ross.fluid_flow.lubricants import lubricants_dict
 
 from plotly import graph_objects as go
 from plotly import figure_factory as ff
@@ -233,15 +233,15 @@ def __init__(
 
         self.oil_flow = self.oil_flow / 60000
 
-        lubricant_properties = lubricant_dict[self.lubricant]
+        lubricant_properties = lubricants_dict[self.lubricant]
 
-        T_muI = Q_(lubricant_properties["temp1"], "degK").m_as("degC")
-        T_muF = Q_(lubricant_properties["temp2"], "degK").m_as("degC")
-        mu_I = lubricant_properties["viscosity1"]
-        mu_F = lubricant_properties["viscosity2"]
-        self.rho = lubricant_properties["lube_density"]
-        self.Cp = lubricant_properties["lube_cp"]
-        self.k_t = lubricant_properties["lube_conduct"]
+        T_muI = Q_(lubricant_properties["temperature1"], "degK").m_as("degC")
+        T_muF = Q_(lubricant_properties["temperature2"], "degK").m_as("degC")
+        mu_I = lubricant_properties["liquid_viscosity1"]
+        mu_F = lubricant_properties["liquid_viscosity2"]
+        self.rho = lubricant_properties["liquid_density"]
+        self.Cp = lubricant_properties["liquid_specific_heat"]
+        self.k_t = lubricant_properties["liquid_thermal_conductivity"]
 
         # Interpolation coefficients
         self.a, self.b = self._interpol(T_muI, T_muF, mu_I, mu_F)

diff --git a/ross/fluid_flow/lubricants.py b/ross/fluid_flow/lubricants.py
@@ -2,41 +2,85 @@
 
 from ross.units import Q_
 
-lubricant_dict = {
+lubricants_dict = {
     "ISOVG32": {
-        "viscosity1": Q_(4.05640e-06, "reyn").to_base_units().m,
-        "temp1": Q_(40.00000, "degC").to_base_units().m,
-        "viscosity2": Q_(6.76911e-07, "reyn").to_base_units().m,
-        "temp2": Q_(100.00000, "degC").to_base_units().m,
-        "lube_density": Q_(873.99629, "kg/m³").to_base_units().m,
-        "lube_cp": Q_(1948.7995685758851, "J/(kg*degK)").to_base_units().m,
-        "lube_conduct": Q_(0.13126, "W/(m*degC)").to_base_units().m,
+        "liquid_viscosity1": Q_(4.05640e-06, "reyn").to_base_units().m,
+        "temperature1": Q_(40.00000, "degC").to_base_units().m,
+        "liquid_viscosity2": Q_(6.76911e-07, "reyn").to_base_units().m,
+        "temperature2": Q_(100.00000, "degC").to_base_units().m,
+        "liquid_density": Q_(873.99629, "kg/m**3").to_base_units().m,
+        "liquid_specific_heat": Q_(1948.7995685758851, "J/(kg*degK)").to_base_units().m,
+        "liquid_thermal_conductivity": Q_(0.13126, "W/(m*degC)").to_base_units().m,
+        "vapour_viscosity": Q_(9e-06, "Pa*s").to_base_units().m,
+        "vapour_density": Q_(0.029, "kg/m**3").to_base_units().m,
+        "vapour_specific_heat": Q_(1000, "J/(kg*degC)").to_base_units().m,
+        "vapour_thermal_conductivity": Q_(0.026, "W/(m*degK)").to_base_units().m,
+        "saturation_pressure": Q_(1e5, "Pa").to_base_units().m,
+        "surface_tension": Q_(0.04, "N/m").to_base_units().m,
+        "molecular_weight": Q_(500, "g/mol").to_base_units().m,
     },
     "ISOVG46": {
-        "viscosity1": Q_(5.757040938820288e-06, "reyn").to_base_units().m,
-        "temp1": Q_(40, "degC").to_base_units().m,
-        "viscosity2": Q_(8.810775697672788e-07, "reyn").to_base_units().m,
-        "temp2": Q_(100, "degC").to_base_units().m,
-        "lube_density": Q_(862.9, "kg/m³").to_base_units().m,
-        "lube_cp": Q_(1950, "J/(kg*degK)").to_base_units().m,
-        "lube_conduct": Q_(0.15, "W/(m*degC)").to_base_units().m,
+        "liquid_viscosity1": Q_(5.757040938820288e-06, "reyn").to_base_units().m,
+        "temperature1": Q_(40, "degC").to_base_units().m,
+        "liquid_viscosity2": Q_(8.810775697672788e-07, "reyn").to_base_units().m,
+        "temperature2": Q_(100, "degC").to_base_units().m,
+        "liquid_density": Q_(862.9, "kg/m**3").to_base_units().m,
+        "liquid_specific_heat": Q_(1950, "J/(kg*degK)").to_base_units().m,
+        "liquid_thermal_conductivity": Q_(0.15, "W/(m*degC)").to_base_units().m,
+        "vapour_viscosity": Q_(9e-06, "Pa*s").to_base_units().m,
+        "vapour_density": Q_(0.029, "kg/m**3").to_base_units().m,
+        "vapour_specific_heat": Q_(1000, "J/(kg*degC)").to_base_units().m,
+        "vapour_thermal_conductivity": Q_(0.026, "W/(m*degK)").to_base_units().m,
+        "saturation_pressure": Q_(1e5, "Pa").to_base_units().m,
+        "surface_tension": Q_(0.04, "N/m").to_base_units().m,
+        "molecular_weight": Q_(500, "g/mol").to_base_units().m,
     },
     "TEST": {
-        "viscosity1": Q_(0.04, "Pa*s").to_base_units().m,
-        "temp1": Q_(40, "degC").to_base_units().m,
-        "viscosity2": Q_(0.01, "Pa*s").to_base_units().m,
-        "temp2": Q_(100, "degC").to_base_units().m,
-        "lube_density": Q_(863.61302696, "kg/m³").to_base_units().m,
-        "lube_cp": Q_(1951.88616, "J/(kg*degK)").to_base_units().m,
-        "lube_conduct": Q_(0.15, "W/(m*degC)").to_base_units().m,
+        "liquid_viscosity1": Q_(0.04, "Pa*s").to_base_units().m,
+        "temperature1": Q_(40, "degC").to_base_units().m,
+        "liquid_viscosity2": Q_(0.01, "Pa*s").to_base_units().m,
+        "temperature2": Q_(100, "degC").to_base_units().m,
+        "liquid_density": Q_(863.61302696, "kg/m**3").to_base_units().m,
+        "liquid_specific_heat": Q_(1951.88616, "J/(kg*degK)").to_base_units().m,
+        "liquid_thermal_conductivity": Q_(0.15, "W/(m*degC)").to_base_units().m,
+        "vapour_viscosity": Q_(9e-06, "Pa*s").to_base_units().m,
+        "vapour_density": Q_(0.029, "kg/m**3").to_base_units().m,
+        "vapour_specific_heat": Q_(1000, "J/(kg*degC)").to_base_units().m,
+        "vapour_thermal_conductivity": Q_(0.026, "W/(m*degK)").to_base_units().m,
+        "saturation_pressure": Q_(1e5, "Pa").to_base_units().m,
+        "surface_tension": Q_(0.04, "N/m").to_base_units().m,
+        "molecular_weight": Q_(500, "g/mol").to_base_units().m,
     },
     "ISOVG68": {
-        "viscosity1": Q_(0.0570722, "Pa*s").to_base_units().m,
-        "temp1": Q_(40, "degC").to_base_units().m,
-        "viscosity2": Q_(0.00766498, "Pa*s").to_base_units().m,
-        "temp2": Q_(100, "degC").to_base_units().m,
-        "lube_density": Q_(867, "kg/m³").to_base_units().m,
-        "lube_cp": Q_(1951, "J/(kg*degK)").to_base_units().m,
-        "lube_conduct": Q_(0.15, "W/(m*degC)").to_base_units().m,
+        "liquid_viscosity1": Q_(0.0570722, "Pa*s").to_base_units().m,
+        "temperature1": Q_(40, "degC").to_base_units().m,
+        "liquid_viscosity2": Q_(0.00766498, "Pa*s").to_base_units().m,
+        "temperature2": Q_(100, "degC").to_base_units().m,
+        "liquid_density": Q_(867, "kg/m**3").to_base_units().m,
+        "liquid_specific_heat": Q_(1951, "J/(kg*degK)").to_base_units().m,
+        "liquid_thermal_conductivity": Q_(0.15, "W/(m*degC)").to_base_units().m,
+        "vapour_viscosity": Q_(9e-06, "Pa*s").to_base_units().m,
+        "vapour_density": Q_(0.029, "kg/m**3").to_base_units().m,
+        "vapour_specific_heat": Q_(1000, "J/(kg*degC)").to_base_units().m,
+        "vapour_thermal_conductivity": Q_(0.026, "W/(m*degK)").to_base_units().m,
+        "saturation_pressure": Q_(1e5, "Pa").to_base_units().m,
+        "surface_tension": Q_(0.04, "N/m").to_base_units().m,
+        "molecular_weight": Q_(500, "g/mol").to_base_units().m,
+    },
+    "HYDRAXP68": {
+        "liquid_viscosity1": Q_(0.0570722, "Pa*s").to_base_units().m,
+        "temperature1": Q_(40, "degC").to_base_units().m,
+        "liquid_viscosity2": Q_(7.66498e-03, "Pa*s").to_base_units().m,
+        "temperature2": Q_(100, "degC").to_base_units().m,
+        "liquid_density": Q_(874, "kg/m**3").to_base_units().m,
+        "liquid_specific_heat": Q_(2000, "J/(kg*degC)").to_base_units().m,
+        "liquid_thermal_conductivity": Q_(0.133, "W/(m*degK)").to_base_units().m,
+        "vapour_viscosity": Q_(9e-06, "Pa*s").to_base_units().m,
+        "vapour_density": Q_(0.029, "kg/m**3").to_base_units().m,
+        "vapour_specific_heat": Q_(1000, "J/(kg*degC)").to_base_units().m,
+        "vapour_thermal_conductivity": Q_(0.026, "W/(m*degK)").to_base_units().m,
+        "saturation_pressure": Q_(81325, "Pa").to_base_units().m,
+        "surface_tension": Q_(0.04, "N/m").to_base_units().m,
+        "molecular_weight": Q_(500, "g/mol").to_base_units().m,
     },
 }
diff --git a/ross/fluid_flow/materials.py b/ross/fluid_flow/materials.py
@@ -0,0 +1,16 @@
+"""This module deals with materials dictionary in the ROSS library."""
+
+from ross.units import Q_
+
+materials_dict = {
+    "steel": {
+        "density": Q_(7850, "kg/m**3").to_base_units().m,
+        "specific_heat": Q_(434, "J/(kg*degC)").to_base_units().m,
+        "thermal_conductivity": Q_(60.5, "W/(m*degK)").to_base_units().m,
+    },
+    "brass": {
+        "density": Q_(8600, "kg/m**3").to_base_units().m,
+        "specific_heat": Q_(380, "J/(kg*degC)").to_base_units().m,
+        "thermal_conductivity": Q_(109, "W/(m*degK)").to_base_units().m,
+    },
+}