Torrefaction Thermogravimetric Analysis and Kinetics of Sorghum Distilled Residue for Sustainable Fuel Production

This study investigated the kinetics of isothermal torrefaction of sorghum distilled residue (SDR), the main byproduct of the sorghum liquor-making process. The samples chosen were torrefied isothermally at five different temperatures under a nitrogen atmosphere in a thermogravimetric analyzer. Afterward, two different kinetic methods, the traditional model-free approach, and a two-step parallel reaction (TPR) kinetic model, were used to obtain the torrefaction kinetics of SDR. With the acquired 92–97% fit quality, which is the degree of similarity between calculated and real torrefaction curves, the traditional method approached using the Arrhenius equation showed a poor ability on kinetics prediction, whereas the TPR kinetic model optimized by the particle swarm optimization (PSO) algorithm showed that all the fit qualities are as high as 99%. The results suggest that PSO can simulate the actual torrefaction kinetics more accurately than the traditional kinetics approach. Moreover, the PSO method can be further employed for simulating the weight changes of reaction intermediates throughout the process. This computational method could be used as a powerful tool for industrial design and optimization in the biochar manufacturing process.


Introduction
Since the mid-20th century, the prosperous growth of human economic activities has caused a significant increase in greenhouse gas (GHG) emissions. The greenhouse effect continues to drive up the global average temperature rapidly, causing global warming and climate change simultaneously. These strong counterattacks from Mother Nature have badly been affecting human civilization and development. Therefore, reducing carbon emissions has become an urgent issue and critical action. Renewable energy, including biochar, is regarded as one of the sustainable and clean energies that can gradually replace traditional fossil fuels.
With the advantages of high heating value and great combustion quality as solid fuel, biochar has been widely used as bioenergy and biofuel [1]. It is an important part of renew- The evolutionary algorithm (EA) is a subset of evolutionary computation (EC) in artificial intelligence (AI). It has successively developed various technologies, to name a few, including genetic algorithm (GA), differential evolution (DE), ant colony optimization (ACO), and particle swarm optimization (PSO). PSO is a well-regarded, simple but powerful technique to optimize the pyrolysis kinetics of microalgae [22]. In this study, to achieve global optimization, PSO was used to simulate the isothermal torrefaction kinetics of SDR. For further evaluating the traditional model-free method and the evolutionary algorithm PSO, the fit qualities of the two methods were calculated to compare the degree of similarity between predicted and actual torrefaction kinetics.

Materials
Sorghum distilled residue (SDR) was the raw material selected for this study. It is a by-product of the liquor production process in a distillery located in southern Taiwan. SDR was ground into powder after drying at 105 • C for 6 h and stored in a dry-keeper until performing the thermogravimetric analysis (TGA). The analyses of various basic properties are detailed in Table 1. The proximate analysis was based on the standard procedure of the American Society for Testing and Materials. The content of crude carbohydrate, crude protein, and crude lipid in the sample was obtained by the phenol-sulfuric acid method, Kjeldahl method, and Soxhlet-extract method, respectively. The calorific value was measured by a bomb calorimeter (IKA C5000). Table 1. Basic properties of sorghum distilled residue (SDR).

Biomass SDR
Photograph Sustainability 2021, 13, x FOR PEER REVIEW 3 o kinetics of SDR for the first time, this research also compared and evaluated the traditio model-free method with the model-based PSO method, which can provide us with m inspiration when studying torrefaction kinetics of SDR. The evolutionary algorithm (EA) is a subset of evolutionary computation (EC) in tificial intelligence (AI). It has successively developed various technologies, to name a fe including genetic algorithm (GA), differential evolution (DE), ant colony optimizat (ACO), and particle swarm optimization (PSO). PSO is a well-regarded, simple but po erful technique to optimize the pyrolysis kinetics of microalgae [22]. In this study, achieve global optimization, PSO was used to simulate the isothermal torrefaction kinet of SDR. For further evaluating the traditional model-free method and the evolution algorithm PSO, the fit qualities of the two methods were calculated to compare the deg of similarity between predicted and actual torrefaction kinetics.

