Shell correlation from Kumar thesis (2005)
z.Shell(pres.pr, temp.pr, tolerance = 1e-13, verbose = FALSE)
pres.pr | pseudo-reduced pressure |
---|---|
temp.pr | pseudo-reduced temperature |
tolerance | controls the iteration accuracy |
verbose | print internal |
# single z point and create a dataframe with info ppr <- 1.5 tpr <- 1.1 z.calc <- z.Shell(pres.pr = ppr, temp.pr = tpr) # From the Standing-Katz chart we obtain a digitized point: z.chart <- getStandingKatzMatrix(tpr_vector = tpr, pprRange = "lp")[1, as.character(ppr)] ape <- abs((z.calc - z.chart) / z.chart) * 100 df <- as.data.frame(list(Ppr = ppr, z.calc =z.calc, z.chart = z.chart, ape=ape)) rownames(df) <- tpr df#> Ppr z.calc z.chart ape #> 1.1 1.5 0.4869976 0.426 14.31868