Next Article in Journal
Application of Extracellular Vesicles in Gynecologic Cancer Treatment
Next Article in Special Issue
Upgrading of Biobased Glycerol to Glycerol Carbonate as a Tool to Reduce the CO2 Emissions of the Biodiesel Fuel Life Cycle
Previous Article in Journal
Establishment of a Serum-Free Hepatocyte Cryopreservation Process for the Development of an “Off-the-Shelf” Bioartificial Liver System
Previous Article in Special Issue
Zr-Based Metal-Organic Frameworks for Green Biodiesel Synthesis: A Minireview
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Kinetics of Biodiesel Production from Microalgae Using Microbubble Interfacial Technology

1
Microfluidics Research Group, Department of Chemical Engineering, COMSATS University Islamabad, Lahore Campus, Lahore 54000, Pakistan
2
Process and Energy Systems Engineering Center-PRESTIGE, Department of Chemical Engineering, COMSATS University Islamabad, Lahore Campus, Lahore 54000, Pakistan
3
Division of Sustainable Development, College of Science and Engineering, Hamad Bin Khalifa University, Qatar Foundation, Doha 34110, Qatar
4
Department of Chemical and Biological Engineering, The University of Sheffield, Sheffield S10 2TN, UK
*
Author to whom correspondence should be addressed.
Submission received: 20 October 2022 / Revised: 18 November 2022 / Accepted: 25 November 2022 / Published: 29 November 2022
(This article belongs to the Special Issue Acceleration of Biodiesel Production)

Abstract

:
As an alternative to fossil fuels, biodiesel can be a source of clean and environmentally friendly energy source. However, its commercial application is limited by expensive feedstock and the slow nature of the pretreatment step-acid catalysis. The conventional approach to carry out this reaction uses stirred tank reactors. Recently, the lab-scale experiments using microbubble mediated mass transfer technology have demonstrated its potential use at commercial scale. However, all the studies conducted so far have been at a lab scale~100 mL of feedstock. To analyze the feasibility of microbubble technology, a larger pilot scale study is required. In this context, a kinetic study of microbubble technology at an intermediate scale is conducted (3 L of oil). Owing to the target for industrial application of the process, a commercial feedstock (Spirulina), microalgae oil (MO) and a commercial catalyst para-toluene sulfonic acid (PTSA) are used. Experiments to characterize the kinetics space (response surface, RSM) required for up-scaling are designed to develop a robust model. The model is compared with that developed by the gated recurrent unit (GRU) method. The maximum biodiesel conversion of 99.45 ± 1.3% is achieved by using these conditions: the molar ratio of MO to MeOH of 1:23.73 ratio, time of 60 min, and a catalyst loading of 3.3 wt% MO with an MO volume of 3 L. Furthermore, predicted models of RSM and GRU show proper fits to the experimental result. It was found that GRU produced a more accurate and robust model with correlation coefficient R2 = 0.9999 and root-mean-squared error (RSME) = 0.0515 in comparison with RSM model with R2 = 0.9844 and RMSE = 3.0832, respectively. Although RSM and GRU are fully empirical representations, they can be used for reactor up-scaling horizontally with microbubbles if the liquid layer height is held constant while the microbubble injection replicates along the floor of the reactor vessel—maintaining the tessellation pattern of the smaller vessel. This scaling approach maintains the local mixing profile, which is the major uncontrolled variable in conventional stirred tank reactor up-scaling.

Graphical Abstract

1. Introduction

The renewable energy generation is an important factor in reducing harmful effects on the environment caused by the excess use of conventional fuels. Biodiesel derived from inexpensive feedstocks, such as microalgae, can be a suitable alternative for replacing fossil fuel. Biodiesel is generally produced from various resources, such as waste cooking oil, animal fats, energy corps, yeast lipids, and microalgae [1,2]. Biodiesel consists of long chains of carboxylic acids of alkyl ester produced from esterification and transesterification of lipids or oils [3,4].
For sustainable production of biodiesel, feedstock cost is the key parameter affecting overall economics of the process. Since biodiesel production through refined oils feedstocks (soybean oil or sunflower oil) is expensive, low-cost unrefined oil (animal or microalgae oil) can be utilized as a cheaper substitute for biodiesel production [5,6,7]. However, the unrefined feedstock’s generally consists of a large quantity of FFA. The presence of these FFA in feedstock are not suitable for transesterification using base catalysts due to the soap formation [8,9]. Therefore, acid catalyzed esterification is employed to reduce FFA content by converting them into FAMEs [10,11]. The acid-catalysis is significantly slower than transesterification due to the low-miscibility of MeOH with oil reducing overall mass transfer; as a result reaction rate decreases. The slower reaction rate directly affects the overall economics of the process [12,13,14]. However, the development of an economical method for commercial biodiesel production would be a revolutionary milestone for the fuel industry. Large scale biodiesel production could address many global challenges, such as waste management, energy supply, and environmental pollution. The major challenges that hinder commercialization of biodiesel include (1) cost-intensive methods of acid catalyzed esterification of biodiesel feedstock, (2) inability of the present technologies to scale-up, (3) expensive two stage production process, i.e., acid catalysis followed by base catalysis, and (4) expensiveness of various feedstocks [15,16,17].
Recently, microbubble mediated mass transfer technology has proven to increase the reaction rate by injecting one of the reactants in the vapor phase [18,19,20]. Microbubbles provide a large interfacial area, low buoyancy force, and high contact time on the bubble surface, which facilitate the rate of reaction of the system [21,22]. Furthermore, a smaller radius causes the increase of pressure inside the microbubble, as stated by Young–Laplace’s law. Hence, the temperature inside the bubble could be predicted to be higher as compared with the boiling point of the alcohol. This also increases the surface energy of the bubble. All of these factors have yielded an unprecedented higher rate and conversion of the esterification reaction [18,20,23]. For example, Fahed et al. (2019) investigated the effect of microbubble on the esterification reaction by producing ethyl acetate and achieved 79.95% in 35 min compared with the conventional method, which achieved 64% conversion of esterification reaction in 350 min for ethyl acetate production [23]. In another study, Naveed et al. (2019) reported a 97% conversion of oleic acid into biodiesel in 30 min using microbubble technology which is a higher conversion than the conventional method, which achieved 80% conversion in 312 min using H2SO4 as a catalyst [20]. The major focus of these studies was to develop microbubble technology for esterification reaction using single component feedstock. In this context, Fahed et al. (2021) validated the effectiveness of microbubble technology in an unrefined feedstock (chicken fat oil) and showed an overall process conversion of 89.90% in 30 min [18]. To further investigate the effect of microbubble technology, Fahed et al. (2022) investigated the effect of microbubble technology by integrating microbubble technology with heterogeneous catalyst using waste cooking oil to further increase the rate of reaction and achieved an 85% conversion in 20 min [24]. The higher conversion was achieved for both chicken fat oil and waste cooking oil in a shorter period of time, indicating that microbubble technology is an economical method for biodiesel production.
However, all previous studies on microbubble technology mainly focused on lab-scale experiments~100 mL of oil and without significant control over vapor pressure of MeOH. The major focus of this study is to scale-up the microbubble mediated mass transfer technology from lab-scale to semi-pilot scale with up to 3 L volume of oil. The experiments were design using response surface methodology (RSM) and compared with gated recurrent unit (GRU). There are many studies in the literature that use RSM for optimization for biodiesel production [25,26]. However, the current study is the first study that implemented both RSM and GRU using microbubble technology. Furthermore, several technologies have been developed to manipulate the reaction equilibrium to achieve a higher conversion. These reactions are usually slow and limited by reaction kinetics and mass transfer, which are key constraints, such as esterification reactions. However, an entirely different technique has been developed using microbubble mediated mass transfer technology: an inherent liquid–liquid reaction is converted into liquid–vapor reaction, which entirely changes the reaction kinetics. The kinetics of a reaction are entirely based on two hypotheses, i.e., (1) by increasing interfacial area, mass transfer of the process increases to which rate of reaction is also enhanced. (2) Simultaneous removal of a reactant could also increase the reaction kinetics in the forward direction. Moreover, this study also investigated reaction kinetics on the semi-pilot scale the first time to understand the feasibility and compatibility of microbubble technology for further scale-up.
Keeping the semi-pilot scale nature of the study, a commercial feedstock (biomass of Spirulina) and a commercial catalyst p-Toluenesulfonic acid (p-TSA) is chosen to investigate microbubble mediated mass transfer. Oil from Spirulina was derived using solvent extraction method. Experiments were designed using response surface methodology (RSM), a robust model derived from RSM compared with another model developed using gated recurrent unit (GRU). This is the first study that demonstrated the potential scale-up of microbubble technology and compared both RSM and GRU models to increase the commercial feasibility of the process.

