725
Views
0
CrossRef citations to date
0
Altmetric
Research Article

Thermochemical behavior of agricultural and industrial sugarcane residues for bioenergy applications

, ORCID Icon, , , , & show all
Article: 2283264 | Received 05 Jul 2023, Accepted 09 Nov 2023, Published online: 20 Nov 2023

ABSTRACT

The Colombian sugarcane industry yields significant residues, categorized as agricultural and industrial. While bagasse, a widely studied industrial residue, is employed for energy recovery through combustion, agricultural residues are often left in fields. This study assesses the combustion behavior of these residues in typical collection scenarios. Additionally, it encompasses the characterization of residues from genetically modified sugarcane varieties in Colombia, potentially exhibiting distinct properties not previously documented. Non-isothermal thermogravimetrical analysis was employed to study the thermal behavior of sugarcane industrial residues (bagasse and pith) alongside agricultural residues from two different sugarcane varieties. This facilitated the determination of combustion reactivity through characteristic combustion process temperatures and technical parameters like ignition and combustion indexes. Proximate, elemental, and biochemical analyses revealed slight compositional differences. Agricultural residues demonstrated higher ash content (up to 34%) due to foreign matter adhering during harvesting, as well as soil and mud attachment during collection. Lignin content also varied, being lower for bagasse and pith, attributed to the juice extraction and milling processes that remove soluble lignin. Thermogravimetric analysis unveiled a two-stage burning process in all samples: devolatilization and char formation (~170°C), followed by char combustion (~310°C). Characteristic temperatures displayed subtle differences, with agricultural residues exhibiting lower temperatures and decomposition rates, resulting in reduced ignition and combustion indexes. This indicates heightened combustion reactivity in industrial residues, attributed to their elevated oxygen percentage, leading to more reactive functional groups and greater combustion stability compared to agricultural residues. This information is pertinent for optimizing sugarcane residues utilization in energy applications.

Research Highlights

  • Weather in collection time affects composition of sugarcane agricultural residues.

  • Combustion of sugarcane residues occurs over similar temperature ranges.

  • Industrial residues are more reactive to combustion than agricultural residues.

  • Overall thermal behavior of sugarcane residues depends on their composition.

1. Introduction

The pressing global concern over CO2 emissions from fossil fuels has spurred a surge in the adoption of renewable and sustainable fuels worldwide [Citation1]. Among these alternatives, biomass stands out as a frontrunner in the battle against global warming, operating within a CO2 fixation and release cycle [Citation2]. Recognized as a pivotal renewable energy resource, biomass plays an increasingly vital role in climate protection efforts [Citation3]. Colombia, for example, annually generates a substantial 71 million tons of agricultural residues [Citation4], showcasing significant potential for bioenergy production. This includes over 5 million tons of sugar cane bagasse (the fibrous residue left after extracting juice from sugar cane in the sugar production process), 457,000 tons of rice straw, and a total of 29 million tons of residual agricultural biomass, all ripe for harnessing in bioenergy applications [Citation5]. In the Valle del Cauca region, a prominent sugarcane-growing area, the residues from sugarcane harvesting, known as Green Harvesting Residues (GHR), represent a significant resource. GHR encompasses green leaves, dried leaves, and buds, which are commonly left in the fields after sugarcane harvesting. It yields approximately 1.82 million tons annually, with a calorific value of 16,965 kJ/kg [Citation6]. The shift from manual to mechanized sugarcane harvesting has led to an abundance of residues, necessitating efficient disposal methods, as burning is no longer a viable option [Citation7].

Coal-fired power generation currently constitutes 8.9% of Colombia’s electricity production [Citation8]. In Europe, the conversion of biomass in coal-fired power stations has yielded significant reductions in carbon emissions [Citation9]. Repurposing coal plants for biomass utilization presents a substantial opportunity to markedly reduce Colombia’s carbon footprint. In 2017, Colombian biomass-based electricity reached 793 MWh, primarily sourced from sugar mills burning bagasse [Citation5]. The deployment of low-pressure boilers in sugar mills for electricity, steam production, and residue management is a common practice [Citation10]. Notably, in 2020, cogeneration with bagasse from sugar mills in specific regions accounted for 99.8% of Colombian biomass energy production [Citation11]. By utilizing agricultural residues for combustion, the sugarcane industry can address the issue of green harvesting residues left in the fields while simultaneously reducing carbon emissions from electricity generation.

Currently, in addition to direct combustion as a strategy for exploiting these residues, there are innovative applications of sugarcane waste, such as anaerobic production of hydrogen from bagasse [Citation12], biogas production [Citation13], and even the extraction of materials extensively required by the modern industry like silica [Citation14]. These demonstrate the diverse potential of using these residues in modern technologies. Additionally, strategies like pyrolysis and gasification have been investigated [Citation15], which, along with combustion, aim to harness the energy potential of these residues. However, the energy utilization of biomass, especially from agricultural residues like GHR, poses challenges due to its sensitivity to initial composition and weather conditions during collection. Limited studies exist on GHR combustion [Citation16–21], also known as sugarcane leaves or straw in some cases. It is worth noting that these studies may not be universally applicable to all sugarcane varieties, given the genetic modifications undertaken based on specific soil conditions and production goals.

This study primarily aims to evaluate the thermal behavior of sugarcane agricultural residues, specifically Green Harvesting Residues (GHR), using non-isothermal Thermogravimetric Analysis (TGA). The focus lies in discerning differences between these agricultural residues and industrial residues like sugarcane bagasse (SB) and sugarcane pith (SP). These distinctions will be assessed in terms of proximate, elemental, and biochemical composition and their impact on thermochemical behavior during combustion. The GHR samples were collected during two predominant Colombian weather conditions, namely the rainy and dry seasons. These variations in weather significantly influence the biochemical composition of the samples, resulting in distinctive characteristics that differentiate agricultural residues from industrial residues such as bagasse. The research endeavors to offer crucial insights into the thermal behavior of different sugarcane residue samples and underscores the importance of ascertaining whether combustion can be integrated into the utilization of agricultural residues to the same extent it has been integrated into the use of industrial residues.

2. Materials and methods

2.1. Sample collection and preparation

