From 7dd18338bddcc26e4ab15c266e553566fff19e5f Mon Sep 17 00:00:00 2001 From: frossi Date: Mon, 3 Jun 2024 13:06:25 -0300 Subject: [PATCH 1/2] Updating properties for CFD simulations --- ross/__init__.py | 3 +- ross/fluid_flow/cylindrical.py | 18 +++---- ross/fluid_flow/lubricants.py | 97 ++++++++++++++++++++++++---------- ross/fluid_flow/materials.py | 16 ++++++ 4 files changed, 95 insertions(+), 39 deletions(-) create mode 100644 ross/fluid_flow/materials.py diff --git a/ross/__init__.py b/ross/__init__.py index b64eb5e8b..c41468be4 100644 --- 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 index 149978a23..825fc717a 100644 --- 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 class THDCylindrical(BearingElement): @@ -230,15 +230,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 index 33173a064..e8820143d 100644 --- a/ross/fluid_flow/lubricants.py +++ b/ross/fluid_flow/lubricants.py @@ -2,41 +2,80 @@ 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_(11325, "Pa").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_(11325, "Pa").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_(11325, "Pa").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_(11325, "Pa").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, + "molecular_weight": Q_(500, "g/mol").to_base_units().m, }, } diff --git a/ross/fluid_flow/materials.py b/ross/fluid_flow/materials.py new file mode 100644 index 000000000..e29b10bca --- /dev/null +++ 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, + }, +} From f30eb1364b21b432375c3385308b7035c62e2a56 Mon Sep 17 00:00:00 2001 From: frossi Date: Mon, 16 Sep 2024 12:39:31 -0300 Subject: [PATCH 2/2] Surface tension coefficient added --- ross/fluid_flow/lubricants.py | 13 +++++++++---- 1 file changed, 9 insertions(+), 4 deletions(-) diff --git a/ross/fluid_flow/lubricants.py b/ross/fluid_flow/lubricants.py index e8820143d..22bdd0975 100644 --- a/ross/fluid_flow/lubricants.py +++ b/ross/fluid_flow/lubricants.py @@ -15,7 +15,8 @@ "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_(11325, "Pa").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": { @@ -30,7 +31,8 @@ "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_(11325, "Pa").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": { @@ -45,7 +47,8 @@ "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_(11325, "Pa").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": { @@ -60,7 +63,8 @@ "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_(11325, "Pa").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": { @@ -76,6 +80,7 @@ "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, }, }