Next Article in Journal
Oily Wastewater Treatment: Methods, Challenges, and Trends
Previous Article in Journal
Influence of Material Microstructure on Machining Characteristics of OFHC Copper C102 in Orthogonal Micro-Turning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Model of a Thermophilic Anaerobic Digestion for Methane Production of Wheat Straw

1
The Stephan Angeloff Institute of Microbiology, Bulgarian Academy of Sciences, Acad. G. Bonchev Str., Bl. 26, 1113 Sofia, Bulgaria
2
Institute of Biophysics and Biomedical Engineering, Bulgarian Academy of Sciences, Acad. G. Bonchev Str., Bl. 105, 1113 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
Processes 2022, 10(4), 742; https://doi.org/10.3390/pr10040742
Submission received: 21 March 2022 / Revised: 5 April 2022 / Accepted: 6 April 2022 / Published: 12 April 2022
(This article belongs to the Section Process Control and Monitoring)

Abstract

:
This paper presents a newly created mathematical model of thermophilic anaerobic digestion of wheat straw carried out in a 2 dm3 bioreactor for methane production. Two batch processes, with 30 mL/dm3 and 35 mL/dm3 organic load, are carried out—one set for parameter identification and one set for model verification. The identification of model parameter values is based on dynamical experiments. It is fulfilled using two different techniques: deterministic sequential quadratic programming algorithm and metaheuristic genetic algorithm. Verification of the developed mathematical models is conducted based on the different data sets of the process. Both models predict the set of the experimental data for all considered process variables well. Genetic algorithm visually fits the data with a higher degree of accuracy, as confirmed by the numerical results for the objective function value.

1. Introduction

Fossil fuels consumption is increasing and has become a major contributor to the increasing concentration of carbon dioxide in the atmosphere, leading to global warming [1]. Attention currently is paid to diverse renewable energy sources; although the development and use of renewable energy is not high yet, it is expected to increase as fossil fuel supplies decline. Various technologies are projected to provide renewable energy production, renewable biomass utilisation being one of them [2]. The development and utilisation of biomass as an energy source is still very small compared to its potential. However, lignocellulosic biomass is a promising alternative for meeting the global sustainability criteria [3]. Biogas is a renewable gas, mainly consisting of methane and carbon dioxide that can be used to produce heat and electricity. Biogas can also be utilised directly as a vehicle’s fuel after an upgrade to biomethane. Increased production of biogas in the world will favour energy supply needs [4]. Utilising agricultural wastes adds another advantage to their use for the production of biofuels [5].
Wheat straw is considered a favourable substrate for biogas production. Lignocellulosic biomass is one of the most abundant among various renewable alternatives. However, due to its complex and rigid structure, its biodegradability during anaerobic digestion is usually low. Challenges associated with the substrate pretreatment are very important for the overall process performance because this is an expensive step that affects the efficiency of all subsequent steps of sustainable valorisation of the biomass components [6]. Temperature and substrate composition are among the main factors affecting the realisation and stability of anaerobic digestion processes. Pretreatments are mainly carried out at high temperatures and appear effective in realising lignin removal. Cellulose and hemicellulose could be hydrolysed into monomeric sugars such as glucose, xylose, and arabinose, which could then be converted to biofuels such as methane [7]. The sought results are in the direction of the appearance of structural changes and improved enzymatic hydrolysis [8]. Liquid hot water pretreatment was studied to explore the feasibility of improving the methane yield of wheat straw in anaerobic digestion, and the results showed that the crosslinking structure of wheat straw had been broken [9].
Improved biomethane yield was also obtained after thermal pretreatment [10]. Conveying the process of anaerobic digestion of wheat straw under thermophilic conditions could strongly affect the performance of biogas digesters, leading to increased hydrolysis rates. Temperature accelerates biochemical processes, enabling higher degradation rates and biogas yields from a wide variety of substrates compared to mesophilic anaerobic digestion processes [11].
The operating physicochemical parameters of anaerobic digesters and the methods of acceleration and optimisation used to improve biogas yields remain of utmost importance [12].
The complicated model structure of anaerobic digestion processes is typically described using nonlinear differential equations. The mathematical modelling of such processes is a complex problem, and the selection of a proper method for optimisation is fundamental for the accurate estimation of model parameters.
There are a few techniques that are appropriate in the considered case problem. Among them are conventional optimisation methods, such as sequential quadratic programming (SQP) [13,14], and metaheuristic techniques, such as genetic algorithms (GA) [15,16]. As a good opportunity for conventional optimisation methods, metaheuristic techniques have been recognised. They are capable of finding a suitable solution for acceptable computational time [17,18,19]. Some of the most powerful metaheuristic algorithms inspired by nature—such as genetic algorithms, artificial bee colony algorithm, Cuckoo search algorithm, etc.—have been proposed and validated for many optimisation problems and have been endorsed as prosperous for the estimation of bioprocess models parameters [20,21,22].
The paper strives to propose a mathematical model of thermophilic anaerobic digestion of wheat straw carried out in a 2 dm3 bioreactor for methane production. For the purposes of modelling, two batch processes with 30 mL/dm3 and 35 mL/dm3 organic load, are carried out—one set for parameter identification and one set for model verification.
The numerous successful GA applications in the field of modelling bioprocesses have inspired the authors to adopt and employ GA for the first time for solving the challenging problem of parameter identification of nonlinear model of the considered anaerobic digestion process. The obtained numerical estimates of model parameters are compared to the results from application of SQP—one of the prominent techniques for solving constrained nonlinear optimisation problems. Both obtained mathematical models are verified with an independent data set.
The paper is organised as follows. After the Introduction, experimental setup and mathematical modelling problem are provided in Section 2. The experimental results are presented and discussed in Section 3. Concluding comments and future research aspects are summarised in Section 4.

