def default_library_gases(self):
try:
from ctREFPROP.ctREFPROP import REFPROPFunctionLibrary
self.list_gases = {}
self.fluid_file_to_final_name = {}
refProp_path = os.environ['RPPREFIX']
if os.path.exists(refProp_path):
self.RefProp = REFPROPFunctionLibrary(refProp_path)
self.RefProp.SETPATHdll(refProp_path)
refProp_fluids_path = get_new_path(refProp_path, "FLUIDS")
list_files = os.listdir(refProp_fluids_path)
for fluid_file in list_files:
if ".BNC" not in fluid_file:
filepath = get_new_path(refProp_fluids_path,
fluid_file)
f = open(filepath, 'r')
line_0 = f.readline()
line_1 = f.readline()
line_2 = f.readline()
f.close()
short_name = line_0.split("!")[0]
full_name = line_2.split("!")[0]
letter = " "
while letter == " ":
short_name = short_name[:-1]
letter = short_name[-1]
letter = " "
while letter == " ":
full_name = full_name[:-1]
letter = full_name[-1]
final_name = short_name if short_name == full_name else f"{short_name} ({full_name})"
self.list_gases[final_name] = [
fluid_file, short_name, full_name
]
self.fluid_file_to_final_name[fluid_file] = final_name
else:
title = "REFPROP installation not detected"
message = "Dear user, REFPROP application files were not found in the computer's default paths. "
message += "Please, install the REFPROP on your computer to enable the set-up of the fluids mixture."
PrintMessageInput([title, message, "ERROR"])
return True
except Exception as log_error:
title = "Error while loading REFPROP"
message = "An error has been reached while trying to load REFPROP data. If the REFPROP module has already been "
message += "installed we recommend running the 'pip install ctREFPROP' command at the terminal to install the "
message += "necessary libraries."
message += f"\n\n{str(log_error)}"
PrintMessageInput([title, message, "ERROR"])
return True
def parallel_evaluate(mixture):
"""
Fully encapsulated for multiprocessing support
Fresh copy of REFPROP will be loaded into each instance, for parallelism
"""
root = os.environ['RPPREFIX']
R = REFPROPFunctionLibrary(root)
R.SETPATHdll(root)
components, compositions = mixture
names = '*'.join([f + '.FLD' for f in components])
r = R.REFPROPdll(names, 'TP', 'D', R.MOLAR_BASE_SI, 0, 0, 300, 10e3,
compositions)
if r.ierr == 0:
return [r.Output[0]]
else:
return r.herr
return out
def parallel_evaluate(mixtures):
"""
Fully encapsulated for multiprocessing support
Fresh copy of REFPROP will be loaded into each instance, for parallelism
"""
root = os.environ['RPPREFIX']
R = REFPROPFunctionLibrary(os.path.join(root, 'REFPRP64.dll')) # Change this as needed
R.SETPATHdll(root)
out = []
names, compositions = mixtures
ierr, herr = R.SETUPdll(2, '|'.join([f+'.FLD' for f in names]), 'HMX.BNC', 'DEF')
if ierr != 0:
print(ierr, herr)
return herr
kphase = 2 # Vapor
kguess = 0 # Don't use guesses
rho, ierr, herr = R.TPRHOdll(300.0, 10.0, compositions, kphase, kguess,-1)
out.append(rho)
return out
def time_REFPROP(n, Nrep):
root = os.environ['RPPREFIX']
RP = REFPROPFunctionLibrary(root)
RP.SETFLUIDSdll('METHANE')
if n == 0:
f = lambda D, T: RP.REFPROP1dll(D, T)
elif n == 1:
f = lambda AS: AS.delta() * AS.dalphar_dDelta()
elif n == 2:
f = lambda AS: AS.delta()**2 * AS.d2alphar_dDelta2()
elif n == 3:
f = lambda AS: AS.delta()**3 * AS.d3alphar_dDelta3()
else:
raise ValueError(n)
T = 300
rho = 3.0
z = [1.0]
tic = timeit.default_timer()
for i in range(Nrep):
RP.PRESSdll(T, rho / 1e3, z)
toc = timeit.default_timer()
elap = (toc - tic) / Nrep
return elap
def RiemFluid_TP(T, P, *, fluid):
from ctREFPROP.ctREFPROP import REFPROPFunctionLibrary
PATH = r'D:\REFPROP'
RP = REFPROPFunctionLibrary(PATH)
MOLAR_BASE_SI = RP.GETENUMdll(0, 'MOLAR BASE SI').iEnum
RP.SETPATHdll(PATH)
RP.SETFLUIDSdll(fluid)
rhomolar, Riem = RP.REFPROPdll('', 'TP', 'D;RIEM', MOLAR_BASE_SI, 0, 0, T,
P, [1.0]).Output[0:2]
rhoN = rhomolar * N_A
z = [1.0] + [0.0] * 19
Ar00, Ar01, Ar10, Ar02, Ar11, Ar20, Ar03, Ar12, Ar21, Ar30 = RP.REFPROPdll(
'', 'TD',
'PHIR00;PHIR01;PHIR10;PHIR02;PHIR11;PHIR20;PHIR03;PHIR12;PHIR21;PHIR30',
MOLAR_BASE_SI, 0, 0, T, rhomolar, z).Output[0:10]
Ao00, Ao01, Ao10, Ao02, Ao11, Ao20, Ao03, Ao12, Ao21, Ao30 = RP.REFPROPdll(
'', 'TD',
'PHIG00;PHIG01;PHIG10;PHIG02;PHIG11;PHIG20;PHIG03;PHIG12;PHIG21;PHIG30',
MOLAR_BASE_SI, 0, 0, T, rhomolar, z).Output[0:10]
Bstar = 1 + 2 * Ar01 + Ar02
cvstar = -Ao20 - Ar20
# First derivatives
TdBdT = -(2 * Ar11 + Ar12)
TdcvdT = ((2 * Ao20 + Ao30) + (2 * Ar20 + Ar30))
rhodBdrho = (2 * (Ar01 + Ar02) + (2 * Ar02 + Ar03))
rhodcvdrho = -Ar21
# Combinations of first derivatives
rhodlnBcvdrho = 1 / Bstar * rhodBdrho + 1 / cvstar * rhodcvdrho
TdlnBcvdT = 1 / Bstar * TdBdT + 1 / cvstar * TdcvdT
rhodlnBovercvdrho = 1 / Bstar * rhodBdrho - 1 / cvstar * rhodcvdrho
# Second derivatives, same treatment as first derivatives,
# Fourth derivatives cancel conveniently
T2dBdT2_plus_rho2dcvdrho2 = (4 * Ar11 + 2 * Ar12 + 2 * Ar21)
Rstar = 1 / (Bstar * cvstar) * (TdBdT *
(1 - 0.5 * TdlnBcvdT) + rhodcvdrho *
(1 - 0.5 * rhodlnBcvdrho) -
0.5 * cvstar * rhodlnBovercvdrho +
T2dBdT2_plus_rho2dcvdrho2)
here = Rstar / rhoN # [m^3]
return here * 1e27, Riem # [nm^3]
class SetFluidCompositionInput(QDialog):
def __init__(self,
project,
opv,
selected_fluid_to_edit=None,
*args,
**kwargs):
super().__init__()
uic.loadUi(
'data/user_input/ui/model/setup/acoustic/setFluidCompositionInput.ui',
self)
self.setWindowFlags(Qt.WindowStaysOnTopHint)
self.setWindowModality(Qt.WindowModal)
icons_path = 'data\\icons\\'
self.icon_pulse = QIcon(icons_path + 'pulse.png')
self.setWindowIcon(self.icon_pulse)
self.icon_animate = QIcon(icons_path + 'play_pause.png')
self.project = project
self.opv = opv
self.opv.setInputObject(self)
self.selected_fluid_to_edit = selected_fluid_to_edit
self.compressor_info = kwargs.get("compressor_info", {})
self.save_path = ""
self.export_file_path = ""
self.userPath = os.path.expanduser('~')
self.fluid_path = project.get_fluid_list_path()
self.map_properties = {
"D": "fluid density",
"CP": "specific heat Cp",
"CV": "specific heat Cv",
"CP/CV": "isentropic exponent",
"W": "speed of sound",
"VIS": "dynamic viscosity",
"TCX": "thermal conductivity",
"PRANDTL": "Prandtl number",
"TD": "thermal diffusivity",
"KV": "kinematic viscosity",
"M": "molar mass"
}
self.selected_fluid = ""
self.str_composition_value = ""
self.unit_temperature = "K"
self.unit_pressure = "Pa"
self.composition_value = 0
self.remaining_composition = 1
self.list_fluids = []
self.fluid_to_composition = {}
self.fluid_states = {}
self.all_fluid_state_properties = {}
self.errors_by_fluid_state = {}
self.complete = False
self.state_index = None
self.label_selected_fluid = self.findChild(QLabel,
'label_selected_fluid')
self.label_title_remaining_fraction = self.findChild(
QLabel, 'label_title_remaining_fraction')
self.label_remaining_composition = self.findChild(
QLabel, 'label_remaining_composition')
self.lineEdit_composition = self.findChild(QLineEdit,
'lineEdit_composition')
self.lineEdit_fluid_name = self.findChild(QLineEdit,
'lineEdit_fluid_name')
self.label_fluid_state_index = self.findChild(
QLabel, 'label_fluid_state_index')
self.lineEdit_temperature = self.findChild(QLineEdit,
'lineEdit_temperature')
self.lineEdit_pressure = self.findChild(QLineEdit, 'lineEdit_pressure')
self.lineEdit_temperature_test = self.findChild(
QLineEdit, 'lineEdit_temperature_test')
self.lineEdit_pressure_test = self.findChild(QLineEdit,
'lineEdit_pressure_test')
self.label_fluid_temperature = self.findChild(
QLabel, 'label_fluid_temperature')
self.label_fluid_pressure = self.findChild(QLabel,
'label_fluid_pressure')
self.label_temperature_unit = self.findChild(QLabel,
'label_temperature_unit')
self.label_pressure_unit = self.findChild(QLabel,
'label_pressure_unit')
self.label_fluid_density = self.findChild(QLabel,
'label_fluid_density')
self.label_fluid_specific_heat_Cp = self.findChild(
QLabel, 'label_fluid_specific_heat_Cp')
self.label_fluid_specific_heat_Cv = self.findChild(
QLabel, 'label_fluid_specific_heat_Cv')
self.label_fluid_isentropic_coefficient = self.findChild(
QLabel, 'label_fluid_isentropic_coefficient')
self.label_fluid_speed_of_sound = self.findChild(
QLabel, 'label_fluid_speed_of_sound')
self.label_fluid_dynamic_viscosity = self.findChild(
QLabel, 'label_fluid_dynamic_viscosity')
self.label_fluid_thermal_conductivity = self.findChild(
QLabel, 'label_fluid_thermal_conductivity')
self.comboBox_temperature_units = self.findChild(
QComboBox, 'comboBox_temperature_units')
self.comboBox_pressure_units = self.findChild(
QComboBox, 'comboBox_pressure_units')
self.comboBox_temperature_units_test = self.findChild(
QComboBox, 'comboBox_temperature_units_test')
self.comboBox_pressure_units_test = self.findChild(
QComboBox, 'comboBox_pressure_units_test')
self.comboBox_temperature_units_test.currentIndexChanged.connect(
self.update_state_treeWidget_info)
self.comboBox_pressure_units_test.currentIndexChanged.connect(
self.update_state_treeWidget_info)
self.pushButton_confirm = self.findChild(QPushButton,
'pushButton_confirm')
self.pushButton_confirm.clicked.connect(self.get_fluid_properties)
self.pushButton_reset_fluid = self.findChild(QPushButton,
'pushButton_reset_fluid')
self.pushButton_reset_fluid.clicked.connect(self.reset_fluid)
self.pushButton_add_gas = self.findChild(QPushButton,
'pushButton_add_gas')
self.pushButton_add_gas.clicked.connect(self.add_selected_gas)
self.pushButton_remove_gas = self.findChild(QPushButton,
'pushButton_remove_gas')
self.pushButton_remove_gas.clicked.connect(self.remove_selected_gas)
self.pushButton_add_fluid_state = self.findChild(
QPushButton, 'pushButton_add_fluid_state')
self.pushButton_add_fluid_state.clicked.connect(self.add_fluid_state)
self.pushButton_remove_fluid_state = self.findChild(
QPushButton, 'pushButton_remove_fluid_state')
self.pushButton_remove_fluid_state.clicked.connect(
self.remove_fluid_state)
self.pushButton_get_fluid_properties_info = self.findChild(
QPushButton, 'pushButton_get_fluid_properties_info')
self.pushButton_get_fluid_properties_info.clicked.connect(
self.get_fluid_properties_by_state)
# self.pushButton_use_remaining_molar_fraction = self.findChild(QPushButton, 'pushButton_use_remaining_molar_fraction')
# self.pushButton_use_remaining_molar_fraction.clicked.connect(self.use_remaining_molar_fraction)
self.tabWidget_general = self.findChild(QTabWidget,
'tabWidget_general')
self.tab_fluid_setup = self.tabWidget_general.findChild(
QWidget, 'tab_main')
self.tab_pretest_analysis = self.tabWidget_general.findChild(
QWidget, 'tab_export')
self.tabWidget_general.currentChanged.connect(
self.update_state_treeWidget_info)
self.treeWidget_reference_gases = self.findChild(
QTreeWidget, 'treeWidget_reference_gases')
self.treeWidget_reference_gases.itemClicked.connect(
self.on_click_item_reference_gases)
self.treeWidget_new_gas = self.findChild(QTreeWidget,
'treeWidget_new_gas')
self.treeWidget_new_gas.itemClicked.connect(self.on_click_item_new_gas)
self.treeWidget_new_gas.setColumnWidth(0, 376)
self.treeWidget_new_gas.itemDoubleClicked.connect(
self.on_double_click_item_new_gas)
self.treeWidget_fluids_states = self.findChild(
QTreeWidget, 'treeWidget_fluids_states')
self.treeWidget_fluids_states.setColumnWidth(0, 60)
self.treeWidget_fluids_states.setColumnWidth(1, 120)
self.treeWidget_fluids_states.setColumnWidth(2, 120)
self.treeWidget_fluids_states.itemClicked.connect(
self.on_click_item_fluid_state)
if self.compressor_info:
self.check_compressor_inputs()
self.update_remainig_composition()
if self.default_library_gases():
return
self.load_default_gases_info()
self.update_selected_fluid()
self.exec()
def keyPressEvent(self, event):
if event.key() == Qt.Key_Enter or event.key() == Qt.Key_Return:
self.get_fluid_properties()
elif event.key() == Qt.Key_Backspace or event.key() == Qt.Key_Delete:
self.remove_selected_gas()
elif event.key() == Qt.Key_Escape:
self.close()
def check_compressor_inputs(self):
font = QFont()
font.setFamily("Arial")
font.setPointSize(12)
font.setBold(True)
font.setWeight(75)
self.lineEdit_temperature.setFont(font)
self.lineEdit_temperature.setStyleSheet("color: rgb(0, 0, 255);")
self.lineEdit_pressure.setFont(font)
self.lineEdit_pressure.setStyleSheet("color: rgb(0, 0, 255);")
self.connection_type_comp = self.compressor_info['connection type']
self.connection_label = "discharge" if self.connection_type_comp else "suction"
self.temperature_comp = self.compressor_info[
f'temperature ({self.connection_label})']
self.pressure_comp = self.compressor_info[
f'pressure ({self.connection_label})']
self.line_id_comp = self.compressor_info['line_id']
self.node_id_comp = self.compressor_info['node_id']
self.isentropic_exponent_comp = self.compressor_info[
'isentropic exponent']
self.pressure_ratio_comp = self.compressor_info['pressure ratio']
self.molar_mass_comp = self.compressor_info['molar mass']
self.k = self.compressor_info['isentropic exponent']
self.p_ratio = self.compressor_info['pressure ratio']
self.T_suc = self.compressor_info[f'temperature (suction)']
self.lineEdit_temperature.setText(str(round(self.temperature_comp, 4)))
self.lineEdit_temperature.setDisabled(True)
self.lineEdit_pressure.setText(str(round(self.pressure_comp, 4)))
self.lineEdit_pressure.setDisabled(True)
def update_selected_fluid(self):
if self.selected_fluid_to_edit:
[fluid_name, temperature, pressure, key_mixture,
molar_fractions] = self.selected_fluid_to_edit
fluid_file_names = key_mixture.split(";")
self.lineEdit_fluid_name.setText(fluid_name)
self.lineEdit_temperature.setText(str(temperature))
self.lineEdit_pressure.setText(str(pressure))
for index, fluid_file_name in enumerate(fluid_file_names):
final_name = self.fluid_file_to_final_name[fluid_file_name]
str_molar_fraction = str(round(molar_fractions[index] * 100,
5))
self.fluid_to_composition[final_name] = [
str_molar_fraction, molar_fractions[index], fluid_file_name
]
self.load_new_gas_composition_info()
self.update_remainig_composition()
def use_remaining_molar_fraction(self):
self.lineEdit_composition.setText(str(self.remaining_composition))
def add_selected_gas(self):
if self.label_selected_fluid.text() != "":
if self.check_composition_input():
return
self.load_new_gas_composition_info()
self.update_remainig_composition()
else:
title = "None 'Fluid' selected"
message = "Dear user, it is necessary to select a fluid in the list to proceed"
PrintMessageInput([title, message, window_title_1])
def update_remainig_composition(self):
self.remaining_composition = 1
for [_, composition_value, _] in self.fluid_to_composition.values():
self.remaining_composition -= composition_value
if round(abs(self.remaining_composition), 5) > 0:
self.label_remaining_composition.setVisible(True)
self.label_title_remaining_fraction.setVisible(True)
_remain = round(self.remaining_composition * 100, 5)
self.label_remaining_composition.setText(str(_remain))
else:
self.label_remaining_composition.setText("")
self.label_remaining_composition.setVisible(False)
self.label_title_remaining_fraction.setVisible(False)
def create_font(self, size):
self.font = QFont()
self.font.setPointSize(size)
self.font.setItalic(False)
# self.font.setBold(True)
# self.font.setWeight(75)
def remove_selected_gas(self):
if self.label_selected_fluid.text() != "":
_fluid = self.label_selected_fluid.text()
if _fluid in self.fluid_to_composition.keys():
self.fluid_to_composition.pop(_fluid)
self.load_new_gas_composition_info()
self.update_remainig_composition()
def reset_fluid(self):
title = f"Resetting of the current 'Fluid Composition'"
message = "Do you really want to reset the current Fluid Composition?\n\n"
message += "\n\nPress the Continue button to proceed with the resetting or press Cancel or "
message += "\nClose buttons to abort the current operation."
read = CallDoubleConfirmationInput(title,
message,
leftButton_label='Cancel',
rightButton_label='Continue')
if read._stop:
return
self.fluid_to_composition.clear()
self.load_new_gas_composition_info()
self.update_remainig_composition()
def load_default_gases_info(self):
self.treeWidget_reference_gases.clear()
self.treeWidget_reference_gases.setGeometry(10, 142, 376, 400)
self.treeWidget_reference_gases.headerItem().setText(
0, "Default fluid library")
for gas in self.list_gases.keys():
new = QTreeWidgetItem([gas])
new.setTextAlignment(0, Qt.AlignCenter)
self.treeWidget_reference_gases.addTopLevelItem(new)
def load_new_gas_composition_info(self):
# if self.selected_fluid != "":
self.treeWidget_new_gas.clear()
self.treeWidget_new_gas.setGeometry(576, 142, 509, 400)
self.treeWidget_new_gas.headerItem().setText(0, "Fluid")
self.treeWidget_new_gas.headerItem().setText(1, "Composition [%]")
for fluid, [str_composition, _,
_] in self.fluid_to_composition.items():
new = QTreeWidgetItem([fluid, str_composition])
new.setTextAlignment(0, Qt.AlignCenter)
new.setTextAlignment(1, Qt.AlignCenter)
self.treeWidget_new_gas.addTopLevelItem(new)
self.label_selected_fluid.setText("")
self.lineEdit_composition.setText("")
def check_composition_input(self):
self.str_composition_value = self.lineEdit_composition.text()
if self.str_composition_value != "":
title = "Invalid input value to the fluid Composition"
message = ""
value = None
try:
value = float(self.str_composition_value)
except Exception as log_error:
message = "Dear user, you have typed an invalid entry at the fluid Composition input. "
message += "\nPlease, check the typed value to proceed with the fluid setup.\n\n"
message += str(log_error)
if value is not None:
if value > 100 or value < 0:
message = "Dear user, you have typed an invalid entry at the fluid Composition input. "
message += "The value must be a positive value less or equals to 100."
message += "\nPlease, check the typed value to proceed with the fluid setup."
if round(value / 100, 5) > round(self.remaining_composition, 5):
_remain = round(self.remaining_composition * 100, 5)
message = "Dear user, you have typed an invalid entry at the Fluid Composition input. "
message += f"The value must be a positive value less or equals to {_remain}%."
message += "\nPlease, check the typed value to proceed with the fluid setup."
if message == "":
self.composition_value = value / 100
fluid_file_name, _, _ = self.list_gases[self.selected_fluid]
self.fluid_to_composition[self.selected_fluid] = [
self.str_composition_value, self.composition_value,
fluid_file_name
]
if self.composition_value == 0:
if self.selected_fluid in self.fluid_to_composition.keys():
self.fluid_to_composition.pop(self.selected_fluid)
return False
else:
PrintMessageInput([title, message, window_title_1])
return True
else:
title = "Empty entry at molar fraction input field"
message = "An Empty entry has been detected at molar fraction input field. "
message += "You should to inform a valid positive float number less than 1 to proceed."
PrintMessageInput([title, message, window_title_1])
return True
def update_label_selected_fluid_font(self):
if len(self.selected_fluid) < 20:
fontsize = 12
elif len(self.selected_fluid) < 40:
fontsize = 11
elif len(self.selected_fluid) < 60:
fontsize = 10
elif len(self.selected_fluid) < 70:
fontsize = 9
else:
fontsize = 8
self.create_font(fontsize)
self.label_selected_fluid.setFont(self.font)
def on_click_item_reference_gases(self, item):
self.selected_fluid = item.text(0)
self.label_selected_fluid.setText(self.selected_fluid)
self.update_label_selected_fluid_font()
def on_double_click_item_new_gas(self, item):
return
def on_click_item_new_gas(self, item):
self.selected_fluid = item.text(0)
self.label_selected_fluid.setText(item.text(0))
self.lineEdit_composition.setText(item.text(1))
self.update_label_selected_fluid_font()
def get_fluid_properties(self):
message = ""
self.fluid_setup = []
self.errors = {}
if round(self.remaining_composition, 5) == 0:
if self.lineEdit_fluid_name.text() != "":
self.fluid_properties = {}
units = self.RefProp.GETENUMdll(0, "MASS BASE SI").iEnum
fluids_string = ""
molar_fractions = []
for _, _fraction, file_name in self.fluid_to_composition.values(
):
fluids_string += file_name + ";"
molar_fractions.append(_fraction)
fluids_string = fluids_string[:-1]
self.unit_temperature_update(self.comboBox_temperature_units)
self.unit_pressure_update(self.comboBox_pressure_units)
values = self.check_input_values_with_units(
self.lineEdit_temperature, self.lineEdit_pressure)
if values is None:
return
else:
[temperature_K, pressure_Pa] = values
self.fluid_properties["temperature"] = temperature_K
self.fluid_properties["pressure"] = pressure_Pa
self.fluid_properties[
"fluid name"] = self.lineEdit_fluid_name.text()
# self.fluid_properties["id"] = ""
# self.fluid_properties["color"] = ""
if self.compressor_info:
if self.connection_label == "discharge":
count = 0
criteria = 100
cache_temperatures = [temperature_K]
while criteria > 0.001 and count <= 100:
for key_prop in [
"D", "CV", "CP", "CP/CV", "W", "VIS",
"TCX", "M"
]: #, "PRANDTL", "TD", "KV"]:
read = self.RefProp.REFPROPdll(
fluids_string, "TP", key_prop, units, 0, 0,
temperature_K, pressure_Pa,
molar_fractions)
if read.herr:
self.errors[self.map_properties[
key_prop]] = read.herr
self.fluid_properties[self.map_properties[
key_prop]] = read.Output[0]
if key_prop == "CP/CV":
k_iter = read.Output[0]
count += 1
temperature_K_iter = self.T_suc * (self.p_ratio**(
(k_iter - 1) / k_iter))
cache_temperatures.append(temperature_K_iter)
criteria = abs(cache_temperatures[-1] -
cache_temperatures[-2]) / (
(cache_temperatures[-1] +
cache_temperatures[-2]) / 2)
temperature_K = temperature_K_iter
self.fluid_properties[
"temperature"] = temperature_K
# print(count, k_iter, cache_temperatures[-1], cache_temperatures[-2], criteria)
elif self.connection_label == "suction":
for key_prop in [
"D", "CV", "CP", "CP/CV", "W", "VIS", "TCX",
"M"
]: #, "PRANDTL", "TD", "KV"]:
read = self.RefProp.REFPROPdll(
fluids_string, "TP", key_prop, units, 0, 0,
temperature_K, pressure_Pa, molar_fractions)
if read.herr:
self.errors[
self.map_properties[key_prop]] = read.herr
self.fluid_properties[
self.map_properties[key_prop]] = read.Output[0]
else:
for key_prop in [
"D", "CV", "CP", "CP/CV", "W", "VIS", "TCX", "M"
]: #, "PRANDTL", "TD", "KV"]:
read = self.RefProp.REFPROPdll(fluids_string, "TP",
key_prop, units, 0, 0,
temperature_K,
pressure_Pa,
molar_fractions)
if read.herr:
self.errors[
self.map_properties[key_prop]] = read.herr
self.fluid_properties[
self.map_properties[key_prop]] = read.Output[0]
self.fluid_properties["impedance"] = round(
self.fluid_properties["fluid density"] *
self.fluid_properties["speed of sound"], 6)
self.fluid_setup = [fluids_string, molar_fractions]
if self.process_errors():
return
self.complete = True
self.close()
# self.actions_to_finalize()
else:
title = "Additional input required"
message = "Define a fluid name at specific input field to proceed."
self.lineEdit_fluid_name.setFocus()
else:
title = "Fluid composition not finished"
message = "Dear user, you should to complete the fluid mixture composition to proceed.\n"
message += "The sum of all fluids molar fractions must be equals to 1. It is recommended "
message += "to check the inserted molar fractions until this requirement is met."
if message != "":
PrintMessageInput([title, message, "ERROR"])
def process_errors(self):
if len(self.errors) != 0:
title = "REFPROP: Error while processing fluid properties"
message = "The following errors were found while processing the fluid properties.\n\n"
for key, _error in self.errors.items():
message += f"{str(key)}: {str(_error)}\n\n"
message += "It is recommended to check the fluid composition and state properties to proceed."
PrintMessageInput([title, message, "ERROR"], fontsizes=[13, 12])
return True
def get_fluid_state_index(self):
index = len(self.fluid_states) + 1
if str(index) in self.fluid_states.keys():
index = 1
while str(index) in self.fluid_states.keys():
index += 1
self.state_index = str(index)
def actions_to_finalize(self):
if self.compressor_info:
if self.compressor_info["connection type"] == 1:
title = "Fluid properties convergence"
message = "The following fluid properties were obtained after completing the iterative updating process:"
message += f"\n\nTemperature (discharge) = {round(self.fluid_properties['temperature'],4)} [K]"
message += f"\nIsentropic exponent = {round(self.fluid_properties['isentropic exponent'],6)} [-]"
message += "\n\nReference fluid properties:"
message += f"\n\nTemperature (suction) = {self.compressor_info['temperature (suction)']} [K]"
message += f"\nPressure (suction) = {self.compressor_info['pressure (suction)']} [Pa]"
message += f"\nPressure (discharge) = {round(self.compressor_info['pressure (discharge)'],4)} [Pa]"
message += f"\nMolar mass = {round(self.fluid_properties['molar mass'],6)} [kg/mol]"
PrintMessageInput([title, message, "WARNING"])
def add_fluid_state(self):
self.get_fluid_state_index()
self.unit_temperature_update(self.comboBox_temperature_units_test)
self.unit_pressure_update(self.comboBox_pressure_units_test)
values = self.check_input_values_with_units(
self.lineEdit_temperature_test, self.lineEdit_pressure_test)
if values is None:
return
else:
if values not in self.fluid_states.values():
self.fluid_states[self.state_index] = values
self.run_pretest_analysis()
def check_temperature_value(self, lineEdit_temperature):
temperature = None
str_value = lineEdit_temperature.text()
if str_value != "":
try:
temperature = float(str_value)
except Exception as log_error:
title = "Invalid entry to the temperature"
message = "Dear user, you have typed an invalid value at the temperature input field."
message += "You should to inform a valid float number to proceed."
window_title = "ERROR"
PrintMessageInput([title, message, window_title])
else:
title = "Empty temperature input field"
message = "Dear user, the temperature input field is empty. Please, inform a valid float number to proceed."
window_title = "ERROR"
PrintMessageInput([title, message, window_title])
lineEdit_temperature.setFocus()
return temperature
def check_pressure_value(self, lineEdit_pressure):
pressure = None
str_value = lineEdit_pressure.text()
if str_value != "":
try:
pressure = float(str_value)
except Exception as log_error:
title = "Invalid entry to the pressure"
message = "Dear user, you have typed an invalid value at the pressure input field."
message += "You should to inform a valid float number to proceed."
window_title = "ERROR"
PrintMessageInput([title, message, window_title])
else:
title = "Empty pressure input field"
message = "Dear user, the pressure input field is empty. Please, inform a valid float number to proceed."
window_title = "ERROR"
PrintMessageInput([title, message, window_title])
lineEdit_pressure.setFocus()
return pressure
def remove_fluid_state(self):
if self.state_index in self.fluid_states.keys():
self.fluid_states.pop(self.state_index)
self.update_state_treeWidget_info()
def unit_temperature_update(self, comboBox_temperature_units):
temperature_unit_labels = ["K", "°C", "°F"]
index_temperature = comboBox_temperature_units.currentIndex()
self.unit_temperature = temperature_unit_labels[index_temperature]
def unit_pressure_update(self, comboBox_pressure_units):
pressure_unit_labels = ["Pa", "kPa", "bar", "psi"]
index_pressure = comboBox_pressure_units.currentIndex()
self.unit_pressure = pressure_unit_labels[index_pressure]
def check_input_values_with_units(self, lineEdit_temperature,
lineEdit_pressure):
_temperature_value = self.check_temperature_value(lineEdit_temperature)
if _temperature_value is None:
return None
if self.unit_temperature == "°C":
_temperature_value += 273.15
elif self.unit_temperature == "°F":
_temperature_value = (_temperature_value - 32) * (5 / 9) + 273.15
if _temperature_value < 0:
title = "Invalid entry to the temperature"
message = "The typed value at temperature input field reaches a negative value in Kelvin scale."
message += "It is necessary to enter a value that maintains the physicall coherence and consistence "
message += "to proceed with the fluid setup."
PrintMessageInput([title, message, window_title_1])
return None
_pressure_value = self.check_pressure_value(lineEdit_pressure)
if _pressure_value is None:
return None
if self.unit_pressure == "kPa":
_pressure_value *= 1e3
elif self.unit_pressure == "bar":
_pressure_value *= 101325
elif self.unit_pressure == "psi":
_pressure_value *= 6894.75729
if _pressure_value < 0:
title = "Invalid entry to the pressure"
message = "The typed value at pressure input field reaches a negative value in Pascal scale."
message += "It is necessary to enter a value that maintains the physicall coherence and consistence "
message += "to proceed with the fluid setup."
PrintMessageInput([title, message, window_title_1])
return None
return [round(_temperature_value, 5), round(_pressure_value, 5)]
def update_fluid_state_header(self):
self.unit_temperature_update(self.comboBox_temperature_units_test)
self.unit_pressure_update(self.comboBox_pressure_units_test)
self.treeWidget_fluids_states.headerItem().setText(0, f"Index")
self.treeWidget_fluids_states.headerItem().setText(
1, f"Temperature [{self.unit_temperature}]")
self.treeWidget_fluids_states.headerItem().setText(
2, f"Pressure [{self.unit_pressure}]")
self.treeWidget_fluids_states.headerItem().setText(
3, "Set of properties")
for i in range(4):
self.treeWidget_fluids_states.headerItem().setTextAlignment(
i, Qt.AlignCenter)
def update_state_treeWidget_info(self):
if self.tabWidget_general.currentIndex() == 1:
self.update_fluid_state_header()
self.treeWidget_fluids_states.clear()
# self.treeWidget_fluids_states.setGeometry(592, 68, 509, 400)
for index, [temperature, pressure] in self.fluid_states.items():
if self.unit_temperature == "°C":
temperature -= 273.15
elif self.unit_temperature == "°F":
temperature = (temperature - 273.15) * (9 / 5) + 32
if self.unit_pressure == "kPa":
pressure *= 1e-3
elif self.unit_pressure == "bar":
pressure *= 1e-5
elif self.unit_pressure == "psi":
pressure /= 6894.75729
if self.errors_by_fluid_state:
if self.errors_by_fluid_state[index] != 0:
status = "invalid"
else:
status = "valid"
else:
status = "--"
new = QTreeWidgetItem([
index,
str(round(temperature, 6)),
str(round(pressure, 6)), status
])
for i in range(5):
new.setTextAlignment(i, Qt.AlignCenter)
self.treeWidget_fluids_states.addTopLevelItem(new)
self.state_index = None
self.lineEdit_temperature.setText("")
self.lineEdit_pressure.setText("")
self.label_fluid_state_index.setText("")
self.pushButton_reset_fluid.setVisible(False)
self.pushButton_confirm.setVisible(False)
self.reset_fluid_properties_labels()
else:
self.pushButton_reset_fluid.setVisible(True)
self.pushButton_confirm.setVisible(True)
def on_click_item_fluid_state(self, item):
str_index = item.text(0)
if str_index != "":
self.label_fluid_state_index.setText(str_index)
self.state_index = str_index
self.fluid_temperature = item.text(1)
self.fluid_pressure = item.text(2)
def run_pretest_analysis(self):
self.errors_by_fluid_state = {}
self.all_fluid_state_properties = {}
for index, [temperature_K, pressure_Pa] in self.fluid_states.items():
# message = ""
if round(self.remaining_composition, 5) == 0:
# self.fluid_properties = {}
units = self.RefProp.GETENUMdll(0, "MASS BASE SI").iEnum
fluids_string = ""
molar_fractions = []
for _, _fraction, file_name in self.fluid_to_composition.values(
):
fluids_string += file_name + ";"
molar_fractions.append(_fraction)
fluids_string = fluids_string[:-1]
fluid_properties_by_state = {}
for key_prop in ["D", "CV", "CP", "CP/CV", "W", "VIS",
"TCX"]: #, "PRANDTL", "TD", "KV"]:
read = self.RefProp.REFPROPdll(fluids_string, "TP",
key_prop, units, 0, 0,
temperature_K, pressure_Pa,
molar_fractions)
if index in self.errors_by_fluid_state.keys():
if read.herr:
self.errors_by_fluid_state[index] += 1
else:
self.errors_by_fluid_state[index] = 0
fluid_properties_by_state[
self.map_properties[key_prop]] = read.Output[0]
self.all_fluid_state_properties[
index] = fluid_properties_by_state
self.update_state_treeWidget_info()
def get_fluid_properties_by_state(self):
if self.all_fluid_state_properties:
fluid_state_properties = self.all_fluid_state_properties[
self.state_index]
self.label_fluid_temperature.setText(str(self.fluid_temperature))
self.label_temperature_unit.setText(f"[{self.unit_temperature}]")
self.label_pressure_unit.setText(f"[{self.unit_pressure}]")
self.label_fluid_pressure.setText(str(self.fluid_pressure))
self.label_fluid_density.setText(
str(round(fluid_state_properties["fluid density"], 6)))
self.label_fluid_specific_heat_Cp.setText(
str(round(fluid_state_properties["specific heat Cp"], 6)))
self.label_fluid_specific_heat_Cv.setText(
str(round(fluid_state_properties["specific heat Cv"], 6)))
self.label_fluid_isentropic_coefficient.setText(
str(round(fluid_state_properties["isentropic exponent"], 6)))
self.label_fluid_speed_of_sound.setText(
str(round(fluid_state_properties["speed of sound"], 6)))
self.label_fluid_dynamic_viscosity.setText(
str(round(fluid_state_properties["dynamic viscosity"], 10)))
self.label_fluid_thermal_conductivity.setText(
str(round(fluid_state_properties["thermal conductivity"], 6)))
def reset_fluid_properties_labels(self):
self.label_fluid_temperature.setText("")
self.label_fluid_pressure.setText("")
self.label_fluid_density.setText("")
self.label_fluid_specific_heat_Cp.setText("")
self.label_fluid_specific_heat_Cv.setText("")
self.label_fluid_isentropic_coefficient.setText("")
self.label_fluid_speed_of_sound.setText("")
self.label_fluid_dynamic_viscosity.setText("")
self.label_fluid_thermal_conductivity.setText("")
def default_library_gases(self):
try:
from ctREFPROP.ctREFPROP import REFPROPFunctionLibrary
self.list_gases = {}
self.fluid_file_to_final_name = {}
refProp_path = os.environ['RPPREFIX']
if os.path.exists(refProp_path):
self.RefProp = REFPROPFunctionLibrary(refProp_path)
self.RefProp.SETPATHdll(refProp_path)
refProp_fluids_path = get_new_path(refProp_path, "FLUIDS")
list_files = os.listdir(refProp_fluids_path)
for fluid_file in list_files:
if ".BNC" not in fluid_file:
filepath = get_new_path(refProp_fluids_path,
fluid_file)
f = open(filepath, 'r')
line_0 = f.readline()
line_1 = f.readline()
line_2 = f.readline()
f.close()
short_name = line_0.split("!")[0]
full_name = line_2.split("!")[0]
letter = " "
while letter == " ":
short_name = short_name[:-1]
letter = short_name[-1]
letter = " "
while letter == " ":
full_name = full_name[:-1]
letter = full_name[-1]
final_name = short_name if short_name == full_name else f"{short_name} ({full_name})"
self.list_gases[final_name] = [
fluid_file, short_name, full_name
]
self.fluid_file_to_final_name[fluid_file] = final_name
else:
title = "REFPROP installation not detected"
message = "Dear user, REFPROP application files were not found in the computer's default paths. "
message += "Please, install the REFPROP on your computer to enable the set-up of the fluids mixture."
PrintMessageInput([title, message, "ERROR"])
return True
except Exception as log_error:
title = "Error while loading REFPROP"
message = "An error has been reached while trying to load REFPROP data. If the REFPROP module has already been "
message += "installed we recommend running the 'pip install ctREFPROP' command at the terminal to install the "
message += "necessary libraries."
message += f"\n\n{str(log_error)}"
PrintMessageInput([title, message, "ERROR"])
return True
import timeit
import json
# https://stackoverflow.com/questions/43356795/can-the-webagg-backend-for-matplotlib-work-with-my-django-site
import matplotlib as mpl
mpl.use('webagg')
# import plotly
# from plotly.offline import init_notebook_mode, plot_mpl
import matplotlib.pyplot as plt
import numpy as np
from ctREFPROP.ctREFPROP import REFPROPFunctionLibrary
root = 'c:/Program Files (x86)/REFPROP'
RP = REFPROPFunctionLibrary(root)
RP.SETPATHdll(root)
def get_PD_data(fluids, Temps, pmax, rhomin, Npts):
RP.SETFLUIDSdll(fluids)
Tt = RP.REFPROPdll("", "TRIP", "T", RP.MOLAR_BASE_SI, 0, 0, 0, 0,
[1.0]).Output[0]
Tc = RP.REFPROPdll("", "CRIT", "T", RP.MOLAR_BASE_SI, 0, 0, 0, 0,
[1.0]).Output[0]
Dmax = RP.REFPROPdll("", "TRIP", "D", RP.MOLAR_BASE_SI, 0, 0, 0, 0,
[1.0]).Output[0]
curves = []
# PEFPROF查找氢气物性模块的封装
# 单位制: kg, m, K, s, mol, J, Pa, W
from ctREFPROP.ctREFPROP import REFPROPFunctionLibrary
RP = REFPROPFunctionLibrary('F:\\1_Refprop\\REFPRP64.DLL', 'dll')
RP.SETPATHdll('F:\\1_Refprop\\REFPROP\\fluids')
r = RP.SETUPdll(1, 'HYDROGEN.FLD', 'HMX.BNC', 'DEF')
assert (r.ierr == 0)
info = RP.INFOdll(1)
def searchHByPD(P, D):
# P: Pa, D: kg/m^3, h: J/kg
sig1 = RP.PDFL1dll(P / 1000, D / info.wmm, [1.0])
return RP.ENTHALdll(sig1.T, D, [1.0]) / info.wmm * 1000
def searchTByDH(D, h):
# D: kg/m^3, h: J/kg, T: K
sig = RP.DHFL1dll(D / info.wmm, h * info.wmm / 1000, [1.0])
return sig.T
def searchPByDH(D, h):
# D: kg/m^3, h: J/kg, p: Pa
sig1 = RP.DHFL1dll(D / info.wmm, h * info.wmm / 1000, [1.0])
return RP.PRESSdll(sig1.T, D / info.wmm, [1.0]) * 1000
def searchDBYPT(P, T):
# D: kg/m^3, p: Pa, T: K
""" A simpified example of calling the PQFLSH function from the legacy version <10.0 API By Ian Bell, 2018, [email protected] """ # Standard library imports import os.sys # pip installable packages from ctREFPROP.ctREFPROP import REFPROPFunctionLibrary root = os.environ['RPPREFIX'] R = REFPROPFunctionLibrary(root) R.SETPATHdll(root) ierr, herr = R.SETUPdll(1, 'WATER.FLD', 'HMX.BNC', 'DEF') if ierr != 0: raise ValueError(str(ierr) + ':' + herr) z = [1.0] p = 101.325 Q = 1.0 o = R.PQFLSHdll(p, Q, z, 1) print(o)
def __init__(self):
RP = REFPROPFunctionLibrary(os.environ['RPPREFIX'])
RP.SETPATHdll(os.environ['RPPREFIX'])
self.UNIT = 7 #RP.GETENUMdll(0,"MKS").iEnum
self.RP = RP