Pyrolytic and Kinetic Analysis of Two Coastal Plant Species: Artemisia annua and Chenopodium glaucum

The large amount of coastal plant species available makes them ideal candidates for energy production. In this study, thermogravimetric analysis was used to evaluate the fuel properties of two coastal plant species, and the distributed activation energy model (DAEM) was employed in kinetic analysis. The major mass loss due to devolatilization started at 154 and 162°C at the heating rate of 10°C min−1 for Artemisia annua and Chenopodium glaucum, respectively. The results showed that the average activation energies of Artemisia annua and Chenopodium glaucum were 169.69 and 170.48 kJ mol−1, respectively. Furthermore, the activation energy changed while the conversion rate increased, and the frequency factor k 0 decreased greatly while the activation energy decreased. The results also indicated that the devolatilization of the two coastal plant species underwent a set of first-order reactions and could be expressed by the DAEM. Additionally, a simplified mathematical model was proposed to facilitate the prediction of devolatilization curves.


Introduction
Biomass is one of the most promising feedstocks for biofuel production because of its merits of renewability and wide distribution. Artemisia annua and Chenopodium glaucum are coastal plant species that can adapt to several stressful environmental conditions, including high saline content, drought, and heavy metal pollution. They are annual halophytes and distributed widely in North America, Europe, Africa, and Asia. Artemisia annua and Chenopodium glaucum obtained our interests because of their high yields, appropriate biomass characteristics, low input demands, and positive environmental impacts, such as helping to improve the soil texture and reduce soil salinity. The large amount of these plant species available makes them ideal candidates for energy production.
Among various energy conversion technologies, pyrolysis is one of the most commonly used techniques, which is characterized by thermal conversion of biomass into useful fuel at high temperature in the absence of oxidizing atmosphere [1]. Furthermore, thermogravimetric analysis is a highly precise method for the study of pyrolysis, and it is shown that each kind of biomass has unique pyrolysis characteristics because of its specific proportions of the components [2]. Pyrolytic and kinetic data from thermogravimetry analysis not only are useful for understanding the processes and mechanisms of the thermal degradation but also can be used as input parameters for a thermal degradation reaction system [3].
Numerous models have been used for the pyrolytic analysis. However, it has been found that distributed activation energy model (DAEM) is more accurate than some pseudomechanistic models, especially when the decomposition is carried out under inert atmosphere [4][5][6]. The present study has evaluated the fuel properties of Artemisia annua and Chenopodium glaucum. The kinetic parameters were determined by the DAEM method, and a simplified mathematical model was proposed to facilitate the prediction of devolatilization curves.

Materials. Artemisia annua and Chenopodium glaucum
were collected from a coastal zone of Yantai, Shandong province, China. The plants were oven-dried at 60 ∘ C to constant weight and then grounded with a Mini-Mill to pass through a 125 m sieve.

Proximate and Ultimate
Analysis. The moisture analysis was conducted according to ASTM E871-82 (2006). The ash content was determined according to ASTM E1755-01 (2007). The volatile matter content was analyzed according to ASTM E872-82 (2006). The fixed carbon was expressed as the 100%ash content-volatile matter-moisture content. The C, H, O, N, and S contents in the samples were measured according to our previous study [7]. All measurements were replicated three times.

Thermogravimetric Analysis.
The powders of Artemisia annua and Chenopodium glaucum were analyzed by a Mettler Toledo TGA/DSC1 STARe thermoanalyzer. The pyrolysis experiments were performed at heating rates of 5, 10, 20, and 30 ∘ C min −1 in a dynamic high purity nitrogen flow of 50 mL min −1 . The temperature of the furnace was programmed to rise from room temperature to 900 ∘ C.

Kinetic Analysis
Using DAEM. DAEM has been widely used in analyzing the thermal reaction system of biomasses [6,8]. It assumes that many irreversible first-order parallel reactions that have different rate parameters occur simultaneously. The model is expressed as where is the mass loss by time , 0 is the total mass loss at the end of pyrolysis, ( ) is the distribution function of activation energy that represents the difference in the activation energies of many first-order parallel reactions, and 0 is the frequency factor corresponding to value. After a series of reduction, including transformation and approximation [9], DAEM can be expressed as Equation (2) establishes a linear relationship between ln( / 2 ) and (1/ ) with the slope of (− / ), where is the heating rate and is the universal gas constant. Activation energy and frequency factor 0 can be determined by the slope and intercept of the Arrhenius plots.

