Abstract

Palm oil empty fruit bunches (POEFBs) can be converted into levulinic acid (LA) and furfural, which are among the top building-block chemicals. The purpose of this study was to investigate separate and simultaneous kinetic model parameters for LA and furfural production from POEFBs, which were pretreated by soaking in aqueous ammonia (SAA). The highest LA yield, which was obtained at a reaction temperature of 170°C after 90 min in an acidic solution with a concentration of 1 M, was 52.1 mol%. The highest furfural yield was 27.94 mol%, which was obtained at a reaction temperature of 170°C after 20 min in an acidic solution with a concentration of 0.5 M. SAA pretreatment affected activation energy in glucose degradation reactions and favoured direct conversion of hemicellulose to furfural. The activation energy of LA production (EakHMF) increases with higher acid catalyst concentration, and the activation energy of furfural production (EakXYN) decreases with higher acid concentration. These trends in the activation energy occurred in both separate and simultaneous kinetic models. Simultaneous kinetic model is better to calculate kinetic parameters of LA and furfural production than separate kinetic models because the simultaneous kinetic model had a lower sum of square error (SSE) when estimating kinetic parameters.

1. Introduction

The palm oil industry continues to grow in response to increased consumption and demand for palm oil. Thus, the palm oil industry represents a major potential source of biomass. During the processing of palm oil, palm oil empty fruit bunches (POEFBs), which are a type of lignocellulosic biomass, are produced as waste. Lignocellulosic biomass consists of cellulose, hemicellulose, lignin, and minor components, such as ash, proteins, and extractives [1].

The largest component of POEFBs is cellulose. Cellulose chains dissociate to produce cellulose fibrils in wood components [2]. Cellulose is surrounded by hemicellulose and lignin, which are the other major fractions from lignocellulosic biomass [3]. Hemicellulose is a branched amorphous polymer compound, of which the main monomer is xylose [4]. Lignin consists of a complex and amorphous three-dimensional network of phenolic polymers; it acts as a support structure for plant cell walls, making these walls resistant to microbial attack [5].

Pretreatment of lignocellulosic biomass is necessary to facilitate hydrolysis of its cellulose and hemicellulose fractions. The pretreatment process is important because it affects production costs. Due to its low boiling point, ammonia can potentially be used for pretreatment of lignocellulosic biomass by the processing industry, and it can easily be recycled back into the process via ventilation [6]. In addition, ammonia is favoured by many studies, including this one, because it is inexpensive and reduces the formation of by-products [7].

Hydrolysis of lignocellulose to produce sugar monomers or other degradation compounds generally uses acid as a catalyst. Acid hydrolysis of lignocellulosic biomass produces solid and liquid fractions. The solid fraction is rich in cellulose and can be used for bioethanol production. The liquid fraction is rich in hemicellulose, hydrolysate (xylose and arabinose), and fermentation inhibitor compounds [8]. The two fermentation inhibitors produced from acid hydrolysis are levulinic acid (LA) and furfural. LA and furfural are two of the top 12 chemicals, with their potential uses as building blocks for a variety of chemicals and derivatives [9].

LA is a short-chain fatty acid with a ketone carbonyl group and an acidic carboxyl group, enabling it to produce a variety of chemical substances [10]. LA production commences with depolymerisation of the biomass cellulose fraction into oligosaccharides and glucose [11]. Next, six-carbon sugars are hydrolysed to 5-hydroxymethylfurfural (HMF), which is then dehydrated into LA and formic acid [12, 13]. LA has applications in biofuels, in which it is converted into γ-valerolactone, 2-methyltetrahydrofuran, and levulinate esters. Both γ-valerolactone (GVL) and 2-methyltetrahydrofuran can be blended directly with gasoline as alternative fuel for vehicles. Levulinate esters can be used in biodiesel because they have the same properties as fatty acid methyl esters (FAME) [14].

Biofine technology refers to a traditional technology used to produce LA from lignocellulosic biomass [15]. The technology consists of a two-stage process. In the first stage, hexose sugar is converted to HMF using a sulfuric acid catalyst (1–4%) at temperatures of 200–230°C and ambient pressure (20–25 bar) for a few seconds [12]. In the second stage, the product produced in the first stage is hydrolysed into LA at temperatures of around 190–220°C and 10–15 bar pressure for 15–30 min [12].

