Process Optimization and Modeling by Response Surface Methodology of Nitrite Electro-Reduction by Ti/RuO2 + IrO2 Electrode ()
1. Introduction
Water is the essential element for all living things life [1] . Its availability and abundance play an important role in the development and evolution of societies [2] . Water pollution is a persistent problem threatening both human health and ecosystems [3] - [11] . Industrial development and excessive use of chemical fertilizers in agriculture contribute to the deterioration of its quality [12] [13] [14] [15] . In recent decades, water pollution caused by anions such as nitrites has become a major problems for researchers [16] [17] [18] . Nitrite is formed during the biodegradation of domestic or industrial nitrogenous wastes as well as some fertilizers. Several detection and treatment methods have been developed, such as spectroscopy, electrocoagulation, electrooxidation and electroreduction [19] . Electrochemical methods are an alternative because of their advantages, for example, environmental compatibility, versatility, high energy efficiency and profitability [20] [21] [22] . This work focuses on the optimization of parameters influencing the yield of water electrochemical denitrification by Ti/RuO2 + IrO2 electrode. To the best of our knowledge, this electrochemical denitrification method has never been applied to a Ti/RuO2 + IrO2 electrode. For this reasons, and with the aim of enhancing the electrochemical method, we examined the nitrite reduction on this choice of electrode at the anode as well as at the cathode. This study will investigate the influence of operating parameters such as electric current intensity, salt concentration and electrolysis time. Moreover, studies of process modeling, supported by the methodology of the response surfaces, lead to a real edification on the efficiency of the process.
2. Materials and Methods
2.1. Chemicals and Apparatuses
All electrochemical manipulations were performed using Ti/RuO2 + IrO2 electrode. Electrolyses were carried out in a standard two-electrode cell, with a Voltalab 40 potentiostat, connected to an interfaced computer that employed Voltamaster 4 software. The stock solution, with a concentration of 2 g/L, was prepared by dissolving NaNO2 solid (>99%, Sigma-Aldrich) in distilled water. The other solutions were obtained by successive dilutions up to the desired concentrations. UV-visible absorption spectroscopy was used to determine the
concentrations before and after electrolysis, using a Varian UV-vis spectrophotometer Cary-60. In order to control the pH, solutions of NaOH or HCl were added before starting the electrolysis. The pH values were measured with a HI 2211 Ph/ORP Meter pH-meter. For the determination of nitrite ions, the Griess reaction was used. Sulfanilamide (>98%, Sigma-Aldrich) and N-(1-naphthyl) ethylenediamine (>98%, Sigma-Aldrich) were used as received.
2.2. Methods
We first established the UV-vis calibration curve which was then used to evaluate the concentration of all solutions before and after electrolysis. For this purpose, we prepared 4 solutions of nitric acid with concentration ranges of 0.5, 1, 2.5 and 5 mg/L. Removal efficiency for nitrite reduction (R) was obtained by Equation (1):
(1)
where Ci and Cf are the initial and final concentration (mg-N/L) of
.
Table 1 represents the experimental factors to be optimized and modeled by response surface methodology using Disign-Expert version 13 software.
The matrix consisting of 4 factors, namely electrolysis time,
concentration, reduction current and pH, is used to draw up a plan of experiments comprising 30 tests including 6 in the center. This was based on Box-Behnken response surface designs by George Box and Donald Behnken.
The number of experiments N to be performed is given by the following relationship [23] :
(2)
where r is the number of replicates in the centre and k is the number of factors. Thus, the mathematical equation below allowed the calculation of the response according to the factors:
(3)
Y being the predicted response (the degree of conversion); b0 the constant coefficient; bi the linear coefficients; bij the interaction coefficients; bii the quadratic coefficients and xi, xj are the coded factors of the operating parameters.
The analysis of the different models studied (linear, 2FI, quadratic and cubic) shows that the quadratic model with a P = 0.0001 less than 0.05 is the most adequate model to predict the response with a correlation coefficient of 0.9664 according to the predicted R2 of 0.9220. For the study parameters, the quadratic model was proposed by Design-Expert 13.
3. Results and Discussion
3.1. Diagnostic Model
The results of the comparison between the values obtained by the model (predicted Y) and the experimental data (experimental Y) are summarized in Figure 1. In this case, the correlation between the theoretical response and the experimental one calculated by the model is satisfactory.
3.2. Numerical Optimization of Parameters
Despite its efficiency in determining the optimal parameters of a study, this optimization has some limitations, as the number of experiments to be carried out, and the impossibility of seeing the interaction of the various parameters affecting the yield. To overcome this, the response surface methodology is used. It allows maximum information with a minimum of experience.
In order to better understand the factors influencing this efficiency, we studied the effect of the following parameters: time, current intensity and nitrite concentration. The results obtained after electrolysis at temperatures between 22˚C and 28˚C were carried out at neutral pH. We first optimized the time, in the range of 10 to 180 min. During the reaction, the other parameters such as nitrite ion concentration (50 mg-N/L), current intensity (−400 mA), pH (7 - 8) and electrode surface area (20 cm2) remained invariant. In Figure 2, an increase in the efficiency is noted up to an electrolysis time of two hours. From this value, increasing the electrolysis time leads to a decrease in the removal efficiency. This may be due to the absorption of other compounds formed in these conditions. A long electrolysis time is also known to favor secondary reactions. In the course of the electrolysis, the solution becomes yellow, in agreement with this hypothesis. This result indicates an optimal time of 120 min. In the following, we will maintain the electrolysis time at two hours.
The effect of the initial concentration on the elimination yield was investigated at room temperature (27˚C) by fixing the electrolysis time at 120 min, the current intensity at −100 mA and the pH at 7 - 8. The results are shown in Figure 3. Although, the increase in the
concentration at high contents leads to a decrease in efficiency, Figure 3 shows that the addition of increasing amounts of nitrites leads to an improvement in efficiency at values below 30 mg-N/L. Beyond this maximum value, any increase in concentration becomes unfavorable to the process. After electrolysis at 30 mg-N/L, we obtained, after analysis by UV-visible spectroscopy, an elimination yield of 97%.
![]()
Table 1. Process factors and their levels.
![]()
Figure 1. Experimental and predicted response for nitrite removal.
![]()
Figure 2. Influence of time on the removal efficiency.
![]()
Figure 3. Influence of concentration on the removal efficiency.
For any wastewater treatment process, the applied current intensity is one of the most important parameter. The current intensity was varied between −50 to −700 mA to test a wide range of reducing current. The effect of the current, as a function of the removal efficiency, is represented in (b). Electrolysis was carried out in the following conditions: current intensity = −100 mA, pH = 7 and electrolysis time = 120 min. From Figure 4, it appears that the nitrite removal efficiency approaches a limiting value close to zero for a wide range of current. This result could be explained by the use of two Ti/RuO2 + IrO2 electrodes in a single compartment, allowing a higher electron production and facilitating nitrite reduction.
Initial pH is an important parameter for wastewater treatment. Starting with 30 mg-N/L of NaNO2, the current was set at −100 mA. Figure 5 shows the degree of conversion in the pH range of 2.0 - 14. In contrast to the situation observed in acid medium, an efficiency of 100% was observed for nitrite reduction. In the literature, it was notably established that the protons are reduced near 0.00 V/SHE in our conditions, all possibly contributing to limiting the faradaic yield of the desired reaction. The results of this work consistently indicate that the ideal pH range for the reduction of nitrite ions is located in the basic range.
3.3. Analysis of Variance
For the choice of model, the Desirability Function Approach (DFA) was used and applied with the help of the Design Expert software to overcome this ambiguity in order to examine the optimal operating conditions for the electro-reduction of nitrite in the chosen experimental field. In order to help with the optimal parameters, a graphical visualization over the entire design space seems more interesting to us than the values of the design parameters that maximize desirability. Table 2 summarizes the parameters for response modeling and optimization by response surface methodology. Analysis of variance is used to describe the significance of the curvature of response at a 95% confidence interval. Indeed, significant process variables are usually decided on the basis of the F-value or P-value [24] [25] . The F-value of the model of 60.6 implies that the model is significant. There is only a 0.01% chance that such a large F-value is due to noise. P-values below 0.05 also indicate that the model terms are significant.
![]()
Table 2. Analysis of variance (ANOVA) for the quadratic response surface model.
DF: Degree of freedom.
![]()
Figure 4. Influence of current on the removal efficiency.
![]()
Figure 5. Influence of pH on the removal efficiency.
Table 2 shows that only the cross-effects AD is significant on the response (the degree of conversion). They suggest that the neutral pH was more favorable for the removal of nitrites.
The result of electrolysis in the same condition shows that the predicted value (100%) is very close to the experimental one (99.32%). This result highlights a valid and applicable model for the prediction of the response.
4. Conclusion
In this work, it was shown that optimization of the operating parameters such as NaNO2 concentration, electrolysis time and current intensity significantly improves the efficiency of the process. The necessity to improve the nitrite removal process led us to parameters modeling (reaction time, pH, NaNO2 concentration and current intensity). In this study, we were able to show the significant interaction between pH and duration of electrolysis. The results indicate a real correlation between the optimal value (100) predicted by the mathematical model and that obtained experimentally (99.32). In this study, the choice of working electrode, Ti/RuO2 + IrO2 is justified on several levels but it would be interesting to locate the trend of other electrodes and guide further research, aiming at the design of an electrode on this process.