Maximum production point tracking method for a solar-boosted biogas energy generation system

Low biogas yield in cold climates has brought great challenges in terms of the flexibility and resilience of biogas energy systems. This paper proposes a maximum production point tracking method for a solar-boosted biogas generation system to enhance the biogas production rate in extreme climates. In the proposed method, a multi-dimensional R–C thermal circuit model is formulated to analyze the digesting thermodynamic effect of the anaerobic digester with solar energy injection, while a hydrodynamic model is formulated to express the fluid dynamic interaction between the hot-water circulation flow and solar energy injection. This comprehensive dynamic model can provide an essential basis for controlling the solar energy for digester heating to optimize anaerobic fermentation and biogas production efficiency in extreme climates. A model predictive control method is developed to accurately track the maximum biogas production rate in varying ambient climate conditions. Comparative results demonstrate that the proposed methodology can effectively control the fermentation temperature and biogas yield in extreme climates, and confirm its capability to enhance the flexibility and resilience of the solar-boosted biogas generation system.


Motivation
Environmental pollution and global warming have prompted the research community to search for renewable energy sources (RESs) as valuable solutions [1]. Large-scale utilization of RESs to replace fossil energy for power generation has become the common goal of countries around the world. Biomass can play a vital role in promoting the use of RESs and is considered as the fourth largest energy resource in the world. Biomass can be converted to biogas through anaerobic fermentation, which can be used to supply local multi-energy demands using a combined heating and power (CHP) unit [2].
For mesophilic digestion, fermentation temperature is the main factor that affects the anaerobic fermentation to promote the degradation and metabolism of anaerobic bacteria [3]. Low ambient temperature can result in extremely low digestion efficiency, and the anaerobic fermentation process is suspended when the temperature falls below around zero [4]. This low-efficiency biogas production is the main obstacle to the promotion of biogas application and development in cold areas [5].
Recent studies have looked at the biogas yield problem in extreme climates. Many physical and chemical approaches have been adopted [6][7][8][9] to improve the fermentation temperature, such as straw and plastic film laying on the periphery of biogas digester, electrothermal film heating, boiler heating, waste heat heating from CHP, etc. The traditional heating approaches rely on the non-renewable fossil fuels, which release greenhouse gas during the utilization processes and have negative impact

Open Access
Protection and Control of Modern Power Systems on the social-economic environment [10]. With the soaring price of fossil fuel and increasing environmental concerns, solar energy technologies are widely harnessed to improve anaerobic fermentation efficiency, including solar collectors and solar greenhouses [11][12][13][14][15][16][17][18]. In addition to electrical generation, CHP waste heat can be used to warm up the biogas digester, which can provide a low-cost and efficient solution for biogas production in cold areas. The inherent intermittency and volatility of solar radiation can lead to uncertainty of the external heating energy injection, thereby accelerating fermentation temperature fluctuation. Solar energy heating can be controlled by the hot-water flow to maintain the optimal temperature to ensure maximum biogas production. Therefore, this paper is devoted to establishing a solarboosted biogas generation system to manage the fermentation temperature in variable climates and rapidly track the maximum biogas production point for flexibility and resilience enhancement. Maximum production point tracking (MPPT) involves not only the thermodynamic effects of solar energy injection but also the hydrodynamic interaction between the water flow and pipe pressure. In such cases, the fermentation temperature requires comprehensive dynamic modelling of water flow, pressure fluctuation, convective heat transfer, and conductive heat transfer in the digester climate. This could lead to high nonlinearity and increase the computational burden throughout the MPPT optimization scheme. Therefore, MPPT is a challenging problem that is difficult to solve with conventional methods because of the uncertainty of solar energy, multi-energy coupling, high-dimensional variables, nonlinear differential thermodynamics and hydrodynamics.

