Skip to main content

Delignification of corncob via combined hydrodynamic cavitation and enzymatic pretreatment: process optimization by response surface methodology



Renewable liquid biofuel production will reduce crude oil import of India. To displace the huge quantity of fossil fuels used for energy production, this research was focused on utilization of unexploited low-cost agricultural residues for biofuel production. Corncobs are a byproduct of corn processing industry, and till now it is not utilized for biofuel production, eventhough it has high lignocellulosic concent. In this study, combined hydrodynamic cavitation and enzymatic (HCE) method was evaluated as a pretreatment method of corncob for biofuel production. The most significant process parameters namely (i) enzyme loading (3–10 U g−1), (ii) biomass loading (2.5–5.0%), and (iii) duration (5–60 min) were optimized and their effects on combined HCE pretreatment of corncob was studied through response surface methodology for lignin reduction, hemicellulose reduction and cellulose increase.


The highest lignin reduction (47.4%) was obtained in orifice plate 1 (OP1) under the optimized conditions namely biomass loading at 5%, enzyme loading at 6.5 U g−1 of biomass, and reaction duration of 60 min. The above tested independent variables had a significant effect on lignin reduction. The cavitational yield and energy consumption under the above-mentioned optimized conditions for OP1 was 3.56 × 10−5 g J−1 and 1.35 MJ kg−1, respectively.


It is evident from the study that HCE is an effective technology and requires less energy (1.35 MJ kg−1) than other biomass pretreatment methods.


Maize (Zea mays L.) is an important cereal grown in India for food and nutritional security, with a total production and productivity of 24 million metric tons and 10.2 t ha−1, respectively. It is estimated that about 180 kg of cobs are produced per ton of shelled maize, leading to huge amounts of corncob accumulation, which are left unutilized. Corncobs are a rich source of lignocellulosic material containing primarily cellulose (36%), hemicelluloses (26%), and lignin (17%). Currently, the corncobs are used for production of paper pulp due to its cellulosic contents. However, it can be a better source for biofuel production due to its constituents. It should be noted that corncob is not a food or feed substrate. Production of biofuels and other value-added industrially important fine chemicals from lignocellulosic biomass will not only lower our dependence on fossil fuels and green gas emissions but also would favor the development of sustainable biorefineries. Hence, utilization of corncob is the best resource for deriving useful industrial chemicals and corncob-based biofuel as one of the viable and promising approaches for sustainable utilization of natural resources and ensuring energy security.

Conversion of lignocellulosic biomass (LCB) into bioalcohols (bioethanol or biobutanol) involves three stages, viz, biomass pretreatment, hydrolysis, and fermentation and among these biomass pretreatment is an inevitable bioprocesses and consumes up to 40% of total production cost [1]. Currently, physical or chemical or biological methods or their combinations are employed for biomass pretreatment. The most common biomass deconstruction methods includes physical (grinding/milling [2, 3]), chemical (acid/alkaline [4, 5], ionic liquids [6], subcritical and supercritical water, organosolv [7]), biological (white rot fungi [8] or laccase enzyme [9]), and combinations of these methods (hot water/autohydrolysis [10, 11], steam explosion [11], supercritical CO2 [12] or ammonia fiber explosion [13]). Generally, combined pretreatment methods are preferred due to its better performance/efficiency in delignification compared to individual methods [14]. The biomass pretreatment methods which are used currently are more complex and energy intensive, eventually requires more efforts during scaling up process.

Recently, hydrodynamic cavitation (HC) technology has been successfully employed in diversified fields namely biodiesel production, wastewater treatment, food processing and different bioprocessing operations. Combined HC pretreatment and other chemical catalysts have been used for production of biofuel from different lignocellulosic feedstocks [15,16,17,18]. During HC pretreatment, highly reactive radicals (OH and H) are formed in the working fluid due to the cavitation effect, which can cause deconstruction of lignin structures. Though HC technology has been considered as a viable pretreatment method for processing diverse feedstocks, it is not optimal for targeting specific end products.

The methodology adopted in HC pretreatment studies creates cavitation by continuous circulation of the working fluid (water containing biocatalyst) by a pump through an orifice plate of the reactor. However, the raw biomass has to be kept within or outside of the HC zone [19]. But in the present investigation, both the biomass and biocatalyst were circulated continuously throughout the reactor to enhance the bioconversion efficiency. The advantages of this approach are less energy requirement and operation at ambient conditions.