2. Materials and Methods

2.1. Experimental Setup

In this study, as substrate for the anaerobic digestion, mechanically pretreated wheat straw was used. For this purpose, native wheat straw was ground to a final particle size of 1–2 mm using a hammer mill, followed by additional milling in a knife mill until a satisfactory particle size was reached.
The experiments were conducted in a bioreactor with a working volume of 2 dm3 for methane production. Initial experiments, as well as inoculum maintenance, were conducted in a bioreactor of 1 dm3. In all experiments, a temperature suitable for growing thermophilic microorganisms (55 °C) was maintained automatically, and continuous stirring was also applied. All experiments were carried out in batch mode of operation. The daily released biogas was kept and measured using a gas holder based on the water displacement method.
Cellulose content was determined by a spectrophotometric method [23]. Cellulose-containing materials were cleared of impurities such as lignin, hemicellulose, xylosans, and other low molecular weight compounds by extraction with an acetate-nitrite mixture. Purified cellulose was dissolved in 67% H2SO4, followed by a colour reaction using an anthrone reagent. Cellulose concentration was determined after measuring the absorbance at 620 nm. Glucose was analysed using the Miller method [24], which is based on the redox reaction between reducing sugars and sodium dinitrosalicylate, resulting in a reddish-brown derivative. Absorption was measured at a wavelength of 530 nm. The concentration of acetate was determined by a Thermo Scientific gas chromatograph (Focus GC model) equipped with a Split/Splitless injector, column: TG-WAXMS A (length 30 m, diameter 0.25 mm, film thickness 0.25 μm), and flame ionisation detector (FID). Prior to injection, the pH of the sample taken from the bioreactor was adjusted to pH 2.0 with 37% H3PO4. After one hour, the sample was centrifuged at 15,000× g min−1 for 10 min, and an aliquot of the supernatant was mixed with an equal volume of 1% 2,2-dimethyl-butyric acid (as internal standard).
The content of the biogas produced from the bioreactor (containing methane) was estimated with a “Dräger, X-am7000” model device (Germany), equipped with infrared sensors for measuring the relative content of CH4 and CO2 in percent by volume and a catalytic sensor for H2S (in ppm).

2.2. Mathematical Modelling Problem