Two samples of sugarcane green harvesting residues (mainly leaves, shoots, and short stems), GHR1 and GHR2, were collected from crop fields, while sugarcane industrial residues, bagasse (SB), and pith (SP), were obtained from a sugar mill in the Valle del Cauca region, southwest Colombia. GHR1 was collected during the dry season, and GHR2 was collected during the rainy season, both originating from indigenous sugarcane varieties developed by Cenicaña, the Colombian Sugarcane Research Center. The collection of GHR samples involved gathering residues left in the field after mechanical sugarcane harvesting, with GHR1 collected before tedding and GHR2 collected after mechanical tedding of the residues.

Subsequently, the samples were dried with atmospheric air until a constant weight was achieved. They were then quartered using a 12-riffle splitter to obtain representative samples. Finally, GHR, SB, and SP samples were milled using a Retsch ZM200 knife mill with a 4 mm screen. All samples were stored in sealed bags until further use.

2.2. Biomass characterization

Samples were reduced to <250 µm for proximate and elemental analysis using knives mill with a 0.1 mm screen. Moisture, volatile matter, ash, and fixed carbon were obtained according to ISO 18,122–2015 (ash) [Citation22], 18123–2015 (volatile matter) [Citation23], 18134-2-2015 (moisture) [Citation24], and the fixed carbon (percent dry basis) is the difference between 100 and the sum of the ash and volatile matter percentage yields, determined on a dry basis. High heating value (HHV) determination was made in a calorimetric pump IKAC5000 using 1.00–1.50 g pellets of 13 mm diameter (Mold SPECAC 13 mm), which was calibrated using benzoic acid from manufacturer, according to the standard test method ASTM D5865–13 [Citation25], originally developed for coal. Elemental composition was carried out according to ISO 16,996–2015 [Citation26] in a Leco CHN628 and Sulfur 628S. The analysis of polymeric composites was performed according to the methodology described by Ayeni et al. [Citation27] for GHR1, and the calculations of Debiagi et al. [Citation28] for GHR2, SB, and SP.

2.3. Thermogravimetric analysis

Thermal behavior of coal and biomass can be studied by thermogravimetrical analysis TGA [Citation29]. Information provided by TGA can be used to preliminarily assess combustion at large scale [Citation30]

TGA can be isothermal or dynamic (non-isothermal). Isothermal TGA is used to determine the thermal stability of the sample, decomposition rate, and gas effect over decomposition rates (such as adsorption-desorption reactions) and kinetic of decomposition reactions. On the other hand, non-isothermal TGA relates weight loss with temperature in oxidative atmosphere.

Thermal behavior was studied using a non-isothermal micro-TGA TA500 in a dried-air atmosphere using a single heating ramp of 10°C/min until 900°C. Combustion reactivity was determined using the maximum mass combustion rate Rmax,average mass combustion rate Rav, and the characteristic temperatures: Ignition temperature (Ti), peak temperature (Tp) and burnout temperature (Tb) obtained from differential thermograms. Ignition index D and the comprehensive combustion index S were calculated according to EquationEquation (1) and EquationEquation (2), respectively [Citation31–33]. Both indexes increase with higher mass combustion rates and lower characteristic temperatures, therefore an increase in S and D is referred to faster ignition and combustion at lower temperatures.

(1) D=RmaxTpTi(1)
(2) S=RmaxRavTi2Tb(2)

3. Results and discussion

In this study, the thermal behavior of industrial and agricultural sugarcane residues was evaluated. The research presents the physicochemical characterization of the residues and the combustion profiles as a starting point to establish differences and similarities in their thermal behavior. This approach allows for an assessment of whether their combustion exhibits distinctive features or shares significant similarities. The aim is to promote the utilization of agricultural residues in line with the established practice for industrial residues.

3.1. Proximate and ultimate analysis

In , a proximate analysis is presented for GHR1, GHR2, SB, and SP. It is evident that GHR2 exhibited the highest ash content. This can be attributed to the rainy conditions during sample collection, as the increased humidity facilitated the adhesion of mud and soil to the GHR surface, which was subsequently collected. The lowest ash content was observed in the bagasse, as it is a residue that undergoes industrial processing before use. Furthermore, the ash contents of GHR1 and bagasse are similar. In general, the fixed carbon content is low for GHR, reaching a maximum of 9.2% for GHR1. Nevertheless, the fixed carbon content of GHR is higher than that of SB and SP, indicating slower combustion. Residues with higher fixed carbon content are also associated with higher activation energy, although this is contingent on the distribution of polymeric compounds such as cellulose and hemicellulose [Citation34]. Regarding high heating value, all analyzed samples were found to be within the same order of magnitude, with minor variations among them. The lowest high heating value coincided with sample GHR2, whose ash concentration was the highest (34.1%).

Table 1. Biomass proximate analysis, %w/w (db)*.

The proximate analysis presented in aligns with findings from other studies that have compared the characterization of industrial and agricultural residues from sugarcane. These studies consistently indicate that industrial residues like bagasse tend to have higher volatile content, resulting in a lower proportion of fixed carbon, leading to an enhanced reactivity in combustion processes. Additionally, they highlight a higher ash content in agricultural residues, referred to as straw or trash in these studies, mentioning the challenges this represents for the integration of agricultural residues into combustion systems, particularly with regard to the slagging and fouling problems [Citation35,Citation36].

In general, it is observed that GHR1 composition presents values very close to those measured for bagasse, so that this residue could be considered as a substitute for the bagasse currently used in co-firing with coal, without incurring a drastic change in the characteristics of the raw material. The main drawback for agricultural residues utilization is the high ash content. In Colombia, ash content in the boilers can be up to 30% so GHR1 ash content would not cause operational issues. GHR1 corresponds to the residues collected during the dry season before tedding so as a first result, only GHR1 could be used for energy production without further treatment. On the other hand, for international energy purposes, according to the ISO 17,225:2014 standard [Citation37], herbaceous biomass intended for biofuel pellets should have an ash content below 10%. In this regard, only SP meets this criterion for biofuel pellet production.

Nonetheless, there are several techniques that could be evaluated to decrease the ash content of the agricultural residues collected during the rainy season, like GHR2. These techniques may involve processes such as washing, leaching, sedimentation, and fractionation, which selectively remove a particular high ash fraction of the biomass, improving its suitability for energy production and other applications [Citation38–40]. To address the high ash content in agricultural residues, this study analyzed the distribution of ash in different particle sizes. illustrates the cumulative ash content and cumulative weight distribution for samples GHR1 and GHR2. By removing particles smaller than 0.250 mm, the ash content of GHR1 decreased from 17% to 10% with an 87.3% yield. Similarly, for GHR2, the ash content decreased from 34.1% to 12.0% with an 81.7% yield. This represents a significant reduction of over 40% in ash content for the GHR collected during the rainy season.

