Modelación matemática del proceso de digestión anaerobia en condiciones de clima frio utilizando
biodigestores tubulares

Mathematical modeling of the anaerobic digestion process under cold weather conditions using tubular digesters

Juan José Sossa¹, René Alvarez^1
1Instituto de Investigación y Desarrollo de Procesos Químicos (IIDEPROQ) UMSA, Campus Universitario Calle N° 30, zona de Cota Cota, La Paz, Bolivia. sljj91@gmail.com ]]> Artículo recibido en: 7 de diciembre 2015 Manejado por: Comité Editorial Aceptado en: 11 de marzo de 2016

Resumen

El propósito principal de esta investigación es modelar y simular matemáticamente el proceso de digestión anaerobia en condiciones de clima frío utilizando biodigestores tubulares, considerando la variación de la temperatura ambiental, el comportamiento hidrodinámico del biodigestor, el mecanismo de reacción, la estequiometría y la cinética de la digestión anaerobia. En el modelo matemático que se propone se ha tomado en cuenta el comportamiento térmico del biodigestor a través de la obtención de correlaciones matemáticas que relacionan el tiempo con la temperatura ambiental y la temperatura de la fase gaseosa del biodigestor, además se asume que la temperatura en la fase líquida se mantiene constante. Por otro lado, se considera que el comportamiento hidrodinámico del biodigestor está representado por cuatro reactores CSTR conectados en serie. Adicionalmente el modelo ADM1 ha sido aplicado para representar el mecanismo de reacción, la estequiometría y la cinética del proceso de digestión anaerobia. El modelo matemático fue calibrado y validado utilizando información experimental recopilada durante el período de trabajo. Como resultado se obtuvo un conjunto de parámetros cinéticos que podrán ser utilizados en posteriores aplicaciones de ADM1 en el caso de la digestión anaerobia de estiércol de vaca. Los resultados de la simulación presentan un grado aceptable de concordancia con los resultados experimentales y las diferencias se deben a variaciones en el comportamiento hidrodinámico del biodigestor.

Palabras clave: Digestión anaerobia, estiércol de vaca, biodigestortubular, modelación matemática, ADM1.

Abstract

The main purpose of this research is to mathematically model and simulate mathematically the process of anaerobic digestion in cold weather conditions using tubular digesters, considering the variation of the air temperature, the hydrodynamic behavior of the digester, the reaction mechanism, stoichiometry and kinetics anaerobic digestion. In the proposed mathematical model the thermal behavior of the digester is taken into account through obtaining mathematical correlations relating time with air temperature and the temperature of the gaseous phase of the digester, in addition it is assumed thatthe temperature in the liquid phase remains constant. Furthermore, it is considered that the hydrodynamic behavior of the digester is represented by four CSTR reactors connected in series. Additionally the ADM1 model has been applied to represent the reaction mechanism, stoichiometry and kinetics of anaerobic digestion process. The mathematical model was calibrated and validated using experimental data collected during the work period. As a result, a set of kinetic parameters, was obtained as they could be used in subsequent applications of ADM1 for anaerobic digestion of cow manure. The simulation results show an acceptable degree of agreement with the experimental results and the differences are due to variations in the hydrodynamic behavior of the digester.

Keywords: Anaerobic digestion, cow manure, tubular digester, mathematical modeling.

INTRODUCCIÓN

La digestión anaerobia es un proceso dependiente de la temperatura donde los fenómenos de transferencia de materia, energía y cantidad de movimiento se ven influenciados porlas condiciones atmosféricas de su entorno y por la configuración del biodigestor en el cual se desarrolla. Es por ello que surge la motivación de estudiar y evaluar los efectos que tiene la temperatura sobre todo el proceso de digestión anaerobia. La modelación y simulación matemática de procesos es una disciplina creada para el desarrollo de un mejor nivel de entendimiento de las partes de un sistema y de todo el sistema como conjunto (Jana 2011). Por lo tanto el abordaje del problema de evaluación de los efectos de la temperatura sobre el proceso de digestión anaerobia se realizará a través de la modelación y simulación matemática.