Proximate and Ultimate Analysis.
Proximate and ultimate analysis facilitates the correlation of the individual composition to the thermal behavior of the biomass materials. Proximate analysis showed that the ash contents of Artemisia annua and Chenopodium glaucum were 4.4 and 5.1 wt%, respectively. Artemisia annua was a little lower in ash content than Chenopodium glaucum (Table 1). Ultimate analysis showed that these coastal plant species consisted of moderately high contents of carbon and oxygen and low amounts of nitrogen, hydrogen, and sulfur (Table 1).

Characterization of Pyrolysis.
The shape of the thermogravimetry (TG) curves did not change with different heating rates. However, TG curves shifted toward the right while the heating rate increased (Figure 1). This was due to the effect that the time to reach a given temperature became shorter by an increased heating rate and caused the entire curve to shift to higher temperatures. The temperatures corresponding to the maximum mass loss rate were also increased with the increasing of heating rate (Table 2). Additionally, the maximum rate of decomposition tended to increase at higher heating rate because there was more thermal energy to facilitate better heat transfer between the surroundings and the insides of the samples (Figure 2). Table 2 showed characteristics of the thermal degradation for the two coastal plant species at different heating rates. The average rate of mass loss at the same heating rate for Artemisia annua was higher than that for Chenopodium glaucum, which indicated that Artemisia annua had a little higher pyrolysis reactivity. Furthermore, the major mass loss due to devolatilization started at lower temperature for Artemisia annua than Chenopodium glaucum. At the heating rate of 10 ∘ C min −1 , for example, the onset of devolatilization occurred at 154 and 162 ∘ C for Artemisia annua and Chenopodium glaucum, respectively.

Kinetic Analysis.
The linear and parallel development for different conversion rates from 0.1 to 0.9 at various heating rates was shown in Figure 3. All plots had fairly high linear correlation coefficients ( 2 ≥ 0.98), indicating that the devolatilization of two plant species underwent a set of firstorder reactions. The results showed that Artemisia annua and Chenopodium glaucum were fairly similar in the average activation energies, which were 169.69 and 170.48 kJ mol −1 , respectively (Table 3).
While the conversion rate increased from 0.1 to 0.9, the value of Artemisia annua decreased from 180.03 to 161.45 kJ mol −1 , whereas that of Chenopodium glaucum decreased from 185.84 to 147.17 kJ mol −1 (Table 3). It can be seen that the frequency factor 0 decreased greatly while values decreased. For Artemisia annua, for example, the 0 value decreased from 2.82 × 10 12 to 1.12 × 10 8 s −1 when the value decreased from 180.03 to 161.45 kJ mol −1 . Figure 3 was available, a simplified mathematical model could be established to predict the devolatilization curves. The full line in Figure 3(a) showed that the data points obtained for each particular devolatilization rate could be linearized. And the dotted line in Figure 3(a) showed that the data points obtained for each heating rate could also be linearized, which could be described in the form of = (1/ )+ +ln (where was the slope of the line and + ln was the intercept of the line) [10]. The simplified process was exhibited in Table 4. Then the temperature at which devolatilization occurred can be determined with the intersection between the linearization    for a heating rate and for each devolatilization rate. can be described in

Prediction of Devolatilization Curves. When
The experimental devolatilization curves at the heating rate of 10 ∘ C min −1 were plotted together with the curves determined by the simplified mathematical model (Figure 4), which matched the experimental data very well for the two coastal plant species. The results were similar for the  heating rate of 5, 20 and 30 ∘ C min −1 . Moreover, the maximum relative errors between the experimental temperature and those obtained by the simplified mathematical model were 0.85% and 0.88% for Artemisia annua and Chenopodium glaucum, respectively ( Figure 5). From this validation, it was evident that the developed simplified mathematical model closely predicted the devolatilization curves of Artemisia annua and Chenopodium glaucum. Therefore, the simplified mathematical model is a useful tool for the prediction of devolatilization curves.

Conclusions
The fuel properties of Artemisia annua and Chenopodium glaucum were evaluated. The major mass loss due to devolatilization started at 154 and 162 ∘ C at the heating rate of 10 ∘ C min −1 for Artemisia annua and Chenopodium glaucum, respectively. The results showed that average activation energies of Artemisia annua and Chenopodium glaucum were 169.69 and 170.48 kJ mol −1 , respectively. The activation energy changed while the conversion rate increased, and the frequency factor 0 decreased greatly while activation energy decreased. The results also indicated that the devolatilization of the two coastal plant species underwent a set of first-order reactions and could be expressed by the DAEM. Additionally, a simplified mathematical model was proved to be a credible tool for the prediction of devolatilization curves.