2. Material and Methods

2.1. Materials

Spirulina-biomass was purchased from Sentron Asia Company in Lahore, Pakistan. p-Toluenesulfonic acid (p-TSA) was purchased from Sigma Aldrich, St. Louis, MO, USA. MeOH of 99% analytical grade and 99% analytical grade n-Hexane were purchased from DAEJUNG chemicals, Siheung-si, South Korea.

2.2. Lipids Extraction from Spirulina Biomass

Microalgae oil/lipids (MO) were extracted from dry biomass using hexane and MeOH in a ratio of 7:3 by vol %. The biomass of Spirulina and solvent were stirred at 1000 rpm for 6 h at room temperature. Afterward, biomass was separated through filtration (whatman filter paper 42), and oil was recovered by evaporating solvent using vacuum evaporation (Buchi R-210, BUCHI Corporation, New Castle, DE, USA) at 60 °C [27]. The gravimetric method determined the oil yield [27]. Physiochemical properties and lipid composition of derived MO are presented in Table 1.

2.3. Pilot-Scale Experimental Setup for Biodiesel Production

The esterification reaction between MO and MeOH was performed using p-TSA as a catalyst. The schematic diagram of the pilot scale process was shown in Figure 1. In the current pilot-scale process, MeOH vapors were formed using a local fabricated digitally controlled vaporizer purchased from EES Technologies, Lahore, Pakistan. The vaporizer provided MeOH vapor at control/desired vapor flowrate, pressure, and temperature. A customized bubble reactor was fabricated using grade 3 sintered borosilicate glass diffuser. The total volume of the bubble reactor was 3.5 L (radius = 43.18 mm and height = 609 mm), and the working volume was up to 3 L. The experiments were designed using RSM (Design Expert 11). The details of the RSM model are given in the next section. The bubble reactor was filled with different volumes of microalgae MO according to the response surface methodology RSM model. The MeOH vapors were produced in a vaporizer then passed through a borosilicate glass diffuser to form microbubbles. The temperature of the reactor was measured through a thermocouple (Digital thermometer, Jiangsu, China). The sample was collected continuously at a regular interval of 10 min. Once, the reaction had been run for the given time, the samples were filtered (whatman filter paper 42) and washed with deionized water. The samples were dried using a vacuum evaporator (Buchi R-210, BUCHI Corporation, New Castle, DE, USA) and stored for further analysis. All the experiments were performed in triplicate, and their average values with standard deviation was reported.

2.4. Modeling and Experimental Design through Response Surface Methodology

RSM is a statistical and mathematical tool that uses multiple variables to design experiments for optimization [28]. In the current study, BBD was used to design experiments. BBD was used to design a process with more than two factors; this method provides fewer experiments than factorial design. Furthermore, BBD follows a cubical design edge using a midpoint with three levels each (−1, 0, +1) [29].
BBD was used with three factors and five center points in the current study. Three factors used in this study are: A (molar ratio of oil to MeOH = 1:5 to 1:25), B (catalyst dosage = 0 to 5 g wt% of MO), and C (TIme-10 to 90 min). According to this design, a total of 17 runs were conducted to evaluate the current process feasibility. The designed experiments and their response with predicted values of RSM and GRU are shown in Table S-I (Supplementary Information) and the RSM final conversion equation with coded value was given in Equation (1) [30,31]. The goodness-of-fit summary provided by RSM shows that quadric model is best suited for current experimental design Table S-II. The suggested model is best suited for current experimental responses, and the current suggested model is also assessed through analysis of variance (ANOVA) as shown in Table S-III.
C o n v e r s i o n = 88.68 + 3.27 A + 26.13 B + 7.97 C + 3.24 A B + 3.49 A C + 3.98 B C 1.74 A 2 29.76 B 2 5.81 C 2

2.5. Modeling through Gated Recurrent Unit (GRU)