In order to find a set of design parameters, x = {x1, x2, ..., xn}, optimisation techniques are used. The resulting outcomes can be defined as optimal ones. Usually, a minimisation or maximisation of a predefined system characteristic, dependent on x, is considered. In a more complex formulation, the objective function f(x), which should be minimised (maximised), might be subject to constraints. The problem can be described as follows:
min x J ( x )
where x is the vector of design parameters, J(x) is the objective function.
The discrepancy between the model and real data is configured as an optimisation criterion:
J = i = 1 k j = 1 m ( Y i j Y ^ i j ) 2
where k is the number of process variables, m is the number of data points, Y i j represents the observed experimental data values, and Y ^ i j represents model predicted values.
The obtained solution depends mainly on the problem complexity, i.e., the number of design variables and constraints. Solution efficiency also depends on constraints characteristics and the objective function.

2.3. Sequential Quadratic Programming

Quadratic programming (QP) considers the minimisation (maximisation) of a quadratic objective function with linear constraints. Trustworthy solution techniques are readily accessible. Constrained minimisation is the problem of finding a vector x that is a local minimum to a scalar function J(x) subject to constraints on the allowable x. The method mimics Newton’s method for constrained optimisation, just as is carried out for unconstrained optimisation. A quasi-Newton updating method is used at each iteration in order to approximate the Hessian of the Lagrangian function. A QP subproblem solution, needed to form a search direction, is then generated. The background of the basic SQP algorithm applied here is described in [25].

2.4. Genetic Algorithms

GA is a metaheuristic technique inspired by genetics and Darwin’s theory of evolution [16].
GA operates on a set of individuals (chromosomes), P o p ( t ) = x 1 t ,   ,   x n t , for a certain generation t. Each chromosome symbolises a potential problem–solution. Each solution is encoded as a finite vector of variables of a certain alphabet, usually a binary one. The individuals of the offspring are modified by the genetic operators’ crossover and mutation [26].
The workflow of a common GA can be roughly presented as follows. The initial set of solutions can be randomly generated. To evaluate the capability of a chromosome compared to the whole population, the fitness function is used. Chromosomes with high values of the fitness function are more likely to be chosen for reproduction. The operators (crossover and mutation) are used in order to reproduce a new offspring, which becomes generation t + 1. The offspring replaces chromosomes existing in the population, and the generation is reiterated until the stopping conditions are fulfilled.

3. Results and Discussion

3.1. Experimental Studies

Experimental studies of batch anaerobic digestion of wheat straw at thermophilic conditions are performed. The mechanically chopped wheat straw with the corresponding quantity (30 mL/dm3 or 35 mL/dm3) is suspended in distilled water, and inoculum from a working biomethane generating bioreactor 30% (vol) is also introduced in the batch process, which continues for 15 days. The inoculum quantity is previously determined as it ensures good buffer capacity at the beginning of the process. Stirring is maintained at 100 min−1. The measurements of cellulose, glucose, acetate, and methane content in the biogas in the bioreactor are available for 30 and 35 mL/dm3 initial substrate organic load. The results are presented in Table 1 for 30 mL/dm3 and in Table 2 for 35 mL/dm3. All concentrations given are mass concentrations, not molar concentrations.

3.2. Mathematical Modelling

3.2.1. Structural Identification

A deterministic nonlinear dynamic system of differential equations was used to describe the thermophilic anaerobic digestion process for methane production of wheat straw as follows:
d S 0 d t = β X 1 S 0
d X 1 d t = μ 1 X 1
d S 1 d t = β X 1 S 0 1 Y 1 μ 1 X 1
d X 2 d t = μ 2 X 2
d S 2 d t = Y b μ 1 X 1 1 Y 2 μ 2 X 2
Q C H 4 = Y C H 4 μ 2 X 2
μ 1 = μ 1 max S 1 k 1 s + S 1 b 1
μ 2 = μ 2 max S 2 k 2 s + S 2 b 2
Equations (3) and (5) present the balance of the effluent substrate (S0)—cellulose and its transformation to glucose (S1) after hydrolysis. Equations (4) and (6) describe the dynamics of the biomass concentrations (X1 and X2). Equation (7) presents the dynamics of the intermediate product—acetate (Ac1), and the algebraic Equation (8)—the flow rate of the methane in the gas phase in the bioreactor. For the specific growth rates of the biomasses, a Monod type function was adopted—Equations (9) and (10).

