2.2. Kinetic Analysis
To clarify the adsorption kinetic process and the involved mechanisms, the experimental data underwent analysis using various kinetic models: pseudo first order, pseudo-second order and intraparticle diffusion (see
Table S1 Supplementary Materials). The determination coefficient obtained in the adjustment (R
2) was the key to understand the mechanism underlying the adsorption process.
For furosemide, the R2 values obtained for pseudo first order model are in a range of 0.987 to 0.834 for β-CDs-EPI. In the case of hydrochlorothiazide, the R2 values are in a similar range (0.991 to 0.805).
The values of experimental
qe and calculated
qe were different, indicating that the adsorption process does not fit the pseudo first order model (
Figure 2A and
Figure 2B), these results are in agreement with those published for other PhACs with CDs polymers such as the case of ciprofloxacin and different organic pollutants where the pseudo first order model did not show a suitable adjustment to the experimental data [
27,
41,
42]; therefore, the pseudo second order model was applied.
The results of the graphical representation (
Figure 2C and 2D) of t/qt versus contact time show straight line for the two PhACs studied and the polymer used. The fit values (R
2), all higher than 0.99, reveal that the experimental data perfectly fit this kinetic model. Furthermore, the values of experimental
qe and
qe calculated are similar (see
Table S1 Supplementary Materials). The adjustment to the PSOM model indicates that the adsorption process is mainly governed by chemisorption. These results are in agreement with those published for other pharmaceuticals such as ciprofloxacin and different dyes [
27].
Different stages are involved in the adsorption process, such as the transport of pharmaceuticals molecules to the surface of the adsorbent material and the subsequent diffusion of PhACs molecules into the polymer. In addition, it was studied whether the intraparticle diffusion is the process that determines the rate of adsorption.
From the graphical representation of the model (
Figure 2E), two zones can be distinguished for furosemide, the first one more curved and later stable zone corresponding to an intraparticle diffusion. On the other hand, for hydrochlorothiazide, a linear rise was observed at the beginning and later a plateau (
Figure 2F), especially at the highest concentrations (15 and 20 mg/L).
The plateau stage means that intraparticle diffusion begins to decrease because there is not enough pollutant available. As could see in
Table S1 (see
Supplementary Materials), the value of
ki increases with the concentration. These result show that intraparticle diffusion is involved in the adsorption processes, but that it is not the only mechanism implied, there being other mechanisms such as the formation of inclusion complexes, adsorption on the external surface of the polymer, ionic exchange and the diffusion into the polymeric network [
41]. The amphiphilic cross-linked EPI-β-CD polymer exhibits both hydrophobic and hydrophilic cavities with a chaotic nature, resulting in a random distribution of shapes and sizes. This aspect is crucial for interpreting the obtained results. The hydrophilic properties of the polymer facilitate interaction with water, enhancing hydration through the potential formation of hydrogen bonds. Additionally, the polymer's hydrophilicity promotes better swelling of the network, thereby increasing its potential for diffusion. Furthermore, the hydrophobic cavities of the CDs and the cage structures formed by the cross-linker, abundant in ethylene oxide groups, attract hydrophobic PhACs [
37].
2.3. Adsorption equilibrium
The distribution in the equilibrium of the PhACs between the polymer and the solution was studied, it was carried out applying the Freundlich, Langmuir and Tempkin isotherm models [
43].
Table 1 shows the parametric values obtained for each model.
For the Freundlich isotherm, which describes the adsorption processes on heterogeneous surfaces, a linear representation was obtained for the two diuretics studied (
Figure 3), reaching K
F constants values of 0.044 (L/g) for furosemide and 0.029 (L/g) for hydrochlorothiazide. The order of K
F values suggests that furosemide exhibits the strongest adsorption and the highest capacity at the β-CDs-EPI polymer surface, followed by hydrochlorothiazide. In addition, the determination coefficients (R
2) of 0.991 and 0.905 for furosemide and hydrochlorothiazide, respectively (
Table 1), indicated that the Freundlich equation fitted the adsorption data better than the Langmuir and Temkin models.
The magnitude of the Freundlich exponent n
F for hydrochlorothiazide and furosemide, that ranged from 0.737 to 0.817 respectively, indicates that that the sorption mechanism is controlled by adsorption and not absorption to the Freundlich model, giving a rational description of the experimental data, involving both chemisorption and physisorption processes, agreeing with other published studies [
27,
44]. The role of the Freundlich isotherm in adsorption lies in its ability to describe multilayer adsorption onto heterogeneous surfaces. It is particularly useful when the adsorption process does not follow a straightforward monolayer adsorption (as assumed in models like Langmuir isotherm) and when multiple layers of solute can be adsorbed onto the surface with varying energies. Moreover, the exponent highlights the variety of energies linked to adsorption of both diuretics on the β-CDs-EPI polymer surface. Moreover, n < 1 for furosemide and hydrochlorothiazide indicates that upon increasing the PhACs concentration/loading the binding energy between the surfaces and both compounds is reduced.
Table 1.
Adsorption isotherm coefficients obtained for β-CDs-EPI polymer by the Freundlich, Langmuir and Tempkin models.
Table 1.
Adsorption isotherm coefficients obtained for β-CDs-EPI polymer by the Freundlich, Langmuir and Tempkin models.
Isotherm |
Parameter |
Furosemide |
Hydrochlorothiazide |
Freundlich |
KF (L/g) |
0.044 |
0.029 |
nF |
0.817 |
0.737 |
R2
|
0.991 |
0.905 |
Langmuir |
qmax(mg/g) |
1.282 |
0.844 |
KL
|
0.050 |
0.038 |
aL
|
0.039 |
0.045 |
∆G |
-16919.810 |
-16730.651 |
R2
|
0.516 |
0.514 |
RL
|
0.838-0.564 |
0.817-0.527 |
Tempkin |
aT
|
0.525 |
0.448 |
bT (kJ/mol) |
6.890 |
6.79 |
R2
|
0.943 |
0.872 |
On the other hand, the Langmuir model proved suitable for scenarios where materials exhibited regularly energetic adsorption sites and monolayer adsorbate coverage. The isotherm postulated that all sites showed uniform surface coverage [
45]. Unlike the Freundlich isotherm, the values of determination coefficients (R
2) (
Table 1) for Langmuir model were lower in all cases, which shows that the data does not fit well to this isotherm.
The most important parameter obtained with the Langmuir model is the maximum monolayer adsorption capacity (qmax), under the studied conditions, being 1.282 mg/g for furosemide and 0.844 mg/g for hydrochlorothiazide, showing that that the order of saturation of the β-CDs-EPI polymer surface with diuretic per mg diuretic was 1.45 times higher for furosemide. This trend fits well with KF values previously described by Freundlich model.
In addition, with the Langmuir isotherm it was analyzed whether the adsorption process is favorable or not, with the separation factor (RL). The adsorption process is considered favorable when it is between 0 and 1. In the case, of furosemide and hydrochlorothiazide, at the concentrations studies, the adsorption process on β-CDs-EPI polymer is favorable.
The Langmuir parameter KL increased in the order 0.038 L/g (hydrochlorothiazide) < 0.050 L/g (furosemide). This constant is mainly related to the adsorption energy and could give information about the β-CDs-EPI polymer–PhACs interaction and binding process strengthens.
Thus, based on Langmuir model, one expects that the β-CDs-EPI polymer–PhACs interaction increases in the order hydrochlorothiazide < furosemide. The same order but with different values 0.448 L/g (hydrochlorothiazide) < 0.525 L/g (furosemide), was obtained from Temkin binding constant aT which is also related to the binding strength.
Lastly, the equilibrium experimental data were analyzed by fitting the results to the Tempkin isotherm (
Figure 3). It assumes the linear decrease in the heat of adsorption of all the molecules found in the external layer, as a consequence of the interactions that occur between the polymer and the PhACs, in this case the binding energies are uniformly distributed [
46].
For ionic exchange, the binding energies are between 8 and 16 kJ/mol and for physisorption they are between -40 kJ/mol. The Temkin b
T values for furosemide (6.890 kJ/mol) and hydrochlorothiazide (6.79 kJ/mol), suggests that the heat of adsorption on the β-CDs-EPI polymer increases in the order hydrochlorothiazide < furosemide as well as that physical and chemical processes are involved in the adsorption. Regardless the Temkin constant b
T, all other parameters from the represented isotherm models in
Figure 3, refer to stronger adsorption and higher capacity for the furosemide diuretic compared to hydrochlorothiazide by considering the number of PhAC mg.
Figure 3.
Isotherm analysis (A,B) Freundlich model; (C,D) Langmuir model; (E,F) Tempkin model for Furosemide (A, C, E) and Hydrochlorothiazide (B, D, F) by β-CDs-EPI polymer. N=3.
Figure 3.
Isotherm analysis (A,B) Freundlich model; (C,D) Langmuir model; (E,F) Tempkin model for Furosemide (A, C, E) and Hydrochlorothiazide (B, D, F) by β-CDs-EPI polymer. N=3.
2.6. Design continuous flow Prototype Adsorption
Once the batch adsorption data of the β-CDs-EPI polymer, using two PhACs as pollutants model, has been determined in distilled water, and the adsorption mechanism established by applying kinetic (pseudo-first order, pseudo-second order and intraparticle diffusion) and equilibrium (Freundlich, Langmuir and Tempkin isotherm) models, where are involved both physical and chemical processes, a laboratory-scale prototype was design to validate in subsequent tests, its behavior in real scenarios working continuously.
For that, a stepwise approach was followed starting with,
i) the characteristics of the pollutant solution and the adsorbent (synthesized as described in
Section 3.2). In our case we used to ascertain the adsorbent performances in batch assays of β-CDs-EPI polymer (
Table 2), working solutions of furosemide and hydrochlorothiazide in distilled water at concentrations between 5 and 20 mg/L;
ii) the establishment of basic processes of adsorption/desorption mechanisms, including work ranges and ratios. As in any engineering process, these two points are key to obtaining good execution and results of the engineering project.
In this sense, we confront a substantial challenge since we must design and build equipment that works in continuous flow in a way that allows us to confirm the results obtained in the laboratory, and also provide sufficient information to validate in a future the operation of the process on an industrial scale.
Once these two points are stated, we must define the adsorption and desorption processes steps, as well as the equipment necessary for their accomplishment. To do this, as a starting point we must prepare a summary (
Table 3) with the recommended design parameters for a standard adsorption system.
As a result of these input parameters, the next step is iii) the definition of the stages of the process, together with the elaboration of a flowchart that helps us identify and understand each of the stages and equipment involved.
Table 2.
Values of physical-chemical parameters obtained for β-CDs-EPI polymer for Furosemide and Hydrochlorothiazide by in batch assay.
Table 2.
Values of physical-chemical parameters obtained for β-CDs-EPI polymer for Furosemide and Hydrochlorothiazide by in batch assay.
Parameter |
Value |
qmax (furosemide) (mg/g) |
1.282 |
qmax (hydrochlorothiazide) (mg/g) |
0.844 |
Density (g/cm³) |
1.06 |
Swelling |
4 ± 1 |
Particle size (mm) |
0.1→0.3 |
Stability range (pH) |
2→11 |
Temperature range (ºC) |
5→35 |
Solubilty in H2O |
Insoluble |
Table 3.
Values of recommended design parameters for a standard adsorption system.
Table 3.
Values of recommended design parameters for a standard adsorption system.
Parameter |
Value |
Adsorbent Volume (L) |
1→3 |
Column diameter |
To define |
Adsorbent bed depth (mm) |
150→550 |
Adsorbent expansion (%) |
Up to 100 |
Contact time (min) |
1→7.5 |
Loading flow rate (BV/h) |
8→40 |
Desorbent flow rate (BV/h) |
2→5 |
Desorbent contact time (min) |
20→60 |
Desorbent displacement (BV of water) |
2→4 |
Final rinse (BV service flow rate) |
2→10 |
Thus, to get started with the design we selected an acceptable adsorbent volume for continuous work in the range of 1 to 3 L. After that, we carry out the necessary calculations to define the necessary column for the piloting. The selection of the column is the point of outmost importance in the design, conditioning its configuration the design of the rest of the prototype components.
For that, adsorbents volumes of 1 L, 2 L and 3 L were used to determine the necessary height and diameter of the column for the continuous process setting for each adsorbent filling volume, a loading flow rates relation of 1:8 and 1:40 (i.e: 1 L adsorbent volume:8 L/h flow; 1 L adsorbent volume:40 L/h flow; 2 L adsorbent volume:16 L/h flow; 2 L adsorbent volume:80 L/h flow; 3 L adsorbent volume:24 L/h flow; 2 L adsorbent volume:120 L/h flow) and contact times of 7.5 min for 1:8 and 1.5 min for 1:40 rates. A column of Ø 90 mm capable of accommodating the β-CDs-EPI adsorbent polymer volumes chosen, with a final adsorbent bed depth do not exceeding 550 mm nor 1,200 mm for backwash expansion column height (the maximum recommended height), was selected.
Transparent polyvinyl chloride with a design pressure of 6 bar, as construction material of the column was selected. This material has chemical compatibility with the PhACs and desorbing (220 mM acetate buffer pH 4.0) solutions proposed for the pilot tests, while allowing us to visualize the behavior of the process at different adsorption/desorption stages.
The columns shall be equipped with two polypropylene nozzles at the top and bottom that will retain the adsorbent within the column. The light passage is set at 100 μm, measured below the particle size of the adsorbent used. The loading and emptying operations of the adsorbent will be carried out through a removable accessory located at both ends of the column.
Taking into account all the design criteria established above, we proceeded to determine the parametric values of the prototype for a column of 90 mm diameter (Ø) and one of 63 mm diameter (
Table 4).
Taking into account the stated criteria 1,200 mm for the maximum recommended height in expansion backwash process and 550 mm of bed depth column, we select Ø90 mm columns because the height for a smaller diameter column Ø63 mm, exceeding the necessary heights in 350 mm for 2 L of adsorbent volume and flow rates of 16 L/h or 80 L/h, and 1,130 mm for 3 L of adsorbent volume and flow rates of 24 L/h or 120 L/h.
Once the column was selected we proceeded with the definition of process stages, iv) adsorption step. As defined in the design parameters, the flow rates established for the adsorption process are in the range from 8 to 120 L/h, ensuring that the contact time between the emergent pollutant solution and the adsorbent β-CDs-EPI polymer will be within the values established in the stated design parameters (1 to 7.5 min).
Table 4.
Column size design calculations of the prototype.
Table 4.
Column size design calculations of the prototype.
|
Column size design calculations |
Parameters |
Ø90 mm |
Ø63 mm |
Flow (L/h) |
8 |
40 |
8 |
40 |
Flow rate (m/h) |
1.43 |
7.17 |
3.10 |
15.51 |
Ad volume (L) |
1 |
1 |
1 |
1 |
BV (BV/h) |
8 |
40 |
8 |
40 |
Area (m²) |
0.0056 |
0.0056 |
0.0026 |
0.0026 |
Bed depth (m) |
0.18 |
0.18 |
0.39 |
0.39 |
Expansion (%) |
100 |
100 |
100 |
100 |
Column height (m) |
0.36 |
0.36 |
0.78 |
0.78 |
Contact time (min) |
7.5 |
1.5 |
7.5 |
1.5 |
Flow (L/h) |
16 |
80 |
16 |
80 |
Flow rate (m/h) |
2.87 |
14.33 |
6.20 |
31.02 |
Ad volume (L) |
2 |
2 |
2 |
2 |
BV (BV/h) |
8 |
40 |
8 |
40 |
Area (m²) |
0.0056 |
0.0056 |
0.0026 |
0.0026 |
Bed depth (m) |
0.36 |
0.36 |
0.78 |
0.78 |
Expansion (%) |
100 |
100 |
100 |
100 |
Column height (m) |
0.72 |
0.72 |
1.55 |
1.55 |
Contact time (min) |
7.5 |
1.5 |
7.5 |
1.5 |
Flow (L/h) |
24 |
120 |
24 |
120 |
Flow rate (m/h) |
4.30 |
21.50 |
9.31 |
46.54 |
Ad volume (L) |
3 |
3 |
3 |
3 |
BV (BV/h) |
8 |
40 |
8 |
40 |
Area (m²) |
0.0056 |
0.0056 |
0.0026 |
0.0026 |
Bed depth (m) |
0.54 |
0.54 |
1.16 |
1.16 |
Expansion (%) |
100 |
100 |
100 |
100 |
Column height (m) |
1.07 |
1.07 |
2.33 |
2.33 |
Contact time (min) |
7.5 |
1.5 |
7.5 |
1.5 |
The adsorption stage begins with the preparation of the PhACs solution in a 50 L tank TK-01 (as described in
Section 3.9). This solution is sent to the adsorption column through a self-priming diaphragm pump (P-01), responsible for propelling the solution towards the column at the determined flow rate and pressure (2 to 4 bar). This flow rate will be measured by an in-line flow meter FI-01 installed in the feed pipe, while the pressure is measured in the inlet and pressure gauges output of column at PI-01 and PI-02, respectively. To regulate the flow, a needle valve RG-01will be used. At the exit of the adsorption column, the treated PhACs solution will be collected in the product tank TK-02 (50 L).
Next, it is necessary to estimate the
v) adsorption cycle. Based on the capacity data obtained from the adsorption isotherms in batch, and once stated the volumes to be used in the columns, we carried out the necessary calculations to estimate the depletion of the β-CDs-EPI polymer used as adsorbent. These theoretical values are based on the data obtained for in batch absorption processes with continuous stirring, but they provide us with the necessary information for a preliminary calculation of the pilot plant design (
Table 5).
Table 5.
Theoretical values for flow adsorption cycle estimation of β-CDs-EPI polymer for Furosemide and Hydrochlorothiazide based in batch assay data.
Table 5.
Theoretical values for flow adsorption cycle estimation of β-CDs-EPI polymer for Furosemide and Hydrochlorothiazide based in batch assay data.
Parameter |
Value |
PhACs concentration (mg/L) |
5→20 |
Tank volume of PhACs solution (L) |
50 |
Amount PhACs concentration (mg) |
250→1000 |
β-CDs-EPI qmax (furosemide) (mg/g) |
1.282 |
β-CDs-EPI qmax (hydrochlorothiazide) (mg/g) |
0.844 |
β-CDs-EPI volume (L) |
1→3 |
β-CDs-EPI weight (g/column) |
1,060→3,180 |
Amount β-CDs-EPI qmax (furosemide) (mg) |
1,358→3,846 |
Amount β-CDs-EPI qmax (hydrochlorothiazide) (mg) |
894→2,683 |
Now, we must propose in the design the two possible flow operation scenarios, co-current and counter-current (as described in
Section 3.9).
To operate co-current, the flow direction is downward. The PhACs solution enters through the upper part of the column, while the effluent exits through the lower part of the column through a polypropylene nozzle with a passage light of 100 µm.
On the other hand, to operate counter-current, the flow direction is upward. The PhACs solution exits through the upper part of the column through a polypropylene nozzle (passage light of 100 µm) located at the upper part of the column.
Table S2 and S3 of
Supplementary Materials details the on/off positions of the valves when working in co-current (
Table S2) and counter-current (
Table S3) way, respectively.
As the prototype is designed to function continuously, the adsorbent polymer must undergo proper regeneration upon reaching its capacity limit by a backwash process.
This involves introduction of a desorbing solution to sponge the adsorbent restoring the capacity of the β-CDs-EPI polymer, ensuring complete contact of the solution with the adsorbent when the operating flow rate is between 55 and 110 L/h. The on/off positions of the valves to carry out the backwash process is the same as described above when working in counter-current flow (see
Table S3 Supplementary Materials).
The
vi) desorption step starts when the β-CDs-EPI polymer lost their adsorbent capacity passing through a desorbing solution (as described in
Section 3.9), to restore its adsorption capacity for the removal of PhACs. This solution is introduced in a 50 L tank TK-03, and the valves HV are placed in the desorption position (see
Table S4 Supplementary Materialswhen working in co-current flow)., turning on the pump P02. The desorbing solution flow rate (2 to 5 BV/h) is regulated with the RG-02 valve and the FI-01 flowmeter, while the pressure is measured in the inlet and pressure gauges output of column at PI-01 and PI-02, respectively. The contact time between the desorbing solution and the adsorbent β-CDs-EPI polymer will be within 20 to 60 min with a displacement of desorbing solution between 2 to 4 BV of water. At the exit of the desorption column, the desorbing solution will be returned to the tank TK-04 (50 L).
Finally, a rinsing process is carried out in order to remove possible traces of desorbing solution that remain in the column before starting a new adsorption process. The volume of effluent used will be between 2 and 10 times higher than the volume of β-CDs-EPI adsorbent polymer filled in the column. The on/off positions of the valves to carry out the rinsing process is the same as described above when working in co-current flow (see
Table S2 Supplementary Materials).
Taking into account the theoretical calculations obtained using previous inputs, a pilot-plant laboratory scale prototype (PPLSP) was built to validate under continuous flows in subsequent tests (
Figure 4), if fulfilled the design criteria regarding the parameters obtained in batch adsorption process.
Figure 4.
Image of the designed pilot-plant laboratory scale prototype (PPLSP).
Figure 4.
Image of the designed pilot-plant laboratory scale prototype (PPLSP).
As can be seen in
Figure 4, we decided to include, in addition to the selected one (Ø90 mm), an additional smaller diameter column Ø63 mm not exceeding the necessary 1,200 mm height (but lower bed depth working), in order to obtain comparative data of the adsorption process in continuous way, as a function of the column size.
Theory results of pharmaceutic retention using this prototype suggest that this simple and inexpensive technological setup could be scaled up to a functional field application to effectively capture emerging pollutants. To confirm this theoretical postulate, the built prototype should be subjected to representative continuous assays with water and wastewater samples enriched with emerging contaminants and the obtained results will be displayed in a succeeding work.