Furthermore, biocatalysts such as enzyme use phenoxy radicals for the removal of recalcitrant fractions of lignocellulosic biomass. The best-characterized lignin depolymerizing enzymes are multi copper oxidases namely laccases (EC, which utilizes O2 instead of H2O2 to mediate substrate oxidation and also not prone to cofactor degradation. To the best of our knowledge, combination of HC and enzyme (HCE) was not attempted so far. In this context, the study highlights the significance of combination of HC with laccase for pretreatment of corncob biomass, optimization of process variables, and their interaction effects on higher recovery of lignin from corncob.

Results and discussion

Optimization of HCE method

Optimization of HCE pretreatment of delignification was performed using orifice plate 1 (OP1) and orifice plate 2 (OP2) at an optimized inlet pressure of 50 and 100 kPa, respectively. To determine the optimum values of variables for HCE pretreatment, response surface methodology (RSM) was employed with all the input variables set in the range, with output responses set at maximum levels. For process optimization, three factors were evaluated for three different responses such as maximum lignin reduction (%), hemicellulose reduction (%) and cellulose increase (%) using Box–Behnken design in RSM, which includes biomass loading (A), enzyme loading (B), and time (C). These variables were statistically optimized using 33 factorial design.

The optimized conditions for HCE pretreatment with orifice plate 1 (HCE-OP1), are biomass loading at 5%, enzyme loading at 6.5 U g−1 of biomass, and a process time of 60 min. The result showed a reduction in lignin (47.4%) and hemicellulose content (3.2%) and an increase in cellulose content (25.3%). Furthermore, it was observed that the biomass loading and duration have resulted in the optimum response with the higher values. However, the enzyme loading resulted in the mid value. The quadratic model of percentage of lignin reduction had a predicted and adjusted R2 values of 0.99 and 0.98, respectively. The predicted R2 value was in agreement with the adjusted R2 value, since the difference between two values was less than 0.02, indicating that the model developed will be able to give a reasonably good estimate of the response of the system [20].

An analysis of variance (ANOVA) was conducted to check the adequacy and the significance of the quadratic model. According to ANOVA (Table 1), the F value determined the significance of each term at the designed level of confidence [21]. The F value of percentage of lignin reduction (1166.54) implies that the model is significant. The p value was used to check the significance of each variable and simultaneously identify the effect of each factor on the response [22, 23]. According to the p values (p < 0.05), the linear model terms (A, B and C), the interactive model terms (AC), and the quadratic model terms (C2) were all significant at 95% confidence level. Among the interactions studied, biomass loading and duration resulted in probability value (p value) of 0.0198, which was highly significant than other interactions. Furthermore, the variables A, B and C have resulted in significant model terms compared to AB and BC. In this model, only significant terms are considered. The coefficient of variance (CV) specifies the degree of precision for the treatments and expressed as percentage (%). The low values of CV% clearly indicated a very high degree of precision and a good reliability of the experimental values [20]. Percentage of lignin reduction resulted in a lower CV value of 1.92%, which showed that the model has high experimental reliability due to closeness of experimental and predicted values [24].

Table 1 Analysis of variance (ANOVA) for quadratic model for lignin reduction of corncob pretreated by HCE-OP1

Adequate precision measures the signal to noise ratio and compares the range of the predicted values at the design points to the average prediction error. This ratio was greater than 4, which is desirable and indicates adequate model discrimination [25]. In the present case, the adequate precision for a percentage of lignin reduction has resulted in a ratio of 103.04, indicating an adequate signal, implying greater predicted response relative associated error.

The response surface plots were generated with one variable kept at its optimum level and other variables varying within the experimental range. The interactions between the model terms were expressed by a three-dimensional surface. A check for interactions was necessary to determine the significance of the model equation [22, 23]. The response surface plots for a percentage of lignin reduction of the OP1 were shown in Fig. 1. An empirical relationship between the response and the variables was expressed by the following quadratic second-order polynomial equation:

figure 1

3D plots for response lignin reduction of corncob pretreated by HCE –OP1. a Enzyme loading vs. biomass loading; b biomass loading vs. time; c enzyme loading vs. time

$$\begin{aligned} {\text{Lignin reduction}} & = - 1.80307 \, + \, 0.505073 \, \times \, A \, + \, 0.037273 \, \times \, B \, + \, 0.185913 \times \, C \, - \, 0.01686 \\ & \quad \times \, AB \, - \, 0.00357 \, \times \, AC \, + \, 0.000169 \, \times \, BC \, + \, 0.000271 \times \, A^{2} + \, 0.003909 \\ & \quad \times \, B^{2} - \, 0.00106 \, \times \, C^{2} \\ \end{aligned}$$

where A is biomass loading in %, B is enzyme loading in U g−1 of biomass, and C is duration in min.

The same experimental designs and conditions were adopted as per OP1 for optimizing variables for OP2. The optimized conditions obtained in HCE pretreatment with plate 2 (HCE-OP2) were the same as plate 1, namely biomass loading (5%), enzyme loading (6.5 U g−1) and a process time (60 min). However, in this case, lignin reduction was 35.91%, which was significantly lower than that observed with OP1. This might be due to the cavitation bubbles collapse and inward propagation of the shock waves. The geometric focusing of this shock wave at the centre of the cluster creates the enhancement of the cavity cluster collapse (higher pressure) and noise associated with the cavity cluster collapse [26]. This kind of shock wave formation has been confirmed by Wang and Brennen [27] and Reisman et al. [28]. The collapse of micro bubble cavities generates enormous destructive forces, which not only brings about the disintegration of biomass via high turbulence and micro-jets [19, 29], but also due to the dissociation of water molecules, which led to the generation of radicals like HO, HOO, and O2. The linear model of lignin reduction (Table 2) showed a predicted and adjusted R-squared (R2) values of 0.94 and 0.96, respectively, revealing a good agreement between experimental and predicted values besides implying that the mathematical model was very reliable. The F value of 140.32 implies that the model was significant. The compositional analysis of corncob biomass samples after HCE with OP1 and OP2 are furnished in Additional file 1: Tables S1, S2.

Table 2 Analysis of variance (ANOVA) for linear model for lignin reduction of corncob pretreated by HCE –OP2

According to the p values (p < 0.05), the linear model terms (A and C) are significant at 95% confidence level. The CV value of 13.43 indicated that the model had reliability of the experimental values. An adequate precision ratio of 130.45 represents an adequate signal. The response surface plots are shown in Fig. 2. The empirical equation fitting the quadratic model is given below:

$$\text{Lignin reduction} = \text{1}{.818182 \times A } + 0{.318588 \times B} + 0{.598915 \times C}$$

where A is biomass loading in %, B is enzyme loading in U g−1 of biomass, and C is duration in min.

figure 2

3D plots for response lignin reduction of corncob pretreated by HCE-OP2—enzyme loading vs. biomass loading

The results of the present study are compared with different HC biomass pretreatment methods which is shown in the Table 3. The present experiment achieved 47% lignin recovery with a lower optimal temperature of operation. More importantly, the energy consumption for the HCE pretreatment process is very low compared to other reported methods.

Table 3 Comparison of different HC biomass pretreatment methods

Thermal behavior of raw and pretreated biomass

Thermogravimetric analysis of lignocellulosic biomass was used to study the thermal degradation profiles of hemicellulose (250–300 °C), cellulose (300–350 °C), and lignin (300–500 °C) [31, 32]. The appearance of the derivative thermogravimetry (DTG) peak at 200–400 °C was pertaining to hemicellulose and cellulose component [31]. This information enables comparison of the changes in composition of biomass due to pretreatment method [33, 34]. Furthermore, it was clear that the two peaks observed between 200 and 400 °C in the DTG curves was associated with the degradation of hemicellulose and cellulose of the sample. Perhaps the ranges of temperature for decomposition of lignin would have overlapped partially with that of hemicelluloses and cellulose. It is obvious that the temperature associated with two DTG peaks of pretreated corncob samples was higher than raw sample and this could be due to pretreatment effect (Fig. 3 and Table 4).

Fig. 3
figure 3

Combined DTG curves for raw and HCE pretreated corncob biomass samples

Table 4 Details of temperature at peak 1 and peak 2 for corncob biomass obtained from DTG curves

Apart from temperature change, the peaks also differed in position and height indicating changes in proximate composition namely reduction in lignin and hemicellulose and increase in cellulose content, that might be happened during pretreatment process.

FT-IR analysis

Pretreated sample was analyzed by Fourier-transform infrared (FT-IR), and the different wave numbers, functional groups and their corresponding polymer were presented in Table 5. It is evident that raw corncob showed clear peaks for cellulose, hemicellulose and lignin at the corresponding wave numbers.

Table 5 Assignment of functional groups and their corresponding polymers in pretreated corncob

The performance of combined HCE pretreatment indicated that there were no major changes in hemicellulose. However, both the orifice plates behaved differently with regard to lignin removal.

A reduction in peak intensity was observed under OP1 than OP2 compared to the raw biomass, at wave numbers 1320, 1424, 1445, and 1509 cm−1 (Fig. 4). This signifies that lignin removal was higher in OP1, and OP2 was less efficient. Probably this difference between the plates could be attributed to the less number of radicals formed in the cavitational zone by the OP2. Moreover, other functional groups also showed similar changes in both OP1 and OP2.

Fig. 4
figure 4

FTIR spectrum for raw and HCE pretreated corncob samples

Scanning electron microscopy (SEM)

The SEM images of raw sample shows a smooth and undisturbed surface (Fig. 5a). By contrast, after the HCE pretreatment, pores appeared on the surface, and it could be due to removal of lignin. According to Grimaldi et al. [35], the biomass pretreatment removes lignin through destruction of cell wall in two stages namely, loss of cohesion between neighboring cell walls, and degradation inside the cell wall by peeling off and forming holes. The breakdown of structure and pinholes in the surface area of treated sample by cavitation could have enhanced the accessibility for hydrolytic enzymes for increasing the saccharification process [36]. HCE-generated potholes in biomass might have occurred due to collapse of the cavities at the surface of the biomass particles, which could lead to a shear of bimolecular structures, while subsequent laccase treatment led to increase in pothole size (1–6 µm). This indicates effective removal of lignin and deconstruction of biomass in OP1 than OP2 (max. 2 µm) (Fig. 5b, c), which could be attributed to the higher intensity of cavities collapse and turbulence generated in OP1 that resulted in a loosening of the biomass for the effective action of laccase than in plate 2.

Fig. 5
figure 5

SEM images a untreated corncob, b corncob treated by HCE-OP1, c corncob treated by HCE-OP2

X-ray diffraction (XRD)

Raw and pretreated samples were estimated for cellulose crystallinity changes. The value of crystallinity index of pretreated samples (56.0%) for both the orifice plates were higher than the raw sample (52.4%), indicating that HCE pretreatment was more efficient in the removal of amorphous fraction as compared to crystalline fractions. No significant change in the crystallinity of cellulose between treatments and type of plates was observed. This observation testifies that laccase removed only lignin and not changed the cellulose crystallinity. Instead, the HCE pretreatment improved the accessibility of cellulose without modifying the cellulose structure.

Cavitational yield and energy consumption

The optimized conditions for biomass treatment in HCE pretreatment with OP1 and OP2 were biomass loading at 5%, enzyme loading at 6.5 U g−1of biomass, and a process time of 60 min. Energy consumption in the HCE pretreatment process for OP1 (1.35 MJ kg−1) was lower than OP2 (3.24 MJ kg−1) with cavitational yields of 3.56 × 10−5 and 2.70 × 10−5 g J−1, respectively, under optimized conditions for OP1 and OP2 for the enzyme pretreatment.


Biomass preparation

The corncob biomass was procured from the farmer fields, Coimbatore, Tamil Nadu, India and dried at ambient temperature (30 °C) for 12 h to reduce the moisture content from 20 to 14%, then the size was reduced in a sequence through shredder (Chudekar Agro Engg Pvt Ltd., Model: 53 H, India), pin mill (Premium Pulman Pvt Ltd., Model: PPM-12, India) and grinder (Aashapura Enterprises, Model: Stylo 750, India) to ≤ 212 µm, which was sieved via ASTM sieve size No. 70.

Enzyme selection

Generally, two families of ligninolytic enzymes are widely considered for enzymatic delignification, which includes phenol oxidase (laccase) and peroxidases (lignin peroxidase, or LiP and manganese peroxidase, or MnP) [37]. Laccase belongs to the copper oxidase enzyme family, similar to other phenol-oxidizing enzymes and it preferably polymerizes lignin by coupling of the phenoxy radicals produced from oxidation of lignin phenolic groups. Laccase enzyme was selected and used as catalyst for this HCE pretreatment [38].

The enzyme laccase used in the study was from microbial source (Trametes versicolor). The laccase was purchased from Sigma-Aldrich, Bangalore, India and used as such. The laccase activity was determined with the reaction mixture contained appropriately diluted enzyme mixed with 1 mM ABTS in sodium phosphate buffer (50 mM, pH 4.5) at 30 °C for 5 min. using 1 mM ABTS by monitoring changes in absorbance at 420 nm (€ max = 3.6 × 104 M−1 cm−1) spectrophotometrically in a Spectramax 360 (Molecular devices, USA). One unit of enzyme activity refers to the amount of enzyme required to oxidize 1 µM min−1 of the ABTS substrate under standard assay conditions.

System description and operating conditions

The main hurdles associated with commercialization of higher biomass loading pretreatment reactors are complexity of design, reactor geometry and upscale, poor mixing characteristics of reactants, and it is an energy-intensive process. For these reasons, low biomass loading rate (< 20%) is largely preferred in most of the biomass pretreatments [39]. The main criteria to be considered for reactor design is rheological properties of different slurries collected from different unit operations in the fermentation process [40]. Mostly, water with catalyst is used as the working fluid in the recent HC biomass pretreatment studies; whereas, in the present study corncob biomass slurry (buffer + biocatalyst + powdered biomass) was used as the working liquid. The process design of this pretreatment is to supply the biomass slurry via holes in the orifice plate to cavitation zone. The rheological study involving different corncob biomass slurries (2.50, 3.75, 5.00, 6.75, 7.50, 8.75 and 10.00%) showed increased viscosity and yield stress, along with increased biomass solid loading and these slurries exhibit pseudo-plastic or shear-thinning behavior [41]. It is also observed that high biomass loading rate (> 6.75%) of corncob slurry have blocked the holes in orifice plate. Based on this result, low biomass loading rates (2.50, 3.75 and 5.00%) were selected for HCE pretreatment. Hydrodynamic cavitation reactor (HCR) consists of circulation tank (6 L capacity), orifice plate, flanges for orifice plate, centrifugal pump, electrical motor, gate valves for priming and bypass, pressure gauges, and pipe accessories (Fig. 6a). Two pressure gauges (P1 and P2) were fixed on both the downstream and upstream sides of the orifice plate to measure the pressure drop. For same area opening in the orifice plate, higher diameter of the hole is recommended for intensive cavitation applications and vice versa. In this reactor, total hole area openings made in the plate were kept as constant (28.26 mm2) and for constant area openings, two configurations of orifice plates were used, viz., orifice plate 1 (OP1:9 holes and 2 mm Ø) and orifice plate 2 (OP2:4 holes and 3 mm Ø) (Fig. 6c). The suction pipe of the pump was connected at the base of the circulation tank. The main pipeline was divided into three sub-pipelines to serve for three purposes namely priming, bypass line, sub mainline to accommodate flanges and orifice plate, and it was connected to the delivery pipe of the centrifugal pump. To make a closed loop circulation, the working fluid was supplied from the circulation tank to the orifice plate with the help of a centrifugal pump and then sent back to circulation tank. The purpose of the bypass valve was to regulate the flow rate/pressure of the working fluid (biomass slurry) by passing through the orifice plate of the reactor.

Fig. 6
figure 6

a Hydrodynamic cavitation reactor, b orifice plate 1, c orifice plate 2

For HCE biomass pretreatment, the biomass slurry was prepared by mixing appropriate quantities of biomass powder and a biocatalyst (laccase enzyme) in an acetate buffer pH 4.5. For example, 200 g of biomass and 6.5 U g−1 of enzyme were added to 4000 mL of acetate buffer to make biomass slurry for 5% biomass loading and 6.5 U g−1 of enzyme. The prepared slurry was added to circulation tank and made to circulate via selected orifice plate continuously, which enabled the biomass to get exposed to cavitation action. Cavitation involves production and aggressive collapse of micro bubbles to generate more hotspots, having higher temperature and pressure. This is sufficient to make the chemical and physical transformations in the lignocellulosic biomass. During the process, decomposition of water molecule causes free radical formation, such as hydroxyl radicals leading to turbulence action of working fluid in the cavitation zone, which eventually helps in rupturing the lignin barrier. Laccase enzyme can oxidize a variety of phenolic subunits of lignin and other aromatic compounds via radical-catalyzed mechanism to yield oxygen-centred free radical and quinine for subsequent reduction in the polymerization reaction [42,43,44]. Hence, HCE-based biomass treatment helps in the generation of highly reactive free radicals in HCR, which can improve the mass transport process rates and enhance the lignin degradation. Based on these results, three biomass loadings (2.50, 3.75 and 5.00%) were selected for HCE pretreatment.

Compositional analysis

After each run, the biomass slurry was collected from the circulation tank, and filtered via filter cloth to separate the pretreated biomass from supernatant. The supernatant was stored at 4 °C for further analysis. The pretreated biomass was washed twice with distilled water to obtain a neutral pH of 7.0 and the biomass samples were dried in a hot air oven at 45 °C for compositional analysis. National Renewable Energy Laboratory procedure was adopted for analyzing biomass composition of raw and pretreated samples [45]. The percentage of lignin reduction by the pretreatment process was calculated using the following equation.

$$\text{Percentage of lignin reduction = }\frac{{\text{lignin in raw biomass{-}lignin in pretreated biomass}}}{{\text{lignin in the raw biomass}}} \times 100$$

Experimental design

Response surface methodology (RSM) is a statistical and mathematical technique widely used for optimization of process parameters and their interactions on output response(s). The added advantages of optimization via RSM approach are (i) only a minimum number of trials are sufficient to find optimum, (ii) to find a correlation between independent inputs and output responses and (iii) reduction in time, materials and cost because of the less number of trials are needed [46]. The optimization of biomass pretreatment process involves studying the influence of operational parameters and their interactions on lignin removal from raw biomass. For enzymatic biomass pretreatment, process parameters such as biomass loading, enzyme loading and reaction time are crucial factors in deconstruction of the lignin structure as well as lead to release of reducing sugars due to solubilization of biomass.

Response surface methodology (RSM) was employed to determine the optimal conditions for HCE pretreatment to attain the maximum percentage of lignin reduction in the corncob biomass samples. The response was assumed to be influenced by three independent variables, catalyst concentration (A), biomass to liquid ratio (B) and reaction time (C). The range of three selected independent variables used for pretreatment process are: biomass loading of 2.5–5.0%, enzyme loading of 3–10 U g−1 of dried biomass, and reaction time of 5–60 min. Based on the results of the preliminary trials, the above range of levels of the three independent variables were fixed. A total of 17 experimental trials of the three variables were designed by Box–Behnken design via Design-Expert software 10.0 (Stat-Ease, Inc., USA) [47].

Structural composition of biomass

FT-IR analysis

The FT-IR spectra of the tested biomass samples were obtained using an FT-IR (FT-IR 6800 JASCO, Japan). Absorbance spectra were recorded between 4000 and 400 cm−1 wave numbers with a spectral resolution of 4 cm−1 and 64 scans per sample.

SEM analysis

The morphology of raw and pretreated corncob biomass was analyzed by scanning electron microscope (SEM; Quanta 250, FEI, Hillsboro, OR, USA) using an Everhart–Thornley Detector (ETD) detector. The SEM was operated in a vacuum, 10 kV, with a spot size of 4 and a pressure of 20 Pa. The sample images were taken at ×4000 magnification.

XRD analysis

The cellulose crystallinity of the biomass samples was measured using an Ultima IV diffractometer (Rigaku, Japan). Copper Kα radiation, 30.0 kV of voltage, 15 mA of current, and a rate of 2.0°min−1 for a 2θ continuous scan from 4.0° to 70.0° were applied. The crystallinity index was obtained from the ratio of the maximum peak intensity 002 (I002, 2θ = 22.0) and minimal depression (Iam 2θ = 16.5) between peaks 001 and 002 [48, 49].

$${\text{Crystallinity index}} = \frac{{I_{002} - I_{\text{am}} }}{{I_{002} }} \times 100$$

where I002 is the diffraction intensity at 2θ = 22.5°, which represents both the crystalline and amorphous regions, and Iam is the diffraction intensity at 2θ = 18.5°, which represents the amorphous regions.

Thermogravimetric analysis (TGA)

The sample size of the corncob used in the experiment was 10 ± 2 mg in a thermogravimetric analyzer (TA instruments, Model: TGA Q50, USA). The test sample was heated at a heating rate of 10 °C min−1 for the temperature range from 50 to 800 °C and nitrogen gas was purged at a flow rate of 30 mL min−1 to create pyrolysis conditions. TGA curve was plotted using TA software (TA Universal Analysis 2000) for both raw and pretreated samples and the results were compared.

Effect of cavitation on temperature of working liquid and residual enzyme activity

During the experimental trials, the temperature of the working liquid was measured at an interval of 5 min by digital thermometer (Multi-thermometer, India). Similarly, the residual enzyme activity during the reaction period was also determined [50].

Working liquid temperature

The acetate buffer was used as working liquid in the HCE biomass pretreatment and there was a rise in the temperature of the working liquid from 30 to 50 °C in 60 min, which could be attributed to heat energy dissipation by sudden collapses of bubbles and cavities. Since laccase enzyme was used, the temperature of working liquid was maintained at 30 °C by circulating cold water to circulation tank during the experiment.

Residual enzyme activity

Enzyme’s protein conformation changes in temperature, pH, ion concentration, and mechanical stress, and microenvironment of a solution. In this study it was observed that, the laccase enzyme was still active (OP1:23.5%, OP2:32.05%) after HCE treatment, this implies that the residual enzyme can be reused for another batch of pretreatment. For OP1, operated with inlet pressure of 50 kPa, the residual enzyme activities at 5, 10, 20 and 30 min were 34.2, 28.6, 23.5, and 23.5%, respectively. The residual enzyme activity was initially reduced and stabilized after 30 min. While OP2 operating with inlet pressure of 100 kPa showed that the residual enzyme activity was gradually reduced over time (34.2, 33.1, 32.1 and 32.05% in 5, 10, 20 and 30 min, respectively). Typically, the residual enzyme activity obtained for OP1 was less than that of OP2, possibly due to variation in higher inlet pressure, coupled with hole numbers and diameter. The frequency of turbulence decreases with an increase in the diameter of the orifice hole [51]. Laccase, being a green catalyst, is widely used for delignification. Although it contributes for the cost factor, in this study, we have shown that laccase after HCE treatment did not lose its activity much. This indicates that the enzyme can be reused for further pretreatment.

Cavitational yield

The cavitational yield is the result of several design parameters optimized for cavitating reactor [52, 53]. This yield can be enhanced by changing flow conditions and reactor geometry. The orifice plates with higher holes results higher cavitational yield due to increased cavitational effects. Cavitational yield can be defined as number of molecules degraded per unit energy dissipated.

$${\text{Cavitational yield}}\, = \,8.834\, \times \,10^{ - 11} \left( {P_{\text{collapse}} } \right)^{1.1633}$$

Collapse pressure can be predicted by

$$P_{\text{collapse}} \, = \,7527 \, \left( F \right)^{ - 2.55} \, \times \,\left( {P_{I} } \right)^{2.46} \left( {R_{0} } \right)^{ - 0.80} \left( {d_{o} } \right)^{2.37}$$

where R0 is the initial cavity size, mm; PI is the inlet pressure, atm; Do is the diameter of the hole in the orifice plate, mm; F is the percentage free area of holes in the total cross-sectional area of the pipe

The energy required per kg of biomass in an HCE pretreatment process

The energy consumed for treating one kg of biomass by the HCE pretreatment process was calculated by the following equation.

$${\text{Energy consumption, MJ/kg}} = \frac{{{\text{Inlet pressure}} \left( {\frac{N}{{m^{2} }}} \right) \times \;{\text{flowrate }}\left( {m^{3} } \right) \times {\text{reaction time}}\left( s \right)}}{{{\text{Weight of biomass used }}\left( {\text{kg}} \right)}}$$


Based on the results it is clear that HCE pretreatment process have decreased the lignin content and increased the cellulose recovery from corncob. The optimum conditions were 5% of biomass loading, 6.5 U g−1 of enzyme loading and the process time of 60 min to get a maximum reduction in lignin content. The percent reduction of lignin by OP1 and OP2 were 47.44 and 35.91%, respectively. This study provides evidence that HCE pretreatment is more efficient in decreasing the lignin content. However, for scaling up, 5% of biomass loading seems to be low and may be increased to achieve an economic output. We hypothesize that it can be accomplished by altering the orifice plate configuration and operating pressure. The conformity studies (TGA, FT-IR, SEM and XRD) have confirmed the deconstruction of biomass, removal of lignin and increase in cellulose accessability by HCE pretreatment. Overall, HCE pretreatment technology is advantageous than other pretreatment method because it requires lesser energy and easy to scale up, and can be used in conjunction with generation of biofuel and other LC-derived bioproducts from biomass as a pretreatment process.



analysis of variance


crystallinity index


coefficient of variance


derivative thermogravimetry


Fourier-transform infrared


hydrodynamic cavitation


hydrodynamic cavitation and enzymatic pretreatment


hydrodynamic cavitation reactor




lignocellulosic biomass


National Renewable Energy Laboratory

OP1 :

orifice plate 1

OP2 :

orifice plate 2


response surface methodology


scanning electron microscopy


thermogravimetric analysis


X-ray diffraction


  1. Sathitsuksanoh N, George A, Zhang Y, Percival H. New lignocellulose pretreatments using cellulose solvents: a review. J Chem Technol Biotechnol. 2013;88(2):169–80.

    Article  CAS  Google Scholar 

  2. Silva GGD, Couturier M, Berrin J-G, Buléon A, Rouau X. Effects of grinding processes on enzymatic degradation of wheat straw. Bioresour Technol. 2012;103(1):192–200.

    Article  CAS  PubMed  Google Scholar 

  3. Zakaria MR, Norrrahim MNF, Hirata S, Hassan MA. Hydrothermal and wet disk milling pretreatment for high conversion of biosugars from oil palm mesocarp fiber. Bioresour Technol. 2015;181:263–9.

    Article  CAS  PubMed  Google Scholar 

  4. Zhou X, Ma J, Ji Z, Zhang X, Ramaswamy S, Xu F. Sun R-c: dilute acid pretreatment differentially affects the compositional and architectural features of Pinus bungeana Zucc. compression and opposite wood tracheid walls. Ind Crops Prod. 2014;62:196–203.

    Article  CAS  Google Scholar 

  5. Ma L, Cui Y, Cai R, Liu X, Zhang C, Xiao D. Optimization and evaluation of alkaline potassium permanganate pretreatment of corncob. Bioresour Technol. 2015;180:1–6.

    Article  CAS  PubMed  Google Scholar 

  6. Qiu Z, Aita GM, Walker MS. Effect of ionic liquid pretreatment on the chemical composition, structure and enzymatic hydrolysis of energy cane bagasse. Bioresour Technol. 2012;117:251–6.

    Article  CAS  PubMed  Google Scholar 

  7. Amiri H, Karimi K, Zilouei H. Organosolv pretreatment of rice straw for efficient acetone, butanol, and ethanol production. Bioresour Technol. 2014;152:450–6.

    Article  CAS  PubMed  Google Scholar 

  8. Potumarthi R, Baadhe RR, Nayak P, Jetty A. Simultaneous pretreatment and saccharification of rice husk by Phanerochaete chrysosporium for improved production of reducing sugars. Bioresour Technol. 2013;128:113–7.

    Article  CAS  PubMed  Google Scholar 

  9. Oliva-Taravilla A, Moreno AD, Demuez M, Ibarra D, Tomás-Pejó E, González-Fernández C, Ballesteros M. Unraveling the effects of laccase treatment on enzymatic hydrolysis of steam-exploded wheat straw. Bioresour Technol. 2015;175:209–15.

    Article  CAS  PubMed  Google Scholar 

  10. Jiang W, Chang S, Li H, Oleskowicz-Popiel P, Xu J. Liquid hot water pretreatment on different parts of cotton stalk to facilitate ethanol production. Bioresour Technol. 2015;176:175–80.

    Article  CAS  PubMed  Google Scholar 

  11. Huang Y, Wei X, Zhou S, Liu M, Tu Y, Li A, Chen P, Wang Y, Zhang X, Tai H. Steam explosion distinctively enhances biomass enzymatic saccharification of cotton stalks by largely reducing cellulose polymerization degree in G. barbadense and G. hirsutum. Bioresour Technol. 2015;181:224–30.

    Article  CAS  PubMed  Google Scholar 

  12. Srinivasan N, Ju L-K. Statistical optimization of operating conditions for supercritical carbon dioxide-based pretreatment of guayule bagasse. Biomass Bioenergy. 2012;47:451–8.

    Article  CAS  Google Scholar 

  13. Zhao C, Ding W, Chen F, Cheng C, Shao Q. Effects of compositional changes of AFEX-treated and H-AFEX-treated corn stover on enzymatic digestibility. Bioresour Technol. 2014;155:34–40.

    Article  CAS  PubMed  Google Scholar 

  14. Alvira P, Tomás-Pejó E, Ballesteros M, Negro M. Pretreatment technologies for an efficient bioethanol production process based on enzymatic hydrolysis: a review. Bioresour Technol. 2010;101(13):4851–61.

    Article  CAS  PubMed  Google Scholar 

  15. Baxi PB, Pandit AB. Using cavitation for delignification of wood. Bioresour Technol. 2012;110:697–700.

    Article  CAS  PubMed  Google Scholar 

  16. Badve MP, Gogate PR, Pandit AB, Csoka L. Hydrodynamic cavitation as a novel approach for delignification of wheat straw for paper manufacturing. Ultrason Sonochem. 2014;21(1):162–8.

    Article  CAS  PubMed  Google Scholar 

  17. Nakashima K, Ebi Y, Shibasaki-Kitakawa N, Soyama H, Yonemoto T. Hydrodynamic cavitation reactor for efficient pretreatment of lignocellulosic biomass. Ind Eng Chem Res. 2016;55(7):1866–71.

    Article  CAS  Google Scholar 

  18. Hilares RT, dos Santos JC, Ahmed MA, Jeon SH, da Silva SS, Han J-I. Hydrodynamic cavitation-assisted alkaline pretreatment as a new approach for sugarcane bagasse biorefineries. Bioresour Technol. 2016;214:609–14.

    Article  CAS  Google Scholar 

  19. Kim I, Lee I, Jeon SH, Hwang T, Han J-I. Hydrodynamic cavitation as a novel pretreatment approach for bioethanol production from reed. Bioresour Technol. 2015;192:335–9.

    Article  CAS  PubMed  Google Scholar 

  20. Maran JP, Manikandan S, Nivetha CV, Dinesh R. Ultrasound assisted extraction of bioactive compounds from Nephelium lappaceum L. fruit peel using central composite face centered response surface design. Arab J Chem J. 2017;10:1145–57.

    Article  CAS  Google Scholar 

  21. Mazaheri H, Lee KT, Bhatia S, Mohamed AR. Subcritical water liquefaction of oil palm fruit press fiber in the presence of sodium hydroxide: an optimisation study using response surface methodology. Bioresour Technol. 2010;101(23):9335–41.

    Article  CAS  PubMed  Google Scholar 

  22. Liu J-Z, Weng L-P, Zhang Q-L, Xu H, Ji L-N. Optimization of glucose oxidase production by Aspergillus niger in a benchtop bioreactor using response surface methodology. World J Microbiol Biotechnol. 2003;19(3):317–23.

    Article  CAS  Google Scholar 

  23. Guo W-Q, Ren N-Q, Wang X-J, Xiang W-S, Ding J, You Y, Liu B-F. Optimization of culture conditions for hydrogen production by Ethanoligenens harbinense B49 using response surface methodology. Bioresour Technol. 2009;100(3):1192–6.

    Article  CAS  PubMed  Google Scholar 

  24. Timung R, Mohan M, Chilukoti B, Sasmal S, Banerjee T, Goud VV. Optimization of dilute acid and hot water pretreatment of different lignocellulosic biomass: a comparative study. Biomass Bioenergy. 2015;81:9–18.

    Article  CAS  Google Scholar 

  25. Maran JP, Sivakumar V, Sridhar R, Thirugnanasambandham K. Development of model for barrier and optical properties of tapioca starch based edible films. Carbohydr Polym. 2013;92(2):1335–47.

    Article  CAS  PubMed  Google Scholar 

  26. Mørch K. Energy considerations on the collapse of cavity clusters. In: Mechanics and physics of bubbles in liquids. New York: Springer; 1982. p. 313–21.

  27. Wang Y-C, Brennen CE. The noise generated by the collapse of a cloud of cavitation bubbles. American Society of Mechanical Engineers; 1995. p. 17–29.

  28. Reisman G, Wang Y-C, Brennen CE. Observations of shock waves in cloud cavitation. J Fluid Mech. 1998;355:255–83.

    Article  Google Scholar 

  29. Nagula KN, Pandit AB. Process intensification of delignification and enzymatic hydrolysis of delignified cellulosic biomass using various process intensification techniques including cavitation. Bioresour Technol. 2016;213:162–8.

    Article  CAS  PubMed  Google Scholar 

  30. Hilares RT, de Almeida GF, Ahmed MA, Antunes Felipe AF, da Silva SS, Han J-I, dos Santos JC. Hydrodynamic cavitation as an efficient pretreatment method for lignocellulosic biomass: a parametric study. Bioresour Technol. 2017;235:301–8.

    Article  CAS  Google Scholar 

  31. Marcilla A, García AN, Pastor M, León M, Sánchez A, Gómez D. Thermal decomposition of the different particles size fractions of almond shells and olive stones. Thermal behaviour changes due to the milling processes. Thermochim Acta. 2013;564:24–33.

    Article  CAS  Google Scholar 

  32. Song C, Hu H, Zhu S, Wang G, Chen G. Nonisothermal catalytic liquefaction of corn stalk in subcritical and supercritical water. Energy Fuels. 2004;18(1):90–6.

    Article  CAS  Google Scholar 

  33. Caballero J, Marcilla A, Conesa J. Thermogravimetric analysis of olive stones with sulphuric acid treatment. J Anal Appl Pyrolysis. 1997;44(1):75–88.

    Article  CAS  Google Scholar 

  34. Chen WH, Tu YJ, Sheen HK. Impact of dilute acid pretreatment on the structure of bagasse for bioethanol production. Int J Energy Res. 2010;34(3):265–74.

    Article  CAS  Google Scholar 

  35. Grimaldi MP, Marques MP, Laluce C, Cilli EM, Sponchiado SRP. Evaluation of lime and hydrothermal pretreatments for efficient enzymatic hydrolysis of raw sugarcane bagasse. Biotechnol Biofuels. 2015;8(1):205.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Kassim M, Potumarthi R, Tanksale A, Srivatsa S, Bhattacharya S. Enzymatic saccharification of dilute alkaline pre-treated microalgal (Tetraselmis suecica) biomass for biobutanol production. Int J Food Agric Vet. 2014;8:967–72.

    Google Scholar 

  37. Krause DO, Denman SE, Mackie RI, Morrison M, Rae AL, Attwood GT, McSweeney CS. Opportunities to improve fiber degradation in the rumen: microbiology, ecology, and genomics. FEMS Microbiol Rev. 2003;27(5):663–93.

    Article  CAS  PubMed  Google Scholar 

  38. Malherbe S, Cloete TE. Lignocellulose biodegradation: fundamentals and applications. Rev Environ Sci Bio. 2002;1(2):105–14.

    Article  CAS  Google Scholar 

  39. Zhang J, Hou W, Bao J: Reactors for high solid loading pretreatment of lignocellulosic biomass. In: Bioreactor engineering research and industrial applications II. Berlin: Springer; 2015. p. 75–90.

  40. Winkler M. Problems in fermenter design and operation. Chemical Engineering Problems in Biotechnology; 1990. p. 215–350.

  41. Kiruthika T, Ramesh D, Sivakumar U. Rheology of different corncob biomass slurries for hydrodynamic cavitation based biomass pretreatment process. Madras Agric J. 2017;104:279–81.

    Article  Google Scholar 

  42. Gonçalves I, Silva C, Cavaco-Paulo A. Ultrasound enhanced laccase applications. Green Chem. 2015;17(3):1362–74.

    Article  Google Scholar 

  43. Christopher LP, Yao B, Ji Y. Lignin biodegradation with laccase-mediator systems. Front Energy Res. 2014;2:12.

    Article  Google Scholar 

  44. Kunamneni A, Plou FJ, Ballesteros A, Alcalde M. Laccases and their applications: a patent review. Recent Pat Biotechnol. 2008;2(1):10–24.

    Article  CAS  PubMed  Google Scholar 

  45. Sluiter A, Hames B, Ruiz R, Scarlata C, Sluiter J, Templeton D, Crocker D. Determination of structural carbohydrates and lignin in biomass. Lab Anal Proc. 2008;1617:1–16.

    Google Scholar 

  46. Morshedi A, Akbarian M. Application of response surface methodology: design of experiments and optimization: a mini review. J Fund Appl Life Sci. 2014;54:2434–9.

    Google Scholar 

  47. Box GE, Behnken DW. Some new three level designs for the study of quantitative variables. Technometrics. 1960;2(4):455–75.

    Article  Google Scholar 

  48. Segal L, Creely J, Martin A Jr, Conrad C. An empirical method for estimating the degree of crystallinity of native cellulose using the X-ray diffractometer. Text Res J. 1959;29(10):786–94.

    Article  CAS  Google Scholar 

  49. Rodrigues Filho G, de Assunção RM, Vieira JG, Meireles CS, Cerqueira DA, Barud SH, Ribeiro SJ, Messaddeq Y. Characterization of methylcellulose produced from sugar cane bagasse cellulose: crystallinity and thermal properties. Polym Degrad Stab. 2007;92(2):205–10.

    Article  CAS  Google Scholar 

  50. Sun J, Peng R-H, Xiong A-S, Tian Y, Zhao W, Xu H, Liu D-T, Chen J-M, Yao Q-H. Secretory expression and characterization of a soluble laccase from the Ganoderma lucidum strain 7071-9 in Pichia pastoris. Mol Biol Rep. 2012;39(4):3807–14.

    Article  CAS  PubMed  Google Scholar 

  51. Barki M, Ganesha T, Math M. CFD analysis and comparison of fluid flow through a single hole and multi hole orifice plate. Int J Res Adv Technol. 2014;2(8):6–15.

    Google Scholar 

  52. Sivakumar M, Pandit AB. Wastewater treatment: a novel energy efficient hydrodynamic cavitational technique. Ultrason Sonochem. 2002;9(3):123–31.

    Article  CAS  PubMed  Google Scholar 

  53. Bashir TA, Soni AG, Mahulkar AV, Pandit AB. The CFD driven optimisation of a modified venturi for cavitational activity. Can J Chem Eng. 2011;89(6):1366–75.

    Article  CAS  Google Scholar 

Download references

Authors’ contributions

SU has received research grants from DBT, Government of India and conceptualized the idea. KT executed most of the experiments, RD supervised the works. OTP is Co-PI of the project. All authors read and approved the final manuscript.


Thanks to Department of Biotechnology, New Delhi and Russian Ministry of Science and Technology, Russia.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

All data generated or analyzed during this study are included in this published article.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.


We thank the Indo-Russian project on “Development of integrated (biotechnological and nanocatalytic) biorefinery for fuels and platform chemicals production from lignocellulosic biomass (crop/wood residues) No. DBT/NRM/CBE/AGM/2014/R016)” for financial support.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Ramesh Desikan or Sivakumar Uthandi.

Additional file

Additional file 1: Table S1.

Results of compositional analysis of pretreated biomass and lignin reduction (%) for HCE with OP1. Table S2. Results of compositional analysis of pretreated biomass and lignin reduction (%) for HCE with OP2.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Thangavelu, K., Desikan, R., Taran, O.P. et al. Delignification of corncob via combined hydrodynamic cavitation and enzymatic pretreatment: process optimization by response surface methodology. Biotechnol Biofuels 11, 203 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: