Kinetic Modeling of the UV/H2O2 Process: Determining the Effective Hydroxyl Radical Concentration

A kinetic model for pollutant degradation by the UV/H2O2 system was developed. The model includes the background matrix effect, the reaction intermediate action, and the pH change during time. It was validated for water containing phenol and three different ways of calculating HO° level time-evolution were assumed (non-pseudo-steady, pseudo-steady and simplified pseudo-steady state; denoted as kinetic models A, B and C, respectively). It was found that the kind of assumption considered was not significant for phenol degradation. On the other hand, taking into account the high levels of HO2° formed in the reaction solution compared to HO° concentration (~10–7 M >>>> ~10–14 M), HO2° action in transforming phenol was considered. For this purpose, phenol-HO2° reaction rate constant was calculated and estimated to be 1.6x103 M-1 s-1, resulting in the range of data reported from literature. It was observed that, although including HO2° action allowed slightly improving the kinetic model degree of fit, HO° developed the major role in phenol conversion, due to their high oxidation potential. In this sense, an effective level of HO° can be determined in order to be maintained throughout the UV/ H2O2 system reaction time for achieving an efficient pollutant degradation.


Introduction
Nowadays, one of the major problems associated with the presence of toxic and persistent pollutants in the aquatic environment is the unfeasibility of conventional treatments for the effective removal of those substances [1][2][3].Hence the application of alternative technologies, such as advanced oxidation processes, is needed [4].Among these techniques, the UV/H 2 O 2 system is included.It consists of the photolysis of hydrogen peroxide (H 2 O 2 ) by applying ultraviolet (UV) radiation resulting in the generation of hydroxyl radicals (HO°) [5,6].This process may be performed at room temperature and pressure, it has no mass transfer problems, it is easy to maintain and operate, no sludge requiring a subsequent treatment and disposal is produced, and it may achieve a complete pollutant mineralization [5].Therefore, the UV/H 2 O 2 system seems to be a promising alternative for the treatment of water containing toxic and recalcitrant substances.However, this kind of technology can be expensive due to the associated electrical and oxidant costs [7].
In order to reduce costs and make the process more feasible for industrial applications, the UV/H 2 O 2 system optimization is required [8] and kinetic models can be considered as functional tools for this purpose.Up to date, several kinetic models have been proposed for describing the UV/H 2 O 2 process and predicting different pollutant removal rates [6][7][8][9][10][11][12][13][14][15].In some of these models [13,15], the proposed set of ordinary differential equations (ODE) defining the studied pollutant degradation rate can be simplified into a pseudo-first-order kinetic expression, whose solution is an exponential one.In such as models, experimental results are fit to that solution.Subsequently, model predictions agree well with laboratory data.In that kind of models the calculated reaction rate constants for the tested pollutant degradation are apparent reaction rate constants (K app ), which include pollutant removal reaction rate constants and the values of parameters such as the quantum yield for the oxidant, the conjugate base (HO 2 − ), and the contaminant photolysis, the initial level of the chemical species involved in pollutant oxidation, the optical path length of the system, the UV-light intensity, and the molar extinction coefficients of H 2 O 2 /HO 2 − and pollutant, among others.Therefore, knowing those parameters is not required.
On the other hand, there are dynamic kinetic models that try to solve the considered ODE set, for which the values of the mentioned variables are required, increasing the complexity of the kinetic model.In order to solve the proposed ODE set, the pseudo-steady state approximation assumption for reactive intermediates, such as HO°, is invoked by arguing that these chemical species are as transient ones as their concentration can be presumed to be at a pseudo-steady state [10,12,14].In other models [6][7][8][9][10][11][12][13][14][15], on the contrary, the non-pseudo-steady state premise in free radical rate expressions for predicting the degradation of the probe compound in a more accurate way is applied.However, although hydroperoxyl radicals (HO 2 °) are involved in those models, none of them, excluding Huang and Shu [11] and Liao and Gurol [12] models, includes HO 2 ° in the target pollutant oxidation final expression.Additionally, some of these models cannot be reproduced unless a conversion factor is included, as demonstrated by Audenaert et al. [8].
In this sense, the aim of this work was to develop a kinetic model based on the main chemical and photochemical reactions for pollutant degradation in water systems by the UV/H 2 O 2 process taking into account the decomposition of the pollutant through direct photolysis, HO° oxidation, and HO 2 ° and superoxide radical (O 2 °−) transformation.Furthermore, HO° scavenging effects of carbonate (CO 3 2− ), bicarbonate (HCO 3 − ), sulfate (SO 4 2− ), and chloride (Cl − ) ions were considered.pH changes in the bulk and the detrimen-Physico-Chemical Wastewater Treatment and Resource Recovery tal action of the organic matter (OM) and the reaction intermediates in shielding UV and quenching HO° were studied.The influence of the pseudo-steady and non-pseudo-steady state hypothesis for determining HO° concentration-evolution with time was also examined and the second-order HO 2 ° reaction rate constant for the studied pollutant was determined.MATLAB software was used to solve the ODE set that characterizes the current model and the results were validated by using experimental data obtained from the literature for phenol (PHE) degradation by the UV/H 2 O 2 process in a completely mixed batch photoreactor., the parent compound and the DOM in terms of dissolved organic carbon (DOC), as well as the formation and disappearance of intermediate products, also considered as DOC.Moreover, it accounts for the solution pH change due to the mineralization of organic compounds and the formation of acids.
In addition to the oxidant photolysis, the target pollutant (C) may interact with the UVradiation, undergoing degradation and producing reaction intermediates.A fraction of those by-products can be dissolved in the solution [16].This fraction is denoted as DOC (Eq.( 10)) [6].The reaction rate for the contaminant direct photolysis is obtained through Eq. (11).
Physico-Chemical Wastewater Treatment and Resource Recovery where c is the quantum yield of pollutant photolysis and , C is the amount of UV-light absorbed by the contaminant, calculated through Eqs. ( 12) and ( 13).
( ) Once HO° are produced, they rapidly react with the pollutant of interest, degrading it to form reaction intermediates (Eq.( 14)), which subsequently can be attacked by HO° and undergo further degradation to produce final products, such as CO 2 , H 2 O, and mineral acids (Eq.( 15)) [7].In this model, intermediate substances were considered as HO° scavengers as well as UVlight absorbers.Additionally, it was assumed that the pH of the solution decreased due to the conversion of the target pollutant, and consequently the reaction intermediates, into carbon dioxide (i.e., H 2 CO 3 * in the aqueous phase); although it must be highlighted that not all the DOC is mineralized, since carboxylic acids are also formed during the oxidation process, making the pH of the bulk decreases as well [8].Under this presumption, Eq. ( 15) is simplified as Eq. ( 16).The mass balances for the evolution of the pollutant, the dissolved organic fraction of the formed by-products, HO° and the H 2 CO 3 * concentrations with time are shown by Eqs. ( 17)-( 20), respectively.14) and ( 15) or ( 16).
In the UV/H 2 O 2 system, recombination of HO° can occur to produce H 2 O 2 .However, these free radicals can also react with H 2 O 2 and HO 2 − , particularly when the oxidant is in excess, to produce HO 2 °.Although HO 2 ° are less reactive than HO° (E° = 0.98 and 2.8 V, respectively) [8], they can also be involved in pollutant degradation, especially if these radicals are produced in high amounts in the system.Furthermore, HO 2 ° can produce O 2 °−, which subsequently can participate in pollutant degradation and mineralization [17].Therefore, the role of these reactive oxygen species was included in the proposed kinetic model.
It is important to note that as the oxidation process develops, the pH of the solution generally goes down, and consequently, some chemical species appear while other species vanish.In order to consider the change of chemical species inside the bulk according to the pH of the solution in the kinetic model, a correction factor (δ R i ) was introduced for the photolysis of H 2 O 2 , HO 2 − , and for the reaction between H 2 O 2 and HO°.This correction factor can adopt two values (0 and 1).When δ R i =1, reaction R i is promoted (i.e., the time-varying concentrations of the chemical species involved in reaction R i are taken into account in the model).When δ R i =0, reaction R i is not considered in the model and, subsequently, the chemical species taking part in reaction R i are neglected.
On the other hand, species commonly present in water, such as DOM and inorganic anions (e.g., CO 3 2− , HCO 3 − , SO 4 , and Cl − , among others) may also have a significant effect because of their ability to absorb UV-light and/or to scavenge HO°.The HO° scavenging effect of matrix constituents drastically limits the oxidation action of HO°, leading to a decrease in the performance of the system [18].
Taking into account all the mentioned processes and in order to give a more realistic view of what happens in the UV/H 2 O 2 process, the kinetic equation describing pollutant degradation can be expressed as Eq. ( 21).
where the terms , C , In order to quantitatively evaluate the contribution of the cited terms in the contaminant removal, parameters d, f, g, and h were introduced in Eq. ( 21), as described by Eq. ( 22).
When d = f = g = h = 0 (i.e., when the initial concentrations of oxidant radical species are equal to zero), the kinetic model only describes the degradation of the pollutant by photolysis.If d = h = 1 and f = g = 0, in addition to the photolysis conversion, in a deionized water, the model can predict pollutant transformation through HO° and CO 3 °− (the latter from the reaction between HO° and HCO 3 − or CO 3
Physico-Chemical Wastewater Treatment and Resource Recovery

Proposed kinetic model
In the developed kinetic model, three different ways of calculating the evolution of HO°c oncentration during time were presumed: (a) a non-pseudo-steady or transient state (i.e., the net formation rate of HO° is different from zero); (b) a pseudo-steady state; and (c) a simplified pseudo-steady state; correspondingly denoted as kinetic model A, B, and C.
In the prediction model A the concentration of HO° can be calculated with Eq. ( 23).From Eq. ( 23), the HO° concentration can be written as Eq. ( 24) (prediction model B).
Considering that the oxidant is in a high level (i.e., where Eq. ( 24) can be simplified to Eq. ( 25) (prediction model C).k i and [X i ] are the reaction rate constants between HO° and species i, and the concentration of species i, respectively.
The initial values of DOC and inorganic anionic species (CO 3 2− , HCO 3 − , SO 4 2− , Cl − , etc.) are set according to the conditions of the water to be treated.

Numerical solution of the proposed kinetic models
The ODE system compiled in Table 1 was solved applying MATLAB software and ODE15S function.For simultaneously solving the ODE set of the proposed kinetic models, it was necessary to define several photochemical parameters such as φ H 2 O 2 , φ pollutant , φ DOC , ε H 2 O 2 , ε HO 2 -, ε pollutant , ε DOC , I 0 and l.Additionally, the initial concentrations of all the species involved in the system and the rate constants of the chemical reactions between those species were required (Table 1).During setting up the model, differential rate equations describing the time dependence of the concentration of the variety of the considered species were defined and plotted.

Results and discussion
In order to validate the proposed kinetic models, experimental data for PHE degradation by the UV/H 2 O 2 process were used from Alnaizy and Akgerman [19] study.These authors conducted a set of experiments in a completely mixed batch cylindrical photoreactor made on Pyrex glass.The used photochemical parameters and the kinetic reaction rate constants of PHE with HO°, O 2 °−, and CO 3 °− are presented in Table 2. Table 2. Values of the parameters used in the kinetic model validation for phenol degradation.

Assumptions taken into consideration
As stated previously, the developed kinetic models A, B, and C employed the non-pseudosteady, the pseudo-steady, and the simplified pseudo-steady state assumption, respectively, to estimate HO° concentration.In the proposed models, the impact of UV radiation individually and/or the combined action of H 2 O 2 and UV light (including the effect of HO°, HO 2 °, O 2 °−, and CO 3 °−) on PHE degradation was studied.
A PHE concentration of 2.23 × 10 −3 M and a H 2 O 2 /PHE ratio of 495 were selected for validating the model.The used PHE solution was prepared by adding the appropriate amount of pollutant solution to deionized water [19].Therefore, the effect of inorganic anions, excluding HCO 3 − and CO 3 2− was not taken into account.In this sense, in the PHE degradation rate expression (ODE 1 ) the contribution of inorganic anion radicals, such as SO 4 °−, H 2 ClO°, HClO°, Cl°, and Cl 2 °−, was not studied.was assumed that the other terms included in ODE 1 contributed to pollutant degradation.
As the treated water was deionized, the presence of OM different from the parent compound in the initial solution was neglected ([DOC] 0 = 0 M).Hence, the DOC in the solution came from PHE photolysis and free radical (HO°, HO 2 °, O 2 °−, and CO 3 °−) oxidation.
On the other hand, several authors agree that OM reduction by direct photolysis in a UV/ H 2 O 2 oxidation process can be neglected [7][8][9].That is the reason why this was not included in the proposed kinetic model.However, it is highlighted that the OM is able to absorb UVlight, preventing UV-penetration into the bulk and avoiding H 2 O 2 /HO 2 − and pollutant direct photolysis.Therefore, the detrimental effect of UV-shielding in PHE degradation was taken into consideration.For including this effect in the kinetic model, OM molar extinction coefficient, referred as DOC molar extinction coefficient (ԑ DOC ), must be previously known.Although this parameter has already been measured [7,8], its value is not a universal one, since DOM is a complex group of aromatic and aliphatic hydrocarbon structures with attached functional groups [20,21].As Alnaizy and Akgerman [19] did not measure this variable, a mean value from Peuravuori and Pihlaja [22] study was presumed.This value corresponded to ԑ DOC(280 nm) = 35 967 M −1 m −1 , which was in the same order of magnitude than ɛ PHE .This number could be acceptable due to the formed aromatic intermediates during PHE conversion [19] conserve structural similarities with the parent compound, like the aromatic ring, responsible for the molecule excitation.
Additionally, it is widely known that the quantum yield of a compound is dependent on the excitation wavelength and the pH of the solution.For 254 nm, PHE quantum yield in an aqueous solution was found to be in the range of 0.02 -0.12 mol Ein −1 at pH 1.6 -3.2 [23].Therefore, an average value (0.07 mol Ein −1 ) was taken as PHE quantum yield.Moreover, it is widely recognized that the solution pH decreases as the process proceeds.This variation in the pH can cause difficulties in modeling studies, since the presence of radical species such as HO 2 °, and O 2 °−, among other chemical species involved in the oxidation system, is significantly dependent on the pH of the medium.Therefore, to give a more realistic view of what happens inside the reaction medium, H 2 O 2 and HO 2 − photolysis reactions (Eqs.( 1) and (2), correspondingly), as well as reactions expressed in Eqs. ( 26) and (27) were discriminated in the model according to the solution pH time evolution, as it is simplified in Figure 3.For selecting the suitable reactions with regard to the pH changes over time, previous information about the evolution of the pH during the performance of the process is required.However, in some occasions this is not provided.In this case, the initial pH of the solution was 6.8 [19].At this pH, one of the predominant species into the bulk was H 2 O 2 , since the pKa of the H 2 O 2 / HO 2 − equilibrium is 11.6.Therefore, the photolysis of HO 2 − was neglected in the model performance (i.e., δ R 3 and δ R 4 ).Additionally, the initial concentrations of HO° and other considered species were assumed to be zero, with the exception of the pollutant and H 2 O 2 .On the other hand, it was found that the pH of the medium rapidly dropped from 6.8 to 4.7-4.2(4.5 as a mean value) within the first 30 min of radiation [19].Therefore, and according to Figure 3, Eq. ( 26) was considered during the first 30 min of the reaction while Eq. ( 27), after that time (i.e., δ R 7 = 1, δ R 8 = 0 and δ R 7 = 0, δ R 8 = 1, before and after the first 30 min of the process, respectively).

Kinetic model validation
Initially, the associated ODE sets with model A, B, and C were solved.Parameters f = g = 0 and d = h = 1 were considered in order to solely investigate the influence of the direct photolysis, HO°, and CO 3 °− on PHE degradation.The measured and simulated results for PHE direct photolysis alone are also presented.It is observed that more than 90% of the initial PHE concentration was removed after 220 min due to both direct photolysis and indirect degradation (primarily due to HO° attack).The effect of CO 3 °− could be seen as marginal because of the reduced number of those radicals, in the range of 10 −15 M, and the low reaction rate constant with PHE compared to HO°.In addition, it is shown that PHE was not completely removed by direct UV photolysis under the tested photochemical conditions (Table 2), since approximately 50% of the total PHE degradation was attributed to UV photolysis, as it was experimentally determined by Alnaizy and Akgerman [19].On the other hand, the figure demonstrates that the prediction kinetic model C was in good agreement with the available experimental data with a relative high correlation factor (R 2 =99.34%).For prediction models A and B, R 2 were 97.41 and 97.37%, respectively.As a result of the subtle differences between kinetic models A and B, the depicted line describing model A overlaps model B line.An acceptable agreement between model predictions considering only the UV radiation and experimental data was also verified (R 2 = 98.25%).Furthermore, Figure 4 clearly presents that the removal rate of the target pollutant was not significantly dependent on the hypothesis assumed to estimate the HO° level (non-pseudosteady, pseudo-steady and simplified pseudo-steady state assumptions).This could be explained from the relatively low concentration of those reactive species in the solution (with a magnitude order of 10 −14 M) when compared to the level of other species involved in the system, such as HO 2 ° and O 2 °−, whose concentrations were in the range of 10 −7 M. The number of HO° remaining in the solution is in concordance with the low final HO° levels found in the literature [24,25] and even higher than those reported by Ray and Tarr [26].
In Figure 5, the evolution of the HO°, HO 2 °, O 2 °−, DOC and H 2 CO 3 * normalized estimated concentrations using the prediction model A is depicted.It is observed that HO 2 ° number increased as the oxidation system proceeded, while O 2 °− level decreased.Typically, the effect of HO 2 ° and O 2 °− radicals are neglected in the UV/H 2 O 2 system [6-8, 15, 17, 27, 28] since they are found to be less reactive than HO°, as stated previously.However, when they are produced in a high level, they could also participate in the contaminant oxidation.That is the current case of HO 2 °, as [HO 2 °] >>>> [HO°].Therefore, the contribution of HO 2 ° to PHE degradation should be studied.On the other hand, in this work the action of O 2 °-in PHE conversion can be omitted since O 2 °-level decreased with the reaction time, as expected, because of the drop in the pH solution.Furthermore, generally, there is a drop in the pH of the medium as the system progresses.This is probably due to acidic compound formation, such as carboxylic acids and H 2 CO 3 * resulting from pollutant degradation and mineralization.In this study the decrease of the pH in the solution was predicted via DOC conversion and the sole generation of H 2 CO 3 * by model A. As presented in Figure 5, DOC generation was progressively increasing as PHE was being degraded up to a certain point (117 min, corresponding to [DOC] max = 2.061 × 10 −3 M, and equivalent to ca. 93% of PHE degradation).From this point, DOC started to decrease until [DOC] f = 1.879 × 10 −3 M. That breakpoint represented the moment at which PHE was almost completely transformed into by-products.In addition, at this point, PHE mineralization began to be more evident, since H 2 CO 3 * level rose approximately in a linear way, up to a final level of 6.493 × 10 −14 M, with the subsequent pH decrease.Similarities between the pattern of this DOC profile and that of the formed intermediate curves reported in Alnaizy and Akgerman [19] research are highlighted.
On the other hand, it is worth noting that HO° level evolution with the reaction time was different when comparing the kinetic prediction model A or B with C. Obviating HO 2 °contribution to PHE degradation, Figure 6 shows that the highest final HO° level was achieved in the prediction model A, with a maximum HO° concentration equal to 9.435 × 10 −14 M.This value was similar to HO° final level in the prediction model B (9.400 × 10 −14 M) and different to that of the prediction model C, where HO° final concentration was 4.486 × 10 −14 M (ca.48% lower than the obtained in the kinetic model A or B).Approximately, in the first 40 min of the process a larger number of HO° in the aqueous medium with the developed kinetic model C was evidenced.Apparently, this amount of HO° was sufficient to degrade about 60% of PHE initial level under the studied experimental conditions.In contrast, models A and B, whose lines are overlapped, produced a lower number of HO° and the theoretical conversion of PHE remained above the experimental data.One possible reason for this discrepancy can be ascribed to DOM and dissolved oxygen positive effects in producing reactive oxygen species (ROS), as HO° [26], which were not considered.After 40 min of reaction, the HO° level was higher in the kinetic models A and B than in model C.This larger amount of HO° might lead to a faster conversion of the pollutant in comparison with the predicted model C, since the hypothetical depletion curve of PHE was below the measured data.Nevertheless, this rapid pollutant degradation did not occur actually.Therefore, there was an amount of HO° produced in excess that was not reacting with the contaminant.This surplus of HO° could be involved in free radical scavenging reactions.As model A and B consider the detrimental effect of HO° consuming reactions, it is suggested that their kinetic rate constants are higher than those ones used in this paper for these reactions to have a larger weight in the system.Additionally, the contradictory outcome between the actual situation and the theoretical one in the first and second stage of the process can also be attributed to the fact that just a fraction of the concentration of the species involved in the whole kinetic equations of the predicted models was actually reacting.Consequently, the real level of the species implicated in each kinetic reaction should be considered.However, it is rather difficult to determine which amount of the chemical species is exactly involving in each reaction for each time step, especially due to the high reactivity of radicals as the oxidation system progresses.In this context, further studies are required to overcome this limitation.
From these findings, it is suggested that there was an effective level of the formed HO°.Below that level, there was a lack of HO° for an efficient pollutant conversion; and above it, an excessive number of HO° was generated.That HO° effective level could be of relevance for industrial applications in order to be maintained throughout the reaction time, allowing an efficient pollutant degradation.

Estimation of PHE-HO 2 ° reaction rate constant
In order to study the action of HO 2 ° for pollutant degradation in the UV/H 2 O 2 system, PHE-HO 2 ° rate constant was calculated.For this purpose, the kinetic model A was used and it was estimated through a non-linear least-square objective function.The objective function for minimizing the error between the predicted and the measured data was defined as Eq. ( 28) [17].
( ) where  of PHE-HO 2 ° rate constant was in a stronger agreement with the experimental data (R 2 = 99.0%)than the previous predicted kinetic model A studied in section 3.2 (R 2 = 97.41%).The same conclusion can be drawn for kinetic model B and C, with a R 2 of 98.78 and 99.57% using the calculated PHE-HO 2 ° rate constant, in comparison with 97.37 and 99.34%, respectively.Therefore, in order to achieve a better fit between experimental and predicted data, HO 2 ° action in pollutant degradation should be considered.Nonetheless, HO 2 ° contribution to PHE degradation is non-significant in comparison with the role developed by HO° attack.

Conclusions
A kinetic model for studying pollutant degradation by the UV/H 2 O 2 system was developed, including the background matrix effect in scavenging free radicals and shielding UV-light and the reaction intermediate action, as well as the change of the pH as the UV/H 2 O 2 process proceeds.Three different ways for calculating HO° level time evolution were assumed (nonpseudo-steady, pseudo-steady and simplified pseudo-steady state; denoted as kinetic models A, B, and C, respectively).It was found that the assumption of pseudo-steady (simplified or not) or transient state for determining the HO° level evolution with time was not significant in PHE degradation rate due to the relatively low HO° level present into the bulk (~10 −14 M).
On the other hand, taking into account the high levels of HO 2 ° formed in the reaction solution compared to HO° concentration (~10 −7 M >>>> ~10 -14 M), HO 2 ° action in transforming PHE was considered.For this purpose, PHE-HO 2 ° reaction rate constant was calculated and estimated to be 1.6 × 10 3 M −1 s −1 , resulting in the range of data reported from literature.It was observed that, although including HO 2 ° action allowed slightly improving the kinetic model degree of fit, HO° developed the major role in PHE conversion, due to their high oxidation potential.
Additionally, it was found that there was an effective level of the HO° formed in solution.
Below that level, there was a lack of HO° for an efficient pollutant conversion; and above it, an excessive number of HO° was generated.That HO° effective level calculated from kinetic model C could be of relevance for industrial applications in order to be maintained throughout the reaction time, allowing an efficient pollutant degradation.
In this study, there was an attempt to contemplate a wide range of the chemical reactions involved in the UV/H 2 O 2 process and although high correlation factors were obtained, it is suggested to include the positive effect of the OM and the dissolved oxygen in generating ROS, as well as the effect of other anions naturally present in water bodies, as phosphate and nitrate, for the model to be a more accurate approximation of reality.

*
3 − and CO 3 2− in the studied water, the kinetic model accounts for the degradation of the pollutant by direct photolysis and oxidation through HO°, HO 2 °, O 2 °−, and CO 3 °−.As the model includes the reactions where Cl − , SO 4 2− , CO 3 2− , and HCO 3 − are involved, it can be used for the treatment of different types of water.Based on the reactions illustrated in Figures 1 and 2, and the different involved parameters, the mass balances and the corresponding ODE of the species of interest (C, DOC, H 2 O 2 , HO 2 -, HO°, HO 2 °, O 2 °−, CO 3 °−, SO 4 °−, HClO°−, H 2 ClO°, HClO°, Cl°, Cl 2 °−, OH − , H + , H 2 CO 3

Figure 1 .
Figure 1.Reaction network of the UV/H 2 O 2 process for the degradation of the parent compound without considering the effect of matrix constituents.The species above and below the arrows are the co-reactants and the co-products, respectively.

Figure 2 .
Figure 2. Reaction network describing the effect of the matrix background in the UV/H 2 O 2 system with co-reactants and co-products above and below the arrows, respectively.

Figure 3 .
Figure 3. (a) Reaction discrimination diagram as a function of the solution pH and (b) pH evolution of the medium during the UV/H 2 O 2 process in a completely mixed batch photoreactor.Operating conditions: H 2 O 2 /PHE = 495; [C] 0 = 2.23 × 10 −3 M; t = 220 min.

Figure 4
Figure4compares the simulation results of the proposed kinetic models A, B, and C with the experimental data for 2.23 × 10 −3 M PHE with a H 2 O 2 /PHE ratio of 495 and a reaction time of 220 min.The measured and simulated results for PHE direct photolysis alone are also presented.It is observed that more than 90% of the initial PHE concentration was removed after 220 min due to both direct photolysis and indirect degradation (primarily due to HO° attack).The effect of CO 3 °− could be seen as marginal because of the reduced number of those radicals, in the range of 10 −15 M, and the low reaction rate constant with PHE compared to HO°.In addition, it is shown that PHE was not completely removed by direct UV photolysis under the tested photochemical conditions (Table2), since approximately 50% of the total PHE degradation was attributed to UV photolysis, as it was experimentally determined by Alnaizy and Akgerman[19].On the other hand, the figure demonstrates that the prediction kinetic model C was in good agreement with the available experimental data with a relative high correlation factor (R 2 =99.34%).For prediction models A and B, R 2 were 97.41 and 97.37%, respectively.As a result of the subtle differences between kinetic models A and B, the depicted line describing model A overlaps model B line.An acceptable agreement between model predictions considering only the UV radiation and experimental data was also verified (R 2 = 98.25%).

7 .
[C] measured and [C] predicted correspond to the evolution of experimental and calculated pollutant concentration, respectively.This expression is a function of PHE-HO 2 ° rate constant.The optimum value for PHE-HO 2 ° second-order rate constant was found to be 1.6 × 10 3 M −1 s −1 , which is consistent with the range of the reported values by Kozmér et al. ((2.7±1.2) × 10 3 M −1 s −1 ) [29].The results of running the new prediction kinetic model A (with and without the contribution of HO 2 ° to pollutant conversion) and the experimental data are presented in Figure The figure shows that the prediction model A with the estimation