Plan
Comptes Rendus

Energy analysis of a micro-cogeneration unit fed by biogas as a function of pyrolysis operating parameters
Comptes Rendus. Chimie, Volume 24 (2021) no. S1, pp. 39-55.

Résumé

In this study, the biomass degradation and the evolution of chemical species during pyrolysis are analysed with the main aim of evaluating the energy performance of a micro-cogeneration unit fed by biogas. The decomposition of the feedstock material is modelled as a two-stage process: firstly, in the reactor, the biomass is decomposed in a residual solid fraction (char) and a gaseous mixture; then, the condensable gases are divided from permanent gases generating the pyro-oil. The mathematical model proposed in this work has been developed considering the dependence of the pyrolysis process from the temperature and within the interval 500–900 °C. The kinetic of the reactions involved during the pyrolysis was also taken into account. Simulations run in AspenPlus exploiting the R-yield reactor supported by a calculator block. Afterwards, the energy recovery line for the valorisation of the pyro-products has been analysed. The gas fraction obtained at the end of the cycle was firstly characterized and then used to feed a micro-CHP system. Results are very promising, with great potential in terms of thermal recovery; more than 60% of the initially fed biogas and about 30% power output can be derived.

Métadonnées
Première publication :
Publié le :
DOI : 10.5802/crchim.77
Mots clés : Pyrolysis, Simulation, Biomass, Agricultural residue, Micro-CHP, Pyro-gas

Antonio Agrifoglio 1 ; Alberto Fichera 1 ; Antonio Gagliano 1 ; Rosaria Volpe 1