Furfural derived from lignocellulose has two functional groups: aldehyde and a furan ring system. Furfural can be utilized in various applications, including the production of chemicals and fuels [16]. Furfuryl alcohol is one of the most widely used furfural derivatives. Furfural production generally commences with the initial hydrolysis of the hemicellulose fraction of lignocellulose into pentose, which is then dehydrated into furfural [17]. In this process, in the batch mode, a sulfuric acid catalyst reacts with the biomass at temperatures of 170–185°C to obtain a furfural yield of approximately 40–50%. Examples of commercial furfural production processes are Quaker Oats, Westpro-modified Huaxia Technology, SupraYield, and Vedernikov [15].

Many parameters, such as temperature, acid concentration, and biomass characteristics, influence the rate of LA and furfural reactions. Previous studies have investigated the kinetic reactions of LA and furfural production in various types of biomass using sulfuric acid catalysts [1626]. Chin et al. [20] performed a kinetic study of POEFB acid hydrolysis, which produced xylose decomposition products and LA. Dussan et al. [21] performed a kinetic study of furfural and LA production from Miscanthus x giganteus. These kinetic studies used separate kinetic models; the LA kinetic model was evaluated using a cellulose degradation kinetic model, and the furfural kinetic model was evaluated using a hemicellulose degradation kinetic model. However, LA and furfural kinetic models should be evaluated using a combined kinetic model because they can be produced simultaneously from biomass in the acid hydrolysis process. Therefore, evaluation of the kinetic model of simultaneous production of furfural and levulinic acid in acid hydrolysis process is important.

There have been no studies evaluating the kinetic parameters of LA and furfural from POEFBs pretreated with ammonia in both separate and simultaneous kinetic models. Therefore, the purpose of this study was to investigate separate and simultaneous kinetic parameters of LA and furfural production using the sulfuric acid hydrolysis process from POEFBs pretreated (soaked) with aqueous ammonia.

2. Materials and Methods

2.1. Biomass and Chemicals

POEFBs were obtained from the palm oil industry PTPN 5 Kertajaya, Banten, Indonesia. The particle size of POEFBs was reduced to 20 mesh, and the POEFBs were washed, dried at temperature 100°C until the moisture content reached 1–5%, and then stored in a plastic bag until further use. Ammonia solution (25%), which was supplied by Merck, Germany, was used as the pretreatment chemical. Sulfuric acid (96.1%) from Mallinckrodt, England, was used as the acid catalyst. Standard analytical grades of glucose, xylose, HMF, furfural, and LA were purchased from Sigma-Aldrich, United States.

2.2. Soaking in Aqueous Ammonia (SAA) Pretreatment

POEFBs were soaked with aqueous ammonia solution (13.13%) for 14 hours at room temperature (27°C) [27]. The solid-to-liquid ratio was 1 : 6. After soaking for 14 hours, the pretreated POEFBs were washed with water to remove ammonia until neutrality and dried until the moisture content reached 1–5%.

2.3. Kinetic Experiments

The kinetic experiments were performed using a pressurized reactor (1 L volume, 16-bar max pressure, and 100 rpm impeller velocity). The reactants were pretreated POEFBs and 0.5 M and 1 M sulfuric acid solution with a mass ratio of 1 : 20. First, pretreated POEFBs and water were added to the reactor. Next, the reactor was tightly sealed and heated. After the reactor reached the desired temperature, a sulfuric acid catalyst was injected into the reactants. The start of the reaction time was when the sulfuric acid catalyst was released into the reactor. Duplicate samples were obtained after 10, 20, 30, 45, 90, and 120 min at 150°C, 160°C, and 170°C (120 min reaction).

2.4. Analytical Methods

