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_rm_rk_delta_1(MODEL_eos, OM, dinputs): Tc, Pc, OM, Vceos = dinputs[0], dinputs[1], dinputs[2], dinputs[3] if MODEL_eos == 'SRK': rm = 0.48 + 1.574 * OM - 0.175 * OM**2 delta_1 = 1.0 return rm, delta_1 elif MODEL_eos == 'PR': rm = 0.37464 + 1.54226 * OM - 0.26992 * OM ** 2 delta_1 = 1.0 + np.sqrt(2.0) return rm, delta_1 elif MODEL_eos == "RKPR": Zc = Pc * Vceos / (RGAS * Tc) del1_init = func_delta_init(Zc) delta_1 = getdel1(Zc, del1_init)[0] # 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) return rk, delta_1
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 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
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]
RT = RGAS * Tc
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 * RT**2 / Pc
b = OMb * RT / Pc
params = [ac, b, rm]
if ICALC == 'parameters_eps':
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 (Tc,Pc,OM,Vceos)
Tc = dinputs[0]
Pc = dinputs[1]
OM = dinputs[2]
Vceos = dinputs[3]
RT = RGAS * Tc
Zc = Pc * Vceos / RT
del1ini = D[0] + D[1] * (D[2] - Zc)**D[3] + D[4] * (D[2] -
Zc)**D[5]
print('del1ini = {0}'.format(del1ini))
del_1 = getdel1(Zc, del1ini)[0]
Zc, OMa, OMb = compressibility_factor_cal(del_1)
print('Zc = {0}'.format(Zc))
ac = OMa * RT**2 / Pc
b = OMb * RT / Pc
# calcular rk
rk = 1
params = [ac, b, del_1, rk]
elif ICALC == 'parameters_eps':
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]
# RKPR EOS (Tc,Pc,OM)
elif ICALC == 'rk_param':
# dinputs = np.array([Tc, Pc, OM, dc, zrat, ac, b, d, rk])
Tc = dinputs[0]
Pc = dinputs[1]
OM = dinputs[2]
dc = dinputs[3]
zrat = dinputs[4]
ac = dinputs[5]
b = dinputs[6]
d = dinputs[7]
delta_1 = d
rk = dinputs[8]
omega = OM
#RT = RGAS * Tc
#print(dinputs)
# THEN ! del1 SPECIFICATION together with Tc,Pc,OM
# READ(NIN,*)del1 ! this line must be active when it is not a list
#rk = dinputs[0]
#delta_1 = dinputs[1]
#print(rk, delta_1, omega)
rk, Pvdat, Tr = initial_data(omega, delta_1, NMODEL, ICALC, Pc)
#print(Tr)
eos_calculation = Parameter_eos()
rk_cal = eos_calculation.resolver_rk_cal(rk, delta_1, Pvdat, Pc,
Tc, Tr)
#rk_cal = Parameter_eos.resolver_rk_cal(rk, delta_1, Pvdat, Pc, Tc, Tr)
#print('rk_cal = '.format(rk_cal))
# RhoLsat SPECIFICATION together with Tc,Pc,OM
elif ICALC == 'density':
# (T, RhoLsat)
# Trho = T / Tc
# del1 = 2.0 #! initial value
# RHOld = 0.0
Tc = dinputs[0]
Pc = dinputs[1]
omega = dinputs[2]
delta_1 = dinputs[3]
T_especific = dinputs[4]
RHOLSat_esp = dinputs[5]
rk, Pvdat, Tr = initial_data(omega, delta_1, NMODEL, ICALC, Pc)
eos_calculation = Parameter_eos()
delta_1_parameter = eos_calculation.resolver_delta_1_cal(
delta_1, rk, Pvdat, RHOLSat_esp, Pc, Tc, Tr)
#print(delta_1_parameter)
print('The NMODEL is eos_{0} and method ICALC is {1}'.format(
NMODEL, ICALC))
if ICALC == 'constants_eps':
print(params)
return params
elif ICALC == 'parameters_eps':
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) is {0}'.format(delta_1_parameter))
return delta_1_parameter