Commonly used artificial neural networks (ANN) involve three main layers, such as (1) input layer, (2) hidden layer (3) output layer and are indicated as x, h, and y. Recurrent neural network, a variant of ANN, was proposed for modeling time series and sequential data [32]. However, the vanishing gradient issue with large sequential data limits the application of RNNs. To solve the vanishing gradient problem, long short-term memory (LSTM), a variant of RNN, was introduced by Hochreiter and Schmidhuber [33]. This new variant of RNN, LSTM, was incorporated with four gates to replace the original hidden state in the memory cell. In 2014, GRU was introduced by adding a gated mechanism to the recurrent neural network [34]. Unlike LSTM, GRU merged the input gate and the forget gate into the update gate, as shown in Figure 2.
Hidden unit activation ( r t ) is processed at a time step using Equation (2):
r t = σ ( W r h t 1 + U r x t )
Here, σ depicts logistic sigmoid function, W r and U r represents weight matrices. After that, by using tanh type layer h t ˜ is calculated using r t :
h t ˜ = tan h ( W ( r t × h t 1 ) + U x t )
Equation (4) is the equation that distinguishes GRU from the LSTM. Here, z t combined the remember gate along with forget gate in LSTM. z t is calculated as follows:
z t = σ ( W z h t 1 + U z x t )
Lastly, the hidden state ( h t ) in GRU is calculated using Equation (5):
h t = ( 1 z t ) ( h t 1 ) + ( z t ) ( h t ˜ )
This study incorporated artificial intelligence modeling for the prediction of response surface, conversion in this case. Gated recurrent unit (GRU) was utilized as a sequence learning deep learning technique that requires suitable value of its hyper-parameters. Suitable architecture of the GRU model was obtained by varying hyper-parameters values (Table 2). To study the comparison between GRU and RSM, different parameters were considered, such as R2, RMSE, MAPE, and mean absolute error MAE.

2.6. Biodiesel Analysis

Biodiesel analysis was performed using gas chromatography (GC) and ASTM method of FFA analysis. Briefly, in GC (Shimadzu GC-2014, Shimadzu Europa, Duisburg, Germany) the system was equipped with a flame ionization detector with column EN14103 (30 m × 0.32 mm id. × 0.25 µm film thickness). Nitrogen was introduced as a carrier with an initial temperature of 523 K and a split ratio of 50:1 [20]. For FFA analysis, the AOCS standard titration method was used [36,37]. To determine the FFA of the solution, the following Equations (6) and (7) were used [27]:
Acid   value ( mgKOH g   biodiesel ) = ( F A F B ) ×   N   × 56.11 W
  Free   fatty   acid   ( FFA )   ( % ) = 1 2 × AV

3. Results and Discussion

3.1. Effect of Different Parameters on Free Fatty Acid Conversion

3.1.1. Effect of Catalyst Loading and Molar Ratio on Free Fatty Acid Conversion

The simultaneous effect of the molar ratio of oil: MeOH and catalyst loading are shown in Figure 3. The graph indicates that increasing catalyst loading from 0 to 5 wt% of MO significantly increases the conversion of the process. Increasing the catalyst loading enhances the protonation of FFA in MO. As the degree of protonation increases, the conversion of FFA and the reaction rate also increases. A further increase in catalyst loading after a certain point conversion of FFA was not increased due to insufficient active sites of MO. Furthermore, Perturbation plot of RSM indicate that most dominate factor in the current study is catalyst loading as show in Figure S-II.
On the other hand, by increasing the molar ratio of oil and MeOH, the FFA conversion also increases due to the contact time of FFA with MeOH vapors increasing. At a lower molar ratio, less volume of MeOH passes through the MO and leaves the system and vice versa. However, increasing the molar ratio after a certain limit, a small amount of MeOH can start to accumulate in the reactor, slightly reducing the reaction rate, as indicated by the results. RSM’s optimized condition was molar ratio: 1:23.73 Oil: (MeOH) and catalyst loading 3.3 wt% of MO.

3.1.2. Effect of Reaction Time and Molar Ratio on Free Fatty Acid Conversion

The effect of both time and molar ratio was investigated and are shown in Figure 4. Time is another important parameter to study the rate of reaction. The results indicate that increasing time conversion and the rate of reaction increase due to an increase in reaction time increases the contact time of oil molecules with MeOH. As a result, conversion of FFA increases. However, a higher molar ratio and less reaction time decrease the reaction rate due to the bubble having less contact time with MeOH, as the flowrate of MeOH vapors is too high. As a result, the result MeOH vapors leave the system unreacted [18,20,23].
Furthermore, an increase in MeOH flowrate also increases the formation of macrobubbles at higher flow rates and tends to produce larger bubbles [38]. Macrobubbles have a large buoyancy-force and less residence time compared to microbubbles. As they rise significantly faster, MeOH does not come into contact with FFA and, as a result, MeOH leaves the system unreacted. Both optimize reaction time and molar ratio can also provide an optimized flow rate to increase FFA conversion. The optimized time and molar ratio according to RSM were molar ratios: 1:23.73 Oil: (MeOH) and Time = 59.79~60 min.

3.1.3. Effect of Reaction Time and Catalyst Loading on Free Fatty Acid Conversion

The simultaneous effect of catalyst loading and reaction time were important parameters for scaling up of biodiesel process. These parameters greatly affect the cost and energy of the biodiesel process. The interactive effect of catalyst loading and reaction time on response surface are shown in Figure 5. FFA conversion tends to increase with the reaction time. However, the effect of reaction time is masked by the effect of catalyst loading. The 3D plot also shows that catalyst loading directly relates to FFA conversion due to an increase in catalyst loading increasing the number of available reaction sides; as a result, higher conversion of FFA was achieved.

3.2. Scale-Up of Microbubble Reactor

The main goal of this work is to investigate the scale-up capacity of microbubble technology at optimized conditions provided by RSM and GRU and verify that provided conditions are suitable for converting FFA into biodiesel. For scale-up of microbubble reactor, different experiments were performed at different volumes of reactor varying from 1 to 3 L (Conversion: 1 L = 99.12%, 2 L = 99.55%, 3 L = 99.45%), as shown in Figure 6. It was observed that an increase in the volume in the reactor has negligible effect on the FFA conversion. An increase in the volume of oil also increased the pressure head of the microbubble reactor. By increasing the volume of oil, microbubbles stay in contact with the oil for a longer period of time. However, it is observed that an equilibrium is achieved at 200 mm and does not change if it is increased any further. This could be explained on the basis of mass transfer occurring across the bubble interface and kinetics of the reaction. At 200 mm, a boundary layer would have formed at the bubble interface stopping the mass transfer of MeOH from inside the bubble to the oil working as a “solid sphere”. Several studies in hot microbubble distillation/stripping show that increasing the liquid layer heat can decrease separation efficiency because the non-equilibrium driving force diminishes [39]. Thus, further upscaling of the reactor should be horizontal by maintaining the microbubble tessellation pattern in the reactor. However, this scaling approach maintains the local mixing profile, which is the major uncontrolled variable in conventional stirred tank reactor upscaling.
Conventionally, biodiesel is produced by mixing both reactants in liquid phase in a batch reactor. The mass transfer is limited by low miscibility of reactants, apart from castor oil/MeOH, which has an OH at C-12, reducing the overall conversion. However, the current process MeOH was injected in the form of vapors (bubbles) and as a result a vapor–liquid system is formed. In vapor–liquid system, conversion of MO is also increased as diffusion of vapor–liquid is higher than the liquid–liquid reaction. Furthermore, microbubbles exhibit less buoyancy force, increasing the residence time of MeOH bubble. FFA is premixed with the catalyst and is already protonated. As the microbubbles of MeOH rise, the reaction between MeOH and protonated MO starts instantaneously. The amount of alcohol available in the interface is in excess as compared to the available MO pushing the reaction in a forward direction. As the bubble rises, the alcohol is transferred into the MO. As the temperature of the reactor is maintained higher than the boiling point of MeOH, it does not condenses in the reactor and leaves as vapor [24]. The vapors of MeOH can be collected, condensed and recycled to improve the process economics. Figure 7 clearly exhibits three different slopes. The highest rate is achieved in the first 10 min of the process. Afterword, the rate slightly slows down and an overall conversion of 97% is achieved in the next 30 min (overall 40 min). Since almost all of the MO has already reacted, only 2% conversion is achieved in the last 20 min. To enhance feasibility and hence the economics, the reaction could be stopped at 10 min and product could be separated using a suitable separation technique, such as distillation or else 99% pure product could be obtained in 60 min. Comparison of current MO study with other conventional acid catalyst based biodiesel production is shown in Table 3.

3.3. Reaction Kinetics of Biodiesel Conversion

Scale-up results already show that by increasing the liquid layer height, the conversion does not change and has been identified; the reaction kinetics for a microbubble mediated esterification system can be upscaled via scaling out the reactor horizontally, maintaining the same aerator pattern along the reactor bottom surface. This is a logical conclusion from the seminal paper of Al-Mashhadani et al. (2015) [41]. In their paper, the geometry of the placement of the internal baffle in an airlift loop reactor is systematically varied, yet the hydrodynamics of the phase distribution is invariant to the baffles position, basically demonstrating that it is insignificant. Only the height above the aerator injection point is shown to matter. This follows logically from the fact that microbubbles that are injected in laminar flow and maintain laminar flow only rise vertically. By tessellating the aerators to provide downcomer regions in between, this configuration replicates the micromixing and bulk mixing profiles. This horizontal tessellation/upscaling approach has been found in all pilot scale studies for microbubble distillation/hot microbubble stripping [39]. Rees-Zimmerman and Chaffin (2021) [42] found that in modeling the hydrodynamics of tall bioreactors with variable bubble size, there is no mass transfer limitation with microbubbles, but the reactor height above a fixed level is immaterial. Hence, studying the kinetics for this critical layer height is sufficient for upscaling horizontally.
The kinetics of this study were investigated under optimized condition of RSM, i.e., the molar ratio of MO to MeOH of 1:23.73 ratio, time of 60 min, and a catalyst loading of 3.3 wt% MO. The results show that the difference between predicted value of RSM (99.10% as shown in Figure S-I) and actual value (99.45%) was less than 0.5% which indicate the validation of the current RSM model. To evaluate the reaction behavior of vapor–liquid, Ha was calculated using Equation (8) to assess the whether a reaction occurs on the bubble’s bulk or surface [43,44].
H a = ( M g / l ) T   k   C b k b l
( M g / l ) 25 = 6.02 × 10 5 ( V l 0.36 μ l 0.61 V g 0.64 )
( M g / l ) T = 4.996 × 10 3 ( M g / l ) 25 exp ( 2539 T )
Mg/l was calculated at 25 and 70 °C by Equations (9) and (10) [45]. For a bubble size less than 2 mm, Equation (11) was used to calculate kbl [46].
k b l = 0.31 ( ( D g / l ) 2 ρ l g μ l ) 1 3
The calculated value of Ha is greater than 1, which indicates that the reaction occurs on the bubble surface, due to which bubble size is a crucial parameter of controlling reaction kinetics. The order of reaction was calculated using E by Equation (12) [43].
E = H a ( 1 H a 1 2 E i )  
Ei was determined by using Equation (13) [43]
E i = 1 + ( M g / l ) T ( C b H b   P g )  
The value of Ha and E show that the current reaction is pseudo first order due to the value of both of them I almost equal. The rate of reaction was calculated by using Equation (14).
    r A = 1 1 k g σ + H a ( D g / l ) T k C b P g  
The values used to calculate for kg  σ , kbl, and H are 5.32 × 10−3 kmol s−1 m−3, 1.24 × 10−4 ms−1 and 43.05 kmol s−1.m−3 Pa, respectively. The final reaction rate was determine using Equation (15);
r A = ( 1.32 × 10 5 ) ( P A × 101 , 325 ) ( C b )
To determine the current activation energy (EA) in the scale-up reactor, the Arrhenius equation was used. The experiments were conducted by varying the reaction temperature (70–90 °C) [43]. The Arrhenius equation was used to develop the relations with rate constants to develop the equation to determine EA using Equation (16) [47]. An Arrhenius plot between ln (k) and 1/T is shown in Figure 8.
  l n   k = E A R T + l n   A °
The EA of the esterification reaction in the current scale-up microbubble reactor is calculated to be 10.01 ± 0.3 kJ mol−1. The current EA was significantly less than conventional processes, as shown in Table 4. This low EA indicates that reactions occur on the surface of bubbles [20]. It also shows that less energy is needed for MO to pass the barrier to form a product. In addition, the latent heat of MeOH was also freely available as free energy, making the reaction nature more exergonic, enhancing the rate of reaction and reducing the EA. The current reaction kinetics and reduced EA successfully show the implementation of current scale-up reactor to industrial level.

3.4. Gated Recurrent Unit and Response Surface Methodology Comparison