Materials
Sorghum distilled residue (SDR) was the raw material selected for this study. It i by-product of the liquor production process in a distillery located in southern Taiw SDR was ground into powder after drying at 105 °C for 6 h and stored in a dry-keep until performing the thermogravimetric analysis (TGA). The analyses of various ba properties are detailed in Table 1. The proximate analysis was based on the standard p cedure of the American Society for Testing and Materials. The content of crude carbo drate, crude protein, and crude lipid in the sample was obtained by the phenol-sulfu acid method, Kjeldahl method, and Soxhlet-extract method, respectively. The calor value was measured by a bomb calorimeter (IKA C5000).

Thermogravimetric Analysis
A thermogravimetric analyzer was used to analyze the SDR pyrolysis and torref tion characteristics and the torrefaction kinetics was conducted accordingly. Each ti about 5 mg of the biomass powder was weighed and loaded inside an aluminum crucib Sample weight was continuously measured under a nitrogen flow rate of 100 mL/m The sample was first heated from room temperature to 105 °C at a heating rate of

Thermogravimetric Analysis
A thermogravimetric analyzer was used to analyze the SDR pyrolysis and torrefaction characteristics and the torrefaction kinetics was conducted accordingly. Each time about 5 mg of the biomass powder was weighed and loaded inside an aluminum crucible. Sample weight was continuously measured under a nitrogen flow rate of 100 mL/min. The sample was first heated from room temperature to 105 • C at a heating rate of 20 • C/min and then held at 105 • C for 10 min to remove moisture. After that, the temperature was further increased to 850 • C at a heating rate of 20 • C/min. In this way, the distribution curves of TGA and DTG were obtained to understand the characteristics of the biomass throughout the whole pyrolysis process.
In the isothermal torrefaction experiment, the samples were heated from room temperature to 105 • C at a heating rate of 20 • C/min and then held up at 105 • C for 10 min for moisture removal. Next, they were further heated to 5 different temperatures at the same heating rate of 20 • C/min; they were 200 • C, 225 • C, 250 • C, 275 • C, and 300 • C and then subjected to isothermal torrefaction at these temperatures for 60 min, followed by heating to 850 • C at a heating rate of 20 • C/min to end the experiments. To ensure a high degree of both accuracy and precision of the experimental data, a thermogravimetric analyzer is regularly calibrated. Each of the same experimental conditions was carried out more than twice to confirm the excellent reproducibility. The relative error was controlled within 3% for the TGA measurement.

Isothermal Torrefaction Kinetics
Equation (1) is the kinetic modeling equation [5,23]. This formula can be used to macroscopically express the rate of thermal degradation to study the isothermal torrefaction kinetics.
where C is the conversion of the sample, k is the reaction rate constant and n is the reaction order. The conversion C can be defined as where W, W i , W f are the instant sample weight, initial sample weight, and final sample weight, respectively. In this study, W i represents the weight of SDR after removing moisture at 105 • C for 10 min and W f represents the weight of SDR at the final temperature of 800 • C.
In the process of isothermal torrefaction, the relationship between conversion and heating time is expressed as [24] ln 1−C 0 where C 0 is the conversion at the onset (t = t 0 ) of isothermal torrefaction. In Equation (3), when n = 1, which means that the thermal degradation is a first-order reaction, by plotting ln(1-C)-1 versus heating time (t-t 0 ), a straight line can be obtained and the slope of this straight line will be the reaction rate constant k, while n = 1, the plot of (1-C)1-n versus torrefaction time t gives a straight line with a slope of (n−1) k. Therefore, regardless of which situation it is in, by processing the data from torrefaction experiments via Equations (2) and (3), the reaction rate constant k can be obtained.
To further discuss the reaction rate constant k, assuming that the torrefaction reaction obeys the Arrhenius equation (Equation (4)), which shows that the chemical rate constant varies as a negative exponential of the reciprocal absolute temperature.
where A is the pre-exponential factor, Ea is the activation energy, R is the universal gas constant and T stands for the absolute temperature (in kelvins). Taking the logarithm of both sides in Equation (4) to get the following Equation (5).
Through the aforementioned approach, five reaction rate constants at different temperatures can be obtained. Furthermore, by plotting ln(k) versus 1/T, a straight line can be obtained with a slope of -Ea/R and an intercept of lnA. Therefore, two important constants, Ea and A, applicable to 200-300 • C can be obtained. Finally, the equation of the thermal degradation conversion in the temperature range of 200-300°C can be obtained and it can be used to describe the actual weight loss curves of SDR in the process of torrefaction.