Figure 1. Cumulative ash and weight percentage for a) GHR1 and b) GHR2.

Figure 1. Cumulative ash and weight percentage for a) GHR1 and b) GHR2.

In , elemental analysis is shown for GHR1, GHR2, SB, and SP. The carbon content was similar for all analyzed biomass samples, at around 35–40%. Coal carbon content is usually around 60–80%, which is related to higher carbonization processes compared with biomass [Citation41]. Moreover, carbon content is directly related to higher heating value as it is normally used as a variable for HHV calculations from elemental analysis [Citation42]. Hydrogen and oxygen contents are related to devolatilization behavior. In primary devolatilization, the biomass releases preferably oxygen and hydrogen by oxidation and hydrogenation reactions [Citation43]. Since all samples have a similar hydrogen content (around 5%), it is likely that the release of volatile matter is consistent among all residues with respect to hydrogenation reactions. In terms of oxygen content, it is higher for SB and SP compared with GHR samples, which indicates a higher reactivity for SB and SP since superficial oxygen compounds promote combustion by reducing stoichiometric oxygen and therefore, accelerating the process by skipping the oxygen adsorption on the surface, which is the slowest stage in combustion [Citation44,Citation45]. Finally, nitrogen and sulfur content are similar for all studied biomass samples at < 1%.

Table 2. Biomass ultimate analysis, %w/w, db×.

Molar relations H/C and O/C are indicators of the type of bonding in biomass structure and can be related with the type of polymeric compounds present. C-H and C-O bonds release a smaller amount of energy when broken than C-C bonds, so the heating value increases as the H/C and O/C ratios decrease [Citation46]. In addition, since breaking C-H and C-O bonds requires less energy than breaking C-C bonds, the higher the H/C and O/C ratios, the higher the reactivity of the fuel. Van Krevelen diagram is displayed in . It is observed that all the studied biomasses exhibit a similar atomic H/C ratio. This suggests that during the devolatilization stage, similar behaviors can be anticipated for all the samples. On the other hand, while the hydrogen content is generally in the range of 5% for most coals [Citation47] – similar to the residues under study – the higher carbon content in coal results in lower H/C ratios. This is exemplified in , where a bituminous coal (from Antioquia-Colombia) is provided as a reference for comparative purposes. This distinction underscores a greater presence of C-H bonds in biomass compared to coal.

Figure 2. Van Krevelen diagram of the biomass samples.

Figure 2. Van Krevelen diagram of the biomass samples.

Industrial residues SB and SP presented the highest values of atomic O/C ratio, followed by GHRs samples. Therefore, it would be expected that sugarcane industrial residues have a higher reactivity than the GHRs. Referenced coal, on the other hand, has lower atomic H/C and O/C ratios, which makes it less reactive but has a higher heating value when subjected to combustion.

Other authors have concurred with a similar analysis, demonstrating that biomass generally exhibits H/C and O/C ratios that can be between two to four times higher than those of coal, depending on the rank or degree of carbonization of the coal. This explains the low calorific value of biomass, but also its high reactivity in primary combustion reactions involving both hydrogen and oxygen, compared to coal [Citation48].

3.2. Biochemical composition

An experimental biochemical composition analysis was performed for GHR1, classifying the sample according to the part of the plant where it comes from. Shoots, stems, and leaves were analyzed by acid hydrolysis; this method is only suitable for biomass with low ash content (<10%); therefore, it could not be used for GHR2. This ash content limitation applies for most of the biochemical composition methods [Citation49,Citation50]. The biochemical composition, reported on a dry ash-free basis, is presented in .

Table 3. GHR1 biochemical composition %(dafb)*.

Similar to the findings of Debiagi et al. [Citation28], GHR1 demonstrates a high extractive content, exceeding 15%, a trait shared with herbaceous biomass like switchgrass, which exhibited a content of 16.99%. The lignocellulosic composition varies depending on the part of the sugarcane plant analyzed. Leaves and shoots are more similar to each other, while stems have higher lignin content and lower hemicellulose. Lignin, being the strongest biochemical compound, provides attachment for the cellulose/hemicellulose fibers of the plants. Therefore, it is more prevalent in the structurally robust parts of the plant [Citation46]. Hemicellulose, an amorphous component, exhibits variable composition depending on the type of biomass, usually falling within the range of 30–40% content. For stems, the lower hemicellulose value indicates lower reactivity compared to leaves and shoots.

The biochemical composition of GHR1 indicates that the majority of the energy derived from GHR comes from the thermochemical decomposition of hemicellulose and cellulose. As noted by Rego et al. [Citation51], lignocellulosic biomass tends to decompose as a mixture of its biochemical compounds. Therefore, it is expected to observe peaks of decomposition around 300°C and 350°C, which are the typical peak decomposition temperatures for hemicellulose and cellulose, respectively [Citation52]. The biochemical composition of the studied biomasses, including GHR samples (GHR1 and GHR2), SB, and SP, was determined using the mathematical algorithm developed by Debiagi et al. [Citation28]. The method consists of an algorithm that employs three theoretical mixtures of cellulose, hemicellulose, and lignin, each with their respective elemental composition of carbon, hydrogen, and oxygen. The resulting polymeric composition of the biomasses is then calculated as a linear combination of these theoretical mixtures, based on the elemental analysis presented in .

Considering that for GHR1 elemental analysis was performed on leaves (main constituent of collected GHR ~ 99%) a comparison between experimental and mathematical biochemical composition determination is shown in . A maximum 10.7% relative error was determined for lignin and a minimum 1.7% for hemicellulose. A thermal method such as the one described by Cano Díaz et al. [Citation53], could be explored for experimental biochemical composition determination of agricultural residues with high ash content.

Table 4. Biochemical composition determination comparison for GHR1 (leaves).

Biochemical composition for GHR, SB, and SP is very similar among them as shown in . Highest lignin content was obtained for GHR1 and GHR2, and highest hemicellulose/cellulose content was obtained for industrial residues SB and SP. These results are attributed to the additional juice extraction and milling processes, which remove soluble lignin from industrial residues [Citation54]. Similar findings have been reported by other authors when comparing the biochemical composition of these two types of sugarcane residues [Citation55]. These differences are correlated with variances in plant tissues and their processing history. This explains the lower lignin content in the bagasse, as it undergoes an extraction process, specifically warm water extraction, which is not the case for the straw.