En el campo de la modelación y simulación matemática de la digestión anaerobia de estiércol de vaca se han desarrollado modelos matemáticos enfocados en distintas partes del proceso. Algunos consideran el mecanismo de reacción, la estequiometría y la cinética, a una determinada temperatura constante (Hill & Barth, 1977; Hill D, 1982; Vavilin et al., 1994; Angelidaki et al., 1999). Existen otros que fueron desarrollados para predecir el comportamiento térmico del biodigestor (Axaopoulos et al., 2001;Perrigault et al., 2012) e incluso llegan a predecir la producción de biogás (Terradas-Ill et al., 2014). Sin embargo aún persiste la necesidad de estudiar todo el proceso en su conjunto. De esta forma surge la intención de desarrollar un modelo matemático que considere los principales aspectos del proceso de digestión anaerobia en reactores tubulares instalados en regiones de clima frío como ser el comportamiento hidrodinámico y térmico del biodigestor, además del mecanismo de reacción, la estequiometría y la cinética del proceso de digestión anaerobia.

Por este motivo el objetivo principal del presente trabajo de investigación es el de desarrollar un modelo matemático capaz de predecir el pH dentro el biodigestor, la reducción de materia orgánica, el flujo de biogás y su composición, considerando la variación de temperatura existente en zonas con climas predominantemente fríos. Como resultado se tendrá una herramienta eficiente para la evaluación de futuras instalaciones de digestores en diferentes zonas de clima predominantemente frío, utilizando la configuración de digestor que se utilizó en el presente trabajo.

MATERIALES Y MÉTODOS

Sistema experimental

El sistema experimental está constituido por un biodigestor tubular instalado a 1 m bajo el nivel del suelo, en una zanja excavada, dentro de un invernadero. Para evitarque el calorse disipe desde el biodigestor hacia el suelo, se aisló el biodigestor con poliestireno expandido (plastoformo), el cual se coloca en el perímetro de la zanja. El biodigestor fue diseñado siguiendo la metodología propuesta por Martí-Herrero & Cipriano (2012), el período de instalación del biodigestor y la construcción del invernadero tomó, aproximadamente, tres semanas. La Figura 1, muestra las dimensiones (en centímetros)del sistema experimental instalado y, en la Tabla 1 se presenta una descripción detallada de las características del biodigestor.

El período de experimentación duró 5 meses (mayo - septiembre de 2015), durante este tiempo el biodigestor fue alimentado con una mezcla de estiércol de vaca y agua, en una relación 1:3, y se trabajó con un tiempo de retención hidráulica (TRH) de 84 d. En este período se monitoreo el sistema, midiendo mensualmente la reducción de materiaorgánica en términos de DQO, el pH a la salida del biodigestor, el flujo de biogás, su concentración en metano y las temperaturas en distintos puntos del sistema experimental (Figura 1b): (i) ambiente, (ii) dentro el invernadero, (iii) en la fase gaseosa y (iv) en la fase líquida del biodigestor.

Métodos analíticos

En los análisis de DQO y pH de las muestras tomadas mensualmente a la entrada y la salida del biodigestor se aplicaron los Standard Methods establecidos por la Asociación Americana de Salud Pública (American Public Health Association APHA, 2005). Para medir el flujo de biogás se utilizó un medidor de flujo de gas natural UGG 2,5 Metrix y la concentración de metano en el biogás fue monitoreada con un analizador de gases GA-21 Plus (MADUR) provisto de un sensor infrarrojo para la medición de metano. En el seguimiento del comportamiento térmico del biodigestor se utilizaron sensores de temperatura HOBO y PT-100, el intervalo de tiempo para el registro de cada dato de temperatura fue de 30 minutos.

Planteamiento del modelo matemático

La tarea de modelación matemática del proceso de digestión anaerobia en el reactor tubular fue realizada tomando en cuenta los principales aspectos del proceso: (a) el comportamiento térmico del sistema experimental, (b) el mecanismo de reacción, la estequiometría y la cinética del proceso de digestión anaerobia y (c) el comportamiento hidrodinámico del biodigestor.

Comportamiento térmico del sistema experimental

Con la información recolectada se obtuvieron funciones matemáticas para establecer una relación entre la temperatura ambiente, el tiempo y la temperatura del invernadero y de las fases gaseosa y líquida del biodigestor, aplicando el método de ajuste de mínimos cuadrados, cuya función objetivo se muestra en la Ecuación 1:

Siendo Z, la función objetivo, Ti_(T) (t), el valor de la temperatura ambiente para un tiempo (t) determinado por la función matemática ajustada; Ti_(t), la temperatura experimental para un tiempo (t) y N, el número total de mediciones experimentales.

Mecanismo de reacción