3.2.2. Parameter Identification

The identification of model parameters values is based on dynamical experiments of thermophilic anaerobic digestion of wheat straw carried out in a 2 dm3 bioreactor with 30 mL/dm3 organic load. It is fulfilled using two different techniques: deterministic SQP algorithm and metaheuristics—GA.
The numerical computations are carried out on a PC/Intel Core i5-2320 [email protected] GHz, 4 GB Memory (RAM), Windows 7 (64 bit) operating system. Modelling and numerical experiments are performed using Matlab R2013a.
To locate solution point algorithms such as SQP, usually involves a few iterations. Such local search methods may fall into local minima. It is known that SQP estimates are very sensitive to the initial solutions. At the same time, GA does not have such a dependence, as the initial solutions (initial population) are generated randomly within the predetermined lower and upper limits of each model parameter.
According to the problems considered here, the GA parameters such as population number, crossover and mutation probability, etc., are tuned previously. A series of tests are conducted in order to reach the best solution quality. The optimal GA parameters are presented in Table 3.
When comparing stochastic algorithms, such as GA, the general approach is to operate with the obtained average results. The GA discussed here has been executed 30 times under the same conditions, e.g., algorithms’ parameters and functions, presented in Table 3. The obtained results are summarised in Table 4. To better distinguish the results from the algorithms applied for the considered problem, in addition to the obtained values of the objective function, the values of Mean Absolute Error (MAE) [27] are calculated.
Primarily, adjustment of the coefficients was carried out with the SQL algorithm, setting wider search limits. The so-found model does not describe experimental data well (SQP1 in Table 4). The observed objective value of the model during the identification is J = 13.91 and MAE = 2.55. For this reason, initially, GA is run for a few iterations, e.g., 100, to give some closure to the optima model parameters estimates. The new lower and upper limits for SQP are defined using these results. Based on the new observed results, a mathematical model with a much higher accuracy is obtained (SQP2 in Table 4)—J = 10.22 and MAE = 0.64.
In the considered identification problem, the SQP and GA optimisation algorithms find two different solutions for the model parameters values. The estimates for μ 1 max , β , μ 2 max , and Y 2 are in similar ranges; the model parameters Y b , Y C H 4 , and b 2 are identified with a small difference; while the estimates of k 1 s , b 1 , Y 1 , and k 2 s are very different. As no information has been found in the literature on the values of the respective coefficients in such processes, it is not possible to conclude which coefficients should be preferred.
The presented comparison (Figure 1 and Figure 2) of the experimental data and model predictions for cellulose, glucose, acetate, and methane obtained by both algorithms demonstrates the better performance of the proposed GA compared to the SQP algorithm.

3.2.3. Model Validation

Verification of the developed mathematical models is carried out based on the different data set of the process, i.e., 35 mL/dm3 initial substrate organic load. The proposed models based on GA and SQP2 algorithms, with the parameters listed in Table 4, are simulated with the new initial substrate organic load (35 mL/dm3 instead of 30 mL/dm3), as well as with the new initial values of cellulose, glucose, acetate, and methane in the new experimental data set. The obtained results (from the comparison between the model and the new data set) are shown in Figure 3 (cellulose and glucose dynamics) and Figure 4 (acetate and methane dynamics). The observed objective value of the model during the verification is JSQP = 18.22 and JGA = 17.91.
As seen in Figure 3 and Figure 4, both GA and SQP models predict well the new set of the real data for variables cellulose, glucose, acetate, and methane process. However, GA visually fits the experimental data with a higher degree of accuracy, as established by the numerical results for the objective function value J.
Simulation models development for the anaerobic digestion process of municipal solid wastes [28], or animal manure and sewage sludge [29,30], had been carried out with the aim of coping with the harmful effects of such disposals on the health of humans, animals, and the environment. In this study, the effect of elevated temperature was revealed as suitable to be explored with the feasibility of improving the methane yield of wheat straw as agricultural waste in anaerobic digestion, as it increased the rate of the process and made unnecessary harsh pretreatment techniques. The created model presented could serve for monitoring and control of such processes.