Table 5. GHR1, GHR2, SB, and SP biochemical composition, %w/w, db×.

3.3. Adiabatic flame temperature

The adiabatic flame temperature refers to the maximum temperature that a fuel would reach if there were no energy losses in the form of heat dissipation or work. It is especially important in combustion processes since, if the adiabatic temperature is higher than the ash fusibility temperature, there is a probability that the ash will melt, causing slagging or fouling problems, which are considered serious problems in the operation of a boiler [Citation32].

An estimate based on simplifications of the complete energy balance of the combustion process was used to calculate the adiabatic flame temperature using EquationEquation (3). Tin is the initial temperature in K, LHV is the lower heating value in J/Kg, f is the mass ratio of fuel/air and Cp,gas refers to the calorific value for produced gas in J/(Kg·K). In this work, the calorific value of the gas was taken from nitrogen calorific value.

(3) Tadiab=Tin+LHVfCp,gas(3)

shows that highest adiabatic temperatures are calculated for SB. Bagasse has the highest H/C and O/C relations, then less stoichiometric oxygen is required in combustion and therefore, less nitrogen is added into the fuel/air mix. Nitrogen acts as an inert gas in the combustion reaction, absorbing energy from the reaction. In general, for all biomass samples, adiabatic flame temperatures are around 2000°C. Lowest temperatures were found for GHR samples, for which LHV is the lowest (14402 KJ/Kg).

Table 6. Calculated adiabatic temperatures (°C) for GHR, SB, and SP.

Ballesteros et al. [Citation56] calculated an adiabatic flame temperature of 772°C using 50% of excess air with sugarcane leaves. In the reported data, a higher calorific value (15600 KJ/kg), lower ash content (3.85%), higher carbon (42.94%), hydrogen (6.26%), and oxygen content (46.65%) are observed in comparison with the data reported in this work. Higher carbon and hydrogen content increase the stoichiometric oxygen required and therefore the stoichiometric air required; so then, the amount of nitrogen present in the flue gas also increases, achieving a lower adiabatic flame temperature. Cobo Barrera [Citation57], calculated the adiabatic flame temperature for agricultural cutting residues at 944°C (1218 K). In this case, the discrepancies in the results are due to the differences between the upper and lower calorific value used, since the author reports a lower calorific value of 3.36 MJ/kg for the residues, while in this work the values are around 15 MJ/kg. Meanwhile, Toscano Morales & Barriga [Citation58] calculated an adiabatic flame temperature of 1215°C for bagasse. The differences in this case are due to differences in the lower heating value used; although the upper heating value reported is similar, for the calculation of the lower heating value, the authors used bagasse with 50% moisture. This value influences the equation for calculating the adiabatic flame temperature by decreasing the energy released in combustion by 50%, so the temperature values decrease drastically compared to the adiabatic temperature obtained in this work.

3.4. Thermal analysis of fuel samples

It was observed from the TGA thermograms in and the DTG curves in that a small weight loss occurs at temperatures between 25 and 105°C for all biomass samples, which is associated with moisture evaporation. The TGA diagrams exhibit similarity across all residues. However, around 310°C, although the thermograms maintain their shape, the weight loss becomes constant first from higher to lower ash content biomass.

Figure 3. Sugarcane residues TG diagram at 10°C/min.

Figure 3. Sugarcane residues TG diagram at 10°C/min.

Figure 4. Sugarcane residues DTG at 10°C/min.

Figure 4. Sugarcane residues DTG at 10°C/min.

In DTG diagrams, three peaks and two non-changing zones can be identified. Firstly, there is a rapid moisture release occurring until 105°C. Then, between temperatures 100–170°C, the biomass absorbs energy without undergoing decomposition. After reaching the ignition temperature around 170°C, a small peak is observed in all samples, mainly associated with hemicellulose devolatilization. This is followed by the maximum decomposition peak attributed to tar and char formation through cellulose consumption [Citation59,Citation60]. Subsequently, another zone with no changes is observed around 350–450°C, during which the biomass absorbs energy. Finally, another peak is observed, corresponding to the combustion of char and tar.

The curves presented in are consistent with studies conducted by other authors who have thermally characterized agricultural and industrial residues from sugarcane [Citation35,Citation61]. These studies demonstrate mass loss profiles that indicate the progressive and differentiated degradation of biomass constituents, over the same temperature ranges. This includes an initial event related to the evaporation of water and low molecular weight components, followed by a distinct shoulder representing hemicellulose decomposition, positioned just prior to the maximum degradation peak associated with cellulose degradation. The final event marks the combustion of the formed char.

3.5. Combustion performance analysis

shows that the highest peak temperatures Tp were obtained for bagasse and pith at 320°C with the highest decomposition rates. The lowest peak temperatures were obtained for the two GHR species at 310 and 314°C with lower decomposition rates. In general, the peak temperatures obtained for the different residues are very similar to each other, with a maximum difference of 10°C between them. The high decomposition rates observed in bagasse and pith are directly correlated with their volatile matter content and biochemical composition. Previous research has demonstrated that hemicellulose and cellulose, both characterized by their high oxygen content, decompose into small molecules in the form of highly volatile gases [Citation62,Citation63]. The degradation of these components induces a structural transformation in the biomass, increasing its porosity, thereby enhancing the combustion performance [Citation64]. This finding aligns with other studies that have shown higher combustion rates and greater reactivity in biomasses with similar characteristics to bagasse and pith, in terms of their volatile matter and oxygen content, as well as their biochemical composition [Citation65,Citation66].

Table 7. Biomass reactivity for combustion.

Reactivity is typically calculated using mathematical relations that involve maximum and average decomposition rates, Rmax and Rav, as well as characteristic temperatures. First, ignition temperature (Ti) is defined as the temperature at which the decomposition rate changes after moisture release [Citation67]. Second, the peak temperature (Tp) is where the maximum decomposition rate (Rmax) occurs, and finally, burnout temperature (Tb) is the point where no further decomposition occurs [Citation68]. While many thermogravimetric analysis studies include these temperatures, there is no standardized methodology for calculating Ti and Tb. In this study, Ti was determined as the inflection point following moisture release. This calculation was done using the second derivative of the percentage weight with respect to temperature. Additionally, Tb was characterized as the temperature at which the magnitude of the combustion rate falls again below 1% per minute after combustion.

