class FluidCharacteristicsView(QWidget, Ui_FluidCharacteristicsForm):
"""
Class to open Molecular Properties Page.
"""
density = """# Density of air
density = 1.293*(273.15 / temperature);
# density for mixtures of gases
#
# Y1 -> mass fraction of component 1
# Y2 -> mass fraction of component 2
rho1 = 1.25051;
rho2 = 1.7832;
A = (Y1 / rho1) + (Y2 /rho2);
density = 1.0 / A;
"""
density_h = """# Density
density = enthalpy / 1040. * 1.29;
# density for mixtures of gases
#
# Y1 -> mass fraction of component 1
# Y2 -> mass fraction of component 2
rho1 = 1.25051;
rho2 = 1.7832;
A = (Y1 / rho1) + (Y2 /rho2);
density = 1.0 / A;
"""
density_wo = """density = 1.25051;
"""
molecular_viscosity = """# Sutherland's Formula
# Gas Cst T0 mu0
# air 120 291.15 18.27e-6
# nitrogen 111 300.55 17.81e-6
# oxygen 127 292.25 20.18e-6
# carbon dioxide 240 293.15 14.8e-6
# carbon monoxide 118 288.15 17.2e-6
# hydrogen 72 293.85 8.76e-6
# ammonia 370 293.15 9.82e-6
# sulfur dioxide 416 293.65 12.54e-6
# helium 79.4 273 19e-6
CST = 120;
T0 = 291.15;
mu_ref = 18.27e-6;
if ( temperature > 0 && temperature < 555) {
molecular_viscosity = mu_ref * ((T0+CST) / (temperature+CST)) * (temperature/T0)^(3./2.);
} else {
molecular_viscosity = -999.0;
}
"""
molecular_viscosity_h = """CST = 120;
T0 = 291.15;
mu_ref = 18.27e-6;
temperature = enthalpy / 1040.;
if ( enthalpy > 0) {
molecular_viscosity = mu_ref * (T0+CST / temperature+CST) * (temperature/T0)^(3./2.);
} else {
molecular_viscosity = -999.0;
}
"""
molecular_viscosity_wo = """CST = 120;
T0 = 291.15;
mu_ref = 18.27e-6;
molecular_viscosity = mu_ref * (T0+CST);
"""
specific_heat = """# specific heat for mixtures of gases
#
# Y1 -> mass fraction of component 1
# Y2 -> mass fraction of component 2
Cp1 = 520.3;
Cp2 = 1040.0;
specific_heat = Y1 * Cp1 + Y2 *Cp2;
"""
volume_viscosity = """# volume_viscosity
"""
thermal_conductivity = """# oxygen
thermal_conductivity = 6.2e-5 * temperature + 8.1e-3;
# nitrogen
thermal_conductivity = 6.784141e-5 * temperature + 5.564317e-3;
# hydrogen
thermal_conductivity = 4.431e-4 * temperature + 5.334e-2;
"""
thermal_conductivity_h = """temperature = enthalpy / 1040.;
thermal_conductivity = 6.2e-5 * temperature + 8.1e-3;
"""
def __init__(self, parent, case):
"""
Constructor
"""
QWidget.__init__(self, parent)
Ui_FluidCharacteristicsForm.__init__(self)
self.setupUi(self)
self.case = case
self.case.undoStopGlobal()
self.mdl = FluidCharacteristicsModel(self.case)
self.notebook = NotebookModel(self.case)
self.m_th = ThermalScalarModel(self.case)
s = self.m_th.getThermalScalarName()
tsm = self.mdl.tsm
# Particular Widget init. taking into account chosen fluid model
mdl_atmo, mdl_joule, mdl_thermal, mdl_gas, mdl_coal, mdl_comp, mdl_hgn=\
self.mdl.getThermoPhysicalModel()
# Combo models
self.modelRho = ComboModel(self.comboBoxRho, 2, 1)
self.modelMu = ComboModel(self.comboBoxMu, 2, 1)
self.modelCp = ComboModel(self.comboBoxCp, 2, 1)
self.modelAl = ComboModel(self.comboBoxAl, 2, 1)
self.modelDiff = ComboModel(self.comboBoxDiff, 2, 1)
self.modelNameDiff = ComboModel(self.comboBoxNameDiff, 1, 1)
self.modelViscv0 = ComboModel(self.comboBoxViscv0, 3, 1)
self.modelDiftl0 = ComboModel(self.comboBoxDiftl0, 3, 1)
self.modelMaterial = ComboModel(self.comboBoxMaterial, 1, 1)
self.modelMethod = ComboModel(self.comboBoxMethod, 1, 1)
self.modelReference = ComboModel(self.comboBoxReference, 1, 1)
if mdl_joule == 'off':
self.modelRho.addItem(self.tr('constant'), 'constant')
self.modelRho.addItem(self.tr('material law'), 'thermal_law')
if mdl_atmo != 'off':
self.modelRho.addItem(self.tr('predefined law'), 'predefined_law')
elif mdl_hgn != "off":
self.modelRho.addItem(self.tr('linear mixture law'),
'predefined_law')
elif mdl_joule != 'off':
self.modelRho.addItem(self.tr('predefined law'), 'predefined_law')
elif mdl_comp != 'off':
self.modelRho.addItem(self.tr('predefined law'), 'predefined_law')
elif mdl_gas != 'off' or mdl_coal != 'off':
self.modelRho.addItem(self.tr('predefined law'), 'predefined_law')
self.modelRho.addItem(self.tr('user law'), 'user_law')
if mdl_joule == 'off':
self.modelMu.addItem(self.tr('constant'), 'constant')
self.modelMu.addItem(self.tr('material law'), 'thermal_law')
if mdl_joule != 'off':
self.modelMu.addItem(self.tr('predefined law'), 'predefined_law')
if mdl_hgn != "off":
self.modelMu.addItem(self.tr('linear mixture law'),
'predefined_law')
self.modelMu.addItem(self.tr('user law'), 'user_law')
if mdl_joule == 'off':
self.modelCp.addItem(self.tr('constant'), 'constant')
self.modelCp.addItem(self.tr('material law'), 'thermal_law')
if mdl_joule != 'off':
self.modelCp.addItem(self.tr('predefined law'), 'predefined_law')
self.modelCp.addItem(self.tr('user law'), 'user_law')
if mdl_joule == 'off':
self.modelAl.addItem(self.tr('constant'), 'constant')
self.modelAl.addItem(self.tr('material law'), 'thermal_law')
if mdl_joule != 'off':
self.modelAl.addItem(self.tr('predefined law'), 'predefined_law')
self.modelAl.addItem(self.tr('user law'), 'user_law')
self.modelDiff.addItem(self.tr('constant'), 'constant')
self.modelDiff.addItem(self.tr('user law'), 'user_law')
self.modelViscv0.addItem(self.tr('constant'), 'constant')
self.modelViscv0.addItem(self.tr('user law'), 'user_law')
self.modelViscv0.addItem(self.tr('material law'), 'thermal_law')
self.modelDiftl0.addItem(self.tr('constant'), 'constant')
self.modelDiftl0.addItem(self.tr('user law'), 'user_law')
self.modelDiftl0.addItem(self.tr('material law'), 'thermal_law')
if mdl_gas != 'off' or mdl_coal != 'off':
self.modelDiftl0.addItem(self.tr('predefined law'),
'predefined_law')
self.scalar = ""
scalar_list = self.mdl.m_sca.getUserScalarNameList()
for s in self.mdl.m_sca.getScalarsVarianceList():
if s in scalar_list: scalar_list.remove(s)
if scalar_list != []:
self.scalar = scalar_list[0]
for scalar in scalar_list:
self.modelNameDiff.addItem(scalar)
# Validators
validatorP0 = DoubleValidator(self.lineEditP0, min=0.0)
self.lineEditP0.setValidator(validatorP0)
validatorT0 = DoubleValidator(self.lineEditT0, min=0.0)
validatorT0.setExclusiveMin(True)
self.lineEditT0.setValidator(validatorT0)
validatorOxydant = DoubleValidator(self.lineEditOxydant, min=0.0)
validatorOxydant.setExclusiveMin(True)
self.lineEditOxydant.setValidator(validatorOxydant)
validatorFuel = DoubleValidator(self.lineEditFuel, min=0.0)
validatorFuel.setExclusiveMin(True)
self.lineEditFuel.setValidator(validatorFuel)
validatorMM = DoubleValidator(self.lineEditMassMolar, min=0.0)
validatorMM.setExclusiveMin(True)
self.lineEditMassMolar.setValidator(validatorMM)
validatorRho = DoubleValidator(self.lineEditRho, min=0.0)
validatorRho1 = DoubleValidator(self.lineEditRho1, min=0.0)
validatorRho2 = DoubleValidator(self.lineEditRho2, min=0.0)
validatorMu = DoubleValidator(self.lineEditMu, min=0.0)
validatorMu1 = DoubleValidator(self.lineEditMu1, min=0.0)
validatorMu2 = DoubleValidator(self.lineEditMu2, min=0.0)
validatorCp = DoubleValidator(self.lineEditCp, min=0.0)
validatorAl = DoubleValidator(self.lineEditAl, min=0.0)
validatorDiff = DoubleValidator(self.lineEditDiff, min=0.0)
validatorViscv0 = DoubleValidator(self.lineEditViscv0, min=0.0)
validatorDiftl0 = DoubleValidator(self.lineEditDiftl0, min=0.0)
validatorRho.setExclusiveMin(True)
validatorRho1.setExclusiveMin(True)
validatorRho2.setExclusiveMin(True)
validatorMu.setExclusiveMin(True)
validatorMu1.setExclusiveMin(True)
validatorMu2.setExclusiveMin(True)
validatorCp.setExclusiveMin(True)
validatorAl.setExclusiveMin(True)
validatorDiff.setExclusiveMin(True)
validatorDiftl0.setExclusiveMin(True)
self.lineEditRho.setValidator(validatorRho)
self.lineEditRho1.setValidator(validatorRho1)
self.lineEditRho2.setValidator(validatorRho2)
self.lineEditMu.setValidator(validatorMu)
self.lineEditMu1.setValidator(validatorMu1)
self.lineEditMu2.setValidator(validatorMu2)
self.lineEditCp.setValidator(validatorCp)
self.lineEditAl.setValidator(validatorAl)
self.lineEditDiff.setValidator(validatorDiff)
self.lineEditViscv0.setValidator(validatorViscv0)
self.lineEditDiftl0.setValidator(validatorDiftl0)
# Connections
self.lineEditP0.textChanged[str].connect(self.slotPressure)
self.lineEditT0.textChanged[str].connect(self.slotTemperature)
self.lineEditOxydant.textChanged[str].connect(self.slotTempOxydant)
self.lineEditFuel.textChanged[str].connect(self.slotTempFuel)
self.lineEditMassMolar.textChanged[str].connect(self.slotMassemol)
self.comboBoxRho.currentIndexChanged[str].connect(self.slotStateRho)
self.comboBoxMu.currentIndexChanged[str].connect(self.slotStateMu)
self.comboBoxCp.currentIndexChanged[str].connect(self.slotStateCp)
self.comboBoxAl.currentIndexChanged[str].connect(self.slotStateAl)
self.comboBoxDiff.currentIndexChanged[str].connect(self.slotStateDiff)
self.comboBoxNameDiff.activated[str].connect(self.slotNameDiff)
self.comboBoxViscv0.currentIndexChanged[str].connect(
self.slotStateViscv0)
self.comboBoxMaterial.activated[str].connect(self.slotMaterial)
self.comboBoxMethod.activated[str].connect(self.slotMethod)
self.comboBoxReference.activated[str].connect(self.slotReference)
self.lineEditRho.textChanged[str].connect(self.slotRho)
self.lineEditRho1.textChanged[str].connect(self.slotRho1)
self.lineEditRho2.textChanged[str].connect(self.slotRho2)
self.lineEditMu.textChanged[str].connect(self.slotMu)
self.lineEditMu1.textChanged[str].connect(self.slotMu1)
self.lineEditMu2.textChanged[str].connect(self.slotMu2)
self.lineEditCp.textChanged[str].connect(self.slotCp)
self.lineEditAl.textChanged[str].connect(self.slotAl)
self.lineEditDiff.textChanged[str].connect(self.slotDiff)
self.lineEditDiftl0.textChanged[str].connect(self.slotDiftl0)
self.lineEditViscv0.textChanged[str].connect(self.slotViscv0)
self.pushButtonRho.clicked.connect(self.slotFormulaRho)
self.pushButtonMu.clicked.connect(self.slotFormulaMu)
self.pushButtonCp.clicked.connect(self.slotFormulaCp)
self.pushButtonAl.clicked.connect(self.slotFormulaAl)
self.pushButtonDiff.clicked.connect(self.slotFormulaDiff)
self.pushButtonViscv0.clicked.connect(self.slotFormulaViscv0)
self.initializeWidget()
self.case.undoStartGlobal()
def initializeWidget(self):
"""
"""
mdls = self.mdl.getThermoPhysicalModel()
mdl_atmo, mdl_joule, mdl_thermal, mdl_gas, mdl_coal, mdl_comp, mdl_hgn = mdls
self.groupBoxMassMolar.hide()
if mdl_atmo != "off":
self.labelInfoT0.hide()
elif mdl_comp != "off" or mdl_coal != "off":
m = self.mdl.getMassemol()
self.lineEditMassMolar.setText(str(m))
self.groupBoxMassMolar.show()
if mdl_thermal != "off" or mdl_gas == 'd3p':
t = self.mdl.getTemperature()
self.lineEditT0.setText(str(t))
if mdl_thermal == "temperature_celsius":
self.labelUnitT0.setText("\xB0 C")
if self.mdl.getMaterials() != "user_material":
self.labelInfoT0.hide()
else:
self.groupBoxTemperature.hide()
if mdl_gas == 'd3p':
self.groupBoxTempd3p.show()
t_oxy = self.mdl.getTempOxydant()
t_fuel = self.mdl.getTempFuel()
self.lineEditOxydant.setText(str(t_oxy))
self.lineEditFuel.setText(str(t_fuel))
else:
self.groupBoxTempd3p.hide()
darc = GroundwaterModel(self.case).getGroundwaterModel()
if darc != 'off':
self.groupBoxPressure.hide()
else:
p = self.mdl.getPressure()
self.lineEditP0.setText(str(p))
if self.mdl.tables == 0 or mdl_joule != 'off' or mdl_comp != 'off':
self.groupBoxTableChoice.hide()
else:
self.groupBoxTableChoice.show()
fluids = self.mdl.getLibPropertiesDict()
o_keys = list(fluids.keys()) # for Python2/3 compatibility
# fluids.keys() ordered and iteratable
# already in Python3
o_keys.sort()
for f in o_keys:
if fluids[f] != 0:
self.modelMaterial.addItem(self.tr(f), f)
else:
self.modelMaterial.addItem(self.tr(f), True)
material = self.mdl.getMaterials()
self.modelMaterial.setItem(str_model=material)
self.updateMethod()
# VoF
if mdl_hgn != 'off':
self.widgetRefRho.hide()
self.widgetVofRho.show()
self.widgetRefMu.hide()
self.widgetVofMu.show()
else:
self.widgetRefRho.show()
self.widgetVofRho.hide()
self.widgetRefMu.show()
self.widgetVofMu.hide()
# compressible
self.groupBoxViscv0.hide()
# combustion
self.groupBoxDiftl0.hide()
if self.scalar == "":
self.groupBoxDiff.hide()
else:
self.groupBoxDiff.show()
self.lineEditDiff.setText(
str(
self.mdl.m_sca.getScalarDiffusivityInitialValue(
self.scalar)))
diff_choice = self.mdl.m_sca.getScalarDiffusivityChoice(
self.scalar)
self.modelDiff.setItem(str_model=diff_choice)
self.modelNameDiff.setItem(str_model=str(self.scalar))
if diff_choice != 'user_law':
self.pushButtonDiff.hide()
self.pushButtonDiff.setEnabled(False)
self.pushButtonDiff.setStyleSheet("background-color: None")
else:
self.pushButtonDiff.show()
self.pushButtonDiff.setEnabled(True)
name = self.mdl.m_sca.getScalarDiffusivityName(self.scalar)
exp = self.mdl.m_sca.getDiffFormula(self.scalar)
if exp:
self.pushButtonDiff.setStyleSheet(
"background-color: green")
self.pushButtonDiff.setToolTip(exp)
else:
self.pushButtonDiff.setStyleSheet("background-color: red")
# Standard Widget initialization
for tag, symbol in self.mdl.lst:
__model = getattr(self, "model" + symbol)
__line = getattr(self, "lineEdit" + symbol)
__button = getattr(self, "pushButton" + symbol)
__label = getattr(self, "label" + symbol)
__labelu = getattr(self, "labelUnit" + symbol)
if tag != 'dynamic_diffusion':
__labelv = getattr(self, "labelVar" + symbol)
c = self.mdl.getPropertyMode(tag)
if c not in __model.getItems():
c = 'constant'
self.mdl.setPropertyMode(tag, c)
__model.setItem(str_model=c)
if c == 'user_law':
__button.show()
__button.setEnabled(True)
__label.setText(self.tr("Reference value"))
else:
__button.hide()
__button.setEnabled(False)
__label.setText(self.tr("Reference value"))
if c == 'thermal_law':
__line.hide()
__label.hide()
__labelu.hide()
__labelv.hide()
else:
__line.show()
__label.show()
__labelu.show()
__labelv.show()
try:
if self.mdl.getMaterials() == "user_material":
__model.disableItem(str_model='thermal_law')
else:
__model.enableItem(str_model='thermal_law')
except Exception:
pass
else:
__label.setText(self.tr("Reference value"))
__line.setText(str(self.mdl.getInitialValue(tag)))
# If no thermal scalar, hide specific heat and thermal conductivity.
if mdl_thermal == 'off':
self.groupBoxCp.hide()
self.groupBoxAl.hide()
if mdl_gas != 'off' or mdl_coal != 'off':
self.groupBoxDiftl0.show()
for tag, symbol in self.mdl.lst:
__model = getattr(self, "model" + symbol)
__line = getattr(self, "lineEdit" + symbol)
__button = getattr(self, "pushButton" + symbol)
__label = getattr(self, "label" + symbol)
__combo = getattr(self, "comboBox" + symbol)
# Gas or coal combustion
if mdl_gas != 'off' or mdl_coal != 'off':
if tag == 'density':
__model.setItem(str_model='predefined_law')
__combo.setEnabled(False)
__button.setEnabled(False)
__button.hide()
self.mdl.setPropertyMode(tag, 'predefined_law')
__label.setText(
self.tr("Calculation by\n perfect gas law"))
__line.setText(str(""))
__line.setEnabled(False)
elif tag == 'dynamic_diffusion':
__model.setItem(str_model='predefined_law')
__combo.setEnabled(False)
__button.setEnabled(False)
__button.hide()
else:
__model.setItem(str_model='constant')
self.mdl.setPropertyMode(tag, 'constant')
# Joule
if mdl_joule == 'arc':
__model.setItem(str_model='predefined_law')
__combo.setEnabled(False)
__button.setEnabled(False)
__button.hide()
self.mdl.setPropertyMode(tag, 'predefined_law')
# Atmospheric Flows
if mdl_atmo != 'off':
if tag == 'density':
__model.disableItem(str_model='constant')
__model.disableItem(str_model='predefined_law')
__model.setItem(str_model='predefined_law')
__combo.setEnabled(False)
__button.setEnabled(False)
__button.hide()
# VoF
if mdl_hgn != 'off':
if tag == 'molecular_viscosity' or tag == 'density':
__model.disableItem(str_model='constant')
__model.disableItem(str_model='predefined_law')
__model.setItem(str_model='predefined_law')
__combo.setEnabled(False)
__button.setEnabled(False)
__button.hide()
# Compressible Flows
if mdl_comp != 'off':
self.groupBoxViscv0.show()
if tag == 'density':
__model.setItem(str_model='predefined_law')
__combo.setEnabled(False)
__button.setEnabled(False)
__button.hide()
self.mdl.setPropertyMode(tag, 'predefined_law')
__line.setEnabled(True)
if tag == 'specific_heat':
__model.setItem(str_model='constant')
__combo.setEnabled(False)
__button.setEnabled(False)
__button.hide()
self.mdl.setPropertyMode(tag, 'constant')
self.groupBoxCp.setTitle('Isobaric specific heat')
if tag == 'volume_viscosity':
__combo.setEnabled(True)
c = self.mdl.getPropertyMode(tag)
if c == 'user_law':
__button.setEnabled(True)
else:
__button.setEnabled(False)
else:
if tag == 'specific_heat':
self.groupBoxCp.setTitle('Specific heat')
if mdl_hgn != 'off':
self.lineEditRho1.setText(str(self.mdl.getVofValueDensity(0)))
self.lineEditRho2.setText(str(self.mdl.getVofValueDensity(1)))
self.lineEditMu1.setText(str(self.mdl.getVofValueViscosity(0)))
self.lineEditMu2.setText(str(self.mdl.getVofValueViscosity(1)))
def updateTypeChoice(self, old_choice):
"""
add/suppress thermo tables for each properties
"""
for tag, symbol in self.mdl.lst:
__model = getattr(self, "model" + symbol)
if self.mdl.getMaterials() == "user_material":
__model.disableItem(str_model='thermal_law')
else:
__model.enableItem(str_model='thermal_law')
if old_choice == "user_material":
self.mdl.setPropertyMode(tag, 'thermal_law')
self.__changeChoice(str("material law"), symbol, tag)
c = self.mdl.getPropertyMode(tag)
__model.setItem(str_model=c)
def updateMethod(self):
"""
update method list with material choice
"""
for nb in range(len(self.modelMethod.getItems())):
self.modelMethod.delItem(0)
material = self.mdl.getMaterials()
method = self.mdl.getMethod()
have_method = False
methods = self.mdl.getLibPropertyMethods(material)
for m in methods:
if method == m[0]:
have_method = True
self.modelMethod.addItem(self.tr(m[0]), m[0], not m[1])
# current method not in list, add it with warning
if not have_method:
self.modelMethod.addItem(self.tr(method), method, True)
# update comboBoxMethod
self.modelMethod.setItem(str_model=method)
self.mdl.setMethod(method)
self.updateReference(material, method)
def updateReference(self, material, method):
"""
update method list with material choice
"""
self.comboBoxReference.hide()
self.labelReference.hide()
for nb in range(len(self.modelReference.getItems())):
self.modelReference.delItem(0)
reference = self.mdl.getReference()
references = self.mdl.getAvailReferences(material, method)
have_reference = False
if references:
for r in references:
if reference == r:
have_reference = True
self.modelReference.addItem(self.tr(r), r)
self.comboBoxReference.show()
self.labelReference.show()
if not have_reference:
reference = references[0]
self.modelReference.setItem(str_model=reference)
self.mdl.setReference(reference)
@pyqtSlot(str)
def slotPressure(self, text):
"""
Input PRESS.
"""
if self.lineEditP0.validator().state == QValidator.Acceptable:
p = from_qvariant(text, float)
self.mdl.setPressure(p)
@pyqtSlot(str)
def slotTemperature(self, text):
"""
Input TEMPERATURE.
"""
if self.lineEditT0.validator().state == QValidator.Acceptable:
t = from_qvariant(text, float)
self.mdl.setTemperature(t)
@pyqtSlot(str)
def slotTempOxydant(self, text):
"""
Input oxydant TEMPERATURE.
"""
if self.lineEditOxydant.validator().state == QValidator.Acceptable:
t = from_qvariant(text, float)
self.mdl.setTempOxydant(t)
@pyqtSlot(str)
def slotTempFuel(self, text):
"""
Input fuel TEMPERATURE.
"""
if self.lineEditFuel.validator().state == QValidator.Acceptable:
t = from_qvariant(text, float)
self.mdl.setTempFuel(t)
@pyqtSlot(str)
def slotMassemol(self, text):
"""
Input Mass molar.
"""
if self.lineEditMassMolar.validator().state == QValidator.Acceptable:
m = from_qvariant(text, float)
self.mdl.setMassemol(m)
@pyqtSlot(str)
def slotMaterial(self, text):
"""
Method to call 'setMaterial'
"""
choice = self.modelMaterial.dicoV2M[str(text)]
old_choice = self.mdl.getMaterials()
self.mdl.setMaterials(choice)
self.updateMethod()
self.updateTypeChoice(old_choice)
@pyqtSlot(str)
def slotMethod(self, text):
"""
Method to call 'setMethod'
"""
choice = self.modelMethod.dicoV2M[str(text)]
self.mdl.setMethod(choice)
self.updateReference(self.mdl.getMaterials(), choice)
@pyqtSlot(str)
def slotReference(self, text):
"""
Method to call 'setReference'
"""
choice = self.modelReference.dicoV2M[str(text)]
self.mdl.setReference(choice)
@pyqtSlot(str)
def slotStateRho(self, text):
"""
Method to call 'getState' with correct arguements for 'rho'
"""
self.__changeChoice(str(text), 'Rho', 'density')
@pyqtSlot(str)
def slotStateMu(self, text):
"""
Method to call 'getState' with correct arguements for 'Mu'
"""
self.__changeChoice(str(text), 'Mu', 'molecular_viscosity')
@pyqtSlot(str)
def slotStateCp(self, text):
"""
Method to call 'getState' with correct arguements for 'Cp'
"""
self.__changeChoice(str(text), 'Cp', 'specific_heat')
@pyqtSlot(str)
def slotStateViscv0(self, text):
"""
Method to call 'getState' with correct arguements for 'Viscv0'
"""
self.__changeChoice(str(text), 'Viscv0', 'volume_viscosity')
@pyqtSlot(str)
def slotStateAl(self, text):
"""
Method to call 'getState' with correct arguements for 'Al'
"""
self.__changeChoice(str(text), 'Al', 'thermal_conductivity')
@pyqtSlot(str)
def slotStateDiff(self, text):
"""
Method to set diffusion choice for the coefficient
"""
choice = self.modelDiff.dicoV2M[str(text)]
log.debug("slotStateDiff -> %s" % (text))
if choice != 'user_law':
self.pushButtonDiff.setEnabled(False)
self.pushButtonDiff.setStyleSheet("background-color: None")
self.pushButtonDiff.hide()
else:
self.pushButtonDiff.show()
self.pushButtonDiff.setEnabled(True)
name = self.mdl.m_sca.getScalarDiffusivityName(self.scalar)
exp = self.mdl.m_sca.getDiffFormula(self.scalar)
if exp:
self.pushButtonDiff.setStyleSheet("background-color: green")
self.pushButtonDiff.setToolTip(exp)
else:
self.pushButtonDiff.setStyleSheet("background-color: red")
self.mdl.m_sca.setScalarDiffusivityChoice(self.scalar, choice)
@pyqtSlot(str)
def slotNameDiff(self, text):
"""
Method to set the variance scalar choosed
"""
choice = self.modelNameDiff.dicoV2M[str(text)]
log.debug("slotStateDiff -> %s" % (text))
self.scalar = str(text)
self.lineEditDiff.setText(
str(self.mdl.m_sca.getScalarDiffusivityInitialValue(self.scalar)))
mdl = self.mdl.m_sca.getScalarDiffusivityChoice(self.scalar)
self.modelDiff.setItem(str_model=mdl)
if mdl != 'user_law':
self.pushButtonDiff.hide()
self.pushButtonDiff.setEnabled(False)
self.pushButtonDiff.setStyleSheet("background-color: None")
else:
self.pushButtonDiff.show()
self.pushButtonDiff.setEnabled(True)
name = self.mdl.m_sca.getScalarDiffusivityName(self.scalar)
exp = self.mdl.m_sca.getDiffFormula(self.scalar)
if exp:
self.pushButtonDiff.setStyleSheet("background-color: green")
self.pushButtonDiff.setToolTip(exp)
else:
self.pushButtonDiff.setStyleSheet("background-color: red")
def __changeChoice(self, text, sym, tag):
"""
Input variable state
"""
__model = getattr(self, "model" + sym)
__line = getattr(self, "lineEdit" + sym)
__combo = getattr(self, "comboBox" + sym)
__label = getattr(self, "label" + sym)
__button = getattr(self, "pushButton" + sym)
__labelu = getattr(self, "labelUnit" + sym)
__labelv = getattr(self, "labelVar" + sym)
choice = __model.dicoV2M[text]
log.debug("__changeChoice -> %s, %s" % (text, choice))
if choice != 'user_law':
__button.hide()
__button.setEnabled(False)
__button.setStyleSheet("background-color: None")
else:
__button.setEnabled(True)
__button.show()
exp = None
if sym == "Rho":
exp = self.mdl.getFormula('density')
elif sym == "Mu":
exp = self.mdl.getFormula('molecular_viscosity')
elif sym == "Cp":
exp = self.mdl.getFormula('specific_heat')
elif sym == "Viscv0":
exp = self.mdl.getFormula('volume_viscosity')
elif sym == "Al":
exp = self.mdl.getFormula('thermal_conductivity')
elif sym == "Diff":
name = self.mdl.m_sca.getScalarDiffusivityName(self.scalar)
exp = self.mdl.m_sca.getDiffFormula(self.scalar)
if exp:
__button.setStyleSheet("background-color: green")
__button.setToolTip(exp)
else:
__button.setStyleSheet("background-color: red")
if choice == 'thermal_law':
__line.hide()
__label.hide()
__labelu.hide()
__labelv.hide()
else:
__line.show()
__label.show()
__labelu.show()
__labelv.show()
self.mdl.setPropertyMode(tag, choice)
@pyqtSlot(str)
def slotRho(self, text):
"""
Update the density
"""
if self.lineEditRho.validator().state == QValidator.Acceptable:
rho = from_qvariant(text, float)
self.mdl.setInitialValueDensity(rho)
@pyqtSlot(str)
def slotRho1(self, text):
"""
Update the density of fluid 1 for VoF module
"""
if self.lineEditRho1.validator().state == QValidator.Acceptable:
rho = from_qvariant(text, float)
self.mdl.setVofValueDensity(0, rho)
@pyqtSlot(str)
def slotRho2(self, text):
"""
Update the density of fluid 2 for VoF module
"""
if self.lineEditRho2.validator().state == QValidator.Acceptable:
rho = from_qvariant(text, float)
self.mdl.setVofValueDensity(1, rho)
@pyqtSlot(str)
def slotMu(self, text):
"""
Update the molecular viscosity
"""
if self.lineEditMu.validator().state == QValidator.Acceptable:
mu = from_qvariant(text, float)
self.mdl.setInitialValueViscosity(mu)
@pyqtSlot(str)
def slotMu1(self, text):
"""
Update the molecular viscosity
"""
if self.lineEditMu1.validator().state == QValidator.Acceptable:
mu = from_qvariant(text, float)
self.mdl.setVofValueViscosity(0, mu)
@pyqtSlot(str)
def slotMu2(self, text):
"""
Update the molecular viscosity
"""
if self.lineEditMu2.validator().state == QValidator.Acceptable:
mu = from_qvariant(text, float)
self.mdl.setVofValueViscosity(1, mu)
@pyqtSlot(str)
def slotCp(self, text):
"""
Update the specific heat
"""
if self.lineEditCp.validator().state == QValidator.Acceptable:
cp = from_qvariant(text, float)
self.mdl.setInitialValueHeat(cp)
@pyqtSlot(str)
def slotViscv0(self, text):
"""
Update the volumic viscosity
"""
if self.lineEditViscv0.validator().state == QValidator.Acceptable:
viscv0 = from_qvariant(text, float)
self.mdl.setInitialValueVolumeViscosity(viscv0)
@pyqtSlot(str)
def slotAl(self, text):
"""
Update the thermal conductivity
"""
if self.lineEditAl.validator().state == QValidator.Acceptable:
al = from_qvariant(text, float)
self.mdl.setInitialValueCond(al)
@pyqtSlot(str)
def slotDiftl0(self, text):
"""
Update the thermal conductivity
"""
if self.lineEditDiftl0.validator().state == QValidator.Acceptable:
diftl0 = from_qvariant(text, float)
self.mdl.setInitialValueDyn(diftl0)
@pyqtSlot(str)
def slotDiff(self, text):
"""
Update the thermal conductivity
"""
if self.lineEditDiff.validator().state == QValidator.Acceptable:
diff = from_qvariant(text, float)
self.mdl.m_sca.setScalarDiffusivityInitialValue(self.scalar, diff)
@pyqtSlot()
def slotFormulaRho(self):
"""
User formula for density
"""
exp, req, sca, symbols_rho = self.mdl.getFormulaRhoComponents()
self.m_th = ThermalScalarModel(self.case)
s = self.m_th.getThermalScalarName()
mdl = self.m_th.getThermalScalarModel()
if mdl == "off":
exa = FluidCharacteristicsView.density_wo
elif mdl == "temperature_celsius":
TempInContext = "(" + s + " + 273.15)"
exa = FluidCharacteristicsView.density.replace(
"temperature", TempInContext)
elif mdl == "enthalpy":
exa = FluidCharacteristicsView.density_h
else:
exa = FluidCharacteristicsView.density
dialog = QMegEditorView(parent=self,
function_type='vol',
zone_name='all_cells',
variable_name='density',
expression=exp,
required=req,
symbols=symbols_rho,
known_fields=sca,
examples=exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaRho -> %s" % str(result))
self.mdl.setFormula('density', str(result))
self.pushButtonRho.setToolTip(result)
self.pushButtonRho.setStyleSheet("background-color: green")
@pyqtSlot()
def slotFormulaMu(self):
"""
User formula for molecular viscosity
"""
exp, req, sca, symbols_mu = self.mdl.getFormulaMuComponents()
self.m_th = ThermalScalarModel(self.case)
s = self.m_th.getThermalScalarName()
mdl = self.m_th.getThermalScalarModel()
if mdl == "off":
exa = FluidCharacteristicsView.molecular_viscosity_wo
elif mdl == "temperature_celsius":
TempInContext = "(" + s + " + 273.15)"
exa = FluidCharacteristicsView.molecular_viscosity.replace(
"temperature", TempInContext)
elif mdl == "enthalpy":
exa = FluidCharacteristicsView.molecular_viscosity_h
else:
exa = FluidCharacteristicsView.molecular_viscosity
dialog = QMegEditorView(parent=self,
function_type='vol',
zone_name='all_cells',
variable_name='molecular_viscosity',
expression=exp,
required=req,
symbols=symbols_mu,
known_fields=sca,
examples=exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaMu -> %s" % str(result))
self.mdl.setFormula('molecular_viscosity', str(result))
self.pushButtonMu.setToolTip(result)
self.pushButtonMu.setStyleSheet("background-color: green")
@pyqtSlot()
def slotFormulaCp(self):
"""
User formula for specific heat
"""
exp, req, sca, symbols_cp = self.mdl.getFormulaCpComponents()
exa = FluidCharacteristicsView.specific_heat
dialog = QMegEditorView(parent=self,
function_type='vol',
zone_name='all_cells',
variable_name='specific_heat',
expression=exp,
required=req,
symbols=symbols_cp,
known_fields=sca,
examples=exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaRho -> %s" % str(result))
self.mdl.setFormula('specific_heat', str(result))
self.pushButtonCp.setToolTip(result)
self.pushButtonCp.setStyleSheet("background-color: green")
@pyqtSlot()
def slotFormulaViscv0(self):
"""
User formula for volumic viscosity
"""
exp, req, sca, symbols_viscv0 = self.mdl.getFormulaViscv0Components()
exa = FluidCharacteristicsView.volume_viscosity
dialog = QMegEditorView(parent=self,
function_type='vol',
zone_name='all_cells',
variable_name='volume_viscosity',
expression=exp,
required=req,
symbols=symbols_viscv0,
known_fields=sca,
examples=exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaViscv0 -> %s" % str(result))
self.mdl.setFormula('volume_viscosity', str(result))
self.pushButtonViscv0.setToolTip(result)
self.pushButtonViscv0.setStyleSheet("background-color: green")
@pyqtSlot()
def slotFormulaAl(self):
"""
User formula for thermal conductivity
"""
exp, req, sca, symbols_al = self.mdl.getFormulaAlComponents()
self.m_th = ThermalScalarModel(self.case)
s = self.m_th.getThermalScalarName()
mdl = self.m_th.getThermalScalarModel()
if mdl == "temperature_celsius":
TempInContext = "(" + s + " + 273.15)"
exa = FluidCharacteristicsView.thermal_conductivity.replace(
"temperature", TempInContext)
elif mdl == "enthalpy":
exa = FluidCharacteristicsView.thermal_conductivity_h
else:
exa = FluidCharacteristicsView.thermal_conductivity
dialog = QMegEditorView(parent=self,
function_type='vol',
zone_name='all_cells',
variable_name='thermal_conductivity',
expression=exp,
required=req,
symbols=symbols_al,
known_fields=sca,
examples=exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaAl -> %s" % str(result))
self.mdl.setFormula('thermal_conductivity', str(result))
self.pushButtonAl.setToolTip(result)
self.pushButtonAl.setStyleSheet("background-color: green")
@pyqtSlot()
def slotFormulaDiff(self):
"""
User formula for the diffusion coefficient
"""
exp, req, sca, sym = self.mdl.getFormulaDiffComponents(self.scalar)
exa = ''
dname = self.mdl.m_sca.getScalarDiffusivityName(self.scalar)
dialog = QMegEditorView(parent=self,
function_type='vol',
zone_name='all_cells',
variable_name=dname,
expression=exp,
required=req,
symbols=sym,
known_fields=sca,
examples=exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaDiff -> %s" % str(result))
self.mdl.m_sca.setDiffFormula(self.scalar, str(result))
self.pushButtonDiff.setToolTip(result)
self.pushButtonDiff.setStyleSheet("background-color: green")
class FluidCharacteristicsView(QWidget, Ui_FluidCharacteristicsForm):
"""
Class to open Molecular Properties Page.
"""
density = """# Density of air
density = 1.293*(273.15 / temperature);
# density for mixtures of gases
#
# Y1 -> mass fraction of component 1
# Y2 -> mass fraction of component 2
rho1 = 1.25051;
rho2 = 1.7832;
A = (Y1 / rho1) + (Y2 /rho2);
density = 1.0 / A;
"""
density_h = """# Density
density = enthalpy / 1040. * 1.29;
# density for mixtures of gases
#
# Y1 -> mass fraction of component 1
# Y2 -> mass fraction of component 2
rho1 = 1.25051;
rho2 = 1.7832;
A = (Y1 / rho1) + (Y2 /rho2);
density = 1.0 / A;
"""
molecular_viscosity="""# Sutherland's Formula
# Gas Cst T0 mu0
# air 120 291.15 18.27e-6
# nitrogen 111 300.55 17.81e-6
# oxygen 127 292.25 20.18e-6
# carbon dioxide 240 293.15 14.8e-6
# carbon monoxide 118 288.15 17.2e-6
# hydrogen 72 293.85 8.76e-6
# ammonia 370 293.15 9.82e-6
# sulfur dioxide 416 293.65 12.54e-6
# helium 79.4 273 19e-6
CST = 120;
T0 = 291.15;
mu_ref = 18.27e-6;
if ( temperature > 0 && temperature < 555) {
molecular_viscosity = mu_ref * ((T0+CST) / (temperature+CST)) * (temperature/T0)^(3./2.);
} else {
molecular_viscosity = -999.0;
}
"""
molecular_viscosity_h="""
CST = 120;
T0 = 291.15;
mu_ref = 18.27e-6;
temperature = enthalpy / 1040.;
if ( enthalpy > 0) {
molecular_viscosity = mu_ref * (T0+CST / temperature+CST) * (temperature/T0)^(3./2.);
} else {
molecular_viscosity = -999.0;
}
"""
specific_heat="""# specific heat for mixtures of gases
#
# Y1 -> mass fraction of component 1
# Y2 -> mass fraction of component 2
Cp1 = 520.3;
Cp2 = 1040.0;
specific_heat = Y1 * Cp1 + Y2 *Cp2;
"""
volume_viscosity="""# volume_viscosity
"""
thermal_conductivity="""# oxygen
thermal_conductivity = 6.2e-5 * temperature + 8.1e-3;
# nitrogen
thermal_conductivity = 6.784141e-5 * temperature + 5.564317e-3;
# hydrogen
thermal_conductivity = 4.431e-4 * temperature + 5.334e-2;
"""
thermal_conductivity_h="""
temperature = enthalpy / 1040.;
thermal_conductivity = 6.2e-5 * temperature + 8.1e-3;
"""
def __init__(self, parent, case):
"""
Constructor
"""
QWidget.__init__(self, parent)
Ui_FluidCharacteristicsForm.__init__(self)
self.setupUi(self)
self.case = case
self.case.undoStopGlobal()
self.mdl = FluidCharacteristicsModel(self.case)
if EOS == 1:
self.ava = eosAva.EosAvailable()
import cs_config
cfg = cs_config.config()
self.freesteam = 0
if cfg.libs['freesteam'].have != "no":
self.freesteam = 1
if CompressibleModel(self.case).getCompressibleModel() != 'off':
self.lst = [('density', 'Rho'),
('molecular_viscosity', 'Mu'),
('specific_heat', 'Cp'),
('thermal_conductivity', 'Al'),
('volume_viscosity', 'Viscv0'),
('dynamic_diffusion', 'Diftl0')]
elif CoalCombustionModel(self.case).getCoalCombustionModel() != 'off' or \
GasCombustionModel(self.case).getGasCombustionModel() != 'off':
self.lst = [('density', 'Rho'),
('molecular_viscosity', 'Mu'),
('specific_heat', 'Cp'),
('dynamic_diffusion', 'Diftl0')]
else:
self.lst = [('density', 'Rho'),
('molecular_viscosity', 'Mu'),
('specific_heat', 'Cp'),
('thermal_conductivity', 'Al')]
self.list_scalars = []
self.m_th = ThermalScalarModel(self.case)
s = self.m_th.getThermalScalarName()
mdl = self.m_th.getThermalScalarModel()
if mdl == "temperature_celsius":
self.list_scalars.append((s, self.tr("Thermal scalar: temperature (C)")))
elif mdl == "temperature_kelvin":
self.list_scalars.append((s, self.tr("Thermal scalar: temperature (K)")))
elif mdl != "off":
self.list_scalars.append((s, self.tr("Thermal scalar")))
self.m_sca = DefineUserScalarsModel(self.case)
for s in self.m_sca.getUserScalarNameList():
self.list_scalars.append((s, self.tr("Additional scalar")))
# Particular Widget initialization taking into account of "Calculation Features"
mdl_atmo, mdl_joule, mdl_thermal, mdl_gas, mdl_coal, mdl_comp = self.mdl.getThermoPhysicalModel()
# Combo models
self.modelRho = ComboModel(self.comboBoxRho, 3, 1)
self.modelMu = ComboModel(self.comboBoxMu, 3, 1)
self.modelCp = ComboModel(self.comboBoxCp, 3, 1)
self.modelAl = ComboModel(self.comboBoxAl, 3, 1)
self.modelDiff = ComboModel(self.comboBoxDiff, 2, 1)
self.modelNameDiff = ComboModel(self.comboBoxNameDiff,1,1)
self.modelViscv0 = ComboModel(self.comboBoxViscv0, 3, 1)
self.modelDiftl0 = ComboModel(self.comboBoxDiftl0, 3, 1)
self.modelMaterial = ComboModel(self.comboBoxMaterial, 1, 1)
self.modelMethod = ComboModel(self.comboBoxMethod, 1, 1)
self.modelPhas = ComboModel(self.comboBoxPhas, 2, 1)
self.modelRho.addItem(self.tr('constant'), 'constant')
self.modelRho.addItem(self.tr('variable'), 'variable')
self.modelRho.addItem(self.tr('thermal law'), 'thermal_law')
if mdl_atmo != 'off':
self.modelRho.addItem(self.tr('defined in atphyv'), 'variable')
elif mdl_joule == 'arc':
self.modelRho.addItem(self.tr('defined in elphyv'), 'variable')
self.modelMu.addItem(self.tr('constant'), 'constant')
self.modelMu.addItem(self.tr('variable'), 'variable')
self.modelMu.addItem(self.tr('thermal law'), 'thermal_law')
if mdl_joule == 'arc':
self.modelMu.addItem(self.tr('defined in elphyv'), 'variable')
self.modelCp.addItem(self.tr('constant'), 'constant')
self.modelCp.addItem(self.tr('variable'), 'variable')
self.modelCp.addItem(self.tr('thermal law'), 'thermal_law')
if mdl_joule == 'arc':
self.modelCp.addItem(self.tr('defined in elphyv'), 'variable')
self.modelAl.addItem(self.tr('constant'), 'constant')
self.modelAl.addItem(self.tr('variable'), 'variable')
self.modelAl.addItem(self.tr('thermal law'), 'thermal_law')
if mdl_joule == 'arc':
self.modelAl.addItem(self.tr('defined in elphyv'), 'variable')
self.modelDiff.addItem(self.tr('constant'), 'constant')
self.modelDiff.addItem(self.tr('variable'), 'variable')
self.modelViscv0.addItem(self.tr('constant'), 'constant')
self.modelViscv0.addItem(self.tr('variable'), 'variable')
self.modelViscv0.addItem(self.tr('thermal law'), 'thermal_law')
self.modelDiftl0.addItem(self.tr('constant'), 'constant')
self.modelDiftl0.addItem(self.tr('variable'), 'variable')
self.modelDiftl0.addItem(self.tr('thermal law'), 'thermal_law')
self.modelPhas.addItem(self.tr('liquid'), 'liquid')
self.modelPhas.addItem(self.tr('gas'), 'gas')
if (self.freesteam == 0 and EOS == 0) or \
self.m_th.getThermalScalarModel() == "off" or \
mdl_joule != 'off' or mdl_comp != 'off':
self.groupBoxTableChoice.hide()
else:
self.groupBoxTableChoice.show()
self.lineEditReference.setEnabled(False)
# suppress perfect gas
self.modelMaterial.addItem(self.tr('user material'), 'user_material')
tmp = ["Argon", "Nitrogen", "Hydrogen", "Oxygen", "Helium", "Air"]
if EOS == 1:
fls = self.ava.whichFluids()
for fli in fls:
if fli not in tmp:
tmp.append(fli)
self.modelMaterial.addItem(self.tr(fli), fli)
if self.freesteam == 1 and EOS == 0:
self.modelMaterial.addItem(self.tr('Water'), 'Water')
material = self.mdl.getMaterials()
self.modelMaterial.setItem(str_model=material)
self.updateMethod()
self.scalar = ""
scalar_list = self.m_sca.getUserScalarNameList()
for s in self.m_sca.getScalarsVarianceList():
if s in scalar_list: scalar_list.remove(s)
if scalar_list != []:
self.scalar = scalar_list[0]
for scalar in scalar_list:
self.modelNameDiff.addItem(scalar)
# Connections
self.connect(self.comboBoxRho, SIGNAL("activated(const QString&)"), self.slotStateRho)
self.connect(self.comboBoxMu, SIGNAL("activated(const QString&)"), self.slotStateMu)
self.connect(self.comboBoxCp, SIGNAL("activated(const QString&)"), self.slotStateCp)
self.connect(self.comboBoxAl, SIGNAL("activated(const QString&)"), self.slotStateAl)
self.connect(self.comboBoxDiff, SIGNAL("activated(const QString&)"), self.slotStateDiff)
self.connect(self.comboBoxNameDiff, SIGNAL("activated(const QString&)"), self.slotNameDiff)
self.connect(self.comboBoxViscv0, SIGNAL("activated(const QString&)"), self.slotStateViscv0)
self.connect(self.comboBoxMaterial, SIGNAL("activated(const QString&)"), self.slotMaterial)
self.connect(self.comboBoxMethod, SIGNAL("activated(const QString&)"), self.slotMethod)
self.connect(self.comboBoxPhas, SIGNAL("activated(const QString&)"), self.slotPhas)
self.connect(self.lineEditRho, SIGNAL("textChanged(const QString &)"), self.slotRho)
self.connect(self.lineEditMu, SIGNAL("textChanged(const QString &)"), self.slotMu)
self.connect(self.lineEditCp, SIGNAL("textChanged(const QString &)"), self.slotCp)
self.connect(self.lineEditAl, SIGNAL("textChanged(const QString &)"), self.slotAl)
self.connect(self.lineEditDiff, SIGNAL("textChanged(const QString &)"), self.slotDiff)
self.connect(self.lineEditDiftl0, SIGNAL("textChanged(const QString &)"), self.slotDiftl0)
self.connect(self.lineEditViscv0, SIGNAL("textChanged(const QString &)"), self.slotViscv0)
self.connect(self.pushButtonRho, SIGNAL("clicked()"), self.slotFormulaRho)
self.connect(self.pushButtonMu, SIGNAL("clicked()"), self.slotFormulaMu)
self.connect(self.pushButtonCp, SIGNAL("clicked()"), self.slotFormulaCp)
self.connect(self.pushButtonAl, SIGNAL("clicked()"), self.slotFormulaAl)
self.connect(self.pushButtonDiff, SIGNAL("clicked()"), self.slotFormulaDiff)
self.connect(self.pushButtonViscv0, SIGNAL("clicked()"), self.slotFormulaViscv0)
# Validators
validatorRho = DoubleValidator(self.lineEditRho, min = 0.0)
validatorMu = DoubleValidator(self.lineEditMu, min = 0.0)
validatorCp = DoubleValidator(self.lineEditCp, min = 0.0)
validatorAl = DoubleValidator(self.lineEditAl, min = 0.0)
validatorDiff = DoubleValidator(self.lineEditDiff, min = 0.0)
validatorViscv0 = DoubleValidator(self.lineEditViscv0, min = 0.0)
validatorDiftl0 = DoubleValidator(self.lineEditDiftl0, min = 0.0)
validatorRho.setExclusiveMin(True)
validatorMu.setExclusiveMin(True)
validatorCp.setExclusiveMin(True)
validatorAl.setExclusiveMin(True)
validatorDiff.setExclusiveMin(True)
validatorDiftl0.setExclusiveMin(True)
self.lineEditRho.setValidator(validatorRho)
self.lineEditMu.setValidator(validatorMu)
self.lineEditCp.setValidator(validatorCp)
self.lineEditAl.setValidator(validatorAl)
self.lineEditDiff.setValidator(validatorDiff)
self.lineEditViscv0.setValidator(validatorViscv0)
self.lineEditDiftl0.setValidator(validatorDiftl0)
if scalar_list == []:
self.groupBoxDiff.hide()
else :
self.groupBoxDiff.show()
self.lineEditDiff.setText(str(self.m_sca.getScalarDiffusivityInitialValue(self.scalar)))
diff_choice = self.m_sca.getScalarDiffusivityChoice(self.scalar)
self.modelDiff.setItem(str_model=diff_choice)
self.modelNameDiff.setItem(str_model=str(self.scalar))
if diff_choice != 'variable':
self.pushButtonDiff.setEnabled(False)
setGreenColor(self.pushButtonDiff, False)
else:
self.pushButtonDiff.setEnabled(True)
setGreenColor(self.pushButtonDiff, True)
#compressible
self.groupBoxViscv0.hide()
# combustion
self.groupBoxDiftl0.hide()
# Standard Widget initialization
for tag, symbol in self.lst:
__model = getattr(self, "model" + symbol)
__line = getattr(self, "lineEdit" + symbol)
__button = getattr(self, "pushButton" + symbol)
__label = getattr(self, "label" + symbol)
__labelu = getattr(self, "labelUnit" + symbol)
if tag != 'dynamic_diffusion':
__labelv = getattr(self, "labelVar" + symbol)
c = self.mdl.getPropertyMode(tag)
__model.setItem(str_model=c)
if c == 'variable':
__button.setEnabled(True)
__label.setText(self.tr("Reference value"))
else:
__button.setEnabled(False)
__label.setText(self.tr("Reference value"))
if c == 'thermal_law':
__line.hide()
__label.hide()
__labelu.hide()
__labelv.hide()
else:
__line.show()
__label.show()
__labelu.show()
__labelv.show()
if self.mdl.getMaterials() == "user_material":
__model.disableItem(str_model='thermal_law')
else:
__model.enableItem(str_model='thermal_law')
else:
__label.setText(self.tr("Reference value"))
self.mdl.getInitialValue(tag)
__line.setText(str(self.mdl.getInitialValue(tag)))
# no 'thermal_conductivity' if not Joule and not Thermal scalar and not
if mdl_joule == 'off' and mdl_thermal == 'off' and mdl_atmo == 'off' and\
CompressibleModel(self.case).getCompressibleModel() == 'off':
self.groupBoxAl.hide()
if mdl_gas != 'off' or mdl_coal != 'off':
self.groupBoxDiftl0.show()
for tag, symbol in self.lst:
__model = getattr(self, "model" + symbol)
__line = getattr(self, "lineEdit" + symbol)
__button = getattr(self, "pushButton" + symbol)
__label = getattr(self, "label" + symbol)
__combo = getattr(self, "comboBox" + symbol)
# Gas or coal combustion
if mdl_gas != 'off' or mdl_coal != 'off':
if tag == 'density':
__model.setItem(str_model='variable')
__combo.setEnabled(False)
__button.setEnabled(False)
self.mdl.setPropertyMode(tag, 'variable')
__label.setText(self.tr("Calculation by\n perfect gas law"))
__line.setText(str(""))
__line.setEnabled(False)
elif tag == 'dynamic_diffusion':
__model.setItem(str_model='variable')
__combo.setEnabled(False)
__button.setEnabled(False)
else:
__model.setItem(str_model='constant')
self.mdl.setPropertyMode(tag, 'constant')
# Joule
if mdl_joule == 'arc':
__model.disableItem(str_model='constant')
__model.disableItem(str_model='variable')
__model.setItem(str_model='variable')
__combo.setEnabled(False)
__button.setEnabled(False)
self.mdl.setPropertyMode(tag, 'variable')
if mdl_joule == 'joule':
__model.setItem(str_model='variable')
__model.disableItem(str_model='constant')
self.mdl.setPropertyMode(tag, 'variable')
# Atmospheric Flows
if mdl_atmo != 'off':
if tag == 'density':
__model.disableItem(str_model='constant')
__model.disableItem(str_model='variable')
__model.setItem(str_model='variable')
__combo.setEnabled(False)
__button.setEnabled(False)
# Compressible Flows
if mdl_comp != 'off':
if tag == 'density':
__model.setItem(str_model='variable')
__combo.setEnabled(False)
__button.setEnabled(False)
__combo.hide()
__button.hide()
self.mdl.setPropertyMode(tag, 'variable')
__line.setEnabled(True)
self.groupBoxViscv0.hide()
if tag == 'specific_heat':
__model.setItem(str_model='constant')
__combo.setEnabled(False)
__button.setEnabled(False)
self.mdl.setPropertyMode(tag, 'constant')
self.groupBoxCp.setTitle('Isobaric specific heat')
if tag == 'volume_viscosity':
__combo.setEnabled(True)
c = self.mdl.getPropertyMode(tag)
if c == 'variable':
__button.setEnabled(True)
else:
__button.setEnabled(False)
self.groupBoxViscv0.show()
else:
if tag == 'specific_heat':
self.groupBoxCp.setTitle('Specific heat')
self.case.undoStartGlobal()
def updateTypeChoice(self):
"""
add/suppress thermo tables for each proprties
"""
for tag, symbol in self.lst:
__model = getattr(self, "model" + symbol)
if self.mdl.getMaterials() == "user_material":
__model.disableItem(str_model='thermal_law')
else:
__model.enableItem(str_model='thermal_law')
c = self.mdl.getPropertyMode(tag)
__model.setItem(str_model=c)
def updateMethod(self):
"""
update method list with material choice
"""
for nb in range(len(self.modelMethod.getItems())):
self.modelMethod.delItem(0)
self.comboBoxPhas.hide()
self.labelPhas.hide()
if self.mdl.getMaterials() == "user_material":
self.modelMethod.addItem(self.tr('user properties'), 'user_properties')
else :
if EOS == 1:
material = self.mdl.getMaterials()
self.ava.setMethods(material)
fls = self.ava.whichMethods()
for fli in fls:
self.modelMethod.addItem(self.tr(fli),fli)
if self.mdl.getMethod() != "freesteam":
self.comboBoxPhas.show()
self.labelPhas.show()
if self.freesteam == 1 and self.mdl.getMaterials() == "Water":
self.modelMethod.addItem(self.tr("freesteam"), "freesteam")
# update comboBoxMethod
method = self.mdl.getMethod()
self.modelMethod.setItem(str_model=method)
self.updateReference()
def updateReference(self):
"""
update Reference with material, method and field nature choice
"""
# update lineEditReference
self.lineEditReference.setText(self.mdl.getReference())
@pyqtSignature("const QString &")
def slotMaterial(self, text):
"""
Method to call 'setMaterial'
"""
choice = self.modelMaterial.dicoV2M[str(text)]
self.mdl.setMaterials(choice)
self.updateMethod()
self.updateTypeChoice()
@pyqtSignature("const QString &")
def slotPhas(self, text):
"""
Method to call 'setFieldNature'
"""
choice = self.modelPhas.dicoV2M[str(text)]
self.mdl.setFieldNature(choice)
self.updateReference()
@pyqtSignature("const QString &")
def slotMethod(self, text):
"""
Method to call 'setMethod'
"""
choice = self.modelMethod.dicoV2M[str(text)]
self.mdl.setMethod(choice)
self.comboBoxPhas.hide()
self.labelPhas.hide()
if self.mdl.getMaterials() != "user_material" and \
self.mdl.getMethod() != "freesteam":
self.comboBoxPhas.show()
self.labelPhas.show()
self.updateReference()
@pyqtSignature("const QString &")
def slotStateRho(self, text):
"""
Method to call 'getState' with correct arguements for 'rho'
"""
self.__changeChoice(str(text), 'Rho', 'density')
@pyqtSignature("const QString &")
def slotStateMu(self, text):
"""
Method to call 'getState' with correct arguements for 'Mu'
"""
self.__changeChoice(str(text), 'Mu', 'molecular_viscosity')
@pyqtSignature("const QString &")
def slotStateCp(self, text):
"""
Method to call 'getState' with correct arguements for 'Cp'
"""
self.__changeChoice(str(text), 'Cp', 'specific_heat')
@pyqtSignature("const QString &")
def slotStateViscv0(self, text):
"""
Method to call 'getState' with correct arguements for 'Viscv0'
"""
self.__changeChoice(str(text), 'Viscv0', 'volume_viscosity')
@pyqtSignature("const QString &")
def slotStateAl(self, text):
"""
Method to call 'getState' with correct arguements for 'Al'
"""
self.__changeChoice(str(text), 'Al', 'thermal_conductivity')
@pyqtSignature("const QString &")
def slotStateDiff(self, text):
"""
Method to set diffusion choice for the coefficient
"""
choice = self.modelDiff.dicoV2M[str(text)]
log.debug("slotStateDiff -> %s" % (text))
if choice != 'variable':
self.pushButtonDiff.setEnabled(False)
setGreenColor(self.pushButtonDiff, False)
else:
self.pushButtonDiff.setEnabled(True)
setGreenColor(self.pushButtonDiff, True)
self.m_sca.setScalarDiffusivityChoice(self.scalar, choice)
@pyqtSignature("const QString &")
def slotNameDiff(self, text):
"""
Method to set the variance scalar choosed
"""
choice = self.modelNameDiff.dicoV2M[str(text)]
log.debug("slotStateDiff -> %s" % (text))
self.scalar = str(text)
self.lineEditDiff.setText(str(self.m_sca.getScalarDiffusivityInitialValue(self.scalar)))
mdl = self.m_sca.getScalarDiffusivityChoice(self.scalar)
self.modelDiff.setItem(str_model=mdl)
if mdl!= 'variable':
self.pushButtonDiff.setEnabled(False)
setGreenColor(self.pushButtonDiff, False)
else:
self.pushButtonDiff.setEnabled(True)
setGreenColor(self.pushButtonDiff, True)
def __changeChoice(self, text, sym, tag):
"""
Input variable state
"""
__model = getattr(self, "model" + sym)
__line = getattr(self, "lineEdit" + sym)
__combo = getattr(self, "comboBox" + sym)
__label = getattr(self, "label" + sym)
__button = getattr(self, "pushButton" + sym)
__labelu = getattr(self, "labelUnit" + sym)
__labelv = getattr(self, "labelVar" + sym)
choice = __model.dicoV2M[text]
log.debug("__changeChoice -> %s, %s" % (text, choice))
if choice != 'variable':
__button.setEnabled(False)
setGreenColor(__button, False)
else:
__button.setEnabled(True)
setGreenColor(__button, True)
if choice == 'thermal_law':
__line.hide()
__label.hide()
__labelu.hide()
__labelv.hide()
else:
__line.show()
__label.show()
__labelu.show()
__labelv.show()
self.mdl.setPropertyMode(tag, choice)
@pyqtSignature("const QString &")
def slotRho(self, text):
"""
Update the density
"""
if self.sender().validator().state == QValidator.Acceptable:
rho = from_qvariant(text, float)
self.mdl.setInitialValueDensity(rho)
@pyqtSignature("const QString &")
def slotMu(self, text):
"""
Update the molecular viscosity
"""
if self.sender().validator().state == QValidator.Acceptable:
mu = from_qvariant(text, float)
self.mdl.setInitialValueViscosity(mu)
@pyqtSignature("const QString &")
def slotCp(self, text):
"""
Update the specific heat
"""
if self.sender().validator().state == QValidator.Acceptable:
cp = from_qvariant(text, float)
self.mdl.setInitialValueHeat(cp)
@pyqtSignature("const QString &")
def slotViscv0(self, text):
"""
Update the volumic viscosity
"""
if self.sender().validator().state == QValidator.Acceptable:
viscv0 = from_qvariant(text, float)
self.mdl.setInitialValueVolumicViscosity(viscv0)
@pyqtSignature("const QString &")
def slotAl(self, text):
"""
Update the thermal conductivity
"""
if self.sender().validator().state == QValidator.Acceptable:
al = from_qvariant(text, float)
self.mdl.setInitialValueCond(al)
@pyqtSignature("const QString &")
def slotDiftl0(self, text):
"""
Update the thermal conductivity
"""
if self.sender().validator().state == QValidator.Acceptable:
diftl0 = from_qvariant(text, float)
self.mdl.setInitialValueDyn(diftl0)
@pyqtSignature("const QString &")
def slotDiff(self, text):
"""
Update the thermal conductivity
"""
if self.sender().validator().state == QValidator.Acceptable:
diff = from_qvariant(text, float)
self.m_sca.setScalarDiffusivityInitialValue(self.scalar, diff)
@pyqtSignature("")
def slotFormulaRho(self):
"""
User formula for density
"""
exp = self.mdl.getFormula('density')
req = [('density', 'Density')]
self.m_th = ThermalScalarModel(self.case)
s = self.m_th.getThermalScalarName()
mdl = self.m_th.getThermalScalarModel()
if mdl == "temperature_celsius":
TempInContext = "("+s+" + 273.15)"
exa = FluidCharacteristicsView.density.replace("temperature", TempInContext)
elif mdl == "enthalpy":
exa = FluidCharacteristicsView.density_h
else:
exa = FluidCharacteristicsView.density
symbols_rho = []
for s in self.list_scalars:
symbols_rho.append(s)
rho0_value = self.mdl.getInitialValueDensity()
ref_pressure = ReferenceValuesModel(self.case).getPressure()
symbols_rho.append(('rho0', 'Density (reference value) = ' + str(rho0_value)))
symbols_rho.append(('p0', 'Reference pressure = ' + str(ref_pressure)))
dialog = QMeiEditorView(self,
check_syntax = self.case['package'].get_check_syntax(),
expression = exp,
required = req,
symbols = symbols_rho,
examples = exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaRho -> %s" % str(result))
self.mdl.setFormula('density', result)
setGreenColor(self.sender(), False)
@pyqtSignature("")
def slotFormulaMu(self):
"""
User formula for molecular viscosity
"""
exp = self.mdl.getFormula('molecular_viscosity')
req = [('molecular_viscosity', 'Molecular Viscosity')]
self.m_th = ThermalScalarModel(self.case)
s = self.m_th.getThermalScalarName()
mdl = self.m_th.getThermalScalarModel()
if mdl == "temperature_celsius":
TempInContext = "("+s+" + 273.15)"
exa = FluidCharacteristicsView.molecular_viscosity.replace("temperature", TempInContext)
elif mdl == "enthalpy":
exa = FluidCharacteristicsView.molecular_viscosity_h
else:
exa = FluidCharacteristicsView.molecular_viscosity
symbols_mu = []
for s in self.list_scalars:
symbols_mu.append(s)
mu0_value = self.mdl.getInitialValueViscosity()
rho0_value = self.mdl.getInitialValueDensity()
ref_pressure = ReferenceValuesModel(self.case).getPressure()
symbols_mu.append(('mu0', 'Viscosity (reference value) = ' + str(mu0_value)))
symbols_mu.append(('rho0', 'Density (reference value) = ' + str(rho0_value)))
symbols_mu.append(('p0', 'Reference pressure = ' + str(ref_pressure)))
symbols_mu.append(('rho', 'Density'))
if CompressibleModel(self.case).getCompressibleModel() == 'on':
symbols_mu.append(('T', 'Temperature'))
ref_temperature = ReferenceValuesModel(self.case).getTemperature()
symbols_mu.append(('t0', 'Reference temperature = '+str(ref_temperature)+' K'))
dialog = QMeiEditorView(self,
check_syntax = self.case['package'].get_check_syntax(),
expression = exp,
required = req,
symbols = symbols_mu,
examples = exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaMu -> %s" % str(result))
self.mdl.setFormula('molecular_viscosity', result)
setGreenColor(self.sender(), False)
@pyqtSignature("")
def slotFormulaCp(self):
"""
User formula for specific heat
"""
exp = self.mdl.getFormula('specific_heat')
req = [('specific_heat', 'Specific heat')]
exa = FluidCharacteristicsView.specific_heat
symbols_cp = []
for s in self.list_scalars:
symbols_cp.append(s)
cp0_value = self.mdl.getInitialValueHeat()
ref_pressure = ReferenceValuesModel(self.case).getPressure()
symbols_cp.append(('cp0', 'Specific heat (reference value) = ' + str(cp0_value)))
symbols_cp.append(('p0', 'Reference pressure = ' + str(ref_pressure)))
dialog = QMeiEditorView(self,
check_syntax = self.case['package'].get_check_syntax(),
expression = exp,
required = req,
symbols = symbols_cp,
examples = exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaRho -> %s" % str(result))
self.mdl.setFormula('specific_heat', result)
setGreenColor(self.sender(), False)
@pyqtSignature("")
def slotFormulaViscv0(self):
"""
User formula for volumic viscosity
"""
exp = self.mdl.getFormula('volume_viscosity')
req = [('volume_viscosity', 'Volumic viscosity')]
exa = FluidCharacteristicsView.volume_viscosity
symbols_viscv0 = []
for s in self.list_scalars:
symbols_viscv0.append(s)
viscv0_value = self.mdl.getInitialValueVolumicViscosity()
ref_pressure = ReferenceValuesModel(self.case).getPressure()
ref_temperature = ReferenceValuesModel(self.case).getTemperature()
symbols_viscv0.append(('viscv0', 'Volumic viscosity (reference value) = '+str(viscv0_value)+' J/kg/K'))
symbols_viscv0.append(('p0', 'Reference pressure = '+str(ref_pressure)+' Pa'))
symbols_viscv0.append(('t0', 'Reference temperature = '+str(ref_temperature)+' K'))
symbols_viscv0.append(('T', 'Temperature'))
dialog = QMeiEditorView(self,
check_syntax = self.case['package'].get_check_syntax(),
expression = exp,
required = req,
symbols = symbols_viscv0,
examples = exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaViscv0 -> %s" % str(result))
self.mdl.setFormula('volume_viscosity', result)
setGreenColor(self.sender(), False)
@pyqtSignature("")
def slotFormulaAl(self):
"""
User formula for thermal conductivity
"""
exp = self.mdl.getFormula('thermal_conductivity')
req = [('thermal_conductivity', 'Thermal conductivity')]
self.m_th = ThermalScalarModel(self.case)
s = self.m_th.getThermalScalarName()
mdl = self.m_th.getThermalScalarModel()
if mdl == "temperature_celsius":
TempInContext = "("+s+" + 273.15)"
exa = FluidCharacteristicsView.thermal_conductivity.replace("temperature", TempInContext)
elif mdl == "enthalpy":
exa = FluidCharacteristicsView.thermal_conductivity_h
else:
exa = FluidCharacteristicsView.thermal_conductivity
symbols_al = []
for s in self.list_scalars:
symbols_al.append(s)
lambda0_value = self.mdl.getInitialValueCond()
ref_pressure = ReferenceValuesModel(self.case).getPressure()
symbols_al.append(('lambda0', 'Thermal conductivity (reference value) = ' + str(lambda0_value)))
symbols_al.append(('p0', 'Reference pressure = ' + str(ref_pressure)))
dialog = QMeiEditorView(self,
check_syntax = self.case['package'].get_check_syntax(),
expression = exp,
required = req,
symbols = symbols_al,
examples = exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaAl -> %s" % str(result))
self.mdl.setFormula('thermal_conductivity', result)
setGreenColor(self.sender(), False)
@pyqtSignature("")
def slotFormulaDiff(self):
"""
User formula for the diffusion coefficient
"""
name = self.m_sca.getScalarDiffusivityName(self.scalar)
exp = self.m_sca.getDiffFormula(self.scalar)
req = [(str(name), str(self.scalar)+'diffusion coefficient')]
exa = ''
sym = [('x','cell center coordinate'),
('y','cell center coordinate'),
('z','cell center coordinate'),]
sym.append((str(self.scalar),str(self.scalar)))
diff0_value = self.m_sca.getScalarDiffusivityInitialValue(self.scalar)
sym.append((str(name)+'_ref', str(self.scalar)+' diffusion coefficient (reference value) = '+str(diff0_value)+' m^2/s'))
dialog = QMeiEditorView(self,
check_syntax = self.case['package'].get_check_syntax(),
expression = exp,
required = req,
symbols = sym,
examples = exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaDiff -> %s" % str(result))
self.m_sca.setDiffFormula(self.scalar, result)
setGreenColor(self.pushButtonDiff, False)
def tr(self, text):
"""
Translation
"""
return text
class FluidCharacteristicsView(QWidget, Ui_FluidCharacteristicsForm):
"""
Class to open Molecular Properties Page.
"""
density = """# Density of air
density = 1.293*(273.15 / temperature);
# density for mixtures of gases
#
# Y1 -> mass fraction of component 1
# Y2 -> mass fraction of component 2
rho1 = 1.25051;
rho2 = 1.7832;
A = (Y1 / rho1) + (Y2 /rho2);
density = 1.0 / A;
"""
density_h = """# Density
density = enthalpy / 1040. * 1.29;
# density for mixtures of gases
#
# Y1 -> mass fraction of component 1
# Y2 -> mass fraction of component 2
rho1 = 1.25051;
rho2 = 1.7832;
A = (Y1 / rho1) + (Y2 /rho2);
density = 1.0 / A;
"""
density_wo = """density = 1.25051;
"""
molecular_viscosity = """# Sutherland's Formula
# Gas Cst T0 mu0
# air 120 291.15 18.27e-6
# nitrogen 111 300.55 17.81e-6
# oxygen 127 292.25 20.18e-6
# carbon dioxide 240 293.15 14.8e-6
# carbon monoxide 118 288.15 17.2e-6
# hydrogen 72 293.85 8.76e-6
# ammonia 370 293.15 9.82e-6
# sulfur dioxide 416 293.65 12.54e-6
# helium 79.4 273 19e-6
CST = 120;
T0 = 291.15;
mu_ref = 18.27e-6;
if ( temperature > 0 && temperature < 555) {
molecular_viscosity = mu_ref * ((T0+CST) / (temperature+CST)) * (temperature/T0)^(3./2.);
} else {
molecular_viscosity = -999.0;
}
"""
molecular_viscosity_h = """CST = 120;
T0 = 291.15;
mu_ref = 18.27e-6;
temperature = enthalpy / 1040.;
if ( enthalpy > 0) {
molecular_viscosity = mu_ref * (T0+CST / temperature+CST) * (temperature/T0)^(3./2.);
} else {
molecular_viscosity = -999.0;
}
"""
molecular_viscosity_wo = """CST = 120;
T0 = 291.15;
mu_ref = 18.27e-6;
molecular_viscosity = mu_ref * (T0+CST);
"""
specific_heat = """# specific heat for mixtures of gases
#
# Y1 -> mass fraction of component 1
# Y2 -> mass fraction of component 2
Cp1 = 520.3;
Cp2 = 1040.0;
specific_heat = Y1 * Cp1 + Y2 *Cp2;
"""
volume_viscosity = """# volume_viscosity
"""
thermal_conductivity = """# oxygen
thermal_conductivity = 6.2e-5 * temperature + 8.1e-3;
# nitrogen
thermal_conductivity = 6.784141e-5 * temperature + 5.564317e-3;
# hydrogen
thermal_conductivity = 4.431e-4 * temperature + 5.334e-2;
"""
thermal_conductivity_h = """temperature = enthalpy / 1040.;
thermal_conductivity = 6.2e-5 * temperature + 8.1e-3;
"""
def __init__(self, parent, case):
"""
Constructor
"""
QWidget.__init__(self, parent)
Ui_FluidCharacteristicsForm.__init__(self)
self.setupUi(self)
self.case = case
self.case.undoStopGlobal()
self.mdl = FluidCharacteristicsModel(self.case)
self.notebook = NotebookModel(self.case)
if EOS == 1:
self.ava = eosAva.EosAvailable()
import cs_config
cfg = cs_config.config()
self.freesteam = 0
if cfg.libs['freesteam'].have != "no":
self.freesteam = 1
self.m_th = ThermalScalarModel(self.case)
s = self.m_th.getThermalScalarName()
tsm = self.mdl.tsm
# Particular Widget initialization taking into account of "Calculation Features"
mdl_atmo, mdl_joule, mdl_thermal, mdl_gas, mdl_coal, mdl_comp = self.mdl.getThermoPhysicalModel(
)
# Combo models
self.modelRho = ComboModel(self.comboBoxRho, 3, 1)
self.modelMu = ComboModel(self.comboBoxMu, 3, 1)
self.modelCp = ComboModel(self.comboBoxCp, 3, 1)
self.modelAl = ComboModel(self.comboBoxAl, 3, 1)
self.modelDiff = ComboModel(self.comboBoxDiff, 2, 1)
self.modelNameDiff = ComboModel(self.comboBoxNameDiff, 1, 1)
self.modelViscv0 = ComboModel(self.comboBoxViscv0, 3, 1)
self.modelDiftl0 = ComboModel(self.comboBoxDiftl0, 3, 1)
self.modelMaterial = ComboModel(self.comboBoxMaterial, 1, 1)
self.modelMethod = ComboModel(self.comboBoxMethod, 1, 1)
self.modelPhas = ComboModel(self.comboBoxPhas, 2, 1)
self.modelRho.addItem(self.tr('constant'), 'constant')
self.modelRho.addItem(self.tr('user law'), 'user_law')
self.modelRho.addItem(self.tr('material law'), 'thermal_law')
if mdl_atmo != 'off':
self.modelRho.addItem(self.tr('defined in atphyv'),
'predefined_law')
elif mdl_joule == 'arc':
self.modelRho.addItem(self.tr('defined in elphyv'),
'predefined_law')
elif mdl_comp != 'off':
self.modelRho.addItem(self.tr('predefined law'), 'predefined_law')
elif mdl_gas != 'off' or mdl_coal != 'off':
self.modelRho.addItem(self.tr('predefined law'), 'predefined_law')
self.modelMu.addItem(self.tr('constant'), 'constant')
self.modelMu.addItem(self.tr('user law'), 'user_law')
self.modelMu.addItem(self.tr('material law'), 'thermal_law')
if mdl_joule == 'arc':
self.modelMu.addItem(self.tr('defined in elphyv'),
'predefined_law')
self.modelCp.addItem(self.tr('constant'), 'constant')
self.modelCp.addItem(self.tr('user law'), 'user_law')
self.modelCp.addItem(self.tr('material law'), 'thermal_law')
if mdl_joule == 'arc':
self.modelCp.addItem(self.tr('defined in elphyv'),
'predefined_law')
self.modelAl.addItem(self.tr('constant'), 'constant')
self.modelAl.addItem(self.tr('user law'), 'user_law')
self.modelAl.addItem(self.tr('material law'), 'thermal_law')
if mdl_joule == 'arc':
self.modelAl.addItem(self.tr('defined in elphyv'),
'predefined_law')
self.modelDiff.addItem(self.tr('constant'), 'constant')
self.modelDiff.addItem(self.tr('user law'), 'user_law')
self.modelViscv0.addItem(self.tr('constant'), 'constant')
self.modelViscv0.addItem(self.tr('user law'), 'user_law')
self.modelViscv0.addItem(self.tr('material law'), 'thermal_law')
self.modelDiftl0.addItem(self.tr('constant'), 'constant')
self.modelDiftl0.addItem(self.tr('user law'), 'user_law')
self.modelDiftl0.addItem(self.tr('material law'), 'thermal_law')
if mdl_gas != 'off' or mdl_coal != 'off':
self.modelDiftl0.addItem(self.tr('predefined law'),
'predefined_law')
self.modelPhas.addItem(self.tr('liquid'), 'liquid')
self.modelPhas.addItem(self.tr('gas'), 'gas')
self.scalar = ""
scalar_list = self.mdl.m_sca.getUserScalarNameList()
for s in self.mdl.m_sca.getScalarsVarianceList():
if s in scalar_list: scalar_list.remove(s)
if scalar_list != []:
self.scalar = scalar_list[0]
for scalar in scalar_list:
self.modelNameDiff.addItem(scalar)
# Validators
validatorP0 = DoubleValidator(self.lineEditP0, min=0.0)
self.lineEditP0.setValidator(validatorP0)
validatorT0 = DoubleValidator(self.lineEditT0, min=0.0)
validatorT0.setExclusiveMin(True)
self.lineEditT0.setValidator(validatorT0)
validatorOxydant = DoubleValidator(self.lineEditOxydant, min=0.0)
validatorOxydant.setExclusiveMin(True)
self.lineEditOxydant.setValidator(validatorOxydant)
validatorFuel = DoubleValidator(self.lineEditFuel, min=0.0)
validatorFuel.setExclusiveMin(True)
self.lineEditFuel.setValidator(validatorFuel)
validatorMM = DoubleValidator(self.lineEditMassMolar, min=0.0)
validatorMM.setExclusiveMin(True)
self.lineEditMassMolar.setValidator(validatorMM)
validatorRho = DoubleValidator(self.lineEditRho, min=0.0)
validatorMu = DoubleValidator(self.lineEditMu, min=0.0)
validatorCp = DoubleValidator(self.lineEditCp, min=0.0)
validatorAl = DoubleValidator(self.lineEditAl, min=0.0)
validatorDiff = DoubleValidator(self.lineEditDiff, min=0.0)
validatorViscv0 = DoubleValidator(self.lineEditViscv0, min=0.0)
validatorDiftl0 = DoubleValidator(self.lineEditDiftl0, min=0.0)
validatorRho.setExclusiveMin(True)
validatorMu.setExclusiveMin(True)
validatorCp.setExclusiveMin(True)
validatorAl.setExclusiveMin(True)
validatorDiff.setExclusiveMin(True)
validatorDiftl0.setExclusiveMin(True)
self.lineEditRho.setValidator(validatorRho)
self.lineEditMu.setValidator(validatorMu)
self.lineEditCp.setValidator(validatorCp)
self.lineEditAl.setValidator(validatorAl)
self.lineEditDiff.setValidator(validatorDiff)
self.lineEditViscv0.setValidator(validatorViscv0)
self.lineEditDiftl0.setValidator(validatorDiftl0)
# Connections
self.lineEditP0.textChanged[str].connect(self.slotPressure)
self.lineEditT0.textChanged[str].connect(self.slotTemperature)
self.lineEditOxydant.textChanged[str].connect(self.slotTempOxydant)
self.lineEditFuel.textChanged[str].connect(self.slotTempFuel)
self.lineEditMassMolar.textChanged[str].connect(self.slotMassemol)
self.comboBoxRho.currentIndexChanged[str].connect(self.slotStateRho)
self.comboBoxMu.currentIndexChanged[str].connect(self.slotStateMu)
self.comboBoxCp.currentIndexChanged[str].connect(self.slotStateCp)
self.comboBoxAl.currentIndexChanged[str].connect(self.slotStateAl)
self.comboBoxDiff.currentIndexChanged[str].connect(self.slotStateDiff)
self.comboBoxNameDiff.activated[str].connect(self.slotNameDiff)
self.comboBoxViscv0.currentIndexChanged[str].connect(
self.slotStateViscv0)
self.comboBoxMaterial.activated[str].connect(self.slotMaterial)
self.comboBoxMethod.activated[str].connect(self.slotMethod)
self.comboBoxPhas.activated[str].connect(self.slotPhas)
self.lineEditRho.textChanged[str].connect(self.slotRho)
self.lineEditMu.textChanged[str].connect(self.slotMu)
self.lineEditCp.textChanged[str].connect(self.slotCp)
self.lineEditAl.textChanged[str].connect(self.slotAl)
self.lineEditDiff.textChanged[str].connect(self.slotDiff)
self.lineEditDiftl0.textChanged[str].connect(self.slotDiftl0)
self.lineEditViscv0.textChanged[str].connect(self.slotViscv0)
self.pushButtonRho.clicked.connect(self.slotFormulaRho)
self.pushButtonMu.clicked.connect(self.slotFormulaMu)
self.pushButtonCp.clicked.connect(self.slotFormulaCp)
self.pushButtonAl.clicked.connect(self.slotFormulaAl)
self.pushButtonDiff.clicked.connect(self.slotFormulaDiff)
self.pushButtonViscv0.clicked.connect(self.slotFormulaViscv0)
self.initializeWidget()
self.case.undoStartGlobal()
def initializeWidget(self):
"""
"""
mdls = self.mdl.getThermoPhysicalModel()
mdl_atmo, mdl_joule, mdl_thermal, mdl_gas, mdl_coal, mdl_comp = mdls
self.groupBoxMassMolar.hide()
if mdl_atmo != "off":
self.labelInfoT0.hide()
elif mdl_comp != "off" or mdl_coal != "off":
m = self.mdl.getMassemol()
self.lineEditMassMolar.setText(str(m))
self.groupBoxMassMolar.show()
if mdl_thermal != "off":
t = self.mdl.getTemperature()
self.lineEditT0.setText(str(t))
if mdl_thermal == "temperature_celsius":
self.labelUnitT0.setText("\xB0 C")
if self.mdl.getMaterials() != "user_material":
self.labelInfoT0.hide()
else:
self.groupBoxTemperature.hide()
if mdl_gas == 'd3p':
self.groupBoxTempd3p.show()
t_oxy = self.mdl.getTempOxydant()
t_fuel = self.mdl.getTempFuel()
self.lineEditOxydant.setText(str(t_oxy))
self.lineEditFuel.setText(str(t_fuel))
else:
self.groupBoxTempd3p.hide()
darc = GroundwaterModel(self.case).getGroundwaterModel()
if darc != 'off':
self.groupBoxPressure.hide()
else:
p = self.mdl.getPressure()
self.lineEditP0.setText(str(p))
if (self.freesteam == 1 or EOS == 1 or coolprop_fluids):
self.tables = True
else:
self.tables = False
if self.tables == False or mdl_joule != 'off' or mdl_comp != 'off':
self.groupBoxTableChoice.hide()
else:
self.groupBoxTableChoice.show()
self.lineEditReference.setEnabled(False)
# suppress perfect gas
self.modelMaterial.addItem(self.tr('user material'),
'user_material')
tmp = ["Argon", "Nitrogen", "Hydrogen", "Oxygen", "Helium", "Air"]
if EOS == 1:
fls = self.ava.whichFluids()
for fli in fls:
if fli not in tmp:
tmp.append(fli)
self.modelMaterial.addItem(self.tr(fli), fli)
if self.freesteam == 1 and EOS == 0:
self.modelMaterial.addItem(self.tr('Water'), 'Water')
if coolprop_fluids and EOS == 0:
have_coolprop = False
if cs_config.config().libs['coolprop'].have != "no":
have_coolprop = True
for fli in coolprop_fluids:
if self.freesteam == 1 and fli == 'Water':
continue
self.modelMaterial.addItem(self.tr(fli), fli,
coolprop_warn)
material = self.mdl.getMaterials()
self.modelMaterial.setItem(str_model=material)
self.updateMethod()
#compressible
self.groupBoxViscv0.hide()
# combustion
self.groupBoxDiftl0.hide()
if self.scalar == "":
self.groupBoxDiff.hide()
else:
self.groupBoxDiff.show()
self.lineEditDiff.setText(
str(
self.mdl.m_sca.getScalarDiffusivityInitialValue(
self.scalar)))
diff_choice = self.mdl.m_sca.getScalarDiffusivityChoice(
self.scalar)
self.modelDiff.setItem(str_model=diff_choice)
self.modelNameDiff.setItem(str_model=str(self.scalar))
if diff_choice != 'user_law':
self.pushButtonDiff.hide()
self.pushButtonDiff.setEnabled(False)
self.pushButtonDiff.setStyleSheet("background-color: None")
else:
self.pushButtonDiff.show()
self.pushButtonDiff.setEnabled(True)
name = self.mdl.m_sca.getScalarDiffusivityName(self.scalar)
exp = self.mdl.m_sca.getDiffFormula(self.scalar)
if exp:
self.pushButtonDiff.setStyleSheet(
"background-color: green")
self.pushButtonDiff.setToolTip(exp)
else:
self.pushButtonDiff.setStyleSheet("background-color: red")
# Standard Widget initialization
for tag, symbol in self.mdl.lst:
__model = getattr(self, "model" + symbol)
__line = getattr(self, "lineEdit" + symbol)
__button = getattr(self, "pushButton" + symbol)
__label = getattr(self, "label" + symbol)
__labelu = getattr(self, "labelUnit" + symbol)
if tag != 'dynamic_diffusion':
__labelv = getattr(self, "labelVar" + symbol)
c = self.mdl.getPropertyMode(tag)
if c not in __model.getItems():
c = 'constant'
self.mdl.setPropertyMode(tag, c)
__model.setItem(str_model=c)
if c == 'user_law':
__button.show()
__button.setEnabled(True)
__label.setText(self.tr("Reference value"))
else:
__button.hide()
__button.setEnabled(False)
__label.setText(self.tr("Reference value"))
if c == 'thermal_law':
__line.hide()
__label.hide()
__labelu.hide()
__labelv.hide()
else:
__line.show()
__label.show()
__labelu.show()
__labelv.show()
if self.mdl.getMaterials() == "user_material":
__model.disableItem(str_model='thermal_law')
else:
__model.enableItem(str_model='thermal_law')
else:
__label.setText(self.tr("Reference value"))
self.mdl.getInitialValue(tag)
__line.setText(str(self.mdl.getInitialValue(tag)))
# If no thermal scalar, hide specific heat and thermal conductivity.
if mdl_thermal == 'off':
self.groupBoxCp.hide()
self.groupBoxAl.hide()
if mdl_gas != 'off' or mdl_coal != 'off':
self.groupBoxDiftl0.show()
for tag, symbol in self.mdl.lst:
__model = getattr(self, "model" + symbol)
__line = getattr(self, "lineEdit" + symbol)
__button = getattr(self, "pushButton" + symbol)
__label = getattr(self, "label" + symbol)
__combo = getattr(self, "comboBox" + symbol)
# Gas or coal combustion
if mdl_gas != 'off' or mdl_coal != 'off':
if tag == 'density':
__model.setItem(str_model='predefined_law')
__combo.setEnabled(False)
__button.setEnabled(False)
__button.hide()
self.mdl.setPropertyMode(tag, 'predefined_law')
__label.setText(
self.tr("Calculation by\n perfect gas law"))
__line.setText(str(""))
__line.setEnabled(False)
elif tag == 'dynamic_diffusion':
__model.setItem(str_model='predefined_law')
__combo.setEnabled(False)
__button.setEnabled(False)
__button.hide()
else:
__model.setItem(str_model='constant')
self.mdl.setPropertyMode(tag, 'constant')
# Joule
if mdl_joule == 'arc':
__model.disableItem(str_model='constant')
__model.disableItem(str_model='predefined_law')
__model.setItem(str_model='predefined_law')
__combo.setEnabled(False)
__button.setEnabled(False)
__button.hide()
self.mdl.setPropertyMode(tag, 'predefined_law')
if mdl_joule == 'joule':
__model.setItem(str_model='user_law')
__model.disableItem(str_model='constant')
self.mdl.setPropertyMode(tag, 'user_law')
# Atmospheric Flows
if mdl_atmo != 'off':
if tag == 'density':
__model.disableItem(str_model='constant')
__model.disableItem(str_model='predefined_law')
__model.setItem(str_model='predefined_law')
__combo.setEnabled(False)
__button.setEnabled(False)
__button.hide()
# Compressible Flows
if mdl_comp != 'off':
self.groupBoxViscv0.show()
if tag == 'density':
__model.setItem(str_model='predefined_law')
__combo.setEnabled(False)
__button.setEnabled(False)
__button.hide()
__combo.hide()
__button.hide()
self.mdl.setPropertyMode(tag, 'predefined_law')
__line.setEnabled(True)
if tag == 'specific_heat':
__model.setItem(str_model='constant')
__combo.setEnabled(False)
__button.setEnabled(False)
__button.hide()
self.mdl.setPropertyMode(tag, 'constant')
self.groupBoxCp.setTitle('Isobaric specific heat')
if tag == 'volume_viscosity':
__combo.setEnabled(True)
c = self.mdl.getPropertyMode(tag)
if c == 'predefined_law':
__button.setEnabled(True)
else:
__button.setEnabled(False)
else:
if tag == 'specific_heat':
self.groupBoxCp.setTitle('Specific heat')
def updateTypeChoice(self, old_choice):
"""
add/suppress thermo tables for each properties
"""
for tag, symbol in self.mdl.lst:
__model = getattr(self, "model" + symbol)
if self.mdl.getMaterials() == "user_material":
__model.disableItem(str_model='thermal_law')
else:
__model.enableItem(str_model='thermal_law')
if old_choice == "user_material":
self.mdl.setPropertyMode(tag, 'thermal_law')
self.__changeChoice(str("material law"), symbol, tag)
c = self.mdl.getPropertyMode(tag)
__model.setItem(str_model=c)
def updateMethod(self):
"""
update method list with material choice
"""
for nb in range(len(self.modelMethod.getItems())):
self.modelMethod.delItem(0)
self.comboBoxPhas.hide()
self.labelPhas.hide()
if self.mdl.getMaterials() == "user_material":
self.modelMethod.addItem(self.tr('user properties'),
'user_properties')
else:
if EOS == 1:
material = self.mdl.getMaterials()
self.ava.setMethods(material)
fls = self.ava.whichMethods()
for fli in fls:
self.modelMethod.addItem(self.tr(fli), fli)
if self.mdl.getMethod() != "freesteam" and self.mdl.getMethod(
) != "CoolProp":
self.comboBoxPhas.show()
self.labelPhas.show()
if self.freesteam == 1 and self.mdl.getMaterials() == "Water":
self.modelMethod.addItem(self.tr("freesteam"), "freesteam")
if self.mdl.getMaterials() in coolprop_fluids:
self.modelMethod.addItem(self.tr("CoolProp"), "CoolProp")
# update comboBoxMethod
method = self.mdl.getMethod()
self.modelMethod.setItem(str_model=method)
self.updateReference()
def updateReference(self):
"""
update Reference with material, method and field nature choice
"""
# update lineEditReference
self.lineEditReference.setText(self.mdl.getReference())
@pyqtSlot(str)
def slotPressure(self, text):
"""
Input PRESS.
"""
if self.lineEditP0.validator().state == QValidator.Acceptable:
p = from_qvariant(text, float)
self.mdl.setPressure(p)
@pyqtSlot(str)
def slotTemperature(self, text):
"""
Input TEMPERATURE.
"""
if self.lineEditT0.validator().state == QValidator.Acceptable:
t = from_qvariant(text, float)
self.mdl.setTemperature(t)
@pyqtSlot(str)
def slotTempOxydant(self, text):
"""
Input oxydant TEMPERATURE.
"""
if self.lineEditOxydant.validator().state == QValidator.Acceptable:
t = from_qvariant(text, float)
self.mdl.setTempOxydant(t)
@pyqtSlot(str)
def slotTempFuel(self, text):
"""
Input fuel TEMPERATURE.
"""
if self.lineEditFuel.validator().state == QValidator.Acceptable:
t = from_qvariant(text, float)
self.mdl.setTempFuel(t)
@pyqtSlot(str)
def slotMassemol(self, text):
"""
Input Mass molar.
"""
if self.lineEditMassMolar.validator().state == QValidator.Acceptable:
m = from_qvariant(text, float)
self.mdl.setMassemol(m)
@pyqtSlot(str)
def slotMaterial(self, text):
"""
Method to call 'setMaterial'
"""
choice = self.modelMaterial.dicoV2M[str(text)]
old_choice = self.mdl.getMaterials()
self.mdl.setMaterials(choice)
self.updateMethod()
self.updateTypeChoice(old_choice)
@pyqtSlot(str)
def slotPhas(self, text):
"""
Method to call 'setFieldNature'
"""
choice = self.modelPhas.dicoV2M[str(text)]
self.mdl.setFieldNature(choice)
self.updateReference()
@pyqtSlot(str)
def slotMethod(self, text):
"""
Method to call 'setMethod'
"""
choice = self.modelMethod.dicoV2M[str(text)]
self.mdl.setMethod(choice)
self.comboBoxPhas.hide()
self.labelPhas.hide()
if self.mdl.getMaterials() != "user_material" and \
self.mdl.getMethod() != "freesteam" and \
self.mdl.getMethod() != "CoolProp":
self.comboBoxPhas.show()
self.labelPhas.show()
self.updateReference()
@pyqtSlot(str)
def slotStateRho(self, text):
"""
Method to call 'getState' with correct arguements for 'rho'
"""
self.__changeChoice(str(text), 'Rho', 'density')
@pyqtSlot(str)
def slotStateMu(self, text):
"""
Method to call 'getState' with correct arguements for 'Mu'
"""
self.__changeChoice(str(text), 'Mu', 'molecular_viscosity')
@pyqtSlot(str)
def slotStateCp(self, text):
"""
Method to call 'getState' with correct arguements for 'Cp'
"""
self.__changeChoice(str(text), 'Cp', 'specific_heat')
@pyqtSlot(str)
def slotStateViscv0(self, text):
"""
Method to call 'getState' with correct arguements for 'Viscv0'
"""
self.__changeChoice(str(text), 'Viscv0', 'volume_viscosity')
@pyqtSlot(str)
def slotStateAl(self, text):
"""
Method to call 'getState' with correct arguements for 'Al'
"""
self.__changeChoice(str(text), 'Al', 'thermal_conductivity')
@pyqtSlot(str)
def slotStateDiff(self, text):
"""
Method to set diffusion choice for the coefficient
"""
choice = self.modelDiff.dicoV2M[str(text)]
log.debug("slotStateDiff -> %s" % (text))
if choice != 'user_law':
self.pushButtonDiff.setEnabled(False)
self.pushButtonDiff.setStyleSheet("background-color: None")
self.pushButtonDiff.hide()
else:
self.pushButtonDiff.show()
self.pushButtonDiff.setEnabled(True)
name = self.mdl.m_sca.getScalarDiffusivityName(self.scalar)
exp = self.mdl.m_sca.getDiffFormula(self.scalar)
if exp:
self.pushButtonDiff.setStyleSheet("background-color: green")
self.pushButtonDiff.setToolTip(exp)
else:
self.pushButtonDiff.setStyleSheet("background-color: red")
self.mdl.m_sca.setScalarDiffusivityChoice(self.scalar, choice)
@pyqtSlot(str)
def slotNameDiff(self, text):
"""
Method to set the variance scalar choosed
"""
choice = self.modelNameDiff.dicoV2M[str(text)]
log.debug("slotStateDiff -> %s" % (text))
self.scalar = str(text)
self.lineEditDiff.setText(
str(self.mdl.m_sca.getScalarDiffusivityInitialValue(self.scalar)))
mdl = self.mdl.m_sca.getScalarDiffusivityChoice(self.scalar)
self.modelDiff.setItem(str_model=mdl)
if mdl != 'user_law':
self.pushButtonDiff.hide()
self.pushButtonDiff.setEnabled(False)
self.pushButtonDiff.setStyleSheet("background-color: None")
else:
self.pushButtonDiff.show()
self.pushButtonDiff.setEnabled(True)
name = self.mdl.m_sca.getScalarDiffusivityName(self.scalar)
exp = self.mdl.m_sca.getDiffFormula(self.scalar)
if exp:
self.pushButtonDiff.setStyleSheet("background-color: green")
self.pushButtonDiff.setToolTip(exp)
else:
self.pushButtonDiff.setStyleSheet("background-color: red")
def __changeChoice(self, text, sym, tag):
"""
Input variable state
"""
__model = getattr(self, "model" + sym)
__line = getattr(self, "lineEdit" + sym)
__combo = getattr(self, "comboBox" + sym)
__label = getattr(self, "label" + sym)
__button = getattr(self, "pushButton" + sym)
__labelu = getattr(self, "labelUnit" + sym)
__labelv = getattr(self, "labelVar" + sym)
choice = __model.dicoV2M[text]
log.debug("__changeChoice -> %s, %s" % (text, choice))
if choice != 'user_law':
__button.hide()
__button.setEnabled(False)
__button.setStyleSheet("background-color: None")
else:
__button.setEnabled(True)
__button.show()
exp = None
if sym == "Rho":
exp = self.mdl.getFormula('density')
elif sym == "Mu":
exp = self.mdl.getFormula('molecular_viscosity')
elif sym == "Cp":
exp = self.mdl.getFormula('specific_heat')
elif sym == "Viscv0":
exp = self.mdl.getFormula('volume_viscosity')
elif sym == "Al":
exp = self.mdl.getFormula('thermal_conductivity')
elif sym == "Diff":
name = self.mdl.m_sca.getScalarDiffusivityName(self.scalar)
exp = self.mdl.m_sca.getDiffFormula(self.scalar)
if exp:
__button.setStyleSheet("background-color: green")
__button.setToolTip(exp)
else:
__button.setStyleSheet("background-color: red")
if choice == 'thermal_law':
__line.hide()
__label.hide()
__labelu.hide()
__labelv.hide()
else:
__line.show()
__label.show()
__labelu.show()
__labelv.show()
self.mdl.setPropertyMode(tag, choice)
@pyqtSlot(str)
def slotRho(self, text):
"""
Update the density
"""
if self.lineEditRho.validator().state == QValidator.Acceptable:
rho = from_qvariant(text, float)
self.mdl.setInitialValueDensity(rho)
@pyqtSlot(str)
def slotMu(self, text):
"""
Update the molecular viscosity
"""
if self.lineEditMu.validator().state == QValidator.Acceptable:
mu = from_qvariant(text, float)
self.mdl.setInitialValueViscosity(mu)
@pyqtSlot(str)
def slotCp(self, text):
"""
Update the specific heat
"""
if self.lineEditCp.validator().state == QValidator.Acceptable:
cp = from_qvariant(text, float)
self.mdl.setInitialValueHeat(cp)
@pyqtSlot(str)
def slotViscv0(self, text):
"""
Update the volumic viscosity
"""
if self.lineEditViscv0.validator().state == QValidator.Acceptable:
viscv0 = from_qvariant(text, float)
self.mdl.setInitialValueVolumeViscosity(viscv0)
@pyqtSlot(str)
def slotAl(self, text):
"""
Update the thermal conductivity
"""
if self.lineEditAl.validator().state == QValidator.Acceptable:
al = from_qvariant(text, float)
self.mdl.setInitialValueCond(al)
@pyqtSlot(str)
def slotDiftl0(self, text):
"""
Update the thermal conductivity
"""
if self.lineEditDiftl0.validator().state == QValidator.Acceptable:
diftl0 = from_qvariant(text, float)
self.mdl.setInitialValueDyn(diftl0)
@pyqtSlot(str)
def slotDiff(self, text):
"""
Update the thermal conductivity
"""
if self.lineEditDiff.validator().state == QValidator.Acceptable:
diff = from_qvariant(text, float)
self.mdl.m_sca.setScalarDiffusivityInitialValue(self.scalar, diff)
@pyqtSlot()
def slotFormulaRho(self):
"""
User formula for density
"""
exp, req, sca, symbols_rho = self.mdl.getFormulaRhoComponents()
self.m_th = ThermalScalarModel(self.case)
s = self.m_th.getThermalScalarName()
mdl = self.m_th.getThermalScalarModel()
if mdl == "off":
exa = FluidCharacteristicsView.density_wo
elif mdl == "temperature_celsius":
TempInContext = "(" + s + " + 273.15)"
exa = FluidCharacteristicsView.density.replace(
"temperature", TempInContext)
elif mdl == "enthalpy":
exa = FluidCharacteristicsView.density_h
else:
exa = FluidCharacteristicsView.density
mci = mei_to_c_interpreter(self.case, False)
mci.init_cell_block(exp, req, symbols_rho, sca, 'density')
dialog = QMegEditorView(self,
mei_to_c=mci,
expression=exp,
required=req,
symbols=symbols_rho,
examples=exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaRho -> %s" % str(result))
self.mdl.setFormula('density', str(result))
self.pushButtonRho.setToolTip(result)
self.pushButtonRho.setStyleSheet("background-color: green")
@pyqtSlot()
def slotFormulaMu(self):
"""
User formula for molecular viscosity
"""
exp, req, sca, symbols_mu = self.mdl.getFormulaMuComponents()
self.m_th = ThermalScalarModel(self.case)
s = self.m_th.getThermalScalarName()
mdl = self.m_th.getThermalScalarModel()
if mdl == "off":
exa = FluidCharacteristicsView.molecular_viscosity_wo
elif mdl == "temperature_celsius":
TempInContext = "(" + s + " + 273.15)"
exa = FluidCharacteristicsView.molecular_viscosity.replace(
"temperature", TempInContext)
elif mdl == "enthalpy":
exa = FluidCharacteristicsView.molecular_viscosity_h
else:
exa = FluidCharacteristicsView.molecular_viscosity
mci = mei_to_c_interpreter(self.case, False)
mci.init_cell_block(exp, req, symbols_mu, sca, 'molecular_viscosity')
dialog = QMegEditorView(self,
mei_to_c=mci,
expression=exp,
required=req,
symbols=symbols_mu,
examples=exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaMu -> %s" % str(result))
self.mdl.setFormula('molecular_viscosity', str(result))
self.pushButtonMu.setToolTip(result)
self.pushButtonMu.setStyleSheet("background-color: green")
@pyqtSlot()
def slotFormulaCp(self):
"""
User formula for specific heat
"""
exp, req, sca, symbols_cp = self.mdl.getFormulaCpComponents()
exa = FluidCharacteristicsView.specific_heat
mci = mei_to_c_interpreter(self.case, False)
mci.init_cell_block(exp, req, symbols_cp, sca, 'specific_heat')
dialog = QMegEditorView(self,
mei_to_c=mci,
expression=exp,
required=req,
symbols=symbols_cp,
examples=exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaRho -> %s" % str(result))
self.mdl.setFormula('specific_heat', str(result))
self.pushButtonCp.setToolTip(result)
self.pushButtonCp.setStyleSheet("background-color: green")
@pyqtSlot()
def slotFormulaViscv0(self):
"""
User formula for volumic viscosity
"""
exp, req, sca, symbols_viscv0 = self.mdl.getFormulaViscv0Components()
exa = FluidCharacteristicsView.volume_viscosity
mci = mei_to_c_interpreter(self.case, False)
mci.init_cell_block(exp, req, symbols_viscv0, sca, 'volume_viscosity')
dialog = QMegEditorView(self,
mei_to_c=mci,
expression=exp,
required=req,
symbols=symbols_viscv0,
examples=exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaViscv0 -> %s" % str(result))
self.mdl.setFormula('volume_viscosity', str(result))
self.pushButtonViscv0.setToolTip(result)
self.pushButtonViscv0.setStyleSheet("background-color: green")
@pyqtSlot()
def slotFormulaAl(self):
"""
User formula for thermal conductivity
"""
exp, req, sca, symbols_al = self.mdl.getFormulaAlComponents()
self.m_th = ThermalScalarModel(self.case)
s = self.m_th.getThermalScalarName()
mdl = self.m_th.getThermalScalarModel()
if mdl == "temperature_celsius":
TempInContext = "(" + s + " + 273.15)"
exa = FluidCharacteristicsView.thermal_conductivity.replace(
"temperature", TempInContext)
elif mdl == "enthalpy":
exa = FluidCharacteristicsView.thermal_conductivity_h
else:
exa = FluidCharacteristicsView.thermal_conductivity
mci = mei_to_c_interpreter(self.case, False)
mci.init_cell_block(exp, req, symbols_al, sca, 'thermal_conductivity')
dialog = QMegEditorView(self,
mei_to_c=mci,
expression=exp,
required=req,
symbols=symbols_al,
examples=exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaAl -> %s" % str(result))
self.mdl.setFormula('thermal_conductivity', str(result))
self.pushButtonAl.setToolTip(result)
self.pushButtonAl.setStyleSheet("background-color: green")
@pyqtSlot()
def slotFormulaDiff(self):
"""
User formula for the diffusion coefficient
"""
exp, req, sca, sym = self.mdl.getFormulaDiffComponents(self.scalar)
exa = ''
dname = self.mdl.m_sca.getScalarDiffusivityName(self.scalar)
mci = mei_to_c_interpreter(self.case, False)
mci.init_cell_block(exp, req, sym, sca, dname)
dialog = QMegEditorView(self,
mei_to_c=mci,
expression=exp,
required=req,
symbols=sym,
examples=exa)
if dialog.exec_():
result = dialog.get_result()
log.debug("slotFormulaDiff -> %s" % str(result))
self.mdl.m_sca.setDiffFormula(self.scalar, str(result))
self.pushButtonDiff.setToolTip(result)
self.pushButtonDiff.setStyleSheet("background-color: green")
def tr(self, text):
"""
Translation
"""
return text
class ThermalScalarView(QWidget, Ui_ThermalScalarForm):
"""
Class to open Thermal Scalar Transport Page.
"""
def __init__(self, parent, case, tree):
"""
Constructor
"""
QWidget.__init__(self, parent)
Ui_ThermalScalarForm.__init__(self)
self.setupUi(self)
self.browser = tree
self.case = case
self.case.undoStopGlobal()
self.thermal = ThermalScalarModel(self.case)
# combo Model
self.modelThermal = ComboModel(self.comboBoxThermal, 4, 1)
self.modelThermal.addItem(self.tr("No thermal scalar"), 'off')
self.modelThermal.addItem(self.tr("Temperature (Celsius)"), 'temperature_celsius')
self.modelThermal.addItem(self.tr("Temperature (Kelvin)"), 'temperature_kelvin')
self.modelThermal.addItem(self.tr("Enthalpy (J/kg)"), 'enthalpy')
self.modelThermal.addItem(self.tr("Potential temperature"), 'potential_temperature')
self.modelThermal.addItem(self.tr("Liquid potential temperature"), 'liquid_potential_temperature')
self.modelThermal.addItem(self.tr("Total energy (J/kg)"), 'total_energy')
self.connect(self.comboBoxThermal, SIGNAL("activated(const QString&)"), self.slotThermalScalar)
# Update the thermal scalar list with the calculation features
for sca in self.thermal.thermalModel:
if sca not in self.thermal.thermalScalarModelsList():
self.modelThermal.disableItem(str_model=sca)
if ElectricalModel(self.case).getElectricalModel() != 'off':
self.comboBoxThermal.setEnabled(False)
if CoalCombustionModel(self.case).getCoalCombustionModel() != 'off':
self.comboBoxThermal.setEnabled(False)
if GasCombustionModel(self.case).getGasCombustionModel() != 'off':
self.comboBoxThermal.setEnabled(False)
if CompressibleModel(self.case).getCompressibleModel() != 'off':
self.comboBoxThermal.setEnabled(False)
else:
self.modelThermal.delItem(6)
if AtmosphericFlowsModel(self.case).getAtmosphericFlowsModel() != 'off':
self.comboBoxThermal.setEnabled(False)
else:
self.modelThermal.delItem(5)
self.modelThermal.delItem(4)
# Select the thermal scalar model
model = self.thermal.getThermalScalarModel()
self.modelThermal.setItem(str_model=model)
self.case.undoStartGlobal()
@pyqtSignature("const QString &")
def slotThermalScalar(self, text):
"""
Update the thermal scalar markup.
"""
th = self.modelThermal.dicoV2M[str(text)]
self.thermal.setThermalModel(th)
self.browser.configureTree(self.case)
def tr(self, text):
"""
Translation
"""
return text
