def func_constans(MODEL_eos, dinputs):
if MODEL_eos == "SRK" and MODEL_eos == "PR":
# SRK and PR
ac, b, rm = dinputs[0], dinputs[1], dinputs[2]
if MODEL_eos == "SRK":
del1, al, be, ga = 1.0, -0.175, 1.574, 0.48 - rm
elif MODEL_eos == "PR":
del1, al, be, ga = 1.0 + np.sqrt(2.0), -0.26992, 1.54226, 0.37464 - rm
Zc, OMa, OMb = compressibility_factor_cal(del1)
elif MODEL_eos == "RKPR":
# RKPR
ac, b, del1, rk = dinputs[0], dinputs[1], dinputs[2], dinputs[3]
Zc, OMa, OMb = compressibility_factor_cal(del1)
if MODEL_eos == "RKPR":
al, be, ga = A1 * Zc + A0, B1 * Zc + B0, C1 * Zc + C0 - rk
Tc = (OMb * ac) / (OMa * RGAS * b)
Pc = OMb * RGAS * Tc / b
OM = acentric_factor_cal(al, be, ga)
Vceos = Zc * RGAS * Tc / Pc
return Tc, Pc, OM, Vceos, Zc
def initial_data(omega, delta_1, MODEL_eos, SPECIFICATION_cal, Pc, dinputs):
Zc, OMa, OMb = compressibility_factor_cal(delta_1)
# initial guess for k parameter
rk = (A1 * Zc + A0) * omega ** 2 + (B1 * Zc + B0) * omega + (C1 * Zc + C0)
# rk = rk * 1.2 # 1.1 #5.2 #3.2
argumen_1 = "constants_eps"
argumen_2 = "parameters_eps"
argumen_3 = "rk_param"
argumen_4 = "density"
condition_1 = SPECIFICATION_cal == argumen_1
condition_2 = SPECIFICATION_cal == argumen_2
condition_3 = SPECIFICATION_cal == argumen_3
condition_4 = SPECIFICATION_cal == argumen_4
if condition_1 or condition_2 or condition_3:
rk *= 1.5
Tr = 0.7
Pvdat = Pc * 10 ** -(1.0 + omega)
elif condition_4:
# 5.2 es otro valor que se puede usar en lugar de 1.5
rk = rk * 1.5
Tr_calculada = dinputs[4] / dinputs[0]
Tr = Tr_calculada
Pvdat = Pc * 10 ** -((1.0 / Tr - 1.0) * 7 * (1.0 + omega) / 3)
return rk, Pvdat, Tr
def call_rkpr_constans_density(MODEL_eos, SPECIFICATION_cal, dinputs):
# SPECIFICATION [Tc, Pc, OM, del1, Tspec, RHOLsat]
# Trho = Tspec / Tc, read initial value of del1
Tc = dinputs[0]
Pc = dinputs[1]
OM = dinputs[2]
delta_1 = dinputs[3]
T_especific = dinputs[4]
RHOLSat_esp = dinputs[5]
Zc, OMa, OMb = compressibility_factor_cal(delta_1)
print('Zc = {0}'.format(Zc))
ac, b = func_ac_b(Tc, Pc, Zc, OMa, OMb)
rk, Pvdat, Tr = initial_data(OM, delta_1, MODEL_eos, SPECIFICATION_cal, Pc, dinputs)
eos_calculation = Parameter_eos()
list_args = [delta_1, rk, Pvdat, RHOLSat_esp, Pc, Tc, Tr]
delta_1_parameter = eos_calculation.resolver_delta_1_cal(list_args)
parameters = [ac, b, rk, delta_1_parameter]
return parameters
def call_rkpr_constans_v_critic(MODEL_eos, SPECIFICATION_cal, dinputs):
# CONSTANTS SPECIFICATION READ [Tc, Pc, OM, Vceos]
Tc = dinputs[0]
Pc = dinputs[1]
OM = dinputs[2]
Vceos = dinputs[3]
Zc = Pc * Vceos / (RGAS * Tc)
del1_init = func_delta_init(Zc)
delta_1 = getdel1(Zc, del1_init)[0]
Zc, OMa, OMb = compressibility_factor_cal(delta_1)
print('Zc = {0}'.format(Zc))
ac, b = func_ac_b(Tc, Pc, Zc, OMa, OMb)
# calcular rk
rk, Pvdat, Tr = initial_data(OM, delta_1, MODEL_eos, SPECIFICATION_cal, Pc, dinputs)
eos_calculation = Parameter_eos()
rk_cal = eos_calculation.resolver_rk_cal(rk, delta_1, Pvdat, Pc, Tc, Tr)
# rk = 1
parameters = [ac, b, rk_cal, delta_1]
return parameters
def func_zc_ac_b(Tc, Pc, Zc, delta_1):
Zc, OMa, OMb = compressibility_factor_cal(delta_1)
# Vceos = (Zc * RGAS * Tc) / Pc
ac, b = func_ac_b(Tc, Pc, OMa, OMb)
return Zc, ac, b
def parameter_ab_cal(self, delta_1_initial, Pc, Tc):
Zc, OMa, OMb = compressibility_factor_cal(delta_1_initial)
dc = Pc / Zc / (RGAS * Tc)
Vceos = 1.0 / dc
self.ac = OMa * (RGAS * Tc)**2 / Pc
self.b = OMb * (RGAS * Tc) / Pc
return self.ac, self.b
def initial_data(omega, delta_1, NMODEL, ICALC, Pc, dinputs):
Zc, OMa, OMb = compressibility_factor_cal(delta_1)
# initial guess for k parameter
rk = (A1 * Zc + A0) * omega**2 + (B1 * Zc + B0) * omega + (C1 * Zc + C0)
#rk = rk * 1.2 # 1.1 #5.2 #3.2
if ICALC == 'constants_eps' or ICALC == 'parameters_eps' or ICALC == 'rk_param':
rk *= 1.5
Tr = 0.7
Pvdat = Pc * 10 ** -(1.0 + omega)
elif ICALC == 'density':
rk = rk * 1.5 #5.2
Tr_calculada = dinputs[4] / dinputs[0]
Tr = Tr_calculada
Pvdat = Pc * 10 ** -((1.0 / Tr - 1.0) * 7 * (1.0 + omega) / 3)
return rk, Pvdat, Tr
def call_rkpr_constans_delta_1(MODEL_eos, SPECIFICATION_cal, dinputs):
# SPECIFICATION [Tc, Pc, OM, del1]
Tc = dinputs[0]
Pc = dinputs[1]
OM = dinputs[2]
delta_1 = dinputs[7]
Zc, OMa, OMb = compressibility_factor_cal(delta_1)
print('Zc = {0}'.format(Zc))
ac, b = func_ac_b(Tc, Pc, Zc, OMa, OMb)
rk, Pvdat, Tr = initial_data(OM, delta_1, MODEL_eos, SPECIFICATION_cal, Pc, dinputs)
eos_calculation = Parameter_eos()
rk_cal = eos_calculation.resolver_rk_cal(rk, delta_1, Pvdat, Pc, Tc, Tr)
parameters = [ac, b, rk_cal, delta_1]
return parameters
def call_rkpr_parameters(MODEL_eos, SPECIFICATION_cal, dinputs):
# SPECIFICATION [ac, b, rk, del1]
ac = dinputs[0]
b = dinputs[1]
del1 = dinputs[2]
rk = dinputs[3]
Zc, OMa, OMb = compressibility_factor_cal(del1)
Tc = OMb * ac / (OMa * RGAS * b)
Pc = OMb * RGAS * Tc / b
Vceos = Zc * RGAS * Tc / Pc
al, be, ga = A1 * Zc + A0, B1 * Zc + B0, C1 * Zc + C0 - rk
OM = acentric_factor_cal(al, be, ga)
# Tc, Pc, OM, Vceos, Zc = func_constans(del1, ac, b, al, be, ga)
constants = [Tc, Pc, OM, Vceos]
return constants
def models_eos_cal(NMODEL, ICALC, dinputs):
if NMODEL == 'SRK' or NMODEL == 'PR':
# CONSTANTS SPECIFICATION READ [Tc, Pc, OM]
if ICALC == 'constants_eps':
Tc = dinputs[0]
Pc = dinputs[1]
OM = dinputs[2]
if NMODEL == 'SRK':
rm = 0.48 + 1.574 * OM - 0.175 * OM**2
del1 = 1.0
elif NMODEL == 'PR':
rm = 0.37464 + 1.54226 * OM - 0.26992 * OM ** 2
del1 = 1.0 + np.sqrt(2.0)
Zc, OMa, OMb = compressibility_factor_cal(del1)
Vceos = (Zc * RGAS * Tc) / Pc
ac = OMa * (RGAS * Tc) ** 2 / Pc
b = OMb * (RGAS * Tc) / Pc
params = [ac, b, rm, del1]
# PARAMETERS SPECIFICATION READ [ac, b, rm]
if ICALC == 'parameters_eps':
ac = dinputs[0]
b = dinputs[1]
rm = dinputs[2]
Tc = (OMb * ac) / (OMa * RGAS * b)
Pc = OMb * RGAS * Tc / b
Vceos = Zc * RGAS * Tc / Pc
if NMODEL == 'SRK':
del1 = 1.0
al = -0.175
be = 1.574
ga = 0.48 - rm
elif NMODEL == 'PR':
del1 = 1.0 + np.sqrt(2.0)
al = -0.26992
be = 1.54226
ga = 0.37464 - rm
OM = acentric_factor_cal(al, be, ga)
constants = [Tc, Pc, OM, Vceos]
elif NMODEL == 'RKPR':
if ICALC == 'constants_eps':
# CONSTANTS SPECIFICATION READ [Tc, Pc, OM, Vceos]
Tc = dinputs[0]
Pc = dinputs[1]
OM = dinputs[2]
Vceos = dinputs[3]
Zc = Pc * Vceos / (RGAS * Tc)
del1ini = D[0] + D[1] * (D[2] - Zc) ** D[3] + D[4] * (D[2] - Zc)** D[5]
print('del1ini = {0}'.format(del1ini))
delta_1 = getdel1(Zc, del1ini)[0]
Zc, OMa, OMb = compressibility_factor_cal(delta_1)
print('Zc = {0}'.format(Zc))
ac = OMa * (RGAS * Tc) ** 2 / Pc
b = OMb * (RGAS * Tc) / Pc
# calcular rk
rk, Pvdat, Tr = initial_data(OM, delta_1, NMODEL, ICALC, Pc, dinputs)
eos_calculation = Parameter_eos()
rk_cal = eos_calculation.resolver_rk_cal(rk, delta_1, Pvdat, Pc, Tc, Tr)
# rk = 1
params = [ac, b, rk, delta_1]
elif ICALC == 'parameters_eps':
# PARAMETERS SPECIFICATION READ [ac, b, rk, del1]
ac = dinputs[0]
b = dinputs[1]
del1 = dinputs[2]
rk = dinputs[3]
Zc, OMa, OMb = compressibility_factor_cal(del1)
Tc = OMb * ac / (OMa * RGAS * b)
Pc = OMb * RGAS * Tc / b
Vceos = Zc * RGAS * Tc / Pc
al = A1 * Zc + A0
be = B1 * Zc + B0
ga = C1 * Zc + C0 - rk
OM = acentric_factor_cal(al, be, ga)
constants = [Tc, Pc, OM, Vceos]
elif ICALC == 'rk_param':
# CONSTANTS SPECIFICATION and del1 READ [Tc, Pc, OM, del1]
Tc = dinputs[0]
Pc = dinputs[1]
OM = dinputs[2]
delta_1 = dinputs[7]
rk, Pvdat, Tr = initial_data(OM, delta_1, NMODEL, ICALC, Pc, dinputs)
eos_calculation = Parameter_eos()
rk_cal = eos_calculation.resolver_rk_cal(rk, delta_1, Pvdat, Pc, Tc, Tr)
elif ICALC == 'density':
# CONSTANTS SPECIFICATION and (T, RhoLsat) READ [Tc, Pc, OM, del1, T, RHOLsat]
# Trho = T / Tc, read initial value of del1
Tc = dinputs[0]
Pc = dinputs[1]
OM = dinputs[2]
delta_1 = dinputs[3]
T_especific = dinputs[4]
RHOLSat_esp = dinputs[5]
rk, Pvdat, Tr = initial_data(OM, delta_1, NMODEL, ICALC, Pc, dinputs)
eos_calculation = Parameter_eos()
delta_1_parameter = eos_calculation.resolver_delta_1_cal(delta_1, rk, Pvdat, RHOLSat_esp, Pc, Tc, Tr)
print('The NMODEL is eos_{0} and method ICALC is {1}'.format(NMODEL, ICALC))
if ICALC == 'constants_eps':
print("params = [ac, b, rm, del1]")
#ac, b, rm, del1 = params
#print("ac = {0} b = {1} rm = {2} del1 = {3}".format(ac, b, rm, del1))
return params
elif ICALC == 'parameters_eps':
print("constants = [Tc, Pc, OM, Vceos]")
print(constants)
return constants
elif ICALC == 'rk_param':
print('The parameter rk_cal is {0}'.format(rk_cal))
return rk_cal
elif ICALC == 'density':
print('The parameter delta1(rho,T) = {0}'.format(delta_1_parameter))
return delta_1_parameter