Ignition and combustion indexes enable the classification of the studied biomass from less to more reactive. SB and SP yielded higher S and D indexes, approximately twice as much as those observed for GHR samples. Biomass reactivity is influenced by physical characteristics such as particle shape and plant structure, as well as elemental composition, O/C and H/C molar ratios, and biochemical composition [Citation69]

In general, the characteristic temperatures for biomass are lower compared to coal [Citation70]. Meanwhile, decomposition rates are higher for coal which leads to higher values for both the ignition index D and the combustion index S. Bagasse characteristic temperatures Ti and Tb are similar to those reported by [Citation71], but the ignition temperature is significantly lower than the 285°C reported in that study. This difference arises from notable composition disparities. While both samples exhibit similar volatile matter content (around 79%), in the referenced study the sample has a 10% higher carbon content and an 8% lower oxygen content, which implies volatile matter with less reactive oxygen compounds. Consequently, this results in a higher ignition temperature compared to the present work. Also, as mentioned before, ignition temperature definition is not standardized. In the study of Chen et al. [Citation70], the use of a higher heating rate (20°C/min) shifts the DTG to the right, hence increasing characteristic temperatures.

The characteristic temperatures are lower for GHRs compared to SB and SP. This discrepancy can be attributed to the structural differences in the plant sources. While industrial residues, such as SB and SP, primarily consist of fibers from the stems post-juice extraction, GHR primarily comprises leaves. The observed differences in reactivity come from these structural distinctions.

Stems, which contribute to SB and SP, possess greater hardness due to the arrangement of multiple layers of robust, thick-walled cells bound together by lignin. This lignin, being less reactive compared to hemicellulose and cellulose, the primary polymeric compounds in leaves, accounts for the variation in reactivity [Citation28,Citation72,Citation73]. Although GHR exhibits higher lignin content than SB and SP, the lignin in the stems, which eventually form SB and SP residues, is more stable. This stability arises because, after juice extraction, the soluble and less stable lignin is removed due to the acidic nature of the sugarcane juice [Citation74].

In terms of ignition and combustion indexes, D and S, it is evident that GHRs exhibit lower reactivity compared to SB and SP. This discrepancy primarily comes from variations in volatile matter content, oxygen content, and fixed carbon, as these factors influence decomposition rates and characteristic temperatures. The differential thermograms of SB and SP closely resemble each other. This is because SP is a residue derived from SB, resulting in very similar characteristics between the two. For GHR1, proximate analysis indicates that the volatile matter and ash contents are very similar to those of SB. However, due to the higher O/C ratio in SB, it experiences a higher maximum decomposition peak. This is attributed to the presence of more oxygen compounds on the surface, making SB more reactive in combustion. Conversely, GHR2 possesses the lowest volatile matter content and the highest ash content among the biomasses analyzed. Consequently, it undergoes lower devolatilization compared to the other biomass samples.

The thermal behavior and composition of sugarcane residues, as discussed earlier, play a pivotal role in their potential utilization for bioenergy production. This study reveals a pivotal factor influencing the utilization of sugarcane agricultural residues for bioenergy production – the weather conditions during the collection period. This insight holds significant implications for the development of effective utilization strategies. Furthermore, the findings demonstrate that despite certain compositional differences, the overall thermal behavior of sugarcane residues remains consistent. This consistency opens up exciting possibilities for the potential integration of agricultural residues as complementary fuels, including the option of blended formulations alongside industrial residues. Through the application of pre-treatment processes to enhance their purity and reactivity, agricultural residues can potentially play a pivotal role in bioenergy production. This research not only sheds light on the seasonal variations and biochemical composition of sugarcane residues but also highlights the viability of specific strategies, such as biomass fractionation, for enhancing their utilization in the bioenergy sector. These findings mark a significant contribution to the field and open avenues for more sustainable and efficient bioenergy production practices.

4. Conclusions

This study assessed the thermochemical behavior of agricultural and industrial sugarcane residues. Results revealed comparable composition, except for agricultural residue GHR2, collected in rainy season, which showed elevated ash content, potentially leading to operational challenges. Agricultural residues displayed lower O/C ratios and higher lignin content than industrial residues, resulting in lower decomposition rates during combustion. In contrast, industrial residues exhibited higher combustion stability. Despite lower decomposition peaks in agricultural residues, combustion occurred over a similar temperature range compared to industrial residues. This provides insights for optimizing sugarcane residue utilization for bioenergy, emphasizing the significance of seasonal variations and biochemical composition.

5. Recommendations for future work

In consideration of potential future studies, several key areas warrant exploration. Firstly, conducting diverse experiments with various sugarcane varieties and employing different collection methods, across various seasons and locations, would yield more representative and reliable data regarding their properties and potential for bioenergy applications. Secondly, investigating the efficacy and feasibility of techniques such as washing and fractionation to reduce the ash content of agricultural residues collected during the rainy season represents important research to deepen. Furthermore, a comprehensive kinetic analysis of the Thermogravimetric Analysis (TGA) data is recommended. This analysis should aim to obtain crucial kinetic parameters, including activation energy, pre-exponential factor, reaction order, or frequency factor for each type of residue. Additionally, it would be valuable to conduct an economic analysis to assess the investment costs associated with collecting, preparing, and utilizing these residues for combustion, as well as exploring the potential for formulating blends with industrial residues. Finally, extending the scope of the study to explore other thermochemical conversion processes, such as pyrolysis or gasification for sugarcane residues, and subsequently comparing their performance with combustion in terms of energy yield, quality, efficiency, emissions, and other relevant factors, could offer valuable insights. These recommendations provide potential directions for further exploration, building upon the foundations laid by the current study.

Abbreviations

GHR=

green harvesting residues

TGA=

thermogravimetric analysis

SB=

sugarcane bagasse

SP=

sugarcane pith

Rmax=

maximum mass combustion rate

Rav=

average mass combustion rate

Ti=

ignition temperature

Tp=

peak temperature

Tb=

burnout temperature

D=

ignition index

S=

combustion index

HHV=

high heating value

LHV=

lower heating value

db=

dry basis

dafb=

dry ash-free basis

Tin=

initial temperature

Cp,gas=

heat capacity for produced gas

Tadiab=

adiabatic temperature

f=

mass ratio fuel/air

DTG=

differential thermogravimetric analysis

Author contribution