Sample compositional analysis of POEFBs before and after the ammonia pretreatment was determined by referring to the method of the National Renewable Energy Laboratory (NREL) [28, 29], and scanning electron microscopy (SEM) of POEFBs was performed to determine the condition of the fibres. The liquid product concentration from kinetic experiments such as glucose, xylose, HMF, furfural, and LA was analysed using high-performance liquid chromatography (HPLC), with an Aminex HPX-87H ion exclusion column (Bio-Rad, Life Science Group Hercules, CA). The eluent was 0.006 N of H2SO4, and the flow rate was 0.6 ml/min. HPLC detector was a refractive index detector, and the temperatures of the detector and column were 60°C.

2.5. Kinetic Modeling

This research used two steps to evaluate kinetic parameters of LA and furfural production which are separate kinetic evaluation and simultaneous kinetic evaluation in the same experimental data. Separate kinetic evaluation assumed that there is no interference between LA and furfural kinetic reaction, so the kinetic parameter calculation of LA and furfural would be done in a separate kinetic model. On the other side, the simultaneous kinetic evaluation calculated LA and furfural kinetic reaction in one kinetic model. The effects of other components in POEFBs such as lignin and ash were negligible from LA and furfural kinetic models, and kinetic parameter calculation since ammonia pretreatment was done before kinetic experiments.

2.6. Separate Kinetic Model of LA and Furfural

In separate kinetic evaluation, LA production was independently evaluated using kinetic model of Girisuta et al. [22], who examined the kinetics of acid hydrolysis of sugarcane bagasse in the production of LA. The reaction model of Girisuta et al. [22] is consecutive reactions from glucan to glucose and humin, then to HMF before finally transforming to LA + FA. The model can be computed with the following kinetic equations:

As the HPLC analysis of xylose and furfural did not reveal the presence of xylose in most of the samples, xylose was not included in this furfural kinetic model. The furfural kinetic model reaction scheme for the first stage is shown in Figure 1. From Figure 1, we can compute the following kinetic equations:

2.7. Simultaneous Kinetic Model of LA and Furfural

Simultaneous kinetic evaluation would calculate LA and furfural production in one kinetic model. The reaction scheme is shown in Figure 2. From Figure 2, we can compute the following kinetic equations:

2.8. Kinetic Parameter Estimation for Separate Kinetic Model of LA and Furfural

Reaction rate constants were optimized by minimizing errors between the experimental data, such as the concentrations of glucose, xylose, HMF, furfural, and LA at various temperatures and concentrations of acid catalyst with kinetic models. The concentrations were transformed into yields to compensate the large spread in concentrations. The sum of square error (SSE) of LA and furfural model prediction was formulated:where is the model yield predicted for every sampling time (mol/mol)%, is the experimental yield for every sampling time (mol/mol)%, and is the data of glucose, xylose, HMF, furfural, and LA.

The yield of LA from the POEFBs on a molar basis () was defined as the ratio of the LA concentration in the acid hydrolysis product () to the initial C6 sugar concentration in the POEFBs ():

The yield of furfural from the POEFBs on a molar basis () was defined as the ratio of the furfural in the acid hydrolysis product () to the initial C5 sugar concentration in the POEFBs ():

Minimization of the sum of square error values was done by fminsearch optimization using MATLAB optimization routine to get LA and furfural reaction rate constants at optimum slope. The value of the reaction rate constants was used to determine the value of the preexponential factor () and activation energy (), using the Arrhenius equation. Systematically, the relationship between the reaction rate constant () and temperature () was expressed by the Arrhenius equation:where is the reaction rate constant, is the preexponential factor, is the activation energy (kJ/mol), is the gas constant (kJ/mol·K), and is the temperature (K).

The effect of acid concentration on LA and furfural production can be predicted using acid reaction order equation that in this equation, there is correlation between reaction rate constant and acid concentration, which can be formulated:where is the reaction rate constant, is the preexponential factor from acid reaction order equation, is the acid concentration (M), and is the acid reaction order. The calculation of acid reaction order was counted based on (8) from the average value of kinetic experiment temperature.

2.9. Kinetic Parameter Estimation for Simultaneous Kinetic Model of LA and Furfural

The calculation of kinetic parameters for simultaneous kinetic model of LA and furfural was same with the separate kinetic model, but the yield from experimental data was obtained from mass basis.