Two-Step Reaction Mechanism
Prins et al. [6] used a two-step reaction mechanism to simulate the kinetics of torrefaction (as given in Equation (6)). First, reactant A is converted into intermediate I and the volatiles V 1 , and intermediate I is further converted into the final product B and the volatiles V 2 .
Assuming that all reactions obey first-order kinetics, the differential rate equations for each of the components are given by Equations (7) where W X is the weight of component X (X = A, I, B, V1, V2).

Particle Swarm Optimization (PSO)
In Section 2.3, the Arrhenius equation was adopted to study the kinetics in the torrefaction process. It assumes that the two important factors in the Arrhenius equation, Ea and A, are temperature-independent constants. In reality, though, Ea and A are dependent on temperature. Therefore, if the "temperature range is large", the Arrhenius equation is not suitable for describing the kinetics of chemical reaction without being modified. In addition to the Arrhenius equation, many other phenomenological formulas about activation energy have been proposed. For example, the deformed Arrhenius equation (DAE), generalized Mott law (GML), modified Arrhenius equation (MAE) and curved Arrhenius plot (CAP) all appropriately describe how Ea varies with temperature [25]. Nevertheless, the Chi-square Test and coefficient of determination found that these corrections still cannot accurately describe the real torrefaction behavior [25]. In summary, in Section 2.3, the set of Ea and A obtained in the temperature range (200-300 • C) may not accurately describe the kinetics at each isothermal torrefaction temperature. Therefore, a more precise approach is needed and the evolutionary algorithm can be used to simulate the kinetics of these five different isothermal torrefaction temperatures. The torrefaction temperatures in the simulation are not a range but five fixed values and the simulated weight loss curve of each group may be more in line with the actual weight loss curve. In this section, momentum-type PSO [26] is applied to achieve global optimization and the equation is given by Equations (12) and (13). (13) where v i is the velocity of the i-th particle, x i is the position of the i-th particle; α C , β C are momentum constants; ϕ 1 and ϕ 2 are cognitive learning rate and social learning rate, respectively; rand() is a random number in the range of [0,1]; pbest i is the best position of the i-th particle; gbest is the global best position; N particle is the number of particles searching for the optimal solution. To predict the weight loss curve, the least-squares method is applied and the target function (TAR) to be calculated is given by Equation (14) TAR where T is the temperature at which the isothermal torrefaction take place; W T j, exp is the experimental sample weight of the j-th datum at a certain reaction temperature; W T j, cal is the calculated sample weight of the j-th datum at a certain reaction temperature.
The flowchart of the PSO algorithm is shown in Figure 1. In addition, the fit quality (F) shown below is used to indicate the similarity between the calculated weight loss curve and the actual weight loss curve (Equation (15)). The PSO algorithm in this study was compiled with the Fortran programming language. First, a set of rough values, Ea and A, must be chosen simply based on relevant literature. If Ea and A were given appropriately at the beginning, it may reduce the number of iterations and have a better chance to converge. After initializing them, the program continues to process them with formulas to obtain the TAR and F values and then update a set of speed and position values of particles based on this result. Afterward, the PSO algorithm repeats the iterative cycle until the optimized criterion is met. In this study, the criterion is that TAR is smaller than 10 −5 . When the criterion is met the iteration stops and Ea and A converge as the globally best values, so the fit quality can also be found. In the process of optimizing Ea and A values, i.e., particles, the velocity and position of particle i determine the position of the next particle i + 1, which is the value of the next set of Ea and A.
where W exp max is the maximum sample weight in the experiment.
where W is the maximum sample weight in the experiment.

Basic Properties and Thermogravimetric Analysis
Excluding other components from the results of the composition analysis in Table 1, the composition could be arranged in descending order according to the amount of content and gives that carbohydrate > protein > lipid. It is speculated that the biggest chunk of the composition is derived from the starch in the endosperm of the sorghum grain. Surprisingly, after the fermentation process of the liquor production, a large amount of starch remains in the SDR without being utilized by the microorganisms. In the approximate analysis, the volatile matter is as high as 68.97%, indicating that the reactivity of the SDR is high [3]. The content of ash in the residue is 9.08% by weight. From the perspective of gasification, the slagging phenomenon in the gasifier can be reduced by blending the residues with coal [3].
The thermogravimetric analysis (TGA) and derivative thermogravimetric (DTG) analysis of SDR are plotted in Figure 2. The whole pyrolysis process can be roughly divided into 4 stages. Dehydration is the first stage (25-200 °C), the second stage is the decomposition of proteins and carbohydrates (200-350 °C) [3,27], the third stage is the decomposition of lipids (350-5 50 °C) [3,8] and, finally, the fourth stage is the decomposition of other components (550-8 00 °C). Hardly any weight change occurs after 800 °C. Obviously, the sample weight drops drastically in the second stage and the DTG curve shows the fastest thermal decomposition rate at 317 °C. It is known from Table 1 that the total content of carbohydrates and proteins exceeds 80%; thus, the depolymerization, decarboxylation, and cracking reactions [28,29] of carbohydrates and proteins [27] dominate the entire pyrolysis process of SDR.

Basic Properties and Thermogravimetric Analysis
Excluding other components from the results of the composition analysis in Table 1, the composition could be arranged in descending order according to the amount of content and gives that carbohydrate > protein > lipid. It is speculated that the biggest chunk of the composition is derived from the starch in the endosperm of the sorghum grain. Surprisingly, after the fermentation process of the liquor production, a large amount of starch remains in the SDR without being utilized by the microorganisms. In the approximate analysis, the volatile matter is as high as 68.97%, indicating that the reactivity of the SDR is high [3]. The content of ash in the residue is 9.08% by weight. From the perspective of gasification, the slagging phenomenon in the gasifier can be reduced by blending the residues with coal [3].
The thermogravimetric analysis (TGA) and derivative thermogravimetric (DTG) analysis of SDR are plotted in Figure 2. The whole pyrolysis process can be roughly divided into 4 stages. Dehydration is the first stage (25-200 • C), the second stage is the decomposition of proteins and carbohydrates (200-350 • C) [3,27], the third stage is the decomposition of lipids (350-550 • C) [3,8] and, finally, the fourth stage is the decomposition of other components (550-800 • C). Hardly any weight change occurs after 800 • C. Obviously, the sample weight drops drastically in the second stage and the DTG curve shows the fastest thermal decomposition rate at 317 • C. It is known from Table 1 that the total content of carbohydrates and proteins exceeds 80%; thus, the depolymerization, decarboxylation, and cracking reactions [28,29] of carbohydrates and proteins [27] dominate the entire pyrolysis process of SDR.

Isothermal Torrefaction
The curves of TGA and heating temperature throughout the entire SDR pyrolysis process are shown in Figure 3. In the temperature holding interval, the two curves of 200 and 225 • C almost overlap, as a consequence of light torrefaction [3] and the heating time almost does not affect the thermal decomposition of SDR. On the contrary, when torrefying at 300 • C, most of the carbohydrates and proteins in SDR decompose severely and cause the sample weight to fall to a large extent at the beginning and then after about 10 min, the sample weight decreases much slightly. In the middle torrefaction temperatures of 250 and 275 • C, the weight loss gradually increases with the torrefaction time.

