Kinetic Model for Anaerobic Digestion of Distillery Spent Wash
Pooja Sharma^{*}, Uttam Kumar Ghosh, Amiya Kumar Ray
Department of Polymer and Process Engineering, Indian Institute of Technology Roorkee, Roorkee, India
Email address:
To cite this article:
Pooja Sharma, Uttam Kumar Ghosh, Amiya Kumar Ray. Kinetic Model for Anaerobic Digestion of Distillery Spent Wash. American Journal of Chemical Engineering. Vol. 4, No. 6, 2016, pp. 139-146. doi: 10.11648/j.ajche.20160406.11
Received: July 21, 2016; Accepted: August 2, 2016; Published: November 15, 2016
Abstract: Kinetic studies were carried out for anaerobic digestion of distillery spent wash considering growth kinetic models due to Monod and Contois. The kinetic parameters were predicted with the help of previously published experimental data for anaerobically treated distillery spent wash. The digester was operated at seven different values of loading rates ranging between 0.954 to 2.459 gCOD/L/day. Maximum specific growth rate and maximum specific substrate utilization rate for Monod model have been predicted as of the order of 0.153 g/g/day and 3.06 g/g/day , respectively, and those for Contois kinetic model were found to be 0.140 g/g/day and 0.421 g/g/day , respectively. However, for both the models the decay rate coefficient, was calculated as 0.0131 g/g/day. The estimated value of empirical constant in Contois model was found to be order of 0.699 (dimensionless). The deviation of model predicted values from experimental results has been compared. Sensitivity analysis of the developed models was carried out by considering 10% increase in the values of both and .
Keywords: Kinetic Modelling, Anaerobic Digester, Distillery Spent Wash, Monod Kinetics, Contois Kinetics and Sensitivity Analysis
1. Introduction
Anaerobic digestion process is used to treat wastewater generated from various municipal and industrial activities and also sludge generated from aerobic processes [1,2]. This process is further classified into suspended growth (contact process, sequencing batch reactor, continuous stirred tank reactor, etc.) and attached growth (upflow packed bed reactor, upflow anaerobic expanded bed reactor, fluidized bed reactor, downflow attached growth anaerobic processes, etc.) biological wastewater treatment processes. Various microbial growth kinetic models such as Monod, Contois, Chen and Hashimoto, Haldane, Mosey, etc., have been developed by various investigators to describe the performance of anaerobic digester [3]. These expressions are required in developing total continuity equations for describing the performance of anaerobic digester and to predict behavior under dynamic conditions. Monod kinetic model, which is based on Michaelis-Menten equation, gives relationship between specific growth rate and growth limiting substrate in dissolved form [3,4] whereas, Contois model is a functional relationship between specific growth rate and concentrations of substrate both at the inlet and the outlet of the digester [3,5]. Chen and Hashimoto kinetic model is concerned with substrate utilization and methane production [5].
Distillery wastewater generated from distillation of the alcohol (rectified spirit from analyzer column) which in fact, vary in composition depending on the kind of raw materials and type of fermentation process. In India fermentation is accomplished with cane molasses using yeast cell and Saccharomyces cerevisiae. The distillery spent wash is the waste generated from distillery based on cane molasses as raw material. The data reported on kinetic parameters of anaerobic treatment of distillery wastes including vinasses (waste generated after fermentation of molasses) are plenty [6-10]. However, microbial kinetic parameters required for design applications and estimation of treatment efficiency for digestion of spent wash based on cane molasses as raw material are scanty. Not any studies reported till date which describes the microbial growth kinetic coefficients for distillery waste based on Contois model for design applications. Large volume of data related to anaerobic degradation kinetics of sewage sludge and other wastes are available in literature [7,8].
The predicted values of the microbial growth kinetic parameters for different kind of wastes vary from waste to waste and depend on the source of waste generation, process conditions, the digester used for treatment and nature of inhibitory compounds present in the wastes. Study carried out for microbial growth kinetics for different types of wastes in complete mixed anaerobic digester operated at mesophilic temperature range. The kinetic parameters were estimated on the basis of mathematical modelling (by developing steady state model) for mixed community of microorganisms based on Monod kinetic model. These are discussed in the following paragraphs:
Bello-Mendoza et al. [1] reported the kinetic coefficient values for sewage sludge consist of maximum specific substrate utilization rate (8.0 g/g/day), yield coefficient (0.3 g/g), decay rate coefficients for acidogens (0.1 day^{-1}) and for methanogens (0.015 day^{-1}). Hu et al. [5] determined the kinetic coefficients for anaerobic contact digester treating ice-cream wastewater which comprised of yield coefficient (0.2116 g/g), empirical constant based on Contois model (0.4818) and half velocity coefficient (0.4028 g/L). The specific bacterial growth rates based on Monod as well as Contois model were found (0.7844 day^{-1}) and (0.9207 day^{-1}), respectively. For cattle waste the above kinetic coefficient was found to be order of 0.37 L/g, 0.30 g/L and 0.25 day^{-1}, respectively. Borja et al. [8] estimated that the yield coefficient of microbial species degrading brewery wastewater is 0.08 g/g and specific substrate utilization rate 0.045 g/g/day. Based on Chen-Hashimoto model, Jimenez et al. [6] investigated on kinetics for untreated and previously treated vinasses where maximum production rate of methane was found to be 1.27 L/day and 0.891 L/day, respectively. The specific growth rate and half velocity coefficient for untreated waste were 0.87 day^{-1} and 0.90 g/L, respectively, and 0.09 day^{-1} and 0.13 g/L, respectively, for treated vinasses.
The aim of the present study is to predict the kinetic coefficient for distillery spent wash using the data reported by Borja et al. [11] for the developed models based on Monod and Contois kinetics and to compare them to describe the performance of complete mix anaerobic digester.
2. Model Development
Steady state model has been developed for the estimation of the kinetic parameters for distillery spent wash based on sugar cane molasses as raw material. The experimental data was published in literature due to Borja et al. [11] was undertaken. The experiments were conducted under mesophilic conditions in continuous complete mix digester of volume 2.35 L operated at seven different loading rates ranging from 0.954 g/L /day to 2.4589 g/L/day. The influent substrate concentration fed to digester in the above experiment had COD level of 4.7 g/L. The laboratory data for spent wash could give the process efficiency in terms of COD removal and rate of methane production. This data fails to describe the inside conditions of the digester such as microbial kinetics, the growth and decay rate terms, specific growth rate and utilization rate of substrate by micro-organisms.
Mass balance equations have been applied across complete mix suspended growth anaerobic digester to derive a new set of equations based on organic loading rate variation. The microbial kinetics has been described by two steady-state models due to Monod as well as Contois. The developed kinetic rate expressions are based on the concepts proposed earlier by Tchobanoglous et al. [4]; Hu et al. [5].
2.1. Kinetic Modelling of Anaerobic Digester Based on Monod Kinetics
2.1.1. Mass Balance Equation for Spent Wash
(1)
Under steady-state condition,
Substrate utilization rate based on Michaelis-Menten equation,
Therefore, Eq. (1) may be written as
(2)
Based on Monod growth kinetics
and
and Eq. (2) can be written as follow
(3)
Now with the help of laboratory data for spent wash and by linear regression method the kinetic parameters and can be evaluated from Eq. (3).
2.1.2. Mass Balance Equation for Biomass
Inflow – Outflow ± Net growth rate= accumulation
(4)
Under steady-state conditions
and assuming the concentration of biomass in the raw spent wash feed to the digester is negligible, Eq. (4) becomes
(5)
Use of Monod kinetic model yields
(6)
With the help of Eq. (6) and also from Eq. (2), the developed equation can be written as
(7)
With the help of laboratory data kinetic parameters can be predicted by linear regression
2.2. Kinetic Modelling of Anaerobic Digester Based on Contois Kinetics
Now, from Eq. (6) using Contois model [3] under steady state conditions the equation can be written as
(8)
From Eq. (8)
(9)
Laboratory data can be used to predict kinetic parameters by regression analysis using Eq. (9).
3. Results and Discussion Based on Distillery Spent Wash
3.1. Evaluation of Kinetic Parameters for Spent Wash
The kinetic parameters for Monod as well as those from Contois model for distillery spent wash have been estimated through the developed model as given in Section 2 and also with the distillery data mentioned in previously published research work due to Borja et al. [11]. The developed model (Eq. (3), Eq. (7) and Eq. (9)) are used to evaluate the microbial kinetic rates.
Figure 1 has been plotted between as a function of for the estimation of kinetic parameters and from Eq. (3) for spent wash based on Monod kinetics. From the slope of the straight line and the intercept the values of and , respectively, can be predicted.
Figure 1. Estimation of saturation constant, , and maximum specific substrate utilization rate, , for distillery spent wash based on Monod growth kinetic model.
Figure 2 has been drawn against for the estimation of and as given in Eq. (7). The values of and can be evaluated by least square fit.
Figure 2. Estimation of decay rate coefficient,, and synthesis yield coefficient, , for distillery spent wash based on Monod growth kinetic model.
For the prediction of kinetic parameters and based on Contois growth kinetics from Eq. (9), plot versus as shown in Figure 3. With the help of linear regression method the values of and are predicted.
Figure 3. Estimation of maximum specific substrate utilization rate, , and empirical constant,, for distillery spent wash based on Contois kinetic model.
The estimated kinetic parameters for spent wash based on model equations are compared with the predicted values for different kinds of wastes [4,5,12-15] as mentioned in Table 1. The evaluated kinetic parameters as given in Table 1 are in close agreement with the range of microbial kinetic parameters based on modelling for different kind of wastes.
Kinetic Coefficient | Distillery Spent wash (present study) | Ice-cream wastewater Hu et al. [5] | Distillery wastewater Prakash et al. [12] | Municipal solid waste Nwabanne et al. [13] | Sewage sludge Tchobanoglous et al. [4] |
(g/g/day) | Monod 0.153 | Monod 0.7844 | 1.56 | 0.53 | 0.35 (0.30-0.38) |
Contois 0.140 | Contois 0.9297 | - | - | - | |
(g/g) | 0.3328 | 0.2116 | 0.5 | 0.367 | 0.08 (0.05-0.1) |
(g/g/day) | 0.0131 | 0.0131 | 0.05 | 0.038 | 0.03 (0.002-0.04) |
(g/L) | 0.4605 | 0.4028 | 0.024 | 0.21 | 0.16 (0.06-0.2) |
(g/g/day) | 3.06 | 3.7 | 3.125 | 0.144 | - |
(g/g/day) | 0.421 | 4.39 | - | - | - |
(Dimensionless) | 0.699 | 0.4818 |
3.2. Deviation of the Model Predicted Values of, , from Experimental Data
The formulae used for prediction of effluent substrate concentration () based on Monod growth kinetic model was obtained by rearranging Eq. (3) and Eq. (9), respectively.
The deviation of model predicted effluent substrate concentration () from the experimental values due to Borja et al. [8] for both Monod and Contois models has been plotted with the help of estimated kinetic parameters. These are shown in Figure 4 based on Monod model whereas Figure 5 is plotted for Contois model. From Figures 4 and Figure 5, the R^{2} values are found of the order of 0.8751 for Monod whereas the same for Contois model 0.9323. The percentage deviation for both the models is found below 12.5%. However, it is observed that Contois model fits the experimental data better for spent wash, indicating its superiority over the other models. Hence, model developed based on Contois model can suitability be used for spent wash.
Figure 4. Comparison between the model predicted values of based on Monod model with the experimental values for distillery spent wash based on Monod kinetics.
Figure 5. Comparison between model predicted values of with the experimental results for distillery spent wash based on Contois kinetic model.
3.3. Sensitivity Analysis for Spent Wash
Sensitivity analysis is essential for analyzing accuracy of the data whether it is experimental or model predicted. Therefore, attempts have been made to carry out the sensitivity analysis by varying the values of maximum specific growth rate coefficient,, and decay rate coefficient,, for both the models (Monod and Contois) and their effects on COD concentration in effluent (). With the help of formulae (based on Eq. (3) and Eq. (9)), the values for for both the models are estimated. Figure 6 through Figure 9 represent the model predicted values as function of HRT for Monod model and organic loading rate (OLR) for Contois model. The values in Contois model are predicted by plotting it against OLR instead of HRT because the effluent substrate concentration in Contois model is dependent upon the influent substrate concentration feed to the digester.
Figure 6. Sensitivity analysis of value and its effect on effluent substrate concentration, as function of HRT for distillery spent wash based on Monod model.
Figure 7. Sensitivity analysis of and its influence on effluent substrate concentration, as function of OLR for distillery spent wash based on Contois model.
The sensitivity of microbial growth kinetic term on process performance in terms of COD removal has been predicted by considering 10% variation in value for both the models. If it is assumed that the process is operated under steady state conditions, there is not any inhibitory effect noticed on microbial growth kinetics [16-18]. At steady state conditions, based on Monod model 10% increase in the value of as shown in Figure 6 is found to decrease with the COD in effluent resulting in higher efficiency. From the Figure 6 it is also observed that if the growth kinetics is assumed to be inhibited under any shock loading condition (10% decrease in ). Under this situation, lower COD removal efficiency is obtained as compared to other two cases (10% increase in and at constant value of ). At the hydraulic retention time of 50 day, COD removal percentage is found to be 94.88%. Increasing by 10% increase of percent COD removal from 94.88% to 96.02% has been observed. Further decreasing by 10% decreases in percentage COD removal from 94.88% to 94.18%. At the retention time of 20.6 day, increase in value of by 10%, COD removal increases from 80.15% to 86.51%, whereas by decreasing by 10%, reduction in percent COD removal occurred from 80.15% to 62.78%.
The same condition has been applied for Contois model and is plotted in Figure 7. At the OLR of 1.606 g/L/day, increasing value by 10% increases percent COD removal from 91.39% to 92.98%, whereas decrease in by 10% decreases the same from 91.39% to 88.91%. From this study it can be concluded that if the microbial growth rate decreases there is significant effect on process performance in terms of COD removal. Again simulation studies have been carried out considering the combined effect of both as well as on COD removal by varying 10% of the values as shown in Figure 8 and Figure 9.
Figure 8 has been plotted between as a function of HRT whereas Figure 9 as a function of OLR. It is well known fact that as the growth rate of microbial species increases there is also simultaneous increase in decay rate of microorganisms. By considering these aspects process performance has been predicted by considering combined effect of both the variables. As shown in Figure 8 based on Monod model if both the kinetic parameters increases by 10%, the percent COD removal increases from 88.61% to 90.46% at the retention time of 25.7 day. Reduction in percent COD removal from 88.61% to 85.13% has been observed at the same retention time on decreasing both the variables by 10%. Figure 9 represents the variation in terms of as function of OLR based on Contois model. If the digester is operated at OLR of 1.856 g/L/day, the combined effect of both the parameters by increasing their values at 10% resulted in increase in percentage COD removal from 88.79% to 90.63%. Further, it is also noticed a reduction in percentage from 90.63% to 86.13% by decreasing both the parameters by 10%.
By performing sensitivity analysis, the importance of the process variables and its effect on process performance can be evaluated. Also methods are to be devised for enhanced performance by controlling the inhibitory substances that effect microbial survival.
Figure 8. Sensitivity analysis of and , and their combined effect on as function of for distillery spent wash based on Monod model.
Figure 9. Sensitivity analysis of and , and their combined effect on as function of for distillery spent wash based on Contois model.
Table 2 represents the confidence interval level of the values predicted from model equations. Student t-distribution method has been adopted for the above estimation [19, 20]. Three different probability levels, 90%, 95% and 99% were selected for estimation of values for both Monod as well as Contois model.
Table 2. Confidence interval for from derived equations for distillery spent wash based on Monod as well as Contois growth kinetic models.
Effluent substrate concentration (g/L) | , | , | , | ||
Monod kinetic model | |||||
(90%) | 6.44±2.5 | 4.57±1.51 | 10.45±5.06 | 5.08±1.72 | 9.47±4.69 |
(95%) | 6.44±3.17 | 4.57±1.90 | 10.45±6.38 | 5.08±2.17 | 9.47±5.89 |
(99%) | 6.44±4.80 | 4.57±2.88 | 10.45±9.66 | 5.08±3.29 | 9.47±8.94 |
Contois kinetic model | |||||
(90%) | 5.92±2.23 | 4.57±1.56 | 8.46±3.68 | 4.81±1.67 | 7.67±3.21 |
(95%) | 5.92±2.81 | 4.57±1.98 | 8.46±4.64 | 4.81±2.11 | 7.67±4.04 |
(99%) | 5.92±4.25 | 4.57±2.99 | 8.46±7.03 | 4.81±3.19 | 7.67±6.12 |
4. Conclusions
From the above analysis the following noteworthy conclusions can be drawn:
• Steady state model equation based on Contois model is found more suitable for analyzing performance of anaerobic digestion of spent wash.
• For spent wash the developed steady state model equations fits well for analyzing performance of anaerobic digestion. Statistical analysis reveals that Contois model is found to be more suitable (with higher R^{2} value of 0.9323) in this type of waste. Confidence interval for at 95% for both Monod as well as Contois based models was found 6.44±3.17 and 5.92±2.81, respectively.
• Sensitivity analysis for spent wash also explains well the impact of kinetic parameters on process performance. Decrease in COD removal efficiency has been observed at lower specific growth rate of microorganisms.
• The predicted values of maximum specific substrate utilization rate for both Monod as well as Contois models for distillery spent wash is lower as compared to that of vinasses. This indicates the presence of inhibitory substances in distillery spent wash that inhibit growth kinetics of microbial species.
Acknowledgement
This work has been supported by Ministry of Human Resources and Development (MHRD), Government of India.
Nomenclature
| Influent substrate concentration (g/L) |
| Effluent substrate concentration (g/L) |
| Substrate utilization rate (g/L/day) |
| Synthesis yield coefficient (g/g) |
| Maximum specific substrate utilization rate for Monod kinetics (g/g/day) |
| Maximum specific substrate utilization rate for Contois Kinetics (g/g/day) |
| Half velocity coefficient for Monod growth kinetics (g/L) |
| Maximum specific bacterial growth rate (g/g/day) |
| Specific bacterial growth rate (g/g/day) |
| Contois kinetic constant (Dimensionless) |
| Influent biomass concentration (g/L) |
| Effluent biomass concentration (g/L) |
| Net biomass growth rate (g/L/day) |
| Decay coefficient (g/g/day) |
| Organic loading rate (g COD/L/day) |
| Hydraulic retention time (day) |
Q | Volumetric flow rate (L/day) |
V | Volume of digester (L) |
References