El mecanismo de reacción utilizado para la representación del proceso de digestión anaerobia fue el propuesto porel ADM1 (Batstone etal., 2002). La estequiometría propuesta en el ADM1 desde su publicación original en 2002 ha ido mejorando constantemente, es por ello que los coeficientes estequiométricos utilizados en el presente trabajo fueron tomados de Rosen & Jeppsson (2006). Asi mismo, los parámetros cinéticos utilizados en la descripción de la cinética del proceso de digestión anaerobia, la caracterización del sustrato y las condicioens iniciales del biodigestor fueron tomados de Page et al.(2008) y Rosen & Jeppsson (2006), a su vez, estos sirvieron para estimar el valor de los parámetros a la temperatura de trabajo.

Efecto de la temperatura sobre la cinética del proceso de digestión anaerobia

La temperatura tiene influencia sobre los procesos fisicoquímicos y bioquímicos dentro la digestión anaerobia, sin embargo existen diferentes expresiones matemáticas derivadas de un análisis termodinámico para expresar esta dependencia. Los procesos fisicoquímicos (equilibrio ácido-base y transferencia de materia líquido-gas) presentan una dependencia respecto a la temperatura de acuerdo a la ecuación de Van't Hoff (Batstone et al., 2002):

Aplicando la Ecuación (1) a los coeficientes cinéticos de los procesos fisicoquímicos se tiene:

Siendo K_a (T) (mol L-1), la constante de disociación ácido-base a una determinada temperatura T (K) ; K_a,0 (mol L^-1), la constante de disociación ácido-base a la temperatura de referencia T_o (K), AH°, el cambio de entalpía en condiciones estándar (J); R, la constante universal de los gases ideales (J mol¹ K^-1); K_H(T) (mol bar^-1), la constante de Henry una determinada temperatura T(K) y K_H,₀ (mol bar^-1), la constante de Henry a la temperatura de referencia T_o (K).

Respecto a los procesos bioquímicos se tiene que la expresión cinética del proceso de crecimiento celulary de consumo de sustrato está representado por la ecuación de Monod:

Siendo j, el subíndice que identifica a cada uno de los procesos de crecimiento celulary de consumo de sustrato;py, la velocidad de reacción (kg DQO m^-3 d^-1) k_m,_j, la velocidad máxima de crecimiento específica de los microorganismos (d^-1); K_s,_j, la constante de velocidad de crecimiento medio de los microorganismos (kg DQO m^-3); Ij, la función de inhibición de cada proceso.

Los coeficientes cinéticos de los procesos de reacción química presentan una dependencia respecto a la temperatura que se representa de acuerdo a la ecuación de Arrhenius:

Siendo k(T), el coeficiente cinético a una temperatura T (°C)k₂₀, el coeficiente a 20°C y 0, el coeficiente de actividad-temperatura.La mayoría de los estudios realizados respecto a la influencia de la temperatura sobre los procesos biológicos tan solo reportan el efecto total sobre el proceso global y no así sobre cada uno de los procesos individuales desarrollados por cada grupo de microorganismos. En el presente trabajo se utilizarán los coeficientes de actividad-temperatura tomados de Grady et al. (2011): 0_km =1.105 y 0_Ks =0.819, donde 0_km y 0_Ks son los coeficientes de actividad-temperatura utilizados para representar el efecto de la temperatura sobre k_m,_j y Ks,j, respectivamente.Finalmente, en el modelo matemático no se tomó en cuenta el aporte energético de las reacciones bioquímicas involucradas en el proceso de digestión anaerobia puesto que se considera que la mayoría de la energía generada en el proceso de desintegración e hidrólisis es utilizada en la producción de acetato y posterior generación de metano, este hecho puede verse reflejado en el análisis termodinámico del proceso de digestión anaerobia desarrollado por Lübken et al., (2007), dónde se demuestra que la energía remanente tiene un aporte muy bajo.

Comportamiento hidrodinámico del biodigestor

El comportamiento hidrodinámico de un reactor tubular puede ser representado a partir de la conexión en serie de varios reactores tanque con agitación continua ideales (CSTR) (Zaher et al., 2006). En este trabajo se asume que el comportamiento hidrodinámico del biodigestor está representado por cuatro CSTR conectados en serie (Pageet al., 2008), de esta forma los balances másicos se realizaron de acuerdo a este comportamiento. Una representación esquemática del biodigestor se muestra en la Figura 2, dónde los flujos de los componentes particulados y solubles son transportados de una etapa a la otra y la sumatoria de los flujos de biogás se constituye en la producción total de biogás.

Estructura del modelo matemático

Finalmente, la estructura del modelo matemático se muestra en el diagrama de flujo de información de la Figura 3. El modelo matemático fue resuelto en el programa computacional de cálculo numérico MatLab.