RSM model has reported 0.9844 correlation coefficient R2 that shows good fitting efficiency, as shown in Figure 9. However, GRU has shown superior efficiency (Figure 10) and reported 0.9999 R2. The values of mean absolute percentage error (MAPE), root mean square error (RMSE) and mean absolute error (MAE) for RSM model and GRU model have been reported in Table 5.
Table 5 confirms the superiority of the GRU model over RSM as its predicted values are in good agreement with actual values, as shown in Figure 10. It can also be confirmed that GRU reported performance criterions used, such as MAPE, RMSE, and MAE, lesser than that of the RSM model. From Figure 9, deviating RSM model prediction values from actual values are clearly visible, causing lower R2 (0.9844) compared to that of GRU (0.9999), as well as a reported larger MAPE, RMSE and MAE compared to that of the GRU model, as show in Figure 10. The MAE for the RSM model is 2.6847, which is approximately 60 times higher compared to that of the GRU model (0.045). Furthermore, MAPE for the GRU model (0.00083) confirms its superiority over the RSM model (0.0465). Furthermore, in terms of RMSE, GRU reported an RMSE of 0.0515 and RSM reported 3.0832.
Furthermore, a 45-degree perfect line was added to the plot to better understand the prediction accuracy of the models. The prediction performance increases as the scatter points approach toward a 45-degree perfect line and hence decreases the error. The plot depicted that the GRU model prediction of conversion% were following the 45-degree perfect line more strictly than that of the RSM model and hence reported lesser RMSE, MAPE, and MAE. Furthermore, a higher R2 value for the GRU model (0.9999) depicts outperforming performance compared to the R2 value for the RSM model (0.9844). This scatter plot comparison also pointed out that the GRU model has performed much better than the RSM model to predict conversion.

4. Conclusions

The current study successfully developed a robust model which shows a high feasibility of microalgae-based biodiesel on a semi-pilot scale. Both models were in line with experimental observation. In addition, the comparison of both RSM and GRU showed that GRU was more accurate than RSM to predict the conversion. Furthermore, predicted models of RSM and GRU show a proper fit to the experimental result. It was found that GRU produced a more accurate and robust model with correlation coefficient R2 = 0.9999 and root-mean-squared error (RSME) = 0.0515 in comparison with the RSM model with R2 = 0.9844 and RMSE = 3.0832, respectively. Furthermore, the kinetics of the pilot scale microbubble reactor revealed that more than 99.45 ± 1.3% conversion of FFA into biodiesel was achieved in 60 min, and the current reaction follows pseudo first order kinetics with respect to MO. Additionally, a lower EA of 10.01 ± 0.3 kJ mol−1 indicates that less energy was required for reactant to jump the barrier to form product. Although RSM and GRU are fully empirical representations, they can be used for reactor upscaling horizontally with microbubbles. Horizontal out-scaling maintains the liquid layer height, while the microbubble injection replicates along the floor of the reactor vessel—keeping the tessellation pattern of the smaller vessel with the same bubble flux per unit area. This scaling approach maintains the local mixing profile, which is the major uncontrolled variable in conventional stirred tank reactor upscaling. This study should prove to be a milestone in future studies for further scale-up of microbubble mass transfer technology. Further research needs to be carried out in terms of reactor development and finding new material for reactor formation, keeping the process cost to a minimum. Additionally, life cycle analysis should be carried out to calculate the environmental impact of the microbubble technology. Life cycle analysis will provide overall insight toward sustainability of the developed approach to analyze manufacturing, ecological effect, and energy expenditures of current technology.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/bioengineering9120739/s1, Table S-I: The designed experiments and their response with predicted values of RSM and GRU; Table S-II: Goodness-of-Fit summary generated through RSM; Table S-III: Statistical ANOVA analysis of current RSM model; Figure S-I: Predicted values from RSM; Figure S-II: Perturbation plot of experimental parameters.

Author Contributions

Conceptualization, F.R. and N.R.; methodology, F.J.; validation, F.R., A.H. and W.B.Z.; formal analysis, F.A.; data curation, M.W.S.-u.-A.; writing—original draft preparation, F.J.; writing—review and editing, F.J.; supervision, F.R.; funding acquisition, F.R. and W.B.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Higher Education Commission Pakistan (HEC), grant number 20-7924/17, 20-7924/17 and Engineering and Physical Sciences Research Council (EPSRC) grant number EP/S031421/1 and EP/N011511/1.The APC was funded by Engineering and Physical Sciences Research Council (EPSRC) grant number EP/S031421/1 and EP/N011511/1.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to confidentiality of work.

Acknowledgments

FJ and FR want to acknowledge National Research Program for Universities (NRPU) provided by higher education commission project number 20-7924/17 and project number 20-3982/14. In addition, WZ would like to acknowledge EPSRC grant numbers EP/S031421/1 and EP/N011511/1.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

WordsAbbreviations
Activation energy EA
Absolute percentage errorMAPE
Acid valueAV
Box–Behnken designBBD
Concentration of MOCb
Enhancement factor E
Free Fatty AcidFFA
Fatty Acid Methyl Esters FAMEs
Gas constantR
Gated recurrent unitGRU
Gas diffusion coefficient Mg/l
Henry constantH
Hatta numberHa
Interfacial area σ
Infinite enhancement factor Ei
Liquid film coefficient kbl
Mass of biodiesel (g)W
Microalgae oilMO
MethanolMeOH
Molar volume of MOvl
Molar volume of MeOHvg
Mean absolute errorMAE
Normality of KOHN
Partial pressure of MeOHpg
Pressure of MeOH (bar) PA
Pre-exponential factor
Rate constantK
Rate of reactionra
Root mean square errorRMSE
Response surface methodologyRSM
Volume of KOH used for titration (mL)FA
Volume of KOH used for blank titration (mL)FB

