The biogas package is an R package for biogas research. It provides
tools for processing biogas data, predicting biogas production, making
conversions, and planning experiments. For example, the
calcBgVol() function can be used for calculating
biochemical methane potential from original measurements. The
calcCOD() function can calculate the chemical oxygen demand
of a compound based on a chemical formula. And the
planBMP() function can suggest substrate or inoculum
quantities for a biochemical methane potential (BMP) experiment based on
multiple constraints.
Much of the biogas package is dedicated to processing of laboratory
measurements, but those functions are difficult to demonstrate in a
short example. So, an alternative is the predBg() microbial
stoichiometry function. With it, the theoretical maximum BMP of, say,
cellulose, could be calculated with
> predBg("C6H10O5")
[1] 413.7274and we could get a lot more information, including an estimate of biogas composition as affected by pH, for a 5:1 mixture of waste paper and waste oil with
> predBg(mcomp = c(C6H10O5 = 5/6, C54H100O7 = 1/6), mass = 1,
fd = 0.8, fs = 0.1, conc.sub = 50, pH = c(6.5, 7, 8.5),
temp = 35, value = "all")
                      form mass mol.mass      moles      COD  fs  fe  fd
1 C1.52618H2.61516O1.00000    1  36.9655 0.02705225 1.454303 0.1 0.9 0.8
2 C1.52618H2.61516O1.00000    1  36.9655 0.02705225 1.454303 0.1 0.9 0.8
3 C1.52618H2.61516O1.00000    1  36.9655 0.02705225 1.454303 0.1 0.9 0.8
  conc.sub temp  pH      hydro      fCH4      xCH4     vCH4      vCO2      vBg
1       50   35 6.5 0.09279292 0.5566161 0.5689900 365.8395 275.91411 641.7536
2       50   35 7.0 0.09279292 0.5566161 0.5864017 365.8395 256.90636 622.7459
3       50   35 8.5 0.09279292 0.5566161 0.8921576 365.8395  44.02904 409.8686
       mCH4      mCO2     mCO2Bg   mCO2.sol       cTIC      m.bio      N.req
1 0.2624617 0.5735654 0.54543300 0.02813242 0.03196140 0.08225028 0.01018521
2 0.2624617 0.5735654 0.50785807 0.06570735 0.07465048 0.08225028 0.01018521
3 0.2624617 0.5735654 0.08703757 0.48652785 0.55274693 0.08225028 0.01018521See the vignettes and help files for more details.
The biogas package is available from CRAN, and so can be installed directly in R with:
install.packages("biogas")To see available vignettes, use
vignette(package = 'biogas')and then, e.g.,
vignette('calcBgVol_function')or
vignette('predBg_function')Alternatively, download the vignettes from the CRAN page.