1 Department of Electrical, Electronics and Computer Engineering, University of Catania, viale Andrea Doria 6, 95125 Catania, Italy
Licence : CC-BY 4.0
Droits d'auteur : Les auteurs conservent leurs droits
@article{CRCHIM_2021__24_S1_39_0,
     author = {Antonio Agrifoglio and Alberto Fichera and Antonio Gagliano and Rosaria Volpe},
     title = {Energy analysis of a micro-cogeneration unit fed by biogas as a function of pyrolysis operating parameters},
     journal = {Comptes Rendus. Chimie},
     pages = {39--55},
     publisher = {Acad\'emie des sciences, Paris},
     volume = {24},
     number = {S1},
     year = {2021},
     doi = {10.5802/crchim.77},
     language = {en},
}
TY  - JOUR
AU  - Antonio Agrifoglio
AU  - Alberto Fichera
AU  - Antonio Gagliano
AU  - Rosaria Volpe
TI  - Energy analysis of a micro-cogeneration unit fed by biogas as a function of pyrolysis operating parameters
JO  - Comptes Rendus. Chimie
PY  - 2021
SP  - 39
EP  - 55
VL  - 24
IS  - S1
PB  - Académie des sciences, Paris
DO  - 10.5802/crchim.77
LA  - en
ID  - CRCHIM_2021__24_S1_39_0
ER  - 
%0 Journal Article
%A Antonio Agrifoglio
%A Alberto Fichera
%A Antonio Gagliano
%A Rosaria Volpe
%T Energy analysis of a micro-cogeneration unit fed by biogas as a function of pyrolysis operating parameters
%J Comptes Rendus. Chimie
%D 2021
%P 39-55
%V 24
%N S1
%I Académie des sciences, Paris
%R 10.5802/crchim.77
%G en
%F CRCHIM_2021__24_S1_39_0
Antonio Agrifoglio; Alberto Fichera; Antonio Gagliano; Rosaria Volpe. Energy analysis of a micro-cogeneration unit fed by biogas as a function of pyrolysis operating parameters. Comptes Rendus. Chimie, Volume 24 (2021) no. S1, pp. 39-55. doi : 10.5802/crchim.77. https://comptes-rendus.academie-sciences.fr/chimie/articles/10.5802/crchim.77/

Version originale du texte intégral (Proposez une traduction )

1. Introduction

Pyrolysis consists of the thermal decomposition of carbonaceous materials (generally wood or agricultural residues) into a solid, liquid and gaseous fraction in absence of oxygen [1]. Products originating from pyrolysis strongly depend on the temperature at which the pyrolytic process occurs as well as on the chemical reactions taking place during the intermediate decomposition of the feedstock material [2]. With respect to the biochemical constituents, cellulose, hemicellulose and lignin have the highest polymeric distribution in plants and, generally, in plant-derived materials. Therefore, the study of their thermal degradation is fundamental to gain awareness on the feedstock chemical composition, on its thermal stability and the decomposition kinetics of the pyrolysis process.

As studied by Ranzi et al. [3], a way to deal with the chemistry of pyrolysis consists in the characterization of the initial biomass and both primary and secondary gasification reactions. In their paper, the authors developed a model to predict the gas composition deriving from biomass gasifiers particularly focusing on primary constituents (cellulose, hemicellulose and lignin) and their degradation. In their subsequent work, they derived the yield of gas, tar and solid residue [4]. A similar contribution can be recognized in the study of Liden et al. [5] focusing on the kinetic prediction model for the liquid fraction from fast pyrolysis. On the same issue, but focussing on slow pyrolysis, Hu et al. [6] compared different kinetic models on the ground of data derived from thermogravimetric analyses (TGA). Carrier et al. [7] based their research on the application of TGA to determine the chemical structure of biomass feedstock. The study of the thermal behaviour of biomass by thermogravimetric analysis is also tackled by Zhou et al. [8] and by Zhang et al. [9]. TGA data are also used by Huang et al. [10] to predict chemical kinetics. Lee et al. [11] developed a methodology based on experimental analyses in order to derive the chemical composition of gaseous products originating from pyrolysis. Differently from the previous experimentally driven contributions, Di Blasi developed a mathematical model to couple the effects deriving from heat, momentum and mass transfer phenomena with the thermal degradation of wood [12]. A numerical study is also presented in the work of Babu and Chaurasia [13], who investigated the optimum set of parameters (temperature, concentration, time and heating rates) for the pyrolysis of biomass.

Thus, in general terms, biomass pyrolysis can be studied considering the initial feedstock composition [14], reaction pathways [15, 4] and final products composition [16]. Apart from the chemical viewpoint, pyrolysis has been extensively studied either as a stand-alone process or as an integration to gasification processes [17]. In particular, simulations can be targeted to specific pyrolysis plants [18] and including sensitivity analysis [19] or can be directed to investigate dedicated energy recovery lines [20, 21], but rarely investigations in the literature include both aspects. To fill this gap, in this paper, simulations aim at increasing awareness on both the decomposition process of biomass (including the chemical evolution of secondary reactions) and the energy performance evaluation of a micro-cogeneration unit coupled with an energy recovery line.

In Table 1, a synthetic framework for prediction-based models tackled from the aforementioned literature is offered. Each paper has been classified as a function of the operative temperature range, the feedstock biomass, modelling or numerical approaches vs. experimental investigation and the kind of analyses.

Table 1.

A synthetic summary of the literature review

PaperTemperature range Feedstock biomass Modelling approach/ experimental investigation Analyses
Ranzi et al. [3] 100–1000 °C Softwood Hardwood Experimental Gas composition prediction
Liden et al. [5] 400–700 °C Wood Experimental Organic liquid yield prediction
Hu et al. [6] 150–650 °C Cellulose Thermogravimetric experiments Biomass decomposition
Di Blasi [12] 100–1000 °C Wood Heat transfer and secondary reactions
Blanco Lopez et al. [16] 100–700 °C Olive stones Experimental Gas mixture composition
Carrier et al. [7] 100–700 °CWood and fern Thermogravimetric experiments Fractionation of biomass
Benanti et al. [19] 100–600 °C Olive pitsCHEMCAD model and experimental measurementsSimulation of char production vs. experimental data
Ducom et al. [20] 100–900 °C Olive residues Experimental Energy recovery

As emerged from the literature review, several experimental studies and models have been proposed by researchers with respect to the chemical and physical behaviours of different types of feedstock material during pyrolysis.

However, due to the intrinsic complexity of the chemical evolution of biomass, it may be useful to develop numerical simulations able to offer a simplified but reliable description of the pyro-products yields and their chemical decomposition. In this perspective and drawing lessons from the existing literature, a novel simulation model for biomass pyrolysis is presented in this study.

The proposed model includes an accurate characterization of the initial feedstock material and predicts the yield and composition of pyrolysis products in order to estimate to what extent the produced pyro-gas can be exploited for energetic purposes, i.e. to recover thermal energy to support plant operations.

In comparison to existing numerical pyrolysis models which provide a simplified composition of the produced pyrogas, the presented model proposes a very detailed analysis of the chemical species originating from the pyrolysis process.

This high granularity is a strategic detail if planning to exploit the produced biogas to fed mechanical engines, for which the exact knowledge about the chemical composition is necessary. The relationships between biomass and its primary constituents are used to accurately correlate the degradation of the feedstock material with the evolution of the chemical species originating from pyrolysis. As a further added value, the decomposition model takes into account both the dependence of degradation from temperature and the kinetics of reactions involved during the process. The model provides information not only on the evolution of the main pyro-products as a function of the temperature but also with respect to the thermal degradation of chemical species.

The paper is structured as follows. Section 2 presents the materials and methods with a focus on the feedstock composition, the modelling of the decomposition process and the characterization of final products. In addition, the simulation scheme and the modelling of the implemented energy recovery line are explained. Section 3 introduces the case study and discusses the main results of the simulations for both the primary and secondary reactions and the energy performances of the micro-cogeneration unit and the thermal recovery line.

2. Materials and methods

2.1. Preliminary consideration about the feedstock composition

Biomass properties are usually determined by taking into consideration the biochemical and elementary composition and the energy content. With respect to the biochemical constituents, as underlined by Chen et al. [15] and Yang et al. [22], the decomposition process for cellulose begins around 300 °C, with mass loss maximum peak around 340 −−360 °C. At higher temperature, typically over 400 °C, almost all cellulose has been pyrolysed leaving a solid residue of circa 7%. Hemicellulose decomposes easier than cellulose, usually between 200 −−300 °C, with mass loss peak around 270 °C. At the end of the pyrolysis process, however, a residue of circa 20% of hemicellulose can still be recognized. Lignin decomposition occurs within a large temperature range (from 100 °C to 900 °C), presenting lower conversion rates and lower mass losses compared to cellulose and hemicellulose; in this case, the solid residue from pyrolysis is around 45% of the initial feedstock.

Going into the detail of intermediate compounds from biomass pyrolysis, at the increase of the temperature, cellulose separates into water and carbon from the initial starting matrix and results in an intermediate compound called “active cellulose”, CELL-A [3, 23]. The further degradation of CELL-A results in the formation of Levoglucosan, a characteristic compound deriving from cellulose pyrolysis, and a mixture of other minor components. Hemicellulose (HCE) firstly decomposes into the two main intermediate compounds: HCE1 and HCE2. Three different sub-reactions, one of them yielding the Xylosan, characterise the HCE1 decomposition, whilst HCE2 decomposed in a mixture of minor components [3, 23]. The thermal degradation of lignin is more complex if compared to cellulose and hemicellulose. Lignin decomposition yields three main intermediate products, each of them characterized by the dominance of either carbon (C), hydrogen (H) or oxygen (O) and indicated, respectively, as LIG–C, LIG–H and LIG–O. The LIG–C decomposes in LIG–CC and a mixture of common compounds and from LIG–CC Coumaryl is obtained. Going forward, LIG–H and LIG–O interact to obtain LIG–OH from which Sinapaldehyde originates. The pyrolysis of lignin is also responsible for the production of Phenol.

To build a reliable feedstock, it is fundamental to consider the elemental composition of biomass along with its biochemical characterization. In this study, the elemental composition is derived from previous studies available in literature [24, 25, 16] and free online databases [26]. Thus, the preliminary biomass structure was obtained by rearranging the information obtained from elemental composition and the stoichiometric balances of complete biomass thermo-chemical decomposition derived from the literature [3]. From this obtained feedstock, the atomic ratios H/C and O/C ratios were calculated.

The energy content of biomass is described by the higher heating value (HHV) calculated using the empirical correlation developed by Ebeling and Jenkins [27], reported in (2.1):

HHV=0.763+0.301[C]+0.525[H]+0.064[O](2.1)
This correlation is valid for any biomass species and correlates the percentages of C, H and O derived from the ultimate analysis.

2.2. Modelling of the decomposition process

The decomposition of biomass can be modelled with respect to just the primary or to both the primary and secondary reactions [12]. The initial feedstock decomposes into solid, liquid and gaseous products; in particular, primary reactions produce intermediates that evolve in the final products after the secondary reactions [28]. In this work, a two-steps decomposition scheme has been adopted, as reported in Figure 1.

Figure 1.

Biomass decomposition scheme.

The initial feedstock decomposes into solid char and a gaseous fraction. The solid fraction (char) is mainly constituted by carbon and contains volatiles and ash; the gaseous fraction includes both condensable and non-condensable gases. Due to cooling and cleaning gas systems, the condensable gaseous fraction undergoes further secondary reactions yielding a liquid fraction (composed of water and pyrolytic oil) and a permanent gaseous fraction. Thus, in compliance with the scheme of Figure 1, the mass balance equation of the entire decomposition process can be written as:

mf=mc+mg1=mc+ml+mg2(2.2)
Being mf the mass of the inlet feedstock in the reactor, mc the mass of the solid char and mg1 the mass of the first-reaction gaseous products (including both condensable and non-condensable gases). This latter is constituted by the mass of the liquid products ml and the non-condensable gaseous fraction mg2.

Figure 2.

Software simulation logic.

To ensure the accuracy of simulations, the kinetic of the chemical decompositions should be reproduced taking into consideration all factors mostly affecting the pyrolysis. The mass changes occurring during the reactions can be calculated from the TG-Analysis through the calculation of the conversion rate 𝛼, expressed as:

𝛼=mimmimf(2.3)
in which mi and mf are the initial and final mass respectively, and m the current mass. The evolution of the species involved during the pyrolysis can be written as:
d𝛼dt=kf(T)(2.4)

The derivative of the conversion rate 𝛼 in (2.4) depends on two factors: the coefficient k, expressing the heating rate and necessary to describe the kinetics of the reaction, and the polynomial function f(T) correlating the species in the compound with the increase of temperature.

The kinetic law of any chemical reaction is governed by the Arrhenius equation as expressed in (2.5):

k=k0eEaRT(2.5)
where the coefficient k is the constant rate and describes the degree of decomposition, k0 is a pre-exponential factor, Ea is the activation energy, R = 8.314 J∕(mol⋅K) is the universal gas constant and T the absolute temperature.

The polynomial function f(T) can be determined through the percentages of each species for a temperature range of 773 < T < 1173 K:

f(T)=a+bT+cT2(2.6)

The coefficients a, b and c were determined through interpolation from the data available in literature.

2.3. Simulation scheme

The understanding of kinetics reactions underlying the pyrolysis process is a multi-faceted problem. To build a reliable simulation model, the chemical species involved during pyrolysis as well as their physical and chemical characteristics should be properly specified as input parameters within any simulation environment. Indeed, the main idea behind the simulations consists of the right choice of a layout ensuring the interaction between compounds and unit blocks [29, 30]. Calculation of mass and energy balances, thermodynamic equilibrium, rate equations, mass-flow rate, temperature, pressure are based on mathematical models and empirical correlations. These are already available in the software library or are developed by users and fed to the calculator block tools. This simulation logic is summarized in Figure 2.

The simulation of biomass decomposition is a complex issue, for which dedicated models are not available within software packages [31, 32]. A way to overcome this difficulty is to build a user-defined material, usually referred to as a nonconventional component, in order to set the composition and characteristics of the feedstock material [29]. Several examples of implemented operation blocks simulating thermochemical conversion processes have been adopted from the current literature. For instance, Benanti et al. simulated a plug flow reactor to study the correlation between temperature and composition of the final products obtained by slow pyrolysis [19]. In their study, they used the decomposition reactions of biomass to obtain the input kinetic parameters. Fonseca et al. chose a Yield Reactor block for simulating wheat straw as nonconventional material and calculated properties and heat capacity of the obtained products from fast pyrolysis [33].

2.4. Characterization of final products

The heat capacity of the gas mixture can be calculated according to (2.8):

cp,g=n=1xyncp,n(2.7)

Being yn and cp,n the mass fraction and heat capacity of each chemical species and x the number of chemical species. In detail, the heat capacities of each component have been calculated according to the NASA polynomial method [34, 35]:

cp,nR=a1+a2T+a3T2+a4T3+a5T4(2.8)
where the coefficients a15 are specifically defined for each compound and are tabulated in literature [34, 35]. Final products differ for their physical states and their chemical compositions. Therefore, to estimate their higher heating values (HHV), it is necessary to implement purpose-built methods. Regarding gaseous compounds and depending on the presence of oxygen, the HHV can be calculated according to (9) and (10), as suggested by Demirbas et al. [36]:
HHV=0.303[C]+1.423[H](2.9)
HHV=0.305[C]+1.423[H]0.154[O](2.10)

The HHV for oil can be expressed by the formula [37]:

HHV=0.325[C]+0.944[H]0.105[O](2.11)

Eventually, for both products, the lower heating value (LHV) can be expressed using (2.12), as recommended in [1]:

LHV=HHVhv9[H]100+M100(2.12)
where hv is the latent heat of vaporization and M the moisture content of gas or pyro-oil.

2.5. Recovering

One way for valorising the gas obtained from pyrolysis is to feed an engine connected to an electrical generator for the simultaneous production of electrical and thermal energy. The heat released during the combustion process is recovered from the cooling water and exhaust gases and redirected to the main plant at the service of drying and/or pyrolysis stages.

The theoretical quantity of oxygen mO 2 necessary to complete the combustion of a fuel can be expressed as [38]:

mO2=YCMO˙2MC+YH4MO2MH+YSMO2MSYO(2.13)

In this equation, Mi is the molecular mass of each compound i, being i = O2, CO2 (kgi∕kmol), Mj the molecular mass of element j, being j = C, H,S,O (kgj∕kmol) and, Yj the mass fraction of element j, with j = C,H,S,O (kgj∕kgfuel).

From these data, and knowing the oxygen percentage in air, the stoichiometric air/fuel ratio AFstoich can be calculated. However, considering complete combustion and high performances for the engine, the air/fuel equivalence ratio 𝜆 is used; in this case, the actual air/fuel ratio AF is obtained by taking into account the excess of air:

AF=𝜆AFstoich(2.14)

The recycling system is schematically represented in Figure 3.

Figure 3.

Cogeneration system cycle.

Table 2.

Biochemical distribution in the final feedstock

Component[C] (wt%)[H] (wt%)[O] (wt%)MW (g/mol)H/C (—)O/C (—)Cell (wt%)HC (wt%)Lig (wt%)HHV (MJ/kg)
Biomass 49.43 6.03 44.544130.761.453 0.67 31.40 41.58 26.02 20.13

In accordance with the specification of the CHP system, thermal energy from cooling water and exhaust gas can be recovered by a heat exchanger and recirculated [39]. The thermal energy from cooling water Qcw is calculated as:

Qcw=mcwcp,cw𝛥Tcw(2.15)

In (2.15), mcw and 𝛥Tcw are the mass of the cooling water and the temperature variation, both derived from the datasheet of the equipment; cp,cw is the heat capacity of water.

Similarly thermal energy Qex is defined as:

Qex=mexcp,ex𝛥Tex(2.16)

Being mex=mair+msyn the total mass of the exhaust gas (considering the mass of air and the mass of syngas, respectively), 𝛥Tex = Ti,exTo,ex the temperature variations and cp,ex the heat capacity of exhaust.

The mechanical power Pm generated by the system can be calculated as:

Pm=QsynQcwQexQl(2.17)
Where Qsyn=mLHVsyn is the heat generated by the syngas combustion and Ql the heat losses due to the combustion chamber, mechanical friction and lubrication. Finally, considering the alternator efficiency 𝜂gen, the produced electrical power Pel is:
Pel=Pm𝜂gen(2.18)

3. Case study and results

In this study, olive pits have been chosen as feedstock material. Their biochemical and elemental compositions have been determined by averaging data and experimental measures available from previous works of the authors [24, 40], other current literature [16] and online databases [26]. The final feedstock has been modelled as a non-conventional solid named “Biomass”, which biochemical and elemental compositions, molecular weight, H/C and O/C ratios and calorific value have been reported in Table 2.

Simulations run considering 30 kg/h of feedstock material with 15% of moisture content. The calculator block processes all components entering the reactor: in this process, the temperature is an import variable, whilst the twenty-one chemical species involved during the thermochemical conversion represent export variables. During the conversion, the evolution of each compound is analysed within the temperature range 500 −−900 °C. It is worth pointing that, in this study, the high granularity of the considered chemical species enhance the validity of the proposed simulation model and ensures a more comprehensive characterization from the chemical evolution viewpoint.

The distribution law is obtained by combining the Arrhenius equation in (2.5) and the polynomial function of the temperature expressed in (2.6). The derived equation for the mass mx of the x compound is reported in (3.1):

mx=kf(T)=k0eEaRTa+bT+cT2(3.1)
Data derived from empirical studies and literature have been adopted to develop the simulation [24, 41, 14, 42, 18]. The coefficients a, b and c used to perform the evolution of the conversion rates of each compound are obtained by comparison among empirical and literature results and are listed in Table 3.

Table 3.

Coefficients for yield distribution

Compoundabc
C2H4O2HAA7.12 − 0.012735.76 × 10−6
C2H2O2GLYOXAL0.9790.00115 − 1.62 × 10−6
C2H4OACETA − 13.980.031 − 1.6 × 10−5
C6H6O3HMFU − 30.650.0688 − 3.7 × 10−5
C3H6OPROPILEN1.61 − 0.003872.32 × 10−6
CO2C-DIOX6.678 − 0.01054.84 × 10−6
CO C-MON − 14.480.03041 − 1.4 × 10−5
CH4METHANE5.204 − 0.01187.00 × 10−6
H2HYDROG − 24.70.0482 − 2.1 × 10−5
H2OWATER10.86 − 0.016157.53 × 10−6
CH2O2FORMIC AC − 45.050.107 − 5.8 × 10−5
C CHAR12.91 − 0.022481.02 × 10−5
C6H10O5LVG − 49191109.97 − 0.05743
CH2OFORMALD − 7.910.017824 − 9.44 × 10−6
C2H6OETHANOL − 4.136130.007301 − 2.91 × 10−6
CH3OHMETHANOL − 3.090.00718 − 3.87 × 10−6
C2H4ETHYLENE4.14 − 0.016.03 × 10−6
C5H8O4XYLOSAN − 43458.31 − 0.00179
C9H10O2COUMARYL − 27.890.03749.93 × 10−6
C6H6OPHENOL − 27.70.02733.83 × 10−5
C11H12O4FE2MACR90610 − 215.70.1285

Simulations were carried out using ASPEN PLUS [29]; the plant design and the flow process used to perform the pyrolysis are reported in Figure 4.

The initial feedstock passes through the inlet (S1). The dryer (B1) removes the moisture from the feedstock and ejects it at (S3). Afterwards, both the daf (dry and ash-free) and pre-heated materials are fed through (S2) to the Yield Reactor (B2). The thermochemical conversion of biomass is simulated by coupling the reactor with a calculator block. More precisely, the chosen reactor does not require exact information on the stoichiometry or kinetic of the material; users just need to specify the production yields of the involved chemical species. It is worth noting that, being these specifications independent from the temperature, the calculator block can be usefully implemented to deal with the formation of chemical species [43]. Outlet products from the reactor are directed in the cyclone (B3) for the gas–solid separation. The gas fraction is firstly cooled in a heat exchanger (B4) and then directed to the cleaning/washing section, simulated as a distillation column (B5). A second calculator is used to set the amount of water needed to enable the condensation process in the scrubber, also taking into consideration that the gas mixture leaving the reactor depends on the temperature. A water flow is then injected from (S8) to decrease the temperature and to separate the condensable products from the permanent gas. The water flow rate w(kgh) is an export-variable of the H2O calculator and is formulated as:

w=gcp,g(Tf,gTi,g)cp,w(Tf,wTi,w)(3.2)

In this equation, Ti,w and Tf,w (°C) are, respectively, the initial and final temperature of water in (S8), i.e. at the inlet and outlet of the scrubber, and cp,w (kJ∕kg⋅K) the heat capacity of water. At the numerator, g and Ti,g are the mass flow rate and the temperature of the gas at the inlet of the line (S7). These values are defined in the calculator block as imported variables.

Figure 4.

Simulation plant design and flow process.

The yields of char, liquid and gaseous fractions as a function of the temperature are reported in Figure 5, whilst the evolution of the chemical species in the range 500 −−900 °C is illustrated in Figure 6.

Figure 5.

Yields of the solid, liquid and gaseous fractions at the increase of the temperature.

Char production has a decreasing trend up to 750 °C beyond which it becomes almost constant. Liquid production decreases slowly and continuously with temperature; this fraction needs some refinery stage to yield chemicals and biofuels.

Gas production shows a rapid increase up to 850 °C and then becomes stable at 900 °C. As said, the gas obtained from this process can be used to feed an engine for the cogeneration of electrical and thermal energy.

The yield of the three main pyrolysis products as a function of the temperature can be described by a polynomial or a linear interpolation. Interpolations provide a simple but effective way to estimate the analytical expression of any process. As can be seen from Figure 5, char and gas productions are well characterized by a quadratic function y = aT2 + bT + c, being respectively:

ychar=0.0052T20.0776T+0.4145(3.3)
ygas=0.005T2+0.0933T+0.3462(3.4)

The sign of the term a indicates either the convexity or concavity of the curve: indeed, in (3.3), as a result of the positive sign of the second-order term, the curve of the char yield is convex, whilst on the other case, in (3.4), the curve of the gas yield is concave.

Oil yield is better described by a simple linear interpolation y = mT + c, such as:

yoil=0.018T+0.2435(3.5)

The analytical expression in (3.5) implies a constant decreasing trend of oil production (also confirmed by the negative sign of the slope m) at increasing the temperature.

A common way to evaluate the effectiveness of the derived analytical expressions consists of calculating the coefficient of determination R2. Values nearby the unit implies a good/optimal prediction capability of the model, i.e. a more than acceptable goodness of fit. From this perspective, as can be read in Figure 5, the derived analytical interpolations (both the quadratic for char and gas and the linear one characterizing the oil) can be considered satisfactory.

As said, the coefficient of determination is a good indicator to confirm the effectiveness of a given analytical law; however, to be more accurate, it could be helpful to include the evaluation of the coefficient of correlation r in the analysis as it allows measuring the strength of the analytical relation between the two variables concerning the simulated dataset.

The coefficient of correlation r can be either positive or negative: if one of the considered variable increases (as in the case of the temperature in the simulation model), the other variable increases (positive correlation) or decreases (negative correlation), as well. It can be calculated as the ratio between the covariance of two variables and the product of their standard deviations.

In this case, the coefficients of linear correlation are rchar = −0.9065, roil = −0.998 and rgas = 0.9675.

It is worth pointing that the coefficient of linear correlation implies some approximations for char and gas, both characterized by a quadratic expression, remaining however highly descriptive for the scope of this study. As can be noted, the coefficient of linear correlation is negative in the case of char and oil, highlighting that the increase of the temperature produces a decrease in the yield of both components. Vice versa for the gas, for which the direction of the relationship is positive, i.e. higher temperatures permit to achieve a higher amount of pyro-gas.

As regards the prediction of the chemical species, Figure 6 depicts the evolution of the main compounds deriving from the pyrolysis of cellulose, hemicellulose and lignin at increasing temperature. The chemical species with higher percentages are reported in Figure 6(a). A drop of carbon (C) can be observed between 450 °C and 700 °C probably due to the decrease of the char production in this temperature range. Similarly, CO2 decreases, whilst H2 and CO rapidly increase becoming the two main components around 750 −−800 °C. The increase of H2 combined with the increase of CO and CH4 and a consequent drop of CO2 significantly changes the composition and properties of pyro-gas. The water content is approximately constant during the considered interval. Cellulose, hemicellulose and lignin contribute to the production of methane CH4, which becomes particularly evident at higher temperatures, i.e. in correspondence with the higher production of pyro-gas.

Figure 6.

Evolution of the chemical species.

Going into the detail of the main biochemical components of biomass, during pyrolysis, cellulose degrades into water and carbon and, as a secondary reaction, in active cellulose CELL-A from which Levoglucosan LVG and other minor components derive. As can be seen, the percentage of CELL-A converted into LVG increases at higher temperatures, reaching a maximum of around 650 −−700 °C. Beyond LVG, cellulose thermal degradation produces other volatiles such as Hydroxyacetaldehyde HAA and Glyoxal, two of the main products deriving from pyrolysis of carbohydrates in biomass and directly affecting the formation of char (coherently with the fact that lower temperatures favour charring reactions and with the trend of HAA). At higher temperatures, typically beyond 600 °C, HAA decomposes into methane CH4, carbon monoxide CO, carbon dioxide CO2, Hydroxymethylfulfural HMFU, formaldehyde CH2O, ketene C2H2O and acetaldehyde C2H4O.

The decomposition scheme for hemicellulose proposed by Ranzi et al. [4] yields the two main activated hemicellulose, HCE1 and HCE2. The first compound, HCE1, decomposes into Xylosan and ethanol C2H5OH, whose evolutions display a slight increase at increasing the temperature. In minor percentages, other components deriving from the hemicellulose sub-components are methane CH4, formaldehyde CH2O and methanol CH3OH. The evolution of these last two compounds, i.e. formaldehyde and methanol, almost imitates the trends of the chemical species originating from cellulose, such as LVG and HMFU. In this sense, from the qualitative viewpoint, the pyrolysis of cellulose and hemicellulose is quite similar and the decomposition schemes for both components with respect to the temperature can be easily compared.

Finally, following the scheme of Ranzi et al. [3], lignin decomposes into LIG–C, LIG–H and LIG–O. The first sub-component, LIG–C, mainly decomposes into LIG–CC, whose evolution produces Coumaryl, phenol C6H6O and, in smaller quantities, methane CH4, carbon monoxide CO and carbon dioxide CO2. The other two sub-components, LIG–H and LIG–O, interact to obtain LIG–OH that decomposes into Sinapaldehyde FE2MACR.

3.1. Energy recovery line

To optimize the process, an energy recovery line consisting of a micro-CHP system with output power up to 100 kWe has been taken into account [44]. Calculations are based on the REC2-40G system of Enerblu Industry [45], working with an electrical output of 43 kWe in full-load and 22 kWe in half-load. The working temperature of 800 °C is chosen to investigate to what extent the heat generated from the combustion of the gaseous fraction in a CHP system can be recovered to feed the process. Outlet products from the reactor (S4), as obtained from the first calculator block using (3.1), are reported in Table 4.

Table 4.

Products at the reactor outlet

Compound kmol/h mol% kg/h wt%
CO20.0973932.998714.28624514.28748
HCOOH0.0053850.16580.2478440.826146
HAA0.0041890.1289770.2515580.838526
H22.1692666.791134.37296714.57656
CH40.0685152.1095621.099173.6639
HMFU0.0026080.080290.3288591.096196
LVG6.07 × 10−30.1870360.9849523.283173
C2H5OH0.0018790.0578420.0865460.288487
C2H40.0097430.2999860.2733290.911095
Coumaryl0.0032850.1011440.4933271.644425
FE2MACR5.61 × 10−40.0172741.17 × 10−10.38937
Phenol0.0079790.245676 0.750942.503134
Xylosan1.03 × 10−30.0318131.37 × 10−10.455022
CH3OH0.0026120.080410.0836810.278936
CH2O0.0260430.801850.7819662.606554
C3H6O0.0045760.1409092.66 × 10−10.886008
C2H4O0.0061880.1905430.2726220.908741
Glyoxal0.0028710.0883970.1666220.555405
CO0.1879015.7854415.263186 17.54395
Char0.2977669.1681563.576465 11.92155
H2O0.34196510.529056.16060220.53534
Total flow3.247826 100 30 100

These products enter the cyclone (B3) that separates about the 97.47% of the produced char, i.e. circa 3.5 kg/h of the initial feedstock. The remaining fraction, i.e. 26.5 kg/h, corresponds to the gas mixture exiting from the cyclone and moving towards the scrubber. In this section, Equations (2.7) and (3.2) are used into the second calculator block to define the heat capacity of the gas mixture and the water mass flow rate needed to cool the gas mixture and to separate the condensable products from permanent gas. In this specific case, these calculations give as results cp,g = 3.498 kJ∕kg⋅K and w=73.5kgh.

Figure 7 reports the products distribution at the end of the cycle. As can be seen, the yield of permanent gas is significantly higher than the other two pyrolysis products and has a production rate of circa 75%; char and oil have comparable yields, around the 12%.

Tables 56 and 7 summarize the composition of the solid, liquid and gaseous fraction, respectively. Here, again, the high granularity of the chemical species can be considered as a benefit for the model, giving that such level of detail also leads to more aware choice with respect to the prime motor of the micro-CHP system.

As can be observed, the resulting chemicals from oil (Table 6) and gas (Table 7) are comparable; this, as explained by [16], can be since both gaseous and liquid products directly derive from secondary reactions of cellulose and hemicellulose, whose degradation mechanism and chemical composition can be considered similar.

Table 5.

Char production at the cyclone outlet (S6)

Compoundkmol/hkg/h
Char 0.293.48
T (°C)800

Table 6.

Oil at the phase separator outlet (S16)

Compound kmol/h mol% kg/h wt%
CO28.26 × 10−60.0214643.63 × 10−40.011451
HCOOH2.45 × 10−40.6362380.0112650.354976
HAA 0.003248.4217750.1945646.130782
H22.44 × 10−60.006354.92 × 10−60.000155
CH41.70 × 10−60.0044062.72 × 10−50.000857
HMFU2.61 × 10−36.7783613.29 × 10−110.36242
LVG6.07 × 10−315.79039.85 × 10−131.03614
C2H5OH3.50 × 10−50.0910780.0016140.050863
C2H49.12 × 10−70.0023722.56 × 10−50.000807
Coumaryl0.0032328.4006424.85 × 10−115.29314
FE2MACR5.61 × 10−41.4582621.17 × 10−13.680663
Phenol0.00748919.466570.704802 22.20854
Xylosan1.03 × 10−32.6857711.37 × 10−14.301363
CH3OH2.66 × 10−50.0691310.0008520.026852
CH2O2.10 × 10−50.0546190.0006310.019881
C3H6O1.96 × 10−50.0508661.14 × 10−30.035813
C2H4O1.60 × 10−50.0414850.0007030.022154
Glyoxal1.84 × 10−50.0479580.0010710.03374
CO2.21 × 10−60.0057356.18 × 10−50.001947
Char0.00754319.606910.0905982.854756
H2O0.00629416.35970.1133823.572702
Total0.038471 1003.173564 100
T (°C)65

At 800 °C, carbon monoxide CO is the predominant component reaching a fraction of 26%, followed by carbon dioxide CO2 and hydrogen H2 with a percentage of 21%. Before investigating the CHP energy performance, it is fundamental to conduct a preliminary analysis on the physical properties and of the energy content of the obtained gaseous fraction at the outlet of the phase separator (S12). The obtained results are summarized in Table 8.

Table 7.

Gas at the phase separator outlet (S12)

Compound kmol/h mol% kg/h wt%
CO29.74 × 10−23.5050024.28588220.59914194
HCOOH 0.005140.1850.2365791.137062126
HAA0.0009490.0341580.0569930.273925665
H22.17 × 10078.074384.37296221.01767266
CH46.85 × 10−22.4658831.0991425.282782416
HMFU5.64 × 10−92.03 × 10−77.12 × 10−73.42072 × 10−6
C2H5OH0.0018440.0663530.0849320.408206834
C2H49.74 × 10−30.3506320.2733031.313570296
Coumaryl5.32 × 10−50.0019150.007990.038401491
FE2MACR1.30 × 10−84.67 × 10−72.70 × 10−61.29981 × 10−5
Phenol 0.000490.0176444.61 × 10−20.221752071
CH3OH0.0025850.0930370.0828290.398096828
CH2O0.0260220.9365550.7813353.755314949
C3H6O4.56 × 10−30.1640090.2646661.27205847
C2H4O0.0061730.2221572.72 × 10−11.306919849
Glyoxal0.0028530.1026661.66 × 10−10.795683709
CO1.88 × 10−16.7627255.26312425.29603902
H2O0.1949887.0178863.51277116.88335526
Total2.778449 10020.80612 100
T (°C)65

Table 8.

Physical properties and energy content of the produced gas

Compoundkmol/hkg/h Volume (m3)𝜌pg (kg∕m3)cp,pg (kJ∕kg⋅K)HHV (MJ/kg)LHV (MJ/kg)
S12–Gas 2.7820.8177.100.274.06 29.04 23.92

Starting from the molar compound, the partial pressure of each compound and the total density of the permanent gas 𝜌pg (i.e. the sum of the partial densities of each gas component) can be calculated. The same approach is used to derive the heat capacity of the gas mixture before condensation, as expressed in (2.6). Finally, the higher and lower heating values have been calculated by (2.9), (2.10) and (2.12).

The available power of the obtained gas can be obtained by multiplying the mass flow rate and the lower heating value LHV [46]. According to the available initial power, the chosen REC2-40G micro-CHP system works at 95.73% of the load. In accordance with the efficiencies specified in the datasheet of the company [45], the theoretical electrical and thermal power are reported in Table 9.

Table 9.

Theoretical electrical and thermal power

Available initial powerEfficiencyFinal power
Pg = 138.27 (kW)𝜂el = 29.5%Pel = 40.79 (kWel)
𝜂th = 62.5%Pth = 86.42 (kWth)

The balances used for power generation are based on the scheme represented in Figure 3 [45]. The oxygen necessary to complete the combustion of the pyro-gas is calculated by (2.13). According to the literature, the Excess Air Coefficient 𝜆 vary between 1.5–2.5. In this study, 𝜆 = 1.8 is assumed as reference value [47, 48, 49]. Thus, the derived values used are listed in Table 10.

Mechanical power is obtained by the difference between the energy content of the pyro-gas, produced thermal power and power losses occurred during the cycle [25]. The power losses due to mechanical friction, heat losses and oil lubrication are assumed to be 5%. Thermal power coming from cooling water and exhaust gases is recovered through the two heat exchangers as depicted in Figure 4. The recovered thermal power values are

 mcwcp,cw(To,cwTi,cw) =mhccp,hc1(T2,hcT1,hc) =60.68(kWth)(3.6)
 mcwcp,cw(To,cwTi,cw) =mhccp,hc1(T2,hcT1,hc) =60.68(kWth)(3.7)

Assuming 𝜂el=91.2% as energy efficiency for the alternator, the resulting electrical power will be Pel = 40.80 (kWel).

The Sankey diagram of Figure 8 shows the normalized input and output energy fluxes resulting from the combustion of the syngas in the micro-CHP system.

Figure 7.

Products distribution at the end of the cycle.

Figure 8.

Sankey diagram for thermal energy fluxes.

Table 10.

Air necessary for a complete combustion

mO˙2 (kgO 2∕kgpyro)mair,stoich (kgair∕kgpyro)mair,p (kgair∕kgpyro)
1.949.2616.66

As can be observed from the diagram, a significant percentage of thermal recovery can be attained with respect to the initial percentage of biogas and being, moreover, the electrical power estimated to be around the 29.5%.

4. Conclusions

This research proposes a reliable and accurate model aiming at simulating the pyrolysis of biomass. The developed methodology permits to gain awareness on:

  • The creation of suitable biomass to feed the process;
  • The prediction of the production yield at the end of the cycle;
  • The composition of primary products;
  • The modelling of a micro-CHP system for the valorisation of the produced biogas.

The choice of a yield-reactor coupled with a calculator block allows operating on the evolution of products with temperature, improving the evaluation and checking stages. Results indicate that, within the range of the considered operative parameters, permanent gas grows at increasing temperature; instead, char and pyro-oil fractions decrease rapidly. Beyond the analysis of the macro-components of pyrolysis, the model is also able to offer a high granularity of the chemical species involved during the first and second reactions. This feature can be considered an added value for the proposed modelling methodology, also in view of the fact that a higher awareness of the composition of chemical species enhances the down-stream decisional process for any energy application purposes.

The proposed case study highlights that the micro-cogeneration unit effectively operates near the full-load, reaching interesting and promising results in terms of both thermal and electrical efficiency.

In particular, about 62.7% of thermal recovery can be attained with respect to the initial percentage of biogas and being, moreover, the power output estimated to be around 30%, which are very promising results.

In addition, the chosen system is also able to work in a wide range of load values (until 50%), allowing the implementation of the process even with lower input power. On the other hand, working at lower temperatures reduces the overall performances.

From the conducted analysis, the choice of the right micro-CHP system depends on the operative parameter of the plant: if it works every time at the same predetermined and standard temperature, it is possible to choose the system concerning that parameter. Otherwise, it is preferable to choose a system able to work in a wide range of input values, even with reduced efficiencies.

Further research will include the utilization of experimental data to improve the model and to highlight the dependence from other working parameters, such as the heating rate, residence time in the reactor and particle grind size, as well as a sensitivity analysis to deepen how biomass quality and types could affect the results.


Bibliographie

[1] P. Basu Biomass Gasification and Pyrolysis, Elsevier Inc., Amsterdam, 2010

[2] C. A. Zaror; D. L. Pyle Proc. Indian Acad. Sci. (Eng. Sci.), 5 (1982), pp. 269-285 | DOI

[3] E. Ranzi; A. Cuoci; T. Faravelli; A. Frassoldati; G. Migliavacca; S. Pierucci; S. Sommariva Energy Fuels, 22 (2008), pp. 4292-4300 | DOI

[4] E. Ranzi; M. Corbetta; F. Manenti; S. Pierucci Chem. Eng. Sci., 110 (2014), pp. 2-12 | DOI

[5] A. G. Liden; F. Berruti; D. S. Scott Chem. Eng. Commun., 65 (1988), pp. 207-221 | DOI

[6] S. Hu; A. Jess; M. Xu Fuel, 86 (2007), pp. 2778-2788 | DOI

[7] M. Carrier; A. Loppinet-Serani; D. Denux; J. M. Lasnier; F. Ham-Pichavant; F. Cansell; C. Aymonier Biomass Bioenergy, 32 (2011), pp. 298-307 | DOI

[8] H. Zhou; Y. Long; A. Meng; Q. Li; Y. Zhang Thermochim. Acta, 566 (2013), pp. 36-43 | DOI

[9] X. Zhang; M. Xu; R. Sun; L. Sun J. Eng. Gas Turb. Power, 128 (2006), pp. 493-496 | DOI

[10] Y. F. Huang; W. H. Kuan; P. T. Chiueh; S. L. Lo Bioresour. Technol., 102 (2011), pp. 9241-9246 | DOI

[11] D. H. Lee; H. Yang; R. Yan; D. T. Liang Fuel, 86 (2007), pp. 410-417 | DOI

[12] C. Di Blasi Combust. Sci. Technol., 90 (2007), pp. 315-340 | DOI

[13] B. V. Babu; A. S. Chaurasia Energy Convers. Manag., 44 (2003), pp. 2135-2158 | DOI

[14] P. A. Brownsort Biomass pyrolysis processes: review of scope, control and variability, 2009 (UKBRC Working paper 5, https://www.biochar.ac.uk/abstract.php?id=16)

[15] T. Chen; J. Wu; J. Zhang; J. Wu; L. Sun Bioresour. Technol., 153 (2014), pp. 223-229 | DOI

[16] M. C. Blanco Lopez; C. G. Blanco; A. Martinez-Alonso; J. M. D. Tascon J. Anal. Appl. Pyrol., 65 (2002), pp. 313-322 | DOI

[17] Y. Kalinci; A. Hepbasli; I. Dincer Int. J. Exergy, 8 (2010), pp. 99-112 | DOI

[18] J. M. Commandrè; H. Lahmidi; S. Salavdor; N. Dupassiuex Fuel Process. Technol., 92 (2011), pp. 837-844 | DOI

[19] E. Benanti; C. Freda; V. Lorefice; G. Braccio; V. Sharma Therm. Sci., 15 (2011), pp. 145-158 | DOI

[20] G. Ducom; M. Gautier; M. Pietraccini; J. P. Tagutchou; D. Lebouil; R. Gourdon Renew. Energy, 145 (2020), pp. 180-189 | DOI

[21] J. Ward; M. G. Rasul; M. M. K. Bhuiya Procedia Eng., 90 (2014), pp. 669-674 | DOI

[22] H. Yang; R. Yan; H. Chen; D. H. Lee; C. Zheng Fuel, 86 (2007), pp. 1781-1788 | DOI

[23] M. B. Gorensek; R. Shukre; C. C. Chen ACS Sustain. Chem. Eng., 7 (2019), pp. 9017-9027 | DOI

[24] A. Gagliano; F. Nocera; F. Patania; M. Bruno; D. G. Castaldo C. R. Chim., 19 (2016), pp. 441-449 | DOI

[25] A. Gagliano; F. Nocera; F. Patania; M. Detommaso; M. Bruno; S. Aneli; A. Agrifoglio Small-scale biomass gasification and micro-CHP plant for the agro-industrial sector, 9°Congresso Internazionale AIGE, Catania, 17–18 September 2015 (2015)

[26] E. B. & E. Efficiency Phyllis2 ([Online]. Available: https://phyllis.nl/Browse/Standard/ECN-Phyllis)

[27] J. M. Ebeling; B. M. Jenkins Trans. Am. Soc. Agric. Eng. (1985), pp. 898-902 | DOI

[28] W. C. Park; A. Atreya; H. R. Baum Combust. Flame, 157 (2010), pp. 481-494 | DOI

[29] AspenTech Version 8.4 of aspenONE Software Expands Activation Capability to Optimize Energy, Economics, and Heat Exchanger Designs, 2013 (available online at https://www.aspentech.com/en/resources/press-releases/version-8-4-of-aspenonert-software-expands-activation-capability-to-optimize-energy-e15032390879)

[30] Chemstations CHEMCAD Version 6. User Guide, 2010 (available online at https://www.chemstations.com/content/documents/chemcad_6_user_guide.pdf)

[31] Chemstations CC - THERM. SHELL & TUBE. User’s Guide and Tutorial, 2006 (version 5.6, available online at https://www.chemstations.com/content/documents/Archive/CCTmanual56.pdf)

[32] AspenTech Aspen Physical Property System, Physical Property Methods and Models 11.1, 2001 (avalaible online at http://web.ist.utl.pt/~ist11061/de/ASPEN/Physical_Property_Methods_and_Models.pdf)

[33] F. Gomes Fonseca; A. Funke; A. Niebel; A. P. Soares Dias; N. Dahmen J. Anal. Appl. Pyrol., 139 (2019), pp. 73-86 | DOI

[34] N. T. Wakelyn; A. G. McLain Polynomial coefficients of thermochemical data for the C–H–O–N System (1975) no. 19750007683 (vol. 72657) (Technical Memorandum (TM))

[35] B. J. McBride; S. Gordon; M. Reno Coefficients for calculating thermodynamic and transport properties of individual species (1993) (vol. 4513) (Technical Memorandum (TM))

[36] A. Demirbas; N. Ak; A. Aslan; N. Sen Pet. Sci. Technol., 36 (2018), pp. 712-717 | DOI

[37] D. Beckman; D. Elliott; B. Gevert; C. Hörnell; B. Kjellström; A. Östman; Y. Solantausta; V. Tulenheimo Techno-economic assessment of selected biomass liquefaction processes, Valtion teknillinen tutkimuskeskus. Tutkimuksia - Research Reports, 697, VTT Technical Research Centre of Finland, 1990

[38] L. M. Dion; M. Lefsrud; V. Orsat; C. Cimon Bioresource, 8 (2013), pp. 1520-1538

[39] O. Balli; H. Aras Energy Explor. Exploit., 25 (2007), pp. 39-62 | DOI

[40] A. Gagliano; F. Nocera; M. Bruno Simulation models of biomass thermochemical conversion processes, gasification and pyrolysis for the prediction of the energetic potential, Advances in Renewable Energies and Power Technologies. Volume 2: Biomass, Fuel Cells, Geothermal Energies, and Smart Grids, Elsevier, 2018, pp. 39-85 | DOI

[41] J. E. Jaimes Figueroa; Y. Camacho Ardila; B. Hoss Lunelli; R. Maciel Filho; M. R. Wolf Maciel Chem. Eng. Trans., 32 (2013), pp. 925-930

[42] V. Tihay; P. Gillard J. Anal. Appl. Pyrol., 88 (2010), pp. 168-174 | DOI

[43] A. Gagliano; F. Nocera; M. Bruno; G. Cardillo Energy Procedia, 111 (2017), pp. 1010-1019 | DOI

[44] A. V. Bridgwater Biomass Bioenergy, 38 (2012), pp. 68-94 | DOI

[45] Cogeneration, E. (http://enerblu-cogeneration.com/en/product_models/301-rec2-40-g.html) | DOI

[46] F. Y. Hagos; R. A. Aziz; S. A. Sulaiman Asian J. Sci. Res., 6 (2013), pp. 187-196 | DOI

[47] J. B. Heywood Internal Combustion Engine Fundamentals, McGraw-Hill, New York, 1988

[48] P. Quaak; H. Knoef; H. Stansen Energy from biomass. A review of combustion adn gasification technologies, World Bank Technical Papers, Energy, 422, World Bank, 1999 | DOI

[49] F. Y. Hagos; A. R. A. Aziz; S. A. Sulaiman Effect of Air-fuel Ratio on the combustion characteristics of syngas (H2:CO) in direct-injection spark_injection engine, The 6th International Conference on Applied Energy - ICAE2014, Taipei, Taiwan (Energy Procedia) (2014), pp. 2567-2571


Commentaires - Politique