4.1. Cost and CO2 Emissions Tool
During the production process in the GM plant, carbon dioxide emissions are due to two processes. The first one is decomposition and the second one is combustion. The decomposition process does not depend on the type of fuel used, but on the raw magnesite ore quantity. The number of emissions due to decomposition equals 47,216.40 tn.
Figure 8 depicts the reduction in CO
2 emissions deriving from fuel combustion when biomass (SHP, OK or WS) substitutes pet coke energy by 30%, 50% or 70%. It is shown that at a 0% substitution (zero biomass), the fuel (pet coke) combustion emissions (47,103.08 tn) correspond to 50% of the total emissions (94,319.48 tn). At maximum substitution rate (70%), the CO
2 emissions from fuel combustion are reduced by the same rate and reach 14,118.74 tn, since biomass is considered carbon neutral, whereas the total CO
2 emissions are reduced by 35% (
Table 15). It is observed that the emissions produced during the combustion process with the co-firing of pet coke and biomass can reach up to 70% in cases of a 70% substitution. This proves that the targets for reducing emissions by 63% in the industry which have been set by the European Commission could be achievable.
The calculation of the total CO
2 emissions is shown in
Table 15 below, where the amount of CO
2 emissions is reduced to 61,335.14 tn for the maximum (70%) substitution rate with biomass. CO
2 emissions from decomposition are unchanged for each scenario examined.
Figure 9,
Figure 10 and
Figure 11 represent the production cost of DBM when using three different biomass types at different substitution rates and different CO
2 prices.
In
Figure 9, all biomass types show similar results at all three substitution rates. With the current price at 350 EUR per ton of pet coke and a price of 92.5 EUR per ton of CO
2 (prices between 21 and 23 December 2022) the best-case scenario is to substitute 70% of the pet coke energy with any kind of biomass, with olive kernel being the optimum choice since it has the highest LHV (17.64 GJ/tn) and the lowest price (65.25 EUR/tn). The production cost at a 70% substitution is reduced by up to 38.71% with this type of biomass. Wood sawdust also shows a similar downward track but with a lower LHV (9.68 GJ/tn); it is obvious that in order to achieve the needed amount of energy for DBM production, a larger amount of wood sawdust is necessary. Thus, production cost is higher when compared to the olive kernel case. In the case of using sunflower husk pellets as the substitution fuel, production cost is at its lowest value (272.82 EUR/ton produced) at a 70% substitution. Based on current prices, the use of OK can reduce production costs by up to 38.7% which, if it proves feasible in practice, is particularly attractive for any company.
Figure 10 demonstrates similar results when replacing a percentage of pet coke energy with biomass but with a decreased CO
2 price at 80 EUR/tn. The price of 80 EUR/tn reflects the average price of CO
2 tariffs for the month of November 2022. Olive kernel and wood sawdust demonstrate a decrease in production cost by 39.04% and 24.97%, respectively, at a 70% substitution, whereas the more expensive sunflower husk, which shows the lowest production cost at a 70% substitution (255.32 EUR/ton produced), demonstrates a decreased production cost by 12.08%. Again, with the current price at 350 EUR per ton of pet coke and a price of 80 EUR per ton CO
2, the best-case scenario is to substitute 70% of the pet coke energy with any kind of biomass, but with olive kernel being the optimum choice since it has the highest LHV (17.64 GJ/tn) and the lowest price (62.25 EUR/tn). Judging based on the application of the optimal proposed solution (substitution with OK by 70%), it is noticed that the reduction in production costs can be particularly high even at lower values of CO
2. In fact, a slightly higher percentage of reduction (39%) is observed here compared to the values presented in
Figure 8 (38.7%).
In
Figure 11, where CO
2 price per ton is increased at 100 EUR/tn, all types of biomasses show a profit from 15.04% when using sunflower husk pellets to 38.53% when using olive kernel at a 70% substitution rate. In the case of wood sawdust, the same trend is observed with around a 26.27% reduction in production cost. This picture can also be characterized as an illustration of future production costs for GM, as according to the European Commission and the targets it has set for the reduction of carbon emissions, the prices of CO
2 allowances will continue to rise and will soon exceed 100 EUR/tn.
Figure 12 illustrates a comparison of the production costs for GM between November 2021 and November 2022. More specifically, the rise in production costs during this year is shown, which is due to both the increase in fuel prices and the increase in the prices of CO
2 tariffs.
Table 16 summarizes the price variation during the year.
The leap in GM production costs is clear. Comparing the costs without the use of biomass, an increase of more than 50% is observed. Corresponding increases are also observed under the optimal substitution scenario (70%) by 41%, 46% and 44% for OK, SHP and WSD, respectively. Indicative is the conclusion that no matter how much fuel prices increase, the use of olive kernels and substitution by 70% will always be the most economically optimal solution. In addition, it should be emphasized that for the scenario that has already been effectively implemented in GM (50% substitution with SHP) the increase is estimated at 49%. Examining the comparison, it is found that in just one year the production cost increased by 35%. This is mainly due to the sharp increase in the price of pet coke, but also to the continuous increase in CO2 tariffs. Obviously, biomass supply prices also increased (by about 45% for each type of biomass); however, the increase is clearly smaller, and the estimated costs are significantly lower. This diagram proves that an increase in prices is expected for the coming years, making the use of biomass to replace part of the pet coke more and more imperative.
Evaluating the conditions, the global political and economic situation, as well as the market trend, it is estimated that both fuel prices and CO
2 tariffs will continue to increase.
Figure 13 shows the cost of production of GM, adopting a worst-case scenario. The following table (
Table 17) summarizes the values estimated to be used in the worst-case scenario as well as their increase compared to the existing ones.
The production cost for the worst-case scenario reaches 386.79 EUR/ton of MgO produced. In this case, the substitution of fossil fuel (pet coke) is considered imperative from an economic point of view. Substituting 70% with olive kernel reduces the production cost by about 40%, estimated it to be about 234.39 EUR/ton produced. Similar, although not as impressive, reductions in costs are also observed from the use of wood sawdust (24%) and sunflower husk pellets (10%). The consecutive increases observed make the worst-case scenario a not-so-improbable case. Taking into account the results, the cost of production can rise to close to 390 EUR/ton produced which can make the production process economically unviable, essentially making the use of biomass a requisite for the sustainability of the system.
4.2. CFD Tool
The CFD tool allows many combinations of fuels, e.g., substituting a given percentage of coke (energy basis) with PSD, or with OK or with a blend of both types of biomasses. Since analyzing all the possibilities is not feasible, we present the results of the replacement using one or the other biomass as way of example.
Figure 14 shows the temperature profiles obtained in the base case (100% coke) in the longitudinal plane of the kiln. The bars shown in the plot are just visual guides to help understand the plots. They represent the distance from the burner. In
Figure 14, there is a comparison of the base case with different substitution rates (SR) using PSD or OK. The maximum temperature reached with coke is above 1700 °C which is in agreement with measurements and simulations by other authors [
27,
28,
29,
30,
31,
32,
36]. When coke is partially replaced with biomass, there is a decrease in the average temperature that is more pronounced the higher the substitution rate is. This is in agreement with the literature on alternative fuels [
33]. The effect is very noticeable for PSD in which the maximum temperature decreases by 500 °C for SR = 30–40% and another 100 °C for SR = 50%. The OK behaves better as the maximum temperature decreases by 200 °C in the three cases studied. Although the substitution is carried out in an energy base, the three fuels have different carbon contents, volatiles contents and different kinetics for devolatilization and subsequent char oxidation. These differences explain the changes in the profiles of temperature and highlight the necessity to perform detailed simulations instead of relying on mass and energy balances.
Comparing the results in
Figure 14 with the proximate analysis of the fuels (
Table 7), the temperature inside the kiln is positively correlated with a fixed carbon and heating value, while it seems to be negatively correlated with the percentage of volatile matter.
An interesting observation is that in the base case the highest temperature between 10 m and 20 m is concentrated near the walls. This is common in rotary kilns, as discussed in [
34]. These can be related to the fly of fuel particles (see Figure below). Instead, in the case of PSD, it is very homogeneous in the transversal direction and, in the case of OK, the maximum is in the center of the kiln.
One advantage of the lower temperatures achieved with biomass is that NO
x production should be lower. In Figures below one can see reductions of up to 50% in NO
x emissions using biomass. This is because thermal NO
x is very sensitive at peak temperatures and increases exponentially for temperatures above 1600 °C. Of course, this is an advantage given that the minimum temperature for magnesite calcination (CCM) is reached. For example, to achieve a reduction of 50% in NO
x, a decrease of 400 °C is required (see Figure below). The behavior of the particles close to the burner can be related to the velocity profiles (
Figure 15)
The values of velocity (
Figure 14) are similar to those reported in the literature [
30] of pulverized coal. They show that the velocities of gases are lower when using biomass. The maximum velocity drops by 30%, approximately. This is because the higher size of the biomass particles, which are more difficult to grind, put an obstacle to gas propagation. On the other hand, the shape of the profile is very similar in all cases. This means that the hydrodynamics of the gases is not affected greatly by the presence of the fuel particles, despite the fact that the fuel dramatically changes the thermal picture.
The oxygen profiles for coal in
Figure 16 are very similar to other CFD studies [
30] which provides further validation for our work. The initial consumption of oxygen is faster in the two cases of biomass than when using only coke. Comparing the base case with PSD, the pattern is very similar between 10 m and 25 m in which a gradient of concentration is clearly visible. This is in contrast with the case of OK where the concentration is homogeneous from 7m onwards. This is in line with the homogeneity of the temperature profiles discussed above being more homogeneous for the OK. Note that the correlation between oxygen concentration and temperature is also tight in the case of PSD, since both quantities decrease fast upon increasing the substitution rate.
The reason for the faster conversion in the two cases of biomass is most likely related to volatile matter. Devolatilization is the first process in fuel combustion, and the carbon content of the volatiles will react with the oxygen. Indeed, there is a negative correlation between the percentage of volatiles in the three fuels and the extension of the red area in
Figure 14. That is to say, the more volatiles, the faster the initial oxygen depletion.
The volume where the heat of reaction is concentrated (see
Figure 17) is correlated with the region where the most oxygen consumption takes place. It is also consistent with the higher reactivity of biomass due to the large number of volatiles, while its calorific value is much lower than that of coke. The higher the amount of biomass, the higher the volume where the heat of reaction reaches significant values.
The heat of reaction of the OK reaches up to 500 W while that of pine is at most 70% of that value.
The precedent results show mostly the behavior of the gas phase. It is also possible to visualize the trajectories of the fuel particles to parse their behavior. In
Figure 17, the trajectories of fuel particles are shown colored according to temperature. Note that only a sample of the particles are represented, otherwise the plot would appear overcrowded and featureless. A common feature of all the results is that the trajectories are very chaotic in the first 5 m from the burner. A close inspection shows that the particles recirculate in that zone several times before starting to travel along the kiln, and finally dragged by the air. From that point onwards the particles follow more ordered trajectories. It is also noticeable that around the 5 m mark, there is a widening of the trajectories of particles that were previously close to the axis of the kiln. In some cases, the trajectories appear braided, while in others they look more or less parallel.
The behavior of the particles close to the burner can be related to the velocity profiles (see above). In all cases, a fraction of the air is ejected directly towards the wall of the kiln. The effect is more apparent in the cases with biomass but is present for coke also. Thus, a fraction of the particles travels directly towards the wall and experience rebounds which gives rise to the trajectories in
Figure 18.
Finally, there is a close correlation between the gas phase temperature (
Figure 14) and the temperatures reported in
Figure 19. In addition, in the latter is clear that the OK gives rise to more homogeneous temperatures than the PSD, and the average temperature is less affected by the substitution rate.
The above results illustrate the use of the CFD simulation to investigate the conditions inside the kiln. Turning to the topic of emissions, the composition of the fumes can be calculated at the exit to obtain an estimation of the compounds of interest.
Figure 19 shows the result of a parametric study with different substitution rates and three blends of biomass: 100% PSD, 100% OK and a blend at 50% of each type (PSD/OK). The reduction in emissions for NO
x and SO
2 along with data for exit temperature and oxygen level at the exit of the fumes are reported.
A reduction of up to 50% can be obtained by sacrificing 400 °C degrees of temperature using pine, which is probably too much to keep the process inside the kiln. Great reductions of NO
x have been reported for similar kilns when using a biomass [
23,
53]. Using olive kernel, it is possible to achieve substantial reductions in NO
x emissions without affecting the temperatures greatly. In fact, the drop of temperature between the base case and case SR 50 for OK is 20 K. Using a blend of both biomasses produces an intermediate result in which substantial reductions of NO
x come at the cost of upsetting the temperature of the process in a substantial way.
Since the set of simulations reported have the same air conditions, the behavior of the oxygen level measured at the exit changes. It goes in different directions for pine and the olive kernel, decreasing with respect to the base case for the former and increasing for the latter. If it is desired to keep the 2% level, it is possible that this is easily achieved by slightly modifying the amount of air supplied to the burner.
The authors of [
54] study the reduction in NO emissions achievable by controlling the operational parameters only. Using an analysis that considers energy balance and kinetic calculations results to a decrease of 15% in NO
x. This is slightly above the minimum saving reported in
Figure 19, so the use of biomass seems a very promising route to obtaining important savings in controlled emissions.
Regarding the accuracy of the results, the NOx emissions are highly dependent on the properties of the biomass nitrogen content and its distribution among volatiles and char. Since in practice is not possible to know these values, the results shown should be taken as an indication of the trends expected. On the other hand, the results concerning temperature, the main control variable, are robust. This is because they depend on the average values of carbon content in the fuel, instead of its distribution, and on the mixing inside the kiln, which is well accounted for by the hydrodynamic calculation. Indeed, comparing different simulations in the literature, the values of temperature are very similar, the differences being due to changes in the process or operating conditions.
Comparative and parametric studies of this kind can be easily and rapidly made with the ROM. Thus, it is possible to determine the most appropriate scenarios from the point of view of emissions and the process conditions.
Both tools simulate the production process and contribute to the understanding of its operational parameters. However, the flexibility tool can also contribute to the decision-making process by providing a solution to issues related to the operation, cost and environmental footprint of the unit. The application of this tool can identify with relative precision and give an initial approximation about the changes that could result from the variation or change of an operational parameter (e.g., fuel change) in costs and emissions (reduction or increase). Therefore, a company using the flexibility tool can test various scenarios at a simulation level, and based on the results decide which of them will be the most beneficial from an economic and environmental point of view so that this can be implemented in a real production process of the demo site.