3.1. Physico-chemical characteristics of materials
The proximate and ultimate analyses of coal, SRF, TDF, and SDF are presented in Fig. 2. The detailed characterization of waste materials can be found in the supplementary information (Table S1-S3). It can be observed that coal comprises 38% volatile matter content, while SRF, TDF, and SDF contain 74%, 65%, and 75.26% volatile matter content, respectively. Thus, the volatile matter content of the co-combustion mixture would improve when these wastes are included with coal co-firing, resulting in better and more stable flame (Mesroghli et al. 2009). Moreover, the lower ash content of SRF (2.6%) and TDF (5.86%) in comparison to coal (8%) suggested a reduction in waste and an increase in the heating value (Buyukada 2016). The ultimate analysis of TDF and SDF suggested higher C content (78.45% and 51.7%, respectively) than coal (55%), indicating a higher heating value. Though SRF has comparatively lower C content (45.16%), higher H and O contents can increase thermal efficiency during combustion. In addition, low S and N content of SRF (0.04%, 0.99%) and SDF (0.56%, 0.64%) suggested that blending of these materials in the co-combustion process could help reduce SOx and NOx emissions.
3.2. Model development and analysis
The experimental design matrix of the mixtures (Table 2) was utilized to develop the empirical models for the co-combustion process. The models were assessed based on the sequential model sum of squares, and the quadratic models containing linear and second order interaction terms (from SRF, TDF, and SDF proportion) were selected to represent the GCV and AFT of the mixed fuels. After eliminating the insignificant terms through stepwise regression and ANOVA analysis, the final model for GCV and AFT are shown in Eqs. (6) and (7), respectively (Myers et al. 2016).
$$\text{G}\text{C}\text{V}= 5296.95\text{A}+5270.42\text{B}+5489.01\text{C}+5292.07\text{D}+215.764\text{C}\text{D} +1110.74 \text{A}\text{C}$$
6
$$\text{A}\text{F}\text{T} = 1211.15\text{A}+1211.17\text{B}+1181.9\text{C}+1190.97\text{D}-42.2634\text{A}\text{B} +33.1815\text{A}\text{D}+102.549\text{B}\text{D}$$
7
Where A, B, C, and D are mass fractions of coal, SRF, TDF, and SDF, respectively, in the mixed-fuel blend. The positive signs of the terms in the equations depict synergistic effects, while the negative ones represent antagonistic effects (Zaib and Kyung 2022). The regression coefficient (R2) and the standard deviation were considered to evaluate the fitness of the models, as shown in Table 3. The correlation coefficients were obtained as 0.86 and 0.81 for GCV and AFT, respectively. Although higher correlation coefficient value is desirable, R² ≥ 0.8 is deemed acceptable (Bashir et al. 2010; Nassar et al. 2016). The smaller standard deviation and good R² demonstrated a reasonable agreement between the experimental and predicted results (Figs. 3 and 4) (Regti et al. 2017). The plots show the reasonable capacity of the model to fit the experimental data representing GCV and AFT as a function of coal mixed with SRF, TDF, and SDF.
Table 3
ANOVA of responses; gross calorific value (GCV) and ash fusion temperature (AFT) of fuel blends.
ANOVA for Reduced Cubic model
|
Response 1: GCV
|
Source
|
SS
|
d.f.
|
Mean Square
|
F-value
|
p-value
|
|
Model
|
69423.29
|
5
|
13884.66
|
18.04
|
< 0.0001
|
significant
|
Linear Mixture
|
64190.51
|
3
|
21396.84
|
27.79
|
< 0.0001
|
|
CD
|
2487.33
|
1
|
2487.33
|
3.23
|
0.0939
|
|
AC(A-C)
|
3124.8
|
1
|
3124.8
|
4.06
|
0.0636
|
|
Residual
|
10778.13
|
14
|
769.87
|
|
|
|
Lack of Fit
|
3296.01
|
9
|
366.22
|
0.2447
|
0.9675
|
not significant
|
Pure Error
|
7482.12
|
5
|
1496.42
|
|
|
|
Cor Total
|
80201.42
|
19
|
|
|
|
|
Std. Dev.
|
27.75
|
|
R²
|
0.8656
|
|
|
Mean
|
5342.35
|
|
Adjusted R²
|
0.8176
|
|
|
CV %
|
0.5194
|
|
Adeq Precision
|
14.3834
|
|
|
ANOVA for Reduced Cubic model
|
Response 2: AFT
|
Source
|
SS
|
df
|
Mean Square
|
F-value
|
p-value
|
|
Model
|
1341.35
|
6
|
223.56
|
9.11
|
0.0005
|
Significant
|
Linear Mixture
|
1042.12
|
3
|
347.37
|
14.16
|
0.0002
|
|
AB
|
125.09
|
1
|
125.09
|
5.1
|
0.0418
|
|
AD
|
78.56
|
1
|
78.56
|
3.2
|
0.0968
|
|
BD(B-D)
|
36.66
|
1
|
36.66
|
1.49
|
0.2432
|
|
Residual
|
318.85
|
13
|
24.53
|
|
|
|
Lack of Fit
|
213.85
|
8
|
26.73
|
1.27
|
0.4123
|
Not significant
|
Pure Error
|
105
|
5
|
21
|
|
|
|
Cor Total
|
1660.2
|
19
|
|
|
|
|
Std. Dev.
|
4.95
|
|
R²
|
0.8079
|
|
|
Mean
|
1198.7
|
|
Adjusted R²
|
0.7193
|
|
|
CV %
|
0.4132
|
|
Adeq Precision
|
9.9904
|
|
|
The analysis of variance (ANOVA) was performed to establish the significance of model terms based on F-values and the p-values (Table 3). The reduced quadratic model representing GCV (F-value = 18.04, p-value < 0.0001) and AFT (p-value < 0.05) remain highly significant. All the individual terms (except for interaction between SRF-SDF in AFT analysis) were significant (p value < 0.1). Moreover, insignificant lack-of-fit test for both models suggested no unusually large residuals from model fitting. The regression coefficient (R2 = 0.8656) and adjusted R2 = 0.8176 were in reasonable agreement, suggesting good model fitness and not overfitting due to additional model terms (Myers et al. 2016).
It is evident from the model that the two-way interaction of TDF with SDF or coal has significant positive interactions on the GCV of the composite mixture (Table 3). Like GCV, the quadratic model for AFT had also been highly significant (p-value = 0.0005), maintaining a non-significant lack-of-fit test (p-value = 0.4123). The difference between R2 (0.8079) and adjusted R2 (0.7193) could be attributed to the smaller sample size in the context of a number of model terms. The antagonistic interaction between coal and SRF can result in decreased AFT of the co-combustion mixture. Studies have suggested that the ash composition of biomass-derived materials can vary, and biomass fired boiler needs redressal (Li et al. 2013; Niu et al. 2016). A significant synergistic interaction between coal and SDF toward AFT was observed from the model. Apart from this, being particularly rich in SiO2 and Al2O3, sewage sludge also contains a higher amount of phosphorous, calcium, and iron, favoring increased AFT of composite fuel material (Zhang et al. 2013).
The quality of the models was further assessed through the coefficient of variance (CV), adequate precision, RMSE, SEP, bias (Bf), accuracy (Af), and MPE (%) (Myers et al. 2016). A coefficient of variance < 10% suggested the reproducibility of the model. The CV of 0.52% and 0.41% for GCV and AFT, respectively, was observed in this study. The adequate precision ratios of the signals with respect to noise for the GCV and AFT models were 14.4. and 9.9, respectively, which are much higher than the minimum desired value of 4.0 (Anderson and Whitcomb 2016). As shown in Table 2, the RMSE and SEP of 5.19 and 0.097% indicated the suitability of chosen mixture design model. Similarly, very low RMSE (0.89) and SEP (0.0744%) indicated excellent modeling ability for AFT of blended fuel. Furthermore, Bf and Af were close to unity for GCV and AFT, indicating a good agreement between the experimental and model predicted values. Besides, the MPE (%) of about 0.3362 and 0.2916 indicated better prediction accuracy for both GCV and AFT of blended fuel.
3.3. Effects of mixture composition on GCV and AFT
The main effects of coal, SRF, TDF, and SDF proportion in mixed fuel composition on GCV and AFT are shown in Fig. 3a and Fig. 4a, respectively. The main effect of the variables is measured as a change in GCV to deviation from the reference mixture composition, i.e., coal (96.25%), SRT (1.25%), TDF (1.25%), and SDF (1.25%). In general, the GCV of the mixture decreased with the increase in coal, SRF, and SDF content and increased with TDF. The GCV of the mixture decreased non-linearly from 5354 kcal/kg to 5285 kcal/kg, with an increase of coal content from 96.25–100%. The decrease in GCV with the increase in SRF and SDF was largely linear. The GCV of the SRF and SDF decreased from 5354 kcal/kg to 5274 kcal/kg (~ 1.49%) and 5292 kcal/kg (~ 1.15%), respectively, with the increase of these constituents by 3.75%. Contrary to other constituents of mixed fuel, the increase in TDF by 3.4% (from 1.6–5%) resulted in an increase in GCV of the mixture from 5354 to 5476 kcal/kg.
The combined effects of the SRF, TDF, and SDF on the GCV of mixed fuel containing 95% coal are shown in Fig. 3b. The GCV of the mixed fuel decreased with the increase of SRT and finally fell below 5300 kcal/kg as the SRT content approached 5% (turquoise contour band). Conversely, the GCV increased beyond 5450 kcal/kg with the increase in TDF composition in a mixed fuel. Together, the trace and contour plots postulate the highest GCV of TDF in the mixture, followed by SDF, coal, and SRT, which were in line with the experimental analyses. The order of GCV of the mixed fuel with different components were observed as: TDF (8723 kcal/kg) > SDF (5795 kcal/kg) > coal (4317 kcal/kg) > SRT (4375 kcal/kg).
Figure 4a shows the trace plots representing the individual effects of different mixed fuel components on AFT. An increase in AFT is observed with an increase in coal and SRF contents. The AFT increased with TDF and SDF content in mixture fuel compared to the reference blend, i.e., 96.25% coal and 1.25% each of SRT, TDF, and SDF. It is observed that the AFT of mixed fuel increases linearly from 1198°C to 1211°C with an increase in coal content from 96.25–100%. A similar but non-linear increase in AFT was observed with a change in SRT content from 1.3–5%. However, a decrease in AFT of the mixed fuel from 1198°C to 1182/1191°C was observed, with an increase in TDF/SDF from 1.3–5%.
Figure 4b shows the gradual increase in AFT with an increase in the relative proportion of SRT (red band towards the top of the cone). The mixture of TDF and SDF resulted in an AFT below 1190°C (blue band in contour plot). Furthermore, an increase in AFT to 1190°C is observed with an increase of SDF content up to 5% in the mixture fuel. Therefore, the trace and contour plots suggested that an increase in AFT is observed with the increase in coal and/or SRT, while decrease in AFT was observed with the increase of TDF and SDF in the mixture fuel. A 3.75% increase in coal and SRT increased the AFT by 13°C, while the same increase in TDF and SDF decreased the AFT by 16°C and 7°C, respectively. Hence, at reference mixture composition (96.25% coal, and 1.25% each of SRT, TDF, and SDF), further addition of coal to the mixed fuel decreased the GCV but increased the AFT; further addition of SRT decreased the GCV but increased the AFT; further addition of TDF increased the GCV but decreased the AFT; and further addition of SDF decreased both GCV and AFT. A desirable fuel mixture should possess high GCV and high AFT.(Adeleke et al. 2020) Therefore, numerical optimization was performed to achieve the mixed fuel blends with optimum GCV and AFT. It is noteworthy that co-combustion of food waste digestate with municipal solid waste was conducive to enhancing the combustion performance (Wei et al. 2021).
3.4. Optimization
The ultimate objective of the study was to propose a mixed fuel blend with reasonably high GCV and acceptable AFT. The constituents' individual and combined combustion patterns indicated that none of the four constituent materials exhibited high GCV and AFT simultaneously. Coal and SRT showed low GCV and high AFT; TDF possessed high GCV but low AFT, while SDF led low GCV and low AFT. Therefore, numerical optimization was performed to maximize both GCV and AFT provided that the coal consumption above 95% is minimized and the waste materials (SRT, TDF, and/or SDF) are utilized between 0 and 5%. The numerical optimization technique utilizes the desirability function to help estimate optimized mixed fuel composition satisfying the desirable characteristics: GCV > 5200 kcal/kg and AFT > 1200°C (Anderson and Whitcomb 2016; Myers et al. 2016). Figure 5 shows the optimized mixed fuel composition and its predicted GCV and AFT. The red dots indicate the factors, while the blue dots represent the responses. The downward ramp of coal (factor A) indicates the direction of the desirability function towards minimization of factor. Similarly, the upward ramp of GCV and AFT represents the maximization of these responses. It was predicted that mixing 95% coal with 3.84% SRT, 0.35% TDF, and 0.81% SDF would result in GCV of 5294 kcal/kg, and the AFT of the residual ash would be 1214°C. In the confirmation experiment, experiments were conducted at the optimized composition, which demonstrated the GCV of the mixed fuel to be 5307 kcal/kg and AFT as 1215°C (Table 4). Similarly, two other optimized mixed fuel compositions were predicted and experimentally validated. The difference between model-predicted and experimentally observed values was less than 1% in all cases, establishing the robustness and applicability of the models (Anderson and Whitcomb 2016).
Table 4
Optimum mixed fuel compositions predicted by numerical optimization using RSM.
|
|
|
|
|
GCV
|
AFT
|
Mixed Fuels
|
Coal
|
SRT
|
TDF
|
SDF
|
RSM
Predicted
|
Exp.
|
Error (%)
|
RSM
Predicted
|
Exp.
|
Error (%)
|
1
|
95.00
|
3.84
|
0.35
|
0.81
|
5294
|
5307
|
0.24
|
1214
|
1225
|
0.93
|
2
|
95.00
|
4.58
|
0.43
|
0.00
|
5291
|
5337
|
0.86
|
1209
|
1219
|
0.85
|
3
|
97.50
|
0.00
|
0.00
|
2.50
|
5289
|
5313
|
0.45
|
1209
|
1220
|
0.87
|