4. Conclusions

In this paper, a new model of the anaerobic digestion of wheat straw is presented. The nonlinear model structure depicted is quite simple but with many coefficients and proven for the anaerobic digestion processes. Parameter identification is made using two different methods—deterministic SQP algorithm and metaheuristic algorithm—GA.
SQP separately has not given satisfactory results with the originally defined initial solutions and parameters limits. GA has been applied here as a competitive metaheuristic technique, proven as an effective optimisation method for solving complex, real problems, such as parameters estimation of yeast and bacterial fermentation models. The best results of SQP are obtained using a technique which combines GA and SQP algorithm. The obtained numerical results regarding the computational efficiency of the algorithms unambiguously show that GA have an advantage over SQP.
Results illustrating the model validity are also presented using a different dataset, which showed good performance with the error even smaller compared with the error from the identification. The results show good model validity.
One of the most difficult questions to answer concerns the interpretation of the coefficients of the model. At the moment in the literature, there is no information about the values of the considered coefficients of the processes models proposed in this paper. This is the reason why it is not possible to draw conclusions about which coefficients should be preferred.
In a future plan, the hybrid SQP-GA technic, for model identification could be developed, aiming to take advantage of different optimisation strategies while avoiding the algorithms’ weaknesses. Designing an appropriate hybridisation of metaheuristics, an algorithm with better performance can be offered.

Author Contributions

Conceptualization, E.C., L.K., I.S. and O.R.; methodology, E.C., I.S. and V.H.; software, E.C. and O.R.; validation, E.C. and O.R.; formal analysis, E.C., L.K. and V.H.; investigation, E.C., L.K., V.H., I.S. and O.R.; data curation, E.C., L.K. and V.H.; writing—original draft preparation, E.C., L.K., V.H. and O.R.; writing— E.C., L.K., V.H. and O.R.; visualization, E.C. and O.R.; supervision, E.C. and O.R.; project administration, E.C.; funding acquisition, E.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research has been supported by the Bulgarian National Science Fund, Grant No KP-06-N-26/5.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

S 0 Cellulose concentration[mL/dm3]
X 1 Acidogenic bacteria concentration[mL/dm3]
S 1 Glucose concentration[mL/dm3]
X 2 Methanogenic bacteria concentration[mL/dm3]
S 2 Acetate concentration[mL/dm3]
Q C H 4 Methane flow rate[mL/dm3]
μ 1 Specific growth rate of acidogenic bacteria[h−1]
μ 2 Specific growth rate of methanogenic bacteria[h−1]
μ1maxMaximum specific growth rate of acidogenic bacteria[h−1]
k 1 s Saturation coefficient for acidogenic bacteria[mL/dm3]
b 1 Maintenance rate of acidogenic bacteria[h−1]
β Coefficient of biodegradability[L/(g·h)]
Y 1 Yield coefficient for acidogenic bacteria[-]
Y b Coefficient[-]
μ2maxMaximum specific growth rate of methanogenic bacteria[h−1]
k 2 s Saturation coefficient for methanogenic bacteria[mL/dm3]
b 2 Maintenance rate of methanogenic bacteria[h−1]
Y 2 Yield coefficient for methanogenic bacteria[-]
Y C H 4 Yield coefficient for methane [mL/dm3]

References

  1. Schneider, S.H.; Easterling, W.E.; Mearms, L.O. Adaptation: Sensitivity to natural variability, agent assumptions, and dynamic climatic changes. Clim. Change 2000, 45, 203–221. [Google Scholar] [CrossRef]
  2. Haryana, A. Biomass Utilization as Renewable Energy for Optimization of National Energy Mix. Bappenas Work. Pap. 2018, 1, 55–65. [Google Scholar] [CrossRef]
  3. Guragain, Y.N.; Vadlani, P.V. Renewable Biomass Utilization: A Way Forward to Establish Sustainable Chemical and Processing Industries. Clean Technol. 2021, 3, 243–259. [Google Scholar] [CrossRef]
  4. Dahlgren, S. Biogas-based fuels as renewable energy in the transport sector: An overview of the potential of using CBG, LBG and other vehicle fuels produced from biogas. Biofuels 2020, 1–13. [Google Scholar] [CrossRef]
  5. Amon, T.; Amon, B.; Kryvoruchko, V.; Machmüller, A.; Hopfner-Sixt, K.; Bodiroza, V.; Hrbek, R.; Friedel, J.; Pötsch, E.; Wagentristl, H.; et al. Methane production through anaerobic digestion of various energy crops grown in sustainable crop rotations. Bioresour. Technol. 2007, 98, 3204–3212. [Google Scholar] [CrossRef]
  6. Vadlani, P.V. Financing strategies for sustainable bioenergy and the commodity chemicals industry. In Green Energy to Sustainability: Strategies for Global Industries, 1st ed.; Vertès, A.A., Qureshi, N., Blaschek, H.P., Yukawa, H., Eds.; John Wiley and Sons Ltd.: Hoboken, NJ, USA, 2020; pp. 569–586. [Google Scholar]
  7. Bakker, R.R.C.; Elbersen, H.W.; Poppens, R.P.; Lesschen, J.P. Rice Straw and Wheat Straw-Potential Feedstocks for the Biobased Economy; NL Agency: Utrecht, The Netherlands, 2013. [Google Scholar]
  8. Zheng, Q.; Zhou, T.; Wang, Y.; Cao, X.; Wu, S.; Zhao, M.; Wang, H.; Xu, M.; Zheng, B.; Zheng, J.; et al. Pretreatment of wheat straw leads to structural changes and improved enzymatic hydrolysis. Sci. Rep. 2018, 8, 321. [Google Scholar] [CrossRef] [Green Version]
  9. Shang, G.; Zhang, C.; Wang, F.; Qiu, L.; Guo, X.; Xu, F. Liquid hot water pretreatment to enhance the anaerobic digestion of wheat straw—Effects of temperature and retention time. Environ. Sci. Pollut. Res. 2019, 26, 29424–29434. [Google Scholar] [CrossRef]
  10. Bolado-Rodriguez, S.; Toquero, C.; Martin-Juarez, J.; Travaini, R.; Garcia-Encina, P.A. Effect of thermal, acid, alkaline and alkaline-peroxide pretreatments on the biochemical methane potential and kinetics of the anaerobic digestion of wheat straw and sugarcane bagasse. Bioresour. Technol. 2016, 201, 182–190. [Google Scholar] [CrossRef]
  11. Moset, V.; Poulsen, M.; Wahid, R.; Hojberg, O.; Moller, H.B. Mesophilic versus thermophilic anaerobic digestion of cattle manure: Methane productivity and microbial ecology. Microb. Biotechnol. 2015, 8, 787–800. [Google Scholar] [CrossRef]
  12. Gómez, X.; Blanco, D.; Lobato, A.; Calleja, A.; Martínez-Núñez, F.; Martin-Villacorta, J. Digestion of cattle manure under mesophilic and thermophilic conditions: Characterization of organic matter applying thermal analysis and 1H NMR. Biodegradation 2011, 22, 623–635. [Google Scholar] [CrossRef]
  13. Byrd, R.H.; Curtis, F.E.; Nocedal, J. Infeasibility Detection and SQP Methods for Nonlinear Optimization. SIAM J. Optim. 2010, 20, 2281–2299. [Google Scholar] [CrossRef] [Green Version]
  14. Gill, P.E.; Wong, E. Sequential Quadratic Programming Methods, UCSD Department of Mathematics; Technical Report NA-10-03; University of California: San Diego, CA, USA, 2010. [Google Scholar]
  15. Holland, J.H. Adaptation in Natural and Artificial Systems, 2nd ed.; MIT Press: Cambridge, UK, 1992. [Google Scholar]
  16. Goldberg, D. Genetic Algorithms in Search, Optimization and Machine Learning; Addison-Wesley Publishing Company: Boston, MA, USA, 1989. [Google Scholar]
  17. Chopard, B.; Tomassini, M. Performance and limitations of metaheuristics. In An Introduction to Metaheuristics for Optimization; Springer International Publishing: Berlin/Heidelberg, Germany, 2018; pp. 191–203. [Google Scholar]
  18. Dillen, W.; Lombaert, G.; Schevenels, M. Performance Assessment of Metaheuristic Algorithms for Structural Optimization Taking Into Account the Influence of Algorithmic Control Parameters. Front. Built Environ. 2021, 7, 618851. [Google Scholar] [CrossRef]
  19. Djebedjian, B.; Abdel-Gawad, H.A.A.; Ezzeldin, R.M. Global performance of metaheuristic optimization tools for water distribution networks. Ain Shams Eng. J. 2021, 12, 223–239. [Google Scholar] [CrossRef]
  20. Roeva, O.; Zoteva, D.; Lyubenova, V. Escherichia coli Cultivation Process Modelling Using ABC-GA Hybrid Algorithm. Processes 2021, 9, 1418. [Google Scholar] [CrossRef]
  21. Roeva, O.; Angelova, M.; Zoteva, D.; Pencheva, T. Water Cycle Algorithm for Modelling of Fermentation Processes. Processes 2020, 8, 920. [Google Scholar] [CrossRef]
  22. Roeva, O.; Atanassova, V. Cuckoo search algorithm for model parameter identification. Int. J. Bioautom. 2016, 20, 483–492. [Google Scholar]
  23. Updegraff, D.M. Semimicrodetermination of cellulose inbiological materials. Anal. Biochem. 1969, 32, 420–424. [Google Scholar] [CrossRef]
  24. Miller, G.L. Use of Dinitrosalicylic Acid Reagent for Determination of Reducing Sugar. Anal. Chem. 1959, 31, 426–428. [Google Scholar] [CrossRef]
  25. Nocedal, J.; Wright, S.J. Numerical Optimization, 2nd ed.; Springer Series in Operations Research; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  26. Roeva, O.; Zoteva, D.; Castillo, O. Joint Set-up of Parameters in Genetic Algorithms and the Artificial Bee Colony Algorithm: An Approach for Cultivation Process Modelling. Soft Comput. 2020, 25, 2015–2038. [Google Scholar] [CrossRef]
  27. Everitt, B.S.; Skrondal, A. The Cambridge Dictionary of Statistics; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  28. Calise, F.; Cappiello, F.L.; Dentice d’Accadia, M.; Infante, A.; Vicidomini, M. Modeling of the Anaerobic Digestion of Organic Wastes: Integration of Heat Transfer and Biochemical Aspects. Energies 2020, 13, 2702. [Google Scholar] [CrossRef]
  29. Xu, F.; Li, Y.; Ge, X.; Yang, L.; Li, Y. Anaerobic digestion of food waste–Challenges and opportunities. Bioresour. Technol. 2018, 247, 1047–1058. [Google Scholar] [CrossRef] [PubMed]
  30. Liu, H.; Wang, X.; Fang, Y.; Lai, W.; Xu, S.; Lichtfouse, E. Enhancing thermophilic anaerobic co-digestion of sewage sludge and food waste with biogas residue biochar. Renew. Energy 2022, 188, 465–475. [Google Scholar] [CrossRef]
