Learn R Programming

chapensk (version 0.4)

ethane_data: Viscosity and second virial coefficient of ethane

Description

Reference values for the second virial coefficient and viscosity of ethane from a intermolecular potential energy surface. The second virial coefficient was calculated semiclassically by means of the Mayer-sampling Monte Carlo technique, while the transport properties were obtained using the classical kinetic theory of polyatomic gases. The computed thermophysical property values are in excellent agreement with the best available experimental data and are recommended as reference values.

Usage

data("ethane_data")

Arguments

Format

A data frame with 107 observations on the following 4 variables.

T

Temperature in K

property

Type of property. B: classically calculated second virial coefficient, BQFH: second virial coefficient calculated by a modification of the pair potential known as the quadratic Feynman-Hibbs (QFH) effective pair potential. viscosity: gas phase viscosity

unit

a factor with levels cm3/mol, uPa.s

value

numerical value of property. For the second virial coefficient the calculated data are supported by experimental data in the temperature range 220 - 623 K. For the viscosity the calculated data are supported by experimental data in the temperature range 90 - 675 K.

Details

Type of fitσε\(/k\)
vs. viscosity data4.38(2)235(6)
vs. second virial coefficient data4.95(6)202(3)
vs. viscosity and second virial coefficient data4.35244

Examples

Run this code
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