Isothermal Torrefaction
The curves of TGA and heating temperature throughout the entire SDR pyrolysis process are shown in Figure 3. In the temperature holding interval, the two curves of 200 and 225 °C almost overlap, as a consequence of light torrefaction [3] and the heating time almost does not affect the thermal decomposition of SDR. On the contrary, when torrefying at 300 °C, most of the carbohydrates and proteins in SDR decompose severely and cause the sample weight to fall to a large extent at the beginning and then after about 10 min, the sample weight decreases much slightly. In the middle torrefaction temperatures of 250 and 275 °C, the weight loss gradually increases with the torrefaction time. In Figure 3, the isothermal torrefaction started at the point when the temperature reached 200, 225, 250, 275, or 300 °C. Then, these five temperatures were held for 60 min

Isothermal Torrefaction
The curves of TGA and heating temperature throughout the entire SDR pyrolysis process are shown in Figure 3. In the temperature holding interval, the two curves of 200 and 225 °C almost overlap, as a consequence of light torrefaction [3] and the heating time almost does not affect the thermal decomposition of SDR. On the contrary, when torrefying at 300 °C, most of the carbohydrates and proteins in SDR decompose severely and cause the sample weight to fall to a large extent at the beginning and then after about 10 min, the sample weight decreases much slightly. In the middle torrefaction temperatures of 250 and 275 °C, the weight loss gradually increases with the torrefaction time. In Figure 3, the isothermal torrefaction started at the point when the temperature reached 200, 225, 250, 275, or 300 °C. Then, these five temperatures were held for 60 min In Figure 3, the isothermal torrefaction started at the point when the temperature reached 200, 225, 250, 275, or 300 • C. Then, these five temperatures were held for 60 min to isothermally torrefy the SDR. Figure 4 graphs only the data for the 60 min isothermal period and plots the weight loss increment against the torrefaction time. For the torrefaction temperature below 225 • C, as the temperature time, the weight loss is not significant. For the medium-temperature torrefaction, i.e., 250 and 275 • C, obviously, the longer the torrefaction time is, the greater the weight loss. For the case of 300 • C, the weight loss increases sharply in the beginning, which indicates that most of the thermal decomposition occurs in the beginning and within a short time. to isothermally torrefy the SDR. Figure 4 graphs only the data for the 60 min isothermal period and plots the weight loss increment against the torrefaction time. For the torrefaction temperature below 225 °C, as the temperature time, the weight loss is not significant. For the medium-temperature torrefaction, i.e., 250 and 275 °C, obviously, the longer the torrefaction time is, the greater the weight loss. For the case of 300 °C, the weight loss increases sharply in the beginning, which indicates that most of the thermal decomposition occurs in the beginning and within a short time. Torrefaction severity index (TSI) is a dimensionless parameter [3], which was developed to account for the degree of torrefaction. The definition of TSI is given by Equation (16).
where ∆ stands for the weight loss increment at a certain temperature and duration; ∆ is the weight loss increment at 300 °C for 1 h. By plotting TSI, torrefaction time, and temperature, a three-dimensional diagram can be obtained as shown in Figure 5. From this three-dimensional profile, it can be pointed out that when torrefaction at the temperature below 225 °C, SDR has only a small degree of thermal decomposition. At medium temperature, TSI is relatively sensitive to the duration, so it is not quite suitable for quality control of biochar. For torrefaction at 300 °C, however, the TSI value could exceed 0.8 within 10 min. The 3D graph of TSI-temperatureduration can be used as a useful tool for SDR upgrading to biochar through torrefaction. For example, if the criterion of TSI for biochar is 0.8, the recommended operating conditions would be 300 °C and 10-15 min. Therefore, the unstable biochar quality caused by the sensitive curve in the first 10 min can be avoided and the energy cost can be reduced by ending up the process in about 15 min. Torrefaction severity index (TSI) is a dimensionless parameter [3], which was developed to account for the degree of torrefaction. The definition of TSI is given by Equation (16). (16) where ∆WI stands for the weight loss increment at a certain temperature and duration; ∆WI max is the weight loss increment at 300 • C for 1 h. By plotting TSI, torrefaction time, and temperature, a three-dimensional diagram can be obtained as shown in Figure 5. From this three-dimensional profile, it can be pointed out that when torrefaction at the temperature below 225 • C, SDR has only a small degree of thermal decomposition. At medium temperature, TSI is relatively sensitive to the duration, so it is not quite suitable for quality control of biochar. For torrefaction at 300 • C, however, the TSI value could exceed 0.8 within 10 min. The 3D graph of TSI-temperature-duration can be used as a useful tool for SDR upgrading to biochar through torrefaction. For example, if the criterion of TSI for biochar is 0.8, the recommended operating conditions would be 300 • C and 10-15 min. Therefore, the unstable biochar quality caused by the sensitive curve in the first 10 min can be avoided and the energy cost can be reduced by ending up the process in about 15 min.

