def input_properties_case_REFPROP___ButaneOctane(comp1, comp2, kij):
this_function_name = sys._getframe().f_code.co_name
print('\n---\n THE DATA USED IN THIS SIMULATION WAS OBTAINED FROM: ',
this_function_name)
print('\n---\n')
'''
SOURCE: REFPROP -- BUTANE AND OCTANE.
'''
critical_pressure = np.array([comp1.pC, comp2.pC])
critical_temperature = np.array([comp1.Tc, comp2.Tc])
acentric_factor = np.array([comp1.AcF, comp2.AcF])
molar_mass = np.array([comp1.MM, comp2.MM])
omega_a = 0.45724 * np.ones_like(molar_mass)
omega_b = 0.07780 * np.ones_like(molar_mass)
binary_interaction = np.array([[0.0000, kij], [kij, 0.0000]])
specific_heat_mass_base = np.array([comp1.Cp, comp2.Cp])
specific_heat = Tools_Convert.convert_specific_heat_massbase_TO_molarbase(
specific_heat_mass_base, molar_mass)
return (critical_pressure, critical_temperature, acentric_factor,
molar_mass, omega_a, omega_b, binary_interaction, specific_heat)
=============================================================================================================
Class Molecule("name", MM, Tc, pC, AcF, Cp)
LEGEND:
name: compost name
MM: molar mass of each component [kg/kmol]
Tc: critical temperature [K]
pC: critical pressure [Pa]
AcF: acentric factor [-]
Cp: specific heat [J / kg K] !<=== HERE IT IS IN MASS BASE, BUT IT IS CONVERTED TO MOLAR BASE INSIDE THE FUNCTION
=============================================================================================================
'''
MM = np.array([58.12, 240.])
specific_heat_mass_base = np.array([1663.0, 1490.0]) #(J/kg K)
specific_heat = Tools_Convert.convert_specific_heat_massbase_TO_molarbase(
specific_heat_mass_base, MM)
#sort as name, molar_mass, Tc, pC, AcF, Cp
comp1 = Molecule("R600a", 58.12, (134.7 + 273.15), 36.4e5, 0.1853,
specific_heat[0])
comp2 = Molecule("AB_ISO5", 240., (675.9 + 273.15), 20.6e5, 0.9012,
specific_heat[1])
kij = -0.02668
'''
===============================================================================================================
Do you have kij above? If you don't have it ==> uncomment the CODE's 4 lines below
===============================================================================================================
'''
# p_C = np.array([comp1.pC, comp2.pC])
# T_C = np.array([comp1.TC, comp2.TC])
# kij_obj = Kij_class(p_C, T_C)