The yield of LA and furfural from the POEFBs on mass basis was defined as the ratio of the LA and furfural mass concentration in the acid hydrolysis product ( and ) to the mass of cellulose and hemicellulose in the POEFBs ():

3. Results and Discussion

3.1. SAA Pretreatment

The purpose of the pretreatment process was to improve the conversion of cellulose and hemicellulose to LA and furfural. The results of the compositional analysis of different types of biomass with SAA pretreatment are shown in Table 1.

According to Table 1, the composition of the POEFBs before and after the SAA pretreatment did not differ greatly, with just a small reduction in lignin composition (from 22.8% to 21.8%). The composition of cellulose and hemicellulose in the untreated POEFBs increased slightly from 39.3% to 29.8%, as compared to the SAA pretreated POEFBs (42% and 32.0%, resp.). As shown in Table 1, previous studies also reported insignificant changes in lignin composition in untreated and SAA pretreated samples at various temperatures, soaking times, and ammonia concentrations [3033].

Lignin can affect LA and furfural production. However, according to Daorattanachai et al. [34], lignin can promote the isomerization reaction of glucose to fructose, which can then be dehydrated to HMF using a phosphoric acid catalyst. As is well known, lignin removal from lignocellulose biomass is difficult and costly. Therefore, in this study, lignin was not removed from the biomass.

In this study, based on Zulkiple et al. [27], SSA pretreatment produced more sugar than raw POEFBs after enzyme hydrolysis. Ammonia pretreatment is an alkali pretreatment. The alkali agent saponifies the ester bonds in the xylan backbone, resulting in the production of carboxyl groups and the breakdown of lignin-hemicellulose bonds [6].

SEM images before and after the SAA pretreatment are presented in Figure 3. As shown in Figure 3(a), prior to the SSA pretreatment, the POEFB fibres looked hard and stiff, with a flat, smooth surface structure. After the SAA pretreatment, the lignin-carbohydrate bonds in the POEFB fibres broke down and the surface contained pores, indicating that the surface area of the POEFBs had increased (Figure 3(b)). Thus, SAA pretreatment appears able to break the bonds between lignin and carbohydrate. As reported earlier, ammonia will also increase the accessibility of cellulose because it acts as a swelling agent [35].

3.2. Evaluation of the Separate Kinetic Parameters of LA and Furfural Production
3.2.1. Evaluation of the Kinetic Parameters of LA Production

LA optimization results are presented in Figure 4. According to Figure 4(f), the highest yield of LA in the kinetic experiments was 52.1 mol%; this was obtained at a reaction temperature of 170°C, after a 90 min reaction, using an acid concentration of 1 M. Girisuta et al. [22] and Dussan et al. [21] examined the kinetics of LA production from sugarcane bagasse and a Miscanthus giganteus cross, respectively, at a temperature of 150–200°C, with a sulfuric acid catalyst (0.1–0.5 M), and obtained yields of around 60–70% mol. The lower yield found in the present study may be due to the operating temperature, which was quite low (150–170°C), although the acid catalyst concentration was high (0.5–1 M) compared to that used by Girisuta et al. [22] and Dussan et al. [21]. This result shows that temperature affects LA production.

Comparison of activation energy values of LA production at different acid concentrations can be seen in Table 2. According to Table 2, activation energy tended to decrease in higher acid concentrations because sulfuric acid functioned as a catalyst to reduce the activation energy for LA production.

The activation energy for glucose formation was lower (132.37 kJ/mol and 108.48 kJ/mol) than the activation energy for HMF formation (212.40 kJ/mol and 119.49 kJ/mol) in 0.5 and 1 M acid concentration, indicating that the glucose formation reaction was faster than the HMF formation reaction. This proved that higher temperature made the rate of the HMF formation reaction faster.

The HMF formation reaction at 0.5 and 1 M acid concentrations had higher activation energies (212.40 kJ/mol and 119.49 kJ/mol) and acid reaction order (2.00) than the activation energies (188.31 kJ/mol and 62.12 kJ/mol) and acid reaction order (1.66) of humin formation. These indicated that more glucose decomposed into HMF rather than into humins at higher temperatures. Likewise, the HMF formation reaction had the highest activation energy and acid reaction order compared to all other reactions. These indicated that both temperature and acid concentration have greater effects on this reaction than on other reactions.