Torrefaction Kinetics from the Traditional Model-Free Approach
The experimental data obtained from five different isothermal torrefaction conditions can be brought into Equation (3) to obtain the rate constant k. Afterward, ln(k) is plotted against 1/T and the coefficient of determination (R 2 ) of each n can be obtained, as shown in Table 2. When n is equal to 3, R 2 reaches a maximum, and thereby the best regression line is obtained. In the case of n = 3, Ea and A can be obtained from the slope and intercept of the straight line in Figure 6. A list of the important kinetics parameters of the torrefaction reaction of SDR in the range of 200-300 • C is detailed in Table 3.
To evaluate the traditional model-free approach, these optimal kinetics parameters were brought into Equation (3) (n = 3) to obtain the optimal weight loss curves, to compare with the actual curves (Figure 7). It shows that the curve predicted by model-free kinetics is not very consistent with the actual TGA curve. In addition, the fit qualities of five different temperatures are in the range between 92 and 97%.

Torrefaction Kinetics from the Traditional Model-Free Approach
The experimental data obtained from five different isothermal torrefaction conditions can be brought into Equation (3) to obtain the rate constant k. Afterward, ln(k) is plotted against 1/T and the coefficient of determination (R 2 ) of each n can be obtained, as shown in Table 2. When n is equal to 3, R 2 reaches a maximum, and thereby the best regression line is obtained. In the case of n = 3, Ea and A can be obtained from the slope and intercept of the straight line in Figure 6. A list of the important kinetics parameters of the torrefaction reaction of SDR in the range of 200-300 °C is detailed in Table 3.         To evaluate the traditional model-free approach, these optimal kinetics parameters were brought into Equation (3) (n = 3) to obtain the optimal weight loss curves, to compare with the actual curves (Figure 7). It shows that the curve predicted by model-free kinetics is not very consistent with the actual TGA curve. In addition, the fit qualities of five different temperatures are in the range between 92 and 97%.

Torrefaction Kinetics from PSO Approach
In this section, a two-step reaction mechanism, shown in Equation (6), was adopted and it was assumed that each step was a first-order reaction. The weight loss curves during the isothermal torrefaction process were simulated with the PSO approach and the results are plotted in Figure 8. From the figure, it is clear that the simulated weight-loss curves and the experimental weight-loss curves are very close to each other. That is, they almost overlap and have very high fit qualities (>99.25). The fit qualities of each temperature are listed in Table 4, while Table 5 details the kinetic parameters at each temperature. PSO approach was further adopted to simulate the weight change of each component in the two-step reaction mechanism and the results are shown in Figure 9, where reactant A represents SDR; reaction intermediate I represents the partially reacted SDR; product B represents the biochar after SDR upgrading; V1 and V2 represent two volatiles in step 1 and step 2 of the mechanism, respectively. The sum of curves A, I and B is equal to the simulation curve.
Comparing the Figure 9a,b, their total solid weight distributions (the sum of A, I, and B) are very similar, where the 225 • C case has a higher final biochar B, yet a smaller amount of the partial-reacted SDR (i.e., I). This could be attributed to a relatively higher extent of reaction due to the higher temperature. Figure 9c indicates that the amount of intermediate I, representing the partial-reacted SDR, is almost unchanged at a high level, which is similar to the curves I at 200 • C and 225 • C. However, the amounts of V1 and V2 are higher than those two low-temperature curves, so they decrease the final amount of biochar B. Comparing Figure 9d,e, the higher temperature of 300 • C make the reaction of step 1 completed, producing a large amount of gas V1 and intermediate solid I. Meanwhile, the reaction in step 2 is gentler and produces less biochar B. In the case of 275 • C, the temperature is milder than 300 • C.
Step 2 fully reacts, so that curve I drops at the initial stage of the reaction and a large amount of intermediate I is continuously converted into biochar B. Therefore, as far as the quality of SDR is concerned, the operating temperature of 275 • C is the best, because biochar B is very high (about 60%) and intermediate I that has not been upgraded is very little (about 2%).