Karen Martinez is responsible for the acquisition, analysis of data and drafting the work. Juan Guerrero is responsible for interpretation of data, manuscript editing and review. Juan Barraza is responsible for design of the work, manuscript editing and review. Edward Lester is responsible for conceptualization and review; Carmen Forero is responsible for conceptualization and review. Orla Williams is responsible for acquisition, analysis of data and review and Nicolas Gil is responsible for conceptualization and review. All authors read, approved the final manuscript and agreed with be accountable for all aspects of the work.

Acknowledgments

The authors gratefully acknowledge the support provided by the University of Nottingham and the Royal Academy of Engineering. Colombian Sugar Cane Research Center, Cenicaña, is thanked for providing all sugarcane residues samples used in this work.

Disclosure statement

No potential conflict of interest was reported by the author(s).

Data availability statement

The data that support the findings of this study are available from the corresponding author, Juan Guerrero, upon reasonable request.

Additional information

Funding

This study was supported by the Colombia Scientific Program within the framework of the call Ecosistema Científico (Contract No. FP44842- 218-2018).

References

  • Wang F, Harindintwali J, Yuan Z. Technologies and perspectives for achieving carbon neutrality. Innovation. 2021;2(4):100180. doi: 10.1016/j.xinn.2021.100180
  • Freitas EN, Salgado JCS, Alnoch RC, et al. Challenges of biomass utilization for bioenergy in a climate change scenario. Biology. 2021;10(12):1277. doi: 10.3390/biology10121277
  • Rather RA, Wani AW, Mumtaz S, et al. Bioenergy: a foundation to environmental sustainability in a changing global climate scenario. J King Saud Univ Sci. 2022;34(1):101734. doi: 10.1016/j.jksus.2021.101734
  • Meza-Sepúlveda DC, Castro AM, Zamora A, et al. Bio-based value chains potential in the management of cacao pod waste in Colombia, a case study. Agronomy. 2021;11(4):1–15. doi: 10.3390/agronomy11040693
  • Eras JJC, Morejón MB, Gutiérrez AS, et al. A look to the electricity generation from non-conventional renewable energy sources in Colombia. Int J Energy Econ Policy. 2019;9(1):15–25. doi: 10.32479/ijeep.7108
  • Gil NJ. La ruta del RAC - Cenicaña. Primer Simposio de Energía y Medio Ambiente, una relación de doble vía. Popayán, Cauca, Colombia: Cenicaña; 2022.
  • Carvalho JLN, Nogueirol RC, Menandro LMS, et al. Agronomic and environmental implications of sugarcane straw removal: a major review. GCB Bioenergy. 2016;9(7):1181–1195. doi: 10.1111/gcbb.12410
  • UPME. Actualización Plan Energético Nacional (PEN) 2022 – 2052. Bogotá: Unidad de Planeación Minero Energética, Ministerio de Minas y Energía; 2023.
  • Beagle E, Belmont E. Comparative life cycle assessment of biomass utilization for electricity generation in the European Union and the United States. Energy Policy. 2019;128(January):267–275. doi: 10.1016/j.enpol.2019.01.006
  • Özüdoğru HMR, Nieder-Heitmann M, Haigh KF, et al. Techno-economic analysis of product biorefineries utilizing sugarcane lignocelluloses: xylitol, citric acid and glutamic acid scenarios annexed to sugar mills with electricity co-production. Ind Crops Prod. 2019;133(March):259–268. doi: 10.1016/j.indcrop.2019.03.015
  • Buriticá-Arboleda CI, Ramírez-Escobar CA, Martínez GAL, et al. Los Recursos Distribuidos De Bioenergía En Colombia. Bogotá: Universidad Distrital Francisco José de Caldas; 2020.
  • Tawfik A, Moanis R, Qyyum MA, et al. Sustainable fermentation approach for biogenic hydrogen productivity from delignified sugarcane bagasse. Int J Hydrogen Energy. 2022;47(88):37343–37358. doi: 10.1016/j.ijhydene.2021.09.200
  • Cevallos-Molina ER, Vélez-Vélez TV, Baquerizo-Crespo RJ, et al. Anaerobic treatment of sugarcane bagasse. An opportunity for sustainability in rural environments? Environ Adv. 2023;13:100427. doi: 10.1016/j.envadv.2023.100427
  • September LA, Kheswa N, Seroka NS, et al. Green synthesis of silica and silicon from agricultural residue sugarcane bagasse ash – a mini review. RSC Adv. 2023;13(2):1370–1380. doi: 10.1039/D2RA07490G
  • Fioranelli A, Bizzo WA. Generation of surplus electricity in sugarcane mills from sugarcane bagasse and straw: challenges, failures and opportunities. Renew Sust Energ Rev. 2023;186:113647. doi: 10.1016/j.rser.2023.113647
  • Colombo G, Ocampo-Duque W, Rinaldi F. Challenges in bioenergy production from sugarcane mills in developing countries: a case study. Energies. 2014;7(9):5874–5898. doi: 10.3390/en7095874
  • Go AW, Conag AT. Utilizing sugarcane leaves/straws as source of bioenergy in the Philippines: a case in the Visayas region. Renewable Energy. 2019;132:1230–1237. doi: 10.1016/j.renene.2018.09.029
  • Imran M, Anwar Khan AR. Characterization of agricultural waste sugarcane bagasse ash at 1100°C with various hours. Mater Today Proc. 2018;5(2):3346–3352. doi: 10.1016/j.matpr.2017.11.577
  • Kumar M, Sabbarwal S, Mishra PK, et al. Thermal degradation kinetics of sugarcane leaves (saccharum officinarum L) using thermo-gravimetric and differential scanning calorimetric studies. Biores Technol. 2019;279:262–270. doi: 10.1016/j.biortech.2019.01.137
  • Smithers J. Review of sugarcane trash recovery systems for energy cogeneration in South Africa. Renew Sust Energ Rev. 2014;32:915–925. doi: 10.1016/j.rser.2014.01.042
  • Zanatta ER, Reinehr TO, Awadallak JA. Kinetic studies of thermal decomposition of sugarcane bagasse and cassava bagasse. J Therm Anal Calorim. 2016;125(1):437–445. doi: 10.1007/s10973-016-5378-x
  • BS EN ISO. Solid biofuels-determination of ash content. London: BSI standards publication; 2017a.
  • BS EN ISO. BSI Standards Publication Solid biofuels — Determination of the content of volatile matter. January. 2017b.
  • BS EN ISO. BSI standards publication Solid biofuels-determination of moisture content-oven dry method part 2: total moisture-simplified method. 2015a.
  • ASTM D5865-13. Standard test method for gross calorific value of coal and coke. 2013. doi: 10.1520/D5865-12
  • BS EN ISO. PD ISO/TS 16996: 2015 BSI Standards publication Solid biofuels—determination of elemental composition by X- ray fluorescence. 2015b.
  • Ayeni AO, Adeeyo OA, Oresegun OM, et al. Compositional analysis of lignocellulosic materials: evaluation of an economically viable method suitable for woody and non-woody biomass. Am J Eng Res. 2015;44:2320–2847.
  • Debiagi PEA, Pecchi C, Gentile G, et al. Extractives extend the applicability of multistep kinetic scheme of biomass pyrolysis. Energy Fuels. 2015;29(10):6544–6555. doi: 10.1021/acs.energyfuels.5b01753
  • Guo F, He Y, Hassanpour A, et al. Thermogravimetric analysis on the co-combustion of biomass pellets with lignite and bituminous coal. Energy. 2020;197:117147. doi: 10.1016/j.energy.2020.117147
  • Alves JLF, da Silva JCG, Sellin N, et al. Upgrading of banana leaf waste to produce solid biofuel by torrefaction: physicochemical properties, combustion behaviors, and potential emissions. Environ Sci Pollut Res. 2022;29(17):25733–25747. doi: 10.1007/s11356-021-17381-x
  • Jayaraman K, Kok MV, Gokalp I. Thermogravimetric and mass spectrometric (TG-MS) analysis and kinetics of coal-biomass blends. Renewable Energy. 2017;101:293–300. doi: 10.1016/j.renene.2016.08.072
  • Protásio TDP, Guimarães Junior M, Mirmehdi S, et al. Combustion of biomass and charcoal made from babassu nutshell. Cerne. 2017;23(1):1–10. doi: 10.1590/01047760201723012202
  • Xie W, Wen S, Liu J, et al. Comparative thermogravimetric analyses of co-combustion of textile dyeing sludge and sugarcane bagasse in carbon dioxide/oxygen and nitrogen/oxygen atmospheres: thermal conversion characteristics, kinetics, and thermodynamics. Biores Technol. 2018;255:88–95. doi: 10.1016/j.biortech.2018.01.110
  • Lane DJ, Truong E, Larizza F, et al. Effect of hydrothermal carbonization on the combustion and gasification behavior of agricultural residues and macroalgae: devolatilization characteristics and Char Reactivity. Energy Fuels. 2018;32(4):4149–4159. doi: 10.1021/acs.energyfuels.7b03125
  • Bizzo WA, Lenço PC, Carvalho DJ, et al. The generation of residual biomass during the production of bio-ethanol from sugarcane, its characterization and its use in energy production. Renew Sust Energ Rev. 2014;29:589–603. doi: 10.1016/j.rser.2013.08.056
  • Camargo JMO, Gallego-Ríos JM, Neto AMP, et al. Characterization of sugarcane straw and bagasse from dry cleaning system of sugarcane for cogeneration system. Renewable Energy. 2020;158:500–508. doi: 10.1016/j.renene.2020.05.107
  • BS EN ISO. Solid biofuels-fuel specifications and classes. London: BSI standards publication; 2014.
  • Kenney KL, Smith WA, Gresham GL, et al. Understanding biomass feedstock variability. Biofuels. 2013;4(1):111–127. doi: 10.4155/bfs.12.83
  • Mantelatto PE, Soares CCSP, Carvalho DJ, et al. Industrial processing. In: Project BRA/10/G31 SUCRE sugarcane renewable electricity. Campinas SP, Brasil: Brazilian Center for Research in Energy and Materials; 2020. p. 63–78.
  • Onwudili JA, Eke JE. Effects of ash removal by agitated aqueous washing and sedimentation on the physico-chemical characteristics and fast pyrolysis of trommel fines. J Energy Inst. 2020;93(1):312–323. doi: 10.1016/j.joei.2019.02.008
  • Vassilev SV, Vassileva CG, Vassilev VS. Advantages and disadvantages of composition and properties of biomass in comparison with coal: an overview. In: Fuel. Vol. 158. Elsevier Ltd; 2015. pp. 330–350. doi: 10.1016/j.fuel.2015.05.050
  • Nhuchhen DR, Afzal MT. HHV predicting correlations for torrefied biomass using proximate and ultimate analyses. Bioengineering. 2017;4(4):7. doi: https://doi.org/10.3390/bioengineering4010007
  • Niksa S. Quantitative interpretations of primary devolatilization behavior. Process Chem Coal Util. 2020;221–270. doi: 10.1016/B978-0-12-818713-5.00006-X
  • Nelson L, Park S, Hubbe MA. Thermal depolymerization of biomass with emphasis on gasifier design and best method for catalytic hot gas conditioning. BioRes. 2018;13(2):4630–4727. doi: 10.15376/BIORES.13.2.NELSON
  • Yang G, Liu Y, Gao L, et al. Investigation of the synergistic effect and kinetic behavior of anthracite and biochar during co-combustion process in pure oxygen atmosphere. J Energy Inst. 2022;101:1–18. doi: 10.1016/j.joei.2021.12.005
  • Basu P. Biomass gasification, pyrolysis and torrefaction: practical design and theory. 3rd ed. San Diego: Academic Press, Elsevier; 2018.
  • Verma M, Loha C, Sinha AN, et al. Drying of biomass for utilising in co-firing with coal and its impact on environment – a review. Renew Sust Energ Rev. 2017;71:732–741. doi: 10.1016/j.rser.2016.12.101
  • Marx S, Laubscher ANE, Bunt JR, et al. Evaluation of sugar cane bagasse hydrothermal liquefaction products for co-gasification with coal as green coal pellet production. Bioresour Technol Rep. 2023;22:101503. doi: 10.1016/j.biteb.2023.101503
  • Karimi K, Taherzadeh MJ. A critical review on analysis in pretreatment of lignocelluloses: degree of polymerization, adsorption/desorption, and accessibility. Biores Technol. 2016;203:348–356. doi: 10.1016/j.biortech.2015.12.035
  • Sluiter A, Sluiter J, Wolfrum EJ. Methods for biomass compositional analysis. In: Behrens M, Datye AK, editors. Catalysis for the Conversion of Biomass and its derivatives. Berlin: Max-Planck-Gesellschaft zur Förderung der Wissenschaften; 2013. p. 223–225.
  • Rego F, Soares Dias AP, Casquilho M, et al. Fast determination of lignocellulosic composition of poplar biomass by thermogravimetry. Biomass Bioenergy. 2019;122(January):375–380. doi: 10.1016/j.biombioe.2019.01.037
  • Waters CL, Janupala RR, Mallinson RG, et al. Staged thermal fractionation for segregation of lignin and cellulose pyrolysis products: an experimental study of residence time and temperature effects. J Anal Appl Pyrolysis. 2017;126:380–389. doi: 10.1016/j.jaap.2017.05.008
  • Cano Díaz GS, Rosas-Aburto A, Vivaldo-Lima E, et al. Determination of the composition of lignocellulosic biomasses from combined analyses of thermal, spectroscopic, and wet chemical methods. Ind Eng Chem Res. 2021;60(9):3502–3515. doi: 10.1021/acs.iecr.0c05243
  • Río C, Lino AG, Colodette JL, et al. Differences in the chemical structure of the lignins from sugarcane bagasse and straw. Biomass Bioenergy. 2015;81:322–338. doi: 10.1016/j.biombioe.2015.07.006
  • Lima CS, Rabelo SC, Ciesielski PN, et al. Multiscale alterations in sugar cane bagasse and straw submitted to Alkaline Deacetylation. ACS Sustain Chem Eng. 2018;6(3):3796–3804. doi: 10.1021/acssuschemeng.7b04158
  • Ballesteros JI, Rueda-Ordóñez Y, Gelvez O, et al. (2018). Modelado y dimensionamiento de una hornilla de biomasa. Hal Open Science, Hal-01700614. HAL Id: hal-01700614 https://hal.archives-ouvertes.fr/hal-01700614
  • Cobo Barrera DF. Pirolisis de residuos de cosecha de caña de azúcar como alternativa de aprovechamiento en procesos de cogeneración. Cali: Universidad del Valle; 2012.
  • Toscano Morales LA, Barriga A. Análisis de los parámetros y Selección de hornos para la combustión de biomasa. (Aplicación a biomasas locales típicas). Revista Tecnológica ESPOL. 2009;1:1–9.
  • Alvarado Flores JJ, Rutiaga Quiñones JG. Estudio de cinética en procesos termogravimétricos de materiales lignocelulósicos. Maderas Cienc Tecnol. 2018;20(ahead):0–0. doi: 10.4067/S0718-221X2018005002601
  • Granados DA, Ruiz RA, Vega LY, et al. Study of reactivity reduction in sugarcane bagasse as consequence of a torrefaction process. Energy. 2017;139:818–827. doi: 10.1016/j.energy.2017.08.013
  • De Aguiar J, Bondancia TJ, Claro PIC, et al. Enzymatic deconstruction of sugarcane bagasse and straw to obtain cellulose nanomaterials. ACS Sustain Chem Eng. 2020;8(5):2287–2299. doi: 10.1021/acssuschemeng.9b06806
  • Dang H, Xu R, Zhang J, et al. Removal of oxygen-containing functional groups during hydrothermal carbonization of biomass: experimental and DFT study. Energy. 2023;276:127436. doi: 10.1016/j.energy.2023.127436
  • Riaza J, Khatami R, Levendis YA, et al. Combustion of single biomass particles in air and in oxy-fuel conditions. Biomass Bioenergy. 2014;64:162–174. doi: 10.1016/j.biombioe.2014.03.018
  • Yuan Y, He Y, Tan J, et al. Co-combustion characteristics of typical biomass and coal blends by Thermogravimetric Analysis. Front Energy Res. 2021;9:753622. doi: 10.3389/fenrg.2021.753622
  • Lu J-J, Chen W-H. Investigation on the ignition and burnout temperatures of bamboo and sugarcane bagasse by thermogravimetric analysis. Appl Energy. 2015;160:49–57. doi: 10.1016/j.apenergy.2015.09.026
  • Toptas A, Yildirim Y, Duman G, et al. Combustion behavior of different kinds of torrefied biomass and their blends with lignite. Biores Technol. 2015;177:328–336. doi: 10.1016/j.biortech.2014.11.072
  • Buratti C, Barberena M, Bartocci P, et al. Thermogravimetric analysis of the behavior of sub-bituminous coal and cellulosic ethanol residue during co-combustion. Biores Technol. 2015;186:154–162. doi: 10.1016/j.biortech.2015.03.041
  • Tang L, Xiao J, Mao Q, et al. Thermogravimetric analysis of the combustion characteristics and combustion kinetics of Coals subjected to different chemical demineralization processes. ACS Omega. 2022;7(16):13998–14008. doi: 10.1021/acsomega.2c00522
  • Gouws SM, Carrier M, Bunt JR, et al. Co-pyrolysis of coal and raw/torrefied biomass: a review on chemistry, kinetics and implementation. Renew Sust Energ Rev. 2021;135:110189. doi: 10.1016/j.rser.2020.110189
  • Chen C, Qin S, Chen F, et al. Co-combustion characteristics study of bagasse, coal and their blends by thermogravimetric analysis. J Energy Inst. 2018;8–13. doi: 10.1016/j.joei.2017.12.008
  • Galina NR, Romero Luna CM, Arce GLAF, et al. Comparative study on combustion and oxy-fuel combustion environments using mixtures of coal with sugarcane bagasse and biomass sorghum bagasse by the thermogravimetric analysis. J Energy Inst. 2019;92(3):741–754. doi: 10.1016/j.joei.2018.02.008
  • Sandhu HS, Singh MP, Gilbert RA, et al. Sugarcane botany : a brief view. In: Sandhu HS editor. Florida sugarcane handbook. Agronomy department UF/IFAS extension; 2016. pp. 1–5. https://edis.ifas.ufl.edu/
  • Wang S, Dai G, Luo Z. Lignocellulosic biomass pyrolysis mechanism: a state-of-the-art review. Prog Energ Combust Sci. 2017;83(4):37–41. doi: 10.1016/j.pecs.2017.05.004
  • Saetear P, Saechua N, Sereenonchai K. Sequential Injection System for analysis of Degree Brix, orthophosphate and pH in raw sugarcane juice applicable to sugar industry. Molecules. 2021;26(21):6484. doi: 10.3390/molecules26216484