References

  1. Malekghasemi, S.; Kariminia, H.-R.; Plechkova, N.K.; Ward, V.C. Direct transesterification of wet microalgae to biodiesel using phosphonium carboxylate ionic liquid catalysts. Biomass Bioenergy 2021, 150, 106126. [Google Scholar] [CrossRef]
  2. Javed, F.; Rehman, F.; Khan, A.U.; Fazal, T.; Hafeez, A.; Rashid, N. Real textile industrial wastewater treatment and biodiesel production using microalgae. Biomass Bioenergy 2022, 165, 106559. [Google Scholar] [CrossRef]
  3. Sitepu, I.R.; Garay, L.A.; Sestric, R.; Levin, D.; Block, D.E.; German, J.B.; Boundy-Mills, K.L. Oleaginous yeasts for biodiesel: Current and future trends in biology and production. Biotechnol. Adv. 2014, 32, 1336–1360. [Google Scholar] [CrossRef] [Green Version]
  4. Soccol, C.R.; Neto, C.J.D.; Soccol, V.T.; Sydney, E.B.; da Costa, E.S.F.; Medeiros, A.B.P.; de Souza Vandenberghe, L.P. Pilot scale biodiesel production from microbial oil of Rhodosporidium toruloides DEBB 5533 using sugarcane juice: Performance in diesel engine and preliminary economic study. Bioresour. Technol. 2017, 223, 259–268. [Google Scholar] [CrossRef] [PubMed]
  5. Banković-Ilić, I.B.; Stojković, I.J.; Stamenković, O.S.; Veljkovic, V.B.; Hung, Y.-T. Waste animal fats as feedstocks for biodiesel production. Renew. Sustain. Energy Rev. 2014, 32, 238–254. [Google Scholar] [CrossRef] [Green Version]
  6. Nath, B.; Das, B.; Kalita, P.; Basumatary, S. Waste to value addition: Utilization of waste Brassica nigra plant derived novel green heterogeneous base catalyst for effective synthesis of biodiesel. J. Clean. Prod. 2019, 239, 118112. [Google Scholar] [CrossRef]
  7. Faisal, A.; Javed, F.; Hassan, M.; Gorji, M.; Akram, S.; Rashid, N.; Rehman, F. Experimental and mathematical nonlinear rheological characterization of chicken fat oil-a sustainable feedstock for biodiesel. Biomass Convers. Biorefinery 2021, 1–8. [Google Scholar] [CrossRef]
  8. Tang, Z.-E.; Lim, S.; Pang, Y.-L.; Shuit, S.-H.; Ong, H.-C. Utilisation of biomass wastes based activated carbon supported heterogeneous acid catalyst for biodiesel production. Renew. Energy 2020, 158, 91–102. [Google Scholar] [CrossRef]
  9. Singh, N.K.; Singh, Y.; Sharma, A. Optimization of biodiesel synthesis from Jojoba oil via supercritical methanol: A response surface methodology approach coupled with genetic algorithm. Biomass Bioenergy 2022, 156, 106332. [Google Scholar] [CrossRef]
  10. Meher, L.C.; Sagar, D.V.; Naik, S. Technical aspects of biodiesel production by transesterification—A review. Renew. Sustain. Energy Rev. 2006, 10, 248–268. [Google Scholar] [CrossRef]
  11. Zimmerman, W.B.; Kokoo, R. Esterification for biodiesel production with a phantom catalyst: Bubble mediated reactive distillation. Appl. Energy 2018, 221, 28–40. [Google Scholar] [CrossRef]
  12. Tran, D.-T.; Chang, J.-S.; Lee, D.-J. Recent insights into continuous-flow biodiesel production via catalytic and non-catalytic transesterification processes. Appl. Energy 2017, 185, 376–409. [Google Scholar] [CrossRef]
  13. Talebian-Kiakalaieh, A.; Amin, N.A.S.; Mazaheri, H. A review on novel processes of biodiesel production from waste cooking oil. Appl. Energy 2013, 104, 683–710. [Google Scholar] [CrossRef]
  14. Li, X.; Wang, Q.; Wu, Y.; Chen, J.; Li, S.; Ye, Y.; Wang, D.; Zheng, Z. Optimization of key parameters using RSM for improving the production of the green biodiesel from FAME by hydrotreatment over Pt/SAPO-11. Biomass Bioenergy 2022, 158, 106379. [Google Scholar] [CrossRef]
  15. Wang, Y.-T.; Fang, Z.; Yang, X.-X. Biodiesel production from high acid value oils with a highly active and stable bifunctional magnetic acid. Appl. Energy 2017, 204, 702–714. [Google Scholar] [CrossRef]
  16. Vardon, D.R.; Moser, B.R.; Zheng, W.; Witkin, K.; Evangelista, R.L.; Strathmann, T.J.; Rajagopalan, K.; Sharma, B.K. Complete utilization of spent coffee grounds to produce biodiesel, bio-oil, and biochar. ACS Sustain. Chem. Eng. 2013, 1, 1286–1294. [Google Scholar] [CrossRef]
  17. Narkhede, N.; Singh, S.; Patel, A. Recent progress on supported polyoxometalates for biodiesel synthesis via esterification and transesterification. Green Chem. 2015, 17, 89–107. [Google Scholar] [CrossRef]
  18. Javed, F.; Shamair, Z.; Hafeez, A.; Fazal, T.; Aslam, R.; Akram, S.; Rashid, N.; Zimmerman, W.B.; Rehman, F. Conversion of Poultry-Fat Waste to a Sustainable Feedstock for Biodiesel Production via Microbubble Injection of Reagent Vapor. J. Clean. Prod. 2021, 311, 127525. [Google Scholar] [CrossRef]
  19. Khan, Z.; Javed, F.; Shamair, Z.; Hafeez, A.; Fazal, T.; Aslam, A.; Zimmerman, W.B.; Rehman, F. Current developments in esterification reaction: A review on process and parameters. J. Ind. Eng. Chem. 2021, 103, 80–101. [Google Scholar] [CrossRef]
  20. Ahmad, N.; Javed, F.; Awan, J.A.; Ali, S.; Fazal, T.; Hafeez, A.; Aslam, R.; Rashid, N.; Rehman, M.S.U.; Zimmerman, W.B. Biodiesel production intensification through microbubble mediated esterification. Fuel 2019, 253, 25–31. [Google Scholar] [CrossRef]
  21. Zimmerman, W.B.; Tesar, V.; Butler, S.; Bandulasena, H.C. Microbubble generation. Recent Pat. Eng. 2008, 2, 1–8. [Google Scholar] [CrossRef] [Green Version]
  22. Zimmerman, W.B.; Al-Mashhadani, M.K.; Bandulasena, H.H. Evaporation dynamics of microbubbles. Chem. Eng. Sci. 2013, 101, 865–877. [Google Scholar] [CrossRef] [Green Version]
  23. Javed, F.; Shamair, Z.; Ali, S.; Ahmad, N.; Hafeez, A.; Fazal, T.; Rehman, M.S.U.; Zimmerman, W.B.; Rehman, F. “Pushing and pulling” the equilibrium through bubble mediated reactive separation for ethyl acetate production. React. Chem. Eng. 2019, 4, 705–714. [Google Scholar] [CrossRef]
  24. Javed, F.; Rizwan, M.; Asif, M.; Ali, S.; Aslam, R.; Akram, M.S.; Zimmerman, W.B.; Rehman, F. Intensification of Biodiesel Processing from Waste Cooking Oil, Exploiting Cooperative Microbubble and Bifunctional Metallic Heterogeneous Catalysis. Bioengineering 2022, 9, 533. [Google Scholar] [CrossRef] [PubMed]
  25. Hariram, V.; Bose, A.; Seralathan, S. Dataset on optimized biodiesel production from seeds of Vitis vinifera using ANN, RSM and ANFIS. Data Brief 2019, 25, 104298. [Google Scholar] [CrossRef] [PubMed]
  26. Manojkumar, N.; Muthukumaran, C.; Sharmila, G. A comprehensive review on the application of response surface methodology for optimization of biodiesel production using different oil sources. J. King Saud Univ.-Eng. Sci. 2020, 34, 198–208. [Google Scholar] [CrossRef]
  27. Park, J.-Y.; Nam, B.; Choi, S.-A.; Oh, Y.-K.; Lee, J.-S. Effects of anionic surfactant on extraction of free fatty acid from Chlorella vulgaris. Bioresour. Technol. 2014, 166, 620–624. [Google Scholar] [CrossRef] [PubMed]
  28. Behera, S.K.; Meena, H.; Chakraborty, S.; Meikap, B. Application of response surface methodology (RSM) for optimization of leaching parameters for ash reduction from low-grade coal. Int. J. Min. Sci. Technol. 2018, 28, 621–629. [Google Scholar] [CrossRef]
  29. Yolmeh, M.; Jafari, S.M. Applications of response surface methodology in the food industry processes. Food Bioprocess Technol. 2017, 10, 413–433. [Google Scholar] [CrossRef]
  30. Jiao, W.; Yu, L.; Feng, Z.; Guo, L.; Wang, Y.; Liu, Y. Optimization of nitrobenzene wastewater treatment with O3/H2O2 in a rotating packed bed using response surface methodology. Desalin. Water Treat. 2016, 57, 19996–20004. [Google Scholar] [CrossRef]
  31. Ruhul, A.; Kalam, M.; Masjuki, H.; Fattah, I.R.; Reham, S.; Rashed, M. State of the art of biodiesel production processes: A review of the heterogeneous catalyst. RSC Adv. 2015, 5, 101023–101044. [Google Scholar] [CrossRef]
  32. Cai, X.; Zhang, N.; Venayagamoorthy, G.K.; Wunsch, D.C. Time series prediction with recurrent neural networks using a hybrid PSO-EA algorithm. In Proceedings of the 2004 IEEE International Joint Conference on Neural Networks (IEEE Cat. No.04CH37541), Budapest, Hungary, 25–29 July 2004; Volume 1642, pp. 1647–1652. [Google Scholar]
  33. Hochreiter, S.; Schmidhuber, J. Long Short-Term Memory. Neural Comput. 1997, 9, 1735–1780. [Google Scholar] [CrossRef]
  34. Cho, K.; Merrienboer, B.v.; Gülçehre, Ç.; Bahdanau, D.; Bougares, F.; Schwenk, H.; Bengio, Y. Learning Phrase Representations Using RNN Encoder Decoder for Statistical Machine Translation. In Proceedings of the EMNLP, Doha, Qatar, 25–29 October 2014. [Google Scholar]
  35. Saif-ul-Allah, M.W.; Qyyum, M.A.; Ul-Haq, N.; Salman, C.A.; Ahmed, F. Gated Recurrent Unit Coupled with Projection to Model Plane Imputation for the PM2.5 Prediction for Guangzhou City, China. Front. Environ. Sci. 2022, 9, 753. [Google Scholar] [CrossRef]
  36. Marchetti, J.; Errazu, A. Comparison of different heterogeneous catalysts and different alcohols for the esterification reaction of oleic acid. Fuel 2008, 87, 3477–3480. [Google Scholar] [CrossRef]
  37. Yin, P.; Chen, L.; Wang, Z.; Qu, R.; Liu, X.; Xu, Q.; Ren, S. Biodiesel production from esterification of oleic acid over aminophosphonic acid resin D418. Fuel 2012, 102, 499–505. [Google Scholar] [CrossRef]
  38. Rehman, F.; Medley, G.J.; Bandulasena, H.; Zimmerman, W.B. Fluidic oscillator-mediated microbubble generation to provide cost effective mass transfer and mixing efficiency to the wastewater treatment plants. Environ. Res. 2015, 137, 32–39. [Google Scholar] [CrossRef] [PubMed]
  39. Desai, P.D.; Turley, M.; Robinson, R.; Zimmerman, W.B. Hot microbubble injection in thin liquid film layers for ammonia separation from ammonia rich-wastewater. Chem. Eng. Process.-Process Intensif. 2021, 180, 108693. [Google Scholar] [CrossRef]
  40. Tashtoush, G.M.; Al-Widyan, M.I.; Al-Jarrah, M.M. Experimental study on evaluation and optimization of conversion of waste animal fat into biodiesel. Energy Convers. Manag. 2004, 45, 2697–2711. [Google Scholar] [CrossRef]
  41. Al-mashhadani, M.K.; Wilkinson, S.J.; Zimmerman, W.B. Laboratory preparation of simulated sludge for anaerobic digestion experimentation. J. Eng. 2015, 21, 131–145. [Google Scholar]
  42. Rees-Zimmerman, C.R.; Chaffin, S.T. Modelling the effect of bioreactor height on stripping fermentation products from the engineered bacterium Geobacillus thermoglucosidasius. Biochem. Eng. J. 2021, 176, 108195. [Google Scholar] [CrossRef]
  43. Levenspiel, O. Chemical Reaction Engineering, 3rd ed.; John Wiley & Sons: Hoboken, NJ, USA, 1998. [Google Scholar]
  44. Kierzkowska-Pawlak, H. Determination of kinetics in gas-liquid reaction systems. An overview. Ecol. Chem. Eng. S 2012, 19, 175–196. [Google Scholar] [CrossRef]
  45. Díaz, M.; Vega, A.; Coca, J. Correlation for the estimation of gas-liquid diffusivity. Chem. Eng. Commun. 1987, 52, 271–281. [Google Scholar] [CrossRef]
  46. Kimweri, H.T.H. Enhancement of Gas-Liquid Mass Transfer in Hydrometallurgical Leaching Systems. Ph.D. Thesis, University of British Columbia, Vancouver, BC, Canada, 2001. [Google Scholar]
  47. Roman, F.F.; Ribeiro, A.E.; Queiroz, A.; Lenzi, G.G.; Chaves, E.S.; Brito, P. Optimization and kinetic study of biodiesel production through esterification of oleic acid applying ionic liquids as catalysts. Fuel 2019, 239, 1231–1239. [Google Scholar] [CrossRef] [Green Version]
  48. Jain, S.; Sharma, M. Kinetics of acid base catalyzed transesterification of Jatropha curcas oil. Bioresour. Technol. 2010, 101, 7701–7706. [Google Scholar] [CrossRef]
  49. Song, E.-S.; Lim, J.-W.; Lee, H.-S.; Lee, Y.-W. Transesterification of RBD palm oil using supercritical methanol. J. Supercrit. Fluids 2008, 44, 356–363. [Google Scholar] [CrossRef]
  50. Ahmad, A.; Yasin, N.M.; Derek, C.; Lim, J. Kinetic studies and thermodynamics of oil extraction and transesterification of Chlorella sp. for biodiesel production. Environ. Technol. 2014, 35, 891–897. [Google Scholar] [CrossRef]
  51. Nautiyal, P.; Subramanian, K.; Dastidar, M. Kinetic and thermodynamic studies on biodiesel production from Spirulina platensis algae biomass using single stage extraction–transesterification process. Fuel 2014, 135, 228–234. [Google Scholar] [CrossRef]