Durante la calibración del modelo matemático se realizó un análisis de sensibilidad, para ello se utilizaron los datos experimentales recolectados durante el período mayo-julio y se calcularon las funciones de sensibilidad absoluta relativa (Jeppson, 1996):

Donde 5_(r,r), es el cambio relativo en una variable de salida x por un 100% en el cambio relativo del parámetro 0. La variable de salida es la reducción de materia orgánica y se utilizaron los parámetros cinéticos de los procesos bioquímicos de la digestión anaerobia. A continuación, se ajustaron los resultados del modelo con los experimentales, modificando los parámetros cinéticos, utilizándose como función objetivo la función de mínimos cuadrados (Donoso-Bravo et al., 2011),

Siendo J la función objetivo, v_exp las mediciones experimentales de la variable seleccionada tomadas durante el período de experimentación, v_sim los valores de la variable de salida seleccionada obtenidos mediante el modelo, 0 el valor de los parámetros que deben ser determinados y N el número de mediciones. Para evaluar la precisión y aplicabilidad de los parámetros calibrados, en la validación del modelo matemático se compararon los resultados simuladosy los datos experimentales correspondientes al período experimental agosto-septiembre.

Comportamiento térmico del biodigestor

Debido a que existe un comportamiento oscilante en las temperaturas del ambiente, el invernadero y la fase gaseosa del biodigestor, fue necesario ajustar los datos de la temperatura ambiente de cada mes a una función sinusoidal. Por otro lado la temperatura ambiente fue correlacionada con la temperatura del invernadero del biodigestor a través de una función exponencial, que fue la que mejor ajuste ofrecía. La temperatura de la fase gaseosa del biodigestor fue relacionada con la temperatura del invernadero a través de una función lineal, finalmente la temperatura de la fase líquida y de la fase gaseosa se encuentra relacionada a través de otra función lineal. La Figura 4 muestra un ejemplo del grado de ajuste entre las correlaciones matemáticas obtenidas y los datos experimentales de temperatura para el mes de junio. La Tabla 2 presenta los coeficientes hallados para las funciones de ajuste de temperatura y los coeficientes de correlación.

De acuerdo a la Tabla 2, se tiene que T_a es la temperatura ambiental (°C); Ti, la temperatura del invernadero (°C), T_fg, de la fase gaseosa del biodigestor; T_fl, la temperatura de la fase líquida del biodigestor y t, el tiempo (d). Para el modelo matemático se utilizaron las constantes promedio calculadas mensualmente para las funciones de ajuste de t vs. T_a, T_a vs. Ti y Ti vs. T_fg ya que existe un aceptable grado de concordancia entre las correlaciones matemáticas ajustadas y los datos experimentales. Por otro lado, como se puede observar en la Figura 4d y la Tabla 2, no existe correlación alguna entre la temperatura de la fase líquida y de la fase gaseosa, puesto que el coeficiente de correlación es muy bajo, por lo tanto en el modelo matemático se utilizó la temperatura promedio de la fase líquida para todo el período de experimentación (20 °C).

Calibración del modelo matemático

La primera tarea en el proceso de calibración del modelo matemático consiste en realizar un análisis de sensibilidad y como resultado se obtienen los parámetros más sensibles del proceso (Tabla 3). La estimación de parámetros se realizó con el objetivo de conseguir un mejor ajuste entre los resultados experimentales y los del modelo matemático, modificando los parámetros más sensibles del proceso. Para este propósito se utilizaron los datos experimentales de reducción de materia orgánica en términos de DQO, pertenecientes al período experimental comprendido entre mayo y julio.

Los parámetros más sensibles son las constantes de velocidad mitad a 20°C para los microorganismos consumidores de valerato y butirato, acetato e hidrógeno. Esto indica que la etapa crítica del mecanismo de reacción es el consumo de los ácidos grasos volátiles e hidrógeno para la producción de metano, (Tabla 3) este hecho también ha sido reportado por Page et al. (2008).

Para Parker (2005) los valores menores de ofrecen resultados más apropiados. Finalmente, el conjunto de parámetros cinéticos, dependientes de la temperatura, que deben aplicarse para la modelación y simulación matemática del proceso de digestión anaerobia de estiércol de vaca en digestores tubulares se muestra en la Tabla 4.

Validación del modelo matemático

El modelo matemático fue validado a partir de la comparación de los resultados obtenidos a través del modelo y los resultados experimentales. La Figura 5, presenta los resultados para: (a) pH a la salida del biodigestor; (b) reducción de materia orgánica (M.O.) en términos de reducción de DQO; (c) Flujo promedio mensual de biogás; (d) Composición promedio mensual del biogás. Respecto al comportamiento del pH a la salida del biodigestor se puede observar una buena correlación de los resultados del modelo respecto a los datos experimentales. Las predicciones del modelo sobre la reducción de materia orgánica en términos de reducción de DQO en el biodigestor para los meses de mayo a junio, no muestran un buen ajuste, esto puede haber sido ocasionado porque el comportamiento hidrodinámico del biodigestor no obedecia exactamente al de 4 CSTR conectados en serie, este hecho puede deberse a la formación de zonas muertas o corto circuitos que incrementan o disminuyen el tiempo de retención hidráulica dentro el biodigestor, disminuyendo e incrementando, respectivamente, la reducción de DQO. Sin embargo para el período agosto-septiembre las predicciones del modelo se ajustan más a los datos experimentales debido a que las predicciones ingresan dentro el rango del error de las mediciones. Respecto a las predicciones del flujo promedio de biogás producido en el biodigestor, se observa que existe una directa relación con la temperatura ambiente, en los meses de menor temperatura se registraron menores flujos de producción de biogás, y en los meses de mayor temperatura, la producción de biogás se incrementa de forma directamente proporcional, además se puede observar que en los meses donde hubo una mayor reducción de materia orgánica se producieron menores volúmenes de biogás. Los resultados predecidos por el modelo no presentan un buen ajuste con los resultados experimentales y estimados obtenidos en el período mayo-septiembre, siendo mayores los resultados del modelo. Los flujos de biogás estimados, fueron calculados en base a la reducción de materia orgánica, considerando que el 10% de la DQO es materia orgánica inerte. Este problema pudo también haber sido originado por las diferencias del comportamiento hidrodinámico del biodigestor, siendo que a mayor tiempo de retención hidráulica, al menos por encima de los 25 d, se tiene una menor producción específica de biogás. La concentración de metano en el biogás presenta una buen grado de concordancia existiendo una diferencia máxima de 8% entre los resultados del modelo y los experimentales para el mes de mayo y de 5% para el mes de julio, el resto de los meses el ajuste es aceptable.

En la Figura 6, se puede observar que el comportamiento del flujo de biogás a través del tiempo es oscilante y está directamente relacionado con la temperatura. El mayor flujo de biogás se da cuando la temperatura de la fase gaseosa del biodigestor alcanza su valor máximo. Esto demuestra que la producción de biogás no es constante durante todo el día.

CONCLUSIONES

Se realizó la modelación y simulación matemática del proceso de digestión anaerobia en un biodigestor tubular instalado en una región de clima predominantemente frío. Se ha desarrollado un modelo matemático que considera el comportamiento térmico e hidrodinámico del biodigestor. Además gracias a la implementación del ADM1, el mecanismo de reacción, la estequiometría y la cinética de la digestión anaerobia fueron incluidas. El modelo matemático fue calibrado utilizando resultados experimentales y fruto de ello se obtuvieron un conjunto de parámetros cinéticos bioquímicos que servirán para la posterior aplicación del ADM1 en el diseño de biodigestores que utilicen como sustrato estiércol de vaca. Finalmente, las simulaciones presentan un grado aceptable de concordancia respecto a las mediciones experimentales y las diferencias existentes se deben a una variabilidad en el comportamiento hidrodinámico del biodigestor.

AGRADECIMIENTOS

Se agradece la colaboración de la Fundación HIVOS que financió la investigación, al equipo técnico de HIVOS que se encargó de la construcción y monitoreo de la temperatura del sistema experimental. A su vez, también se agradece la colaboración del personal del IIAREN, y la Estación Experimental Choquenaira por realizar el mantenimiento y operación del biodigestor bajo estudio.

APHA; AWWA and WEF. 2005. Standard Methods for the examination of water & wastewater. Washington, DC, U.S.A.:APHA.        [ Links ]

Angelidaki, I, L. Ellegaard and B. K.Ahring. 1999. A comprehensive model of anaerobic bioconversion of complex substrates to biogas. Biotechnology and bioengineering. Vol. 63 (3): 363-372.        [ Links ]

Axaopoulos, P., P. Panagakis, A. Tsavdaris and D. Georgakakis. 2001. Simulation and experimental performance of a solar-heated anaerobic digester.. Solar Energy. Vol. 70 (2): 155-164.        [ Links ]