Literature review
Traditional physical insulation measures, such as covering the ground around the digester with charcoal [9], laying straw and plastic film on the periphery of the biogas digester, are commonly used to prevent heat loss and maintain fermentation temperature. However, by adding insulation materials, the fermentation temperature can only be passively sustained rather than proactively adjusted. Thus, several studies have looked at solving the low biogas yield problem with external energy injection. In [4], engine waste heat is employed to provide the steam required for thermal hydrolysis pre-treatment and the heat for anaerobic digestion. Electricity from natural gas and fossil fuel to heat the anaerobic digestion is also effective in improving the biogas yield [7], though waste of primary energy and negative impact on the environment cannot be avoided.
The focus of recent research is to use RESs to facilitate relatively high fermentation temperature, including the ground source heat pump [19], biomass boiler [20], groundwater source heat pump [6], methane liquid heat recovery pump etc. However, considering availability and economy, these methods are not suitable for biogas production in rural areas. Compared with other renewable technologies, solar energy is independent of the influence of the geographic environment and offers eco-friendly and least expensive heating sources [4,21]. It has been proved that solar heating is an efficient way to facilitate biogas production. In [15], a solar heating biogas fermentation system is tested for five months and the heat collected by the 2m 2 solar collector after eight hours in October can meet the heat demands of 6m 3 biogas digester's fermentation at Hohhot. A solar heated biogas system integrated with a phase change thermal storage device is introduced in [16] and the preliminary results show that the system could be a promising solution to promote the application of biogas technology in cold areas. The solar-heated biogas digester is modeled in [17] to simulate and predict its thermal behavior. Two solar greenhouses are used to heat the digester in [18] and the simulation indicates that the proposed system can increase biogas production and revenue compared to the system with no heating resource. Nevertheless, existing work has mainly focused on the enhancement of steadystate production performance while trying to verify the economical and technical feasibilities of the solar-heated biogas system. The system thermodynamic and hydrodynamic behaviors, which are essential for the development of advanced control strategies to track the maximum biogas production, are ignored.

Contribution and organization
In this paper, a model predictive control (MPC)-based biogas MPPT method is proposed to enhance the biogas production rate of the solar-boosted biogas generation system in extreme climates. The contributions of the paper are: 1. A dynamic model of a solar-boosted biogas energy generation system is formulated to describe the thermodynamic and hydrodynamic interactions during the biogas production process with solar energy injection. An R-C thermal network is modeled to analyze the digesting thermodynamic effect with solar energy injection, while a hydrodynamic model is formulated to control the solar hot-water circulation flow. Compared with the previous solar-heated biogas system with uncontrollable dynamic modelling, this model allows accurate control of thermodynamic and hydrodynamic behaviors, which can provide an essential basis for controlling the solar energy for digester heating.

A MPC-based biogas MPPT method is proposed
to maximize the biogas production rate in variable climates. The strategy systematically integrates the rolling optimization and feedback mechanism based on short-term prediction of the solar radiation and ambient temperature, to optimize the digesting fermentation temperature, so as to track the maximum biogas production rate and enhance the resilience of the solar-boosted biogas generation system to varying ambient climate conditions.
The rest of the paper is organized as follows. Section 2 illustrates the low biogas yield problem in extreme climates and shows the proposed biogas generation system. Section 3 describes the proposed MPC-based MPPT method, while Sect. 4 analyzes the performance of the proposed method through simulation. Section 5 presents the conclusions.

Problem description
Among all the factors that affect the anaerobic fermentation, temperature is of greatest importance. For mesophilic anaerobic fermentation, the increase in temperature of the anaerobic digester will promote the degradation and metabolism of anaerobic bacteria, thus improving the biogas production [22]. Therefore, solar energy is utilized for biogas digester heating to enhance the fermentation temperature in this paper. The hotwater circulation for biogas digester heating is illustrated in Fig. 1. As can be seen, the cold water in the tank is pumped to the solar thermal collector. Then solar radiation is absorbed and converted to thermal energy for water heating in the pipe, which is transmitted to the digester tank to heat the biogas fermentation liquid. However, with further increase of fermentation temperature, such overtemperature can reduce the activity of microorganisms in the fermentation process, resulting in inefficient biogas production. Therefore, it is important to develop control strategies that can track the maximum biogas production by optimizing the fermentation temperature. This paper aims to maintain the fermentation temperature at the optimal level so that the system can output maximum production in real time, i.e., to achieve MPPT. With the digesting thermodynamic effects on biogas production, the fermentation temperature can be influenced by external solar energy injection. The heat transmitted to the digester can be controlled by the hot-water flow, while water flow is highly nonlinear to the pressure drop along pipelines based on hydrodynamic analysis. In addition, the inherent intermittency and volatility of solar radiation make the analysis of the associated optimal control problem challenging. In order to maintain optimal biogas production, the primary objectives of this paper are to propose: a biogas generation system based on thermodynamic and hydrodynamic analysis, and a MPC-based MPPT to enable the biogas generation system to rapidly track the desired maximum point during dynamic changes of ambient temperature.

