3.1. Pilot Scale Fed-Batch Fermentation
The results obtained show that the CG fermentation process withstood the scale change, being possible to reach final values of ethanol and biomass concentration higher than those obtained in R5, as described in Cofré et al [
15]. This experience was divided chronologically into three stages: a batch stage (stage I, 63.5 hour), a fed-batch stage (stage II, 96 hour) and a post-batch stage (stage III, 92 hour) and a summary of them are presented in
Table 1.
During stage II, the biomass yield with respect to the CG consumed YX/S was equal to 0.03 g g-1 and for the case of ethanol YP/S reached a value of 0.43 g g-1 for a concentration of ethanol on the order of 5.4 g L-1. This final ethanol concentration value is within the order of magnitude of the results obtained in R5, being 14% lower. At the end of stage II, a CG consumption of 52.7% had been achieved, reaching a concentration value in the medium of 10.6 g L-1. On the other hand, the biomass reached the steady state after 155 h of fermentation and remained at a constant concentration value equal to 0.65 g L-1 until 250 h. At the end of stage III, all the CG had been consumed, reaching the maximum concentration of ethanol at the same time. Indeed, during stage III the remaining CG was completely consumed after 68 h (≈75% of the total time), a time that coincided with the maximum ethanol concentration that was reached and that corresponded to 8.15 g L-1.
The global behavior of PR was like that obtained in previously experiments at smaller scales (R5-5.5L batch fermentation), regarding the consumption of CG and production of ethanol during the post batch (stage III) the maximum concentration of ethanol reached in was equal to 7.58 g LO while in RP the maximum concentration was higher and equal to 8.51 g L-1, CG accumulation profile was remarkably similar except for the decrease registered in R5 at 68.5 h of fermentation time.
3.2. Mathematical Modeling
The considerations and assumptions made for the development of the mathematical model used correspond to an exponential variation of the feed flow rate over time and, in addition, it was considered that the specific population growth rate could be described according to the Monod model. (Monod, 1949). Thus, we have then:
where:
µm: maximum specific growth rate (h-1)
KS: saturation constant (g L-1)
The kinetic parameters for the Monod equation were obtained from two sources. First, the value of the maximum specific growth speed (µ
m) was assumed equal to 0.044 h
-1 according to the value indicated by Murarka et al. (2008) [
7]. Second, the value of the mean saturation constant (K
S) was considered equal to 6.5645 g L
-1. In the literature [
17] it is possible to find a K
S value for
E. coli, under anaerobic conditions with glycerol as substrate, equal to 14.6464 g L
-1, this value was evaluated during the adjustment of the experimental data to the theoretical model for fed-batch and the result obtained will be discussed later. The time interval chosen for the numerical integration, using the 4th order Runge-Kutta method, required to solve the system of equations was between 63.5 and 159.5 h, corresponding to the beginning and end of the feeding stage, respectively. Therefore, the initial values of reaction volume (V), substrate concentration (S) and biomass concentration (X) were chosen according to the experimental values for a fermentation time t = 63.5 h, which then corresponded to V
0 = 50 L, S
0 = 0.06 g L
-1 and X
0 = 1.19 g L
-1. It must be considered that PR fermentation was carried out considering a design value for F
0 equal to 0.25 L h
-1 which corresponds to a feeding time of 97 h. The fit of the theoretical model against the experimental data was carried out based on three scenarios, described in
Table 2. In a preliminary way, other alternative scenarios were considered, considering as a restriction, when discriminating between one result and another, the physical meaning of each parameter. The result of the adjustment for each of them was not satisfactory since the numerical values found did not comply with the proposed restriction.
The resolution of the system of equations was carried out by means of the 4th order Runge-Kutta method and the result of the parameter adjustment with respect to the theoretical model is presented in
Table 3, where it is possible to observe that the E1 scenario predicts a value for F
0 that corresponds to the one furthest from the design value.
For both E2 and E3, the F
0 value obtained from the numerical adjustment is closer to the design value. In the last case (E3), the parameter Ks was left free to obtain a value from the numerical adjustment. The volume variation profiles (V), CG consumption (SV) and biomass production (XV) were graphically represented with respect to the fermentation time corresponding to the experimental data and were compared with the simulated profiles from the numerical resolution of the model.
Figure 1a, 1b and 1c show the experimental data for PR contrasted with the curves obtained for the simulation of E2 and E3.
The existing discrepancy for biomass (XV) may be related to the µm value that was used with the same value (0.044 h-1) for all the cases presented. Apparently, the value of Ks does not present a significant influence for XV since the curves obtained for E2 and E3 are practically the same. Moreover, during the actual process, despite the accumulation of CG within the system, there is substrate consumption, but the biomass grows at a slower rate than the rate of volume increase and ethanol production, a fact that is evident when there is an observed negative variation of the biomass concentration compared to the positive variation of the ethanol concentration during fermentation. Taken to cell mass and product mass, a rate of increase for XV can be estimated approximately equal to 0.897 g h-1 and in the case of ethanol this rate could be set at 14.25 g h-1 so it can be deduced that the observed CG consumption is more related to ethanol production than to biomass growth. However, this condition is not considered in the proposed model and the accumulation of CG translated as an increase in the value of S, provides a value of µ greater than the real one, which results in an oversized prediction of the cell mass within the reactor.
With respect to the adjustment corresponding to the substrate (SV), as described in Cofré et al. [
15], for the experimental points corresponding to 130 h and 140 h it is where the greatest discrepancy occurs with respect to the E2 simulation, while the initial zone of points included between 63.5 and 115 hours is reasonably described by the model. When the comparison with the E3 simulation is made, a contrary behavior can be observed since the curve corresponding to the simulation passes precisely through the experimental value. In fact, the value provided by the model corresponds to 496.88 g while the experimental value is 497.35 g.
Figure 2 compares the residuals obtained for E2 and E3, defined as the difference between the real response value (y) and the response value predicted by the model (ŷ) [
18].
The influence of the value that Ks takes is notorious when analyzing the variation of the substrate where the final curvature of the SV profile is better described with an estimated Ks equal to 10.7703 g L-1, which is 64.1% greater than the value obtained in this work, while the value of F0 corresponding to the E3 adjustment is 12.18% greater than that obtained from the E2 adjustment. On the other hand, an extra scenario was evaluated leaving the corresponding values for F0, µm and Ks free. However, the values obtained for µm and Ks (0.0242 h-1 and 1.7217 g L-1, respectively) did not make sense according to what was observed throughout the development of this work and in the consulted literature, the explanation can be based on the compensation that must be created, between the values assigned to each parameter, when solving and minimizing the difference between the values delivered by the model and the real data.
From the above, E2 is the scenario chosen to describe the fed-batch zone of RP, because the value of F0 obtained by means of the numerical adjustment, and equal to 0.302 L h-1, is the closest to the value design (0.25 L h-1) with a difference of 20%. Although it is true, this difference can be debatable and a lower value for it would be desirable, the model under E2 allows us to adequately describe RP as a fermentation in fed-batch mode with exponential feeding. In addition, it should be noted that under E2 the weight of the adjustment fell only on F0 while µm was obtained from the literature and agreed with the µ values observed during the development of this work (approximately 0.03 h-1); while the value of Ks was obtained from the adjustment of real data generated in batch conditions at a smaller scale (reactor R5), which in a certain way would allow a validation of the value of this parameter.
3.3. Olive Pomace Oil Projections
In recent years, global interest in the production of renewable energy from biomass has increased significantly because it can be considered as a tool to help reduce global warming, due to its potential to replace the use of fossil fuels [
19]. However, for it to be considered as a competitive alternative to energy obtained from oil, the energy required for biomass production must not exceed the amount of energy that can be extracted from biomass. Edible vegetable oils such as canola, soybeans and corn have been used to produce biodiesel as a substitute for diesel, impacting other productive areas such as food production, agriculture, and land use, among others. Additionally, in the literature [
20], it has been reported that a major obstacle in the commercialization of biodiesel production from edible vegetable oil is its high production cost, which is due to the high cost of edible oil. In this sense, the above has boosted the search for alternative sources to obtain oil for biodiesel production [
21,
22], where it is possible to detach recovered cooking oil, microalgae, or agro-industrial waste such as olive pomace.
Rajaeifar et al. [
23] compared, through LCA (Life Cycle Assessment), the environmental impacts of biodiesel derived from Olive pomace oil (OPO) for B100, B20 and petroleum diesel No. 2, identifying the existence of environmental trade-offs between the categories of climate change and damage to resources. Furthermore, it was found that biodiesel required some remedial practices (in view of the agricultural and combustion stages) to be more environmentally friendly in all the mentioned damage categories. The results obtained showed that the B100 mixture was the worst in the Human Health and Ecosystem Quality damage categories due to the significant amount of emissions in the agricultural stage and the increase, of approximately 15%, in the amount of emissions of NOx during combustion compared to diesel. In contrast, in the Climate Change and Resources damage categories, the B100 mixture was the best and the No.2 petroleum diesel was the worst. Considering the results, authors suggested to adopt changes in agricultural practices to optimize the consumption of fertilizers and manure, as well as in the stage of use of biodiesel derived from OPO. An economic analysis presented in the same study showed a benefit / cost ratio of 1.45 revealing the economic viability of the production of biodiesel derived from OPO.
According to literature [
24], olive pomace has a proximal composition corresponding to: moisture (% in wet base) 6.86, ashes (% in dry base) 5.55, fixed carbon (% in dry base) 17.28 and volatile matter (% in dry base) 77.18, in addition to a higher heating value (HHV) of 22,032.95 kJ kg
-1 (dry base) 33.8% lower compared to petroleum coke 33,280.37 kJ kg
-1 (dry basis) [
25]. Olive pomace contains 2 to 12% oil, depending on the olive oil production technology[
26]. Regarding the use of OPO as raw material to produce biodiesel, it has been reported [
27] that it is possible to obtain a product with the same yields as other primary sources for the production of biodiesel, such as soybean oil or rapeseed oil [
23].
Each year, one million tons of olives are used in the production of olive oil and about 450,000 t of olive pomace are obtained [
26]. OPO is extracted from solid waste called two-phase pomace or “alperujo”. Alperujo is readily available (Spain generates about 4-6 million t per year) and is a low-cost raw material, but other important aspects must be considered such as its high humidity (50-70%) and low content of oil [
27,
28]. Specifically, with respect to OPO, as an example, in Greece the annual production of OPO is equivalent to 40,000 t [
29] and in Spain, for the 2005-2006 season, this value reached 56,000 t [
27]. Ten to twenty-seven kilograms of olive oil and 35 to 45 kg of olive pomace are obtained on average from 100 kg of olives; From 6 to 7.5 kg of olive pomace oil and 60 to 70 kg of dry olive pomace are obtained from 100 kg of olive pomace [
26].
From the results obtained in this work, it is possible to present a projection of the existing theoretical potential for the use of the generated OPO, taking Spain as an example, from the conversion of OPO to biodiesel, the obtaining of glycerin as a by-product of the previous process to the biotransformation of the latter through a fermentation process to produce ethanol and hydrogen. Taking into consideration the parameters contained in
Table 4, it is possible to estimate a production, expressed in (t year
-1), equivalent to: 53312 (biodiesel), 2968 (glycerol), 979.44 (ethanol) and 36.13 (hydrogen).