The effect of acid concentration also had an important role in humin formation, which can be seen from the second highest acid reaction order value after HMF formation. The higher acid concentration triggered the formation of humins. Humin formation is undesirable in the production of LA because this reaction competes with HMF formation when using glucose as the raw material. Therefore, temperature and acid concentration become important factors to prevent humin formation.

The lowest activation energy was recorded for LA formation (42.55 kJ/mol and 56.08 kJ/mol), suggesting that HMF was quickly converted into LA. The HMF concentration during the kinetic experiments was low due to rapid HMF degradation and LA formation.

The comparison of kinetic parameters in LA production (0.5 M acid concentration) with other research is shown in Table 3. The results from this study were compared with Girisuta et al. [22] and Dussan et al. [21], who examined the kinetics of LA production from sugarcane bagasse and a Miscanthus giganteus cross, respectively, at a temperature of 150–200°C with a lower sulfuric acid catalyst (0.1–0.5 M).

The activation energy of glucose formation in the present study (132.37 kJ/mol) was not too different from the studies by Girisuta et al. [22] and Dussan et al. [21], who reported figures of 144.85 kJ/mol and 188.9 kJ/mol, respectively. The differences could be due to the different characteristics of the biomass used.

In the present study, the lowest activation energy for LA formation was the same as that recorded in other research. Girisuta et al. [22] and Dussan et al. [21] reported low activation energies of about 101.63 kJ/mol and 121.3 kJ/mol, respectively. In this study, the activation energy conversion of LA formation was much smaller (42.55 kJ/mol), indicating that LA readily forms at a lower activation energy.

SAA pretreatment was the most likely cause of the higher activation energy in HMF formation (212.40 kJ/mol) as compared to humin formation (188.31 kJ/mol) in this study. In contrast, Girisuta et al. [22] and Dussan et al. [21] reported that the activation energy of HMF formation was lower than the activation energy of humin formation. The ammonia pretreatment likely improved the accessibility of the glucose, thereby resulting in increased production of HMF rather than humins. The results suggested that, at higher temperatures, the reaction pathway will tend to lead HMF formation rather than humins. The humins in this study may be products from cellulose conversions other than HMF.

The acid concentration had a larger effect on LA formation, which is consistent with the data in earlier studies. Girisuta et al. [22] and Dussan et al. [21] reported acid reaction orders of approximately 1.32 and 1.95, respectively, which were relatively similar to the acid reaction order in the present study (1.15). Therefore, higher acid concentration will enhance LA formation.

3.2.2. Evaluation of the Kinetic Parameters of Furfural Production

Furfural optimization results are presented in Figure 5. As shown in Figure 5(e), the largest furfural yield in this study (27.94 mol%) was obtained at a temperature of 170°C, with a 20 min reaction and a 0.5 M acid concentration. The furfural yield in this study was low. This was likely due to the high temperature (170°C), resulting in the formation of degradation products (formic acid and tar), as reported previously by Danon et al. [36]. In addition, the presence of lignin may have affected furfural production. Lamminpaa et al. [37] showed that lignin can increase the pH of a reactant solution, leading to low conversion of xylose into furfural.

In the present study, xylose was detected only at a temperature of 150°C, after a 60 min reaction time in a 0.5 M acid concentration. Xylose was undetectable at higher temperatures, higher acid concentrations, and longer reaction times, in contrast to the findings of Dussan et al. [21] and Chin et al. [20]. Both studies detected xylose at higher reaction temperatures of about 200°C and 180°C, respectively. The kinetic model by Dussan et al. [21], which incorporated xylose, was not in line with the experimental data. Therefore, xylose was not included in the furfural kinetic model.

