3.1 Surface Morphology of Different Structures
As we can see from Fig. 1a ~ d, it can be clearly observed that CPS and CPSG have a 3D porous structure. In contrast, PS only forms a 2D film, and CP forms a 3D structure without porous structure. At first, as a water-soluble, soft, and long-chain polymer, PVA has widely been applied in the preparation of permeable membrane, aerogel, and hydrogel due to its excellent film-forming and excellent chemical properties (Sonker et al. 2018). However, even if we added any foaming agent in this research, there was still no porous structure. This may be because we utilized freezing instead of freeze-drying in this experiment. So, in the subsequent thaw process, its 3D structure collapsed. In addition, bubbles generated by the foaming agent could not form porous structure in the final solidification due to the collapse of the overall structure.
In contrast to CP, CPS and CPSG formed a good 3D structure. The main reason is the addition of cellulose. It is well-known that cellulose consists of a well-organized crystallinity region that contributes to the strength, high stiffness, and amorphous region, which makes fiber more flexible (Kalia et al. 2011; Lavoine et al. 2012). Cellulose is added to the system as a reinforcing component filler, forming a uniform network and stable 3D structure in the mixed solution with PVA through hydrogen bonding and supramolecular interaction. Song et al. prepared a cellulose-PVA hydrogel material. They found that after adding cellulose, the mechanical strength of the composite material was significantly enhanced, and the structural stability of the material was improved (Song et al. 2020).
Foam molding technology has developed quickly in the past few decades. As an anionic surfactant, SDS has been widely used to foam cellulosic suspensions (Alimadadi and Uesaka 2016; Hou and Wang 2017; Viitala et al. 2020). Lee et al. successfully prepared porous 3D cellulose foam using 3D printing and foaming technology, with good mechanical properties and interesting structure (Lee et al. 2021). In this research, cellulose/PVA/EG 3D porous foam was physically cross-linked. Then the PEG and SDS were washed away in the subsequent washing process (Fig. S4) and initially utilized SDS as a foaming agent and surfactant to prepare foam material. Compared with CPS and CPSG, SDS is the key to porous forming. Foam is a complex gas/liquid dispersion system. In the effect of high shear force, massive bubbles are generated and accumulated in a trace amount of surfactant solution to form a liquid foam slurry (Weaire and Phelan 1996). As shown in Fig. 1g, those irregular bubbles fix the cellulose and EG particles between the gas and liquid phases (Hjelt et al. 2011). In addition, the liquid foam can be prepared by rearranging SDS molecules at the surface of air bubbles entrapped in a suspension (Liu et al. 2018). The hydrophobic area of SDS extends out of the bubble surface, and the hydrophilic area combines with water to form a surfactant layer. The negative charge of cellulose and EG caused electrostatic repulsion with SDS. So, they are repelled between the bubbles to form a water layer containing cellulose and EG. Although it has the same charge as cellulose and EG, they still adsorb to hydrophobic sites to form aggregates as SDS is a strong surfactant(Tucker et al. 2012).
3.2 Influence of Different Factors on Adsorption Performance
Adsorption is a surface phenomenon in which the solid surface of the adsorbent is bonded by adsorptive molecules (gas/liquid/solid phase). The whole adsorption process is mainly divided into three steps: (1) Film diffusion is also called external diffusion, in which the adsorbate is transported to the external surface of the adsorbent in the space; (2) Intraparticle diffusion (IPD), means that the adsorbate diffuses from the external surface of the adsorbent into the pores; (3) Surface reaction, means that the adsorbate is adsorbed on the inner surface of the adsorbent (Ho et al. 2000; Tan and Hameed 2017). In this case, the mass transfer effect needs to be considered. In the adsorption process, the first two steps are mass transfer, and the last step is the reaction step. The adsorption efficiency is determined by the adsorption resistance. Any step of the adsorption resistance in the above steps will have a huge impact on the total adsorption efficiency (Amanullah et al. 2000).
At first, transmission resistance is affected by many factors, including the type and structure of the adsorbent. We investigated the adsorption performance of different samples at relatively lower dye concentrations (25 mg/L). As seen in Fig. 2a, CPSG has the highest removal rate of dye (96.21%) and adsorption capacity (12.03 mg/g) compared to CP (16.51%, 2.06 mg/g) and CPS (65.02%, 8.13 mg/g). This is mainly due to the porous structure of CPSG, and MB molecules easily diffuse from the external surface to the internal surface. Moreover, EG provides more active sites and significantly improves the adsorption performance.
In addition, the effect of different initial dye concentrations (25 mg/L – 250 mg/L) on the adsorption performance was investigated. It can be seen from Fig. 2b that as the initial dye concentration increases, the adsorption capacity increases significantly from 12.03 mg/g to 110.81 mg/g, which is possibly caused by concentration polarization. At higher dye concentration, concentration polarization accelerates the diffusion rate of MB molecules on the surface of the adsorbent, and more dye molecules are adsorbed (Eltaweil et al. 2020a).
The pH value is also one of the important factors influencing adsorption performance. The research results show that the removal rate and adsorption capacity of MB increase with the increase of pH value, which is consistent with the relative researches (Fig. 2c) (Wang et al. 2018; Abd-Elhamid et al. 2019; Eltaweil et al. 2020a). At low pH, due to protonation, the positive charge on the surface of the material increases, which causes repulsion between cationic dyes and positively composite, resulting in low adsorption efficiency (Premarathna et al. 2019).
Moreover, the influence of different temperatures on the adsorption performance was investigated. It can be seen from Fig. 2d that as the temperature increases, the dye removal rate and adsorption capacity of the adsorbent decrease significantly, indicating that temperature is also an important influencing factor of the adsorption performance. As the temperature increases, the thermal mobility of dye molecules increases, which may accelerate their transport to the surface of the adsorbent. However, at the same time, it may also accelerate the desorption of dye molecules from the surface of the adsorbent (Sabarish and Unnikrishnan 2018).
3.3 Adsorption Kinetic of Adsorbent in MB
There are two possible interactions between adsorbate and adsorbent, namely chemical and physical interaction. Among them, chemical interaction is also called chemical adsorption in that there are chemical or covalent bonds between adsorbate and adsorbent by sharing or transferring electrons. On the other hand, the main force that dominates physical adsorption are van der Waals forces mainly(Dąbrowski 2001).
It is very important to figure out whether an adsorption phenomenon is a chemical or physical adsorption, which can help us understand its adsorption mechanism. Adsorption kinetics reveals the relationship between adsorption efficiency and time and is a good tool to analyze the adsorption rate. This research chooses the three most commonly used kinetic models (Pseudo-1st order, Pseudo-2nd order and IPD model) to study the adsorption rate. The equation of Linearized form is shown in Table S1.
Batch adsorption experiments were performed with different concentrations of dye solutions at 25 ℃, and the samples were tested every 30 minutes. The adsorption kinetics curves and parameters are shown in Fig. 3 and Table 1. It is obvious that the Pseudo-2nd order is more suitable than the Pseudo-1st order. It reveals that the adsorption process depends on the number of active sites (Wang and Guo 2020a).
However, some research has reported that the kinetic data near equilibrium may have serious deviations if a Pseudo-2nd order kinetic model is applied (Simonin 2016). Therefore, other models are needed for further study on the adsorption rate.
The IPD model is used to determine the rate-limiting steps of the entire adsorption process. As shown in Fig. 3c and Table 1, \({k}_{i-1 stage}>{k}_{i-2 stage}>{k}_{i-3 stage}\) indicates changes in diffusion rate, revealing that the adsorption process mainly follows three steps: (1) The first stage is the diffusion of MB molecules from the body to the external surface of the adsorbent; (2) Then, MB molecules diffuse through the pores and diffuse into the internal surface; (3) The last stage is near equilibrium stage. In this stage, the balance of adsorption and desorption is reached.
Table 1
Parameters of adsorption kinetic
Models
|
Parameters
|
25 mg/L
|
50mg/L
|
100 mg/L
|
150 mg/L
|
200 mg/L
|
250 mg/L
|
Pseudo-1st order
|
qe,exp(mg/L)
|
12.03
|
24.91
|
48.94
|
73.49
|
90.90
|
110.81
|
qe,cal(mg/L)
|
4.49
|
5.55
|
12.40
|
15.92
|
42.84
|
98.50
|
k1(min− 1)
|
-0.8963
|
-0.8271
|
-0.4998
|
-0.4887
|
-0.5565
|
-1.083
|
R2
|
0.902
|
0.741
|
0.688
|
0.629
|
0.888
|
0.944
|
Pseudo-2nd order
|
qe,exp(mg/L)
|
12.03
|
24.91
|
48.94
|
73.49
|
90.90
|
110.81
|
qe,cal(mg/L)
|
12.56
|
25.06
|
48.19
|
72.20
|
92.17
|
123.03
|
k2(min− 1)
|
0.4240
|
0.6130
|
0.2210
|
0.2014
|
0.0370
|
0.0161
|
R2
|
0.999
|
0.999
|
0.999
|
0.999
|
0.998
|
0.998
|
IPD 1 Stage
|
ki
|
1.86
|
3.64
|
10.07
|
13.97
|
17.00
|
63.95
|
C
|
20.77
|
6.84
|
33.72
|
53.03
|
54.24
|
16.26
|
R2
|
0.998
|
0.976
|
0.872
|
0.911
|
0.768
|
0.992
|
IPD 2 Stage
|
ki
|
0.85
|
1.16
|
3.53
|
3.26
|
14.08
|
25.81
|
C
|
22.95
|
9.76
|
40.59
|
64.94
|
57.40
|
60.64
|
R2
|
0.715
|
0.939
|
0.915
|
0.912
|
0.797
|
0.649
|
IPD 3 Stage
|
ki
|
0.72
|
1.02
|
1.36
|
2.61
|
13.33
|
8.17
|
C
|
23.01
|
9.89
|
44.31
|
65.49
|
58.29
|
93.19
|
R2
|
0.986
|
0.882
|
0.984
|
0.878
|
0.941
|
0.999
|
3.4 Adsorption Isotherm of Adsorbent in MB
The adsorption isotherm is to study the relationship between the concentration of adsorbate in the liquid phase and the adsorbent at a specific temperature at equilibrium. Modeling the equilibrium adsorption data can help us study the adsorption mechanism, the interaction between the adsorbate and the adsorbent, the maximum adsorption capacity, and so on (Wang and Guo 2020b). In this study, several models (Langmuir, Freundlich and Temkin) were selected to analyze the adsorption mechanism (Table S2).
It can be seen from Fig. 4 and Table 2 that the Freundlich model (R2 = 0.954, 0.991, 0.979) is more suitable than Langmuir (R2 = 0.983, 0.863, 0.771) and Temkin model (R2 = 0.935, 0.871, 0.919) for describing the adsorption of MB molecules on cellulose/PVA/EG porous materials. It reveals the multilayer adsorption of MB on heterogeneous surfaces. Multilayer adsorption is due to mutual attraction between molecules, and additional molecules are superimposed on the first adsorption layer to form multilayer adsorption, which is essentially physical adsorption. In addition, \(n\) represents the adsorption strength of the adsorbent. When \(n>1\), it indicates that the adsorption behavior is favorable (Eltaweil et al. 2020a). At the same time, porous adsorbent has a greater adsorption strength at 25 ℃ as the lowest temperature in the test, indicating that high temperature is unfavorable for the adsorption behavior, which further supports its physical adsorption behavior.
Table 2
Key parameters of adsorption isotherm model
|
Langmuir
|
Freundlich
|
Temkin
|
\({q}_{m}\) (mg/g)
|
\(b\) (L/mg)
|
R2
|
\({k}_{f}\) (L/mg)
|
\(n\)
|
R2
|
\({B}_{1}\) (Jmol− 1)
|
\(A\) (L/mg)
|
R2
|
25 ℃
|
112.38
|
0.5161
|
0.983
|
43.1121
|
3.516
|
0.954
|
16.12
|
21.4108
|
0.935
|
50℃
|
30.30
|
0.0113
|
0.863
|
13.2766
|
2.1404
|
0.991
|
18.06
|
1.5371
|
0.871
|
75℃
|
70.72
|
0.0048
|
0.771
|
9.8355
|
2.4114
|
0.979
|
13.19
|
1.0879
|
0.919
|
3.5 Thermodynamic of Adsorption
Adsorption thermodynamic is used to survey the trend, degree, and driving force of the adsorption process and plays an important role in explaining the characteristics, laws, and mechanism of adsorption. The changes of parameters of thermodynamic of the adsorption process are generally calculated by equations of Gibbs equation and Vant’ Hoffs (Lakouraj et al. 2015):
$$\varDelta G=-RTln{K}_{c}$$
3
$${K}_{c}=\frac{{q}_{e}}{{C}_{e}}$$
4
$${lnK}_{c}=\frac{\varDelta S}{R}-\frac{\varDelta H}{T}$$
5
Where \(\varDelta G\) (kJ/mol) is Gibbs Free Energy, R and T are gas constant (8.314 J/mol/K) and absolute temperature (K), \({K}_{c}\) (L/mol) is the equilibrium constant of adsorption, \(\varDelta S\) and \(\varDelta H\) are entropy and enthalpy, respectively.
As shown in Table 3, the free energy \(\varDelta G\) is negative, indicating that the porous material adsorbs MB as a spontaneous process. As the temperature increases, the absolute value of \(\varDelta G\) decreases, indicating that the spontaneous tendency of adsorption decreases, which is not conducive to adsorption. The negative value of the adsorption enthalpy \(\varDelta H\) reveals that the adsorption process is an exothermic reaction. The negative value of the adsorption entropy \(\varDelta S\) indicates that the adsorption reduces the degree of freedom of the adsorbate molecules.
Table 3
Parameters of the thermodynamic model
|
\(\varDelta G\) (kJ/mol)
|
\(\varDelta H\) (kJ/mol)
|
\(\varDelta S\) (kJ/mol/K)
|
25 ℃
|
-7.9121
|
-55.7537
|
-0.1622
|
50℃
|
-2.0405
|
|
|
75℃
|
-0.0663
|
|
|
3.6 Reusability of the adsorbent
For industrial applications, the regeneration and recycling of adsorbents are very important. In this study, the adsorbent after adsorption was immersed in an ethanol solution to regenerate it. The experimental result shows that after five cycles, the MB dye removal rate of the composite material drops from 97.98–71.93%, possibly for a combination of chemistry and physics adsorption of cellulose/PVA/EG foam (Fig. 5). As is known to all, chemical adsorption is irreversible, but physical adsorption is reversible. Therefore, the adsorption process was mainly led by physical adsorption in the subsequent cycle. To sum up, the removal rate of adsorbents after regeneration remains high, which is appropriate for the treatment of MB wastewater.
3.7 Adsorption mechanism of Cellulose/PVA/EG Foam in MB
The adsorption mechanism of the adsorbate on the adsorbent is affected by many factors, such as the size and structure of the adsorbent. Based on the above experimental results, the adsorption mechanism scheme lists in Fig. 6. According to analyzing adsorption kinetics and isotherm, the adsorption process is mainly physical adsorption affected by the number of active sites and adsorbent structure. So, the possible interactions that determine the effect of MB adsorption are electrostatic adsorption, hydrogen bonding, and van der Waals force.
3.8 Comparison to other adsorbents of Adsorption Performance in MB
As shown in Table 4, the adsorbent has significant advantages. The main components of the materials are cellulose, PVA, and EG. These materials are low-cost, renewable, easy to recycle, and degraded. Although the materials developed by Boukhalfa, N., and Eltaweil, A.S. et al. have high adsorption capacity (Boukhalfa et al. 2019; Eltaweil et al. 2020b), the cost of their adsorption materials, such as carbon nanotubes or graphene oxide, is very high.
Table 4
Adsorption performance of different adsorbents on MB molecules
Adsorbent
|
Concentration of MB
|
Removal percentage (%)
|
qe (mg/g)
|
Ref.
|
Cellulose based porous foam composites
|
250 mg/L
|
90.9
|
110.8
|
This research
|
Maghemite/alginate/functionalized multiwalled carbon nanotubes
|
230 mg/L
|
-
|
905.5
|
(Boukhalfa et al. 2019)
|
CMC/GOCOOH composite microbeads
|
250 mg/L
|
72.09
|
180.32
|
(Eltaweil et al. 2020b)
|
CMC-Alg/GO hydrogel beads
|
15 mg/L
|
96.2
|
78.5
|
(Allouss et al. 2019)
|
Salecan-g-PAI
|
500 mg/L
|
22.1
|
107.1
|
(Qi et al. 2018)
|
Magnetic chitosan/clay beads
|
100 mg/L
|
-
|
82
|
(Bée et al. 2017)
|
activated carbon/cellulose biocomposite films
|
100 mg/L
|
-
|
103.66
|
(Somsesta et al. 2020)
|