Torrefaction Kinetics from PSO Approach
In this section, a two-step reaction mechanism, shown in Equation (6), was adopted and it was assumed that each step was a first-order reaction. The weight loss curves during the isothermal torrefaction process were simulated with the PSO approach and the results are plotted in Figure 8. From the figure, it is clear that the simulated weight-loss curves and the experimental weight-loss curves are very close to each other. That is, they almost overlap and have very high fit qualities (> 99. 25). The fit qualities of each temperature are listed in Table 4, while Table 5 details the kinetic parameters at each temperature.    PSO approach was further adopted to simulate the weight change of each component in the two-step reaction mechanism and the results are shown in Figure 9, where reactant A represents SDR; reaction intermediate I represents the partially reacted SDR; product B represents the biochar after SDR upgrading; V1 and V2 represent two volatiles in step 1 and step 2 of the mechanism, respectively. The sum of curves A, I and B is equal to the simulation curve. Comparing the Figure 9a,b, their total solid weight distributions (the sum of A, I, and B) are very similar, where the 225 °C case has a higher final biochar B, yet a smaller amount of the partial-reacted SDR (i.e., I). This could be attributed to a relatively higher extent of reaction due to the higher temperature. Figure 9c indicates that the amount of intermediate I, representing the partial-reacted SDR, is almost unchanged at a high level, which is similar to the curves I at 200 °C and 225 °C. However, the amounts of V1 and V2 are higher than those two low-temperature curves, so they decrease the final amount of biochar B. Comparing Figure 9d,e, the higher temperature of 300 °C make the reaction of step 1 completed, producing a large amount of gas V1 and intermediate solid I. Meanwhile, the re-

Conclusions
The results of composition analysis show that SDR remains a lot of carbohydrates and proteins that are not utilized by microorganisms. Therefore, in the TGA diagram, a violent thermal decomposition occurs around 317 • C. This could be attributed to the reactions of depolymerization, decarboxylation, and cracking of carbohydrates and proteins. The 3D TSI is a useful reference tool for upgrading SDR to biochar. The analysis reveals that the torrefaction should be held at 300 • C for 10-15 min to avoid the difficulty of quality control of biochar caused by the sensitive curve in the first 10 min and to reduce the energy cost by ending up the process in 15 min. In the Arrhenius equation, Ea and A can be regarded as constants only on an important assumption that the reaction takes place in a situation where the "reaction temperature is a constant". The optimal parameters obtained by the conventional kinetic approach can not accurately predict the torrefaction behavior and the fit qualities are at a relatively low level of 92-97%. However, the PSO algorithm can accurately simulate the real torrefaction behavior, all the fit qualities are higher than 99%. PSO algorithm is not only superior to the traditional Arrhenius equation in terms of kinetics behavior approach, it can also be further adopted to simulate the weight change distribution curves of each component in the two-step reaction mechanism. For example, the results from Table 5 show that k V1 increases with increasing temperature. Therefore, Figure 9a-e shows that V1 increases significantly with increasing temperature. When the torrefaction temperature reaches 300 • C, half of the SDR converted into gas phase V1, so that product B is at a low level of 22%. Another example is that the results in Table 5 show that the maximum value of k 2 appears at 275 • C, so a large amount of impure intermediate I is converted into product B. Figure 9d shows that I is only 2% at 275 • C. If higher biochar quality is required, the operating temperature of 275 • C is better than 300 • C, because the reaction in step 2 is relatively complete, the final SDR-biochar has a high product B yield and a low partial-reacted intermediate I content. This study shows that the PSO algorithm is a promising tool that can be used to accurately predict the torrefaction kinetics and even predict reaction intermediates that affect product quality.