tests/gas_class.R

In chapensk: Estimation of Gas Properties from the Lennard-Jones Potential

library(chapensk)

CH4 <- Gas("methane")
Ar <- Gas("argon")
# virial coefficient
stopifnot(round(1E6 * CH4$B(T = 300), 0) == -41)
# viscosity
stopifnot(round(1E6 * CH4$viscosity(T = 300), 1) == 11.2)
# density
stopifnot(round(CH4$density(T = 300), 3) == 0.652)
# diffusion coefficient
stopifnot(round(1E4 * CH4$diffusion(T = 300), 3) == 0.199)
# thermal conductivity
stopifnot(round(Ar$thermal_conductivity(T = 300), 3) == 0.018)
# polar molecule
HCl <- Gas("hydrogen chloride")
HCl$viscosity(T = 300)