1. Phys-chemical properties of the OMWW
The physicochemical analysis of the studied OMWW (Table 2) shows these effluents have an acidic pH (4.81), indicating that OMWW is an acidic effluent. The same pH value was observed by Bouknana et al. (2014), Lee et al. (2019), Dehmani et al. (2020) and Lissaneddine et al. (2021). Table 2 shows a lower electrical conductivity (12 mS cm−1) than that found by Lissaneddine et al. (2021), Elayadi et al. (2021) and El Ghadraoui et al. (2021). The Chemical Oxygen Demand (80 g O2. L−1) is high and is also characterized by a predominance of toxic substances and the presence of PCs (Vuppala et al. 2021). Inorganic loadings such as potassium (4.8 g L−1), and sodium (0.3 g L−1).
Table 2
Physicochemical parameters of the OMWW
Parameters (Unit)
|
OMWW
|
pH
|
4.81
|
Electrical Conductivity (mS cm−1)
|
12
|
COD (g-O2 L−1)
|
84.1
|
Total phosphorus (mg L−1)
|
2.54
|
Orthophosphate (mg L−1)
|
1.69
|
Nitrite (mg L−1)
|
0.06
|
Sodium (g L−1)
|
0.3
|
Potassium (g L−1)
|
4.8
|
Ammoniacal nitrogen (mg L−1)
|
0.96
|
Kjeldahl nitrogen (NTK) (mg L−1)
|
215
|
Sulfate (mg L−1)
|
2769.23
|
Nitrate (mg L−1)
|
0.331
|
2. Characterization of the bio-sorbent
Figure 1(a and b) shows the FTIR spectrum of OP before and after treatment by hydrogen peroxide. FTIR spectrums were obtained in the wavelength range of 480-4000 cm−1. The FTIR spectrum of OP waste (Fig. 1a), the peak at 3457 cm−1, can be attributed to the -OH stretching vibration. The peak at 2925 cm−1 is due to C-H stretching vibrations in aliphatic CH, CH2, and CH3. The peaks at 1461 and 1645 cm−1 were attributed to the symmetrical stretching vibrations of CH2 and C=O groups, respectively (Wang et al. 2019). The FTIR spectrum of bio-sorbent (Fig. 2b) at 3416 cm−1 corresponds to the -OH of hydroxyl functional groups. The peak at 1628 cm−1 is attributed to the stretching vibration of carboxylic groups -C=O (Martín-Lara et al. 2008), the bond at 1380 cm−1 is attributed to the symmetrical stretching vibration CH2 (Allwar et al. 2020). The peak at 1033 cm−1 can be attributed to the bending vibration of the -OH group (Nasrullah et al. 2018). The FTIR spectra (Fig. 2c) between the bio-sorbent before and after PCs adsorption showed no significant difference.
Biosorption can be explained by the interaction between the functional groups present on the bio-sorbent and the CPs, considering two types of chemical and physical interactions. Biosorption involves different mechanisms such as coordination, complexation, microprecipitation or electrostatic attraction (Veglio and Beolchini 1994). Identifying these processes and the characterization of the active sites of the bio-sorbent are essential steps for optimizing the operating conditions in the development of the adsorption capacity of PCs.
SEM analysis allows visualization of the bio-sorbent pores. This analysis was performed to determine the elements present on the surface of the samples by EDX. The results obtained from olive pomace before and after chemical activation at the same magnifications are presented in Fig. 2.
Figure 2a shows that the surface of the olive pomace does not contain pores or voids. However, after activation, relatively more homogeneous pores with constant diameters appeared on the external surface of the bio-sorbent, improving the porosity; consequently, a high number of pores on the surface of the bio-sorbent.
The EDX analysis of the OP before and after activation allowed the elemental analysis of the elements present on the external surface. The analysis shows that both materials are composed mainly of carbon atoms and contain oxygen (Table 3a). After activation (Table 3b), Mg and Cu appear on the surface.
Şirazi et al. (2021) observed similar texture characteristics using OP activated by KOH. This shows that the OP's SEM surface before activation has a porosity less than that of the adsorbent after activation, which is similar to the present study results.
Table 3 EDX spectra of OP before (a) and after chemical activation (b)
(a)
Element
|
Atomic %
|
C
|
74.13
|
O
|
24.46
|
Si
|
0.26
|
P
|
0.26
|
S
|
0.17
|
K
|
0.54
|
Ca
|
0.18
|
Total
|
100
|
(b)
Element
|
Atomic %
|
C
|
70.50
|
O
|
29.10
|
Cu
|
0.22
|
Mg
|
0.18
|
Totals
|
100
|
The XRD patterns for pure and activated OP are shown in Fig. 3. The XRD patterns of OP before and after activation are very similar, suggesting that no alteration of the OP structure occurred during the preparation of the bio-sorbent. This is normal because we only did a chemical activation without carbonization. The two diffractograms showed peaks located at 2θ = 22.31° (Trache et al. 2014), 39.16°, 45.43°, 61.58°, 77.65°, 84. 47°, which is attributed to the native cellulose, while the other constituents of OP are mainly amorphous, which is explained by the appearance of the spectrum Schneider and Brenbner (1985) and Zhang et al. (2015).
The pH at which the zero surface charge is called the point of zero charges (pzc). The surface's electrokinetic properties are defined by pHPZC (Yagub et al. 2014). The plot pHf = f (pHi) is shown in Fig. 4. Results show that pHPZC of the bio-sorbent was around 8.64.
This value, superior to 7, indicates a concordance between this result, the acid's contents, and the surface's basic functions. Thus, the bio-sorbent surface is positively charged for pH values lower than pHPZC. As a result, the solution will become less acidic (consumption of protons from the solution), while for higher values, the surface is negatively charged, and the solution will become more acidic (release of surface protons). This is consistent with the result obtained by Şirazi et al. (2021) since the KOH-activated OP, pHPZC for adsorbent, was 8.0.
3. Process optimization
3.1. Optimization of preparation conditions for bio-sorbent
The preparation variables are presented in a design matrix; their ranges and responses (adsorption capacity (mg g−1)) are displayed in Table 4. The statistical software MINITAB is applied to compare and correlate the independent variables in Table 5. It was applied to develop polynomial regression equations, representing all quadratic expressions suggested by MINITAB. The equation expression was selected according to the sequential sum of the square model, which is based on the highest order of the polynomial where the model was not aliased and the extra terms were significant (Sahu et al. 2010; Garba et al. 2014; Salman 2014). There was a clear correlation between the experimental and predicted data, as indicated by the model R2 values 0.9745 for adsorption capacity. The final empirical model equation for the Yq response of the adsorption capacity is given by equation (eq. (8)):
Yq= 686 – 3.72X1 – 13.9X2 + 85X3 + 0.00583X12 + 0.190X22 + 58.3X32 + 0.0292X1X2 + 0.238X1X3 – 2.98X2X3 (Eq. 8)
Table 4
Minitab experimental design
Run
|
Level
|
Variables
|
Repones
q (mg g−1)
|
X1
|
X2
|
X3
|
1
|
-1
|
-1
|
-1
|
240
|
60
|
3.1: 1
|
273.21
|
2
|
+1
|
+1
|
-1
|
180
|
60
|
4.6: 1
|
123.21
|
3
|
-1
|
-1
|
-1
|
240
|
60
|
6.2: 1
|
326.78
|
4
|
+1
|
+1
|
-1
|
180
|
80
|
4.6: 1
|
416.35
|
5
|
-1
|
-1
|
+1
|
180
|
70
|
3.1: 1
|
340.71
|
6
|
+1
|
+1
|
+1
|
240
|
80
|
6.2: 1
|
486.07
|
7
|
-1
|
-1
|
+1
|
180
|
70
|
6.2: 1
|
305.35
|
8
|
+1
|
+1
|
+1
|
300
|
60
|
4.6: 1
|
292.85
|
9
|
-1.682
|
0
|
0
|
240
|
80
|
3.1: 1
|
551.78
|
10
|
+1.682
|
0
|
0
|
300
|
70
|
6.2: 1
|
510.71
|
11
|
0
|
-1.682
|
0
|
240
|
70
|
4.6: 1
|
332.14
|
12
|
0
|
+1.682
|
0
|
240
|
70
|
4.6: 1
|
332.14
|
13
|
0
|
0
|
-1.682
|
300
|
70
|
3.1: 1
|
489.00
|
14
|
0
|
0
|
+1.682
|
240
|
70
|
4.6: 1
|
332.14
|
15
|
0
|
0
|
0
|
300
|
80
|
4.6: 1
|
656.00
|
Table 5
Regression analysis (ANOVA) for response surface methodology
Source
|
df
|
Sum of Squares
|
Mean Square
|
F-Value
|
p-Value
|
Model
|
9
|
242316
|
26924
|
21.20
|
0.002
|
A-Time
|
1
|
72758
|
72758
|
57.28
|
0.001
|
B-Temperature
|
1
|
149644
|
149644
|
117.81
|
0.000
|
C-Ratio
|
1
|
83
|
83
|
0.07
|
0.808
|
A2
|
1
|
1624
|
1624
|
1.28
|
0.309
|
B2
|
1
|
1332
|
1332
|
1.05
|
0.353
|
C2
|
1
|
12563
|
12563
|
9.89
|
0.026
|
AB
|
1
|
1225
|
1225
|
0.96
|
0.371
|
AC
|
1
|
814
|
814
|
0.64
|
0.460
|
BC
|
1
|
3557
|
3557
|
2.80
|
0.155
|
Table 4 represents the variance (ANOVA) results of the full factional design using the statistical software MINITAB. The model F-value was recorded at 21.20, enlightening that the model was significant. On the other hand, high values of R2 (0.9745) and Adj-R2 (0.9285) were estimated, indicating that the selected model can describe 92.85% of the total variation on adsorption capacity data. The parameters having a p-value probability value of less than 0.05 are significant. In this case, A, B and C2 are the significant model terms. Values superior to 0.1 are considered non-significant terms.
3.2. The Pareto chart
The presentation of the relative importance of the main effects and their interaction was achieved using the Pareto chart, as shown in Fig. 5. The individual factors and the likely combinations obtained are represented as bars. To show whether the effects studied differ significantly from zero, the t-test was used (Abdel-ghani et al. 2016). For a 97.45% confidence level and nine degrees of freedom, the t-test value was found to be equal to 2.57. The values on the horizontal axes are the t-test values for each effect and their interaction (Carmona et al. 2005; Rathinam et al. 2011; Saadat and Karimi-Jashni 2011). Variables A, B and C2 have an absolute value greater than 2.57, which places them on the right side of the vertical line and makes them significant. While all other factors have an absolute value below the reference line, making them insignificant. By analyzing the Pareto chart, it can be seen that B, A and C2 have the greatest influence on the adsorption of PCs in the bio-sorbent. On the other hand, the effect of ratio (C) made the smallest contribution to the adsorption of PCs into bio-sorbent because it is far from the reference line of the Pareto chart.
3.3. Contour plots and response surface
The 3D contour plots and response surface show the interactive effects of several parameters on the response. The contour plots and response surface were acquired by Minitab software to indicate the adsorption efficiency of PCs considering independent factors such as activation time, activation temperature and ratio were shown in Fig. 6. The effect of activation time, ratio and activation temperature was investigated using response surface methodology (RSM) by Minitab software. From Table 5, the F-value of the model is 21.20 and the p-value is 0.023. The model terms were considered significant when the p-value was less than 0.05. The contour plot (Fig. 6a) for the interaction between activation time and ratio shows at time 300 min and ratio 6.2:1 (3) the interaction was declared significant. The graph of the interaction between activation temperature and ratio (Fig. 6b) shows that at ratio 3.2:1(1) and temperature 80°C, the adsorption of PCs is significant. The interaction between activation temperature and activation time (Fig. 6c) indicates at time 300 min and temperature 80°C, the biosorption is significant.
The optimal conditions for bio-sorbent preparation using response optimization were obtained at ratio = 6.2:1 (3), activation time = 300 min and activation temperature = 80°C.
4. Adsorption efficiency
4.1. Batch adsorption
4.1.1. Effect of initial concentration
The driving force is very important by the initial concentration of PCs to overcome all the mass transfer limitations between the two phases, liquid and solid (Liu et al. 2019).
The plot of adsorption capacity versus time at different concentrations was presented in Fig. 7. The maximum adsorption capacity is 789.28 mg g-1, at the initial phenol concentration of 4000 mg L-1 and time 120 min. This increase in adsorption capacity with increasing concentration was attributed to the fact that this concentration was the driving force in order to overcome the resistance in the mass transfer of PCs between the liquid phase (OMWW) and the solid phase (bio-sorbent surface), resulting in the increase of their transfer rate into the adsorbent (Lissaneddine et al. 2021). Comparing the results with different studies, the adsorption of PCs on activated carbon by Achak et al. (2009) reported a decrease in adsorption efficiency with an increase in the initial concentration. This phenomenon was attributed to activated carbon adsorbing more phenolic compounds at low initial concentrations than higher initial concentrations. On the other hand, in the study of Lissaneddine et al. (2021) on the use of SA-AC beads for the adsorption of PCs, an increase in adsorption capacity was observed with the increase in initial concentration. This is due to the increasing amount of PCs in the solution, which causes a rise in the mass transfer driving force (Liu et al. 2019).
4.1.2. Adsorption isotherm
When the adsorption process reached equilibrium, the adsorption isotherms were analyzed to describe the distribution of PCs between the solid and liquid phases (Liu et al. 2019). Several models have been proposed to describe this relationship during the adsorption of PCs OMWW, the main ones being the isothermal models of Langmuir and Freundlich. Some assumptions are applicable according to the Langmuir theory. First of all, all sites are energetically equivalent, i.e., the energy of the adsorption sites is considered equal, while they can contain at most one adsorbate molecule, forming a monolayer on the adsorbent surface. In addition, the surface of the adsorbent is assumed to be homogeneous, with the identification of the active sites. Moreover, there is no interaction between the adsorbed molecules (Swenson and Stadie 2019). On the contrary, Freundlich's theory considers the heterogeneous adsorption surface, while several layers of adsorbate molecules can be formed. Furthermore, this theory considers that interactions between adsorbed molecules occur (Dada et al. 2012; Osagie and Owabor 2015). Figure 8 shows the adsorption isotherm of PCs with the bio-sorbent. The adsorption capacity of PCs increased with the increase of PCs concentration until it reached 789.28 mg g−1 at the concentration of 4000 mg L −1. The driving force of mass transfer is the increase of the PCs concentration in the solution (Li et al. 2019).
From the linear regression data (Table 6) in Fig. 9(a and b), it is clear that the model with a correlation coefficient (R2= 0.9898), indicating that the active adsorption sites of PCs on the bio-sorbent are energetically homogeneous and a significant monolayer coverage of PCs on the outer surface of the bio-sorbent is formed without any interaction between the PCs (Mojoudi et al. 2019). Therefore, the Langmuir model represents the best model for the adsorption of PCs onto the bio-sorbent with a maximum adsorption capacity is 1125 mg g-1.
Table 6
Parameters of adsorption isotherms of PCs on bio-sorbent
Freundlich
|
Langmuir
|
KF (mg 1−(1/n) L 1/n g−1)
|
n
|
R2
|
KL (L mg−1)
|
qm (mg g−1)
|
R2
|
0.934
|
1.22
|
0.967
|
16.10−3
|
1250
|
0.9898
|
4.1.3. Effect of contact time
The effect of contact time allowed the essential parameters during adsorption. Adsorption time can vary depending on the nature of adsorbent and adsorbate, while it is affected by different parameters, such as concentration of PCs in the solution, temperature, and pH of the OMWW. In all cases, the equilibrium is reached at a given time; this is the state in which the increase of the contact time does not significantly increase PC adsorption efficiency. The equilibrium time varies with the studied adsorbents, such as SA-AC beads (120 min) (Lissaneddine et al. 2021), pomegranate seed (20 min) (Papaoikonomou et al. 2019), the raw clay (10 min) (Chaari et al. 2020), MgCl2-impregnated activated carbons (60 min) (Hamadneh et al. 2020) and hydroxyapatite-sodium alginate composite (Benaddi et al. 2021).
Figure 10 shows the sorption equilibrium obtained after about 150 min, then remains stable with an adsorption capacity value of 789.28 mg g-1. Adsorption can be considered a two-step process, a rapid initial adsorption step, attributed to many adsorption sites available on the bio-sorbent for PCs adsorption. This rapid step leads to a rapid increase in the amount of PCs accumulated on the bio-sorbent surface. Existing residual sites are then difficult to occupy by PCs, as repulsive forces develop between the PCs on the solid surface of the bio-sorbent and the bulk phase (Lissaneddine et al. 2021). In addition, PCs are small molecules and can diffuse into the internal pores (Papaoikonomou et al. 2019), thus reducing the driving force for mass transfer (Din et al. 2009). Obviously, during the slow step, the molecules encounter greater resistance to enter the adsorption sites (Lissaneddine et al. 2021; Achak et al. 2009).
4.1.4. Adsorption kinetic
The adsorption kinetic study provides important information to predict the adsorption process. In this work, two linear kinetic models studied the adsorption of PCs from OMWW by the bio-sorbent as a function of time, including pseudo-first-order and pseudo-second-order.
From Fig. 11, both models show clear linearity of the equation with a correlation value of R²=0.9456 for pseudo-first-order kinetics and a correlation of R²=0.992 for pseudo-second-order kinetics. Therefore, according to Table 7, the model with the experimentally calculated maximum equilibrium sorption capacity value approximately equal to that determined theoretically is the pseudo-first-order. Therefore, the adsorption kinetics of PCs on the bio-sorbent is of pseudo-second-order type. Furthermore, this model assumes that the sorption of PCs is mainly controlled by chemisorption, which includes valence changes through electron exchange or partitioning between the bio-sorbent and PCs. Therefore, external adsorption occurs more often than micropore adsorption (Abdelhay et al. 2017; Tao et al. 2019).
In a batch test, the bio-sorbent removed 789.28 mg g−1 of PCs. Compared to olive pomace treated with H3PO4+ 2MHNO3 (28.57 mg g−1) (Soudani et al. 2013), HCl + ZnCl2 (78.74 mg g−1) (Temdrara et al. 2015) and H3PO4 (110.30 mg g−1) (Soudani et al. 2017). The bio-sorbent of this study showed better removal efficiency of PCs.
Table 7
Values of characteristics constants for two models of adsorption kinetic
Pseudo-first-order
|
Pseudo-second-order
|
qethero
(mg g−1)
|
qeexp
(mg g−1)
|
Klag
(min−1)
|
R2
|
qethero
(mg g−1)
|
qeexp
(mg g−1)
|
Kb
(mg g−1 min−1)
|
R2
|
446
|
408.92
|
0.024
|
0.9456
|
446
|
423.21
|
60.64×10−4
|
0.9929
|
4.1.5. Effect pH
pH is a more important parameter in the adsorption process. The pH affects the sorption mechanisms onto the bio-sorbent surface and the nature of physicochemical interactions between the bio-sorbent adsorption sites and the PCs (Achak et al. 2009). In addition, pH can also affect the surface charge of the adsorbent (Aksu and Gönen 2004).
The adsorption capacity as a function of time at different pH values was shown in Fig. 12. A significant effect of pH on the adsorption of PCs was observed. It can be explained by the π-π electron donor-acceptor interactions between the aromatic ring of the PCs (electron acceptor) and the free oxygen of the surface basic sites (electron donor) (Moreno-Castilla 2004; Dabrowski et al. 2005; Hamdaoui and Naffrechoux 2007). The adsorption capacity of PCs on the bio-sorbent increases from 150 to 473 mg g−1 for pH values between 2 and 4, the maximum adsorption of PCs is reached at pH equal to 4.0. If the pH value is exceeded 4.0, there is a decrease in the adsorption capacity of PCs. This decrease can be attributed to the phenol ionization to phenolate ions form since the absorption of this later is prevented by hydroxyl ions presented on the adsorbent (Halhouli and Darwish 1995). The adsorption capacity is higher at acidic pH than at basic pH. The pH of the OMWW is acidic, which is advantageous in this study, so pH adjustment would not be necessary during effluent treatment. Therefore, OMWW can be used directly without pH adjustment in large-scale treatment systems, which minimizes costs during treatment. A study by Alwan (2008) showed that pH adjustment in industrial wastewater treatment plants influences the treatment efficiency, which leads to an increase in wastewater treatment cost.
4.1.6. Effect of temperature
Temperature is one of the most important parameters during the sorption process. Figure 13 shows the plot of adsorption capacity versus time at different temperatures. The maximum adsorption capacity is found to be 1307 mg g−1 at a temperature of 60°C. The increase in adsorption capacity with increasing temperature can be attributed to the increase in temperature promoting the polymerization phenomenon between the PCs, increasing adsorption efficiency. Furthermore, this increase in adsorption capacity can be associated with an increase in adsorbent swelling, allowing more active sites of the bio-sorbent to become available to the adsorbate (Ververi and Goula 2019). However, the temperature is a parameter that depends on other characteristics of the adsorbent and the solution (Moreno-Castilla 2004).
4.1.7. Thermodynamic studies
The thermodynamic parameters (ΔG0), (ΔH0) and (ΔS0) were calculated by using the following equations 10 and 11:
\(\varDelta {G}^{^\circ }=\varDelta {H}^{^\circ }-T\times \varDelta {S}^{^\circ }\)(Eq. 10)
\(ln{K}_{c}=-\frac{\varDelta {H}^{^\circ }}{R\times T}+\frac{\varDelta {S}^{^\circ }}{R}\)(Eq. 11)
Plotting ln(Kc) as a function of 1/T (Fig. 14) allows the calculated ΔH0, ΔG0 and ΔS0. The thermodynamic parameters for the adsorption of PCs onto the bio-sorbent are given in Table 8. The value of ΔG0 for all three temperatures was obtained negative, confirming the spontaneous and feasible adsorption of PCs on the bio-sorbent. In addition, the ΔG0 values for physisorption are less than -20 kJ mol−1 and chemisorption is greater than -80 kJ mol−1 (Turco et al. 2019). In this work, ΔG0 was -4229.38, -4355.49, -4732.51 and -5036.33 J mol−1, which indicated the adsorption of PCs with the bio-sorbent is chemisorption. In addition, ΔH0 and ΔS0 are other critical thermodynamic parameters that provide information about the adsorption process of PCs. In this study, the positive value of ΔH0 was 30.82 kJ mol−1, which demonstrates the endothermic nature of the adsorption process of PCs with the bio-sorbent. As for ΔS0, a positive value was obtained 49.26 J mol−1 K−1, so the increased randomness can explain it at the liquid/solid interface during the adsorption process (Sun et al. 2005). These outcomes are identical to those of Lissaneddine et al. (2021), who have found that ΔS° = 74 J K−1 mol−1 and ΔH° = 29 kJ mol−1.
Table 8
Thermodynamic parameters of adsorption of PCs onto bio-sorbent
ΔG0 (J mol−1)
|
ΔH0 (KJ mol−1)
|
ΔS0 (J mol−1 K−1)
|
293.5K
|
298.5K
|
313.5K
|
333.5K
|
-4229.38
|
-4355.49
|
-4732.51
|
-5036.33
|
30.82
|
49.26
|
5. Fixed-bed column
A continuous fixed-bed column adsorption process using a bio-sorbent as a sorbent was considered an alternative method to separate PCs from OMWW. Figure 15 shows the breakthrough curves for the adsorption of PCs on bio-sorbent in a fixed-bed column adsorption system at 60°C.
In Fig. 15, C0 is the initial concentration of PCs entering the column and Ct is the concentration of PCs leaving the column every 15 minutes. All PCs were adsorbed in the first 40 minutes, resulting in a low concentration of PCs in OMWW output. During the adsorption process, the concentration of PCs in the treated effluent gradually increases. The bio-sorbent sites are saturated and the adsorption zone progresses vertically in the column (Kundu and Gupta, 2005). In our work, the qe value reached 643.92 mg g−1 for Ci= 4251 mg L−1 and Q= 0.5 mL min −1, while the removal rate was 64% (Table 9).
Using Yoon-Nelson and Thomas's empirical models, the mathematical analysis of PCs adsorption study data on fixed-bed columns were analyzed. The model proposed by Yoon-Nelson considers the adsorption probability of each molecule to be directly proportional to the breakthrough probability of the adsorbate on the adsorbent and the adsorption probability of the adsorbate (Ajmani et al. 2020). The Thomas model is a theoretical model commonly applied in analyzing column adsorption data. It is based on the assumptions of the Langmuir isotherm model and the second-order kinetic model of the batch adsorption process (Ajmani et al. 2020). The Yoon-Nelson and Thomas model parameters for the PCs adsorption on the bio-sorbent in a fixed-bed column are presented in Table 10.
From Table 10, the Yoon-Nelson model correlates R2= 0.9572, demonstrating that the Yoon-Nelson model was approachable and could be used to describe the adsorption of PCs in continuous mode on the bio-sorbent. Table 10 shows that the mean absolute error (MAE) values were 0.0125 and 0.0104 for the Yoon-Nelson and Thomas model, respectively. In contrast, the two models' root means square error (RMSE) values were calculated as 0.12 and 0.10, respectively. The RMSE and MAE values of the Thomas model were lower than those of the Yoon-Nelson model. The Thomas model fit the experimental data; the results followed Langmuir kinetics (Chue 2010).
Table 9
Parameters of breakthrough curves of the packed bed column for PCs adsorption onto the bio-sorbent
Ci (mg L−1)
|
Q (mL min−1)
|
Vef (mL)
|
mtotal (mg)
|
R (%)
|
qe (mg g−1)
|
4251
|
0.5
|
60
|
500
|
64
|
643.92
|
Table 10
Parameters of Thomas and Yoon-Nelson model for the adsorption of PCs into bio-sorbent in a fixed-bed column
Model
|
Parameters
|
Value
|
Thomas
|
KTh (mL mg− 1 min− 1)
|
0.00144
|
q0 cal (mg g− 1)
|
547.78
|
R2
|
0.94571
|
MAE
|
0.0104
|
RMSE
|
0.10
|
Yoon-Nelson
|
KYN (min− 1)
|
0.3411
|
tCal (min)
|
57
|
R2
|
0.9571
|
MAE
|
0.0125
|
RMSE
|
0.12
|
6. Desorption
The desorption study aims to elucidate the adsorption process, the recovery of PCs, and the adsorbents' reuse. In the first adsorption cycle, the bio-sorbent was saturated with PCs with a removal rate of 64%. Then, the bio-sorbent was washed with distilled water and treated with 0.1M HCl solution for 120 min at room temperature. The adsorption process was duplicated several times, and the efficiency of the bio-sorbent in removing PCs was analyzed. At the second reuse cycle, the adsorption rate of the bio-sorbent is decreased from 64 to 30% (Fig. 16), which shows that the bio-sorbent is a good sorbent with excellent reusability and stability. Several studies on the desorption of PCs for adsorbent regeneration have been performed. According to Lissaneddine et al. (2021), who used Sodium-alginate-active carbon beads, the percentage of desorption of PCs was 58.5%. Yuney et al. (2020) used CuO-coated OP, reaching 77% of the PCs concentration desorbed.