The SAA pretreatment explains the differences between xylose production in this study as compared to those observed in the kinetic experiments conducted by Dussan et al. [21] and Chin et al. [20]. Hemicellulose is an amorphous polymer, which is more easily degraded than cellulose. In the presence of ammonia pretreatment, the increased accessibility of hemicellulose means it is likely to be converted into a variety of products. Therefore, during acid hydrolysis, xylose located in hemicellulose amorphous fibres will be directly converted into furfural.

In this study, according to Table 4, the activation energy of furfural formation (76.76 kJ/mol) was higher than the activation energy of furfural decomposition (9.89 kJ/mol) in 0.5 M acid concentration. This means that, at higher temperatures, reactions can produce more furfural. The largest furfural yield in this study (27.94 mol%) was obtained at a temperature of 170°C in 0.5 M acid concentration.

In contrast, the higher acid concentration (1 M) made the activation energy of furfural formation (59.22 kJ/mol) lower than the activation energy of furfural decomposition (77.08 kJ/mol). This means that, at higher temperatures, reactions can degrade furfural, as indicated by the fact that the furfural yield was lower at a temperature of 170°C in 1 M acid concentration in the present study.

3.3. Evaluation of Simultaneous Kinetic Parameters of LA and Furfural Production

The results of the evaluation of simultaneous kinetic parameters of LA and furfural are presented in Figure 6, and the comparison of separate and simultaneous LA and furfural kinetic parameters at different acid concentrations can be seen in Table 5, which shows that the activation energy of LA production was reduced because of the sulfuric acid catalyst. The LA activation energy trends are the same in separate and simultaneous LA production, except for EakGLN.

As shown in Table 5, the activation energy of glucose formation (EakGLN) in the simultaneous kinetic model tended to increase when the acid concentration was higher. The higher activation energy of the glucose formation (EakGLN) indicated that there were complex hydrolysis interaction reactions of pretreated POEFBs when LA and furfural kinetic parameters were calculated simultaneously. Glucose was the main monomer derived from the pretreated POEFB cellulose fraction, which was the highest fraction (42%) in pretreated POEFBs; therefore, glucose was the largest monomer composition in pretreated POEFBs. The degradation process from cellulose fraction into glucose was a complex reaction in biomass hydrolysis. In contrast, LA was the product of a series of reaction steps in the acid hydrolysis of biomass. LA formation was a complex reaction produced from C6 sugars (glucose and fructose) in lignocellulose biomass [12]. Therefore, the activation energy of glucose formation (EakGLN) depends on reaction conditions and the interaction between compounds that affect its degradation into other products (HMF and humins).

According to Table 5, the activation energy of furfural formation (171.74 kJ/mol) was higher than the activation energy of furfural decomposition (98.21 kJ/mol) in 0.5 M acid concentration, and the activation energy of furfural formation (88.81 kJ/mol) was lower than the activation energy of furfural decomposition (175.14 kJ/mol) in 1 M acid concentration. This means that, at higher temperatures, reactions can produce more furfural in 0.5 acid concentration, but furfural would be degraded in 1 M acid concentration. Therefore, it was shown that the interaction effect on biomass hydrolysis is very low for furfural formation because hemicellulose is comprised of amorphous fibres that are easily degraded in the biomass hydrolysis process.

The acid reaction order in the simultaneous kinetic model of LA and furfural production showed higher values than that in the separate kinetic model. This suggested that acid had an important role as a catalyst in the hydrolysis reaction of pretreated POEFBs and LA-furfural production.

According to Table 5, the activation energy of LA and furfural acid production from simultaneous kinetic model calculation had the same trend as the separate kinetic model calculation. In separate and simultaneous kinetic models, the activation energy of LA production (EakHMF) increased with higher acid catalyst concentration and the activation energy of furfural production (EakXYN) decreased with higher acid concentration.

Simultaneous kinetic model is better than separate kinetic models for evaluating kinetic parameters in LA and furfural production. Glucan and xylan in separate kinetic models of LA and furfural were combined into one kinetic model (Figure 2) because glucan and xylan are fractions in POEFBs; this is the main reason simultaneous kinetic model was developed in this study. The sum of square error (SSE) of the optimization result between the experimental data and the models in the kinetic parameter estimation was lower in the simultaneous kinetic model (14.32) compared to those in the separate kinetic models (88.23 for LA and 44.04 for furfural). Therefore, the simultaneous kinetic model is better than separate kinetic models.