Digesting thermodynamic analysis
For biogas digestion, temperature is the primary environmental factor that affects the digestion rate and process. In general, the increase of fermentation temperature will lead to the increase of the substrate use and microbial metabolic rate, thereby improving the biogas production. However, if the fermentation temperature exceeds the optimum temperature, the activity of anaerobic bacteria and the efficiency of the digestion process will be reduced. According to the field measured data in [23], biogas production rate f bio is a function of fermentation temperature and can be modeled as: where B o is the biochemical methane potential; S o is the influent volatile solids concentration; HRT is the hydraulic retention time; µ m is the maximum specific growth rate of microorganisms; T d is the fermentation temperature [24]. K is a dimensionless kinetic parameter related to the concentration of influent volatile solids [25], and is given as: The relationship between the maximum specific growth rate of microorganisms and the fermentation temperature is given as [25]: The available solar energy can be used for digester heating, thereby promoting the anaerobic fermentation process for biogas production enhancement [26]. A multi-dimensional R-C thermal circuit model of the underground biogas digester is established to analyze the digesting thermodynamic effects under external energy injection. Here, the heat transfer and storage are represented by thermal resistances and capacitances. The inside, internal wall, peripheral wall and ambient temperatures of the digestion reactor can be characterized by nodes with thermal potentials. Using the lumped model analysis of the structure and environmental conditions of the cylindrical digester in [24], three thermal circuital branches are established to analyze the nodal digesting thermodynamic effects on different surfaces. For the i th branch, the internal wall is the insulation layer and can be denoted as a node connected to the node inside the digester by two resistances Z ini /2 and Z di , and the peripheral wall node is connected to the external node by two series resistances Z wi /2 and Z outi . Two resistances Z ini /2 and Z wi /2 are used to represent the thermal interaction and transfer of the two-layer walls. Thus, the digesting thermodynamic effects can be modeled based on the nodal analysis of thermal circuits, as: where C d , C ini , and C wi represent the thermal capacitances of the digester inside, the i th surface of the insulation layer and the peripheral wall, respectively; Z di , Z outi , Z ini , Z wi and Z s represent the thermal resistances that are used for the convective heat transfer inside the digester and through the digester cover, conductive heat transfer inside, peripheral walls and soil, respectively; T d , T out , T s , T ini , and T wi are the temperatures inside and outside the digester, the surrounding soil, the first and second layer walls, respectively; k represents the operating time; ρ water is the water density; E sf represents the controllable amount of thermal energy injection from renewable harvesting; T hw,k is the temperature of hot-water flow injected into the digester at time k.

Hydrodynamic modelling of solar water circulation
The heat transmitted to the digester can be controlled by the hot-water flow. A hydrodynamic networking model is established to describe the interaction between the water flow and the pressure drop along the pipelines. Solar collector, heat storage tank, throttle valve and biogas fermentation tank can be represented by nodes with specific pressure as shown in Fig. 2. Between adjacent nodes, the flow coefficient of the pipeline is used to describe the flow resistance, and there is a nonlinear relationship between water flow and pressure drop along the transmission pipeline [27]. According to the node analysis of the hydrodynamic network loop, the water flow Q d,k from the throttle valve to the digester will change with the pressure at the throttle valve outlet, which can be described as: where α bp is the flow coefficient of the pipeline; U d,k is the pressure of the biogas digester at time k; U val,out,k is the flow pressure at the outlet.
Other water flows as well as their temperatures in the model are described as: where Q res,out,k is the water flow from the heat storage tank to the throttle valve; Q PTv,k is the total water flow from the collector; Q PTv,k , Q PTr,k are the water flows from the collector to the throttle valve and the tank, respectively; U res,in,k and U res,out,k are the respective inlet and outlet pressures of the tank; U PT,k and U val,in,k are the respective pressures of the solar collector and inlet of the throttle valve; T PT,k and T res,k are the temperatures of the collector and tank at time k, respectively. The throttle valve controls the water flow by changing the throttling area. Assuming that the throttling process is an isenthalpic process, the water flow through the throttle valve at time k Q val,k has the following relationship with the inlet and outlet pressures: where β val is the flow coefficient; A val,k is the throttling area at time k; U val,in,k is the flow pressure at the inlet; A val,min and A val,max are the respective lower and upper limits of the allowable throttle area; U val,in,min , U val,in,max , U val,out,min and U val,out,max are the allowable pressure limits at the inlet and outlet of the throttle valve. The throttling area of the valve can be changed by adjusting the opening of the valve, so as to control the water flow. It is noted from (17) that the change of throttling area will lead to a pressure change at the outlet of the throttle valve, as: When the light and heat are insufficient, the heat storage tank provides heat so the water flow of the throttle valve is equal to the outlet water flow of the storage tank and the collector. In contrast, when the light and heat are sufficient, the excess thermal energy is stored in the tank so the water flow of the throttle valve is equal to the water flow of the collector.

Solar-biogas generation system for multi-energy supplies
The use of solar energy to heat the biogas digester is efficient in overcoming the obstacle of low biogas yield in a cold climate and improving solar energy penetration. Figure 3 illustrates a solar-biogas generation system with multi-energy carriers, where solar and biomass energy are converted into electricity and heat through the solar thermal collector and biogas digester. Given the thermodynamic effects of fermentation temperature on biogas production, the thermal energy from the solar collector is used to heat the biogas digester, thereby promoting anaerobic fermentation to enhance the biogas production rate. The biogas can be stored in the biogas storage tank or converted into electricity and heat via a CHP unit to supply local multi-energy demands [28]. The hot-water storage tank and biogas storage tank can provide sufficient storage capacity for the available heat and biogas. Therefore, a variety of energy conversion and storage devices are integrated to combine and couple these carriers to form a centerinternal energy redundant connection, and can provide flexibility and synergy for multi-energy supplies.  Figure 4 shows the MPC-based MPPT framework. The predicted solar radiation and ambient temperature can be obtained from the historical observation data using a nonlinear regression method [3]. Monte Carlo simulations are implemented to generate a series of scenarios to capture the prediction uncertainties of solar radiation and ambient temperature [5]. Thus the fermentation temperature and biogas production rate in each time slot can be obtained from the multi-dimensional R-C thermodynamic model of the biogas digester and the hydrodynamic model of the hot-water circulation flow based on the prediction scenarios and input parameters. During the optimization process, the control decision is obtained by solving the optimization problem at the current time while taking into account the uncertainties of solar radiation and ambient temperature of future horizons [29], while only the optimal control decision for the current time k is implemented as the actual control input [30].
Considering the inaccuracy of the predicted outputs obtained from the proposed model, a feedback mechanism is adopted based on the actual outputs of the solarboosted biogas generation system at the current time k to tune the multi-dimensional R-C thermodynamic model of the biogas digester and the hydrodynamic model of the hot-water circulation flow, as well as the prediction of solar radiation and ambient temperature.

Objective function
The optimization objective aims to minimize the difference between the maximum biogas production and the biogas production at the current time k, as well as the expected biogas production of all remaining future scenarios, which can be formulated as: where K end represents the last time slot of the rolling horizon; N s represents the total number of scenarios; π s is the probability of scenario s; f bio,k,s is the calculated biogas production at time k; and f bio,max is the maximum biogas production of the mesophilic fermentation digestion.

Solar thermal collector
The flat-plate solar collectors have simple mechanical structure and stable performance, and are widely used in household and industrial applications. They usually include a black absorber with a glass cover and insulation. While the glass cover with low iron content can increase light transmittance and reduce radiation loss, the insulation on the back and sides can reduce conduction losses. The plate can absorb and convert solar radiation into useable heat to heat the water up to 90 °C while the thermal efficiency can reach 70%. The useful energy E PT,k is calculated based on the energy balance as [31]: where G PT,k is the solar radiation at time k; A s is the lighting area of the collector; η s is the efficiency of the solar collector; C water is the specific heat capacity of the hotwater in the solar thermal collector; V PT,max is the rated volume of the collector; U PT,min and U PT,max are the allowable pressure limits of the collector.

Stratified heat storage tank
Solar energy has significant instability and discontinuity due to the changes of weather and other factors. Storage of solar energy is necessary when the fermentation temperature exceeds the optimal temperature, while the hot-water stored in the tank when the solar radiation intensity is sufficient can provide thermal feedback for the biogas digester in an underlit scene, such as a cloudy day or night.
A natural layered heat storage tank is adopted, and its specific structure is shown in Fig. 1. The tank is directly connected to the solar thermal collector where the outlet and return pipes are located above and below the tank, respectively. The inlet ports of the solar collector and the outlet ports of the biogas tank are located at the high temperature top of the tank whereas the outlet ports of (25) the solar collector and the inlet ports of the biogas tank are located at the low temperature bottom of the tank. The density of water varies with temperature, so the returning water tends to flow up or down to the storage tank because its temperature is significantly different from the temperature at the top and bottom of the tank [32].
Assuming that the heat losses in and out of the water pipe are negligible, the energy balance of the heat storage tank is given as: The operating constraints of the tank are: where E res,total,k and E res,loss,k are the total heat and heat loss in the tank at time k, respectively; m res,k is the mass of the hot-water in the tank at time k;A res is the surface area of the tank; E res,in,k and E res,out,k are the inlet and outlet heat in the tank at time k, respectively; T out,k is the temperature of the outside of the tank at time k; θ hl is the heat loss coefficient; U res,in,min , U res,in,max , U res,out,min and U res.out,max are the allowable pressure limits at the inlet and outlet of the tank; V res,min and V res,max are the lower and upper limits of the hot-water volume, respectively.

Problem reformulation
MPPT is a nonlinear programming optimization problem which is difficult to solve directly using existing solvers because of the nonlinear differential thermodynamics (29) E res,total,k+1 = E res,total,k + E res,in,k − E res,out,k − E res,loss,k = E res,total,k + m res,k C water (T res,k+1 − T res,k ) E res,out,k = C water Q res,out,k ρ water (T res,k − T d,k+1 ) (33) m res,k+1 = m res,k + (Q PTr,k − Q res,out,k )ρ water �k (34) E res,total,min ≤ E res,total,k ≤ E res,total,max (35) U res,out,min ≤ U res,out,k ≤ U res,out,max (36) U res,in,min ≤ U res,in,k ≤ U res,in,max (37) U res,out,min ≤ U res,out,k ≤ U res,out,max (38) V res,min ρ water ≤ m res,k ≤ V res,max ρ water and hydrodynamics. In such cases, the thermodynamics-based R-C thermal network and the hydrodynamic networking model are reformulated to decrease the nonlinearity and computational burden.
As the multi-dimensional R-C thermal circuit model in (4)-(9) is highly nonlinear, the linearization method of the general nonlinear differential system in [33] is adopted to decrease the nonlinearity of the digesting thermodynamic equations. It has been demonstrated in [33] that this linearization will not result in noticeable error. This is because of the small fermentation temperature range of the biogas digester. Thus, the linearized state space form of the multi-dimensional R-C thermal model is discretized as: where T k,i = T in,k,i , T w,k,i , T s,k,i is the state vector representing the nodal temperatures of the R-C thermal model at time k; u k denotes the input vector including solar energy injection and environmental temperature; A and B are the coefficient matrices which can be obtained from the capacitances and resistances of thermal circuital branches.
In order to decrease the nonlinearity of water flows in (10)- (14), a linearization method in [5] is applied to linearize the hydrodynamic model of hot-water circulation flow. In Fig. 5  In addition, the capacitance and resistance parameters of the dynamic models can be derived from [12,22,27], and have been validated against the historical data (e.g., digestion temperature, outside temperature, thermal energy for digester heating, etc.) for parameter tuning [3]. Figure 6 illustrates the flowchart of the proposed MPPT method procedure. In this paper, an MPPT method is proposed to improve the fermentation temperature and biogas yield in extreme climates considering the thermodynamic and hydrodynamic effects, thereby enhancing the flexibility and resilience of the solar-boosted biogas generation system. Specifically, for the multi-dimensional R-C thermal circuit model, the fermentation temperature is determined by all heat changes inside the digester, including the hot water injected into the digester, outside temperature, etc. For the hydrodynamic networking model, the solar water circulation is determined by all pressure changes inside the digester, including the pressure at the throttle valve outlet, heating water from renewables, etc. Based on the reformulated thermodynamics-based R-C thermal network and the hydrodynamic networking model in Sect. 3.4, the MPC-based MPPT is adopted. The optimization problem is solved every 15 min in a 24-h rolling horizon and only the first control step is implemented at the current time k. Considering the inaccuracy of the parameters and prediction data, a feedback mechanism is executed to adjust the control decisions of the next time slot. All the tests are performed via the commercial Matlab platform on a laptop with 1.8-GHz Intel Core i7 CPU and 8 GB RAM, and solved with Yalmip solver CPLEX with their default settings. The optimized results of 24 h can be obtained and further analyzed to calculated the system performance index for comparative study.

Experimental data and settings
In this case study, the proposed MPPT method is tested on a demonstration project in Yiyang, Hunan, China.
The project is equipped with a 600m 3 cylindrical biogas digester and a 600m 3 gas tank for biogas production and storage. The actual parameters of the biogas digester, solar thermal collector, and stratified heat storage tank from the rural microgrid demonstration are shown in Table 1.
The complete scheduling cycle k and time interval ∆k of the biogas generation system are set to 24 h and 1 h respectively. Monte Carlo simulations are performed to generate a series of scenarios to capture the prediction uncertainties of solar radiation and ambient temperature. The initial scenario tree consists of 2000 scenarios, and each scenario has a 1/2000 probability of occurrence. In order to decrease the computational effort, only 10 scenarios are reserved for rolling optimization in each time slot, as shown in Fig. 7. Because of the seasonal differences in meteorological conditions, typical daily meteorological data in summer and winter in Hunan are adopted, as shown in Figs. 8 and 9.

Comparative results and analysis
In order to validate the effectiveness and superior performance of the proposed methodology, comparative studies are conducted with three schemes: 1. Scheme 1 is the proposed MPC-based MPPT method in Sects. 2 and 3. 2. Scheme 2 is the solar-boosted biogas generation system without adopting the MPPT method. 3. Scheme 3 is the biogas generation system without using the solar energy. Figures 10 and 11 depict the dynamic changes in water temperature and the volume of heated water inside the solar thermal collector in summer and winter, respectively. As can be seen, the solar collector stops working during the nighttime because of zero solar radiation, and thus the water temperature of the collector is equal to the environmental temperature and the volume of heated water is 0. During the daytime, the collector absorbs solar radiation so both the water temperature and the volume of heated water increase. It can be seen that the heated water volume follows the temperature, such as during 6 h to 14 h in winter, for full absorption of solar radiation under the protection temperature limit.
The dynamic changes of heat storage/release and the remaining heat of the stratified heat storage tank in summer and winter are illustrated in Figs. 12, 13, 14 and 15. The positive and negative values in Figs. 12 and 13 represent heat storage and heat release, respectively. It should be noted that the change of heat storage in the tank is mainly affected by the ambient temperature and solar Table 1 Technical specifications of the solar-boosted biogas generation system [5,23] Device

Power factors Power generation efficiency
Digester  Fig. 7 The daily solar radiation and ambient temperature of 10 scenarios  radiation. For example, the increment of heat storage at 13 h decreases in summer time as shown in Fig. 12 as the magnitude of the solar radiation increase is less than the temperature decrease. As the solar radiation in winter is much less than that in summer, no excess heat can be used for storage in winter. Thus, the dynamic change of the heat storage/release and the remaining heat of the storage tank in Scheme 1 are completely consistent with those in Scheme 2, as shown in Fig. 13. Furthermore, as can be seen from Figs. 12 and 14, excess thermal energy in scheme 1 is automatically stored from 12 to 16 h and released from 16 to 21 h to improve the fermentation temperature because of the nodal pressure regulation by the MPPT method. The surplus heat in scheme 2 cannot be stored in the storage tank, and thus no heat can be used for promoting the biogas yield at night, such as during 19 h to 23 h in summer. As shown in Figs. 16 and 17, the contribution of solar energy feedback to heat the digester in schemes 1 and 2, and thus to the fermentation temperature, increases biogas production compared with scheme 3. Since the MPPT method is not considered in scheme 2, the surplus thermal energy cannot be stored in the tank, and has to be directly injected into the digester. As a result, the fermentation temperature from 11 to 19 h in scheme 2 is above the optimal temperature, which leads to the decrease of biogas yield as the activity of anaerobic microorganisms is affected by the high temperature. The proposed methodology outperforms scheme 2 by maintaining the optimal fermentation temperature and biogas yield from 11 to 22 h, with the excess solar energy fully utilized and stored in hot-water. Figures 18 and 19 depict the curves of biogas fermentation temperature and yield in winter. The average daily fermentation temperature in winter is 5.37℃ in both schemes 1 and 2, which is more than twice as high as that in scheme 3. Since the fermentation temperature in scheme 3 is below 8 °C all day and below 0 °C for the first 7 h, the anaerobic bacteria can barely survive in winter, and thus the total production of biogas in scheme 3 is 0. It is clear that schemes 1 and 2 reduce the low-temperature     time. This can maintain the activity of the anaerobic bacteria. Compared with scheme 3, sufficient solar thermal energy is used to maintain the optimal fermentation temperature in schemes 1 and 2, and this enhances the biogas production even in cold climates. The 24-h comparative performance results on the average fermentation temperature, biogas yield, electricity procurement cost, gas procurement cost, and biogas supply break time with schemes 1-3 on typical summer and winter days are summarized in Tables 2  and 3, respectively. Considering that electricity is the main energy requirement of local consumers for the incomplete district heating network in the south of China, the CHP unit is running in the following-electric-load (FEL) mode with fixed electrical generation efficiency of 0.4 and thermal generation efficiency of 0.45 [3]. As shown in Figs. 20, 21, 22 and 23, the thermal and electrical outputs cannot completely meet the local demand in schemes 1-3. Hence, the CHP unit has to consume all generated biogas at every moment in the FEL mode for energy supply without excess biogas reserved, and will be shut down during biogas supply break time. Furthermore, considering a better coordination among the collector and heat storage tank under the proposed MPPT method in scheme 1, it can be seen that scheme 1 in summer can reduce the procurement cost of electricity and gas by 11.25% compared with scheme 2. The biogas supply break time in schemes 1 and 2 can also be reduced by 9 h compared to scheme 3 in winter, because of the solar energy feedback from the collector. The comparative results demonstrate the superior performance of the proposed scheme 1 for the improvements in electricity and gas procurement cost and biogas supply break time. This confirms its capability to enhance the flexibility and resilience of the solarboosted biogas generation system.

Conclusions
In this paper, a solar-boosted biogas energy generation system considering thermodynamic and hydrodynamic effects is established to control the solar energy for digester heating. A MPC-based MPPT method is developed to enhance the biogas production rate in extreme climates. The effectiveness and validity of the proposed methodology have been extensively tested on a demonstration project equipped with a 600m 3 cylindrical biogas digester in Hunan, China. The conclusions of this investigation are summarized as follows: 1. The proposed method can effectively enhance the flexibility and resilience of the solar-boosted biogas generation system, especially for the low biogas yield problem in extreme climates. 2. The proposed method can make full use of the collected solar energy and prevent excessive fermentation temperature in summer by storing the excess heat in the stratified heat storage tank. 3. The proposed method can maintain the optimal fermentation temperature and track the maximum biogas production in varying climate conditions.