Batstone, D. J., Keller, J., Angelidaki, I., Kalyuzhnyi, S. V., Pavlostathis, S. G., Rozzi, A., Sanders, W. T. M., Siegrist, H. and V. A. Vavilin. 2002. The IWA Anaerobic Digestion Model No 1 (ADM 1). London, UK: IWApublishing        [ Links ]

Characklis, W. G. and W. Gujer. 1978. Temperature dependency of microbial reactions (Eds. Jenkins, S.H.) Kinetics of Wastewater Treatment: Proceedings of a Post-Conference Seminar Held at theTechnical University of Denmark, Copenhagen. Technical University of Denmark. Copenhagen, Denmark. pp. 111-130        [ Links ]

Donoso-Bravo,A.,J.Mailier, C.Mailier, J. Rodríguez, C. A.Aceves-Lara &Wouwe, A. V. 2011. Model selection, identification and validation in anaerobic digestion: A review. Water Research. Vol. 45(1): 5347-5364.        [ Links ]

Grady, C. P. L., G. T. Daigger, N. G. Love and C.D.M. Filipe. 2011 .Biological Wastewater Treatment. Florida, USA: CRC Press        [ Links ]

Hill, D. T. and C. L. Barth. 1977. A dynamic model for simulation of animal waste digestion. Journal (Water Pollution Control Federation). Vol. 49 (10): 2129-2143.        [ Links ]

Hill, D.T. 1982. A comprehensive dynamic model for animal waste methanogenesis..Transactions of the ASAF (American Society of Agricultural Engineers). Vol. 25 (): 1374-1380.        [ Links ]

Jana, A. K. 2011. Chemical Process Modelling and Computer Simulation. Delhi: PHI Learning Private Limited        [ Links ]

Jeppson, U. 1996. ModellingAspectsofWastewater Treatment Processes.Doctoral Dissertation.Lund, Sweden: Lund University.        [ Links ]

Lübken, M., M. Wichern, M. Schlattmann, A. Gronauer and H. Horn. 2007. Modelling the energy balance of an anaerobic digester fed with cattle manure and renewable energy crops. Water research. Vol. 41 (18): 4085-4096.        [ Links ]

Martí-Herrero, J. and J. Cipriano. 2012. Design methodology for low cost tubular digesters. Bioresource Technology. Vol. 108 (): 21-27.        [ Links ]

Page, D. I., K. L. Hickey, R. Naruta, A. L. Main and S. J. Grimberg. 2008. Modeling anaerobic digestion of diary manure using the IWA Anerobic Digestion Model No. 1 (ADM1). Water Science & Technology. Vol. 58 (3): 689-695.        [ Links ]

Parker, W. J. 2005. Application of theADM1 model to advanced anaerobic digestion. Bioresource Technology. Vol. 96 (1): 1832-1842.        [ Links ]

Perrigault, T., V., Martí-Herrero, J. Weatherford and D. Poggio. 2012. Towards thermal design optimization of tubular digesters in cold climates: A heat transfer model. Bioresource technology. Vol. 124 ():259-268.        [ Links ]

Rosen, C. and U. Jeppsson. 2006.Aspects on ADM1 Implementation within BSM2 Framework. Lund, Sweden: Lund University        [ Links ]

Tabatabaei, M., A. Sulaiman, A. M. Nikbakht, G. Najafpour and N. Yusof. 2011 Influential parameters on biomethane generation in anaerobic wastewatertreatment plants (Eds. Manzanera, M.) Alternative Fuel.INTECH Open Access Publisher. Rijeka, Croatia. pp. 227-262        [ Links ]

Terradas-Ill, G., C. H. Pham, J. M. Triolo, J. Martí-Herrero and S. G. Sommer. 2014. Thermic model to predict biogas production in unheated fixed-dome digesters buried in the ground. Environmental science & technology. Vol. 48 (6): 3253-2362.        [ Links ]

Vavilin, V. A., V. B. Vasiliev, A. V. Ponomarev and S. V. Rytow. 1994. Simulation model 'methane'as a tool for effective biogas production during anaerobic conversion of complex organic matter. Bioresource technology. Vol. 48 (1): 1-8.        [ Links ]

Zaher, U. and S. Chen. 2006. Interfacing the IWA Anaerobic Digestion Model No. 1 (ADM1) With Manure and Solid Waste Characteristics. Proceedings of theWaterEnvironmentFederation. Vol. 2006(9): 3162-3175.        [ Links ]