4. Conclusion

This study evaluated separate and simultaneous kinetic models of LA and furfural production from POEFBs that were pretreated with ammonia. A kinetic experiment was performed using a pressurized reactor at a temperature of 150–170°C, with a sulfuric acid catalyst at concentrations of 0.5 M and 1 M. In the kinetic experiments, the greatest LA yield was 52.1 mol%, which was obtained at a reaction temperature of 170°C, after a 90 min reaction, using an acid concentration of 1 M. The highest furfural yield was 27.94 mol%, which was obtained at a temperature of 170°C, after a 20 min reaction, using an acid concentration of 0.5 M. SAA pretreatment affected activation energy in glucose degradation reactions and favoured direct conversion of hemicellulose to furfural. Based on the evaluation of the kinetic parameters, the simultaneous kinetic model has been shown to have the same trends as the separate kinetic models in LA and furfural production: the activation energy of LA production (EakHMF) increases with higher acid catalyst concentration, and the activation energy of furfural production (EakXYN) decreases with higher acid concentration. Higher reaction temperature and acid concentration will increase LA production. However, higher acid concentration can reduce furfural production because of furfural decomposition. Based on the lower sum of square error (SSE) of the optimization result between the experimental data and the models in the kinetic parameter estimation, simultaneous kinetic model is better to calculate kinetic parameters of LA and furfural production than separate kinetic models.

Nomenclature

A:Preexponential factor (s−1)
C5:C5 sugar concentration (mol/L)
C6:C6 sugar concentration (mol/L)
Cfurfural:Furfural concentration (mol/L)
Cglucan:Glucan concentration (mol/L)
Cglucose:Glucose concentration (mol/L)
CHMF:5-Hydroxylmethylfurfural concentration (mol/L)
CHUM:Humin concentration (mol/L)
CLA:Levulinic acid concentration (mol/L)
CRES:Resinification product concentration (mol/L)
Cxylan:Xylan concentration (mol/L)
Ea:Activation energy (kJ/mol)
FA:Formic acid
FAME:Fatty acid methyl esters
GVL:γ-Valerolactone
H2SO4:Sulfuric acid
HCl:Hydrochloric acid
HMF:5-Hydroxylmethylfurfural
HPLC:High-performance liquid chromatography
k:Reaction rate constant (s−1)
kFUR:Furfural reaction rate constant (s−1)
kGLC1:Glucose reaction rate constant 1 (s−1)
kGLC2:Glucose reaction rate constant 2 (s−1)
kGLN:Glucan reaction rate constant (s−1)
kHMF:5-Hydroxylmethylfurfural reaction rate constant (s−1)
kXLN:Xylan reaction rate constant (s−1)
LA:Levulinic acid
n:Acid reaction order
NREL:National Renewable Energy Laboratory
POEFBs:Palm oil empty fruit bunches
R:Gas constant (kJ/mol·°C)
SAA:Soaking in aqueous ammonia
SEM:Scanning electron microscopy
SSE:Sum of square error
T:Temperature (°C)
YF:Furfural yield (mol%)
YLA:Levulinic acid yield (mol%)
Yn,experiment:Experimental yield for every sampling time (mol/mol)%
Yn,predicted:Model yield predicted for every sampling time (mol/mol)%.

Data Availability

The authors confirm that all data underlying the findings in the paper are fully available without restriction. The calculations and data used to support the findings of this study are included within the paper and supplementary information file.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors acknowledge financial support from Indonesia Estate Crop Fund for Palm Oil (BPDPKS; Research Grant no. Peng-01/DPKS.4/2015); from the Ministry of Research, Technology and Higher Education of the Republic of Indonesia through the World Class Professor (WCP) Program (no. 168.A10/D2/KP/2017); and from the USAID through the SHERA program—Centre for Development of Sustainable Region (CDSR).

Supplementary Materials

LA and furfural reaction rate constant data at optimum slope after MATLAB optimization can be seen in supplementary materials. (Supplementary Materials)