Characterization of RHAC1 E RHAC2, RHAP
Immediate analysis
From the immediate analysis of the ash presented in Table 2, it can be seen that the two ash samples from the combustion in a mobile grate system (RHAC1 and RHAC2) have a low proportion of volatile material and fixed carbon, in contrast to the pyrolysis ash (RHAP), which showed high content of volatile material and fixed carbon. These results reflect the thermochemical processes where these ashes were generated. It can be seen that RHAC2, when compared to RHAC1, underwent a more intense combustion process, reducing both fixed and volatile carbon levels and presenting a higher ash content. For the RHAP samples, the fixed, volatile, and ash carbon levels are explained by the low process temperature and residence time (Dhyani & Bhaskar, 2019). The RHAP sample was obtained in the slow pyrolysis process, where it is observed the highest biochar formation (RHAP). Thus, the chemical bonds of the lignocellulosic material are broken slowly without significantly affecting the polymeric structure, favoring rearrangement reactions on the solid surface, secondary reactions associated with a longer residence time in the reactor, which benefits the formation of a carbon rich solid residue (Valadão, et al. 2021). The results described in Table 2 are in agreement with those determined by Fernandes et al. (2016), Flores et al. (2021), and Hessien et al. (2009) for RHAC and, Conz, (2015), Vieira, (2018) for RHAP.
Immediate analysis ( % p.p)
|
RHAC1
|
RHAC2
|
RHAP
|
Table 2
General characteristics of rice husks (% p.p ± SD)
Moisture
|
2.0 ± 0.1
|
1.5 ± 0.2
|
3.5 ± 0.3
|
Volatile Material
|
4.6 ± 0.1
|
2.7 ± 0.1
|
8.4 ± 0.1
|
Fixed Carbon *
|
5.5
|
2.0
|
44.5
|
Ash
|
87.9 ± 0.4
|
93.8 ± 0.2
|
43.6 ± 0.2
|
* Fixed carbon is obtained mathematically: 100 minus % moisture, % Volatile Material, and % ash. |
Where: SD: standard deviation. %(p.p): % Weight
X-ray diffraction (XRD)
Figure 2 shows the comparison of diffractograms for the three samples used in this study. It can be seen that the RHAP diffractogram shows only the amorphous halo characteristic of silica with a centralized peak, located between the 2θ positions of 15° and 30°, confirming that the silica present in this sample is amorphous. On the other hand, the RHAC1 and RHAC2 ash diffractograms show a decrease in the silica amorphism halo and show the presence of crystalline phases of silica. In the RHAC1 diffractogram, where the halo of amorphism decreases, it is possible to observe a peak of Cristobalite in the 2θ position of 22° with an intensity of up to 600 cps. In the RHAC2 diffractogram, where the halo of amorphism is significantly smaller, the cristobalite peak reaches the intensity of 2000 cps.
Table 3 lists the amorphism indices of RHAC1 and RHAC2, RHAP, and their respective crystallite sizes. When comparing these data with the immediate analysis (Table 2), it can be seen that the ash amorphism index increases in the following order: RHAP < RHAC1 < RHAC2, which is exactly the reverse order of volatile material and fixed carbon in the ash, which follows the following ascending order: RHAC2 < RHAC1 < RHAP. Considering the literature, an increase in temperature for the thermochemical process is one of the deciding factors for the silica transformation from amorphous to crystalline. In addition, when comparing the crystallite size increase for RHAC1 and RHAC2 ash samples, which present silica crystalline phases, the crystallite size is greater for the RHAC2 sample, as it has a lower volatile material and fixed carbon content. Possibly, these characteristics are a result of the burning conditions of the rice husks, mainly due to the temperature used. Some authors point out that crystallization starts at temperatures above 700°C (Cordeiro, 2009; Fernandes et al. 2016), while others suggest temperatures above 800°C (Chauhan & Kumar, 2013; Ghorbani et al. 2013; Soltani et al. 2015). Cordeiro (2014) proved the correlation between the amorphism index and the firing temperature, where it is possible to verify that the increase in process temperature causes a drop in the amorphism index. It should also be noted that, in addition to the firing temperature, there are other variables that directly affect the reactivity of these ashes, such as residence time, type of process, cooling, granulometry, and material composition. Finally, the XRD agrees with the immediate analysis of the samples, where RHAC1 and RHAC2 (combustion) showed lower levels of fixed and volatile carbon, higher amounts of ash, and also showed an increase in their crystallinity, as a reflection of the thermochemical conditions that were generated. These characteristics decrease the silica extraction efficiency. This situation does not apply to RHAP, generated by the Pyrolysis process (600°C), which presented 100% amorphous silica.
Table 3
Crystallinity index, amorphism content, and crystallite size of RHA
|
RHAC1
|
RHAC2
|
RHAP
|
Crystallinity index (%)
|
5
|
15
|
0
|
Amorphism content (%)
|
95
|
85
|
100
|
Crystallite size (nm)
|
25
|
30
|
-
|
RHAP showed no crystalline phase, so it is not possible to calculate crystallite size. |
X-ray fluorescence analysis (XRF)
Table 4 presents the composition of the oxides determined by the XRF for the RHAC1, RHAC2, and RHAP samples used in this study, and compares with other authors results. Similar composition is observed between RHAC1 and RHAC2 with Silica contents of 88.2% and 88.9%. Ash similar to the ones in this study were used in a work by Flores (2021), which were obtained from the combustion of the same type of RH. The results are also in agreement with those described by Fernandes (2016), for RHA obtained by combustion in mobile grate systems similar to the RHA of this study.
As observed in the immediate analysis, RHAC2 has a lower loss of ignition (L.O.I) content and higher levels of K2O. The potassium contained in RHA accelerates both the fusion of the particles and the crystallization of amorphous silica into cristobalite by lowering the melting point of the material (Ugheoke & Mamat, 2012). The RHAP ash has 71.3% in silica content, it is obtained at a lower temperature, thus leaving a greater amount of compounds that volatilize at higher temperatures (combustion), as indicated here by Chloride and L.O.I (Dhyani & Bhaskar, 2019). These results are in agreement with the Immediate analysis and with the XRD.
Table 4 X-ray fluorescence of RHA samples and comparative results
Oxide
|
RHAC1
|
RHAC2
|
RHAP
|
Fernandes et al. (2016) RHA MG
|
Flores et al. (2021)
|
SiO2
|
88.21
|
89.45
|
71.31
|
90.02
|
87.90
|
Na2O
|
<0.01
|
<0.01
|
<0.01
|
N.A.
|
N.A.
|
MgO
|
2.98
|
2.26
|
1.75
|
<0.01
|
1.26
|
Al2O3
|
<0.01
|
<0.01
|
<0.01
|
0.07
|
0.06
|
P2O5
|
1.13
|
1.09
|
0.93
|
0.34
|
<0.01
|
Cl
|
0.35
|
0.29
|
0.92
|
0.03
|
<0.01
|
CaO
|
0.80
|
1.03
|
0.87
|
<0.01
|
3.40
|
MnO
|
0.21
|
0.26
|
0.25
|
<0.01
|
N.A.
|
Fe2O3
|
0.06
|
0.08
|
0.08
|
0.01
|
0.33
|
K2O
|
1.69
|
2.78
|
2.20
|
0.81
|
5.58
|
SO3
|
0.16
|
0.12
|
0.17
|
0.07
|
<0.01
|
L.O.I
|
4.35
|
2.58
|
21.44
|
9.88
|
N.A.
|
L.O.I = loss on ignition; N.A.=Not analyzed; RHA MG = RHA obtained in combustion process with mobile grate.
Extraction of amorphic silica
Assessment of cavitations zones
Based on the tests with the aluminum foil strips, the regions of greatest wear on the paper were identified (Fig.3), which shows that the bath presents differences in ultrasonic intensity depending on the location in the device's bowl. To overcome these variations during the extractions the Erlenmeyers were attached to a support, which alternated positions every 10 minutes to ensure a uniform effect of the ultrasonic bath during the process.
Model development and statistical analysis
Table 5 shows the experimental design with the applied conditions and the results obtained in the extraction. The best silica yields were observed in tests 4, 6, and 8, ranging from 36.45–38.46%. The conditions of test nº 4 (molar ratio of 4.4 and time of 107 minutes) were chosen to replicate in the other samples and compare with the traditional method. Samples 6 and 8 did not differ from 4 with a significance level of 5%, wherein sample 6 used the highest test value for molar ratio and sample 8 used the highest test value for extraction time, demonstrating the combined effect of these two variables. Experiments 1 and 7 showed the lowest yields and were also the experiments with the shortest extraction times, which demonstrates the strong influence of time, because even increasing the molar ratio did not lead to a significant difference in yields. This influence is also observed when comparing samples 1 and 2, where the increase in time caused a significant improvement in yield.
Using the values of the dependent variables (income), the Shapiro-Wilk normality test was performed, obtaining a value of p − value = 0.08044 greater than 0.05, there is no evidence to reject the null hypothesis, so the data follows the normal distribution.
Table 5
Planning matrix with coded dependent variables, natural, and responses
Test nº
|
Coded variables
|
Natural variables
|
Yield (%)
|
X1
|
X2
|
X1
|
X2
|
Predictive
|
Observed
|
1
|
-1
|
-1
|
1.6
|
43
|
16.53
|
16.46f
|
2
|
-1
|
1
|
1.6
|
107
|
32.74
|
33.13b
|
3
|
1
|
-1
|
4.4
|
43
|
27.42
|
26.92d
|
4
|
1
|
1
|
4.4
|
107
|
38.50
|
38.46a
|
5
|
-1.4142
|
0
|
1
|
75
|
24.19
|
23.94e
|
6
|
1.4142
|
0
|
5
|
75
|
36.08
|
36.45a
|
7
|
0
|
-1.4142
|
3
|
30
|
17.90
|
18.29f
|
8
|
0
|
1.4142
|
3
|
120
|
37.09
|
36.82a
|
9
|
0
|
0
|
3
|
75
|
30.13
|
30.58c
|
10
|
0
|
0
|
3
|
75
|
30.13
|
30.71c
|
11
|
0
|
0
|
3
|
75
|
30.13
|
29.35c
|
12
|
0
|
0
|
3
|
75
|
30.13
|
29.89c
|
X1 = Molar ratio; X2 = Extraction time.The observed yields are the average results of each point of the response variable, in terms of the silica yield obtained. The mean values followed by the same lowercase letter in the column did not differ statistically by Tukey's test.
The statistical analysis of the effects, considering the yield as a response, confirms that the linear time is the most significant (1.6 times greater than the linear molar ratio effect) and positive variable. The linear molar ratio, interaction, and quadratic time effects were also significant, with the last two causing a decrease in silica yield. The quadratic molar ratio effect was not significant, with significance below p = 0.05, being taken out of the regression calculation. The analysis of standardized variables using the Pareto diagram (Fig. 4) demonstrates that the effects that are above the threshold value for the confidence interval are significant.
Multiple regression analyzes were performed to correlate the silica production responses (Y%) with the selected independent variables, using the polynomial illustrated by Eq. 7. The quadratic regression equation that leads to optimal silica production yield for this study, after eliminating the statistically non-significant coefficients, is shown as Eq. 8.
$$Y= -8.53984+5.11938{X}_{1 }+0.49439{X}_{2 }-0.0013{X}_{2}^{2}-0.0286{X}_{1}{X}_{2}$$
8
The statistical model presented in Eq. 8 must be predictive and meaningful to adequately represent the experimental data. The prediction and significance of the statistical model were evaluated using analysis of variance and Fischer's F test. Fischer's F test showed that the Fcalculated (Fcalculated =537.48) was more than 119 times the value of Ftabulated (Ftabulated=4.5). This implies that the use of the quadratic regression described in Eq. 8 is highly significant, with a significance level of 99%.
The calculation of the coefficient of determination (R2 = 0.993) shows that 99.30% of the variability in the response is compatible with the data predicted by the selected model. The predicted yield values of amorphous silica production are similar to the experimentally observed ones, showing that the experimental data follows a predictive model.
Through the normal probability graph of the residuals (Fig. 5), the assumption of normality of the errors was verified. They were normally distributed, and the interval was between − 1.2 and + 1.0, demonstrating the normality of the errors.
From the validated model, the response surface was constructed (Fig. 6). It can be observed that the increase in the reaction time and the molar ratio (n = nNaOH / nSilica) increases the silica yield, as discussed earlier. It is possible to verify that, at low molar ratios and high reaction times, yields close to those with high molar ratios are obtained, a fact that represents an economy both in the obtaining process and in the generation of residues. It is worth noting that the minimum ratio used in the process must respect the reaction stoichiometry for silica with NaOH (Eq. 1), which is 1, to avoid NaOH becoming the limiting reagent.
Silica extraction yield: traditional and modified method assisted ultrasound
The yield of silica extraction using ultrasound was evaluated in comparison with extraction by the traditional method. The results, presented in Table 6, show that the ultrasonic extraction method proved to be reproducible, with coefficients of variation (CV) lower than 5% for the different ash samples analyzed. It is also possible to observe the dependence of the matrix to be extracted, where samples with higher crystallinity showed lower yields. Effect also observed using the traditional method. A more complex matrix, such as that of RHAC1 and RHAC2, which present differences both in their compositions and in the variables of the thermochemical process applied to AC, such as operating temperature, residence time, and cooling, can result in different morphologies, structures, and reactivity (Fernandes, 2020), which ends up influencing the extraction yield.
The silica yields calculated considering the percentage of silica present in the ashes, determined by XRF, and the percentage of amorphism, determined by XRD, were higher than those initially calculated. They better represent the extraction efficiency by considering the real influence of the sample on the extraction, as it is known that only amorphous silica is extracted and not all ash is silica. It can be seen that the highest percentage of increase in yield occurred in RHAC2 due to its lower degree of amorphism, proving the influence of the rice husk combustion process conditions on the obtained ash and, consequently, on the extraction yield.
Table 6 Yields on silica (% ± SD), comparing the traditional and ultrasound method
Sample
|
Yield on sílica (%) Traditional
|
Yield on sílica (%) Traditional *
|
Yield on sílica (%) Ultrasound
|
Yield on sílica (%) Ultrasound *
|
RHAC2
|
33.9 ± 1.5%
|
43.4% *
|
24.3 ± 1.2%
|
31.1% *
|
RHAC1
|
63.3 ± 0.8 %
|
72.2% *
|
38.5 ± 0.8 %
|
43.9% *
|
RHAP
|
77.8 ± 0.9 %
|
85.7% *
|
66.5 ± 2.5%
|
73.3% *
|
Commercial Silica
|
98.2±0.5%
|
|
98.8 ± 0.6 %
|
|
SD: standard deviation* Yields were calculated considering the amount of silica determined by XRF and ash amorphism determined by XRD.
When applied to a commercial amorphous silica sample, the yields reached the recovery value of 98.8%, indicating minimal losses during the silica extraction process. The influence of the morphology in the extractive process was confirmed when the dissolution of commercial silica was reached within 20 minutes. The situation is also confirmed when compared with the traditional method. Samples with a higher degree of amorphisms, such as RHAP and commercial silica, showed less difference in yields between the methods, but this behavior was not so evident among combustion ash samples (RHAC), which suggests the influence of other components of the matrix in addition to the morphology of the silica to be extracted.
Considering the rice production of the South region, 2020/2021 harvest of 1,521,934 t/yearof which 20% are rice husk, which corresponds to 304,387 t/year.I In the conversion via the thermochemical process of CA, the average generation of 20% of RHA occurs (Fernandes, 2020), which would result in the generation of 60,877 t/year of environmental residue. Subjecting these ashes to the extraction process carried out in this work, considering yields of 73.3%, there would be the production of 44,622 t/year (3,719 t/month) of silica with purity above 95%.
The reaction temperature in this work remained at 60°C due to limitations of the ultrasound equipment, results by the traditional method with 90°C and bibliographic data suggest that better yields can be achieved with higher temperatures ( Mourhly et al. 2019; Zaky et al. 2008).
Characterization of the silica obtained
X-ray Diffraction Analysis (XRD)
Figure 7 shows the comparison of the XRD for the silicas obtained by the traditional and ultrasound methods. It is observed the characteristic amorphism halo of the silica with a centralized peak, located between the 2θ positions of 15° and 30°, and the absence of a crystalline peak, suggesting that the silica is in its amorphous form. Differences between the tested silicas were not evidenced by XRD, regardless of the origin of the ash or the extraction method adopted. The amorphism of the extracted silica denotes that only amorphous silica was extracted from the ashes, which may explain, in part, the lower extraction yields for the ashes that showed crystalline phases of silica (RHAC1 and RHAC2).
Scanning Electron Microscopy (SEM) and Energy Dispersive X-ray Spectroscopy (EDX)
Figure 8 shows the SEM micrographs for the silicas obtained from the SRHAC1 sample by the two methods, traditional and ultrasound. Both samples show small particles of different sizes and large agglomerates. It is not possible to identify significant differences between the silicas obtained by the traditional process (Fig.8b) and the silica obtained by ultrasound (Fig.8a). The images prove the presence of particles smaller than 10 µm, which can reach nanometric dimensions.
EDX analysis confirms the purity of the silica product, with silicon and oxygen representing 100% of the composition of the silicas obtained by both processes. The gold peak (2.10 keV) is due to the preparation of the samples, which were coated with gold for analysis, while the carbon peak (0.277 keV) is due to the carbon tape used to mount the sample. The elemental composition of the samples was confirmed by EDX, proving they are majorly composed of silicon and oxygen, as shown in Table 7 and Figure 9.
Table 7 EDX spectrum of silica extracted from RHAC1 Ultrasound and RHAC1 Traditional
Element
|
RHAC1 Ultrasound
|
RHAC1 Traditional
|
|
|
Weight%
|
Weight%
|
O
|
k
|
43.37
|
56.71
|
Si
|
k
|
56.63
|
43.29
|
|
|
100
|
100
|
K = electronic layer.
X-ray Fluorescence Analysis (XRF)
Table 8 presents the results of the chemical analysis by XRF. It can be seen that the silicas obtained by ultrasound and by the traditional method are similar to each other. The SRHAC1 obtained by the traditional method still showed an accentuated amount of Na2O remaining from the process. According to Kalapahty et al. (2001), metal ions (Na) are trapped in the silica gel network during gel formation, making it difficult to remove them when washing this gel. Sodium impurities can be removed by repeated washing with deionized water. In the work presented by Eliana dos Santos (2017), where the results are also listed in Table 9, it was verified that the amount of Na2O and SO3 was very high, and a subsequent washing was necessary after drying the silica. This washing with water avoids acid washing, leading to less consumption of chemical products.
The Fire Loss values of the ultrasonic and traditional SRHAC1 samples can be explained by the higher concentration of silane groups, which are located on the silica surface. These groups transform into siloxanes upon heating above 600°C, releasing H2O and transforming into a hydrophobic surface. To confirm this procedure, the SRHAC2 samples were heated to 600 °C in Muffle to remove the OH groups, as observed by the reduction of the LOI value and increase of oxides in XRF analysis.
In the work of Steven et al. (2021), a series of industrial applications were presented for silicas with purity above 95%, including civil construction (high-performance cement, concrete, pozzolan additive); chemical industry (cellulose and paper, detergent, petrochemicals, fertilizer, carrier); special materials (polymer, glass, ceramic, solar panel, graphene, zeolite, catalyst support); electronics (battery, capacitor, semiconductor); insulators (thermal insulators, refractory materials, flame retardants); adsorbent materials (moisture adsorbent in the form of silica gel, activated carbon); and pharmaceutical and health industry (as an active ingredient vehicle in formulations). This demonstrates that the silicas obtained in this work, except for the sodium contaminated traditional SRHAC1, are suitable for use. Considering only the oxides, all the silicas obtained in this work, except for the traditional SRHAC1, have SiO2 values above 95%.
Table 8 X-ray fluorescence of the silica obtained and comparative results
Oxides
|
Silica Ultrasound SRHAC1
|
Traditional Silica SRHAC1
|
Silica Ultrasound SRHAC2
|
Traditional Silica SRHAC2
|
Silica Ultrasound SRHAP
|
Silica before washing.Dos Santos (2017)
|
Silica after washing. Dos Santos (2017)
|
SiO2 *
|
92.97
|
83.97
|
94.87
|
94.99
|
95.59
|
76.50
|
95.52
|
(96.98)
|
(86.45)
|
(95.69)
|
(95.81)
|
(96.43)
|
|
|
Na2O
|
<0.01
|
9.90
|
<0.01
|
<0.01
|
<0.01
|
10.20
|
0.34
|
MgO
|
0.85
|
1.09
|
1.13
|
1.24
|
1.05
|
0.02
|
<0.05
|
Al2O3
|
1.48
|
1.42
|
1.65
|
1.49
|
1.60
|
0.12
|
<0.09
|
P2O5
|
0.21
|
0.39
|
0.46
|
0.43
|
0.46
|
0.19
|
<0.05
|
Cl
|
0.20
|
0.23
|
0.70
|
0.81
|
0.20
|
N.A.
|
N.A
|
CaO
|
0.08
|
0.07
|
0.16
|
0.17
|
0.16
|
0.03
|
<0.05
|
MnO
|
0.03
|
0.01
|
0.03
|
0.01
|
0.01
|
N.A.
|
N.A
|
Fe2O3
|
0.02
|
0.02
|
0.02
|
0.02
|
0.02
|
0.03
|
<0.05
|
K2O
|
<0.01
|
0.02
|
0.08
|
0.06
|
0.03
|
0.41
|
<0.05
|
SO3
|
<0.01
|
<0.01
|
<0.01
|
<0.01
|
<0.01
|
10.40
|
N.A
|
ZnO
|
0.02
|
<0.01
|
0.01
|
<0.01
|
0.01
|
N.A.
|
N.A
|
Eu2O3
|
0.02
|
0.01
|
0.02
|
0.01
|
0.01
|
N.A.
|
N.A
|
L.O.I
|
4.12
|
2.86
|
0.85
|
0.74
|
0.84
|
1.20
|
4.02
|
L.O.I = loss on ignition; N.A.=Not analyzed; * Values in parentheses were calculated considering LOI = 0%.
BET (Brunauer, Emmett, Teller) and BJH (Barrett-Joyner-Halenda) analyzes
Figure 10 shows the N2 adsorption-desorption isotherms for the silica samples SRHAC1, SRHAC2, and SRHAP obtained by the Ultrasound method, and the SRHAC1 obtained by the Traditional method. All curves exhibited type IV isotherms with hysteresis, as shown in Figure 10 (according to IUPAC classification). This indicates the mesoporous nature of the obtained silica products.
The BET and BJH results obtained for the SRHAC1, SRHAC2 and SRHAP samples are shown in Table 9. It is verified that the silica surface area for the SRHAC1 sample is 1.65 times greater when using the ultrasound method, instead of the traditional method, and the pore size is approximately half. This confirms the effect of ultrasound on the improvement in the textural properties of silica, which provides greater adsorption power and greater reactivity, in agreement with the studies by Peres (2018) and Ortiz (2016). The surface areas of the silicas are in agreement with the values found in the literature using the sol-gel method, such as Fernandes et al. (2016) and Yuvakkumar et al. (2014). It can also be seen that there are differences in the values of surface area and pore size for the silicas of the three ashes extracted by ultrasound. The smallest surface area observed was precisely for the silica obtained from RHAC2, indicating once again that the thermochemical method of ash generation and the ash matrix itself can influence not only the amount of silica extracted but also the textural properties of the silica obtained.
Table 9 Surface area, pore volume, and pore diameter
Sample
|
BET (m2.g-1)
|
Vt(cm3.g-1)
|
DBJH (nm)
|
SRHAC1 Traditional
|
275.89
|
0.84
|
8.94
|
SRHAC1 Ultrasound
|
456.92
|
0.77
|
4.69
|
SRHACP Ultrasound
|
347.41
|
0.83
|
7.23
|
SRHAC2 Ultrasound
|
265.81
|
0.73
|
7.98
|
BET - surface area determined by BET calculation, Vt= total pore volume determined by BJH, DBJH= Pore diameter determined by BJH.