Figure 1. Comparison of the experimental data and model predictions for substrates cellulose and glucose obtained by GA and SQP models—identification results.
Figure 1. Comparison of the experimental data and model predictions for substrates cellulose and glucose obtained by GA and SQP models—identification results.
Processes 10 00742 g001
Figure 2. Comparison of the experimental data and model predictions for acetate and methane obtained by GA and SQP models—identification results.
Figure 2. Comparison of the experimental data and model predictions for acetate and methane obtained by GA and SQP models—identification results.
Processes 10 00742 g002
Figure 3. Comparison of the experimental data and model predictions for substrates cellulose and glucose obtained by GA and SQP models—validation results.
Figure 3. Comparison of the experimental data and model predictions for substrates cellulose and glucose obtained by GA and SQP models—validation results.
Processes 10 00742 g003
Figure 4. Comparison of the experimental data and model predictions for acetate and methane obtained by GA and SQP models—validation results.
Figure 4. Comparison of the experimental data and model predictions for acetate and methane obtained by GA and SQP models—validation results.
Processes 10 00742 g004
Table 1. Batch experiments with 30 mL/dm3 organic load.
Table 1. Batch experiments with 30 mL/dm3 organic load.
Duration,
day
Cellulose,
mL/dm3
Glucose,
mL/dm3
Acetate,
mL/dm3
Methane,
mL/dm3
010.730.0130.130.0000
19.690.0531.360.0099
28.840.0821.690.0771
38.160.1011.000.1224
47.620.1130.110.1313
57.200.1190.140.0301
66.900.1220.130.0122
76.680.1220.110.0136
86.540.1220.090.0151
96.450.1220.090.0130
106.390.1240.090.0120
116.350.1260.090.0110
126.300.1260.080.0100
136.230.1240.080.0090
146.130.1220.080.0090
155.960.1200.080.0090
Table 2. Batch experiments with 35 mL/dm3 organic load.
Table 2. Batch experiments with 35 mL/dm3 organic load.
Duration,
day
Cellulose,
mL/dm3
Glucose,
mL/dm3
Acetate,
mL/dm3
Methane,
mL/dm3
011.040.0190.150.000
19.830.0601.400.010
28.850.0921.700.012
38.070.1161.270.140
47.480.1320.960.150
57.040.1430.430.100
66.730.1490.280.047
76.540.1520.120.030
86.440.1530.080.030
96.410.1520.080.030
106.420.1520.080.030
116.450.1530.080.030
126.480.1520.080.030
136.490.1500.080.020
146.450.1500.080.020
156.340.1500.080.020
Table 3. GA parameters.
Table 3. GA parameters.
ParameterGA
Population number (Npop)200
Generation gap (ggap)0.97
Maximal number of iterations (iter_max)100
Crossover probability (pc)0.75
Mutation probability (pm)0.01
Table 4. Results from model parameter identification based on SQP and GA.
Table 4. Results from model parameter identification based on SQP and GA.
ParameterDeterministic
Algorithm
Metaheuristic
Algorithm
Deterministic
Algorithm
SQP1GASQP2
ValueBestValue
J13.919.8710.22
MAE2.550.620.64
μ 1 max 0.290.740.98
k 1 s 1.6719.1610.96
b 1 0.100.00030.004
β 0.340.510.52
Y 1 0.040.0010.001
Y b 2.0319982476
μ 2 max 0.53.694.85
k 2 s 0.311.0051.83
Y 2 0.187.225.74
Y C H 4 3.650.0050.01
b 2 0.10.0460.029
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chorukova, E.; Kabaivanova, L.; Hubenov, V.; Simeonov, I.; Roeva, O. Mathematical Model of a Thermophilic Anaerobic Digestion for Methane Production of Wheat Straw. Processes 2022, 10, 742. https://doi.org/10.3390/pr10040742

AMA Style

Chorukova E, Kabaivanova L, Hubenov V, Simeonov I, Roeva O. Mathematical Model of a Thermophilic Anaerobic Digestion for Methane Production of Wheat Straw. Processes. 2022; 10(4):742. https://doi.org/10.3390/pr10040742

Chicago/Turabian Style

Chorukova, Elena, Lyudmila Kabaivanova, Venelin Hubenov, Ivan Simeonov, and Olympia Roeva. 2022. "Mathematical Model of a Thermophilic Anaerobic Digestion for Methane Production of Wheat Straw" Processes 10, no. 4: 742. https://doi.org/10.3390/pr10040742

APA Style

Chorukova, E., Kabaivanova, L., Hubenov, V., Simeonov, I., & Roeva, O. (2022). Mathematical Model of a Thermophilic Anaerobic Digestion for Methane Production of Wheat Straw. Processes, 10(4), 742. https://doi.org/10.3390/pr10040742

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop