c2h6 <- Gas("ethane")
# estimate LJ-coefficients from viscosity data
ethane_viscosity <- subset(ethane_data,
(property=="viscosity") & (T>=90) & (T<=675), select=c(T, value)
)
c2h6$zeta <- 0
ethane_viscosity$value <- 1E-6*ethane_viscosity$value
c2h6$fit_viscosity_data(ethane_viscosity)
print(c2h6$sigma)
print(c2h6$epsk)
# estimate LJ-coefficients from second virial coefficient
ethane_B <- subset(ethane_data, (property=="BQFH") & (T>= 220) & (T<=623))
ethane_B$value <- 1E-6*ethane_B$value
c2h6$fit_B_data(ethane_B)
print(c2h6$sigma)
print(c2h6$epsk)
Run the code above in your browser using DataLab