Figure 1. Process flow diagram of scale-up microbubble reactor.
Figure 1. Process flow diagram of scale-up microbubble reactor.
Bioengineering 09 00739 g001
Figure 2. Comparison parameter between (a) LSTM and (b) GRU framework [35].
Figure 2. Comparison parameter between (a) LSTM and (b) GRU framework [35].
Bioengineering 09 00739 g002
Figure 3. Effect of molar ratio (oil: MeOH) and catalyst loading on biodiesel production.
Figure 3. Effect of molar ratio (oil: MeOH) and catalyst loading on biodiesel production.
Bioengineering 09 00739 g003
Figure 4. Interactive effect of molar ratio and reaction time.
Figure 4. Interactive effect of molar ratio and reaction time.
Bioengineering 09 00739 g004
Figure 5. Interactive effect of reaction time and catalyst loading.
Figure 5. Interactive effect of reaction time and catalyst loading.
Bioengineering 09 00739 g005
Figure 6. Effect on reaction conversion by increasing microalgae oil volume.
Figure 6. Effect on reaction conversion by increasing microalgae oil volume.
Bioengineering 09 00739 g006
Figure 7. Effect microalgae oil conversion with respect to reaction time.
Figure 7. Effect microalgae oil conversion with respect to reaction time.
Bioengineering 09 00739 g007
Figure 8. Arrhenius plot for microalgae oil based biodiesel reaction.
Figure 8. Arrhenius plot for microalgae oil based biodiesel reaction.
Bioengineering 09 00739 g008
Figure 9. Response surface methodology model performance prediction (actual versus predicted).
Figure 9. Response surface methodology model performance prediction (actual versus predicted).
Bioengineering 09 00739 g009
Figure 10. Gated recurrent unit model prediction.
Figure 10. Gated recurrent unit model prediction.
Bioengineering 09 00739 g010
Table 1. Physicochemical properties of microalgae oil.
Table 1. Physicochemical properties of microalgae oil.
ParametersUnitsValue
FFA content%32.5 ± 2
Density (25 °C)Kg m3920 ± 5
Kinematic viscosity (40 °C)mm2 s−130.06 ± 3
MO composition
Mystic acid (C14:0)%1.90 ± 0.5
Palmitic acid (C16:0)%35.67 ± 3
Palmitoleic acid (C16:1)%6.11 ± 2
Linoleic acid (C18:2)%48.55 ± 2
Linolenic acid (C18:3)%2.17 ± 0.5
Stearic acid (C18:0)%5.60 ± 2
Table 2. Hyperparameters for gated recurrent unit model.
Table 2. Hyperparameters for gated recurrent unit model.
HyperparametersBoundsSet Values
Number of hidden unitsPositive integers50
Gradient threshold0–10.1
Initial learning rate0–10.01
Learn rate drop factor0–10.2
Learn rate drop periodPositive integers100
Training EpochsPositive integers150
Table 3. The comparison of current up-scaling study with other conventional catalyst based biodiesel production.
Table 3. The comparison of current up-scaling study with other conventional catalyst based biodiesel production.
FeedstockCatalystReaction Time
(min)
Conversion
(%)
Reference
Conventional methodH2SO412078[40]
Microbubble TechnologyH2SO43098[20]
Microbubble Technologyp-TSA3097[20]
Microbubble Technologyp-TSA3089.90[18]
Microbubble TechnologySr/ZrO22085[24]
Microbubble Technology (semi pilot-scale)p-TSA6099.45 ± 1.3This study
Table 4. Comparison of activation energy with different biodiesel feedstocks.
Table 4. Comparison of activation energy with different biodiesel feedstocks.
FeedstockMethodScale of ExperimentsCatalystEA
(kJ mol−1)
Reference
JatrophaConventional methodLab-scale1% H2SO4 and 1% NaOH87.808[48]
MicroalgaeSupercritical methodLab-scaleNo catalyst105[49]
ChlorellaConventional methodLab-scaleHCl38.892[50]
Spirulina platensisSingle stage extraction–transesterification processLab-scaleH2SO414.518[51]
Oleic acidMicrobubble technologyLab-scale7% H2SO426.37[20]
Chicken fat oilMicrobubble technologyLab-scale7% PTSA24.9[18]
SpirulinaMicrobubble technologySemi pilot scale3.3% PTSA10.01 ± 0.3This study
Table 5. Performance criteria comparison (response surface methodology vs. gated recurrent unit).
Table 5. Performance criteria comparison (response surface methodology vs. gated recurrent unit).
CriteriaConversion % Prediction Performance
RSM ModelGRU Model
R20.98440.9999
MAPE0.04650.00083
RMSE3.08320.0515
MAE2.68470.045
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Javed, F.; Saif-ul-Allah, M.W.; Ahmed, F.; Rashid, N.; Hussain, A.; Zimmerman, W.B.; Rehman, F. Kinetics of Biodiesel Production from Microalgae Using Microbubble Interfacial Technology. Bioengineering 2022, 9, 739. https://0-doi-org.brum.beds.ac.uk/10.3390/bioengineering9120739

AMA Style

Javed F, Saif-ul-Allah MW, Ahmed F, Rashid N, Hussain A, Zimmerman WB, Rehman F. Kinetics of Biodiesel Production from Microalgae Using Microbubble Interfacial Technology. Bioengineering. 2022; 9(12):739. https://0-doi-org.brum.beds.ac.uk/10.3390/bioengineering9120739

Chicago/Turabian Style

Javed, Fahed, Muhammad Waqas Saif-ul-Allah, Faisal Ahmed, Naim Rashid, Arif Hussain, William B. Zimmerman, and Fahad Rehman. 2022. "Kinetics of Biodiesel Production from Microalgae Using Microbubble Interfacial Technology" Bioengineering 9, no. 12: 739. https://0-doi-org.brum.beds.ac.uk/10.3390/bioengineering9120739

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop